First Order Phase Transition Material

First Order Phase Transition

This exploratory MPM material models a first order phase transition between two materials (it is currently only available in OSParticulas). It has no material response itself. Instead, it acts as a "parent" material to two "child" materials, where one child is a "solid phase" and the other is a "liquid" phase. The properties of this material control transition from the solid to liquid phase, where solid is the low-temperature phase, and the liquid is the high temperature phase.

A first order phase transition occurs at a melting temperature, T_melt. The meaning of a first order transition is that thermodynamic energy functions (e.g., Gibbs free energy) are continuous at T_melt, but the phase which is the lowest energy changes. Although free energy is constant, the first derivatives of free energy undergo a discrete change (hence the origin of the term first-order transition). Namely, there is a change in enthalpy (ΔH_f), entropy (ΔS_f), and volume (ΔV_f), where subscript "f" means fusion. Because free energies of the two phases are equal at T_melt, the first two are related by

      [math]\displaystyle{ \Delta G(T_m) = 0 = \Delta H_f - T_m \Delta S_f \qquad{\rm or}\qquad \Delta S_f = {\Delta H_f\over T_m} }[/math]

Enthalpy and Volume Changes in the Transition

Most materials do not have a sharp transition, but rather undergo a transition over a range of temperatures, such as from T_i to T_f. The enthalpy of fusion is experimentally related to the heat capacity during this transition by:

      [math]\displaystyle{ \Delta H_f = \int_{T_i}^{T_f} \Delta C\ dT }[/math]

where ΔC is the excess heat capacity during the transition compared to the heat capacity of the material in the absence of a transition. Likewise, the density of the liquid will differ from the density of the solid, which means there is a volume change of fusion. The volume change of fusion can be expressed from thermal expansion (or shrinkage) during the transition of

      [math]\displaystyle{ \Delta V_f = \int_{T_i}^{T_f} \Delta \alpha\ dT }[/math]

where Δα is the excess thermal expansion during the transition compared to the thermal expansion of the material in the absence of a transition.

Modeling with a Finite Width Transition

One approach to modeling a first order transition is to spread out the transition region and allow [math]\displaystyle{ \Delta C }[/math] and [math]\displaystyle{ \Delta \alpha }[/math] to depend on temperature within that region. We can define [math]\displaystyle{ f_L }[/math] as the fracture of the material that is liquid where [math]\displaystyle{ f_L }[/math] transitions from 0 in the solid phase to 1 in the liquid phase. The excess heat capacity and thermal expansion can then be written as:

      [math]\displaystyle{ \Delta C = \Delta H_f \frac{df_L}{dT} \qquad{\rm and}\qquad \Delta \alpha = \Delta V_f \frac{df_L}{dT} }[/math]

This material can use several different f_L functions:

• Hat (shape 0): This option uses a piece-wise quadratic function for f_L such that ΔC is a triangular "hat" function from T_melt-ΔT to T_melt+ΔT with peak ΔH_f/ΔT at T_melt.
• Box (shape 1): This option sets f_L to a linear function that varies from 0 to 1 over the range T_melt-ΔT to T_melt+ΔT. The resulting ΔC is a constant value of the same range of ΔH_f/(2ΔT).
• Sech^2 (shape 2): This option sets f_L to an exponential step function:
        [math]\displaystyle{ f_L = \frac{1}{2}\left[1-\tanh\left(\frac{T-T_m}{\Delta T}\right)\right] }[/math]
  The resulting ΔC is a sech^2 function with peak equal to ΔH_f/(2ΔT) at T_melt. Here ΔT is a parameter characterizing the width of the transition
• Smooth Step (shape 3): This option sets f_L to an cubic step function from 0 to 1 with zero slopes at T_melt-ΔT to T_melt+ΔT.

For numerical modeling, the you enter the total heat of fusion, f_L shape style, and transition width parameter ΔT. The same transition function is used for Δα where ΔV_f is found from density difference between solid and liquid phases. Because implementation of heat of fusion is done through heat capacity, all simulations with phase transition materials should activate coupled conduction calculations (conduction is also needed to see changes in temperature that can lead to phase changes). The specific values ΔT and shape style may need to varied for numerical stability in specific problems.

Transition Temperature Range

The temperature range for which a material undergoes a transition is commonly caused by spatial heterogeneity in the object. For example, polymer materials (which have a particularly large temperature range for first order transitions) have wide variation in crystal perfection and therefore a wide temperature range for melting. The ΔT property mentioned above is one way to model a temperature range, but that parameter is more associated with numerical implementation of implementing an infinite heat capacity or infinite thermal expansion that is seen for a sharp transitions. A second approach to modeling a transition temperature range is to allow T_melt to be a random variable representing heterogeneity in the crystal structures. This form of variation is modeled by entering the standard deviation for the melting temperature. It can be picked based on observations of the width of the melting transition.

The current implementation randomly assigns a normal distribution of melting points to all particles with zero spatial correlation. In the future, spatial correlation will be added.

Material Properties

The unusual task for this material is to use multiple Gk and tauk properties (all with the same property name) to enter a material with multiple relaxation times.

Property	Description	Units	Default
SolidPhase	Enter the solid phase material by ID (scripted files only) or by number	none	none
LiquidPhase	Enter the liquid phase material by ID (scripted files only) or by number	none	none
Tmelt	The transition temperature (K). The solid phase is stable below Tmelt and liquid phase above it	K	none
Tsigma	The standard deviation for the transition temperature (K). The temperature range will be limited to +/- three standard deviations	K	0
DeltaT	A temperature range for application of increased heat capacity and increased thermal expansion (must be > 0)	K	10
Hfusion	The heat of fusion for the transition	heat of fusion units	0

Note that the materials that are capable of being phase materials are limited to those that have been revised to work with this phase transition material. The currently allowed materials are Mooney Material, Neo-Hookean Material, Isotropic, Hyperelastic-Plastic Material, Isotropic, Hyperelastic-Plastic Mie-Grüneisen Material, Clamped Neo-Hookean Material, and Tait Liquid Material.

History Variables

This material tracks two history variables:

1. The currently active phase with 0 for solid phase and 1 for liquid phase.
2. The melting point of the particle. It is constant throughout the analysis, but can be visualized as random variable if Tsigma>0.

For history variable n>2, the output is history variable n-2 of the solid phase material. The liquid phase may also have history variables, but currently those variables can not be archived.