A time-domain finite element formulation of the equivalent fluid model for the acoustic wave equation

Sound-absorptive materials such as foam can be described by the equivalent fluid (EF) model. The homogenized fluid’s acoustic behavior is thereby described by complex-valued, frequency-dependent acoustic material parameters. When transforming the acoustic wave equation for the EF model from the frequency domain to the time domain, convolution integrals arise. The auxiliary differential equation (ADE) method is used to circumvent the direct calculation of these convolution integrals. The wave equation and the coupled set of ordinary ADEs are solved in the time domain using the finite element (FE) method. The approach relies on approximating the complex-valued frequency response functions of the inverse equivalent bulk modulus and density by a sum of rational functions consisting of real and complex poles. The order of the rational function approximation defines the number of additionally introduced auxiliary variables per nodal degree of freedom. The presented FE formulation includes a narrow-band non-reflecting boundary condition (NRBC) for normal incidence. The implementation in openCFS shows optimal temporal and spatial convergence for a semi-infinite duct based on the analytic plane wave solution for harmonic excitation. The simulation of a pressure pulse propagating in an infinite EF domain with a scatterer demonstrates the capability for multidimensional, actual transient problems.