Skip to content

Commit

Permalink
Add inputs file demonstrating inconsistent EBFactory data.
Browse files Browse the repository at this point in the history
  • Loading branch information
cgilet committed Sep 25, 2023
1 parent 6a44ae9 commit deedaf4
Showing 1 changed file with 92 additions and 0 deletions.
92 changes: 92 additions & 0 deletions test_2d/inputs_cyl_right
Original file line number Diff line number Diff line change
@@ -0,0 +1,92 @@
amrex.fpe_trap_invalid = 1

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# SIMULATION STOP #
#.......................................#
max_step = 0 # Max umber of time steps
steady_state = 0 # Steady-state solver?

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# TIME STEP COMPUTATION #
#.......................................#
incflo.fixed_dt = 1.e-3 # Use this constant dt if > 0
#incflo.cfl = 0.45

incflo.do_initial_proj = 1
incflo.initial_iterations = 0

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# INPUT AND OUTPUT #
#.......................................#
amr.plot_int = 1 # Steps between plot files
amr.plt_regtest = 1
amr.check_int =1

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# PHYSICS #
#.......................................#
incflo.gravity = 0. 0. # Gravitational force (3D)
incflo.ro_0 = 1. # Reference density

incflo.fluid_model = "newtonian" # Fluid model (rheology)
incflo.mu = 0.0 # Dynamic viscosity coefficient
incflo.constant_density = false #

incflo.advect_momentum = true
incflo.advection_type = "MOL"
incflo.diffusion_type = 0
incflo.redistribution_type = "StateRedist"

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# ADAPTIVE MESH REFINEMENT #
#.......................................#
amr.n_cell = 34 16 # 544 256 #1088 512 #272 128 # 136 64 # Grid cells at coarsest AMRlevel
amr.max_level = 0 # Max AMR level in hierarchy
amr.max_grid_size_x = 1024
amr.max_grid_size_y = 1024
amr.max_grid_size_z = 1024
amr.blocking_factor = 2

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# GEOMETRY #
#.......................................#
geometry.prob_lo = 0. 0. # Lo corner coordinates
geometry.prob_hi = 3.4 1.6 # Hi corner coordinates
geometry.is_periodic = 0 1 # Periodicity x y z (0/1)

# Boundary conditions
xlo.type = "mi"
xlo.velocity = 1.0 0.0 0.0
xlo.tracer = 1
xlo.density = 1.

xhi.type = "po"
xhi.pressure = 0.0

# Add box
incflo.geometry = "sphere"

sphere.internal_flow = false
sphere.radius = 0.4
sphere.center = 1.1 0.800 0.0

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# INITIAL CONDITIONS #
#.......................................#
incflo.ic_u = 1.0 #
incflo.ic_v = 0.0 #
incflo.ic_p = 0.0 #

#¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨¨#
# VERBOSITY #
#.......................................#
incflo.verbose = 1 # incflo itself
mac_proj.verbose = 1 # MAC Projector
nodal_proj.verbose = 1 # Nodal Projector

#.......................................#
# EB FLOW #
#.......................................#
eb_flow.velocity = 1.0 0.0 0.0

#eb2.small_volfrac = 1.0e-6

0 comments on commit deedaf4

Please sign in to comment.