Next Article in Journal
Symmetry GAN Detection Network: An Automatic One-Stage High-Accuracy Detection Network for Various Types of Lesions on CT Images
Next Article in Special Issue
j-Dimensional Integral Involving the Logarithmic and Exponential Functions: Derivation and Evaluation
Previous Article in Journal
Equiauxetic Hinged Archimedean Tilings
Previous Article in Special Issue
Double Diffusive Magneto-Free-Convection Flow of Oldroyd-B Fluid over a Vertical Plate with Heat and Mass Flux
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transient Propagation of Spherical Waves in Porous Material: Application of Fractional Calculus

1
Laboratory of Mechanics and Acoustics, French National Centre for Scientific Research LMA, CNRS, UMR 7031, Centrale Marseille, Aix-Marseille University CEDEX 20, F-13402 Marseille, France
2
Laboratory of Theoretical Physics, Faculty of Physics, University of Science and Technology Houari-Boumediene BP 32 El Alia, Bab Ezzouar 16111, Algeria
3
French Aerospace Lab, ONERA/Multi-Physics Department for Energy, Toulouse University, F-31055 Toulouse, France
4
Department of Physics and Materials Science, Maseno University, Maseno 40105, Kenya
5
Acoustics Laboratory of the University of Le Mans, French National Centre for Scientific Research, CNRS UMR 6613, UFR STS Avenue O. Messiaen, CEDEX 09, F-72085 Le Mans, France
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(2), 233; https://doi.org/10.3390/sym14020233
Submission received: 18 December 2021 / Revised: 18 January 2022 / Accepted: 19 January 2022 / Published: 25 January 2022

Abstract

:
A fractional-order wave equation is established and solved for a space of three dimensions using spherical coordinates. An equivalent fluid model is used in which the acoustic wave propagates only in the fluid saturating the porous medium; this model is a special case of Biot’s theory obtained by the symmetry of the Lagrangian (invariance by translation and rotation). The basic solution of the wave equation is obtained in the time domain by analytically calculating Green’s function of the porous medium and using the properties of the Laplace transforms. Fractional derivatives are used to describe, in the time domain, the fluid–structure interactions, which are of the inertial, viscous, and thermal kind. The solution to the fractional-order wave equation represents the radiation field in the porous medium emitted by a point source. An important result obtained in this study is that the solution of the fractional equation is expressed by recurrence relations that are the consequence of the modified Bessel function of the third kind, which represents a physical solution of the wave equation. This theoretical work with analytical results opens up prospects for the resolution of forward and inverse problems allowing the characterization of a porous medium using spherical waves.

1. Introduction

Acoustic absorbers, such as plastic foams or fibrous materials [1,2,3,4], are widely used for sound insulation, passive sound control [1], and noise reduction in the aeronautical and building industries [5]. These materials have the particularity of being porous [6,7] and consisting of a solid phase and a fluid phase, the fluid part generally being saturated with air. When an acoustic (sound) wave propagates in this type of medium, the solid generally remains motionless (does not vibrate), and the wave essentially propagates in the fluid that saturates the porous material [1,6,7]. The equivalent fluid model is used, which is a special case of the Biot theory obtained by the symmetry of the Lagrangian (by translation and rotation) [1].
Acoustic propagation in porous media involves geometrical parameters that are different depending on the frequency domain. For example, in high frequencies, there is the tortuosity and the viscous and thermal characteristic lengths [8,9,10,11,12], but also the viscous and thermal surfaces for the most absorbent media [13,14]. At low frequencies, the most relevant parameters are the viscous and thermal permeabilities [15,16,17], the low-frequency inertial parameter, and a viscous parameter that was recently experimentally measured [18]. It is therefore important to model acoustic propagation in porous media in these two frequency regimes.
The modeling of acoustic propagation in porous materials allows us to predict sound absorption and acoustic attenuation [1]. Several methods have been developed for the resolution of the direct and inverse problems allowing the characterization of porous materials [2,4,9,10,14,16,17,18,19,20,21,22,23,24]. The direct problem consists in solving the propagation equation and deducing the reflected and transmitted acoustic fields [11,12,15,19,25,26,27,28]. Knowing the experimental incident signal and the solution of the propagation equation (Green function of the medium) makes it possible to simulate the acoustic field propagating in the porous material, as well as to predict the waves reflected and transmitted by a slab of material [19]. The solution of the propagation equation is therefore a key step in directly solving the problem [11,12]. To characterize the porous material using sound waves, the solution of the inverse problem becomes necessary. The inverse problem is posed as follows: Knowing the experimental incident signal, the reflected and/or transmitted signal, we go back by inversion to the geometric parameters describing the porous medium [2,4,14,16,19,20,21,24].
In the study of wave propagation phenomena, the search for the elementary solution to the equation of propagation is a fundamental problem [1,11,12]. In spherical coordinates [29], this solution represents the field radiating from a point source. The solution is also essential in perturbation theory [30] when one is interested in propagation in a weakly inhomogeneous (with punctual defects) medium. In this work, we establish and resolve a fractional propagation equation in the time domain for a three-dimensional space by using the properties of the Laplace transform. The problem of the propagation of spherical waves in a homogeneous and isotropic porous medium in the presence of losses is studied. The visco-thermal losses are modeled by fractional derivatives in the time domain. This work generalizes our previous approach, which was to solve the wave equation in one space dimension [11,12]. We show in this paper, by using analytical computational developments, that Green’s function of the porous medium (in spherical coordinates) and the solution of its fractional wave equation can be expressed by recurrence relations and depend on Green’s function of a unidirectional propagation problem and its fractional wave equation solution. These results open up new perspectives for the prediction of sound absorption in porous media and their ultrasonic characterization, since the resolution of the propagation equation is the basis for solving the direct and inverse problems.

