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

Journal of Applied Mathematics and Physics
Vol.06 No.05(2018), Article ID:84679,9 pages
10.4236/jamp.2018.65086

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

Victor Polenov1, Anatoli Chigarev2

1Department of Mathematics, Voronezh State University, Voronezh, Russia

2Department of Theoretical Mechanics, Belarusian National Technical University, Minsk, Belarus    Received: March 30, 2018; Accepted: May 19, 2018; Published: May 22, 2018

ABSTRACT

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.

Keywords:

Wave Propagation, Porous Medium, Dynamic Quantities, Elasticity, Viscoelastic, Fluid 1. Introduction

The propagation of shock waves in an isotropic viscoelastic homogeneous and inhomogeneous medium was considered in    . 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    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   .

2. 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

${T}_{ij}={\Lambda }^{\ast }{e}_{kk}^{\left(1\right)}{\delta }_{ij}+2{M}^{\ast }{e}_{ij}^{\left(1\right)}+{A}_{1}{e}_{kk}^{\left(2\right)}{\delta }_{ij}$ , $N={A}_{1}{e}_{kk}^{\left(1\right)}+{A}_{2}{e}_{kk}^{\left(2\right)}$ (1)

Here ${\Lambda }^{\ast }$ , ${M}^{\ast }$ are linear integral operators whose relaxation kernels depend continuously on the spatial coordinates:

$\begin{array}{l}{\Lambda }^{\ast }=\lambda \left(1+\stackrel{˜}{\Lambda }\right),\text{\hspace{0.17em}}\stackrel{˜}{\Lambda }{e}^{\left(1\right)}=\underset{0}{\overset{\infty }{\int }}\Lambda \left(t,{x}_{i}\right){e}^{\left(1\right)}\left(t-{t}_{1}\right)\text{d}{t}_{1}\\ {M}^{\ast }=\mu \left(1+\stackrel{˜}{M}\right),\text{\hspace{0.17em}}\stackrel{˜}{M}{e}^{\left(1\right)}=\underset{0}{\overset{\infty }{\int }}M\left(t,{x}_{i}\right){e}^{\left(1\right)}\left(t-{t}_{1}\right)\text{d}{t}_{1}\end{array}$ (2)

here $\lambda =\lambda \left({x}_{i}\right)$ , $\mu =\mu \left({x}_{i}\right)$ , ${A}_{1}\left({x}_{i}\right)=\left(1-m\right){R}_{0}$ , ${A}_{2}\left({x}_{i}\right)=m{R}_{0}$ ―continuous functions of coordinates, $m=m\left({x}_{i}\right)$ porosity, ${R}_{0}={R}_{0}\left({x}_{i}\right)$ ―modulus of fluid compressibility, ${Т}_{ij}$ ―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}_{ij}^{\left(1\right)}$ , ${e}_{kk}^{\left(2\right)}$ ―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

$\begin{array}{l}{\rho }_{11}{\stackrel{˙}{V}}_{i}^{\left(1\right)}+{\rho }_{12}{\stackrel{˙}{V}}_{i}^{\left(2\right)}={T}_{ij,j},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\rho }_{12}{\stackrel{˙}{V}}_{i}^{\left(1\right)}+{\rho }_{22}{\stackrel{˙}{V}}_{i}^{\left(2\right)}={N}_{,i}\\ {\rho }_{11}={\rho }_{1}-{\rho }_{12},\text{\hspace{0.17em}}\text{\hspace{0.17em}}{\rho }_{22}={\rho }_{2}-{\rho }_{12}\end{array}$ (3)

and Cauchy formulas

$2{e}_{ij}^{\left(1\right)}={u}_{i,j}^{\left(1\right)}+{u}_{j,i}^{\left(1\right)}$ , ${e}_{kk}^{\left(2\right)}={u}_{k,k}^{\left(2\right)}$ (4)

represent a closed system for describing the dynamic behavior of a fluid-saturated inhomogeneous two-component medium.

In formulas (3) and (4), ${\rho }_{11}={\rho }_{11}\left({x}_{i}\right)$ ―the effective density of the viscoelastic component, the ${\rho }_{22}={\rho }_{22}\left({x}_{i}\right)$ ―effective density of the fluid, ${\rho }_{12}={\rho }_{12}\left({x}_{i}\right)$ ―the coefficient of dynamic coupling of the viscoelastic component and the liquid, ${\rho }_{1}\left({x}_{i}\right)$ , ${\rho }_{2}\left({x}_{i}\right)$ ―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.

