FlowNoDiffuseInterface

central-pressure-wave

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

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

#@ [central-pressure-wave]
#@ exe = hydro
#@ check-file = reference/central-pressure-wave.dat
#@ dim    = 2
#@ nprocs = 1

alamo.program = hydro

### OUTPUT ###

plot_file = tests/FlowNoDiffuseInterface/output

### MESHING ###

amr.plot_dt = 0.1
#mr.plot_int = 1
amr.max_level = 0
amr.max_grid_size = 500000
amr.blocking_factor = 1
amr.grid_eff = 0.8
amr.n_cell = 64 8

### DIMENSIONS ###

geometry.prob_lo = 0.0 -0.5 0.0 # [ m ] 
geometry.prob_hi = 8.0  0.5 0.0 # [ m ]

geometry.is_periodic = 0 0 0

### TIME STEP ###

timestep  = 1.0e-2
stop_time = 2.5

### ETA IC ###

Ldot_0 = 0.0

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

deltapInterface.ic.type = expression
deltapInterface.ic.expression.region0 = "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.region0        = "1.0"

### HYDRO PARAM ###

gamma = 1.4
cfl   = 0.6
mu    = 0.0
pref = 0.0

rho_solid = 1.0 #kg/m^3
rho_fluid = 1.0 #kg/m^3

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

pressure.ic.type = expression
pressure.ic.expression.constant.gamma = 1.4
pressure.ic.expression.region0 = "(0.5*tanh((x-4.0)/0.5)+0.5)*0.1 + 1.0"

density.ic.type = constant
density.ic.constant.value=1.0


solid.velocity.ic.type=constant
solid.velocity.constant.value=1.0
solid.momentum.ic.constant.value=1.0
solid.energy.ic.constant.value=1.0
solid.density.ic.constant.value=1.0
m0.ic.constant.value=1.0 # todo - remove or rename
u0.ic.constant.value=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 = neumann
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 = 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.xhi = neumann neumann
momentum.bc.constant.type.xlo = neumann neumann
momentum.bc.constant.type.ylo = neumann neumann
momentum.bc.constant.type.yhi = neumann neumann
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 = 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

#TODO
v_solid = 0.0