How to combine Multi-Fluid VOF with Geo-Reconstruct and energy source in a phase change process?
my goal is to solve 1D Stefan problem on a mesh consisted of 50 cells (3 cells of superheated water-vapor adjacent to 10 K superheated wall and the rest is saturated water-liquid; linear temperature profile in the vapor cells) using Multi-fluid VOF model (sharp interface modeling) available in ANSYS Fluent 2019 R3. The most promising results is obtained with explicit VOF formulation, but the problem is that the slope of interface (alpha_g = 0.5) evolution is not in a good agreement with the analytical solution as shown in the attached figure.
The mass transfer is defined in accordance to Mr. Amine's recommendation (ttps://studentcommunity.ansys.com/thread/zone-of-applicability--define_mass_transfer/?order=all#comment-7a51709b-3091-4657-8157-aa4d00d4cf4d) in a sense that involves volume fraction gradient and reads
m_lg = (Kg * grad(T_g) * grad(alpha_g))/h_lg,
where "Kg" is thermal conductivity of vapor, "grad(T_g)" is temperature gradient of the vapor, "grad(alpha_g)" is gradient of vapor volume fraction and "h_lg" is the latent heat.
This equation is included in "DEFINE_MASS_TRANSFER" macro, while the gradient multiplication is realized using UDM inside "DEFINE_ADJUST" macro.
Also, the energy equation of each phase is supplemented with necessary energy source term (S_e.v = m_lg * h_lg; S_e.l = - S_e.v).
Together with the aforementioned mass transfer model, the available phase interaction terms are used: drag ("anisotropic-drag"), heat (Nu = 1) and interfacial area ("ia-gradient").
The applied discretisation schemes are as follows: Phase Coupled SIMPLE (pressure-velocity coupling), Least Squares Cell Based (gradient), PRESTO! (pressure), First Order Upwind (momentum and energy), Geo-Reconstruct (volume fraction) and First Order Implicit (transient).
I tried with other interface reconstruction approaches, but Geo-Reconstruct produced the best result.
Finally, my question is: "How can I improve the solution process in order to have smooth distribution of the interface position?"
Many thanks in advance!