3. Velocities of Shock Waves

A shock wave in an inhomogeneous two-component medium saturated with a liquid is determined by an isolated surface $\sum \left(t\right)$ on which the ${u}_{i}^{\left(\alpha \right)}\left(\alpha =1,2\right)$ displacements are continuous in time and coordinates, and the stresses ${T}_{ij}$ for ces N, and velocities of phase ${V}_{i}^{\left(\alpha \right)}\left(\alpha =1,2\right)$ discontinuities undergo. The parameters of the porous medium and their gradients are continuous.

Dynamic relationships must be satisfied on the wave surface $\sum \left(t\right)$ of a two-component medium 

$\left[{T}_{ij}\right]{\nu }_{j}=-{\rho }_{11}с\left[{V}_{i}^{\left(1\right)}\right]-{\rho }_{12}с\left[{V}_{i}^{\left(2\right)}\right],\text{\hspace{0.17em}}\left[N\right]{\nu }_{i}=-{\rho }_{12}с\left[{V}_{i}^{\left(1\right)}\right]-{\rho }_{22}с\left[{V}_{i}^{\left(2\right)}\right]$ (5)

where $с\left({x}_{i},t\right)$ ―the normal velocity of propagation of the surface, the ${\nu }_{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

$\begin{array}{l}\left\{\left({\rho }_{11}{c}^{2}-\mu \right){\delta }_{ij}-\left(\lambda +\mu \right){\nu }_{i}{\nu }_{j}\right\}\left[{V}_{j}^{\left(1\right)}\right]+\left({\rho }_{12}{c}^{2}{\delta }_{ij}-{A}_{1}{\nu }_{i}{\nu }_{j}\right)\left[{V}_{j}^{\left(2\right)}\right]=0\\ \left({\rho }_{12}{с}^{2}{\delta }_{ij}-{A}_{1}{\nu }_{i}{\nu }_{j}\right)\left[{V}_{j}^{\left(1\right)}\right]+\left({\rho }_{22}{c}^{2}{\delta }_{ij}-{A}_{2}{\nu }_{i}{\nu }_{j}\right)\left[{V}_{j}^{\left(2\right)}\right]=0\end{array}$ (6)

Assuming that $\left[{V}_{j}^{\left(1\right)}\right]={S}_{j}^{\left(1\right)}{\nu }_{j}={S}_{1}\ne 0$ , $\left[{V}_{j}^{\left(2\right)}\right]={S}_{j}^{\left(2\right)}{\nu }_{j}={S}_{2}\ne 0$ on the surface $\sum \left(t\right)$ where ${S}_{j}^{\left(\alpha \right)}\left(\alpha =1,2\right)$ the values of the discontinuities of the

displacements of the components, multiply (6) by ${\nu }_{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 $\left(с={с}_{p}\right)$

$k{c}_{p}^{4}-{k}_{1}{c}_{p}^{2}+{k}_{2}=0$ , $k={\rho }_{11}{\rho }_{22}-{\rho }_{12}^{2}$

${k}_{1}={\rho }_{11}{A}_{2}+{\rho }_{22}{\Lambda }_{1}-2{\rho }_{12}{A}_{1}$ , ${k}_{2}={\Lambda }_{1}{A}_{2}-{A}_{1}^{2}$ , ${\Lambda }_{1}=\lambda +2\mu$

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}}=\sqrt{\frac{{k}_{1}±\sqrt{{k}_{1}^{2}-4k{k}_{2}}}{2k}}$ (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 $\left[{V}_{j}^{\left(1\right)}\right]{\nu }_{j}=0$ , $\left[{V}_{j}^{\left(2\right)}\right]{\nu }_{j}=0$ , then from (6) it follows that the propagation velocity of the transverse wave $\left(c={c}_{t}\right)$ is given by

${c}_{t}=\sqrt{\frac{\mu {\rho }_{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, $\left[{V}_{i}^{\left(\alpha \right)}\right]{\nu }_{i}={S}_{\alpha }\ne 0$ , $\left[{V}_{i}^{\left(\alpha \right)}\right]{\nu }_{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 $\left({\rho }_{12}=0,{A}_{1}=0\right)$ , then from (7) and (8)

${c}_{{p}_{1}}=\sqrt{\frac{{\Lambda }_{1}}{{\rho }_{11}}}=\sqrt{\frac{\lambda +2\mu }{{\rho }_{1}}}$ , ${c}_{{p}_{2}}=\sqrt{\frac{{A}_{2}}{{\rho }_{2}}}$ , ${с}_{t}=\sqrt{\frac{\mu }{{\rho }_{11}}}=\sqrt{\frac{\mu }{{\rho }_{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.

4. 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

$\begin{array}{l}\left[{\stackrel{˙}{T}}_{ij}\right]=\lambda \left[{V}_{k,k}^{\left(1\right)}\right]{\delta }_{ij}+\mu \left(\left[{V}_{i,j}^{\left(1\right)}\right]+\left[{V}_{j,i}^{\left(1\right)}\right]\right)-\lambda \Lambda \left(0,{x}_{i}\right)\left[{u}_{k,k}^{\left(1\right)}\right]{\delta }_{ij}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}-\mu M\left(0,{x}_{i}\right)\left(\left[{u}_{i,j}^{\left(1\right)}\right]+\left[{u}_{j,i}^{\left(1\right)}\right]\right)+{A}_{1}\left[{V}_{k,k}^{\left(2\right)}\right]{\delta }_{ij},\\ \left[\stackrel{˙}{N}\right]={A}_{1}\left[{V}_{k,k}^{\left(1\right)}\right]+{A}_{2}\left[{V}_{k,k}^{\left(2\right)}\right]\end{array}$ (10)

Taking into account the compatibility conditions of the first order  and the equations of motion written in the discontinuities, (10) take the following form

$\begin{array}{c}\frac{\delta \left[{T}_{ij}\right]}{\delta t}-{M}_{ij}c=\lambda \left({L}_{k}^{\left(1\right)}{\nu }_{k}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(1\right)}\right]}_{,\alpha }{x}_{\beta }^{k}\right){\delta }_{ij}+\mu \left\{{L}_{i}^{\left(1\right)}{\nu }_{j}+{L}_{j}^{\left(1\right)}{\nu }_{i}\right\}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\mu {g}^{\alpha \beta }\left({\left[{V}_{i}^{\left(1\right)}\right]}_{\alpha }{x}_{\beta }^{j}+{\left[{V}_{j}^{\left(1\right)}\right]}_{,\alpha }{x}_{\beta }^{i}\right)+\left\{\lambda \Lambda \left(0,{x}_{i}\right)\left[{V}_{k}^{\left(1\right)}\right]{\nu }_{k}{\delta }_{ij}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+\mu M\left(0,{x}_{i}\right)\left(\left[{V}_{i}^{\left(1\right)}\right]{\nu }_{j}+\left[{V}_{j}^{\left(1\right)}\right]{\nu }_{i}\right)\right\}{c}^{-1}\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{A}_{1}\left({L}_{k}^{\left(2\right)}{\nu }_{k}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(2\right)}\right]}_{,\alpha }{x}_{\beta }^{k}\right){\delta }_{ij}\end{array}$

$\begin{array}{l}\frac{\delta \left[N\right]}{\delta t}-\chi c={A}_{1}\left({L}_{k}^{\left(1\right)}{\nu }_{k}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(1\right)}\right]}_{,\alpha }{x}_{\beta }^{k}\right)+{A}_{2}\left({L}_{k}^{\left(2\right)}{\nu }_{k}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(2\right)}\right]}_{,\alpha }{x}_{\beta }^{k}\right)\\ -{\rho }_{11}с{L}_{i}^{\left(1\right)}+{\rho }_{11}\frac{\delta \left[{V}_{i}^{\left(1\right)}\right]}{\delta t}-{\rho }_{12}с{L}_{i}^{\left(2\right)}+{\rho }_{12}\frac{\delta \left[{V}_{i}^{\left(2\right)}\right]}{\delta t}={M}_{ij}{\nu }_{i}+{g}^{\alpha \beta }{\left[N\right]}_{,\alpha }{x}_{\beta }^{j}\\ -{\rho }_{12}с{L}_{i}^{\left(1\right)}+{\rho }_{12}\frac{\delta \left[{V}_{i}^{\left(1\right)}\right]}{\delta t}-{\rho }_{22}с{L}_{i}^{\left(2\right)}+{\rho }_{22}\frac{\delta \left[{V}_{i}^{\left(2\right)}\right]}{\delta t}=\chi {\nu }_{i}+{g}^{\alpha \beta }{\left[N\right]}_{,\alpha }{x}_{\beta }^{i}\end{array}$ (11)

