Slab of explosive in a stick form
In this test case, a slab of explosive in a stick form is confined by an inert material before detonating. The non-dimensional computational domain is taken to be (x, y) ∈ [0, 20] × [−6, 0] and the unreacted explosive initially takes up the region (x, y) ∈ [4, 20] × [−2, 0]. The explosive is initiated by a booster of the same material, placed at (x, y) ∈ [2, 4] × [−2, 0], and the entire stick is confined by a simple confiner material taking up the remainder of the domain ((x, y) ∈ ([0, 20]×[−6, 0])\ ([2, 20] × [−2, 0])). The simulation is cylindrically symmetric, with the centerline of the rate stick acting as axis of symmetry. This example comes from reference [1] and tests one of the multiphase flow solvers with hybrid formulation.
References:
[1] L.Michael, N.Nikiforakis, A hybrid formulation for the numerical simulation of condensed phase explosives, J. Comp. Phys. 316 (2016) 193-217.
Resulting pressure fringe:
*DUALCESE_MODEL
*DUALCESE_CONTROL_SOLVER
*DUALCESE_CONTROL_TIMESTEP
*DUALCESE_CONTROL_LIMITER
*DUALCESE_INCLUDE_MODEL
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
*DUALCESE_INITIAL_TWO-HYBRID
*DUALCESE_INITIAL_TWO-HYBRID_SET
*DUALCESE_PART_MULTIPHASE
*DUALCESE_EOS_SET
*DUALCESE_EOS_JWL
*DUALCESE_REACTION_RATE_IG
*DUALCESE_D3PLOT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3DUMP
*KEYWORD
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DUALCESE_CONTROL_SOLVER
$ icese igeom iframe MixSelect
$ Euler axis fixed hybrid
Euler 2D fixed hybrid
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
1 .7 1.e-3
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
2 2. 1.0 .0
$
*DUALCESE_INCLUDE_MODEL
fluid_m2h.k
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup the boundary conditions for fluid
$
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
$ ssid
5
2
3
$
$ Solid BCs (up & down)
$
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
$*DUALCESE_BOUNDARY_AXISYMMETRIC_SEGMENT_SET
$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET
$ ssid
1
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup the initial conditions for fluid
$
*DUALCESE_INITIAL_HYBRID
$ z1 ra uic vic wic rho_1 rho_a rho_b
.9999999 0.0 0.0 0.0 0.0 1.0 1.0 1.0
$ pic tic
0.001
$
*DUALCESE_INITIAL_HYBRID_SET
$ ssid idef_func
62
$ z1 ra uic vic wic rho_1 rho_a rho_b
1.0e-8 1.0 0.0 0.0 0.0 1.0 1.0 1.0
$ pic tic
0.3
*DUALCESE_INITIAL_HYBRID_SET
$ ssid idef_func
63
$ z1 ra uic vic wic rho_1 rho_a rho_b
1.0e-8 1.0 0.0 0.0 0.0 1.0 1.0 1.0
$ pic tic
0.001
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup fluid properties
$
*DUALCESE_PART_MULTIPHASE
$ pid ReactinID EosSet_id mid FSI-type MOVMSHALG
1 3 11
*DUALCESE_EOS_SET
$ Setid inert_eosid react_eosid product_eosid
11 5 6 7
*DUALCESE_EOS_JWL
$ eosid A B R1 R2 Gama0 rho0 e0
5 100.4 -0.04 10.0 1.5 0.8 1.0 0.0
$ Cv
0.006569
$ 12.514
*DUALCESE_EOS_JWL
$ eosid A B R1 R2 Gama0 rho0 e0
$ 6 692.51 -0.045 11.3 1.13 0.8938 1.0 0.0
6 692.51 -0.045 11.3 1.13 0.8938 1.0 0.5446e-2
$ Cv
0.0066
$ 12.573
*DUALCESE_EOS_JWL
$ eosid A B R1 R2 Gama0 rho0 e0
$ 7 13.18 0.57 6.2 2.2 0.5 1.0 0.0
7 13.18 0.57 6.2 2.2 0.5 1.0 0.06141
$ Cv
0.0027
$ 5.1435
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_REACTION_RATE_IG
$ react_id Ign aa bb xx Grow1 cc dd
3 4.0e+6 0.22 0.667 7.0 6389.32 0.667 1.0
$ 3 13.18 0.22 0.667 7.0 6389.32 0.667 1.0
$ yy Grow2 ee gg zz IGmax G1max G2min
3.0 33.71 0.667 0.667 1.0 0.02 0.8 0.8
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Handle output of state and restart data
$
*DUALCESE_D3PLOT
density
pressure
velocity
total_energy
internal_energy
temperature
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
