FlowCouette

eps0.20

	Two-dimensional
	Serial
	Validated using check script
	./bin/hydro-2d-g++ tests/FlowCouette/input

eps0.10

	Two-dimensional
	Serial
	Validated using check script
	./bin/hydro-2d-g++ tests/FlowCouette/input eta.ic.expression.constant.eps="0.1"

eps0.05

	Two-dimensional
	Serial
	Validated using check script
	./bin/hydro-2d-g++ tests/FlowCouette/input eta.ic.expression.constant.eps="0.05" amr.max_level="2"

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

#@
#@ [eps0.20]
#@ exe=hydro
#@ dim=2
#@ check-file = reference/eps0.20.dat
#@
#@ [eps0.10]
#@ exe=hydro
#@ dim=2
#@ args = eta.ic.expression.constant.eps=0.1
#@ check-file = reference/eps0.10.dat
#@
#@ [eps0.05]
#@ exe=hydro
#@ dim=2
#@ args = eta.ic.expression.constant.eps=0.05
#@ args = amr.max_level=2
#@ check-file = reference/eps0.05.dat
#@


alamo.program = hydro

### OUTPUT ###

plot_file = ./tests/FlowCouette/output

### MESHING ###

amr.plot_dt = 0.5
amr.max_level = 1
amr.max_grid_size = 500000
amr.blocking_factor = 2
#amr.regrid_int = 100
amr.grid_eff = 0.8
amr.n_cell = 4 64 0

### DIMENSIONS ###
geometry.prob_lo = 0.0 0.0 0.0 # [ m ] 
geometry.prob_hi = 0.1875 3.0 0.0 # [ m ]
geometry.is_periodic = 1 0 0

### TIME STEP ###

timestep  = 1.0e-3
stop_time = 30.0

### HYDRO PARAM ###

gamma = 3.0
cfl   = 0.4
mu    = 10.0
pref = 100.0

### ETA IC ###

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

m0.ic.type = constant
m0.ic.constant.value = "0.0"

u0.ic.type = expression
u0.ic.expression.region0 = "0.001*(y>1.5)"
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.region0 = "0.0"
velocity.ic.expression.region1 = "0.0"

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

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

solid.momentum.ic.type = expression
solid.momentum.ic.expression.region0 = "0.1*(y>1.5)" #"0.5*(tanh((y-1.5)/.05) + 1)" # 0 at bottom plate and 1 at top plate
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 = periodic
pf.eta.bc.constant.type.xhi = periodic
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 = periodic
density.bc.constant.type.xlo = periodic
density.bc.constant.type.ylo = neumann
density.bc.constant.type.yhi = neumann
density.bc.constant.val.xhi = 0.0
density.bc.constant.val.xlo = 0.0
density.bc.constant.val.ylo = 0.0
density.bc.constant.val.yhi = 0.0

energy.bc.constant.type.xhi = periodic
energy.bc.constant.type.xlo = periodic
energy.bc.constant.type.ylo = neumann
energy.bc.constant.type.yhi = neumann
energy.bc.constant.val.xhi = 0.0
energy.bc.constant.val.xlo = 0.0
energy.bc.constant.val.ylo = 0.0 
energy.bc.constant.val.yhi = 0.0

momentum.bc.constant.type.xlo = periodic periodic
momentum.bc.constant.type.xhi = periodic periodic
momentum.bc.constant.type.ylo = neumann neumann
momentum.bc.constant.type.yhi = neumann neumann
momentum.bc.constant.val.xlo = 0.0 0.0
momentum.bc.constant.val.xhi = 0.0 0.0
momentum.bc.constant.val.ylo = 0.0 0.0
momentum.bc.constant.val.yhi = 0.0 0.0

### HYDRO REFINEMENT CRITERIA ###

r_refinement_criterion   = 0.1
e_refinement_criterion   = 0.1
m_refinement_criterion   = 0.1
eta_refinement_criterion = 0.1
omega_refinement_criterion = 0.1
gradu_refinement_criterion   = 0.1
p_refinement_criterion   = 0.1
rho_refinement_criterion   = 0.1

#TODO