The quantities, ${M}_{ij}$ , ${L}_{i}^{\left(\alpha \right)}\left(\alpha =1,2\right)$ , χ and are determined on the surface $\sum \left(t\right)$ 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}^{\alpha \beta }$ ―the coefficients of the first fundamental quadratic form of the surface, ${x}_{\beta }^{i}$ ―the derivatives of the Cartesian coordinates ${x}_{i}$ along the curvilinear coordinates ${y}_{\beta }$ and $\delta /\delta t-\delta$ ―the time derivative.

By the standard method, we exclude the quantities ${M}_{ij}$ , χ from (11). As a result, we get

$\begin{array}{l}{\rho }_{11}c\frac{\delta \left[{V}_{i}^{\left(1\right)}\right]}{\delta t}+{\rho }_{12}c\frac{\delta \left[{V}_{i}^{\left(2\right)}\right]}{\delta t}-{\rho }_{11}{c}^{2}{L}_{i}^{\left(1\right)}+\lambda \left({L}_{k}^{\left(1\right)}{\nu }_{k}{\nu }_{i}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(1\right)}\right]}_{,\alpha }{x}_{\beta }^{k}{\nu }_{i}\right)\\ +\mu \left\{{L}_{i}^{\left(1\right)}+{L}_{k}^{\left(1\right)}{\nu }_{k}{\nu }_{i}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(1\right)}\right]}_{,\alpha }{x}_{\beta }^{i}{\nu }_{k}\right\}+\lambda \Lambda \left(0,{x}_{i}\right){c}^{-1}\left[{V}_{k}^{\left(1\right)}\right]{\nu }_{k}{\nu }_{i}\\ +\mu M\left(0,{x}_{i}\right){c}^{-1}\left(\left[{V}_{i}^{\left(1\right)}\right]+\left[{V}_{k}^{\left(1\right)}\right]{\nu }_{k}{\nu }_{i}\right)+{A}_{1}\left({L}_{k}^{\left(2\right)}{\nu }_{k}{\nu }_{i}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(2\right)}\right]}_{,\alpha }{x}_{\beta }^{k}{\nu }_{i}\right)\\ -{g}^{\alpha \beta }{\left[{T}_{ik}\right]}_{,\alpha }c{x}_{\beta }^{k}-\frac{\delta \left[{T}_{ik}\right]}{\delta t}{\nu }_{k}=0\end{array}$ (12)

