Phase Field Model Units

While phase field models are often nondimensionalized, the inherent equations are dimensional. Thus, it is critical to understand and correctly implement the units, especially when coupling multiple physics.

Units Analysis of Phase Field Equations

In the unit analysis, we break the equations down by units, with the following conventions:

  • - length

  • - time

  • - energy

  • - moles

Free Energy Functional

The driving force for evolution in the phase field method is the minimization of a free energy functional. The free energy functional is written as

where defines the local free energy density with the units and is a function of concentration variables (typically a molar fraction) and order parameters . The gradient energy density

where and are gradient energy coefficients with the units . Finally, describes any additional sources of energy in the system, such as deformation or electrostatic energy, with units of .

Allen-Cahn Equation

The Allen-Cahn equation after the variational derivative takes the form

The units of this equation are

where the units of are .

Cahn-Hilliard Equation

The Cahn-Hilliard equation after the variational derivative takes the form

where is the molar volume of the reference state of the material with units of . The units of this equation are

where the units of are . Note that some models include the in the mobility term, such that it has units of .

CALPHAD Local Free Energies

A common approach for phase field models is to use free energies for the individual phases taken from the CALPHAD approach. CALPHAD free energies typically have units of , therefore they must be converted to free energy densities using the molar volume according to

Note that when using a multiphase model, the CALPHAD free energies of each phase must be converted with the same reference molar volume , with changes in volume between phases handled by stress-free strains.

Handling Units in MOOSE

There is no inherent unit system in MOOSE. Thus, the units of the phase field equations are set by the user when they define a model. Specifically, the units are set by the local free energy density and the and mobility parameters ( and ). The units in all these terms must be consistent. Additional energy sources, such as the elastic energy, also must have consistent units. In the phase field module, all of these values are created using Material objects. Thus, the units of your system are not set by the kernels but rather by the materials.

One useful practice is to create your material objects to take SI units as input parameters. Then use length_scale, time_scale, and energy_scale input parameters to convert the actual units of the problem. As an example of this, see the PFParamsPolyFreeEnergy material, where the input file block looks like


[./Copper]
  type = PFParamsPolyFreeEnergy
  block = 0
  c = c
  T = 1000 # K
  int_width = 30.0
  length_scale = 1.0e-9
  time_scale = 1.0e-9
  D0 = 3.1e-5 # m^2/s, from Brown1980
  Em = 0.71 # in eV, from Balluffi1978 Table 2
  Ef = 1.28 # in eV, from Balluffi1978 Table 2
  surface_energy = 0.7 # J/m^2
[../]

Taken from phase_field/tests/PolynomialFreeEnergy/split_order4_test.i