!* 2D example showing the effect of lattice parameter induced strain on particle growth


SET_NUMBER_OF_THREADS
4!Number of threads /1/

USER_CONSTANTS
!*User constants may not be altered once entered
!*Constant names must start with UC_
!*Names will be converted to upper case
uc_deltaE!Enter name of constant or EXIT to exit
0.107!Enter constant value
uc_gradeps!Enter name of constant or EXIT to exit
-2.625e-10!Enter constant value
uc_Mphi!Enter name of constant or EXIT to exit
0.05!Enter constant value
uc_Lalpha!Enter name of constant or EXIT to exit
-10e3!Enter constant value
uc_Lbeta!Enter name of constant or EXIT to exit
-20e3!Enter constant value
uc_mobab!Enter name of constant or EXIT to exit
1e-16!Enter constant value
uc_youngs!Enter name of constant or EXIT to exit
200e9!Enter constant value
uc_poissons!Enter name of constant or EXIT to exit
0.3!Enter constant value
uc_a_alpha!Enter name of constant or EXIT to exit
1.015!Enter constant value
uc_a_beta!Enter name of constant or EXIT to exit
1.0!Enter constant value
exit!Enter name of constant or EXIT to exit


MANDATORY_DATA
analytical!GES file / ANALYTICAL [<analytic init file>]
!Number of elements (>0) /2/
a!Element name
b!Element name
2!Number of phases (>0) /1/
alpha!Phase name
beta!Phase name
!(Initial) Temperature [K] /1000.0000/
!Pressure [Pa] /100000.00/
!Substitutional (A) /YES/
!Substitutional (B) /YES/
!Initially homogeneous composition /YES/
0.55!Initial mole fraction (A)
0.45!Initial mole fraction (B)
!State of phase (ALPHA) /1/
!State of phase (BETA) /1/
2!Dimensionality /1/
128!Number of grid points /10/
128!Number of grid points /10/
1.28e-5!Domain width [m] /1.0000000/
1.28e-5!Domain width [m] /1.0000000/
!Periodic boundary /NO/
!Periodic boundary /NO/
!Grid type /1/
!Grid type /1/
15!Simulation time [s] /1.0000000/
!Maximum timestep [s] /1.5000000/
!Minimum timestep [s] /0.10000000E-09/
!Initial timestep [s] /0.10000000E-09/
!Save interval /1/
0.1!Minimum elapsed time between saves [s] /1.0000000/

PHASE_FIELD_DATA
!Use phase-field variables /YES/
!Total nr of phase-field variables /2/
!File containing initial phase-field variable values
1!Index of associated phase of phase-field variable # 1
2!Index of associated phase of phase-field variable # 2
gs11!Gm expression for g_1
uc_deltaE*gs21!Gm expression for g_2
uc_gradeps!E0102
uc_Mphi!Enter phase-field mobility expression

ANALYTICAL_EXPRESSIONS
10E3*XA+20E3*XB+R*T*(XA*LOG(XA)+XB*LOG(XB))+uc_Lalpha*XA*XB!Gm for ALPHA
20E3*XA+10E3*XB+R*T*(XA*LOG(XA)+XB*LOG(XB))+uc_Lbeta*XA*XB!Gm for BETA
uc_mobab!Mobility for A in ALPHA
uc_mobab!Mobility for B in ALPHA
uc_mobab!Mobility for A in BETA
uc_mobab!Mobility for B in BETA
!Will temperature vary spatially /NO/


NUMERICAL_SETTINGS
!Degree of implicity /1.0000000/
!Relative error tolerance /0.50000000E-01/
!Absolute error tolerance /0.10000000E-07/
!Jacobian perturbation parameter /0.10000000E-02/
!Clean Jacobian /NO/
!Minimum mole fraction /0.10000000E-14/
!Minimum mobility /0.10000000E-39/
!Maximum temperature /6000.0000/
!Maximum timestep change factor /2.0000000/
!Enforce a maximum \Delta T per time-step /NO/
!Temporally piece-wise locally constant coefficients /YES/
y!Enable linearized time-steps /NO/
y!Assume a maximum width of interfaces /NO/
1.0e-6!Enter width of interface (nr of grid points, allowed values 3,5,7,..) /5/
!Set non-active phase-field variables away from interfaces to zero /YES/


