[mikulo]

Originally Posted by janghel

Anyways I got it working for mushy zone ( both linear and schiel law case)

Hello, would you mind if you can share your solver and a test case? I will greatly appreciate it. Thanks!

[gtarang]

Originally Posted by mikulo

Hi, thank you for your shared case files and solver. May I ask in which solver will it work? It would have been nice if it will work on OF 7. I will try it later..

I ran it with Openfoam 6, so I guess it will work with openfoam 7 too. But I doubt whether it will run with openfoam8. Feel free to post any errors you get.


-
gtarang

[gtarang]

Originally Posted by mikulo

Hi, thank you for your shared case files and solver. May I ask in which solver will it work? It would have been nice if it will work on OF 7. I will try it later..

Hi,

I just saw the message. I have used this with Openfoam 6.

Thanks

[mikulo]

Originally Posted by blebon

If you update melt_fraction using the if statement, energy is not conserved, because solution of the energy balance equation is based on previous values of melt fraction and hence inconsistent with the current temperature. If I am not mistaken, the correction process (i.e. what is implemented in solidificationMeltingSource) is based under the assumption that enthalpies is the same both for solving the energy balance equation and updating the melt fraction. Upon convergence of the iterative process, both energy is conserved and melt_fraction residuals should be very small.

When you have a mushy zone, you need to replace Tmelt_ but the temperature corresponding to alpha1_[celli] (either using a formula -- lever rule, Scheil or whatever --) or a lookup table.

Hello, the way I understood the code, the alpha is just computed once for every time step. How could you say that the alpha residuals are computed and should be very small?

-Mike