FlowDendrite

serial

	Two-dimensional
	Serial
	Validated using check script
	./bin/sfi-2d-g++ tests/FlowDendrite/input stop_time="0.6"

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

#@ [serial]
#@ exe=sfi
#@ dim=2
#@ args = stop_time=0.6

# Simulation length
timestep = 1E-3
stop_time = 1.0

# AMR parameters
amr.plot_dt = 0.1
amr.max_level = 2
amr.n_cell = 48 32 32
amr.blocking_factor = 8
amr.regrid_int = 10
amr.grid_eff = 1.0

# Specify geometry and unrefined mesh
geometry.prob_lo = -4.0 0 0
geometry.prob_hi = 4.0 6.0 6.0
geometry.is_periodic= 0 0 0

alamo.program = dendrite
plot_file = ./output.dendrite


#amr.n_cell = 64 64
#geometry.prob_lo = -24.0 -24.0 0.0 # [ m ] 
#geometry.prob_hi =  24.0  24.0 0.0 # [ m ]
#geometry.is_periodic = 0 0 0
#stop_time = 1000.0


tstart = 000.1
invert = 1

dendrite.ic.phi.region0  = "(x-1.0)^2 + (y-1.0)^2 + z^2 < 0.1"
dendrite.ic.temp.constant.DeltaT = 0.25
dendrite.ic.temp.region0 = "-DeltaT*(x^2 + y^2 + z^2 > 0.1)"
dendrite.tau = 3E-4
dendrite.diffusion = 1.0
dendrite.eps = 1E-2
dendrite.alpha = 0.9
dendrite.gamma = 55.0
dendrite.delta = 0.5
dendrite.theta = 15.0
dendrite.bc.temp.constant.type.xhi = neumann
dendrite.bc.temp.constant.type.xlo = neumann
dendrite.bc.temp.constant.type.yhi = neumann
dendrite.bc.temp.constant.type.ylo = neumann
dendrite.bc.temp.constant.type.zhi = neumann
dendrite.bc.temp.constant.type.zlo = neumann
dendrite.bc.phi.constant.type.xhi = neumann
dendrite.bc.phi.constant.type.xlo = neumann
dendrite.bc.phi.constant.type.yhi = neumann
dendrite.bc.phi.constant.type.ylo = neumann
dendrite.bc.phi.constant.type.zhi = neumann
dendrite.bc.phi.constant.type.zlo = neumann




### ETA IC ###
#hydro.rhoInterface.ic.type = expression
#hydro.rhoInterface.ic.expression.region0 = "20.0"
#hydro.vInjected.ic.type = expression
#hydro.vInjected.ic.expression.region0 = "0.0"
#hydro.vInjected.ic.expression.region1 = "0.0"
#hydro.deltapInterface.ic.type = expression
#hydro.deltapInterface.ic.expression.region0 = "0.0"
### HYDRO PARAM ###
hydro.gamma = 1.4
hydro.mu    = 2.0
hydro.pref    = 10000
#hydro.rho_solid = 30.0 #kg/m^3
#hydro.rho_fluid = 20.0 #kg/m^3
#hydro.E_solid   = 5.0
hydro.density.ic.constant.value = 40.0
hydro.velocity.ic.type = expression
hydro.velocity.ic.expression.constant.eps   = 0.1
hydro.velocity.ic.expression.constant.x0    = -2.0
hydro.velocity.ic.expression.region0 = "0.0"
hydro.velocity.ic.expression.region1 = "0.0"
hydro.pressure.ic.type = expression
hydro.pressure.ic.expression.constant.gamma = 1.4
hydro.pressure.ic.expression.constant.E_solid = 5.0
hydro.pressure.ic.expression.region0 = "(gamma - 1.0) * E_solid"
### ETA BC  ###
hydro.pf.eta.bc.constant.type.xlo = neumann
hydro.pf.eta.bc.constant.type.xhi = neumann
hydro.pf.eta.bc.constant.type.ylo = neumann
hydro.pf.eta.bc.constant.type.yhi = neumann
hydro.pf.eta.bc.constant.type.zlo = neumann
hydro.pf.eta.bc.constant.type.zhi = neumann
hydro.pf.eta.bc.constant.val.xlo = 0.0
hydro.pf.eta.bc.constant.val.xhi = 0.0
hydro.pf.eta.bc.constant.val.ylo = 0.0
hydro.pf.eta.bc.constant.val.yhi = 0.0
hydro.pf.eta.bc.constant.val.zlo = 0.0
hydro.pf.eta.bc.constant.val.zhi = 0.0
### HYDRO BC ###
hydro.density.bc.constant.type.xhi = neumann
hydro.density.bc.constant.type.xlo = dirichlet
hydro.density.bc.constant.type.ylo = neumann
hydro.density.bc.constant.type.yhi = neumann
hydro.density.bc.constant.val.xhi = 0.0
hydro.density.bc.constant.val.xlo = 1.0
hydro.density.bc.constant.val.ylo = 0.0 
hydro.density.bc.constant.val.yhi = 0.0
hydro.energy.bc.constant.type.xhi = neumann
hydro.energy.bc.constant.type.xlo = dirichlet #neumann
hydro.energy.bc.constant.type.ylo = neumann
hydro.energy.bc.constant.type.yhi = neumann
hydro.energy.bc.constant.val.xhi = 0.0
hydro.energy.bc.constant.val.xlo = 1.0 #0.0
hydro.energy.bc.constant.val.ylo = 0.0 
hydro.energy.bc.constant.val.yhi = 0.0
hydro.momentum.bc.constant.type.xhi = neumann neumann
hydro.momentum.bc.constant.type.xlo = neumann neumann #dirichlet neumann
hydro.momentum.bc.constant.type.ylo = neumann dirichlet
hydro.momentum.bc.constant.type.yhi = neumann dirichlet
hydro.momentum.bc.constant.val.xhi = 0.0 0.0
hydro.momentum.bc.constant.val.xlo = 0.0  0.0 #10. 0.0
hydro.momentum.bc.constant.val.ylo = 0.0 0.0
hydro.momentum.bc.constant.val.yhi = 0.0 0.0
hydro.eta.ic.constant.value = 1.0
### HYDRO REFINEMENT CRITERIA ###
#hydro.r_refinement_criterion   = 0.1
#hydro.e_refinement_criterion   = 0.1
#hydro.m_refinement_criterion   = 0.1
hydro.eta_refinement_criterion = 0.1
hydro.omega_refinement_criterion = 10000000
hydro.gradu_refinement_criterion = 10000000

hydro.solid.momentum.ic.constant.value = 0.0
hydro.solid.density.ic.constant.value = 40.0
hydro.solid.energy.ic.constant.value = 0.0

hydro.m0.ic.constant.value = 0.0
hydro.u0.ic.constant.value = 0.0
hydro.q.ic.constant.value = 0.0

hydro.small=1E-8
hydro.cutoff=0.01
hydro.lagrange=200.0

#hydro.Lfactor=-10
#hydro.Pfactor=10



amr.abort_on_nan = 0

dynamictimestep.on = 1
dynamictimestep.verbose = 1
dynamictimestep.max  = 2E-4
dynamictimestep.min  = 1E-8
hydro.cfl   = 10.0
hydro.cfl_v   = 1.0
#hydro.small=1E-4