7.17.1.10.15. Melting

Syntax

Source For Cvfem Energy ([{of} SpeciesName] | [{in} MaterialPhaseName] | [{ls} {a | b | c}]) {@ | at | for | on | over} MeshPartIdentifier = Melting ([Ts = Ts] | [Tl = Tl] | [Latent_Heat = latent_heat] | [Profile_Shape = profile_shape] | [Add_Advection = add_advection])

Summary

Energy source term to represent the latent heat of melting of a material.

Description

Uses a Gaussian (default) or Square distribution based on solidus (TS) and liquidus (TL) temperatures to compute a liquid fraction $f_{L}$ . The Gaussian option is set so that only one percent of the energy gets applied outside of this range. The source term is computed as:

$\rho L \frac{f_{L}^{n+1}-f_L{^n}}{\Delta t}$

Uses either the LATENT_HEAT ( $L$ ) expression or a provided “latent_heat = X” argument for each source. If both are present, the provided argument is used. This also assumes that the LATENT_HEAT is provided per mass (this expression multiplies it by the density). The solidus and liquidus temperatures can either be provided inline or in the material model. If both are provided, the inline values are used. The ADD_ADVECTION option adds an additional term that tracks the advection of the solidification front, which is currently implemented as:

$\rho L\frac{\partial f_{L}}{\partial T}\boldsymbol{u}\cdot\nabla T$

If this option is enabled, users should be aware that this term may be subject to temporal errors where large enough time steps can advance the solidification front over and incorrectly withhold energy over many cells.

Parameter	Value	Default
{of}	{of \| species \| subindex}	–
SpeciesName	string	–
{in}	{in \| material_phase}	–
MaterialPhaseName	string	–
{ls}	{levelset_phase \| ls}	–
MeshPartIdentifier	string	–
Ts	real	0
Tl	real	0
latent_heat	real	0
profile_shape	“string”	gaussian
add_advection	integer	–