underworldcode / underworld2

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

Compressible mantle convection with tala approximation #688

Open gp37 opened 4 months ago

gp37 commented 4 months ago

Hi,

I saw that Underworld 3 can do compressible mantle convection using the TALA approximation: https://www.underworldcode.org/compressible-convection-in-cartesian-coordinates-in-underworld3/

Does this mean that Underworld 2 cannot do compressible convection because it cannot change the conservation of mass assumption that div(v) = 0 (incompressible) to div(v) + 1/rho u div(rho) = 0 (compressible)? where v is velocity and rho is density.

Are there any workarounds for this in underworld 2?

NTIaN7 commented 4 months ago

As far as I know, it seems that ALA or TALA approximation cannot be well implemented in underworld2. For similar problems, you can refer to the previous discussion.