Inciter: Example set up for the four material gas impact problem
This example runs the Multi-material DG hydrodynamics algorithm for the four material gas impact problem from Pandare et al. (Comp. Fluids 2021). It uses a second-order finite volume method for spatial discretization, with THINC material interface reconstruction. The problem is configured using the user-defined problem class with box-type initial conditions, and symmetric boundary conditions on all domain boundaries. For further details about the numerical method see Pandare et al. (Comp. Fluids 2021)
Initial conditionsControl file
# vim: filetype=sh: # This is a comment # Keywords are case-sensitive title "Gas impacting two slabs problem" inciter term 10.0 cfl 0.75 ttyi 25 # TTY output interval scheme p0p1 limiter vertexbasedp1 partitioning algorithm mj end multimat physics veleq depvar u prelax 1 # finite-rate pressure relaxation prelax_timescale 0.1 # time scale for pressure relaxation intsharp 1 # THINC interface reconstruction intsharp_param 2.5 # interface steepness parameter nmat 4 material gamma 50 1.66666667 1.249 1.4 end # ratio of specific heats cv 0.32539682 28.7 200.0 717.5 end # specific heats at constant volume pstiff 0.0 0.0 0.0 0.0 end # stiffness parameter end ic # background materialid 4 end pressure 1.0 end temperature 0.0034843206 end velocity 0.0 0.0 0.0 end # slab-1 box materialid 2 xmin 1.0 xmax 1.1 ymin -0.5 ymax 2.0 zmin -1.0 zmax 1.0 pressure 1.0 temperature 0.0034843206 end # slab-2 box materialid 3 xmin 1.1 xmax 1.3 ymin -0.5 ymax 2.0 zmin -1.0 zmax 1.0 pressure 1.0 temperature 0.0034843206 end # impactor box materialid 1 xmin 0.25 xmax 0.75 ymin -0.5 ymax 1.1 zmin -1.0 zmax 1.0 pressure 2.0 temperature 0.0062717771 velocity 0.2 0.0 0.0 end end end bc_sym sideset 1 2 3 4 5 6 end end end diagnostics interval 250 format scientific error l2 end field_output interval 2500 var elem material_indicator # material indicator function density # bulk density pressure # bulk pressure specific_total_energy # bulk specific total energy x-velocity # bulk x-velocity y-velocity # bulk y-velocity z-velocity # bulk z-velocity end end end
Example run on 16 CPUs
./charmrun +p16 Main/inciter -c mmimpact4mat.q -i gasimpacthalf_ss6_struc24k.exo -v
Setup screen output
< FACTORY >
---- o ----
* Unknowns data layout (CMake: FIELD_DATA_LAYOUT):
--------------------------------------------------
unknown-major
* PDE factory legend, policy codes:
-----------------------------------
Physics configuration (policy code: p)
A - Advection
D - Advection + diffusion
E - Euler
N - Navier-Stokes
V - Velocity equilibrium
Test problem (policy code: t)
U - User-defined
S - Shear-diffusion
V - Vortical flow
N - Nonlinear energy growth
R - Rayleigh-Taylor
T - Taylor-Green
Z - Zalesak's slotted cylinder
G - Advection of 2D Gaussian hump
C - Advection of cylinder
X - Deformation of cylinder in a vortex
F - Shedding flow over triangular wedge
H - Sod shock-tube
O - Rotated Sod shock-tube
B - Sedov blast-wave
I - Interface advection
A - Advection of 2D Gaussian hump for Euler equations
W - Water-air shock-tube
E - Shock He-bubble problem
D - Underwater explosion problem
* Registered PDEs using continuous Galerkin (CG) methods:
---------------------------------------------------------
Unique equation types : 2
With all policy combinations : 30
Compressible single-material flow : p:EN, t:ABFHNORTUV
Transport : p:AD, t:CGSXZ
* Registered PDEs using discontinuous Galerkin (DG) methods:
------------------------------------------------------------
Unique equation types : 3
With all policy combinations : 21
Compressible multi-material flow : p:V, t:DEHIUW
Compressible single-material flow : p:E, t:ABFHNORTUV
Transport : p:A, t:CGSXZ
< PROBLEM >
---- o ----
* Title: Gas impacting two slabs problem
-----------------------------------------
* Discretization parameters:
----------------------------
Discretization scheme : P0P1+RK
Limiter function : VERTEXBASEDP1
PE-locality mesh reordering : false
Operator-access mesh reordering : false
Local time stepping : false
Number of time steps : 18446744073709551615
Start time : 0
Terminate time : 10
CFL coefficient : 0.75
* Partial differential equations integrated (0):
------------------------------------------------
< Compressible multi-material flow >
dependent variable : u
physics : Velocity equilibrium
problem : User-defined
flux : AUSM
number of materials : 4
finite pressure relaxation : 1
pressure relaxation time-scale : 0.100000
interface sharpening : 1
interface sharpening parameter : 2.500000
number of components : 15
start offset in unknowns array : 0
ratio of specific heats : { 50 1.66667 1.249 1.4 }
specific heat at constant volume : { 0.325397 28.7 200 717.5 }
material stiffness : { 0 0 0 0 }
IC background material id : 4
IC background velocity : { 0 0 0 }
IC background pressure : 1
IC background temperature : 0.00348432
IC box 0 : { 1 1.1 -0.5 2 -1 1 }
IC box 0 material id : 2
IC box 0 density : 0
IC box 0 velocity : { }
IC box 0 pressure : 1
IC box 0 internal energy per unit mass : 0
IC box 0 mass : 0
IC box 0 internal energy per unit volume : 0
IC box 0 temperature : 0.00348432
IC box 1 : { 1.1 1.3 -0.5 2 -1 1 }
IC box 1 material id : 3
IC box 1 density : 0
IC box 1 velocity : { }
IC box 1 pressure : 1
IC box 1 internal energy per unit mass : 0
IC box 1 mass : 0
IC box 1 internal energy per unit volume : 0
IC box 1 temperature : 0.00348432
IC box 2 : { 0.25 0.75 -0.5 1.1 -1 1 }
IC box 2 material id : 1
IC box 2 density : 0
IC box 2 velocity : { 0.2 0 0 }
IC box 2 pressure : 2
IC box 2 internal energy per unit mass : 0
IC box 2 mass : 0
IC box 2 internal energy per unit volume : 0
IC box 2 temperature : 0.00627178
* Input/Output filenames and directories:
-----------------------------------------
Input mesh(es) : { gasimpacthalf_ss6_struc24k.exo }
Volume field output file(s) : out.e-s.<meshid>.<numchares>.<chareid>
Surface field output file(s) : out-surf.<surfid>.e-s.<meshid>.<numchares>.<chareid>
History output file(s) : out.hist.{pointid}
Diagnostics file : diag
Checkpoint/restart directory : restart/
* Output intervals:
-------------------
TTY : 25
Field and surface : 2500
Diagnostics : 250
Checkpoint/restart : 1000
* Output fields:
----------------
Elem field(s) : { density material_indicator pressure specific_total_energy x-velocity y-velocity z-velocity }
Solution contours
Bulk density (top-half) and material indicator (bottom-half) contours for the four material gas impact test:
t = 2.0
t = 5.5
t = 8.5