2. Fractional Model for Porous Media

When porous media are saturated with air, solid-frame vibrations can be neglected in the absence of direct contact with the sound source [1]. It is then considered that the waves propagate only in the fluid saturating the porous medium [6,7]. We will use the equivalent fluid model, which is a special case of the Biot model in which the exchanges between fluid and structure are described by two relaxation functions: the dynamic tortuosity α ( ω ) given by Johnson et al. [6] and the dynamic compressibility β ( ω ) given by Allard [1]. Biot’s equations [31,32] are Euler’s equations deduced from a Lagrangian taking into account the symmetry due to invariances by translation and rotation.
The dynamic tortuosity describes the inertial and viscous effects while the dynamic compressibility describes the thermal effects. In the frequency domain, these functions correct for the density of the fluid, as well as its compressibility, and they represent the deviation from free-space fluid behavior as the frequency changes. In this context, the equations for acoustics in a porous material are given by:
ρ α ( ω ) i ω V ( r , ω ) = P ( r , ω ) , β ( ω ) K i ω P ( r , ω ) = . V ( r , ω ) .
In Equation (1): i 2 = 1 , ρ is the saturating fluid density and K is the compressibility modulus of the fluid. P ( r , ω ) and V ( r , ω ) are the Fourier transforms of v ( r , t ) and p ( r , t ) . The expression of the dynamic compressibility is given by:
β ( ω ) = γ γ 1 α ( ω ) ,
where α ( ω ) is the dynamic thermal tortuosity and γ is the heat capacity ratio (equal to 1.4 for air). In the high-frequency domain, which corresponds to viscous and thermal boundary layers ( δ ( ω ) = 2 η ρ ω and δ ( ω ) = 2 η ρ ω P r , respectively) that are small compared with the characteristic pore radius, we have [13,14]:
α ( ω ) = α + 2 α Λ η i ω ρ 1 / 2 + 3 α Σ η i ω ρ + O 1 i ω 3 / 2 ,
α ( ω ) = 1 + 2 Λ η i ω ρ P r 1 / 2 + 3 Σ η i ω ρ P r + O 1 i ω ρ P r 3 / 2 ,
where: η is the fluid viscosity, ω is the frequency pulsation, α is the tortuosity, Λ is the viscous characteristic length, Λ is the thermal characteristic length, Σ is the characteristic viscous surface, Σ is the characteristic thermal surface, and P r is the Prandtl number. From Equations (2) and (4), we can deduce the expression of β ( ω ) :
β ( ω ) = 1 + ( γ 1 ) 2 Λ η i ω ρ P r 1 / 2 + 3 Σ η i ω ρ P r + O 1 i ω ρ P r 3 / 2 ,
In the time domain, α and β become time operators whose expressions depend on fractional derivatives; Equations (3) and (5) are then in the time domain:
α ˜ ( t ) = α δ ( t ) + 2 α Λ η ρ 1 / 2 1 / 2 t 1 / 2 + 3 α Σ η ρ 1 t 1 + O 3 / 2 t 3 / 2 ,
β ˜ ( t ) = δ ( t ) + 2 ( γ 1 ) Λ η P r ρ 1 / 2 1 / 2 t 1 / 2 + 3 ( γ 1 ) Σ η ρ P r 1 t 1 + O 3 / 2 t 3 / 2 ,
In Equations (6) and (7), δ ( t ) is the Dirac function. In this model, 1 t 1 represents a temporal integral, and the term 1 / 2 t 1 / 2 is interpreted as a semi-derivative operator following the definition of the fractional derivative of order ν given in Samko et al. [33]:
ν t ν [ x ( t ) ] = 1 Γ ( ν ) 0 t ( t u ) ν 1 x ( u ) d u ,
where Γ ( x ) is the gamma function. In the time domain, the constitutive relations (1) in a porous medium are given by
0 t ρ α ˜ ( t t ) v ( r , t ) t d t = p ( r , t ) ,
0 t 1 K β ˜ ( t t ) p ( r , t ) t d t = K . v ( r , t ) .
Using the expressions (6) and (7) of the tortuosity and compressibility operators and the constitutive relations (9) and (10), we obtain the following relations:
ρ α v ( r , t ) t + 2 α ρ η Λ 1 / 2 v ( r , t ) t 1 / 2 + 3 α Σ η ρ v ( r , t ) = p ( r , t ) ,
1 K p ( r , t ) t + 2 γ 1 K Λ η P r ρ 1 / 2 1 / 2 p ( r , t ) t 1 / 2 + 3 ( γ 1 ) Σ η ρ P r p ( r , t ) = . v ( r , t ) .
Equations (11) and (12) generalize those obtained in previous work for a unidirectional propagation problem [11,12,19]. Moreover, here, the characteristic viscous and thermal surfaces are taken into account, which is not the case in the previous developments. The use of the fractional derivative in the time domain makes it possible to describe the viscous and thermal exchanges between fluid and structure in the constitutive equations. These losses are at the origin of the attenuation and the acoustic absorption in a porous medium.
Using the general constitutive relations (11) and (12), the propagation equation in the time domain is given by:
Δ p ( r , t ) 1 c 2 2 p ( r , t ) t 2 b 3 / 2 p ( r , t ) t 3 / 2 a p ( r , t ) t = 0 ,
where c = K / ρ α is the front wave velocity, a = 1 c 2 η ρ 3 ( γ 1 ) Σ P r + 4 ( γ 1 ) Λ Λ P r + 3 Σ and b = 2 α ρ η K 1 Λ + ( γ 1 ) P r Λ are terms of attenuation and dispersion, which depend mainly on the tortuosity α and viscous and thermal characteristic lengths Λ and Λ . In Equation (13), we have stopped at the term / t , which corresponds to the term ( 1 / i ω ) in the expansion of Equations (4) and (5). This fractional wave equation describes the acoustic propagation in a porous material with a rigid frame. The viscous and thermal losses are modeled using the fractional derivatives in the time domain of order 3 / 2 . Using the definition of the fractional derivative (8), the fractional propagation equation can be written as
Δ p ( r , t ) 1 c 2 2 p ( r , t ) t 2 a p ( r , t ) t b 0 t 2 p ( r , t ) t 2 d t t t = 0 .
The fractional differential Equation (14) describes the propagation of a spherical wave generated by a point-like source, with the spherical wavefronts (connected parts of the wave that all share the same displacement, such as the crests of water waves) spreading radially outward from the point source. The porous medium is isotropic and homogeneous. This propagation equation was solved for the case of unidirectional propagation [11,12]. We will try in this paper to solve it as a three-dimensional problem using spherical coordinates.

