Open access peer-reviewed chapter

# Mathematical Fundamentals of a Diagnostic Method by Long Nonlinear Waves for the Structured Media

By Vyacheslav Vakhnenko, Dmitri Vengrovich and Alexandre Michtchenko

Submitted: November 26th 2019Reviewed: February 1st 2020Published: March 12th 2020

DOI: 10.5772/intechopen.91462

## Abstract

We have proven that the long wave with finite amplitude responds to the structure of the medium. The heterogeneity in a medium structure always introduces additional nonlinearity in comparison with the homogeneous medium. At the same time, a question appears on the inverse problem, namely, is there sufficient information in the wave field to reconstruct the structure of the medium? It turns out that the knowledge on the evolution of nonlinear waves enables us to form the theoretical fundamentals of the diagnostic method to define the characteristics of a heterogeneous medium using the long waves of finite amplitudes (inverse problem). The mass contents of the particular components can be denoted with specified accuracy by this diagnostic method.

### Keywords

• diagnostics
• nonlinear waves
• inverse problem
• asymptotic model
• structured medium

## 1. Introduction

Natural media, in the general case, should not be treated as structureless. The experiments have shown that the intrinsic structure of a medium influences the wave motions [1, 2, 3, 4, 5, 6, 7, 8]. Existing inhomogeneities complicate the problem and, at the same time, are fully manifested under the propagation of nonlinear waves. The principal part of the problem is associated with the phenomena caused by the nonlinear behavior of natural media, in particular, more substantial increase of nonlinear effects in structured media than that in homogeneous ones [1, 2, 3, 4, 5, 7].

## 2. The model notion of the medium with structure

The wave processes in heterogeneous media are usually described in terms of more or less complicated models. Under the conditions of local equilibrium, the media are traditionally modelled irrespective of their structure. In the framework of continuum mechanics, the known idealization of a real medium as a homogeneous one has been relatively successful in the description of wave processes (see, e.g., [9, 10, 11]). The continuum models are commonly applied to the mixtures whose dispersive dissipative properties are treated with regard to the interactions between the components [12, 13, 14, 15]. On this level, the media are modelled in the framework of a homogeneous elastic, viscous elastic, and elastic-plastic media. In this case, the features of the medium structure are taken into account indirectly through the kinetic parameters (relaxation time, viscous coefficients, etc.) [4, 5, 6, 10, 12, 13, 14, 15].

The model of multi-velocity interpenetratable continua was developed in terms of classical continuum mechanics [16] and statistical physics in order to describe the dynamical behavior of multicomponent media. A fundamental assumption in the theory of mixtures [15] reproduces the assumption in the model of multi-velocity interpenetratable continua [16], namely, that the particles of the existing components occupy each microvolume. The equations of motion for each component involve the terms describing the mass, force, and energy interactions between the components. The problem is complicated by the necessity to employ, in the general case, the experimental data for establishing theoretical relations between the macro parameters at the component interaction level. Moreover, if the component interaction is determined, then these models would be indispensable in the theory of multicomponent media.

In all the models mentioned, the formalism of continuum mechanics is based on the principle of local action as well as on the generalization of the mechanics laws relating the point mass to the continuum [11].

When going from the integral equations to differential balance equations, the existence of a differentially small microvolume dvis assumed. On the one hand, this volume is so small that the mechanics laws of the point can be extended to the whole microvolume. On the other hand, the volume contains so many structural elements of the medium that, in this sense, it can be regarded as macroscopic one despite its smallness as compared to the entire volume occupied by the medium. So, the passage to the differential balance equations is based on the assumption that microstructural scales εare small as compared to the characteristic macroscopic scale of the λand the passage should be made to the limiting case ε/λ0. Contraction of the volume dvto the point is in the general case correct for continuous functions [11, 15]. This means that all points within the differentially small volume are equivalent. Hence, for the case of a mixture, the equivalence of the points implies that field characteristics should be averaged over dv. Hence, it is assumed that the equations of motion can be written in terms of average density, mass velocity, and pressure of each component. We note that these models do not contain the exact sizes of components.

The application of the models of a homogeneous medium to the description of the dynamical wave processes in a structured natural medium is associated with specific fundamental difficulties [3, 4, 6, 8]. In what follows, we treat the medium structure at the macrolevel. We abandon the assumption that the differentially small volume dvcontains all the components of the medium. Nevertheless, we consider the long-wave approach with the wavelength λmuch higher than the characteristic length of the medium structure ε. We consider a structured medium (Figure 1) in which separated components are considered as a homogeneous medium (the differentially small volume dvis much smaller than the characteristic size of a particular component ε).

Within continuum mechanics [17], the known idealization of a real medium as homogeneous is widely used for modelling their dynamic behavior. In these models, the effect of heterogeneity is taken into account indirectly throughout the kinetic parameters such as a viscous coefficient and relaxation time. The inner processes, in this case, are manifested through dispersive dissipative properties of a medium. Traditionally, it was considered that in heterogeneous media with wavelength appreciably exceeding the size of the structural heterogeneities, the perturbations propagate in the same way as inhomogeneous media [9, 10, 15]. However, this statement should be proven, and we shall show that this approximation is not universally correct. In general case, the wave propagation cannot be described in terms of the average characteristics (continuum model).

