Fluid Cavity Definition

A surface-based fluid cavity:

  • can be used to model a liquid-filled or gas-filled structure;

  • is associated with a node known as the cavity reference node;

  • is defined by specifying a surface that fully encloses the cavity;

  • is applicable only for situations where the pressure and temperature of the fluid in a particular cavity are uniform at any point in time;

  • can be used to model an airbag using the assumptions of an ideal gas mixture under adiabatic conditions; and

  • has a name that can be used to identify history output associated with the cavity.

This page discusses:

Defining the Fluid Cavity

You must associate a name with each fluid cavity.

Specifying the Cavity Reference Node

Every fluid cavity must have an associated cavity reference node. Along with the fluid cavity name, the reference node is used to identify the fluid cavity. In addition, it may be referenced by fluid exchange and inflator definitions. The reference node should not be connected to any elements in the model.

Specifying the Boundary of the Fluid Cavity

The fluid cavity must be completely enclosed by finite elements unless symmetry planes are modeled (see About Surface-Based Fluid Cavities). Surface elements can be used for portions of the cavity surface that are not structural. The boundary of the cavity is specified using an element-based surface covering the elements that surround the cavity with surface normals pointing inward. By default, an error message is issued if the underlying elements of the surface do not have consistent normals. Alternatively, you can skip the consistency checking for the surface normals.

Specifying Additional Volume in a Fluid Cavity

An additional volume can be specified for a fluid cavity. The additional volume will be added to the actual volume when the boundary of the cavity is defined by a specified surface. If you do not specify a surface forming the boundary of the fluid cavity, the fluid cavity is assumed to have a fixed volume that is equal to the added volume. In Abaqus/Standard, along with the added volume, the surface forming the boundary of the fluid cavity must be specified.

Specifying the Minimum Volume

When the volume of a fluid cavity is extremely small, transients in an explicit dynamic procedure can cause the volume to go to zero or even negative causing the effective cavity stiffness values to tend to infinity. To avoid numerical problems, you can specify a minimum volume for the fluid in Abaqus/Explicit. If the volume of the cavity (which is equal to the actual volume plus the added volume) drops below the minimum, the minimum value is used to evaluate the fluid pressure.

You can specify the minimum volume either directly or as the initial volume of the fluid cavity. If the latter method is used and the initial volume of the fluid cavity is a negative value, the minimum volume is set equal to zero.

Defining the Fluid Cavity Behavior

The fluid cavity behavior governs the relationship between cavity pressure, volume, and temperature. A fluid cavity in Abaqus/Standard can contain only a single fluid. In Abaqus/Explicit a cavity can contain a single fluid or a mixture of ideal gases.

Fluid Behavior with a Homogeneous Fluid

To define a fluid cavity behavior made of a single fluid, specify a single fluid behavior to define the fluid properties. You must associate the fluid behavior with a name. This name can then be used to associate a certain behavior with a fluid cavity definition.

Fluid Behavior with a Mixture of Ideal Gases in Abaqus/Explicit

In Abaqus/Explicit you can define a fluid cavity behavior made of multiple gas species. To define a fluid cavity behavior made of multiple gas species, you specify multiple fluid behaviors to define the fluid properties. Specify the names of the fluid behaviors and the initial mass or molar fractions defining the mixture to associate a certain group of behaviors with a fluid cavity definition.

User-Defined Fluid Behavior in Abaqus/Standard

In Abaqus/Standard the fluid behavior can be defined in user subroutine UFLUID.

Defining the Ambient Pressure for a Fluid Cavity

For pneumatic fluids the equilibrium problem is generally expressed in terms of the “gauge” pressure in the fluid cavity (that is, ambient atmospheric pressure does not contribute to the loading of the solid and structural parts of the system). You can choose to convert gauge pressure to absolute pressure p~ as used in the constitutive law. For hydraulic fluids you can define the ambient pressure, which can be used to calculate the pressure difference in the fluid exchange between a fluid cavity and its environment. The pressure value given as degree of freedom 8 at the cavity reference node is the value of the gauge pressure. The ambient pressure, pA, is assumed to be zero if you do not specify it.

Isothermal Process

For hydraulic fluids and pneumatic fluids in problems of long time duration, it is reasonable to assume that the temperature is constant or a known function of the environment surrounding the cavity. In this case the temperature of the fluid can be defined by specifying initial conditions (see Defining Initial Temperatures) and predefined temperature fields (see Predefined Temperature) at the cavity reference node. For a pneumatic fluid the pressure and density of the gas are calculated from the ideal gas law, conservation of mass, and the predefined temperature field.

