FlowCouette
eps0.20
Two-dimensional |
|
Serial |
|
Validated using check script |
|
./bin/hydro-2d-g++ tests/FlowCouette/input eta.ic.expression.constant.eps="0.2"
|
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 |
|
Parallel (4 procs) |
|
Validated using check script |
|
mpiexec -np 4 ./bin/hydro-2d-g++ tests/FlowCouette/input eta.ic.expression.constant.eps="0.05" timestep="5e-4"
|
Input file (../../tests/FlowCouette/input)
#@
#@ [eps0.20]
#@ exe=hydro
#@ dim=2
#@ args = eta.ic.expression.constant.eps=0.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
#@ nprocs=4
#@ args = eta.ic.expression.constant.eps=0.05
#@ args = timestep=5e-4
#@ check-file = reference/eps0.05.dat
#@
##### This should work if trying to debug a high accuracy sim,
##### but it is too expensive to keep as a test case
#####
##### #@ [eps0.03]
##### #@ exe=hydro
##### #@ dim=2
##### #@ nprocs=4
##### #@ args = eta.ic.expression.constant.eps=0.03
##### #@ args = amr.n_cell=8 124 0
##### #@ args = timestep=2e-4
##### #@ check = false
### 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 = 10.0_Pa*s
gas.transport.mixture_averaged.k = 0.023_W/m/K
gas.eos.type = cpg
### OUTPUT ###
plot_file = output.couette
### MESHING ###
amr.plot_dt = 0.5
amr.max_level = 1
amr.max_grid_size = 500000
amr.blocking_factor = 2
amr.regrid_int = 10
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 = 10.0
### HYDRO PARAM ###
cfl = 0.4
small = 1e-16
### ETA IC ###
eta.ic.type = expression
eta.ic.expression.constant.eps = 0.10
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.000*(y>1.5)"
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.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 ###
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