search

geodynamics / aspect

A parallel, extensible finite element code to simulate convection in both 2D and 3D models.
https://aspect.geodynamics.org/
Other
217 stars 232 forks source link

Add fluid velocity postprocessor #5902

Open danieldouglas92 opened 2 weeks ago

danieldouglas92 commented 2 weeks ago

Currently, the fluid velocity can only be output when Include melt transport = true. However, a fluid phase can exist and be advected with the darcy velocity if the compositional field method is set to darcy field. This PR allows the user to output the fluid velocity when simulating fluid migrating using the darcy field method, and not just when Include melt transport = true. This PR builds off of #5856, and so right now it fails because MeltOutputs are not being filled and created by the Reactive Fluid Transport Model, and additionally the test made here is essentially the same test as well, so #5856 should be merged first.

This addresses one of the points in #4748.

danieldouglas92 commented 2 weeks ago

I think this is good for another run through when you get a chance @jdannberg

danieldouglas92 commented 23 hours ago

This is ready for another look through @gassmoeller! I added the test for this functionality into the existing test for the darcy_convection_timestep