Defining the Ambient Temperature for a Fluid Cavity

For pneumatic fluids with adiabatic behavior the ambient temperature is needed when the heat energy flow is defined between a single cavity and its environment and the flow definition is based on analysis conditions. The ambient temperature, θA, is assumed to be zero if you do not specify it.

Hydraulic Fluids

The hydraulic fluid model is used to model nearly incompressible fluid behavior and fully incompressible fluid behavior in Abaqus/Standard. Compressibility is introduced by assuming a linear pressure-volume relationship. The required parameters for compressible behavior are the bulk modulus and the reference density. You omit the bulk modulus to specify fully incompressible behavior in Abaqus/Standard. Specifying a high bulk modulus may affect the stable time increment in Abaqus/Explicit. Temperature dependence of the density can be modeled as a thermal expansion of the fluid.

Defining the Reference Fluid Density

The reference fluid density, ρR, is specified at zero pressure and the initial temperature, θI:

ρR=ρ(0,θI).

Defining the Fluid Bulk Modulus for Compressibility

The compressibility is described by the bulk modulus of the fluid:

p=-K(V(p,θ)-V0(θ)V0(θI))=-KρR(ρ-1(p,θ)-ρ0-1(θ)),

where

p

is the current pressure,

θ

is the current temperature,

K

is the fluid bulk modulus,

V(p,θ)

is the current fluid volume,

ρ(p,θ)

is the density at current pressure and temperature,

V0(θ)

is the fluid volume at zero pressure and current temperature,

V0(θI)

is the fluid volume at zero pressure and initial temperature, and

ρ0(θ)

is the density at zero pressure and current temperature.

It is assumed that the bulk modulus is independent of the change in fluid density. However, the bulk modulus can be specified as a function of temperature or predefined field variables.

Defining the Fluid Expansion

The thermal expansion coefficients are interpreted as total expansion coefficients from a reference temperature, which can be specified as a function of temperature or predefined field variables. The change in fluid volume due to thermal expansion is determined as follows:

V0(θ)=V0(θI)[1+3α(θ)(θ-θ0)-3α(θI)(θI-θ0)],

where θ0 is the reference temperature for the coefficient of thermal expansion and α(θ) is the mean (secant) coefficient of thermal expansion.

If the coefficient of thermal expansion is not a function of temperature or field variables, the value of θ0 is not needed.

Thermal expansion can also be expressed in terms of the fluid density:

ρ0(θ)=ρR/[1+3α(θ)(θ-θ0)-3α(θI)(θI-θ0)].

Pneumatic Fluids

Compressible or pneumatic fluids are modeled as an ideal gas (see Equation of State). The equation of state for an ideal gas (or the ideal gas law) is given as

p~=ρR(θ-θZ),

where the absolute (or total) pressure p~ is defined as

p~=p+pA,

and pA is the ambient pressure, p is the gauge pressure, R is the gas constant, θ is the current temperature, and θZ is absolute zero on the temperature scale being used. The gas constant, R, can also be determined from the universal gas constant, R~, and the molecular weight, MW, as follows:

R=R~MW.

Conservation of mass gives the change of mass in the fluid cavity as

m˙=m˙in-m˙out,

where m is the mass of the fluid, m˙in is the mass flow rate into the fluid cavity, and m˙out is the mass flow rate out of the fluid cavity.

Defining the Molecular Weight

You must specify the value of the molecular weight of the ideal gas, MW.

Specifying the Value of the Universal Gas Constant

You can specify the value of the universal gas constant, R~.

Specifying the Value of Absolute Zero

You can specify the value of absolute zero temperature, θZ.

Adiabatic Process

By default, the fluid temperature is defined by the predefined temperature field at the cavity reference node. However, for rapid events the fluid temperature in Abaqus/Explicit can be determined from the conservation of energy assumed in an adiabatic process. With this assumption, no heat is added or removed from the cavity except by transport through fluid exchange definitions or inflators. An adiabatic process is typically well suited for modeling the deployment of an airbag. During deployment, the gas jets out of the inflator at high pressure and cools as it expands at atmospheric pressure. The expansion is so quick that no significant amount of heat can diffuse out of the cavity.

The energy equation can be obtained from the first law of thermodynamics. By neglecting the kinetic and potential energy, the energy equation for a fluid cavity is given by

d(mE)dt=m˙inHin-m˙outHout-W˙-Q˙,

where the work done by the fluid cavity expansion is given as

