Propagation of acoustic waves in a one-dimensional macroscopically inhomogeneous poroelastic material

Bayesian inference of human bone sample properties using ultrasonic reflected signals

The non-intrusiveness and low cost of ultrasonic interrogation is motivating the development of new means of detection of osteoporosis and other bone deficiencies. Bone is a porous media saturated with a viscous fluid and could thus be well characterized by the Biot model. The main purpose of this work is to present an in vitro methodology for the identification of the properties and structural parameters of the bone, adopting a statistical Bayesian inference technique using ultrasonic reflected signals at normal incidence. It is, in this respect, a companion paper to a previous work [J. Acoust. Soc. Am. 146, 3 (2019), pp. 1629-1640], where ultrasonic transmitted signals were considered. This approach allows the retrieval of some important parameters that characterize the bone structure and associated uncertainties. The method was applied to seven samples of bone extracted from femoral heads, immersed in water, and exposed to ultrasonic signals with a center frequency of ≈500 kHz. For all seven samples, signals at different sites were acquired to check the method robustness. The porosity, pore mean size and standard deviation, and the porous frame bulk modulus were all successfully identified using only ultrasonic reflected signals.

Graded and Anisotropic Porous Materials for Broadband and Angular Maximal Acoustic Absorption

The design of graded and anisotropic materials has been of significant interest, especially for sound absorption purposes. Together with the rise of additive manufacturing techniques, new possibilities are emerging from engineered porous micro-structures. In this work, we present a theoretical and numerical study of graded and anisotropic porous materials, for optimal broadband and angular absorption. Through a parametric study, the effective acoustic and geometric parameters of homogenized anisotropic unit cells constitute a database in which the optimal anisotropic and graded material will be searched for. We develop an optimization technique based on the simplex method that is relying on this database. The concepts of average absorption and diffuse field absorption coefficients are introduced and used to maximize angular acoustic absorption. Numerical results present the optimized absorption of the designed anisotropic and graded porous materials for different acoustic targets. The designed materials have anisotropic and graded effective properties, which enhance its sound absorption capabilities. While the anisotropy largely enhances the diffuse field absorbing when optimized at a single frequency, graded properties appear to be crucial for optimal broadband diffuse field absorption.

Acoustic wave propagation in effective graded fully anisotropic fluid layers

This work deals with the sound wave propagation modeling in anisotropic and heterogeneous media. The considered scattering problem involves an infinite layer of finite thickness containing an anisotropic fluid whose properties can vary along the layer depth. The specular transmission and reflection of an acoustic plane wave by such a layer is modeled through the state vector formalism for the acoustic fields. This is solved using three different numerical techniques, namely, the transfer matrix method, Peano series, and transfer Green's function. These three methods are compared to demonstrate the convergence of the numerical solutions. Moreover, the implemented numerical procedures allow the authors to retrieve the internal acoustic fields and show their dependency along with the fluid anisotropic properties. Results are presented to illustrate the changes in absorption that can be achieved by tuning the fluid anisotropy as well as the variation of these properties across the depth of the layer. The results presented are in very good agreement across the different methods. Given that many porous materials can be modeled as equivalent fluids, the results presented show the potential offered by such numerical techniques, and can further give more insight into inhomogeneous anisotropic porous materials.

Estimating the material parameters of an inhomogeneous poroelastic plate from ultrasonic measurements in water

The estimation of poroelastic material parameters based on ultrasound measurements is considered. The acoustical characterisation of poroelastic materials based on various measurements is typically carried out by minimising a cost functional of model residuals, such as the least squares functional. With a limited number of unknown parameters, least squares type approaches can provide both reliable parameter and error estimates. With an increasing number of parameters, both the least squares parameter estimates and, in particular, the error estimates often become unreliable. In this paper, the estimation of the material parameters of an inhomogeneous poroelastic (Biot) plate in the Bayesian framework for inverse problems is considered. Reflection and transmission measurements are performed and 11 poroelastic parameters, as well as 4 measurement setup-related nuisance parameters, are estimated. A Markov chain Monte Carlo algorithm is employed for the computational inference to assess the actual uncertainty of the estimated parameters. The results suggest that the proposed approach for poroelastic material characterisation can reveal the heterogeneities in the object, and yield reliable parameter and uncertainty estimates.

