FlowDrivenCavity
This tests the hydrodynamics solver for the classic driven cavity problem.
2d-serial
Two-dimensional |
|
Serial |
|
Validated using check script |
|
./bin/hydro-2d-g++ tests/FlowDrivenCavity/input
|
Input file (../../tests/FlowDrivenCavity/input)
#@
#@ [2d-serial]
#@ exe=hydro
#@ dim=2
#@
alamo.program = hydro
### OUTPUT ###
plot_file = tests/FlowDrivenCavity/output
### MESHING ###
amr.plot_dt = 0.1
amr.max_grid_size = 500000
amr.blocking_factor = 2
amr.regrid_int = 100
amr.grid_eff = 0.8
amr.max_level = 2
timestep = 0.5e-3
amr.n_cell = 32 32
### DIMENSIONS ###
geometry.prob_lo = -12.0 -12.0 0.0 # [ m ]
geometry.prob_hi = 12.0 12.0 0.0 # [ m ]
geometry.is_periodic = 0 0 0
### TIME STEP ###
stop_time = 1.0
q.ic.type = expression
q.ic.expression.region0 = "0.0"
q.ic.expression.region1 = "0.0"
eta.ic.type = expression
eta.ic.expression.region0 = "1.0"
### HYDRO PARAM ###
gamma = 1.4
cfl = 0.4
mu = 20.0
density.ic.type = expression
density.ic.expression.region0 = "1.2"
solid.density.ic.type = expression
solid.density.ic.expression.region0 = "1.2"
velocity.ic.type = expression
velocity.ic.expression.region0 = "0.0"
velocity.ic.expression.region1 = "0.0"
pressure.ic.type = expression
pressure.ic.expression.region0 = "103320.0"
pref = 103320.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
### HYDRO BC ###
density.bc.constant.type.xhi = neumann
density.bc.constant.type.xlo = neumann
density.bc.constant.type.ylo = neumann
density.bc.constant.type.yhi = neumann
energy.bc.constant.type.xhi = neumann
energy.bc.constant.type.xlo = neumann
energy.bc.constant.type.ylo = neumann
energy.bc.constant.type.yhi = neumann
momentum.bc.constant.type.xhi = dirichlet dirichlet
momentum.bc.constant.type.xlo = dirichlet dirichlet
momentum.bc.constant.type.ylo = dirichlet dirichlet
momentum.bc.constant.type.yhi = dirichlet dirichlet
momentum.bc.constant.val.xhi = 0.0 0.0
momentum.bc.constant.val.xlo = 0.0 0.0
momentum.bc.constant.val.ylo = 0.0 0.0
momentum.bc.constant.val.yhi = 500.0 0.0
### HYDRO REFINEMENT CRITERIA ###
eta_refinement_criterion = 1.1
omega_refinement_criterion = 10.55
#TODO
solid.momentum.ic.constant.value = 0.0
solid.energy.ic.constant.value = 0.0
m0.ic.constant.value = 1.0
u0.ic.constant.value = 0.0