Mathematical Modeling of Shock Waves in Inhomogeneous Viscoelastic Two-Component Media

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.


Introduction
The propagation of shock waves in an isotropic viscoelastic homogeneous and inhomogeneous medium was considered in [1] [2] [3].The nature of the damping of waves in such media is determined by the geometry of the discontinuity surface and the rheological model described by the corresponding relaxation nuclei.
It is of interest to consider shock waves in a fluid-saturated inhomogeneous two-component medium when one of the components is viscoelastic.Porosity [4] [5] [6] is understood to be effective porosity, taking into account only interconnected pores.Isolated pores are considered as elements of the viscoelastic part of the porous skeleton.It is assumed that the pore sizes are small in comparison with the distance at which the kinematic and geometric characteristics of the motion change.This allows us to assume that the viscoelastic and liquid components of the medium are continuous and at each point of space in this case there will be two displacement vectors, the displacement vector of the viscoelastic component (the skeleton of the porous medium) and the displacement vector of the liquid.We shall assume that the fluid is compressible.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 [7] [8].

Formulation Problem
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 Here * Λ , M * are linear integral operators whose relaxation kernels depend continuously on the spatial coordinates: here Т -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, ( ) , , and Cauchy formulas ( ) -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.

Velocities of Shock Waves
A shock wave in an inhomogeneous two-component medium saturated with a liquid is determined by an isolated surface displacements are continuous in time and coordinates, and the stresses ij T for ces N, and velocities of phase ( ) ( ) The parameters of the porous medium and their gradients are continuous.
Dynamic relationships must be satisfied on the wave surface where ( ) , i с x 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 Assuming that ( ) ( ) the values of the discontinuities of the displacements of the components, multiply (6) by i ν and sum over the repeat- ing index i, we obtain a homogeneous system of equations for 1 S , 2 S .For sys- tem (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 ( ) 7), generally speaking, has two distinct roots, each of which determines the propagation velocity of the longitudinal wave 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 , then from (6) it follows that the propagation velocity of the transverse wave ( ) given by 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, ( ) 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 ( ) , then from ( 7) and ( 8) 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.

Determination of Shock Wave Amplitudes
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 Taking into account the compatibility conditions of the first order [6] and the equations of motion written in the discontinuities, (10) take the following form [ ] The quantities, ij M , ( ) ( ) , χ 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, i x β -the derivatives of the Cartesian coordinates i x along the cur- vilinear coordinates y β and t δ δ δ − -the time derivative.
By the standard method, we exclude the quantities ij M , χ from (11).As a re- sult, we get From the rheological relations (1) written in the discontinuities, the dynamic (5) and kinematic conditions of first-order compatibility, we find, [ ] , 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 { } 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 ( ( ) ( ) For longitudinal waves, we multiply ( 14) and ( 15) by i ν and sum over the repeated indices.As a result, we We eliminate the quantity ( ) 16) and (17).To do this, we multiply (16) by, When Equation (18) was obtained, the first Equation ( 6) and Formula (7) were used.

Conclusions
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.
system for describing the dynamic behavior of a fluid-saturated inhomogeneous two-component medium.In formulas (3) and (4), of dynamic coupling of the viscoelastic component and the liquid, 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.
-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.Journal of Applied Mathematics and Physics