Bayesian inference of a human bone and biomaterials using ultrasonic transmitted signals

Ultrasonic techniques could be good candidates to aid the assessment of osteoporosis detection, due to their non-intrusiveness and low cost. While earlier studies made use of the measured ultrasonic phase velocity and attenuation inside the bone, very few have considered an inverse identification of both the intrinsic pore microstructure and the mechanical properties of the bone, based on Biot's model. The main purpose of this work is to present an in vitro methodology for bone identification, adopting a statistical Bayesian inference technique using ultrasonic transmitted signals, which allows the retrieval of the identified parameters and their uncertainty. In addition to the bone density, Young's modulus and Poisson's ratio, the bone pore microstructure parameters (porosity, tortuosity, and viscous length) are identified. These additional microstructural terms could improve the knowledge on the correlations between bone microstructure and bone diseases, since they provide more information on the trabecular structure. In general, the exact properties of the saturating fluid are unknown (bone marrow and blood in the case of bone study) so in this work, the fluid properties (water) are identified during the inference as a proof of concept.

A transverse isotropic equivalent fluid model combining both limp and rigid frame behaviors for fibrous materials

Due to the manufacturing process, some fibrous materials like glasswool may be transversely isotropic (TI): fibers are mostly parallel to a plane of isotropy within which material properties are identical in all directions whereas properties are different along the transverse direction. The behavior of TI fibrous material is well described by the TI Biot's model, but it requires one to measure several mechanical parameters and to solve the TI Biot's equations. This paper presents an equivalent fluid model that can be suitable for TI materials under certain assumptions. It takes the form of a classical wave equation for the pressure involving an effective density tensor combining both limp and rigid frame behaviors of the material. This scalar wave equation is easily amenable to analytical and numerical treatments with a finite element method. Numerical results, based on the proposed model, are compared with experimental results obtained for two configurations with a fibrous material. The first concerns the absorption of an incident plane wave impinging on a fibrous slab and the second corresponds to the transmission loss of a splitter-type silencer in a duct. Both configurations highlight the effect of the sample orientation and give an illustration of the unusual TI behavior for fluids.

A mixture approach to the acoustic properties of a macroscopically inhomogeneous porous aluminum in the equivalent fluid approximation

The acoustic properties of an air-saturated macroscopically inhomogeneous aluminum foam in the equivalent fluid approximation are studied. A reference sample built by forcing a highly compressible melamine foam with conical shape inside a constant diameter rigid tube is studied first. In this process, a radial compression varying with depth is applied. With the help of an assumption on the compressed pore geometry, properties of the reference sample can be modelled everywhere in the thickness and it is possible to use the classical transfer matrix method as theoretical reference. In the mixture approach, the material is viewed as a mixture of two known materials placed in a patchwork configuration and with proportions of each varying with depth. The properties are derived from the use of a mixing law. For the reference sample, the classical transfer matrix method is used to validate the experimental results. These results are used to validate the mixture approach. The mixture approach is then used to characterize a porous aluminium for which only the properties of the external faces are known. A porosity profile is needed and is obtained from the simulated annealing optimization process.

Acoustic behavior of a rigidly backed poroelastic layer with periodic resonant inclusions by a multiple scattering approach

The acoustic response of a rigidly backed poroelastic layer with a periodic set of elastic cylindrical inclusions embedded is studied. A semi-analytical approach is presented, based on Biot's 1956 theory to account for the deformation of the skeleton, coupling mode matching technique, Bloch wave representation, and multiple scattering theory. This model is validated by comparing the derived absorption coefficients to finite element simulations. Numerical results are further exposed to investigate the influence of the properties of the inclusions (type, material properties, size) of this structure, while a modal analysis is performed to characterize the dynamic behaviors leading to high acoustic absorption. Particularly, in the case of thin viscoelastic membranes, an absorption coefficient larger than 0.8 is observed on a wide frequency band. This property is found to be due to the coupling between the first volume mode of the inclusion and the trapped mode induced by the periodic array and the rigid backing, for a wavelength in the air smaller than 11 times the material thickness.

