Added sharp interface case for channel flow.

input.ChannelFlow_eps0.00

@@ -0,0 +1,138 @@
+alamo.program = hydro
+
+### OUTPUT ###
+
+plot_file = ./output..ChannelFlow_eps0.00
+
+### MESHING ###
+
+amr.plot_dt = 0.25
+amr.plot_int = 10
+amr.max_level = 0
+amr.max_grid_size = 500000
+amr.blocking_factor = 2
+amr.regrid_int = 100
+amr.grid_eff = 0.8
+amr.n_cell = 400 80 0
+
+### DIMENSIONS ###
+
+geometry.prob_lo = 0.0 0.5 0.0 # [ m ] 
+geometry.prob_hi = 10.0 2.5 0.0 # [ m ]
+
+
+geometry.is_periodic = 0 0 0
+
+### TIME STEP ###
+
+timestep  = 1.0e-3
+stop_time = 70.0
+
+### HYDRO PARAM ###
+
+gamma = 3.0
+cfl   = 0.0001
+mu    = 10.0
+pref = 100.0
+
+dynamictimestep.on = 1
+dynamictimestep.verbose = 1
+dynamictimestep.max  = 1E-3
+dynamictimestep.min  = 1E-8
+
+
+
+### ETA IC ###
+
+eta.ic.type = expression
+eta.ic.expression.constant.eps   = 0.05
+eta.ic.expression.region0        = "(0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
+
+rho_injected.ic.type = expression
+rho_injected.ic.expression.region0 = "100.0"
+
+mdot.ic.type = expression
+mdot.ic.expression.region0 = "0.0"
+mdot.ic.expression.region1 = "0.0"
+
+q.ic.type = expression
+q.ic.expression.region0 = "0.0"
+q.ic.expression.region1 = "0.0"
+
+density.ic.type = expression
+density.ic.expression.region0 = "100.0"
+
+velocity.ic.type = expression
+velocity.ic.expression.constant.eps = 0.03
+velocity.ic.expression.region0 = "0.01"
+velocity.ic.expression.region1 = "0.0"
+
+pressure.ic.type = expression
+pressure.ic.expression.region0 = "0.0" #"(0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
+
+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.type.xlo = neumann
+pf.eta.bc.type.xhi = neumann
+pf.eta.bc.type.ylo = neumann
+pf.eta.bc.type.yhi = neumann
+pf.eta.bc.type.zlo = neumann
+pf.eta.bc.type.zhi = neumann
+
+pf.eta.bc.val.xlo = 0.0
+pf.eta.bc.val.xhi = 0.0
+pf.eta.bc.val.ylo = 0.0
+pf.eta.bc.val.yhi = 0.0
+pf.eta.bc.val.zlo = 0.0
+pf.eta.bc.val.zhi = 0.0
+
+### HYDRO BC ###
+
+density.bc.type.xhi = neumann
+density.bc.type.xlo = neumann
+density.bc.type.ylo = neumann
+density.bc.type.yhi = neumann
+density.bc.val.xhi = 0.0
+density.bc.val.xlo = 0.0
+density.bc.val.ylo = 0.0
+density.bc.val.yhi = 0.0
+
+energy.bc.type.xhi = dirichlet
+energy.bc.type.xlo = neumann
+energy.bc.type.ylo = neumann
+energy.bc.type.yhi = neumann
+energy.bc.val.xhi = 0.0
+energy.bc.val.xlo = 0.0
+energy.bc.val.ylo = 0.0 
+energy.bc.val.yhi = 0.0
+
+momentum.bc.type.xlo = dirichlet dirichlet
+momentum.bc.type.xhi = neumann neumann
+momentum.bc.type.ylo = dirichlet dirichlet
+momentum.bc.type.yhi = dirichlet dirichlet
+momentum.bc.val.xlo = "1.0" "0.0"
+momentum.bc.val.xhi = "0.0" "0.0"
+momentum.bc.val.ylo = "0.0" "0.0"
+momentum.bc.val.yhi = "0.0" "0.0"
+
+### HYDRO REFINEMENT CRITERIA ###
+
+r_refinement_criterion   = 1000
+e_refinement_criterion   = 1000
+m_refinement_criterion   = 1000
+eta_refinement_criterion = 1000
+omega_refinement_criterion = 1000
+p_refinement_criterion   = 1000
+rho_refinement_criterion   = 1000
+
+#TODO

input.VortexShed

@@ -2,7 +2,7 @@ alamo.program = hydro
 
 ### OUTPUT ###
 
-plot_file = ./output
+plot_file = ./output.VortexShed
 
 ### MESHING ###
 
@@ -59,7 +59,7 @@ eta.ic.expression.region0        = "0.5*(1 + tanh((sqrt((x-x0)**2 + (y)**2) - 3.
 ### HYDRO PARAM ###
 
 gamma = 1.4
-cfl   = 0.1
+cfl   = 0.01
 mu    = 10
 #pref = 1000.0
 pref = 1000.0
@@ -135,7 +135,7 @@ momentum.bc.val.yhi = 0.0 0.0
 
 eta_refinement_criterion = 0.1
 omega_refinement_criterion = 0.01
-gradu_refinement_criterion = 0.05
+gradu_refinement_criterion = 0.1
 
 #TODO