Mathematical modeling of two-component media with a saturated liquid began over 90 years with studies of the consolidation of soils. Two-component must be taken into account when solving a significant number of applied problems arising in various areas of human activity (soils, foams, various cement mortars, sand, porous ceramics, porous sintered composite materials, etc.). Two-component media are widely used in the national economy. For example, in the construction of new airfields and the restoration of destroyed, where the building materials used contain a significant number of voids. The study of wave processes is also very important for the development of new diagnostic methods, new technologies for creating two-component environments that could be applied in the field of engineering, construction, instrumentation, metallurgy, nuclear power and the defense capability of the country. However, the complexity of describing the effects of the interaction of components, heat transfer, and other related processes has led to the fact that until now the generally accepted models (elastic medium-liquid) for a fluid-saturated two-component medium have not been fully developed. Therefore, it is of interest to develop a mathematical two-component model when one of the components represents an inhomogeneous viscoelastic medium and the other is a compressible fluid. The presence and degree of porosity in materials is accounted for by a porosity coefficient equal to the ratio of the pore volume to the total volume occupied by the medium.
The propagation of shock waves in an isotropic viscoelastic homogeneous and inhomogeneous medium was considered in [
It is of interest to consider shock waves in a fluid-saturated inhomogeneous two-component medium when one of the components is viscoelastic. Porosity [
Differential equations are obtained for determining the amplitude of longitudinal and transverse shock waves in a fluid-saturated inhomogeneous two-component viscoelastic medium using the mathematical theory of discontinuities [
The interpenetrating motion of a viscoelastic component and a liquid is considered as a fluid motion in a deformable porous medium.
Rheological relations in such a medium through deformations will be written as
T i j = Λ ∗ e k k ( 1 ) δ i j + 2 M ∗ e i j ( 1 ) + A 1 e k k ( 2 ) δ i j , N = A 1 e k k ( 1 ) + A 2 e k k ( 2 ) (1)
Here Λ ∗ , M ∗ are linear integral operators whose relaxation kernels depend continuously on the spatial coordinates:
Λ ∗ = λ ( 1 + Λ ˜ ) , Λ ˜ e ( 1 ) = ∫ 0 ∞ Λ ( t , x i ) e ( 1 ) ( t − t 1 ) d t 1 M ∗ = μ ( 1 + M ˜ ) , M ˜ e ( 1 ) = ∫ 0 ∞ M ( t , x i ) e ( 1 ) ( t − t 1 ) d t 1 (2)
here λ = λ ( x i ) , μ = μ ( x i ) , A 1 ( x i ) = ( 1 − m ) R 0 , A 2 ( x i ) = m R 0 ―continuous functions of coordinates, m = m ( x i ) porosity, R 0 = R 0 ( x i ) ―modulus of fluid compressibility, Т i j ―of the stress tensor of a porous medium, N―force acting on the liquid per unit area of the cross section of the porous medium, e i j ( 1 ) , e k k ( 2 ) ―respectively, the components of the strain tensors of components, δ―the Kronecker symbol. Index 1 above after the letter refers to the viscoelastic component, 2 to the liquid.
Relations (1), (2) together with the equations of motion
ρ 11 V ˙ i ( 1 ) + ρ 12 V ˙ i ( 2 ) = T i j , j , ρ 12 V ˙ i ( 1 ) + ρ 22 V ˙ i ( 2 ) = N , i ρ 11 = ρ 1 − ρ 12 , ρ 22 = ρ 2 − ρ 12 (3)
and Cauchy formulas
2 e i j ( 1 ) = u i , j ( 1 ) + u j , i ( 1 ) , e k k ( 2 ) = u k , k ( 2 ) (4)
represent a closed system for describing the dynamic behavior of a fluid-saturated inhomogeneous two-component medium.
In formulas (3) and (4), ρ 11 = ρ 11 ( x i ) ―the effective density of the viscoelastic component, the ρ 22 = ρ 22 ( x i ) ―effective density of the fluid, ρ 12 = ρ 12 ( x i ) ―the coefficient of dynamic coupling of the viscoelastic component and the liquid, ρ 1 ( x i ) , ρ 2 ( x i ) ―and the density of the viscoelastic component and liquid, respectively.
The repeated Latin indices assume a summation of one to three, in Greek―from one to two. The dot above the letter indicates the time derivative, and the comma below the letter is the derivative of the coordinate.
A shock wave in an inhomogeneous two-component medium saturated with a liquid is determined by an isolated surface ∑ ( t ) on which the u i ( α ) ( α = 1 , 2 ) displacements are continuous in time and coordinates, and the stresses T i j for ces N, and velocities of phase V i ( α ) ( α = 1 , 2 ) discontinuities undergo. The parameters of the porous medium and their gradients are continuous.
Dynamic relationships must be satisfied on the wave surface ∑ ( t ) of a two-component medium [
[ T i j ] ν j = − ρ 11 с [ V i ( 1 ) ] − ρ 12 с [ V i ( 2 ) ] , [ N ] ν i = − ρ 12 с [ V i ( 1 ) ] − ρ 22 с [ V i ( 2 ) ] (5)
where с ( x i , t ) ―the normal velocity of propagation of the surface, the ν i ―components of the unit normal to the surface, directed into the unperturbed region, which depend on the spatial coordinates.
The sign [ ] denotes the difference in the values of a certain quantity on different sides of the rupture surface.
From the rheological relations (1) and (2) written in the discontinuities, dynamic (5) and kinematic conditions of first-order compatibility for the phases, we obtain the system of equations
{ ( ρ 11 c 2 − μ ) δ i j − ( λ + μ ) ν i ν j } [ V j ( 1 ) ] + ( ρ 12 c 2 δ i j − A 1 ν i ν j ) [ V j ( 2 ) ] = 0 ( ρ 12 с 2 δ i j − A 1 ν i ν j ) [ V j ( 1 ) ] + ( ρ 22 c 2 δ i j − A 2 ν i ν j ) [ V j ( 2 ) ] = 0 (6)
Assuming that [ V j ( 1 ) ] = S j ( 1 ) ν j = S 1 ≠ 0 , [ V j ( 2 ) ] = S j ( 2 ) ν j = S 2 ≠ 0 on the surface ∑ ( t ) where S j ( α ) ( α = 1 , 2 ) the values of the discontinuities of the
displacements of the components, multiply (6) by ν i and sum over the repeating index i, we obtain a homogeneous system of equations for S 1 , S 2 . For system (6) to have a nonzero solution, its determinant must be zero. This condition leads to an equation for determining the propagation velocity of a longitudinal wave ( с = с p )
k c p 4 − k 1 c p 2 + k 2 = 0 , k = ρ 11 ρ 22 − ρ 12 2
k 1 = ρ 11 A 2 + ρ 22 Λ 1 − 2 ρ 12 A 1 , k 2 = Λ 1 A 2 − A 1 2 , Λ 1 = λ + 2 μ
Equation (7), generally speaking, has two distinct roots, each of which determines the propagation velocity of the longitudinal wave c p 1 and c p 2
c p 1 , 2 = k 1 ± k 1 2 − 4 k k 2 2 k (7) (I)
Thus, in an inhomogeneous viscoelastic two-component medium, there are two types of longitudinal waves, which are determined by Formula (I).
On the other hand, if on the surface [ V j ( 1 ) ] ν j = 0 , [ V j ( 2 ) ] ν j = 0 , then from (6) it follows that the propagation velocity of the transverse wave ( c = c t ) is given by
c t = μ ρ 22 k (8) (II)
Thus, it is shown that in a two-component inhomogeneous viscoelastic medium, there are three types of shock waves: two types of longitudinal waves and one transverse, for which, [ V i ( α ) ] ν i = S α ≠ 0 , [ V i ( α ) ] ν i = 0 respectively, and the velocities of these waves are determined by Formulas (7) and (8).
If the connection between the viscoelastic component and the liquid is weak ( ρ 12 = 0 , A 1 = 0 ) , then from (7) and (8)
c p 1 = Λ 1 ρ 11 = λ + 2 μ ρ 1 , c p 2 = A 2 ρ 2 , с t = μ ρ 11 = μ ρ 1 (9)
It follows from (9) that the propagation velocities of waves in a viscoelastic two-component medium are equal to the velocities of waves propagating separately in a continuous viscoelastic medium and liquid.
Let us determine the change in the amplitude of the waves. To do this, we differentiate (1) with respect to u and take their difference on different sides of the discontinuity surface, we obtain
[ T ˙ i j ] = λ [ V k , k ( 1 ) ] δ i j + μ ( [ V i , j ( 1 ) ] + [ V j , i ( 1 ) ] ) − λ Λ ( 0 , x i ) [ u k , k ( 1 ) ] δ i j − μ M ( 0 , x i ) ( [ u i , j ( 1 ) ] + [ u j , i ( 1 ) ] ) + A 1 [ V k , k ( 2 ) ] δ i j , [ N ˙ ] = A 1 [ V k , k ( 1 ) ] + A 2 [ V k , k ( 2 ) ] (10)
Taking into account the compatibility conditions of the first order [
δ [ T i j ] δ t − M i j c = λ ( L k ( 1 ) ν k + g α β [ V k ( 1 ) ] , α x β k ) δ i j + μ { L i ( 1 ) ν j + L j ( 1 ) ν i } + μ g α β ( [ V i ( 1 ) ] α x β j + [ V j ( 1 ) ] , α x β i ) + { λ Λ ( 0 , x i ) [ V k ( 1 ) ] ν k δ i j + μ M ( 0 , x i ) ( [ V i ( 1 ) ] ν j + [ V j ( 1 ) ] ν i ) } c − 1 + A 1 ( L k ( 2 ) ν k + g α β [ V k ( 2 ) ] , α x β k ) δ i j
δ [ N ] δ t − χ c = A 1 ( L k ( 1 ) ν k + g α β [ V k ( 1 ) ] , α x β k ) + A 2 ( L k ( 2 ) ν k + g α β [ V k ( 2 ) ] , α x β k ) − ρ 11 с L i ( 1 ) + ρ 11 δ [ V i ( 1 ) ] δ t − ρ 12 с L i ( 2 ) + ρ 12 δ [ V i ( 2 ) ] δ t = M i j ν i + g α β [ N ] , α x β j − ρ 12 с L i ( 1 ) + ρ 12 δ [ V i ( 1 ) ] δ t − ρ 22 с L i ( 2 ) + ρ 22 δ [ V i ( 2 ) ] δ t = χ ν i + g α β [ N ] , α x β i (11)
The quantities, M i j , L i ( α ) ( α = 1 , 2 ) , χ and are determined on the surface ∑ ( t ) and characterize the jumps of the first derivative of the stresses, the forces acting on the fluid and the velocity of displacement of the components, respectively, g α β ―the coefficients of the first fundamental quadratic form of the surface, x β i ―the derivatives of the Cartesian coordinates x i along the curvilinear coordinates y β and δ / δ t − δ ―the time derivative.
By the standard method, we exclude the quantities M i j , χ from (11). As a result, we get
ρ 11 c δ [ V i ( 1 ) ] δ t + ρ 12 c δ [ V i ( 2 ) ] δ t − ρ 11 c 2 L i ( 1 ) + λ ( L k ( 1 ) ν k ν i + g α β [ V k ( 1 ) ] , α x β k ν i ) + μ { L i ( 1 ) + L k ( 1 ) ν k ν i + g α β [ V k ( 1 ) ] , α x β i ν k } + λ Λ ( 0 , x i ) c − 1 [ V k ( 1 ) ] ν k ν i + μ M ( 0 , x i ) c − 1 ( [ V i ( 1 ) ] + [ V k ( 1 ) ] ν k ν i ) + A 1 ( L k ( 2 ) ν k ν i + g α β [ V k ( 2 ) ] , α x β k ν i ) − g α β [ T i k ] , α c x β k − δ [ T i k ] δ t ν k = 0 (12)
A 1 ( L k ( 1 ) ν k ν i + g α β [ V k ( 1 ) ] , α x β k ν i ) + A 2 ( L k ( 2 ) ν k ν i + g α β [ V k ( 2 ) ] , α x β k ν i ) − ρ 12 c 2 L i ( 1 ) + ρ 12 c δ [ V i ( 1 ) ] δ t − ρ 22 c 2 L i ( 2 ) + ρ 22 c δ [ V i ( 2 ) ] δ t − g α β [ N ] , α c x β i − δ [ N ] δ t ν i = 0 (13)
From the rheological relations (1) written in the discontinuities, the dynamic (5) and kinematic conditions of first-order compatibility, we find, [ T i k ] , α ,
[ N ] , α , δ [ T i k ] δ t ν k , δ [ N ] δ t ν i , and substitute in (12)-(13). After the transformations,
we obtain a system of equations determining the change in the components of
the wave amplitude vector ( [ V i ( α ) ] = S i ( α ) , α = 1 , 2 )
{ ( λ + μ ) ν k ν I + ( μ − ρ 11 c 2 ) δ i k } L k ( 1 ) + ( A 1 ν k ν i − ρ 12 c 2 δ i k ) L k ( 2 ) + 2 ρ 11 c δ S i ( 1 ) δ t + 2 ρ 12 c δ S i ( 2 ) δ t + ρ 11 δ c δ t S i ( 1 ) + ρ 12 δ c δ t S i ( 2 ) + 1 c { ( λ Λ ( 0 , x i ) + μ M ( 0 , x i ) ) ν k ν i + μ M ( 0 , x i ) δ i k } S k ( 1 ) + ( λ + μ ) g α β S k , α ( 1 ) x β i ν k + ( λ + μ ) g α β S k , α ( 1 ) x β k ν i − 2 μ Ω S i ( 1 ) + λ , α g α β x β i S β ( 1 ) ν k + λ g α β x β i S k ( 1 ) ν k , α + g α β x β k ( μ S k ( 1 ) ν i , α + μ , α S k ( 1 ) ν i ) + A 1 g α β S k , α ( 2 ) x β k ν i + g α β x β i ( A 1 , α S k ( 2 ) + A 1 S k , α ( 2 ) ν k + A 1 S k ( 2 ) ν k , α ) = 0 (14)
( A 1 ν k ν i − ρ 12 c 2 δ i k ) L k ( 1 ) + 2 ρ 12 c δ S i ( 1 ) δ t + ρ 12 δ c δ t S i ( 1 ) + A 1 g α β S k , α x β k ν i + ( A 2 ν k ν i − ρ 22 c 2 δ i k ) L k ( 2 ) + 2 ρ 22 c δ S i ( 2 ) δ t + ρ 22 δ c δ t S i ( 2 ) − g α β x β i { ρ 12 c c , α S k ( 1 ) ν k + ρ 22 c c , α S k ( 2 ) ν k − ( A 1 S k ( 1 ) ν k ) , α − ( A 2 S k ( 2 ) ν k ) , α } = 0 (15)
where Ω―the average curvature of the wave surface.
From the relations (14) and (15) we obtain differential equations for determining the change in the amplitude of the longitudinal and transverse waves ( [ V i ( α ) ] ν i = S i ( α ) ν i ≠ 0 , S p ( α ) = S i ( α ) S i ( α ) , [ V i ( α ) ] ν i = 0 , α = 1 , 2 ).
For longitudinal waves, we multiply (14) and (15) by ν i and sum over the repeated indices. As a result, we get
( Λ 1 − ρ 11 с p 2 ) L i ( 1 ) ν i + ( A 1 − ρ 12 c p 2 ) L i ( 2 ) ν i + ρ 11 δ c δ t S p ( 1 ) + 2 ρ 11 c p δ S p ( 1 ) δ t − 2 Λ 1 Ω p S p ( 1 ) + 1 c p ( λ Λ ( 0 , x i ) + 2 μ M ( 0 , x i ) ) S p ( 1 ) + ρ 12 δ c δ t S p ( 2 ) + 2 ρ 12 c p δ S p ( 2 ) δ t − 2 A 1 Ω p S p ( 2 ) = 0 , Λ 1 = λ + 2 μ (16)
( A 1 − ρ 12 c p 2 ) L i ( 1 ) ν i + ( A 2 − ρ 22 c p 2 ) L i ( 2 ) ν i + 2 c p ( ρ 12 δ S p ( 1 ) δ t + ρ 22 δ S p ( 2 ) δ t ) + δ c p δ t ( ρ 12 S p ( 1 ) + ρ 22 S p ( 2 ) ) − 2 Ω p ( A 1 S p ( 1 ) + A 2 S p ( 2 ) ) = 0 (17)
We eliminate the quantity L i ( 2 ) from (16) and (17). To do this, we multiply (16) by, A 2 − ρ 22 c p 2 and (17) by A 1 − ρ 12 c p 2 and subtract, we obtain the following equation of the first component of the two-component medium
δ S p ( 1 ) δ t = { Ω p c p − β 2 2 c p ⋅ δ c p δ y − A 1 − ρ 12 c p 2 2 c p 2 a 1 ( λ Λ ( 0 , x i ) + 2 μ M ( 0 , x i ) ) } S p ( 1 ) a 1 = ( ρ 11 ρ 12 A 2 + ρ 12 ρ 22 Λ 1 − 2 ρ 11 ρ 22 A 1 ) c p 2 − 2 ρ 12 A 2 Λ 1 + ρ 22 A 1 Λ 1 + ρ 11 A 1 A 2 β 1 = ( ρ 11 A 1 − ρ 12 Λ 1 ) ( ρ 12 A 2 − ρ 22 A 1 ) , β 2 = A 1 a 1 + ( 4 β 1 − ρ 12 a 1 ) c p 2 a 1 ( A 1 − ρ 12 c p 2 ) (18)
When Equation (18) was obtained, the first Equation (6) and Formula (7) were used.
For the transverse wave from (14) and (15), taking into account the second Equation (6) and Formula (8), we similarly obtain
δ S t ( 1 ) δ t = ( Ω t c t − 1 2 c t ⋅ δ c t δ t − M ( 0 , x i ) 2 ) S t ( 1 ) (19)
Passing to the variable σ ≥ 0 denoting the distance along the normals to the surface ∑ ( t 0 ) in (18) and (19), we obtain equations for changing the amplitude of the longitudinal and transverse waves during their propagation in the first component
d S l ( 1 ) d σ = { Ω l − χ 1 d ln c l 2 d σ − χ 2 } S ( 1 ) , l = p , t (20) (III)
where χ 1 = β 2 ―for longitudinal waves; χ 1 = 1 for a transverse wave,
χ 2 = A 1 − ρ 12 c l 2 c l 2 a 1 ( λ Λ ( 0 , σ ) + 2 μ M ( 0 , σ ) ) ―for longitudinal waves; χ 2 = M ( 0 , σ )
―for a transverse wave.
The change in the amplitude of the waves in the second component is obtained from (6)
S l ( 2 ) = Γ α S l ( 1 ) ( α = 1 , 2 ; l = p , t ) , Γ 1 = ρ 11 с p 2 − Λ 1 A 1 − ρ 12 c p 2 , Γ 2 = − ρ 12 ρ 22 (21) (IV)
Then the change in the amplitude of the waves during their propagation in a porous medium will be written in the form
S l = S l ( 1 ) + S l ( 2 ) = ( 1 + Γ α ) S l ( 1 ) (22) (V)
Equations (20) contain mean curvature Ω l , which is a function of σ.
The average curvature Ω l ( σ ) is related to the first g α β and second b α β quadratic forms and the Gaussian curvature K l ( σ ) of the wave surface by the equations [
The amplitude S l ( 1 ) ( l = p , t ) level satisfying Equation (20) is found by the method of successive approximations under the initial conditions, S ( 0 ) ( 0 ) = S 0 ( 0 ) , S ( i ) ( 0 ) = 0 ( i = 1 , 2 , ⋯ ) ―approximation. A homogeneous medium corresponds to the zeroth approximation. In the first approximation, the rate of change in the inhomogeneity of the viscoelastic two-component medium along the beam is taken into account, and in the second, across the beam.
A more detailed solution of equations of the form (20) by the method of successive approximations was considered in [
We consider an inhomogeneous viscoelastic two-component medium saturated with a liquid, characterized by parameters and relaxation nuclei (a special case)
Λ 1 ( x ) = 0.7 x 2 , ρ 11 ( x ) = 0 , 6 x , ρ 12 ( x ) = − 0.1 x , ρ 22 ( x ) = 0.6 x , m ( x ) = 0.5 x , R 0 ( x ) = 0.2 x , Λ ( 0 , x ) + 2 M ( 0 , x ) = 0.7 x (23)
At the moment of time t = 0 in the plane xy along the axis OX, the wave front propagates with velocity c p and c t , which is determined from (7) (8) with x i = x .
c p 1 = 0.7 x , c p 2 = 0.45 x
c t = 0.87 x
Since by the assumption of the problem the wave surfaces form a family of parallel planes, then the average curvature Ω l = 0 at all times.
Then from (20) for, we obtain the dependence of the level of the wave amplitude on the velocity, the physico-mechanical characteristics and the relaxation kernels of the viscoelastic two-component medium for the first phase
S p ( 1 ) = S 0 p ( 1 ) exp ( − ∫ 0 x [ β 2 ( x ) d ln c ( x ) d x + ( 1 − 0.5 x ) 0.2 + 0.1 ⋅ c 2 ( x ) 2 a 1 ( x ) c 3 ( x ) ⋅ 0.17 x 3 ] d x ) (24)
where S 0 p ( 1 ) the value of the function S p ( 1 ) for x = 0 , β 2 ( x ) , и a 1 ( x ) and are from (18).
By the formula (22) we obtain
S p = S 0 p ( 1 + Γ 1 ( x ) ) S p ( 1 ) ( x ) , Γ 1 ( x ) = ρ 11 ( x ) c 2 ( x ) − ( λ ( x ) + 2 μ ( x ) ) ( 1 − m ( x ) ) R 0 ( x ) − ρ 12 ( x ) c 2 ( x ) (25)
Expression (25) determines the change in the level of the wave amplitude in a fluid-saturated inhomogeneous viscoelastic two-component medium.
If we specify in (25) a specific form of the physico-mechanical characteristics of the medium and the relaxation nuclei, we obtain an expression for determining the change in the level of the wave velositi in a fluid-saturated inhomogeneous viscoelastic two-component medium.
Applying the mathematical theory of discontinuities to the basic equations of a two-component viscoelastic porous medium, formulas were first obtained for determining the velocities and changing the level of the amplitude of longitudinal and transverse waves.
In the absence of dynamic coupling between the propagation velocity components of the first and second types of longitudinal waves, the wave velocities propagating separately in a continuous viscoelastic component and a continuous liquid component are equal.
In the medium under consideration, in the absence of a connection between the viscoelastic component and the liquid, the transverse wave propagates only in a continuous viscoelastic component.
Using geometric and dynamic conditions of second-order compatibility of the theory of discontinuities, equations are obtained for determining shock wave amplitudes in a saturated inhomogeneous viscoelastic two-component medium.
The amplitude of the waves in an inhomogeneous viscoelastic two-component medium depends on the porosity of the medium, the viscosity, the initial mean curvature, and the Gaussian curvature of the wave surface.
Polenov, V. and Chigarev, A. (2018) Mathematical Modeling of Shock Waves in Inhomogeneous Viscoelastic Two-Component Media. Journal of Applied Mathematics and Physics, 6, 997-1005. https://doi.org/10.4236/jamp.2018.65086
Т i j : stress tensor of a porous medium
N: force acting on the liquid per unit area of the cross section of the porous medium
e i j ( 1 ) : components of the strain tensors of the elastic component
e k k ( 2 ) : components of fluid deformations
u i ( 1 ) : displacement of the elastic component
u 1 ( 2 ) : fluid movement
Λ ∗ , M ∗ : linear elastic integral operators
S P ( α ) : wave amplitudes
Ω : average curvature of the wave surface
g α β : the coefficients of the first fundamental quadratic form of the surface