SummaryThe equation of state model
Ideal GasWrite an ideal gas equation of state in the form of where is the ambient pressure, is the gas constant, is the current temperature, and is the absolute zero on the temperature scale used. It is an idealization to real gas behavior and can model any gases approximately under appropriate conditions (for example, low pressure and high temperature).
JWLThe Jones-Wilkins-Lee (or JWL) equation of state models the pressure generated by the release of chemical energy in an explosive. This model is implemented in a form referred to as a programmed burn, which means that the reaction and initiation of the explosive is not determined by shock in the material. Instead, the initiation time is determined by a geometric construction using the detonation wave speed and the distance of the material point from the detonation points. You can write the JWL equation of state in terms of the internal energy per unit mass, , as where , where , where , , and are user-defined material constants; is the user-defined density of the explosive; and is the density of the detonation products. Explosive materials generally have some nominal volumetric stiffness before detonation. It may be useful to incorporate this stiffness when elements modeled with a JWL equation of state are subjected to stress before initiation of detonation by the arriving detonation wave. You can define the predetonation bulk modulus, . The pressure is computed from the volumetric strain and until detonation, at which time the pressure is determined by the procedure outlined above. The initial relative density ( ) used in the JWL equation is assumed to be unity. The initial specific energy is assumed to be equal to the user-defined detonation energy . If you specify a nonzero value of , you can also define an initial stress state for the explosive materials.
Currently, no support for specifying coordinates of detonation point is available for this model. Ignition and growthThe ignition and growth description of the equation of state models shock initiation and detonation wave propagation of solid explosives that react to form gaseous products. It requires the specification of the reaction rate to convert the unreacted solid explosive to reacted gas and the specific heat properties of the reacted gas. This model assumes the heterogenous explosive to be made of a homogeneous mixture of unreacted solid explosive and the reacted gaseous products. Separate JWL equations of state can be prescribed for each phase at each material point: where and The subscript refers to the unreacted solid explosive, and refers to the reacted gas products. and are user-defined material constants used in the JWL equations; is the detonation energy; is the user-defined reference density of the explosive, and is the density of the unreacted explosive or the reacted products.
The specific heat of the gas can be specified as a function of temperature and field variables.
The conversion of unreacted solid explosive to reacted gaseous products is governed by the reaction rate. The reaction rate equation in the ignition and growth model is a pressure-driven rule, which includes three terms: These three terms are defined as follows: where , and are reaction rate constants; and is a reference pressure. The first term, , describes hot spot ignition by igniting some of the material relatively quickly but limiting it to a small proportion of the total solid . The second term, , represents the growth of reaction from the hot spot sites into the material and describes the inward and outward grain burning phenomena; this term is limited to a proportion of the total solid . The third term, , is used to describe the rapid transition to detonation observed in some energetic materials.
Us - UpThe Us - Up options allow you to define the Mie-Grüneisen equation of state: where and define the linear relationship between the shock velocity, , and the particle velocity, , as follows: The reference density, , is the material Density, which must be entered separately from the EOS behavior.
TabularThe tabulated equation of state provides flexibility in modeling the hydrodynamic response of materials that exhibit sharp transitions in the pressure-density relationship, such as those induced by phase transformations. The tabulated equation of state is linear in energy and assumes the form where and are functions of the logarithmic volumetric strain only, where , and is the reference density. You can specify the functions and directly in tabular form. The tabular entries must be given in descending values of the volumetric strain; that is, from the most tensile to the most compressive states. The app will use a piecewise linear relationship between data points. Outside the range of specified values of volumetric strains, the functions are extrapolated based on the last slope computed from the data. Specify at least two rows of data, and the values you specify must conform to the following criteria:
Plastic compactionThe Us - Up and the tabular equation of state models also support irreversible compaction behavior during volumetric response for porous solid materials. Typical applications include modeling compaction of granular materials like partially saturated sand during under water explosion or porous metals like aluminum and iron under dynamic loads. ![]() The model requires , an elastic limit for pressure below which compaction is assumed to be reversible and , a maximum pressure for full compaction. The degree of compaction is expressed in terms of distension ( at full compaction), with its initial value in terms of initial porosity . The pressure response as shown in Figure 1 can be separated into a plastic response and several branches of elastic response depending on the deformation history: Here , is the minimum value of representing the maximum compaction achieved by the material based on the deformation history experienced so far whenever . In this sense, is a state variable and determines the relevant branch of the elastic response in Figure 1 on elastic unloading or reloading. The plastic portion of the pressure relation is given by Or, equivalently in terms of the inverse relation Similarly, the elastic portion of the pressure response is given by: Again, equivalently the relation is The slope of the elastic pressure response is given by the expression: where is the elastic bulk modulus of the solid material at small nominal strains; is the reference density of the solid; and and are the reference sound speeds in the solid and virgin (porous) materials, respectively. In these equations, represents the degree of compaction at which the material first experiences irreversible behavior. This parameter is internally calculated by specifying the following nonlinear equation in . If the solid is modeled with the Us - Up model, is equal to the reference sound speed, , which is user specified. On the other hand, if the solid phase is modeled using the tabulated equation of state, is computed from the initial bulk modulus and reference density of the solid material, . In this case the reference density is required to be constant; it cannot be a function of temperature or field variables.
UserThe user option provides a general capability to specify pressure as a function of current density and internal energy per unit mass through the user subroutine VUEOS .
|