The properties of a medium deviate from the equilibrium state under the propagation of intensive waves. Moreover, an unperturbed medium can be in one of the unstable stationary states. So, a geophysical medium, within a current physical concept, is an open thermodynamic system, which substantially influences the exchanges of energy and mass. Thus, a description of open systems should take into account the peculiarities of their inner structure, dynamical processes occurring on the level of structural elements. What is more, the state of media under the action of high-frequency wave perturbations departs from equilibrium, and, thus, the behavior of media cannot be described in the framework of equilibrium thermodynamics. Consequently, there is a necessity to develop new mathematical models to take into account the nonlinear wave perturbations and irreversible inner exchange processes.

## 3. Asymptotic averaged model for structured medium

We describe the wave processes in nonequilibrium heterogeneous media in terms of an asymptotic averaged model [18, 19, 20, 21, 22]. The obtained integral differential system of equations cannot be reduced to the average terms (pressure, mass velocity, specific volume) and contains the terms with characteristic sizes of individual components.

### 3.1 Background and initial equations

The most straightforward heterogeneous media for which the effect of the structure can be analyzed are media with a regular structure. Features of the propagation of long-wave perturbations will be investigated by using, as an example, a periodic medium under conditions of equality of stresses and mass velocities on the boundaries of neighboring components. It is supposed that the microstructure elements of medium dv(see Figure 1) are large enough that it is possible to submit to the laws of classical continuum mechanics for each component. At the same time, the inner processes in each component will be considered within a relaxation approach. The notions based on the relaxation nature of a phenomenon are regarded to be promising and fruitful. We consider that the properties of the medium, such as density, sound velocity, and relaxation time, vary periodically (although this assumption is unessential in the final result).

#### 3.1.1 Motion equations for individual component

The analysis of wave motions is based on the hydrodynamic approach. This restriction can be imposed for the modelling of nonlinear waves in water-saturated soils, bubble media, aerosols, etc. [13]. The set of acceptable media could be extended to solid media where the powerful loads are studied in the condition that the strength and plasticity of the material can be neglected [23]. In the hydrodynamic approach, we have considered the media without tangential stresses, while there are equalities of the stresses as well as of mass velocities on boundaries of neighboring components. Also, we assume that the medium is barotropic. The individual components of the medium are considered to be described by the classical equations of hydrodynamics. In the Lagrangian coordinate system lt, the equations of one-dimensional motion for each component have the form

rνlν=VV0,u=rt,ut+V0rlν1pl=0.E1

The equation of continuity can also be used in the alternative form

VtνV0rν1ulν=0.E2

Here V=ρ1is the specific volume; νis a parameter of symmetry, where ν=1is planar symmetry, ν=2is cylindrical one, and ν=3is spherical one; and index 0 relates to the initial state. The other notations are those that are generally accepted.

Conditions for matching are the equality of mass velocities and pressures on the boundaries of the components:

u=0,p=0.E3

#### 3.1.2 Dynamic state equation

Considering the models of a relaxing medium as more general than the equilibrium models for describing the evolution of high-gradient waves, we will take into account the relaxing processes for each component. Thermodynamic equilibrium is disturbed owing to the propagation of fast perturbations in a medium. There are processes of the interaction that tend to return the equilibrium. The parameters characterizing this interaction are referred to as the inner variables, unlike the macro parameters such as the pressure p, mass velocity u, and density ρ. In essence, the change of macro parameters caused by the changes of inner parameters is a relaxation process. From the nonequilibrium thermodynamics standpoint, the models of a relaxing medium are more general than the equilibrium models for describing the wave propagation.

An equilibrium state equation of a barotropic medium is one-parameter equation. As a result of relaxation, an additional variable ξ(inner parameter) appears in the state equation. It defines the completeness of the relaxation process:

p=pρξ.E4

There are two limiting cases:

1. The lack of relaxation (inner interaction processes are frozen) ξ=1,

p=pρ1=pfρ.E5

• The relaxation complete (there is the local thermodynamic equilibrium) ξ=0,

p=pρ0=peρ.E6

• The equations of state (5) and (6) are considered to be known. These relationships enable us to introduce the sound velocities for fast processes:

cf2=dpf/E7

and for slow processes

ce2=dpe/.E8

The slow and fast processes are compared, utilizing relaxation time τp. The dynamic state equation is written down in the form of the first-order differential equation:

τpdtcf2dpdt+ρρe=0.E9

The equilibrium equations of state are considered to be known:

ρeρ0=p0pce2dp.E10

Clearly, for the fast processes ωτp1, we have the relation (5), and for the slow ones ωτp1, we obtain (6).

The substantiation of Eq. (9) within the framework of the thermodynamics of irreversible processes has been given in [17, 24, 25, 26]. As far as we know, the first work in this field was paper by Mandelshtam and Leontovich (see Section 81 in [17]). We note that the mechanisms of the exchange processes are not explicitly defined when deriving Eq. (9), and the thermodynamic and kinetic parameters appear only in this equation. These characteristics can be found experimentally.