$\begin{array}{l}{A}_{1}\left({L}_{k}^{\left(1\right)}{\nu }_{k}{\nu }_{i}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(1\right)}\right]}_{,\alpha }{x}_{\beta }^{k}{\nu }_{i}\right)+{A}_{2}\left({L}_{k}^{\left(2\right)}{\nu }_{k}{\nu }_{i}+{g}^{\alpha \beta }{\left[{V}_{k}^{\left(2\right)}\right]}_{,\alpha }{x}_{\beta }^{k}{\nu }_{i}\right)\\ \text{ }-{\rho }_{12}{c}^{2}{L}_{i}^{\left(1\right)}+{\rho }_{12}c\frac{\delta \left[{V}_{i}^{\left(1\right)}\right]}{\delta t}-{\rho }_{22}{c}^{2}{L}_{i}^{\left(2\right)}+{\rho }_{22}c\frac{\delta \left[{V}_{i}^{\left(2\right)}\right]}{\delta t}\\ \text{ }-{g}^{\alpha \beta }{\left[N\right]}_{,\alpha }c{x}_{\beta }^{i}-\frac{\delta \left[N\right]}{\delta t}{\nu }_{i}=0\end{array}$ (13)

From the rheological relations (1) written in the discontinuities, the dynamic (5) and kinematic conditions of first-order compatibility, we find, ${\left[{T}_{ik}\right]}_{,\alpha }$ ,