3. Solution of the Fractional Spherical-Wave Equation

We propose in this section the solution the fractional integro-differential equation describing the propagation of transient ultrasonic waves in porous materials with a rigid frame given by Equation (13) or by:
r = r ( x , y , z ) ; Δ = 2 x 2 + 2 y 2 + 2 z 2 .
We will solve this fractional wave equation with the method of the Laplace transform. We note P ( r , s ) , the Laplace transform of p ( r , t ) , defined by [34]
P ( r , s ) = L p ( r , t )   = 0 exp ( s t ) p ( r , t ) d t .
The Laplace transform of Equation (14) yields
Δ P ( r , s ) s 2 1 c 2 + a s + b π s P ( r , s ) = 1 c 2 + b π s × s p ( r , 0 ) + p t ( r , 0 ) a p ( r , 0 ) .
The initial conditions are assumed as:
p ( r , 0 ) = p ( r , 0 ) t = 0 .
Equation (16) becomes:
Δ P ( r , s ) s 2 1 c 2 + a s + b π s P ( r , s ) = 0 .
By putting:
k 2 ( s ) = 1 c 2 s 2 + b s s + c s , b = b c 2 π , c = a c 2 ,
we obtain the following Helmholtz equation:
Δ k 2 ( s ) P ( r , s ) = 0 .
Using the spherical coordinates r = r ( r , θ , φ ) , the Laplacian is given by:
Δ = P r 2 + 2 r 2 ,
where
P r = 1 r r r = r + 1 r ,
is an operator depending only on the radial variable r, and is an operator depending only on the angular variables θ and φ and is given by:
2 = 1 sin 2 θ sin θ θ sin θ θ + 2 φ 2 = 2 θ 2 + cos θ sin θ θ + 1 sin 2 θ 2 φ 2 .
It is also in Cartesian coordinates:
= r , where t = x , y , z and   by   noting t = l x , l y , l z ,
we have:
x = y z z y = sin φ θ + cot θ cos φ φ , y = z x x z = cos φ θ cot θ sin φ φ , z = x y y x = φ .
Equation (20) is, in fact, an eigenvalue equation, which can be written as
H P ( r , s ) = 0 ,
where H is an operator given by
H = Δ k 2 ( s ) = P r 2 + 2 r 2 k 2 ( s ) .
It can be seen from Equations (22) and (24) that H and 2 commute because H commutes with l x , l y , and l z , and thus commutes with any function of these operators and, in particular, 2 . The operators H, 2 , and l z have at least one complete set of eigenfunctions in common. We can therefore find the eigenfunctions of H from those of 2 and l z . However, these are none other than the spherical harmonics Y m ( θ , φ ) and are such that
2 Y m ( θ , φ ) = ( + 1 ) Y l m ( θ , φ ) , l z Y l m ( θ , φ ) = i m Y m ( θ , φ ) , i 2 = 1 ,
where l is all positive or zero integers, = 0 , 1 , 2 , . . , and m is all integer values in the range ( , + )   m + .
We look for a solution in the form
P ( r , s ) = P ( r , θ , φ , s ) = R ( r , s ) Y m ( θ , φ ) ; l = 0 , 1 , 2 , , and m + .
By using the two Equations (23) and (26), we find that taking (25) into account satisfies the differential equation:
d 2 d r 2 + 2 r d d r ( + 1 ) r 2 k 2 ( s ) R ( r , s ) = 0 ,
which is a Bessel differential equation. Its special solutions are modified spherical Bessel functions of the first kind, π 2 k r I + 1 / 2 ( k r ) , of the second kind, π 2 k r I 1 / 2 ( k r ) , and of the third kind, π 2 k r K + 1 / 2 ( k r ) [35].
Two linearly independent solutions are formed by the two functions π 2 k r I + 1 / 2 ( k r ) and π 2 k r I 1 / 2 ( k r ) , or by both functions: π 2 k r I + 1 / 2 ( k r ) and π 2 k r K + 1 / 2 ( k r ) . Keeping only the solutions that are finite at infinity, the general solution of (27) is given by:
U ( r , s ) = φ ( s ) R ( r , s ) = φ ( s ) K + 1 / 2 ( r , s ) ; l = 0 , 1 , 2 ,
where φ ( s ) is a constant of integration determined by the initial conditions. It is known that the Bessel function can be calculated, for example, from the recurrence relation [35]:
R ( r , s ) = R 2 ( r , s ) + ( 2 1 ) k r R 1 ( r , s ) .
where R 0 ( r , s ) and R 1 ( r , s ) are given by:
R 0 ( r , s ) = π 2 r k ( s ) K 1 / 2 ( r k ( s ) ) = π 2 r k ( s ) exp r k ( s ) ,
R 1 ( r , s ) = π 2 r k ( s ) K 3 / 2 ( r k ( s ) ) = π 2 r k ( s ) 1 + 1 r k ( s ) exp r k ( s ) .