INITIAL_SINGLE_PHASE
y!Set whole domain to consist of a single phase /NO/
2!Enter corresponding phase-field variable index

INITIAL_PRECIPITATE
1!Enter corresponding phase-field variable index
1!Precipitate option
0.64e-5!r0 spatial dimension #1
0.64e-5!r0 spatial dimension #2
2e-6!Enter precipitate radius
0.79!Element A
0.21!Element B
10!Relative grid fineness around precipitate (only relevant with adaptive meshing) /0/


ADAPTIVE_MESHING
y!Use adaptive meshing /NO/
!Initial allocation factor (>=1) /1.0000000/
!Adapt mesh based on composition gradients /YES/
!Max composition Delta for A /0.10000000/
!Min composition Delta for A /0.10000000E-01/
!Max composition Delta for B /0.10000000/
!Min composition Delta for B /0.10000000E-01/
!Adapt mesh based on phase-field variable gradients /YES/
0.01!Max phase-field variable Delta /0.10000000/
0.001!Min phase-field variable Delta /0.10000000E-01/
n!Adapt mesh based on temperature gradients /YES/
1!Integer aspect ratio, 1st spatial dimension
1!Integer aspect ratio, 2nd spatial dimension
!Input integer (1) or real valued (2) parameters /1/
7!Initial depth (>0)
7!Minimum depth (>0)
7!Maximum depth (>0)
!Base interval to adapt mesh on cumulative change /NO/


STRESS_EFFECTS
y!Include effect from stresses /NO/
!*For linear thermal expansion coefficient, choose option:
!*1.Enter manually
!Enter option /1/
0.0!Expression for linear thermal expansion coefficient
!*For Youngs modulus, choose option:
!*1.Enter manually
!Enter option /1/
uc_youngs!Expression for Youngs modulus
!*For Poissons ratio, choose option:
!*1.Enter manually
!Enter option /1/
uc_poissons!Expression for Poissons ratio
!*For lattice parameter, choose option:
!*1.Enter manually
!Enter option /1/
uc_a_alpha*p01**3*(10-15*p01+6*p01**2)+uc_a_beta*p02**3*(10-15*p02+6*p02**2)!Expression for lattice parameter
!*Enter expression for linear thermal strain
!*For example: ALPHA_THERM*(T-T0_STRESS)
0.0!Expression for linear thermal strain /ALPHA_THERM*(T-T0_STRESS)/
!*Enter expression for lattice parameter induced strain
!*For example:
!*(A_STRESS-A0_STRESS)/A0_STRESS
!Expression for lattice parameter induced linear strain /(A_STRESS-A0_STRESS)/A0_STRESS/
!*Give boundary conditions.
!*Note: vertex at x=y=0 is always prescribed to have zero displacement
!*In the x-direction, select option:
!*1.Unclamped (zero force at boundary)
!*2.Clamped (zero displacement at boundary)
!*3.Applied tension at x=xmax, zero displacement at x=xmin
2!Select boundary condition option (x) /1/
!*In the y-direction, select option:
!*1.Unclamped (zero force at boundary)
!*2.Clamped (zero displacement at boundary)
!*3.Applied tension at y=ymax, zero displacement at y=ymin
2!Select boundary condition option (y) /1/
!*Enter option for 2D stress state:
!*1. Plane stress
!*2. Plane strain
!2D stress state /1/
!Evaluate stress state every time-step /YES/




SOLVER
!*Solvers:
!*1.Jacobian solver
!*2.Implicit solver
2!Choose solver /1/

