search

underworldcode / underworld2

underworld2: A parallel, particle-in-cell, finite element code for Geodynamics.
http://www.underworldcode.org/
Other
174 stars 59 forks source link

Uncoupled temperature and compositional field in UW Geodynamics #679

Open bush2014 opened 1 year ago

bush2014 commented 1 year ago

I am currently mentoring a high school student using UW Geodynamics, and I must say it's been great for educational purposes; we really appreciate the effort put into it.

In this project, where we've imposed a velocity (1 cm/year) from one side of the wall. However, we've noticed that the temperature of the slabs and the composition don't seem to match well when the slabs are being subducted. In my experience with UW, we usually address this issue with the following comment: dt = min(advDiff.get_max_dt(), advector.get_max_dt()). Is there a similar approach or workaround in UW Geodynamics to address this temperature and composition mismatch? We are using the version underworld2_2.13.1b, which combines UW and UW Geodynamics.

Any information or guidance on this matter would be greatly appreciated. Thanks!

julesghub commented 1 year ago

Hi @bush2014,

Thanks for the great feedback - it's good to know UWGeodynamics is being used for educational purposes. Are you using a model from the repository or do you have picture to illustrate the issue. That would help lots in discussing what's occurring.

UWGeodynamics by default applies a CFL condition as is used generally in Underworld but you can override it using this idea https://underworld2.readthedocs.io/en/latest/UWGeoUserGuide.html#specify-a-timestep

Does this help?

bush2014 commented 1 year ago

Hello Julian, I really appreciate your prompt response! This is incredibly helpful. I believe this is exactly what I was seeking, and I will inform you if it doesn't work out. Thanks!