Slab of explosive confined in 3D cylinder
In this test case, a slab of explosive in a stick form, is confined inside a container. The explosive material is initiated by the booster, which forms a detonation wave that propagates to the right in reactive material (explosives). The high pressure behind the detonation front forces the structure confiner to deformation. The explosive considered in this scenario is LX-17 and the Ignition and Growth (I&G) reaction rate law is assumed.
References:
[1] L.Michael, N.Nikiforakis, A hybrid formulation for the numerical simulation of condensed phase explosives, J. Comp. Phys. 316 (2016) 193-217.
*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_HYBRID
*DUALCESE_INITIAL_HYBRID_SET
*DUALCESE_PART_MULTIPHASE
*DUALCESE_EOS_SET
*DUALCESE_EOS_JWL
*DUALCESE_EOS_IDEAL_GAS
*DUALCESE_REACTION_RATE_IG
*DUALCESE_D3PLOT
*CONTROL_TERMINATION
*DATABASE_BINARY_D3PLOT
*DATABASE_BINARY_D3DUMP
*INCLUDE
*PART
*SECTION_SOLID
*MAT_ELASTIC
*HOURGLASS
*CONTROL_CONTACT
*CONTROL_STRUCTURED
*KEYWORD
*TITLE
Rate Stick in 3D container
$
$===========================================================================
$ Units:
$ rho (density): Kg/m^3
$ V (velocity): m/s
$ P (pressure): Pa
$ T (temparature): K
$============================================================================
*CONTROL_TERMINATION
$ endtim endcyc dtmin endeng endmas
15.e-6 20000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Handle output of state and restart data
$
*DATABASE_BINARY_D3PLOT
$ dt/cycl lcdt beam npltc
5.0e-7 0
*DATABASE_BINARY_D3DUMP
$ dt/cycl lcdt beam npltc
10000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*INCLUDE
struct_setup.k
$
*DUALCESE_MODEL
dcese_setup.k
*END
$---------------------------------------------------------------------------------------------------------------------------------------------
*KEYWORD
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DUALCESE_CONTROL_SOLVER
$ icese igeom iframe MixSelect
Euler 3D fixed hybrid
*DUALCESE_CONTROL_TIMESTEP
$ iddt cfl dtint
2 .9 9.e-9
*DUALCESE_CONTROL_LIMITER
$ idlmt alfa beta epsr
2 1.0 1.0 .0
$
*DUALCESE_INCLUDE_MODEL
fluid_mesh1.k
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup the boundary conditions for fluid
$
*DUALCESE_BOUNDARY_NON_REFLECTIVE_SEGMENT_SET
$ ssid
1
2
$
$ Solid BCs (up & down)
$
*DUALCESE_BOUNDARY_REFLECTIVE_SEGMENT_SET
$*DUALCESE_BOUNDARY_SOLID_WALL_SEGMENT_SET
$ ssid
5
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup the initial conditions for fluid
$
*DUALCESE_INITIAL_HYBRID
$ z1 ra uic vic wic rho_1 rho_a rho_b
0.9999999 0.0 0.0 0.0 0.0 1.225 1905.0 1905.0
$ pic tic
1.01325e5
$ 1.1236e8
$
*DUALCESE_INITIAL_HYBRID_SET
$ ssid idef_func
62
$ z1 ra uic vic wic rho_1 rho_a rho_b
1.0e-7 1.0 0.0 0.0 0.0 1905.0 1905.0 1905.0
$ pic tic
0.33708e11
*DUALCESE_INITIAL_HYBRID_SET
$ ssid idef_func
63
$ z1 ra uic vic wic rho_1 rho_a rho_b
1.0e-7 1.0 0.0 0.0 0.0 1905.0 1905.0 1905.0
$ pic tic
1.1236e8
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$ Setup fluid properties
$
*DUALCESE_PART_MULTIPHASE
$ pid ReactinID EosSet_id mid FSI-type MOVMSHALG
$ 1 3 11
1 3 11 ibm
*DUALCESE_EOS_SET
$ Setid inert_eosid react_eosid product_eosid
11 5 6 7
*DUALCESE_EOS_IDEAL_GAS
$ eosid cv cp t0 v0
5 713.5 1001.5
*DUALCESE_EOS_JWL
$ eosid A B R1 R2 Gama0 rho0 e0
6 778.1e11 -5.031e9 11.3 1.13 0.8938 1905.0 0.2552e6
$ Cv
1.3055e3
*DUALCESE_EOS_JWL
$ eosid A B R1 R2 Gama0 rho0 e0
714.8105e11 0.6379e11 6.2 2.2 0.5 1905.0 0.3622e7
$ Cv
.524934e3
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*DUALCESE_REACTION_RATE_IG
$ react_id Ign aa bb xx Grow1 cc dd
3 4.0e12 0.22 0.667 7.0 4.5e-24 0.667 1.0
$ yy Grow2 ee gg zz IGmax G1max G2min
3.0 3.0e-4 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
volume_fraction
reactant_mass_fraction
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*END
$---------------------------------------------------------------------------------------------------------------------------------------------
*KEYWORD
*TITLE
struct_shell
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*PART
solid
$ PID SECID MID EOSID HGID GRAV ADPOPT TMID
11 6 5 0 1
*SECTION_SOLID
6 1
*MAT_ELASTIC
5,1.0e+4,2.0e+11,0.1,0.0,0.0,0.0
*HOURGLASS
1 4
$-------------------------------------------------------------------------------
*CONTROL_CONTACT
0.10000
2
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*INCLUDE
struct_mesh1.k
*CONTROL_STRUCTURED
*END
