Inclusion

This demonstrates the evolution of an initially spherical inclusion, with anisotropic boundary energy, and with a volume constraint enforced using a Lagrange multiplier. See Integrator::PhaseFieldMicrostructure for additional information.

2d-serial

	Two-dimensional
	Serial
	Validated using check script
	29.49s (beaker) 20.54s (statler) 15.39s (waldorf)
	./bin/alamo-2d-g++ tests/Inclusion/input

2d-serial-coverage

	Two-dimensional
	Serial
	Not validated
	./bin/alamo-2d-g++ tests/Inclusion/input stop_time="0.02" anisotropy.tstart="0"

Input file (../../tests/Inclusion/input)

#@
#@  [2d-serial]
#@  dim = 2
#@  check = true
#@  benchmark-beaker = 29.49
#@  benchmark-statler = 20.54
#@  benchmark-waldorf = 15.39
#@  
#@  [2d-serial-coverage]
#@  dim = 2
#@  coverage=true
#@  check=false
#@  args = stop_time=0.02
#@  args = anisotropy.tstart=0
#@

alamo.program = microstructure

timestep = 0.01
stop_time = 2.0

plot_file = tests/Inclusion/output

amr.plot_dt = 0.1
amr.max_level = 2
amr.n_cell = 64 64 64
amr.blocking_factor = 8
amr.regrid_int = 10
amr.grid_eff = 1.0
amr.max_grid_size = 8

amr.thermo.plot_dt = 0.1

ic.type=sphere
ic.sphere.center = 4 4 4
ic.sphere.radius = 1

geometry.prob_lo = 0 0 0
geometry.prob_hi = 8 8 8 
geometry.is_periodic= 0 0 0

bc.eta.type = constant
bc.eta.constant.type.xlo = neumann
bc.eta.constant.type.xhi = neumann
bc.eta.constant.type.ylo = neumann
bc.eta.constant.type.yhi = neumann

pf.number_of_grains = 2
pf.M = 1.0 
pf.mu = 10.0
pf.gamma = 1.0
pf.l_gb=0.1
pf.sigma0=0.075

anisotropy.on=1
anisotropy.timestep=0.0005
anisotropy.tstart= 1.
anisotropy.type=sin
anisotropy.sin.sigma0=0.075
anisotropy.sin.sigma1=0.07
anisotropy.sin.theta0=45
anisotropy.beta= 0.00001

lagrange.on = 1
lagrange.lambda = 2.0
lagrange.tstart = 1.0
lagrange.vol0 = 4.0