SCPSandwich

serial

	Two-dimensional
	Serial
	Not validated
	./bin/alamo-2d-g++ tests/SCPSandwich/input stop_time="0.01"

serial-coverage

	Two-dimensional
	Serial
	Not validated
	./bin/alamo-2d-g++ tests/SCPSandwich/input stop_time="0.0001"

parallel

	Two-dimensional
	Parallel (4 procs)
	Validated using check script
	83.9s (waldorf)
	mpiexec -np 4 ./bin/alamo-2d-g++ tests/SCPSandwich/input

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

#@
#@  [serial]
#@  dim = 2
#@  check = false
#@  args = stop_time=0.01
#@  
#@  [serial-coverage]
#@  dim = 2
#@  check = false
#@  args = stop_time=0.0001
#@  coverage = true
#@
#@  [parallel]
#@  dim = 2
#@  nprocs = 4
#@  check = true
#@  benchmark-waldorf = 83.9
#@  

alamo.program = flame

# AMR parameters
plot_file = tests/SCPSandwich/output
amr.plot_dt = 0.01
amr.max_level = 7 
amr.n_cell = 16 2 2
amr.blocking_factor = 2
amr.regrid_int = 100
amr.base_regrid_int = 20
amr.grid_eff = 0.7

# Geometry
geometry.prob_lo = 0.0 -0.25 -0.1 # [mm]
geometry.prob_hi = 4.0 0.25 0.1 # [mm]
geometry.is_periodic = 0 1 1

# Timestep and duration
timestep = 0.00005 # [s]
stop_time = 0.2 # [s]

# Phase field params
pf.eps = 0.0005 # [mm]
pf.lambda = 0.001
pf.kappa = 1.0
pf.gamma = 0.02726
pf.w1 = 1.0
pf.w12 = 2.0
pf.w0 = 0.0

# AP / HTPB regression rate params 
pressure.r_ap = 1.222 # [ mm / s ]
pressure.n_ap = 1.042 # [ 1 ]
pressure.r_htpb = 0.1 # [ mm / s]
pressure.r_comb = 10.0 # [ mm / s]
pressure.n_htpb = 0.0 # [ 1 ]
pressure.n_comb = 0.0 # [ 1 ]
# Pressure used in regression rate calc
pressure.P = 0.2

# Phase field IC
eta.ic.type = constant
eta.ic.constant.value = 1

# Species field params
phi.ic.type = laminate
phi.ic.laminate.number_of_inclusions = 1
phi.ic.laminate.center  = 0.0 0.0 0.0
phi.ic.laminate.thickness = 0.1
phi.ic.laminate.orientation = 0 1 
phi.ic.laminate.eps = 0.015
phi.ic.laminate.singlefab = 1
phi.ic.laminate.invert = 1

# Initial conditions on eta
pf.eta.ic.type = constant
pf.eta.ic.constant.value = 1.0

# Boundary conditions on eta
pf.eta.bc.type = constant
pf.eta.bc.constant.type.xlo = dirichlet
pf.eta.bc.constant.type.xhi = dirichlet
pf.eta.bc.constant.type.ylo = neumann
pf.eta.bc.constant.type.yhi = neumann
pf.eta.bc.constant.val.xlo = 0.0
pf.eta.bc.constant.val.xhi = 1.0
pf.eta.bc.constant.val.ylo = 0.0
pf.eta.bc.constant.val.yhi = 0.0

pf.eta.ic.type = constant
pf.eta.ic.constant.value = 1.0

thermal.on = 0 
elastic.type = disable
#amr.node.any = 0