/MAT/LAW41 (LEE_TARVER)

Format

Definitions

Example (LX17)

#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#-  2. MATERIALS:
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW41/1
LX17 (unit Mg-mm-s)
#              RHO_I               RHO_0
                1900                   0
#    Ireac
         2
#                 Ar                  Br                 R1r                 R2r                 R3r
       4930000000000       -166000000000                7.44                3.72           3.3337E-5
#                 Ap                  Bp                 R1p                 R2p                 R3p
        696000000000          2500000000                 4.4                 .94              4.3E-6
#                Cvr                 Cvp                  Eq
                2781                1000                .088
#     iter                           EPS               check
         0                             0                   0
#                rki                  ex                  ri
           100000000                   1                   4
#                rkg                  yg                  zg                 ex1
          1000000000                .371                   3                .191
#                  K                   X                 tol
                   0                   0                   0
#              grow2                 ex2                 yg2                 zg2
                   0                   1                   1                   1
#              ccrit               fmxig               fmxgr               fmngr
                   0                 .25                   1                 100
#                  G                  Ti
            75000000                 298
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
#ENDDATA
/END
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|

Reference

E.L. Lee and C.M. Tarver "Phenomenological model of shock initiation in heterogeneous explosives" Phy. Fluids Vol. 23, No. 12, December 1980.

Comments

1. If f is the mass fraction of the products and p is the reduced pressure:
   Ireact = 1: "Ignition and growth" according to Lee/Tarver(1)

   $${\frac{\mathrm{df}}{\mathrm{dt}})}_{\mathrm{i}}={\mathit{r}}_{\mathrm{ki}}*{\left(1-\mathrm{f}\right)}^{\mathrm{ex}}*{\left(\frac{\mathrm{\rho }}{{\mathrm{\rho }}_0}-1\right)}^{{\mathit{r}}_{\mathrm{i}}}$$

   (2)

   $${\frac{\mathrm{df}}{\mathrm{dt}})}_{\mathrm{g}}={\mathit{r}}_{\mathrm{kg}}*{\left(1-\mathrm{f}\right)}^{\mathrm{ex}}*{\mathrm{f}}^{{\mathrm{y}}_{\mathrm{g}}}*{\mathrm{f}}^{{\mathrm{z}}_{\mathrm{g}}}$$
   Ireac = 2: "Ignition and growth" according to the formulation introduced in Dyna-2D(3)

   $${\frac{\mathrm{df}}{\mathrm{dt}})}_{\mathrm{i}}={\mathit{r}}_{\mathrm{ik}}*{\left(\mathrm{fmxig}-\mathrm{f}\right)}^{\mathrm{ex}}*{\left(\frac{\mathrm{\rho }}{{\mathrm{\rho }}_0}-1-\mathrm{ccrit}\right)}^{{\mathrm{r}}_{\mathrm{i}}}$$

   (4)

   $${\frac{\mathrm{df}}{\mathrm{dt}})}_{\mathrm{g}1}={\mathrm{grow}}_1*{\left(1-\mathrm{f}\right)}^{{\mathrm{ex}}_1}*{\mathrm{f}}^{{\mathrm{y}}_{\mathrm{g}1}}*{\mathrm{P}}^{{\mathrm{z}}_{\mathrm{g}1}}$$

   (5)

   $${\frac{\mathrm{df}}{\mathrm{dt}})}_{\mathrm{g}2}={\mathrm{grow}}_2*{\left(1-\mathrm{f}\right)}^{{\mathrm{ex}}_2}*{\mathrm{f}}^{{\mathrm{y}}_{\mathrm{g}2}}*{\mathrm{P}}^{{\mathrm{z}}_{\mathrm{g}2}}$$
2. Coefficient grow₁ is initialized by r_kg
3. Coefficients y_g1 and z_g1 are respectively initialized by y_g and z_g.
4. Coefficients R₃ and $\omega$ are linked by the relation: ${\mathrm{R}}_3=\omega {\mathrm{C}}_{\mathrm{\nu }}$
5. Coefficients fmxgr and fmngr are the limiters of the growth rate according to the mass fraction of products.
6. This material law is not compatible with ALE.
7. Heat reaction energy E_Q is supposed to be constant.
8. Reagent pressure P_r and detonation products pressure P_r are computed using a modified Jones-Wilkins-Lee equation of state:
   In term of relative volume $v$ :(6)

   $$P\left(v,T\right)={\mathit{Ae}}^{-R_{1}v}+{\mathit{Be}}^{-R_{2}v}+R_{3}Tv$$

   Where, $v=\frac{V}{V_0}$

   In terms of $\mu$ :(7)

   $$P\left(\mu ,T\right)={\mathit{Ae}}^{-R_1/\left(1+\mu \right)}+{\mathit{Be}}^{-R_2/\left(1+\mu \right)}+R_{3}T/\left(1+\mu \right)$$

   Where, $\mu =\frac{\rho }{{\rho }_0}-1=\frac{1}{v}-1$ .