4. Green Function of the Porous Medium in the Time Domain

Determining Green’s function is an important step [11,12,19] in finding the solution to the fractional wave equation. Indeed, Green’s function, sometimes called a “propagator", describes the propagation of a wave in the porous medium (velocity, attenuation, and dispersion) and depends on the wave number k ( s ) (Equation (19)). The solution of the wave equation is thus obtained by carrying out a time convolution product of the incident signal with Green’s function of the medium.
The Green function of the porous medium in the time domain is given by taking the inverse Laplace transform to solution (27):
G ( r , t ) = L 1 R ( r , s ) .
The inverse Laplace transforms of the relations (29)–(31) are used to calculate G ( r , t ) from G 0 ( r , t ) and G 1 ( r , t ) . Indeed, the Green function of the porous medium G ( r , t ) satisfies the recurrence relation:
G ( r , t ) = G 2 ( r , t ) + ( 2 1 ) r L 1 1 k ( s ) G 1 ( r , t ) ,
with G 0 ( r , t ) and G 1 ( r , t ) such that:
G 0 ( r , t ) = π 2 r L 1 1 k ( s ) L 1 exp r k ( s ) ,
and
G 1 ( r , t ) = π 2 r L 1 1 k ( s ) + 1 r L 1 1 k 2 ( s ) L 1 exp r k ( s ) .
In Equations (33)–(35), * denotes the time convolution product.
The term L 1 exp r k ( s )   = F ( r , t ) corresponds to the Green function of the porous medium for a one-dimensional problem, which has already been calculated analytically in our previous work [11,12] and is given by:
F ( r , τ ) = 0 , if 0 τ r / c , r c b 4 π τ r / c 3 / 2 exp b 2 r 2 16 c 2 ( τ r / c ) + Δ 0 τ r / c h ( ξ ) d ξ , if t r / c ,
where h ( ξ ) has the following form:
h ( ξ ) = 1 π 3 / 2 1 ( τ ξ ) 2 r 2 / c 2 1 ξ 0 1 e x p χ ( μ , τ , ξ ) 2 ( χ ( μ , τ , ξ ) 1 ) μ d μ 1 μ 2 ,
and where χ ( μ , τ , ξ ) = Δ μ ( τ ξ ) 2 r 2 / c 2 + b ( τ ξ ) 2 / 8 ξ , Δ = b 2 4 c .

5. Analytical Solution in the Time Domain