W˙=pV˙

and Q˙ is the heat energy flow rate due to the heat transfer through the surface of the fluid cavity. A positive value for Q˙ will generate the heat energy flow out of the primary fluid cavity. The specific energy is given by

E=EI+θIθcv(T)dT,

where EI is the initial specific energy at the initial temperature θI, cv is the specific heat at constant volume (or the constant volume heat capacity), which depends only upon temperature for an ideal gas, H is the specific enthalpy, and V is the volume occupied by the gas. By definition, the specific enthalpy is

H=HI+θIθcp(T)dT,

where HI is the initial specific enthalpy at the initial (or reference) temperature θI and cp is the specific heat at constant pressure, which depends only upon temperature for an ideal gas. The pressure, temperature, and density of the gas are obtained by solving the ideal gas law, the energy balance, and mass conservation.

Adiabatic behavior will always be used for the fluid cavity if an adiabatic or coupled procedure is used.

Defining the Heat Capacity at Constant Pressure

You must define the heat capacity at constant pressure for the ideal gas. It can be defined either in polynomial or tabular form. The polynomial form is based on the Shomate equation according to the National Institute of Standards and Technology. The constant pressure molar heat capacity can be expressed as

c~p=a~+b~(θ-θZ)+c~(θ-θZ)2+d~(θ-θZ)3+e~(θ-θZ)2,

where the coefficients a~, b~, c~, d~, and e~ are gas constants. These gas constants together with molecular weight are listed in Table 1 for some gases that are often used in airbag simulations. The constant pressure heat capacity can then be obtained by

cp=c~pMW.

The constant volume heat capacity, cv, can be determined by

cv=cp-R.
Table 1. Properties of some commonly used gases (SI units).
Gas MW a~ b~ c~ d~ e~ θ
(× 10−3) (× 10−6) (× 10−9) (× 106) (kelvin)
Air 0.0289 28.110 1.967 4.802 −1.966 0.0 273–1800
Nitrogen 0.028 26.092 8.218 –1.976 0.1592 0.0444 298–6000
Oxygen 0.032 29.659 6.137 –1.186 0.0957 –0.219 298–6000
Hydrogen 0.00202 33.066 −11.36 11.432 –2.772 –0.158 273–1000
Carbon monoxide 0.028 25.567 6.096 4.054 −2.671 0.131 298–1300
Carbon dioxide 0.044 24.997 55.186 −33.691 7.948 –0.136 298–1200
Water vapor 0.0180 32.240 1.923 0.105 −3.595 0.0 273–1800

You can use the polynomial form for specifying the heat capacity at constant pressure, in which case you enter the coefficients a~, b~, c~, d~, and e~. Alternatively, you can define a table of constant pressure heat capacity versus temperature and any predefined field variables.

A Mixture of Ideal Gases

Abaqus/Explicit can model a mixture of ideal gases in the fluid cavity. For ideal gas mixtures the Amagat-Leduc rule of partial volumes is used to obtain an estimate of the molar-averaged thermal properties (Van Wylen and Sonntag, 1985). Let each species have constant pressure and volume heat capacities, cpi and cvi; molecular weight, MWi; and mass fraction, fi. The constant pressure and volume heat capacities for the mixed gas are then given by

cp=ificpi,
cv=ificvi,

and the molecular weight is given by

MW=1/ifiMWi.

The specific energy and enthalpy for the mixed gas are then given by

E=ifiEi,
H=ifiHi.

The energy flow entering the fluid cavity is given by

m˙inHin=im˙iniHini,

and the energy flow out of the fluid cavity is given by

m˙outHout=im˙outiHouti.

Using the properties of a mixture of ideal gases as given above, the pressure and temperature can be obtained from the ideal gas law and the energy equation.

Averaged Properties for Multiple Fluid Cavities

If the output of the state of the fluid inside the cavity is requested for a node set that contains more than one fluid cavity, the averaged properties of the multiple fluid cavities will also be output automatically. The average pressure is calculated by volume weighting cavity pressure contributions. The average temperature for an adiabatic ideal gas is obtained by mass weighting cavity temperature contributions. Let each fluid cavity have pressure pk, temperature θk, volume Vk, gas constant Rk, and mass mk. The average pressure of the fluid cavity cluster is then defined as

pavg=kpkVk/kVk,

and the average temperature is

θavg=kRkθkmk/kRkmk.

References

  1. Van Wylen G. J. and RESonntag, Fundamentals of Classical Thermodynamics, Wiley, New York, 1985.