# Presentation: 036325

## Surface Energy and Surface Stress in Phase-Field Models of Elasticity

- J. Slutsker, G. McFadden, J. Warren, W. Boettinger, (NIST)

- K. Thornton, A. Roytburd, P. Voorhees, (U Mich, U Md, NWU)

- Surface Energy and Surface Stress in Phase-Field Models of Elasticity

- Surface excess quantities and phase-field models
- 1-D Elastic equilibrium – axial stress & biaxial strain
- 3-D Equilibrium of two-phase spherical systems

- Goal: illuminate phase-field description of surface energy and surface strain by simple examples

##     Surface Excess Quantities (Gibbs)

-     _Surface Excess Quantities (Gibbs)_

##     Kramer’s Potential (fluid system)

-     _Kramer’s Potential (fluid system)_

- (surface energy)

##     1-D Elastic System (single component)

- z

- Solid

- “Liquid”

-     _1-D Elastic System (single component)_

##     “Kramer’s Potential” (elastic system)

-     “_Kramer’s Potential” (elastic system)_

## Planar Geometry

- Solid and “liquid” separated by an interface
- Planar geometry
- No dynamics
- Applied uniaxial stress or biaxial strain
- 1D problem

***0*

- z

- Solid

- Liquid

- Examine
-  Equilibrium temperature (*T**0*)
-  Surface energy and surface stress (Gibbs adsorption)
- Analytical results and numerical results are compared

*e**S*

##     Phase-Field Model of Elasticity

-     Phase-Field Model of Elasticity

##     1-D Phase-Field Solution

-     _1-D Phase-Field Solution_

##     1-D Stress and Strain Fields

-     _1-D Stress and Strain Fields_

## Analytical Results: Melting Temperature

- Here is the analytical results for the melting temperature.  Using the first integral of the system, we obtain the constant quantity that gives equilibrium.  Here, f is the free energy including the temperature effect, fel is the elastic energy density, this is the gradient energy, and the last term is the work term due to applied stress, where tau0 is the applied stress, U’ is the strain ezz.  Using analytic solution for the elastic equilibrium, we obtain this expression for the melting temperature shift due to the presence of applied stress or applied strain.  Note that the shift depends only on the jump of the quantities, and does not depend on the interfacial thickness.

## Numerical Simulation: Melting Temperature

-  “Physical” parameters for Aluminum eutectic is used
-  Variables are non-dimensionalized using the latent heat per unit volume and the system length
-  Here, we focus on applied stress with no misfit: * *

- For numerical simulation, we take a physical parameter for Aluminum eutectic (or an estimate for it.  The equations are non-dimensionalized using the latent heat per unit volume and the system size.  We now focus on the case with applied stress only, that is tau0 is non zero, and msifit and applied strain is zero.  Then the expression for the melting temperautre change becomes this.  The change is proportional to the applied stress^2.

## Simulation and analytics agree

- Indeed the simulation results agree with the analytical results.

## Analytical Results: Surface Energy

## Numerical and analytical results agree

- Here, I plot dgama dtau0 agaist tau0.  The numerical solution agrees with the analytical solution for dgamma dtau0.

## Elastic Equilibrium of a Spherical Inclusion

- Bulk modulus, KL=KS=K

- Shear modulus, =0 in “liquid”

- VS<VL
- Self-strain: jk  in liquid
-                     0      in solid

- R1

- R

- f=fS-fL= LV (T-T0)/T0

- (1)

- (2)

- Compare phase-field & sharp interface results for Claussius-Clapyron/Gibbs-Thomson effects [numerics & asymptotics]
-                                                                                      [Johnson (2001)]

- _Elastic Equilibrium of a Spherical Inclusion_

## Phase-Field Model

## Sharp-Interface Model

## Interface Conditions

## Solid Inclusion

- L

- S

## Liquid Inclusion

- L

- S

## Phase-Field Calculations

- T/T0

- Liquid fraction

- S

- L

- _Phase-Field Calculations_

- Liquid-Solid volume mismatch produces stress and alters equilibrium temperature (Claussius-Clapyron)

## Phase Field vs Sharp Interface (no surface energy)

- _Phase Field vs Sharp Interface (no surface energy)_

- Liquid fraction

## Phase Field vs Sharp Interface (surface energy fit)

- _Phase Field vs Sharp Interface (surface energy fit)_

- Liquid fraction

## Conclusions

- Future Work

-  Phase-field models provide natural surface excess quantities
-  Surface stress is included – but sensitive to interpolation through the interface
-  Surface energy and Clausius-Clapyron effects included

-  More detailed numerical evaluation of surface stress in 3-D
-  Derive formal sharp-interface limit of phase-field model

## (End)