${\left[N\right]}_{,\alpha }$ , $\frac{\delta \left[{T}_{ik}\right]}{\delta t}{\nu }_{k}$ , $\frac{\delta \left[N\right]}{\delta t}{\nu }_{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 $\left(\left[{V}_{i}^{\left(\alpha \right)}\right]={S}_{i}^{\left(\alpha \right)},\alpha =1,2\right)$

$\begin{array}{l}\left\{\left(\lambda +\mu \right){\nu }_{k}{\nu }_{I}+\left(\mu -{\rho }_{11}{c}^{2}\right){\delta }_{ik}\right\}{L}_{k}^{\left(1\right)}+\left({A}_{1}{\nu }_{k}{\nu }_{i}-{\rho }_{12}{c}^{2}{\delta }_{ik}\right){L}_{k}^{\left(2\right)}+2{\rho }_{11}c\frac{\delta {S}_{i}^{\left(1\right)}}{\delta t}\\ \text{ }+2{\rho }_{12}c\frac{\delta {S}_{i}^{\left(2\right)}}{\delta t}+{\rho }_{11}\frac{\delta c}{\delta t}{S}_{i}^{\left(1\right)}+{\rho }_{12}\frac{\delta c}{\delta t}{S}_{i}^{\left(2\right)}+\frac{1}{c}\left\{\left(\lambda \Lambda \left(0,{x}_{i}\right)+\mu M\left(0,{x}_{i}\right)\right){\nu }_{k}{\nu }_{i}\\ \text{ }+\mu M\left(0,{x}_{i}\right){\delta }_{ik}\right\}{S}_{k}^{\left(1\right)}+\left(\lambda +\mu \right){g}^{\alpha \beta }{S}_{k,\alpha }^{\left(1\right)}{x}_{\beta }^{i}{\nu }_{k}+\left(\lambda +\mu \right){g}^{\alpha \beta }{S}_{k,\alpha }^{\left(1\right)}{x}_{\beta }^{k}{\nu }_{i}-2\mu \text{ }\Omega {S}_{i}^{\left(1\right)}\\ \text{ }+{\lambda }_{,\alpha }{g}^{\alpha \beta }{x}_{\beta }^{i}{S}_{\beta }^{\left(1\right)}{\nu }_{k}+\lambda {g}^{\alpha \beta }{x}_{\beta }^{i}{S}_{k}^{\left(1\right)}{\nu }_{k,\alpha }+{g}^{\alpha \beta }{x}_{\beta }^{k}\left(\mu {S}_{k}^{\left(1\right)}{\nu }_{i,\alpha }+{\mu }_{,\alpha }{S}_{k}^{\left(1\right)}{\nu }_{i}\right)\\ \text{ }+{A}_{1}{g}^{\alpha \beta }{S}_{k,\alpha }^{\left(2\right)}{x}_{\beta }^{k}{\nu }_{i}+{g}^{\alpha \beta }{x}_{\beta }^{i}\left({A}_{1,\alpha }{S}_{k}^{\left(2\right)}+{A}_{1}{S}_{k,\alpha }^{\left(2\right)}{\nu }_{k}+{A}_{1}{S}_{k}^{\left(2\right)}{\nu }_{k,\alpha }\right)=0\end{array}$ (14)

$\begin{array}{l}\left({A}_{1}{\nu }_{k}{\nu }_{i}-{\rho }_{12}{c}^{2}{\delta }_{ik}\right){L}_{k}^{\left(1\right)}+2{\rho }_{12}c\frac{\delta {S}_{i}^{\left(1\right)}}{\delta t}+{\rho }_{12}\frac{\delta c}{\delta t}{S}_{i}^{\left(1\right)}+{A}_{1}{g}^{\alpha \beta }{S}_{k,\alpha }{x}_{\beta }^{k}{\nu }_{i}\\ +\left({A}_{2}{\nu }_{k}{\nu }_{i}-{\rho }_{22}{c}^{2}{\delta }_{ik}\right){L}_{k}^{\left(2\right)}+2{\rho }_{22}c\frac{\delta {S}_{i}^{\left(2\right)}}{\delta t}+{\rho }_{22}\frac{\delta c}{\delta t}{S}_{i}^{\left(2\right)}\\ \text{ }-{g}^{\alpha \beta }{x}_{\beta }^{i}\left\{{\rho }_{12}c{c}_{,\alpha }{S}_{k}^{\left(1\right)}{\nu }_{k}+{\rho }_{22}c{c}_{,\alpha }{S}_{k}^{\left(2\right)}{\nu }_{k}-{\left({A}_{1}{S}_{k}^{\left(1\right)}{\nu }_{k}\right)}_{,\alpha }-{\left({A}_{2}{S}_{k}^{\left(2\right)}{\nu }_{k}\right)}_{,\alpha }\right\}=0\end{array}$ (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 ( $\left[{V}_{i}^{\left(\alpha \right)}\right]{\nu }_{i}={S}_{i}^{\left(\alpha \right)}{\nu }_{i}\ne 0$ , ${S}_{p}^{\left(\alpha \right)}=\sqrt{{S}_{i}^{\left(\alpha \right)}{S}_{i}^{\left(\alpha \right)}}$ , $\left[{V}_{i}^{\left(\alpha \right)}\right]{\nu }_{i}=0$ , $\alpha =1,2$ ).

For longitudinal waves, we multiply (14) and (15) by ${\nu }_{i}$ and sum over the repeated indices. As a result, we get

$\begin{array}{l}\left({\Lambda }_{1}-{\rho }_{11}{с}_{p}^{2}\right){L}_{i}^{\left(1\right)}{\nu }_{i}+\left({A}_{1}-{\rho }_{12}{c}_{p}^{2}\right){L}_{i}^{\left(2\right)}{\nu }_{i}+{\rho }_{11}\frac{\delta c}{\delta t}{S}_{p}^{\left(1\right)}+2{\rho }_{11}{c}_{p}\frac{\delta {S}_{p}^{\left(1\right)}}{\delta t}\\ \text{ }-2{\Lambda }_{1}{\Omega }_{p}{S}_{p}^{\left(1\right)}+\frac{1}{{c}_{p}}\left(\lambda \Lambda \left(0,{x}_{i}\right)+2\mu M\left(0,{x}_{i}\right)\right){S}_{p}^{\left(1\right)}+{\rho }_{12}\frac{\delta c}{\delta t}{S}_{p}^{\left(2\right)}\\ \text{ }+2{\rho }_{12}{c}_{p}\frac{\delta {S}_{p}^{\left(2\right)}}{\delta t}-2{A}_{1}{\Omega }_{p}{S}_{p}^{\left(2\right)}=0,\text{\hspace{0.17em}}{\Lambda }_{1}=\lambda +2\mu \end{array}$ (16)

$\begin{array}{l}\left({A}_{1}-{\rho }_{12}{c}_{p}^{2}\right){L}_{i}^{\left(1\right)}{\nu }_{i}+\left({A}_{2}-{\rho }_{22}{c}_{p}^{2}\right){L}_{i}^{\left(2\right)}{\nu }_{i}+2{c}_{p}\left({\rho }_{12}\frac{\delta {S}_{p}^{\left(1\right)}}{\delta t}+{\rho }_{22}\frac{\delta {S}_{p}^{\left(2\right)}}{\delta t}\right)\\ +\frac{\delta {c}_{p}}{\delta t}\left({\rho }_{12}{S}_{p}^{\left(1\right)}+{\rho }_{22}{S}_{p}^{\left(2\right)}\right)-2{\Omega }_{p}\left({A}_{1}{S}_{p}^{\left(1\right)}+{A}_{2}{S}_{p}^{\left(2\right)}\right)=0\end{array}$ (17)

We eliminate the quantity ${L}_{i}^{\left(2\right)}$ from (16) and (17). To do this, we multiply (16) by, ${A}_{2}-{\rho }_{22}{c}_{p}^{2}$ and (17) by ${A}_{1}-{\rho }_{12}{c}_{p}^{2}$ and subtract, we obtain the following equation of the first component of the two-component medium

$\begin{array}{l}\frac{\delta {S}_{p}^{\left(1\right)}}{\delta t}=\left\{{\Omega }_{p}{c}_{p}-\frac{{\beta }_{2}}{2{c}_{p}}\cdot \frac{\delta {c}_{p}}{\delta y}-\frac{{A}_{1}-{\rho }_{12}{c}_{p}^{2}}{2{c}_{p}^{2}{a}_{1}}\left(\lambda \Lambda \left(0,{x}_{i}\right)+2\mu M\left(0,{x}_{i}\right)\right)\right\}{S}_{p}^{\left(1\right)}\\ {a}_{1}=\left({\rho }_{11}{\rho }_{12}{A}_{2}+{\rho }_{12}{\rho }_{22}{\Lambda }_{1}-2{\rho }_{11}{\rho }_{22}{A}_{1}\right){c}_{p}^{2}-2{\rho }_{12}{A}_{2}{\Lambda }_{1}+{\rho }_{22}{A}_{1}{\Lambda }_{1}+{\rho }_{11}{A}_{1}{A}_{2}\\ {\beta }_{1}=\left({\rho }_{11}{A}_{1}-{\rho }_{12}{\Lambda }_{1}\right)\left({\rho }_{12}{A}_{2}-{\rho }_{22}{A}_{1}\right),\text{\hspace{0.17em}}{\beta }_{2}=\frac{{A}_{1}{a}_{1}+\left(4{\beta }_{1}-{\rho }_{12}{a}_{1}\right){c}_{p}^{2}}{{a}_{1}\left({A}_{1}-{\rho }_{12}{c}_{p}^{2}\right)}\end{array}$ (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

$\frac{\delta {S}_{t}^{\left(1\right)}}{\delta t}=\left({\Omega }_{t}{c}_{t}-\frac{1}{2{c}_{t}}\cdot \frac{\delta {c}_{t}}{\delta t}-\frac{M\left(0,{x}_{i}\right)}{2}\right){S}_{t}^{\left(1\right)}$ (19)

Passing to the variable $\sigma \ge 0$ denoting the distance along the normals to the surface $\sum \left({t}_{0}\right)$ in (18) and (19), we obtain equations for changing the amplitude of the longitudinal and transverse waves during their propagation in the first component

$\frac{\text{d}{S}_{l}^{\left(1\right)}}{\text{d}\sigma }=\left\{{\Omega }_{l}-{\chi }_{1}\frac{\text{d}\mathrm{ln}{c}_{l}}{2\text{d}\sigma }-{\chi }_{2}\right\}{S}^{\left(1\right)}$ , $l=p,t$ (20) (III)

where ${\chi }_{1}={\beta }_{2}$ ―for longitudinal waves; ${\chi }_{1}=1$ for a transverse wave,

${\chi }_{2}=\frac{{A}_{1}-{\rho }_{12}{c}_{l}^{2}}{{c}_{l}^{2}{a}_{1}}\left(\lambda \Lambda \left(0,\sigma \right)+2\mu M\left(0,\sigma \right)\right)$ ―for longitudinal waves; ${\chi }_{2}=M\left(0,\sigma \right)$

―for a transverse wave.

The change in the amplitude of the waves in the second component is obtained from (6)

${S}_{l}^{\left(2\right)}={\Gamma }_{\alpha }{S}_{l}^{\left(1\right)}\left(\alpha =1,2;l=p,t\right)$ , ${\Gamma }_{1}=\frac{{\rho }_{11}{с}_{p}^{2}-{\Lambda }_{1}}{{A}_{1}-{\rho }_{12}{c}_{p}^{2}}$ , ${\Gamma }_{2}=-\frac{{\rho }_{12}}{{\rho }_{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}^{\left(1\right)}+{S}_{l}^{\left(2\right)}=\left(1+{\Gamma }_{\alpha }\right){S}_{l}^{\left(1\right)}$ (22) (V)

Equations (20) contain mean curvature ${\Omega }_{l}$ , which is a function of σ.

The average curvature ${\Omega }_{l}\left(\sigma \right)$ is related to the first ${g}^{\alpha \beta }$ and second ${b}^{\alpha \beta }$ quadratic forms and the Gaussian curvature ${K}_{l}\left(\sigma \right)$ of the wave surface by the equations  , where its method of finding is given.

The amplitude ${S}_{l}^{\left(1\right)}\left(l=p,t\right)$ level satisfying Equation (20) is found by the method of successive approximations under the initial conditions, ${S}^{\left(0\right)}\left(0\right)={S}_{0}^{\left(0\right)}$ , ${S}^{\left(i\right)}\left(0\right)=0$ $\left(i=1,2,\cdots \right)$ ―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   .

5. Example

We consider an inhomogeneous viscoelastic two-component medium saturated with a liquid, characterized by parameters and relaxation nuclei (a special case)

$\begin{array}{l}{\Lambda }_{1}\left(x\right)=0.7{x}^{2},\text{\hspace{0.17em}}{\rho }_{11}\left(x\right)=0,6x,\text{\hspace{0.17em}}{\rho }_{12}\left(x\right)=-0.1x,\text{\hspace{0.17em}}{\rho }_{22}\left(x\right)=0.6x,\\ m\left(x\right)=0.5x,\text{\hspace{0.17em}}{R}_{0}\left(x\right)=0.2x,\text{\hspace{0.17em}}\Lambda \left(0,x\right)+2M\left(0,x\right)=0.7x\end{array}$ (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\sqrt{x}$ , ${c}_{{p}_{2}}=0.45\sqrt{x}$

${c}_{t}=0.87\sqrt{x}$

Since by the assumption of the problem the wave surfaces form a family of parallel planes, then the average curvature ${\Omega }_{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}^{\left(1\right)}={S}_{0p}^{\left(1\right)}\mathrm{exp}\left(-\underset{0}{\overset{x}{\int }}\left[{\beta }_{2}\left(x\right)\frac{\text{d}\mathrm{ln}\sqrt{c\left(x\right)}}{\text{d}x}+\frac{\left(1-0.5x\right)0.2+0.1\cdot {c}^{2}\left(x\right)}{2{a}_{1}\left(x\right){c}^{3}\left(x\right)}\cdot 0.17{x}^{3}\right]\text{d}x\right)$ (24)

where ${S}_{0p}^{\left(1\right)}$ the value of the function ${S}_{p}^{\left(1\right)}$ for $x=0$ , ${\beta }_{2}\left(x\right)$ , и ${a}_{1}\left(x\right)$ and are from (18).

By the formula (22) we obtain

${S}_{p}={S}_{0p}\left(1+{\Gamma }_{1}\left(x\right)\right){S}_{p}^{\left(1\right)}\left(x\right)$ , ${\Gamma }_{1}\left(x\right)=\frac{{\rho }_{11}\left(x\right){c}^{2}\left(x\right)-\left(\lambda \left(x\right)+2\mu \left(x\right)\right)}{\left(1-m\left(x\right)\right){R}_{0}\left(x\right)-{\rho }_{12}\left(x\right){c}^{2}\left(x\right)}$ (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.

6. 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.

Cite this paper

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

References

1. 1. Rossikhin, Yu.A. and Meshkov, S.I. (1969) Propagation of Shock Waves in Linear Hereditary Media. Sat. Some Problems of Solid State Physics. News of the State Pedagogical University, 1, 49-53.

2. 2. Blitshtein, Yu.M., Meshkov, S.I. and Chigarev, A.V. (1972) Propagation of Waves in a Linear Viscoelastic Inhomogeneous Medium. Izvestiya AN SSSR. Mechanics of a Solid Body, № 3, 40-47.

3. 3. Polenov, V.S. (2014) Propagation of Elastic Waves in a Porous Medium Saturated with a Viscous Liquid. Applied Mathematics and Mechanics, 78, 501-507.

4. 4. Biot, M.A. (1956) Theory Propagation of Elastic Waves in a Fluid-Saturated Porous Solid I. Low-Frequency Range. The Journal of the Acoustical Society of America, 28, 168-178. https://doi.org/10.1121/1.1908239

5. 5. Biot, M.A. (1956) Theory Propagation of Elastic Waves in a Fluid-Saturated Porous Solid. II. Higher Frequency Range. The Journal of the Acoustical Society of America, 28. 179-191. https://doi.org/10.1121/1.1908241

6. 6. Kosachevsky, L.Ya. (1959) On the Propagation of Elastic Waves in Two-Component Media. Applied Mathematics and Mechanics, 23, 1115-1123.

7. 7. Bykovtsev, G.I. and Verveiko, N.D. (1966) On the Propagation of Waves by an Elastic-Viscous-Plastic Medium. Mechanics of a Solid, № 4, 111-123.

8. 8. Thomas, T.Y. (1961) Plastic Flow and the Fracture in Solids. Academic Press, New York, 267 p.

9. 9. Polenov, V.S. and Chigarev, A.V. (2010) Propagation of Waves in an Inhomogeneous Porous Medium Saturated with a Liquid. Applied Mathematics and Mechanics. 74, 276-284.

10. 10. Chigarev, A.V. (1980) On the Geometry of Wave Fronts in Inhomogeneous Media, Acoustic Journal, 26, 905-912.

A List of Symbols

${Т}_{ij}$ : stress tensor of a porous medium

N: force acting on the liquid per unit area of the cross section of the porous medium

${e}_{ij}^{\left(1\right)}$ : components of the strain tensors of the elastic component

${e}_{kk}^{\left(2\right)}$ : components of fluid deformations

${u}_{i}^{\left(1\right)}$ : displacement of the elastic component

${u}_{1}^{\left(2\right)}$ : fluid movement

${\Lambda }^{\ast },{M}^{\ast }$ : linear elastic integral operators

${S}_{P}^{\left(\alpha \right)}$ : wave amplitudes

$\Omega$ : average curvature of the wave surface

${g}^{\alpha \beta }$ : the coefficients of the first fundamental quadratic form of the surface