Corrections to AMR criterions. Channel flow and Couette flow input fi…

…les.

input.ChannelFlow_eps0.00

@@ -2,25 +2,24 @@ alamo.program = hydro
 
 ### OUTPUT ###
 
-plot_file = ./output..ChannelFlow_eps0.00
+plot_file = ./output.channelflow_eps0.00
 
 ### MESHING ###
 
-amr.plot_dt = 0.25
-amr.plot_int = 10
-amr.max_level = 0
+amr.plot_dt = 1.0
+#amr.plot_int = 10
+amr.max_level = 1
 amr.max_grid_size = 500000
 amr.blocking_factor = 2
-amr.regrid_int = 100
+#amr.regrid_int = 100
 amr.grid_eff = 0.8
-amr.n_cell = 400 80 0
+amr.n_cell = 200 60 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 ###
@@ -31,22 +30,15 @@ stop_time = 70.0
 ### HYDRO PARAM ###
 
 gamma = 3.0
-cfl   = 0.0001
+cfl   = 0.4
 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))"
+eta.ic.expression.region0        = "1.0" #"(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"
@@ -63,8 +55,8 @@ 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.constant.eps = 0.05
+velocity.ic.expression.region0 = "0.01" # * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
 velocity.ic.expression.region1 = "0.0"
 
 pressure.ic.type = expression
@@ -120,19 +112,20 @@ 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.xlo = "0.1" "0.0" #*(y>0.5)*(y<2.5)" "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
+r_refinement_criterion   = 0.1
+e_refinement_criterion   = 0.1
+m_refinement_criterion   = 0.1
+eta_refinement_criterion = 0.1
+omega_refinement_criterion = 0.1
+gradu_refinement_criterion = 0.1
+p_refinement_criterion   = 0.1
+rho_refinement_criterion   = 0.1
 
 #TODO

input.ChannelFlow_eps0.03

@@ -2,18 +2,18 @@ alamo.program = hydro
 
 ### OUTPUT ###
 
-plot_file = ./output
+plot_file = ./output.channelflow_eps0.03
 
 ### MESHING ###
 
 amr.plot_dt = 1.0
 #amr.plot_int = 10
-amr.max_level = 1
+amr.max_level = 3
 amr.max_grid_size = 500000
 amr.blocking_factor = 2
 #amr.regrid_int = 100
 amr.grid_eff = 0.8
-amr.n_cell = 800 240 0
+amr.n_cell = 200 60 0
 
 ### DIMENSIONS ###
 
@@ -24,8 +24,8 @@ geometry.is_periodic = 0 0 0
 
 ### TIME STEP ###
 
-timestep  = 1.0e-4
-stop_time = 1000.0
+timestep  = 5.0e-4
+stop_time = 70.0
 
 ### HYDRO PARAM ###
 
@@ -56,7 +56,7 @@ density.ic.expression.region0 = "100.0"
 
 velocity.ic.type = expression
 velocity.ic.expression.constant.eps = 0.03
-velocity.ic.expression.region0 = "0.01 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
+velocity.ic.expression.region0 = "0.001 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
 velocity.ic.expression.region1 = "0.0"
 
 pressure.ic.type = expression
@@ -99,23 +99,24 @@ density.bc.val.xlo = 0.0
 density.bc.val.ylo = 0.0
 density.bc.val.yhi = 0.0
 
-pressure.bc.type.xhi = dirichlet
-pressure.bc.type.xlo = neumann
-pressure.bc.type.ylo = neumann
-pressure.bc.type.yhi = neumann
-pressure.bc.val.xhi = 0.0
-pressure.bc.val.xlo = 0.0
-pressure.bc.val.ylo = 0.0 
-pressure.bc.val.yhi = 0.0
-
-velocity.bc.type.xlo = dirichlet dirichlet
-velocity.bc.type.xhi = neumann neumann
-velocity.bc.type.ylo = neumann neumann
-velocity.bc.type.yhi = neumann neumann
-velocity.bc.val.xlo = "0.1*(y>0.5)*(y<2.5)" "0.0"
-velocity.bc.val.xhi = "0.0" "0.0"
-velocity.bc.val.ylo = "0.0" "0.0"
-velocity.bc.val.yhi = "0.0" "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 = neumann neumann
+momentum.bc.type.yhi = neumann neumann
+#momentum.bc.val.xlo = "0.1*(y>0.5)*(y<2.5)" "0.0"
+momentum.bc.val.xlo = "0.1 * (0.5*tanh((y - 0.5)/0.03) - 0.5*tanh((y - 2.5)/0.03))" "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 ###
 