The phenomenological approach for describing the relaxation processes in hydrodynamics has been developed in many publications [13, 17, 26]. The dynamic equation of state was used (a) for describing the propagation of sound waves in a relaxing medium [17], (b) for taking into account the exchange processes within media (gas-solid particles) [26], and (c) for studying wave fields in gas-liquid media and in soil [13]. In most works, the equation of state has been derived from the concept of precise mechanism for the inner process. Within the context of mixture theory, Biot [12] attempted to account for the nonequilibrium in velocities between components directly in the equations of motion in the form of dissipative terms.

We assume that the relaxation time and sound velocities do not depend on time, but they are functions of pressure and the individual properties of the components. This means that in the process of a relaxation interaction, we can take into account the exchange of moment and heat but not that of mass. The dynamic equation of state determines the features of the intrastructural interaction for each component.

The equations of motion (1) have been written in the Lagrangian coordinate system. The necessity of such a description stems from the fact that the dynamic equation of state (9) has been written to the mass element of a medium. Besides, the use of the Lagrangian coordinates is essential for the application of the method of asymptotic averaging, since in these coordinates, the structure is independent of a wave process.

### 3.2 Asymptotic averaged system of equations

A regularity of structure and a nonlinearity of long-wave processes investigated here specify the choice of mathematical methods. One way of studying this heterogeneous medium is based on a method of asymptotic averaging of equations with high-oscillating coefficients [20, 27, 28, 29, 30]. The essence of this method consists in the application of a multiscale method in combination with a space averaging. In accordance with this method, the mass space coordinate m=lν/V0is divided into two independent coordinates: slow coordinate sand fast one ξ, wherein

m=s+εξ,m=s+ε1ξ.E11

The slow coordinate scorresponds to a global change of the wave field and sis a constant value during a period, while the fast coordinate ξtraces the variations of a field in the structure period. The dependent functions are presented as a degree series over the structure period ε

Vmt=V0stξ+εV1stξ+ε2V2stξ+pmt=p0stξ+εp1stξ+ε2p2stξ+umt=u0stξ+εu1stξ+ε2u2stξ+.rνmt=rν0stξ+εrν1stξ+ε2rν2stξ+E12

where pi, ui, Vi, and riare defined as the one-period functions of ξ. In the Lagrangian mass coordinates the period is a constant which allows the averaging procedure to be performed.

We now will prove that p0=p0st, p1=p1st, u0=u0st, and rν0=rν0stare independent of the fast variable ξ. Indeed, after substitution of Eqs. (11) and (12) into the initial equations of motion, we obtain

ε1rν0ξ+ε0rν0srν1ξV0+=0,ε0u0r0t+=0,ε1νrν10p0ξ+ε0u0t+νrν10p0s+νrν11p0ξ+νrν10p1ξ+=0,ε1νrν10u0ξ+ε0V0t+νrν10u0sνrν11u0ξνrν10u1ξ+=0,E13

According to the general theory of the asymptotic method, the terms of equal powers of εshould vanish independently of each other. Thus, p0/ξ=0, u0/ξ=0, and rν10/ξ=0, i.e., p0=p0st, u0=u0st, and r0=r0st, are independent of ξ. Furthermore

rν0s+rν1ξ=V0,u0=r0t,u0t+νrν10p0s+νrν10p1ξ=0,V0tνrν10u0sνrν11u0ξνrν10u1ξ=0.E14

Thus, we can average the equations during the period ξ. We define =01and perform the normalization 01=1. Since p1, u1, and r1are periodic, the integrals can be calculated as p1/ξ=0, u1/ξ=0, and r1/ξ=0. Moreover, as u0=u0, p0=p0than p1/ξ=0. This means that p1does not also depend on ξ. After integrating over the structure period the equations containing the value of zero order of ε, we obtain the averaged system

rν0s=V0,u0=r0t,u0t+νrν10p0s=0,V0tνrν10u0s=0E15

with the averaged equation of state

dV0=V02cf2dpV0τpVep0V0Vep0dt.E16

Unlike the values u0, p0, p1, and r0, the specific volume V0is a function of ξ. Hereafter, we will consider only the zero approximation of the equations, and, therefore, the upper index 0 is omitted. Choosing the wavelength λto be large enough, we can always reduce the effect to zero from other approximation terms.

The averaged system of Eqs. (15) and (16) is an integrodifferential one and, in the general case, is not reduced to the averaged variables p, u, and V0. The derivation of Eqs. (15) and (16) relates to a rigorous periodic medium. However, it may be shown that Eqs. (15) and (16) are also relevant to media with a quasiperiodic structure. Indeed, the pressure pand the mass velocity uare independent of the fast variable ξ. Hence on a microscale ξ, the action is statically uniform (waveless) over the whole period of the medium structure, while on the slow scale s, the action of perturbation is manifested by the wave motion of the medium. On a microlevel, the behavior of medium adheres only to the thermodynamic laws. There is a mechanical equilibrium. On a macrolevel, the motion of the medium is described by the wave dynamics laws for averaged variables. Mathematically, in the zero-order case of ε, the size of the period is infinitesimal ε0. This signifies that the location of particular components in the period is irrelevant. Eqs. (15) and (16) do not change their form if the components are broken and/or change their location in an elementary cell. This means that Eqs. (15) and (16) describe the motion of any quasiperiodic (statistical heterogeneous) medium which has a constant mass content of components on the microlevel and the location of these components within the cell is not important.