The solution to Equations (26) and (27) in the time domain is obtained by the time convolution product of the Green functions given by relations (33)–(35) with the incident signal p ( 0 , t ) = L 1 φ ( s ) ,
U ( r , t ) = L 1 U ( r , s ) = p ( 0 , t ) G ( r , t ) = 0 t G ( r , τ ) p ( 0 , t τ ) d τ .
We can therefore express the solution U ( r , t ) with the recurrence relations
U ( r , t ) = U 2 ( r , t ) + ( 2 1 ) r L 1 1 k ( s ) U 1 ( r , t ) ,
and
U 0 ( r , t ) = π 2 r L 1 1 k ( s ) 0 t F ( r , τ ) p ( 0 , t τ ) d τ ,
U 1 ( r , t ) = π 2 r L 1 1 k ( s ) + 1 r L 1 1 k 2 ( s ) 0 t F ( r , τ ) p ( 0 , t τ ) d τ .
The convolution product used in Equations (40) and (41) of the Green function F ( r , t ) with the incident signal p ( 0 , t ) gives the solution u ( r , t ) to the propagation equation in a 1D porous medium [11,12]:
u ( r , t ) = 0 t F ( r , τ ) p ( 0 , t τ ) d τ .
It is interesting to see here how the solution of the propagation equation for a spherical wave depends on the solution of the unidirectional propagation problem.
To express U 0 ( r , t ) and U 1 ( r , t ) given by Equations (40) and (41), we need to calculate the inverse Laplace transforms of 1 k ( s ) and 1 k 2 ( s ) , where the wavenumber k ( s ) is given by Equation (19).
Let D ( t ) be the inverse Laplace transform of 1 k ( s ) ; we have (Appendix A):
D ( t ) = L 1 1 k ( s ) = c exp b 2 t J 0 t 2 4 c b 2 , if 4 c b 2 > 0 , c π t 0 exp u 2 4 t exp b u 2 J 0 u 2 4 c b 2 d u , 4 c b 2 < 0 ,
when 4 c b 2 < 0 , we use the relation
J 0 U 2 4 c b 2 = J 0 i U 2 b 2 4 c = I 0 U 2 b 2 4 c ,
where J 0 is the Bessel function and I 0 is the modified Bessel function [34].
To calculate U 1 ( r , t ) given by Equation (48), we need to calculate the inverse Laplace transform of 1 k 2 ( s ) . Let us put S ( t ) = L 1 1 k 2 ( s ) ; according to the calculations made in Appendix B, one has
S ( t ) = L 1 1 k 2 ( s ) = 2 c b 4 c 0 Erfc u 2 t exp b 2 u sinh u 2 b 2 4 c d u if b 2 4 c > 0 , 2 c 4 c b 2 0 Erfc y 2 t exp b 2 y sin y 2 4 c b 2 d y if 4 c b 2 > 0 .
The solution to Equation (13) is obtained by taking the inverse Laplace transform of Equation (26):
p , m ( r , t ) = L 1 P ( r , θ , φ , s ) = U ( r , t ) Y m ( θ , φ ) , = 0 , 1 , 2 , , m ,
where the radial solution U ( r , t ) is given by Equation (28) and satisfies the recurrence relation:
U ( r , t ) = U 2 ( r , t ) + ( 2 1 ) 0 t D ( τ ) U 1 ( r , t τ ) d τ , = 2 , 3 , .
The expression of D ( t ) is given by (43). The functions U 0 ( r , t ) and U 1 ( r , t ) are given by:
U 0 ( r , t ) = π c 2 r 0 t D ( τ ) u ( r , t τ ) d τ ,
U 1 ( r , t ) = π c 2 r 0 t D ( τ ) + c r S ( τ ) u ( r , t τ ) d τ ,
where S ( t ) is given by Equation (44) and u ( r , t ) by Equation (42).
Using solution (45), it is easy to calculate the first terms for = 0 , which induces m = 0 . Knowing that Y 0 0 ( θ , φ ) = 1 4 π , one has
p 0 , 0 ( r , t ) = U 0 ( r , t ) Y 0 0 ( θ , φ ) = 1 4 π U 0 ( r , t ) .
When = 1 , m = 1 , 0 , 1 , and knowing that
Y 1 0 ( θ , φ ) = 3 4 π cos θ and Y 1 ± 1 ( θ , φ ) = 3 8 π sin θ exp ± i φ ,
we have the following solutions:
p 1 , 0 ( r , t ) = 3 4 π cos θ . U 1 ( r , t ) ,
p 1 , ± 1 ( r , t ) = 3 8 π sin θ . exp ± i φ U 1 ( r , t ) ,
where U 1 ( r , t ) is given by Equation (48).
The higher-order terms of the solution p , m ( r , t ) (Equation (45)) can be calculated from the first two solutions p 0 , 0 ( r , t ) and p 1 , m ( r , t ) given by relations (50) and (49) and the recurrence relation (46), which originates from the properties of the Bessel function (28). The analytical calculation of the inverse Laplace transforms of expressions (43) and (44) makes it possible to obtain an analytical form in the time domain of the solution to the fractional spherical-wave equation. This solution allows us to predict the ultrasonic wave in spherical coordinates propagating in a porous medium by knowing the wave p ( 0 , t ) incident in the medium given by the initial conditions and the properties of the material represented by the wave number k ( s ) , which essentially depends on the porosity, tortuosity, viscous and thermal characteristic lengths, and viscous and thermal surfaces.

