Defining the Liquid Phase for the Melting-Solidification Model
You define the liquid phase and apply the Melting-Solidification model.
You model the freezing water as a Eulerian phase in Simcenter STAR-CCM+. The properties of water come from the material database, with some modifications to suit the problem. For this tutorial, the value of latent heat is reduced so that freezing takes place sooner than it would with the default value. For a real case, however, leave the value at its original setting.
To define the liquid phase for the Melting-Solidification model:
- In the Physics 1 continuum, right-click the node and create a new phase.
- Rename the Phase 1 node to H2O.
-
For the
H2O phase, select the following models in order:
Group Box
Model
Enabled Models
Energy (selected automatically)
Turbulent (selected automatically)
Material
Liquid
Equation of State
Constant Density
Optional Models
Melting-Solidification
- In the H2O Model Selection dialog, in the left pane, select the Physics 1 node.
- In the Physics 1 Model Selection dialog, in the Optional Models group box, select Melting-Solidification Flow Stop.
- Click Close.
- Select the 0.1. node and set Under-Relaxation Factor to
You set the solid volume fraction threshold above which flow is stopped in the fluid phase.
- Select the 0.999. node and set Value to
The thermal conductivity of H2O varies with phase and temperature. For this simulation, assume that the thermal conductivity is constant for each phase.
The thermal conductivity of ice is 2.33 W / m K while the conductivity of H2O is 0.620271 W/m K. Use a conditional expression to set the conductivity of the H2O phase. Set the conductivity to 2.33 W/m K when the temperature is less than or equal to 273 K, and 0.620271 W / m K for higher temperatures.
Create a step function to set the conductivity of H2O:
- Right-click the node and select from the pop-up menu.
- Rename the Variable Thermal Conductivity. to
-
Select the
Variable Thermal Conductivity node and set the following properties:
Property Setting Function Name Variable Thermal Conductivity Definition (${Temperature} <= 273.) ? 2.33 : 0.620271 -
Set
Dimensions to
Power/Length-Temperature by filling the
Dimensions dialog with the following values:
The Properties window appears as shown.
-
Edit the
node and set the following properties:
Node Property Setting Thermal Conductivity Method Field Function Scalar Function Variable Thermal Conductivity Value 100000.0 J/kg Value 273.0 K Value 273.0 K - Save the simulation.