In the case of nonlinear wave propagation, the individual components suffer different compressions. The structure of the medium is changed, with the result that the averaged specific volume Vis changed. This change differs from the change of the specific volume for homogeneous medium under the same loading. Thus, the structure of medium is manifested in the wave motion, despite the fact that the equations of motion (15) (but not the equation of state) are written down for the averaged values u, p, and Vonly.

### 3.3 System of equations in Eulerian coordinates

In certain cases of theoretical analysis, it is more convenient to use the Eulerian coordinate system. The immediate employment of a method of the asymptotic averaging in Eulerian variables is impossible because of the variability of the microstructure sizes. However, from the zero approximation in the equations of motion (11), which are presented by the averaged values p, u, and V, the equations can be rewritten in the Eulerian system of coordinates rtEutilizing a transformation from the Lagrangian system st[18, 19, 20, 21, 22]:

r=rst,tE=t.E17

There is an essential presumption that the velocity of the particle in the zero approximation is constant throughout of the structure, and, consequently, we can describe an averaged trajectory for the particle:

rstts=ust.E18

From the physical point of view, it is clear that the position of the particle is unambiguously defined by its coordinate and time:

From the mathematical point of view, this means that in the transformation (19) the value drνis a total differential. Therefore, we must have

At=νrν1us.E20

This condition is satisfied if A=V, because the equation converts into the continuity Eq. (15). We obtain the following transformation between Lagrangian and Eulerian systems of coordinates:

drν=Vds+νrν1udt,tE=t.E21

It is reasonable to define the slow Lagrangian coordinate (non-mass one) as

Rν=sV.E22

Eq. (15) in the Eulerian system of coordinates then take the form

V1tE+rν1uV1r=0,utE+uur+Vpr=0.E23

It is convenient to determine the fast Eulerian coordinate ζas

ζξt=ρ˜ρξ.E24

It should be noted that the average density ρ˜in the Eulerian coordinates is a value usually used for density. A chain of identities

V=01Vξ=01Vρρ˜=ρ˜1E25

proves that V1is the average density of the medium in the Eulerian coordinates. Note that ρ˜ρ. The value ρ˜is a real density. The value Vis the specific volume averaged in units of mass over the period, and it is expressed as the ratio of the volume to the mass inside this volume. This value can be determined experimentally. At the same time, the averaged values pand ucoincide in both Lagrangian and Eulerian systems of coordinates. Now the equations of motion (23) can be written in the usual form of the averaged density ρ˜.

The notation of the equations of motion in the averaged values enables us to suggest the method of the computer solution for the system of equations, where the integration step is restricted by the perturbation wavelength and not by the period of the structure [22]. Then the main computational problem associated with the smallness of the integration step can be avoided, and the equations of motion can be solved at a significant distance of wave propagation within a reasonable time.

### 3.4 Nonlinear waves

We will analyze the propagation of nonlinear waves in a structured medium. To make the results more clear, we will restrict our consideration to a nonrelaxation media c=cf=ce. The averaged equation of state in this case is simplified to the form

dV=V2c2dp,E26

and we can introduce an effective sound velocity by the formula

ceff=V2/V2c2.E27

We obtain a traditional representation of the system of Eqs. (15) and (26).

The system of the equations is concerned in the hyperbolic type of a system. Now we restrict ourselves to the plane symmetry ν=1. Substituting the equation of state (26) into the equation of the continuity (the last equation in (15)), we get

V2c2pt+us=0.E28

The combination of this equation with the third equation in (15) ν=1leads to the relationships

ut±V2c21/2pt±V2c21/2us±V2c21/2ps=0.E29

From this relationship, it is seen that the averaged system of the equations pertains to the hyperbolic system. The equations for the characteristic in Lagrangian coordinates (mass space coordinate) have the forms

dsdt=±V2c21/2.E30

In characteristic, the relations are the following:

I±=u±V2c21/2dp.E31

Analogously to the homogeneous medium, we call these relations as the Riemann invariants. The value (30) has the physical meaning, namely, it is the averaged velocity of the wave propagation in the Lagrangian coordinates. This velocity depends on pressure and integrally on a structure. Note the particular case. It is known that in vacuum the wave does not propagate. This result also follows formally from Eq. (30). The hyperbolism of a system points up that this system can describe the shock wave. The equations for the characteristic (30) and the Riemann invariants (31) are the integrodifferential equations, since they retain the variable V2/c2, which depends on the properties of the structural elements in medium.

It should be noted that ceffis not an averaged value, i.e., ceff2c2. Evidently, the structure of the medium introduces a certain contribution to the nonlinearity. In fact, even if cffp, then in the general case, the value of ceffis a function of pressure.

