SCPThermalSandwich

2d-serial

	Two-dimensional
	Serial
	Validated using check script
	./bin/alamo-2d-g++ tests/SCPThermalSandwich/input stop_time="0.001_s" amr.plot_dt="0.001_s"

2d-serial-coverage

	Two-dimensional
	Serial
	Not validated
	./bin/alamo-2d-g++ tests/SCPThermalSandwich/input stop_time="1.0e-4"

2d-parallel

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

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

#@ [2d-serial]
#@ dim = 2
#@ check = true
#@ check-file = reference/serial.dat
#@ args = stop_time=0.001_s
#@ args = amr.plot_dt=0.001_s
#@ 
#@ [2d-serial-coverage]
#@ dim = 2
#@ check = false
#@ args = stop_time=1.0e-4
#@ coverage = true
#@ 
#@ [2d-parallel]
#@ dim = 2
#@ nprocs = 4
#@ check = true
#@ check-file = reference/reference.csv

alamo.program = flame

plot_file = tests/SCPThermalSandwich/output

system.length = m
system.time = s

# AMR and output parameters
amr.plot_dt = 0.001_s
amr.max_level = 7
amr.max_grid_size = 50
amr.blocking_factor = 2
amr.base_regrid_int = 10
amr.grid_eff = 0.8
amr.refinement_criterion = 0.1
amr.refinement_criterion_temp = 10.0_K

# Geometry and base grid size
amr.n_cell = 16 2 2
geometry.prob_lo =    0.0_um -250._um -250._um
geometry.prob_hi = 4000.0_um  250._um  250._um
geometry.is_periodic = 0 1 1

## Timestep
timestep = 1.0e-5_s # 0.00005 # [s]
stop_time = 0.03_s

small = 1.0e-8
pf.eps = 1.0_um #0.0005 # [m]
pf.lambda = 0.001_J/m^2

## Eta initial condition
pf.eta.ic.type = constant 
pf.eta.ic.constant.value = 1.0 
#pf.eta.ic.type = expression
#pf.eta.ic.expression.constant.eps = 1.0_um
#pf.eta.ic.expression.constant.x0 = 40.0_um
#pf.eta.ic.expression.region0 = "0.5 + 0.5*erf((x-x0)/eps)"

## Phi field initial conditions
phi.ic.type = laminate
phi.ic.laminate.center = 0.0 0.0 0.0 
phi.ic.laminate.thickness = 100.0_um
phi.ic.laminate.orientation = 0 1 
phi.ic.laminate.singlefab = 1 
phi.ic.laminate.invert = 1
phi.ic.laminate.eps = 20._um

# Phase field parameters
pf.kappa      = 1.0_J/m^2
pf.w1         = 1.0_1
pf.w12        = 2.0_1
pf.w0         = 0.0_1

# Eta boundary conditions
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 = periodic
pf.eta.bc.constant.type.yhi = periodic
pf.eta.bc.constant.type.zlo = periodic
pf.eta.bc.constant.type.zhi = periodic
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.bc.constant.val.zlo = 0.0
pf.eta.bc.constant.val.zhi = 0.0

# Propellant model - requires thermal model
propellant.type = fullfeedback
# arrhenius rate law parameters
propellant.fullfeedback.m_ap      = 1.45e5_m/s
propellant.fullfeedback.m_htpb    = 1.4e1_m/s
propellant.fullfeedback.E_ap      = 11000.0_K
propellant.fullfeedback.E_htpb    = 7500.0_K
# thermal transport parameters
propellant.fullfeedback.rho_ap      =  1950.0_kg/m^3
propellant.fullfeedback.rho_htpb    =   920.0_kg/m^3
propellant.fullfeedback.k_ap        =  0.4186_W/m/K
propellant.fullfeedback.k_htpb      =  0.1300_W/m/K 
propellant.fullfeedback.cp_ap       = 1297.90_J/kg/K 
propellant.fullfeedback.cp_htpb     = 2418.29_J/kg/K
propellant.fullfeedback.mlocal_ap   =  1000.0_kg/m^3/s
propellant.fullfeedback.mlocal_comb =     0.0_kg/m^3/s
# Mass-to-heat flux fit parameters
propellant.fullfeedback.a1       =    1.114_1
propellant.fullfeedback.a2       =    0.460_1
propellant.fullfeedback.a3       =    2.797_1
propellant.fullfeedback.b1       =    0.323_1
propellant.fullfeedback.b2       =    0.420_1
propellant.fullfeedback.b3       =   0.3225_1
propellant.fullfeedback.c1       = -0.09906_1
propellant.fullfeedback.phi.zeta =     20.0_um
propellant.fullfeedback.Pref     =      1.0_MPa

# Thermal transport parameters
thermal.on = 1  
#thermal.Tfluid = 300.0_K
thermal.Tref = 300.0_K

# Temperature field initial condition
temp.ic.type = constant
temp.ic.constant.value = 300.0_K

# Temperature field boundary conditions
thermal.temp.bc.type = constant
thermal.temp.bc.constant.type.xlo = neumann
thermal.temp.bc.constant.type.xhi = neumann
thermal.temp.bc.constant.type.ylo = neumann
thermal.temp.bc.constant.type.yhi = neumann
thermal.temp.bc.constant.type.zlo = neumann
thermal.temp.bc.constant.type.zhi = neumann
thermal.temp.bc.constant.val.xlo = 0.0
thermal.temp.bc.constant.val.xhi = 0.0
thermal.temp.bc.constant.val.ylo = 0.0 
thermal.temp.bc.constant.val.yhi = 0.0
thermal.temp.bc.constant.val.zlo = 0.0 
thermal.temp.bc.constant.val.zhi = 0.0

# volumetric heat flux multiplier
thermal.hc = 1.0e7_W/m^2

# Laser initial condition
laser.ic.type = constant
laser.ic.constant.value = 1.0e6_W/m^2

# Chamber pressure (used for regression model)
chamber.pressure = 4.0_MPa  ## [MPa]

# Disable elasticity
elastic.type = disable

# Disable plotting all extra fields
plot_field = 1

# Disable NAN checking
amr.abort_on_nan = 0