6. Discussion and Conclusions

In this paper, we have solved the fractional propagation equation in spherical coordinates. This equation describes ultrasonic propagation in a porous material obeying the equivalent fluid model, in which the acoustic wave propagates only in the fluid saturating the porous medium, with the structure remaining motionless. This model is obtained through symmetry considerations (invariance by translation and rotation) of the Lagrangian of the Biot theory. We have shown that Green’s function of the porous medium (in spherical coordinates) and the solution of its fractional wave equation are, respectively, expressed by recurrence relations and depend on Green’s function of a unidirectional propagation problem and its fractional wave equation solution. Calculations are performed analytically in the Laplace domain, thus obtaining analytical expressions in the time domain of the response of the medium to an acoustic excitation. The time-domain analysis has the advantage of being the most suitable for transient signals, such as pulses or bursts, which are often used in experimental reality for the characterization and resolution of direct and inverse problems.
This essentially theoretical work with analytical results constitutes an important first step for the resolution of the direct problem allowing us to predict the ultrasonic propagation of an acoustic wave in a porous medium in spherical coordinates while knowing the incident wave and the physical properties of the porous medium (porosity ϕ , tortuosity α , viscous Λ and thermal lengths Λ , and viscous Σ and thermal surfaces Σ ).
In the future, we plan to perform numerical simulations to validate and compare the proposed approaches and consider the problem of reflection and transmission of an acoustic wave by a porous medium by associating the appropriate boundary conditions. Our objective is to propose a model in spherical coordinates that allows us to solve the inverse problem based on experimental data of reflected and/or transmitted waves in order to recover the mechanical and acoustic properties as what currently exists for a unidirectional problem.

Author Contributions

Writing—original draft preparation, Z.E.A.F., M.F., E.O., N.O.O. and C.D.; analysis, analytical calculations, and processing of results, Z.E.A.F., R.R., M.F., E.O. and C.D.; writing—review and editing, Z.E.A.F., M.F., E.O. and C.D.; resources, Z.E.A.F., M.F., E.O. and C.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Calculus of D ( t )   =   L 1 1 k ( s )

Let
d ( s ) = L D ( t ) = c s 2 + b s + c ,
we have:
D ( t ) = L 1 d ( s ) = c exp b 2 t J 0 t 2 4 c b 2 if 4 c b 2 > 0 .
Using the property
L 1 d ( s ) s = 1 π t 0 exp u 2 4 t d u ,
we obtain:
D ( t ) = L 1 1 k ( s ) = L 1 c s 2 + b s + c = c π t 0 exp U 2 4 t exp b U 2 J 0 U 2 4 c b 2 d u ,
where 4 c b 2 < 0 , and we take into account the fact that
J 0 U 2 4 c b 2 = J 0 i U 2 b 2 4 c = I 0 U 2 b 2 4 c .

Appendix B. Calculus of S ( t )   =   L 1 1 k 2 ( s )

We have:
L 1 c s 2 + b s + c = 2 c b 4 c exp b 2 t sinh t 2 b 2 4 c , if b 2 4 c > 0 , 2 c 4 c b 2 exp b 2 t sin 4 c b 2 , if 4 c b 2 > 0 .
Let be E ( s ) the inverse Laplace transform of a function e ( t ) ; we have the property:
L 1 E ( s ) s = 0 Erfc U 2 t e ( u ) d u ,
where Erfc ( x ) is the complementary error function given by:
Erfc ( x ) = 2 π x + exp t 2 d t .
Tacking
E ( s ) = c s 2 + b s + c ,
and we have:
E ( s ) s = 1 k 2 ( s ) = c z 2 + b z z + c z ,
so that
S ( t ) = L 1 1 k 2 ( s ) = 2 c b 4 c 0 Erfc U 2 t exp b 2 U sinh U 2 b 2 4 c d U if b 2 4 c > 0 , 2 c 4 c b 2 0 Erfc U 2 t exp b 2 U sin U 2 4 c b 2 d U if 4 c b 2 > 0 .