The system of Eq. (15) is hyperbolic ones, and this specifies the breaking solutions, which are shock waves. For the analysis of such solutions, it is necessary to present Eq. (15) in the form of integral conservation laws:

Vds+udt=0,udspdt=0.E32

Now we can quickly formulate the conditions on the shock front, when there is conservation of the fluxes of mass and impulse through the shock front:

V1V0D+u1u0=0,u1u0Dp1+p0=0,E33

where indexes 0 and 1 relate to the parameters of the flow before and after the front, respectively. Hence, the formula for the averaged velocity of the shock front in terms of the Lagrangian variable D(dimension Dis kg/s) and the mass velocity ufollow from the following relations:

D=p1p0/V0V1,u1u0=p1p0V0V1.E34

These systems are not expressed in the average hydrodynamical terms; hence the dynamical behavior of the medium cannot be modelled by a homogeneous medium even for long waves, if they are nonlinear. The structure of the medium influences the nonlinear wave propagation.

In the next section, it will be proven that the heterogeneity of the medium structure always introduces additional nonlinearity that does not arise in a homogeneous medium. This effect makes it possible to formulate the theoretical grounds of a new diagnostic method that determines the characteristics of a heterogeneous medium with the use of finite-amplitude long waves (inverse problem). This diagnostic method can also be employed to find the mass contents of individual components.

## 4. Diagnostics of a medium by long nonlinear waves

In the previous section, we proved that the long wave with finite amplitude responds to the structure of medium. At the same time, a question appears, namely, is there sufficient information in the wave field to reconstruct the structure of medium? It turns out that the knowledge on the evolution of nonlinear waves enables one to define with certain accuracy the concentrations of medium components.

### 4.1 The increase of nonlinearity in medium with structure

In this section, we shall prove the statement that the structure of medium always exalts the nonlinear effects under the propagation of long waves. At first, let us consider the sound velocity in homogeneous chomand heterogeneous ceffmedia. Now we will show that in the general case with pressure increase the velocity of the sound becomes higher in a structured medium than in a homogeneous one:

ceffchom.E35

For the sake of clarity, we consider a medium in which the sound velocities of individual components are independent of the pressure:

cfp,dc/dp=0.E36

The equality sign is fulfilled (a) for an initial pressure, by virtue of the normalization, and also (b) for a special structured medium in which the relation Vξ/c2ξis not a function on the fast variable ξ. We must prove which case results in equality and which gives the inequality.

Let us write the relations (36) for homogeneous medium consisting only one component:

chomfp,dchom/dp=0.E37

For multicomponent medium, the derivative dceff/dpis defined from the relationship

dceffdp=2VV2/c2VV3c4V2c220.E38

This last inequality follows from the well-known Cauchy-Schwarz inequality (see, e.g., [31]). Therefore, with the increase of pressure, the sound velocity ceffincreases. Consequently, we have the inequality (35) at pp0.

Moreover, at p>p0the shock adiabatic curve for the medium with a structure always lies above that for the homogeneous medium (they touch only at the initial point p=p0):

d2pdV2d2pdV2hom.E39

Indeed, a ratio of these derivatives is equal to

d2pdV2/d2pdV2hom=V3/c4V2/c23chom2V3=V3/c4Vceff2chom2V2/c22V3/c4VV2/c221.E40

Hence, a long wave with a finite amplitude responds to the structure of the medium, and the nonlinear effects increase as compared with those in the homogeneous medium. The nonlinearity takes place even if individual components are described by the linear evolution equation (i.e., at condition (36)).

The exception, as it was noted already, is a medium with the properties of structure Vξ/c2ξfξ. For this medium, only the equality sign is correct in the inequalities (35) and (39). Particular elements of the structure respond to the pressure variations, but the relative structure does not change, i.e., the ratio Vξp/Vξp0does not depend on ξ. In this case, the value ceff=c2is an averaged characteristic (see Eq. (27)). Therefore, the system of equations may be presented using the averaged variables p, u, V, and ceff=c2. Heterogeneity does not introduce an additional nonlinearity for this medium, and the structure of medium does not affect the wave motion.

In addition to the analysis of the sound velocity in homogeneous and heterogeneous media, we consider now the evolution equations with the nonlinear term and compare the coefficients of nonlinearity in these media. Let us derive the evolution equation with weak nonlinearity. First of all, we have to note that the mass velocity uis related to the pressure pby means of [22]

u=popV2/c2dp.E41

Functional dependence of an average specific value on the pressure increment p=pp0with the accuracy Op2can be presented as a series:

Vp=V0+dVdpp=p0p+12d2Vdp2p=p0p2.E42

In this case, the system of Eq. (23) for planar symmetry ν=1can be written as

V0ux+V2c20pt12d2Vdp2p=p0p2t=0,E43
ut+V0px=0.E44

The relationship upx=puxfollows from Eq. (41) with the assumed accuracy Op2and was used for derivation of the first equation. The evolution equation for one variable assumes the form

V022px2V2c202pt2+12d2Vdp2p=p02p2t2=0.E45

Now let us consider the waves propagating in one direction, and then with the indicated accuracy, we can write (hereinafter index 0 is omitted):

