In this test case, a slab of explosive in a stick form is confined by an inert material before detonating. This example tests one of the multiphase flow solvers with hybrid formulation and makes use of a JWL equation of state .

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