FlowDrivenCavity

This tests the hydrodynamics solver for the classic driven cavity problem.

• Integrator::Hydro

2d-serial

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

2d-serial-rk4

	Two-dimensional
	Serial
	Validated using check script
	./bin/hydro-2d-g++ tests/FlowDrivenCavity/input integration.type="RungeKutta" integration.rk.type="4" timestep="1.0e-3"

2d-serial-ssprk3

	Two-dimensional
	Serial
	Validated using check script
	./bin/hydro-2d-g++ tests/FlowDrivenCavity/input integration.type="RungeKutta" integration.rk.type="3" timestep="1.0e-3"

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

#@
#@ [2d-serial]
#@ exe=hydro
#@ dim=2
#@
#@ [2d-serial-rk4]
#@ exe=hydro
#@ dim=2
#@ args = integration.type=RungeKutta
#@ args = integration.rk.type=4
#@ args = timestep=1.0e-3
#@
#@ [2d-serial-ssprk3]
#@ exe=hydro
#@ dim=2
#@ args = integration.type=RungeKutta
#@ args = integration.rk.type=3
#@ args = timestep=1.0e-3


### 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 = 1
timestep  = 5.0e-4
amr.n_cell = 64 64

### 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"

### 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 = 29.1_J/mol/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

### HYDRO PARAM ###

cfl   = 0.4

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

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

pressure.ic.type = expression
pressure.ic.expression.region0 = "1.0e5"

### 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 = dirichlet
density.bc.constant.val.yhi =  5.0

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.type = expression
momentum.bc.expression.constant.lidvelocity = 300.0
momentum.bc.expression.type.xhi = dirichlet dirichlet
momentum.bc.expression.type.xlo = dirichlet dirichlet
momentum.bc.expression.type.ylo = dirichlet dirichlet
momentum.bc.expression.type.yhi = dirichlet dirichlet
momentum.bc.expression.val.xhi = "0.0" "0.0"
momentum.bc.expression.val.xlo = "0.0" "0.0"
momentum.bc.expression.val.ylo = "0.0" "0.0"
momentum.bc.expression.val.yhi = "lidvelocity" "0.0"

### HYDRO REFINEMENT CRITERIA ###

eta_refinement_criterion = 1.1
omega_refinement_criterion = 10.55

#TODO


solid.density.ic.constant.value = 1.2
solid.momentum.ic.constant.value = 0.0
solid.energy.ic.constant.value = 0.0
m0.ic.constant.value = 0.0
u0.ic.constant.value = 0.0