FlowChannel

serial

	Two-dimensional
	Serial
	Validated using check script
	./bin/hydro-2d-g++ tests/FlowChannel/input amr.plot_dt="0.5" stop_time="2.5"

fully-developed

	Two-dimensional
	Parallel (4 procs)
	Validated using check script
	mpiexec -np 4 ./bin/hydro-2d-g++ tests/FlowChannel/input amr.plot_dt="0.5" stop_time="10" velocity.ic.expression.region0=""0.0"" pressure.ic.expression.region0=""100 - 0.06*x"" density.bc.constant.type.xhi="neumann" density.bc.constant.val.xhi="0.0" momentum.bc.expression.type.xlo="neumann neumann" momentum.bc.expression.type.xhi="neumann neumann" momentum.bc.expression.val.xlo=""0" "0"" momentum.bc.expression.val.xhi=""0" "0"" energy.bc.expression.val.xhi=""49.91 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))""

parallel

	Two-dimensional
	Parallel (4 procs)
	Validated using check script
	mpiexec -np 4 ./bin/hydro-2d-g++ tests/FlowChannel/input amr.plot_dt="1.0" stop_time="10"

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

#@ [serial]
#@ exe = hydro
#@ dim = 2
#@ args = amr.plot_dt = 0.5
#@ args = stop_time=2.5
#@ check-file = reference/reference-2d.csv

#@ [fully-developed]
#@ exe = hydro
#@ dim = 2
#@ nprocs = 4
#@ args = amr.plot_dt = 0.5
#@ args = stop_time=10
#@ args = velocity.ic.expression.region0 = "0.0"
#@ args = pressure.ic.expression.region0 = "100 - 0.06*x"
#@ args = density.bc.constant.type.xhi=neumann
#@ args = density.bc.constant.val.xhi=0.0
#@ args = momentum.bc.expression.type.xlo=neumann neumann
#@ args = momentum.bc.expression.type.xhi=neumann neumann
#@ args = momentum.bc.expression.val.xlo="0" "0"
#@ args = momentum.bc.expression.val.xhi="0" "0"
#@ args = energy.bc.expression.val.xhi = "49.91 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
#@ check-file = reference/reference-2d-full.csv

#@ [parallel]
#@ exe = hydro
#@ dim = 2
#@ nprocs = 4
#@ args = amr.plot_dt = 1.0
#@ args = stop_time=10
#@ check-file = reference/reference-2d-long.csv

# Re = 0.01, Le ~= 1.26 (Chen, et al.)

### RIEMANN SOLVER ###
solver.type = hllc

### GAS MODEL ###
system.amount = kmol
gas.mw = 28.0_g/mol

gas.thermo.type = cpconstant
gas.thermo.cpconstant.cp_moles = 12471.69_J/kmol/K
gas.thermo.cpconstant.h0 = 0.0
gas.thermo.cpconstant.s0 = 0.0
gas.thermo.cpconstant.Tref = 0.0_K

gas.transport.type = mixture_averaged
gas.transport.mixture_averaged.type = constant
gas.transport.mixture_averaged.mu = 20.0_Pa*s
gas.transport.mixture_averaged.k = 0.023_W/m/K

gas.eos.type = cpg

### OUTPUT ###

plot_file = ./output

### MESHING ###

amr.plot_dt = 1.0
#amr.plot_int = 10
amr.max_level = 0
amr.max_grid_size = 500000
amr.blocking_factor = 2
amr.regrid_int = 10
amr.grid_eff = 0.8
amr.n_cell = 120 40 0
#amr.abort_on_nan=0

### DIMENSIONS ###

geometry.prob_lo = 0.0 0.0 0.0 # [ m ] 
geometry.prob_hi = 3.0 3.0 0.0 # [ m ]

geometry.is_periodic = 0 0 0

### TIME STEP ###

timestep  = 2.5e-4
stop_time = 10.0

### HYDRO PARAM ###

#gamma = 3.0
cfl   = 0.4

small=1E-8
#cutoff=1E-5
#lagrange=1E4

### ETA IC ###

eta.ic.type = expression
eta.ic.expression.constant.eps   = 0.1
eta.ic.expression.region0        = "(0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"

m0.ic.type = expression
m0.ic.expression.region0 = "0.0"

u0.ic.type = expression
u0.ic.expression.region0 = "0.0"
u0.ic.expression.region1 = "0.0"

q.ic.type = expression
q.ic.expression.region0 = "0.0"
q.ic.expression.region1 = "0.0"

density.ic.type = expression
density.ic.expression.region0 = "100.0"

velocity.ic.type = expression
velocity.ic.expression.constant.eps = 0.1
velocity.ic.expression.region0 = "0.001 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
velocity.ic.expression.region1 = "0.0"

pressure.ic.type = expression
pressure.ic.expression.region0 = "100.0"

solid.density.ic.type = expression
solid.density.ic.expression.region0 = "100.0"

solid.momentum.ic.type = expression
solid.momentum.ic.expression.region0 = "0.0"
solid.momentum.ic.expression.region1 = "0.0"

solid.energy.ic.type = expression
solid.energy.ic.expression.region0 = "0.0"

### ETA BC  ###

pf.eta.bc.constant.type.xlo = neumann
pf.eta.bc.constant.type.xhi = neumann
pf.eta.bc.constant.type.ylo = neumann
pf.eta.bc.constant.type.yhi = neumann
pf.eta.bc.constant.type.zlo = neumann
pf.eta.bc.constant.type.zhi = neumann

pf.eta.bc.constant.val.xlo = 0.0
pf.eta.bc.constant.val.xhi = 0.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

### HYDRO BC ###

density.bc.constant.type.xhi = neumann
density.bc.constant.type.xlo = dirichlet
density.bc.constant.type.ylo = neumann
density.bc.constant.type.yhi = neumann
density.bc.constant.val.xhi = 0.0
density.bc.constant.val.xlo = 100.0
density.bc.constant.val.ylo = 0.0
density.bc.constant.val.yhi = 0.0

energy.bc.type = expression
energy.bc.expression.constant.eps = 0.1
energy.bc.expression.type.xhi = dirichlet
energy.bc.expression.type.xlo = dirichlet
energy.bc.expression.type.ylo = neumann
energy.bc.expression.type.yhi = neumann
energy.bc.expression.val.xhi = "49.91 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
energy.bc.expression.val.xlo = "50.0 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
energy.bc.expression.val.ylo = "0.0 "
energy.bc.expression.val.yhi = "0.0"

momentum.bc.type = expression
momentum.bc.expression.constant.eps = 0.1
momentum.bc.expression.type.xlo = dirichlet dirichlet
momentum.bc.expression.type.xhi = neumann neumann
momentum.bc.expression.type.ylo = neumann neumann
momentum.bc.expression.type.yhi = neumann neumann
momentum.bc.expression.val.xlo = "0.1 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))" "0.0"
momentum.bc.expression.val.xhi = "0.0" "0.0"
momentum.bc.expression.val.ylo = "0.0" "0.0"
momentum.bc.expression.val.yhi = "0.0" "0.0"

### HYDRO REFINEMENT CRITERIA ###

eta_refinement_criterion = 0.1
omega_refinement_criterion = 0.2
gradu_refinement_criterion = 0.2
p_refinement_criterion   = 0.5
rho_refinement_criterion   = 0.5

#TODO
#
#dynamictimestep.on = 1
#dynamictimestep.verbose = 1
#dynamictimestep.cfl = 0.0004
#dynamictimestep.max  = 1.0
#dynamictimestep.min  = 0.0