FlowVortexShed

serial

	Two-dimensional
	Serial
	Validated using check script
	./bin/hydro-2d-g++ tests/FlowVortexShed/input stop_time="25"

parallel

	Two-dimensional
	Parallel (4 procs)
	Validated using check script
	mpiexec -np 4 ./bin/hydro-2d-g++ tests/FlowVortexShed/input stop_time="100"

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

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

#@ [parallel]
#@ exe=hydro
#@ dim=2
#@ nprocs = 4
#@ args = stop_time = 100
#@ check-file = reference/parallel.csv

alamo.program = hydro

### OUTPUT ###

plot_file = ./output.VortexShed

### MESHING ###

amr.plot_dt = 5.0
#amr.plot_int = 1000
amr.max_grid_size = 500000
amr.blocking_factor = 4
amr.regrid_int = 100
amr.grid_eff = 0.8


amr.max_level = 1
timestep  = 2e-2
amr.n_cell = 60 32


### DIMENSIONS ###

geometry.prob_lo = -6.0 -5.0 0.0 # [ m ] 
geometry.prob_hi =  14.0 5.0 0.0 # [ m ]

geometry.is_periodic = 0 0 0

### TIME STEP ###

stop_time = 100.0

### ETA IC ###

m0.ic.type = constant
m0.ic.constant.value = 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"

eta.ic.type = expression
eta.ic.expression.constant.eps   = 0.2
eta.ic.expression.constant.vel_x = 0.0 
eta.ic.expression.constant.vel_y = 0.0 
eta.ic.expression.constant.x0    = 0.0 
eta.ic.expression.constant.y0    = 0.0
eta.ic.expression.constant.R1    = 1.5
eta.ic.expression.constant.R2    = 1.6
eta.ic.expression.constant.ycut  = -1.5
eta.ic.expression.region0        = "(y>=ycut)*(0.5*(1 + tanh((sqrt((x-x0)**2 + (y-y0)**2) - R1)/eps))) + (y<ycut)*(0.5*(1 + tanh((sqrt((x-x0)**2 + (y-y0)**2) - R2)/eps)))"

### HYDRO PARAM ###

gamma = 1.4
cfl   = 0.01
mu    = 1e-1
pref = 1000.0

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

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

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

### SOLID ###

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.constant.type.xhi = neumann
energy.bc.constant.type.xlo = neumann
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 = dirichlet neumann
momentum.bc.constant.type.xhi = neumann neumann
momentum.bc.constant.type.ylo = neumann dirichlet
momentum.bc.constant.type.yhi = neumann dirichlet
momentum.bc.constant.val.xlo = 100.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.2
gradu_refinement_criterion = 1.0

#TODO

amr.abort_on_nan = 0

dynamictimestep.on = 1
dynamictimestep.verbose = 1
dynamictimestep.cfl = 0.01
dynamictimestep.max  = 1.0
dynamictimestep.min  = 0.0


small=1E-8
cutoff=0.01
lagrange=2000.0