@@ -124,7 +125,8 @@ e_refinement_criterion   = 0.1
 m_refinement_criterion   = 0.1
 eta_refinement_criterion = 0.1
 omega_refinement_criterion = 0.1
-p_refinement_criterion   = 1.0
-rho_refinement_criterion   = 0.0
+gradu_refinement_criterion = 0.1
+p_refinement_criterion   = 0.1
+rho_refinement_criterion   = 0.1
 
 #TODO

input.ChannelFlow_eps0.05

@@ -2,18 +2,18 @@ alamo.program = hydro
 
 ### OUTPUT ###
 
-plot_file = ./output.ChannelFlow_eps0.05
+plot_file = ./output.channelflow_eps0.05
 
 ### MESHING ###
 
-amr.plot_dt = 0.25
+amr.plot_dt = 1.0
 #amr.plot_int = 10
-amr.max_level = 2
+amr.max_level = 1
 amr.max_grid_size = 500000
 amr.blocking_factor = 2
 amr.regrid_int = 100
 amr.grid_eff = 0.8
-amr.n_cell = 100 30 0
+amr.n_cell = 200 60 0
 #amr.abort_on_nan=0
 
 ### DIMENSIONS ###
@@ -25,7 +25,7 @@ geometry.is_periodic = 0 0 0
 
 ### TIME STEP ###
 
-timestep  = 1.0e-3
+timestep  = 2.0e-3
 stop_time = 70.0
 
 ### HYDRO PARAM ###
@@ -57,7 +57,7 @@ density.ic.expression.region0 = "100.0"
 
 velocity.ic.type = expression
 velocity.ic.expression.constant.eps = 0.05
-velocity.ic.expression.region0 = "0.01 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
+velocity.ic.expression.region0 = "0.001 * (0.5*tanh((y - 0.5)/eps) - 0.5*tanh((y - 2.5)/eps))"
 velocity.ic.expression.region1 = "0.0"
 
 pressure.ic.type = expression
@@ -113,28 +113,28 @@ momentum.bc.type.xlo = dirichlet dirichlet
 momentum.bc.type.xhi = neumann neumann
 momentum.bc.type.ylo = neumann neumann
 momentum.bc.type.yhi = neumann neumann
-#momentum.bc.val.xlo = "1.0*(y>0.5)*(y<2.5)" "0.0"
+#momentum.bc.val.xlo = "0.1*(y>0.5)*(y<2.5)" "0.0"
 
-momentum.bc.val.xlo = "1.0 * (0.5*tanh((y - 0.5)/0.05) - 0.5*tanh((y - 2.5)/0.05))" "0.0"
+momentum.bc.val.xlo = "0.1 * (0.5*tanh((y - 0.5)/0.05) - 0.5*tanh((y - 2.5)/0.05))" "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   = 1E100
-e_refinement_criterion   = 1E100
-m_refinement_criterion   = 1E100
-eta_refinement_criterion = 0.001
-omega_refinement_criterion = 1E100
-p_refinement_criterion   = 1E100
-rho_refinement_criterion   = 1E100
+r_refinement_criterion   = 0.1
+e_refinement_criterion   = 0.1
+m_refinement_criterion   = 0.1
+eta_refinement_criterion = 0.1
+omega_refinement_criterion = 0.1
+gradu_refinement_criterion = 0.1
+p_refinement_criterion   = 0.1
+rho_refinement_criterion   = 0.1
 
 #TODO
-
-dynamictimestep.on = 0
-dynamictimestep.verbose = 1
-dynamictimestep.cfl = 0.001
-dynamictimestep.max  = 1.0
-dynamictimestep.min  = 1E-8
-cfl = 0.001
+#
+#dynamictimestep.on = 1
+#dynamictimestep.verbose = 1
+#dynamictimestep.cfl = 0.0004
+#dynamictimestep.max  = 1.0
+#dynamictimestep.min  = 0.0