!* Simple example of a 1D diffusion controlled phase transformation


SET_NUMBER_OF_THREADS
4!Number of threads /1/

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/
!Dimensionality /1/
128!Number of grid points /10/
1.28e-5!Domain width [m] /1.0000000/
!Periodic boundary /NO/
!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
(0.856/8)*gs21!Gm expression for g_2
-2.625e-10!E0102
0.05!Enter phase-field mobility expression

ANALYTICAL_EXPRESSIONS
10E3*XA+20E3*XB+R*T*(XA*LOG(XA)+XB*LOG(XB))-10E3*XA*XB!Gm for ALPHA
20E3*XA+10E3*XB+R*T*(XA*LOG(XA)+XB*LOG(XB))-20E3*XA*XB!Gm for BETA
1e-16!Mobility for A in ALPHA
1e-16!Mobility for B in ALPHA
1e-16!Mobility for A in BETA
1e-16!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/
8e-7!Enter width of interface [m]
!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
2!Precipitate option
1e-6!x-axis cut-off
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/
!Input integer (1) or real valued (2) parameters /1/
7!Initial depth (>0)
5!Minimum depth (>0)
7!Maximum depth (>0)
!Base interval to adapt mesh on cumulative change /NO/




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

MATRIX_SOLVER
!*Matrix solvers
!*1.QMR
!*2.BiConjugate Gradient Stabilized
!*3.PARDISO
3!Matrix solver /1/


SET_SCREEN_OUTPUT
3!Screen output. 1=Normal, 2=Minimal, 3=Less than minimal, 4=Almost nothing /2/