A time-domain numerical modeling of two-dimensional wave propagation in porous media with frequency-dependent dynamic permeability

An explicit finite-difference scheme is presented for solving the two-dimensional Biot equations of poroelasticity across the full range of frequencies. The key difficulty is to discretize the Johnson-Koplik-Dashen (JKD) model which describes the viscous dissipations in the pores. Indeed, the time-domain version of Biot-JKD model involves order 1/2 fractional derivatives which amount to a time convolution product. To avoid storing the past values of the solution, a diffusive representation of fractional derivatives is used: The convolution kernel is replaced by a finite number of memory variables that satisfy local-in-time ordinary differential equations. The coefficients of the diffusive representation follow from an optimization procedure of the dispersion relation. Then, various methods of scientific computing are applied: The propagative part of the equations is discretized using a fourth-order finite-difference scheme, whereas the diffusive part is solved exactly. An immersed interface method is implemented to discretize the geometry on a Cartesian grid, and also to discretize the jump conditions at interfaces. Numerical experiments are proposed in various realistic configurations.

Acoustics of porous materials with partially opened porosity

A theoretical and experimental study of the acoustic properties of porous materials containing dead-end (or partially opened) porosity was recently proposed by Dupont, Leclaire, Sicot, Gong, and Panneton [J. Appl. Phys. 110, 094903 (2011)]. The present article provides a description of partially opened porosity systems and their numerous potential applications in the general context of the study of porous materials, the classical models describing them, and the characterization techniques. It is shown that the dead-end pore effect can be treated independently and that the description of this effect can be associated with any acoustic model of porous media. Different theoretical developments describing the dead-end porosity effect are proposed. In particular, a model involving the average effective length of the dead-end pores is presented. It is also shown that if the dead-end effect can be treated separately, the transfer matrix method is particularly well suited for the description of single or multilayer systems with dead-end porosity.

Scattering of acoustic waves by macroscopically inhomogeneous poroelastic tubes

Wave propagation in macroscopically inhomogeneous porous materials has received much attention in recent years. For planar configurations, the wave equation, derived from the alternative formulation of Biot's theory of 1962, was reduced and solved recently: first in the case of rigid frame inhomogeneous porous materials and then in the case of inhomogeneous poroelastic materials in the framework of Biot's theory. This paper focuses on the solution of the full wave equation in cylindrical coordinates for poroelastic tubes in which the acoustic and elastic properties of the poroelastic tube vary in the radial direction. The reflection coefficient is obtained numerically using the state vector (or the so-called Stroh) formalism and Peano series. This coefficient can then be used to straightforwardly calculate the scattered field. To validate the method of resolution, results obtained by the present method are compared to those calculated by the classical transfer matrix method in the case of a two-layer poroelastic tube. As an example, a long bone excited in the sagittal plane is considered. Finally, a discussion is given of ultrasonic time domain scattered field for various inhomogeneity profiles, which could lead to the prospect of long bone characterization.

An application of the Peano series expansion to predict sound propagation in materials with continuous pore stratification

This work reports on an application of the state vector (Stroh) formalism and Peano series expansion to solve the problem of sound propagation in a material with continuous pore stratification. An alternative Biot formulation is used to link the equivalent velocity in the oscillatory flow in the material pores with the acoustic pressure gradient. In this formulation, the complex dynamic density and bulk modulus are predicted using the equivalent fluid flow model developed by Horoshenkov and Swift [J. Acoust. Soc. Am. 110(5), 2371-2378 (2001)] under the rigid frame approximation. This model is validated against experimental data obtained for a 140 mm thick material specimen with continuous pore size stratification and relatively constant porosity. This material has been produced from polyurethane binder solution placed in a container with a vented top and sealed bottom to achieve a gradient in the reaction time which caused a pore size stratification to develop as a function of depth [Mahasaranon et al., J. Appl. Phys. 111, 084901 (2012)]. It is shown that the acoustical properties of this class of materials can be accurately predicted with the adopted theoretical model.