[wuxuetianzi]

Problem description：
This is a melting simulation with an ice hemisphere siting in the center of a cylinder container (aluminum), the rest of container filled with air. the geometry picture is attached. The second picture indicates that I will heat the wall of container, including bottom, up and keep temperature at 350K.
Here is my setting:
General:
1. Mesh seems to be ok with positive values for “minimum volume” and “minimum face area”.
2. Solver Pressure Based, Time unsteady, Gravity checked (-9.81 in Y direction), 2D-space Planner.
Models:
1. Multiphase model (VOF) with 2 eulerian phase. I have activated the check mark of “Explicit” for volume fraction parameters and “implicit body force” for body force formulation.
2.Solidification/melting is ON (so as the energy equation)
Phases:
air -primary, and ice- secondary.
Material:
add ice from fluent database and set solidus temperature at 273K; liquidus temperature at 274K; pure solvent melting heat at 330000 j/Kg, Cp at 2100 j/kg-k and keep the rest as default.
Cell zone condition:
Zones are air, ice and wall and wall type is solid
Boundary condition:
1.only “heat” part set with temperature at 350K, which is shown in the picture
2.Both fluid zones are automatically set as “interior”.
Solution method:
Using default setting.
Solution control: Solution control is left untouched.
initialization: Standard. Initial temperature is 273K, patched ice at 272K, air at 350K, I did not patch volume fraction value for ice because once I did that, everything went wrong, I don't know why.
Then, I run the simulation by time step size of 0.05s, number of time steps is 100

Finally, I got liquid fraction result as shown in the picture, which is not I expect. What I want is a shape changing ice like real meilting.


Thanks a lot!

[wuxuetianzi]

could anyone help me？please

[wuxuetianzi]

could anyone help me？Thank you so much

[wuxuetianzi]

could anyone help me？Thank you so much