References

  1. Allard, J.F. Propagation of Sound in Porous Media: Modeling Sound Absorbing Materials; Chapman and Hall: London, UK, 1993. [Google Scholar]
  2. Fellah, Z.E.A.; Fellah, M.; Motri, F.G.; Sebaa, N.; Depollier, C.; Lauriks, W. Measuring permeability of porous materials at low frequency range via acoustic transmitted waves. Rev. Sci. Instrum. 2007, 78, 114902. [Google Scholar] [CrossRef] [PubMed]
  3. Zhao, D.; Li, X.Y. A review of acoustic dampers applied to combustion chambers in aerospace industry. Prog. Aerosp. Sci. 2015, 74, 114–130. [Google Scholar] [CrossRef]
  4. Ogam, E.; Fellah, Z.E.A.; Sebaa, N.; Groby, J.P. Non-ambiguous recovery of Biot poroelastic parameters of cellular panels using ultrasonicwaves. J. Sound Vib. 2011, 330, 1074–1090. [Google Scholar] [CrossRef]
  5. Garay, R.; Nidia Pino, N.G. Acoustic Behavior in Three Types of Housing: Brick Social Housing, Structural Insulated Panel (SIP) Emergency Housing and Mediagua Emergency Housing. Environ. Sci. Rev. Constr. 2019, 18, 96–110. [Google Scholar] [CrossRef]
  6. Johnson, D.L.; Koplik, J.; Dashen, R. Theory of dynamic permeability and tortuosity in fluid-saturated porous media. J. Fluid Mech. 1987, 176, 379–402. [Google Scholar] [CrossRef]
  7. Lafarge, D.; Lemarinier, P.; Allard, J.-F.; Tarnow, V. Dynamic compressibility of air in porous structures at audible frequencies. J. Acoust. Soc. Am. 1997, 102, 1995–2006. [Google Scholar] [CrossRef] [Green Version]
  8. Brown, N.; Melon, M.; Montembault, V.; Castagnède, B.; Lauriks, W.; Leclaire, P. Evaluation of the viscous characteristic length of air-saturated porous materials from the ultrasonic dispersion curve. Comptes R. Acad. Sci. 1996, 322, 122–127. [Google Scholar]
  9. Leclaire, P.; Kelders, L.; Lauriks, W.; Melon, M.; Brown, N.; Castagnede, B. Determination of the viscous and thermal characteristic lengths of plastic foams by ultrasonic measurements in helium and air. J. Appl. Phys. 1996, 80, 2009–2012. [Google Scholar] [CrossRef] [Green Version]
  10. Roncen, R.; Fellah, Z.E.A.; Simon, F.; Piot, E.; Fellah, M.; Ogam, E.; Depollier, C. Bayesian inference for the ultrasonic characterization of rigid porous materials using reflected waves by the first interface. J. Acoust. Soc. Am. 2018, 144, 210–221. [Google Scholar] [CrossRef] [Green Version]
  11. Fellah, Z.E.A.; Depollier, C.; Fellah, M.; Lauriks, W.; Chapeelon, J.-Y. Influence of dynamic tortuosity and compressibility on the propagation of transient waves in porous media. Wave Motion 2005, 41, 145–161. [Google Scholar] [CrossRef]
  12. Fellah, Z.E.A.; Fellah, M.; Lauriks, W.; Depollier, C.; Chapelon, J.Y.; Angel, Y.C. Solution in time domain of ultrasonic propagation equation in a porous material. Wave Motion 2003, 38, 151–163. [Google Scholar] [CrossRef]
  13. Kergomard, J.; Lafarge, D.; Gilbert, J. Transients in porous media: Exact and modelled time-domain Green’s functions. Acta Acust. Acust. 2013, 99, 557–571. [Google Scholar] [CrossRef] [Green Version]
  14. Roncen, R.; Fellah, Z.E.A.; Piot, E.; Simon, F.; Ogam, E.; Fellah, M.; Depollier, C. Inverse identification of a higher order viscous parameter of rigid porous media in the high frequency domain. J. Acoust. Soc. Am. 2019, 145, 1629–1639. [Google Scholar] [CrossRef] [PubMed]
  15. Fellah, Z.E.A.; Fellah, M.; Mitri, F.G.; Sebaa, N.; Lauriks, W.; Depollier, C. Transient acoustic wave propagation in air-saturated porous media at low frequencies. J. Appl. Phys. 2019, 102, 084906. [Google Scholar] [CrossRef]
  16. Fellah, Z.E.A.; Fellah, M.; Sebaa, N.; Lauriks, W.; Depollier, C. Measuring flow resistivity of porous materials at low frequencies range via acoustic transmitted waves. J. Acoust. Soc. Am. 2006, 119, 1926–1928. [Google Scholar] [CrossRef] [PubMed]
  17. Sadouki, M.; Fellah, M.; Fellah, Z.E.A.; Ogam, E.; Sebaa, N.; Mitri, F.G.; Depollier, C. Measuring static thermal permeability and inertial factor of rigid porous materials. J. Acoust. Soc. Am. 2011, 130, 2627–2630. [Google Scholar] [CrossRef]
  18. Roncen, R.; Fellah, Z.E.; Lafarge, D.; Piot, E.; Simon, F.; Ogam, E.; Fellah, M.; Depollier, C. Acoustical modeling and Bayesian inference for rigid porous media in the low-mid frequency regime. J. Acoust. Soc. Am. 2018, 144, 3084–3101. [Google Scholar] [CrossRef]
  19. Fellah, Z.E.A.; Fellah, M.; Lauriks, W. Direct and inverse scattering of transient acoustic waves by a slab of rigid porous material. J. Acoust. Soc. Am. 2003, 113, 61–72. [Google Scholar] [CrossRef]
  20. Cuenca, J.; Goransson, P.; De Ryck, L.; Lahivaara, T. Deterministic and statistical methods for the characterisation of poroelastic media from multi-observation sound absorption measurements. Mech. Syst. Signal Process. 2022, 163, 108186. [Google Scholar] [CrossRef]
  21. Fellah, Z.E.A.; Sadouki, M.; Fellah, M.; Mitri, F.G.; Ogam, E.; Depollier, C. Simultaneous determination of porosity, tortuosity, viscous and thermal characteristic lengths of rigid porous materials. J. Appl. Phys. 2013, 114, 204902. [Google Scholar] [CrossRef]
  22. Jaouen, L.; Gourdon, E.; Gle, P. Estimation of all six parameters of Johnson-Champoux-Allard-Lafarge model for acoustical porous materials from impedance tube measurements. J. Acoust. Soc. Am. 2020, 148, 1998–2005. [Google Scholar] [CrossRef] [PubMed]
  23. Fellah, Z.E.A.; Berger, S.; Lauriks, W.; Depollier, C.; Chapelon, J.-Y. Inverse problem in air-saturated porous media via reflected waves. Rev. Sci. Instrum. 2003, 74, 2871–2879. [Google Scholar] [CrossRef]
  24. Zielinski, T.G. Normalized inverse characterization of sound absorbing rigid porous media. J. Acoust. Soc. Am. 2015, 137, 3232–3243. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Alruwaili, A.D.; Seadawy, A.R.; Ali, A.; Beinane, S.A.O. Novel Analytical Approach for the Space-Time Fractional (2+1)-Dimensional Breaking Soliton Equation via Mathematical Methods. Mathematics 2021, 9, 3253. [Google Scholar] [CrossRef]
  26. Shokri, A. A new eight-order symmetric two-step multiderivative method for the numerical solution of second-order IVPs with oscillating solutions. Numer. Algorithm 2018, 77, 95–109. [Google Scholar] [CrossRef]
  27. Bockstal, K.V. Existence of a Unique Weak Solution to a Nonlinear Non-Autonomous Time-Fractional Wave Equation (of Distributed-Order). Mathematics 2020, 8, 1283. [Google Scholar] [CrossRef]
  28. De Rosa, S.; Polimeni, A.; De Velli, G.; Conte, M.; Sorrentino, S.; Spaccarotella, C.; Mongiardo, A.; Sabatino, J.; Contarini, M.; Todaro, D.; et al. Reliability of Instantaneous Wave-Free Ratio (iFR) for the Evaluation of Left Main Coronary Artery Lesions. J. Clin. Med. 2019, 8, 1143. [Google Scholar] [CrossRef] [Green Version]
  29. Skudrzyk, E. Solution of the Wave Equation in General Spherical Coordinates. In The Foundations of Acoustics; Springer: Vienna, Austria, 1971; pp. 378–391. [Google Scholar]
  30. Keller, J. Stochastic equations and wave propagation in random media. In Stochastic Processess in Mathematical Physics and Engineering; Richard Ernest Bellman, American Mathematical Society: Providence, RI, USA, 1964; pp. 145–170. [Google Scholar]
  31. Biot, M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. I. Low-frequency range. J. Acoust. Soc. Am. 1956, 28, 168–178. [Google Scholar] [CrossRef]
  32. Biot, M.A. Theory of propagation of elastic waves in a fluid-saturated porous solid. II. Higher frequency range. J. Acoust. Soc. Am. 1956, 28, 179–191. [Google Scholar] [CrossRef]
  33. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1993; pp. 93–112. [Google Scholar]
  34. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products Academic Press; Jeffrey, A., Zwillinger, D., Eds.; Elsevier Academic Press: Amsterdam, The Netherlands, 2007. [Google Scholar]
  35. Abramowitz, M.; Stegun, I.A. Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables; United States Department of Commerce: Washington, DC, USA, 1964.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fellah, Z.E.A.; Fellah, M.; Roncen, R.; Ongwen, N.O.; Ogam, E.; Depollier, C. Transient Propagation of Spherical Waves in Porous Material: Application of Fractional Calculus. Symmetry 2022, 14, 233. https://doi.org/10.3390/sym14020233

AMA Style

Fellah ZEA, Fellah M, Roncen R, Ongwen NO, Ogam E, Depollier C. Transient Propagation of Spherical Waves in Porous Material: Application of Fractional Calculus. Symmetry. 2022; 14(2):233. https://doi.org/10.3390/sym14020233

Chicago/Turabian Style

Fellah, Zine El Abiddine, Mohamed Fellah, Rémi Roncen, Nicholas O. Ongwen, Erick Ogam, and Claude Depollier. 2022. "Transient Propagation of Spherical Waves in Porous Material: Application of Fractional Calculus" Symmetry 14, no. 2: 233. https://doi.org/10.3390/sym14020233

APA Style

Fellah, Z. E. A., Fellah, M., Roncen, R., Ongwen, N. O., Ogam, E., & Depollier, C. (2022). Transient Propagation of Spherical Waves in Porous Material: Application of Fractional Calculus. Symmetry, 14(2), 233. https://doi.org/10.3390/sym14020233

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop