7.16.1.20.43. Melting

Syntax

Source For Energy [{of} SpeciesName | {in} MaterialPhaseName | {ls} {a | b | c}] {@ | at | for | on | over} MeshPartIdentifier = Melting [Power_Output = power_output | Toggle = toggle | 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.

The following options are present in most boundary conditions/sources. The meaning of each option is as follows:

Power_Output

If specified, the total (integrated) power of the given BC/Src is calculated and stored into the global variable name specified as the power_output argument.

Toggle

If given, specifies which toggle block controls whether this BC is active or not.

Parameter

Value

Default

{of}

{of | species | subindex}

–

SpeciesName

string

–

{in}

{in | material_phase}

–

MaterialPhaseName

string

–

{ls}

{levelset_phase | ls}

–

MeshPartIdentifier

string

–

power_output

“string”

–

toggle

“string”

–

Ts

real

0

Tl

real

0

latent_heat

real

0

profile_shape

“string”

gaussian

add_advection

integer

–