3.7.1 Volumetric drag

Products: ABAQUS/Standard  ABAQUS/Explicit  

Elements tested

AC1D2    AC1D3   

AC2D3    AC2D4    AC2D4R    AC2D6    AC2D8   

AC3D4    AC3D6    AC3D8    AC3D8R    AC3D10    AC3D15    AC3D20   

ACAX3    ACAX4    ACAX4R    ACAX6    ACAX8   

Features tested

Acoustic analysis in steady-state (direct and subspace-based) and transient analyses with high discontinuity in volumetric drag.

Problem description

The model consists of a tube of fluid 4 m long with a constant cross-sectional area. The tube lies horizontally (along the -axis) and has a sound source at = 0 m, which is given in the form of an inward volume acceleration. From = 0 m to = 3 m, the acoustic material in the tube is air with a bulk modulus of 1.424 × 105 N/m2 and a density of 1.21 kg/m3. The region from = 3 m to = 4 m is filled with a dissipative material with the same bulk modulus and density of air but with a volumetric drag of 10,000 Ns/m4. The condition at = 4 m is a closed end. The tube is modeled using 400 first-order or 200 second-order acoustic elements.

The speed of sound for these air constants is = 343 m/s. At the highest frequency of 1100 Hz the wavelength is 0.312 m. The internodal interval (distance between nodes) for the meshes is always .01 m; therefore, at this frequency there are 30 first-order elements per wavelength or 15 second-order elements.

Both direct-solution and subspace-based steady-state dynamic analyses are performed in ABAQUS/Standard over 3 frequencies ranging from 100 to 1100 Hz. The transient simulations are performed in ABAQUS/Explicit using an excitation frequency of 100 Hz. Different excitation frequencies can be tested by changing the parameters defined in the input files. The transient analysis is also performed in ABAQUS/Standard using the AC2D4 element for the purpose of providing a reference solution for ABAQUS/Explicit.

Results and discussion

For ABAQUS/Standard at the highest frequency the results with the second-order meshes lie within 0.1% of the analytical solution for the pressure and the phase in the air region. With the first-order meshes the results lie within 7%. As is to be expected, the second-order elements perform considerably better than first-order elements for the same number of degrees of freedom. Results for both types of mesh improve at lower frequencies (where there are more elements per wavelength).

The results from the transient analyses in ABAQUS/Explicit agree very well with those obtained from ABAQUS/Standard.

Input files

ABAQUS/Standard input files

ec12afad.inp

AC1D2 elements.

ec13afad.inp

AC1D3 elements.

ec23afad.inp

AC2D3 elements.

ec24afad.inp

AC2D4 elements.

ec26afad.inp

AC2D6 elements.

ec28afad.inp

AC2D8 elements.

ec34afad.inp

AC3D4 elements.

ec36afad.inp

AC3D6 elements.

ec38afad.inp

AC3D8 elements.

ec3aafad.inp

AC3D10 elements.

ec3fafad.inp

AC3D15 elements.

ec3kafad.inp

AC3D20 elements.

eca3afad.inp

ACAX3 elements.

eca4afad.inp

ACAX4 elements.

eca6afad.inp

ACAX6 elements.

eca8afad.inp

ACAX8 elements.

ec12afaf.f

FORTRAN: analytical solution.

ec24afad_trans.inp

AC2D4 elements.

ABAQUS/Explicit input files

eca3afad_trans_xpl.inp

ACAX3 elements.

eca4arad_trans_xpl.inp

ACAX4R elements.

ec23afad_trans_xpl.inp

AC2D3 elements.

ec24arad_trans_xpl.inp

AC2D4R elements.

ec34afad_trans_xpl.inp

AC3D4 elements.

ec36afad_trans_xpl.inp

AC3D6 elements.

ec38arad_trans_xpl.inp

AC3D8R elements.