cimlib\tests\modele\Binary_DirectionalSolidification\Enthalpie_gsTabule

Note:
Changes to do to transform the test case to an Exercice and obtain outputs corresponding to curves below:
In Maillage.mtc: change the value of TemperatureInitialeC from 622 to 800
In Increment.mtc: change the value of TempsFin from 3 to 1000
Results of reference for the long simulation are available in the test folder

Table of caracteristics of the model

In this model, the solid fraction is linked to the temperature through a tabulation. For temperatures is the mushy zone, this tabulation corresponds to the Gulliver-Scheil model ($f_S = 1- \left(\frac{T_M-T}{T_M-T_L}\right)^{1/(k-1)} $ where $k$ is the segregation coefficient; $T_M$ is the melting temperature and $T_L$ is the liquidus temperature corresponding to the nominal composition). For temperatures under the eutectic temperature, the solid fraction is equal to $1$.

In the simulation, the temperature field and the solid fraction field are recorded on sensors which are placed every 2 cm along the z direction. In the following figures, solid lines correspond to outputs of the simulation on these sensors and dotted lines correspond to results of the fornt tracking model presented on page Front tracking 1D (for interfaces with low undercoolings).

Outputs on the temperature field on sensors

Outputs on the solid fraction field on sensors

Results of the model are in good agreement with the ones of the front tracking model. For the last two sensors, the transition from mushy to eutectic is longer than for other sensors. This is because the decrease of temperature is very slow in this range of time and in this part of the domain.