Custom mometum source term / modifying solidification drag eqn
 

Hi,
In Flow3D the solidification flow attenuation is undertaken by ramping the viscosity and applying a drag term. Details in the mannual.

How can I modify the equation used for the drag term or create my own momentum source term to replace the drag function?

I'd be happy to turn off the solidification model and create my own by controlling viscosity, latent heat, and applying a momentum source term. I just can't work out how to create a momentum source term dependent on solid fraction/temperature.


- I can modify the viscosity ramp via the mucal.F90 user function
- I can implement the latent heat via modifying the Cp curve
- I can't currently work out how to modify the drag term, or recreate it as a sub-routine
- I see there are mass source user sub-routine, but I don't see a momentum one


Any thoughts or suggestions?