!* 2D simulation of grain growth using the obsolescent GRAIN_GROWTH_MODEL


!*YAPFI macro
!*-----------
SET_NUMBER_OF_THREADS
4!Number of threads /1/

PRE_PROCESSOR

RANDOM_NUMBER_SETTINGS
y!Use a single seed for the pseudorandom number generator /NO/
-1!Enter integer seed (1 digit) (<0 for internal vector) /1/

VORONOI_STRUCTURE
32!Nr of Voronoi seeds (>1)
2!Dimensionality (1-3)
64!Number of grid points
64!Number of grid points
!Domain width [a.u.] /1.0000000/
!Domain width [a.u.] /1.0000000/
1!Random coordinates of seeds (=1), or manually entered (=2)
y!Periodic boundary
y!Periodic boundary
n!Varying growth rate of seeds
0.0!Minimum radial distance between seeds (>=0)
!*Output format options:
!*1.Suitable for GRAIN_GROWTH_MODEL
!*2.Suitable for PHI_SWITCHING
1!Select output option
!Output file name (leave blank to keep structure in memory)

RETURN



MANDATORY_DATA
analytical!GES file / ANALYTICAL [<analytic init file>]
1!Number of elements (>0) /2/
a!Element name
!Number of phases (>0) /1/
alpha!Phase name
!(Initial) Temperature [K] /1000.0000/
!Pressure [Pa] /100000.00/
!Element A substitutional /YES/
!Initially homogeneous composition /YES/
1.0!Initial mole fraction of element A
!*Enter state of phases, 1:Entered, 2: Suspended
!State of phase (ALPHA) /1/
2!Dimensionality /1/
64!Nr of grid points in dimension #1 /10/
64!Nr of grid points in dimension #2 /10/
64e-7!Domain width [m] in dimension #1 /1.0000000/
64e-7!Domain width [m] in dimension #2 /1.0000000/
y!Periodic boundary in dimension #1 /NO/
y!Periodic boundary in dimension #2 /NO/
!*Grid type. 1:Linear, 2:Geometric, 3: Double geometric, 4: Read from file
!Grid type in dimension #1 /1/
!Grid type in dimension #2 /1/
!*Info: Initial time is always zero
5!Simulation time [s] /1.0000000/
!Maximum timestep [s] /0.50000000/
!Minimum timestep [s] /0.10000000E-09/
!Initial timestep [s] /0.10000000E-09/
!Save interval /1/
0.1!Minimum elapsed time between saves [f(time)], [s] /1.0/

GRAIN_GROWTH_MODEL
!Enable the grain growth model /YES/
32!Initial nr of grains
!File containing corresponding initial phase-field variable values
!*Each grain should be associated with one phase
!*Phase indices:
!*1: ALPHA
1!Interactive phase association input (=1) or read from file (=2)
!*Simulation is single phase. All phi will be associated with same phase.
4!Initial nr of spare phase-field variables (>=0)
7!Give the width of the grain boundaries in finite volumes (1,3,5,..)
y!Use hard coded g2 expression
!*Info: The (unitless) height of the barrier at a boundary is
!*unity sans the energy prefactor
468.75!Grain/phase boundary energy prefactor (>0)

PHASE_FIELD_DATA
!Use phase-field variables /YES/
!*Info: The grain growth model has set the nr of phase-field variables to 12
!*Enter g_1
!*Shorthand notations:
!*gs11=\sum_\alpha G_{m\alpha}(\Phi_\alpha^3(10-15\Phi_\alpha+6\phi_\alpha^2))/
!*\sum_\alpha\Phi_alpha^3(10-15\Phi_\alpha+6\Phi_alpha^2)
!*where \Phi_alpha=\sum_{i\in\alpha}\phi_i
!*gs12=\sum_iG_{mi}\phi_i
!*variable GXX correspond to Gm of phase
0.0!Gm expression for g_1
!*Enter g_2
!*Shorthand notations:
!*gs21=8\sum_i\phi_i^2(1-\phi_i)^2
!*gs22=0.626\sum_i(\frac{\tanh(2)^2-1}{\tanh(2)^2}\log[1-\tanh(2)^2(2\phi_i-1)^2]-(2\phi_i-1)^2)
!*gs23=10.67(0.25-0.5\sum\phi_i^2+0.25(\sum\phi_i^2)^2+0.5\sum_i\sum_{j>i}\phi_i^2\phi_j^2
0.0!Gm expression for g_2
!*Enter epsilon_ij expressions (<=0 is enforced)
!*Enter REPEAT before the first expression to repeat it
repeat -0.2*10000*1e-7**2!E0102
0.05!Enter phase-field mobility expression

ANALYTICAL_EXPRESSIONS
0.0!Gm expression for ALPHA
0.0!Mobility expression for A in ALPHA
!Will temperature vary spatially /NO/

NUMERICAL_SETTINGS
0.0!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/
7e-7!Enter width of interface as real value (>0.0) or nr of grid points (3,5,7,..)
!Set non-active phase-field variables away from interfaces to zero /YES/


SAVE_SETTINGS
!Filename prefix /yapfisim/
!*Possible additional properties to be saved:
!*1.Gm
!*2.\varepsilon_{ij} (only for phase-field)
!*3.d\phi_i/dt (only for phase-field)
!*4.dg_{12}/d\phi_i (only for phase-field)
!*5.dg_\nabla/d\phi_i (only for phase-field)
!*6.\nabla\xi_i (only for phase-field)
!*7.active \phi (only for phase-field)
!*8.g_1 (only for phase-field)
!*9.g_2 (only for phase-field)
!*10.\sum\nabla\phi_i\sum\nabla\phi_j\varepsilon_{ij} (only for phase-field)
!*11.\sum\nabla c_i\sum\nabla c_j\kappa_{ij} (only if gradient energy contribution)
!*12.User defined functions
!*13.Save coordinates and variable values only in csv-like separate file
y!Add additional properties to be saved /NO/
10!Enter code for additional save property (-1*code to cancel, zero to exit)
0!Enter code for additional save property (-1*code to cancel, zero to exit)