V2/c2Vt+x2xE46

(see, e.g., Section 93 in Ref. [17]). Thus, after factorization of Eq. (45) we get

pt+ceffpx+12VV2c23/2d2Vdp2ppx=0.E47

The coefficient of nonlinearity αpfor the structured medium, when the sound velocities in the individual components are independent of the pressure cfp, can be presented as

αp12VV2c23/2d2Vdp2=du+ceffdp=VV3c4V2c23/2.E48

For all cases we take αp>0. For a homogeneous medium with dc/dp=0, we have αphom=V/c.

In certain media the value V/c2does not change within the period. The individual elements of the structure respond to the pressure variations so that a relative structure does not change, i.e., the ratio Vξp/Vξp0does not depend on ξ. In this case, the value ceff=c2derived from Eq. (27) is the averaged characteristic. Consequently, the system of equations may be presented in the averaged variables p, u, V, and ceff=c2. Heterogeneity does not introduce the additional nonlinearity for these media. Such media behave like the homogeneous media under the action of the nonlinear wave perturbations.

For media, when the sound velocity is independent of the pressure cfp, it is possible to show that heterogeneity of the medium, in the general case, introduces the additional nonlinearity. Let us consider the ratio of the nonlinearity coefficients for heterogeneous and homogeneous media. In the space of dimensionless normalized variables, this implies that at p=p0we have V0=1as well as V2/c20=1for the compared media.

Using the conditions (27) we can obtain

αpαphom=VV3c4V2c221.E49

This inequality is the well-known Cauchy-Schwarz inequality (see formula (15.2-3) in Ref. [31]). Since V0and V/c20, we prove

VV3/c4VdξV3c4=V2c2Vc21V2c2Vc2V2c2Vc21V2c2Vc22=V2c22V2/c22.E50

It only remains to find the condition for the equality sign in (49). For this purpose, we apply the Cauchy-Schwarz inequality in vector form (see formula (15.2-5) in Ref. [31]):

ab2aabb.E51

However, the equality sign is realized if and only if the vectors aand bare linearly dependent, i.e., a=kb(k=const). By designating aaV/c2and bbV2/c2, it is easy to notice that the equality sign is realized if and only if

V2c2Vc21/V2c2Vc2=const.E52

(see sections 14.2-6 in Ref. [31]), i.e., when the value V/c2=constdoes not vary within the period (Vξ/cξ2fξ). This heterogeneous medium has been considered above. For all other heterogeneous media for which the value V/c2changes within period, the inequality is realized in Eq. (49). So, in a heterogeneous medium, the value αpis always greater than αphomin a homogeneous medium. Thus, it is proved that, in the general case, the heterogeneities in a medium introduce the additional nonlinearity. This effect provides the basis for a new method of diagnostics to define the properties of multicomponent media using the propagation of long nonlinear waves in such media.

### 4.2 Fundamentals of new diagnostic method

The structure of the medium affects the wave field. There are different methods which allow the detection of gas bubbles and/or cracks in liquid [32], concrete [33], and ice cover [34] employing the nonlinear effects.

In this section, we describe our new diagnostic method for the properties of the medium. The features of the motion of finite-amplitude long waves and the effect of the increase of nonlinearity in the heterogeneous medium in comparison with homogeneous medium form the basis for the development of theoretical fundamentals of the diagnostic method. In this method, the properties of individual components are defined by long waves of finite amplitudes; more specifically, the dependence V/c2=V/c2ζon the fast Eulerian coordinate ζ(see Eq. (24)) is defined.

Thus, the nonlinear wave evolution allows one to obtain the structure of the medium with an inherent accuracy. As a final result, the mass concentrations of the individual components can be found using this method.

It should be kept in mind that the period of the structure of medium is infinitely small in the long-wave model, so it is not always possible to indicate the location of the structure elements inside the period reliably. Hence, the media with the different structures plotted in Figure 2, for example, affect identically on wave fields. These two media are indistinguishable in the framework of the suggested method. Taking into account this indefiniteness, we consider the function V/c2=V/c2ζthat is to be the decreasing, integrable, mutually one-valued function on the interval ζ01, and equal to zero outside of this interval.

Now we represent the theoretical fundamentals for new method of diagnostics of medium by means of the long nonlinear waves. Let us prove the principal relation which enables us to obtain the inverse function ζ=ζV/c2for the desired function V/c2=V/c2ζthrough the inverse Fourier transformation [19, 35, 36, 37]:

ζVc2=F1n=0VVc2n+1n+1!VinqnVc2.E53

It is known from theory of probability that the distribution function fx(any one-valued, integrable, positive function) can be expressed by its central moments:

αn=xnfxdx.E54

Indeed, by using the characteristic function

χq=Ffxq,E55

any positive integrable function fxcan be written as follows:

fx=F1χqx,E56

where Fis the Fourier transformation and F1is the inverse Fourier transformation.

We take into account the important fact from the theory of probability: the characteristic function χqis uniquely determined by the central moments αn:

χq=n=0αninqnn!.E57

Hence, the function fxcan be found by means of the inverse Fourier transform:

fx=F1n=0αninqnn!x,E58

if series n=0αnsn/n!converges absolutely for some value s>0(see Section 18.3.7 in Ref. [31]).

These facts from the theory of probability are used to prove such a statement: if V/c2=V/c2ζis a decreasing positive integrable function on the interval ζ01and equals to zero outside of it, then the inverse function ζ=ζV/c2for the required function V/c2=V/c2ζcan be written as (53) in the averaged values

VV/c2nVV/c2n.E59

Indeed, for the monotonic one-valued function V/c2=V/c2ζ, we find the integral in (59) by integrating the inverse function ζ=ζV/c2, since the transformation Jacobian is not equal to zero. We have the chain of identifies

VV/c2n=01VξVc2n=V01VVc2nρdζ=VVc2ndV/c2dV/c2.E60

In the geometric sense, this relation signifies that the integral (in our case, it is an area between the curve V/c2=V/c2ζand axes and OV/c2) can be calculated either over ζor over V/c2(see Figure 2). Whereas, the inequality is realized for the monotonic decreasing function V/c2=V/c2ζ.

For a function defined on a finite interval, if this function is positive and bounded above, we have

VV/c2n=VVc2ndV/c2dV/c2=nVVc2n1ζdV/c2.E61

This relation provides the connection between the central moment αnand the value VV/c2n

VV/c2n=nVαn1.E62

Then the characteristic function χqfor the inverse function ζ=ζV/c2is expressed through VV/c2n. By applying the inverse Fourier transformation, finally, we find the required relationship (53).

The physical value Vc2is bounded by some constant M, hence

αn=Vc2n1ζdV/c20MVc2ndV/c2=Mn+1n+1.E63

The series n=0αnsn/n!n=0Mn+1sn/n+1!converge at s<M1. Consequently, the power series (53) also converges.

The coefficients VV/c2n(n=3,4,) in Eq. (53) can be easily calculated, if we know the functional dependence Vpor V2/c2p. Indeed, they can be successively defined by the recurrence relation

dVVc2ndp=n+1VVc2n+1,E64

that follows directly from the equation of state. With mentioned accuracy, it is possible to diagnose the structural properties of the medium.

We have proven the principal relation (53) for the method of diagnostics that allows one to find the properties of the individual components in structured media by means of the long nonlinear waves.

### 4.3 Approximation of diagnosed medium by layer medium

Diagnostics of the structured medium properties by the long nonlinear waves is connected with the definition of values VV/c2n. As indicated above, there is a problem related to the accuracy of the description of the structure by finite series (53).

Now, we shall show that the partial sum of series (53) is a step-function and approximates the desired function ζ=ζV/c2with certain accuracy, namely, the diagnosed medium can be approximated by a layer medium. Let us write down the chain of the identities for any integrable function:

2πfx=FFfxqx=Fn=0inqnn!αn=n=0inαnn!2πinδnx.E65

Here we used the known relationships for the Fourier transform [31]):

FFfxqx=2πfx,E66
Fqnx=2πinδnx.E67

Hence, any integrable function can be represented by a series:

fx=n=0αnn!δnx.E68

We will prove that the finite series (53) approximates the desired function fxby step-function. Consider the step-function f1xconsisting of Nsteps:

f1x=φ1,0<xb1,φ2,b1<xb2,φN,bN1<xbNE69

in order to approximate the desired function fx. The relation (69) can be written down through the Heavyside functions as follows:

f1x=φ1ΘxΘxb1+φ2Θxb1Θxb2++φNΘxbN1ΘxbN,E70

Evidently, by increasing the number of steps Nand choosing the values φiand bi, any integrable function fxcan be approximated by the step-function f1x. It is convenient to use a notation

f1(x)=φ1[Θ(x+b1)Θ(x)]+φ2[Θ(x+b2)Θ(x+b1)]++φN[Θ(x+bN)Θ(x+bN1)],E71

that follows immediately from (70) after substitution:

Θx=1Θx.E72

The Heavyside function Θx+bcan be expanded into a Taylor series in the neighborhood of point x:

Θx+b=Θx+n=1bnn!Θnx.E73

It is well-known that the derivative of Heavyside function Θxis δx-function, then Θn+1x=δnx. We equate functions (68) and (71) and consider that the number of steps for function f1xis infinitely larger, and in this case we obtain

φ1n=0b1n+1n+1!δnx+φ2n=0b2n+1b1n+1n+1!δnx+++φNn=0bNn+1bN1n+1n+1!δnx+=n=0αnn!δnx.E74

This relationship shows that when we use the partial sum of series on the right-hand side of Eq. (74) n=02N1αnn!δnxand also the Nleading terms on the left-hand side, then the desired function fxis approximated by the step-function f1xwith Nsteps. In other words, if it is necessary to restore the structure of medium by means of Nperiodic repeated layers, then we need to know the 2N1moments αn, i.e., the values VVc2n.

For the sake of convenience, we write down the relation (74) in the expanded form. For this purpose, we multiply it by xnand integrate over x. We obtain the nonlinear system of the equations in the unknowns b1, b2, ,bN, φ2, φ3, and ,φN(variable φ1=1owing to normalization):

φ1b1+φ2b2b1+φ3b3b2++φNbNbN1=α0,φ1b12+φ2b22b12+φ3b32b22++φNbN2bN12=2α1,φ1b12N1+φ2b22N1b12N1+φ3b32N1b22N1+++φNbN2N1bN12N1=2N1α2N2.E75

Now, if biimplies the partition of V/c2iand φiimplies the partition of ζi, we can obtain the system of Eq. (75) to define the structure of medium. Solution of these equations gives the information about the component properties of the medium, namely, the value V/c2on the structure period ζ01is found in the form of the step-function.

Let us note the special case of a periodic medium for which the value V/c2is constant within the period. This medium, as we already know, does not differ from a homogeneous one for the propagation of the long nonlinear waves. The same result follows from a system (75). Indeed, for homogeneous media the moments αnare equal to

αn=VVc2n+1n+1V=bn+1n+1.E76

Here, the conditions of normalization V2/c20=V2/c20=1and V0=V0=1have been used as before. Therefore, the values in the right-hand side of Eq. (75) are equal to bVc2=const. It is easy to see that the solution of system is b1=b2==bN=b=1and φ1=1(where φiis any value for i2). This corresponds to the layer medium, for which V/c2fζ, in particular, this medium can be a homogeneous one.

According to the asymptotic averaged model of a structured medium, the period of the structure is infinitely small, and this diagnostic method cannot give the exact location of the structure elements inside the period. Hence, using this method, only the mass contents of the particular components can be determined.

We present, as an example, the results of the calculation to define the structure of layer media, which can properly approximate the diagnosed medium. The structure of the diagnosed medium is V/c2=0.2+0.81ζ2in Figure 3. In order to approximate the diagnosed medium by layer periodic medium, which has Nlayers within the period, it is necessary to know 2N1values VVc2nfor finite series (53). If we regard that the 2N1averaged characteristics VVc2ncoincide for the diagnosed medium and the layer medium, these averaged values at n2N1can be calculated from the known distributions V/c2=0.2+0.81ζ2. At n>2N1the values VVc2nfor diagnosed medium and for approximated layer medium are different. The distributions of V/c2ζwithin the period for diagnosed medium and for approximated media with Ncomponents are shown in Figure 3. On the one hand, the calculated distributions for layered media are the best approximation for the medium we test. On the other hand, we have illustrated the accuracy of the approximation of the diagnosed medium by the finite series (53).

Thus, the new method for the diagnostics of the medium characteristics by long nonlinear waves is suggested on the basis of the asymptotic averaged model of the structured medium. The mass contents of the particular components can be denoted by the abovementioned diagnostic method.

## 5. Conclusion

The asymptotic averaged model is suggested for the description of the wave processes in nonequilibrium heterogeneous media. The obtained integral differential system of equations cannot be reduced to the average terms (pressure, mass velocity, specific volume) and contains the terms with characteristic sizes of individual components.

On the microstructure level of the medium, the dynamical behavior is governed only by the laws of thermodynamics. On the macrolevel, the motion of the medium can be described by the wave-dynamical laws for the averaged variables with the integrodifferential equation of state containing the characteristics of the medium microstructure. A rigorous mathematical proof is given to show that finite-amplitude long waves respond to the structure of the medium in such a way that the homogeneous medium model is insufficient for the description of the behavior of the structured medium. An important result that follows from this model is that, for a finite-amplitude wave, the medium structure (in particular, existence of microcracks) produces nonlinear effects even if the individual components of the medium are described by a linear law. Finding the wave fields in the structured medium is the direct problem, on the one hand.

On the other hand, the system analyzed here is not expressed in the average hydrodynamical terms; hence the dynamical behavior of the medium cannot be modelled by a homogeneous medium even for long waves, if they are nonlinear. The heterogeneity of the medium structure always introduces additional nonlinearity that does not arise in a homogeneous medium. This effect enabled one to formulate the theoretical grounds of a new diagnostic method that determines the characteristics of a heterogeneous medium with the use of finite-amplitude long waves (inverse problem). This diagnostic method can also be employed to find the mass contents of individual components.

Physics and Astronomy Classification Scheme (PACS): 43.25.Zx, 43.25.Ba, 43.25.Ed, 43.25.Gf

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2020 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Vyacheslav Vakhnenko, Dmitri Vengrovich and Alexandre Michtchenko (March 12th 2020). Mathematical Fundamentals of a Diagnostic Method by Long Nonlinear Waves for the Structured Media, Advances in Complex Analysis and Applications, Francisco Bulnes and Olga Hachay, IntechOpen, DOI: 10.5772/intechopen.91462. Available from:

### Related Content

#### Advances in Complex Analysis and Applications

Edited by Francisco Bulnes

Next chapter

#### Introductory Chapter: Frontiers and Future Developments of the Complex Analysis

By Francisco Bulnes

#### Recent Advances in Integral Equations

Edited by Francisco Bulnes

First chapter

#### Introductory Chapter: Frontier Research on Integral Equations and Recent Results

By Francisco Bulnes

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

## Book Subject Areas

View all Books