On the Propagation of Longitudinal Stress Waves in Solids and Fluids by Unifying the Navier-Lame and Navier-Stokes Equations

Propagation of mechanical waves’ phenomenon is the result of infinitely small displacements of integrated individual particles in the materials. These displacements are governed by Navier-Lame and Navier-Stokes equations in solids and fluids, respectively. In the present work, a generalized Kelvin-Voigt model of viscoelasticity has been proposed with the aim of bridging the gap between solids and fluids leading to a new concept of viscoelasticity which unifies the Navier-Lame and the Navier-Stokes equations. On solving this equation in one dimension, propagation of stress disturbance in the so-called “Kelvin-Voigt materials” will be studied. The model of these materials involves all the elastic and viscoelastic solids, as well as fluids and soft materials.


Introduction
The subject matter of mechanics is the study of motion, in how a physical object changes position with time and why [1].In continuum mechanics we are concerned with the mechanical behavior and shape of materials under load.The physical reasons for this behavior can be quite different for different materials.Solid media will deform when forces are applied on them.These materials are called elastic, if the object will return to its initial shape and size when these forces are removed.Hence elasticity is the tendency of solid materials to return to their original shape after being deformed [2].In continuum mechanics we consider the basic equations describing the physical effects created by external forces acting upon solids and fluids.In addition to the basic equations that are applicable to all continua, there are other equations, called constitutive equations which are constructed to take into account material characteristics.In the study of solids the constitutive equations for a linear elastic material are a set of relations between stress and strain.Hook's law represents the material behavior and relates the Cauchy stress tensor  = (  ) and infinitesimal strain tensor  = (  ).The general form of Hook's law in components is where   is the fourth-order stiffness tensor [3].If the body is isotropic and homogenous, then this law is simply written as [3]  =  + 2 (2) in which  is the trace of ,  is the unit tensor, and  and  are the Lame constants.Perfect elasticity is an approximation of the real world and few materials remain purely elastic even after very small deformations.When an elastic material is not stressed in tension or compression beyond its elastic limit, its individual particles perform elastic movement.The displacement of the particles mass center, denoted by vector field  = ( 1 ,  2 ,  3 ), is related to the strain tensor by the relation [4] 2

Mathematical Problems in Engineering
Newton's second law of motion is based on the conservation of momentum and expressed as [5] div in which  is the mass density of the body.Now, in an isotropic, homogenous elastic solid, by combining (2), (3), and (4) we have which is the vector form of the Navier-Lame equations and governs the infinitesimal movements of the body's integrated individual particles [5].
Viscous materials resist shear flow and strain linearly with time, when stress is applied [6].If a material exhibits a linear response, it is categorized as a Newtonian material [6].In this case the stress is linearly proportional to the strain rate.If the material exhibits a nonlinear response to the strain rate, it is categorized as non-Newtonian fluid.To get the equation that governs the small movement of integrated particles of fluids, let us first consider the continuity equation, governing the continuity of the integrated individual particles which are derived from the principle of conservation of mass and is given as where the vector field  represents the velocity of the particles and  is the mass density of the fluid [7].In the absence of body forces, conservation of momentum, on which Newton's second law is based, is expressed as [7] div  =    .
In the study of fluids, the constitutive equations consist of a set of relations between stress and rate of strain; these equations in tensor form are given as [8] where the pressure  is induced by tension, that is, the difference between the dynamical and thermodynamical pressures,  and  are independent parameters characterizing viscosity,  is the unit tensor, and  is the rate of strain which is a symmetric tensor of order two.
Finally the kinematical relation of fluids describing the relation between tensor  and velocity vector field  is given as Moreover, the pressure can be written as [8]  = −( + 2 3 ) tr().(10) It can be deduced from ( 7), ( 8), (9), and (10) that However, in terms of displacement vector field , which is the set of Navier-Stokes equations.
In fluid continuum, the motion of substances is described by the Navier-Stokes equations.These equations are strictly the statement of conservation of momentum and are based on the assumption that the fluid, at the scale of interest, is a continuum.In other words, it is not made up of discrete particles but rather a continuous substance.Another necessary assumption is that all the fields of interest like pressure, velocity, density, and temperature are differentiable weakly at least [9].
There are materials for which a suddenly applied and maintained state of uniform stress induces an instantaneous deformation followed by a flow process which may or may not be limited in magnitude as time grows [10].These materials exhibit both solid and fluid characteristics.Behavior of these materials clearly cannot be described by either elasticity or viscosity theories alone, as it combines features of each and is called viscoelastic.Viscoelasticity is a generalization of elasticity and viscosity [11].
All materials exhibit some viscoelastic response.In common metals such as steel or aluminum, as well as in quartz, at room temperature, and at small strain, the behavior does not deviate much from linear elasticity, which is the simplest response of a viscoelastic material.Synthetic polymers, wood, and human tissue as well as metals at high temperature display significant viscoelastic effects.Some phenomena in viscoelastic materials are as follows.
(i) If the stress is held constant, the strain increases with time (creep).
(ii) If the strain is held constant, the stress decreases with time (relaxation).
The material creeps, that gives the prefix visco-, and the material fully recovers, which gives the suffix-elasticity [12].The viscoelastic materials play an important role in many engineering structures.Those materials such as polymers are being used, for example, to dissipate and to insulate vibration caused by rotating or reciprocal movements [13].They also have potential application in a new Hopkinson pressure bar testing apparatus [14][15][16][17][18]. Therefore, having the knowledge of the behavior of these materials, particularly related to their mechanical parameters, is essential.The theory which illustrates this behavior is the viscoelasticity theory.This theory is used in many fields, such as solid mechanics, seismology, exploration geophysics, acoustics, and engineering [19].
Modeling and model parameter estimation are of great importance for a correct prediction of the foundation behavior [19].In many cases, elastic constitutive models work well when time dependent effects can be neglected.However in those cases when time dependent effects cannot be neglected, we will need to utilize different constitutive models.Basically, time dependent effects indicate that the stress-strain behavior of material will change with time.The elastic material model for time dependent effects is viscoelasticity.Many researchers as Alfrey [20], Barberán and Herrera [21], Achenbach and Reddy [22], Bhattacharya and Sengupta [23], and Acharya et al. [24] formulated and developed this theory [19].Further, Bert and Egle [25], Abd-Alla and Ahmed [26], and Batra [27] successfully applied this theory to wave propagation in homogenous, elastic media [19].Murayama and Shibata [28] and Schiffman et al. [29] have proposed higher order viscoelastic models of five and seven parameters to represent the soil behavior.Jankowski et al. [30] discussed the linear viscoelastic model and the nonlinear viscoelastic model [19].
One of the most important classic models which is focused on, in this paper is Kelvin-Voigt model.The constitutive relation of this model is expressed as a linear first-order differential equation, which can be derived as below.
The most simple one-dimensional stress and strain tensors are of the form which, according to the generalized Hook's law, are related by where  is Young's modulus of elasticity [3].In the Kelvin-Voigt model for viscoelastic materials, ( 14) is augmented to include viscosity, leading to the generalized equation where  is the constant of viscosity [31,32].This model represents a solid undergoing reversible, viscoelastic strain and is extremely good with modeling creep in materials, but with regard to relaxation the model is much less accurate [33].
All material substances are comprised of particles.When a material is not stressed in tension or compression beyond its elastic limit, its individual particles may be forced into vibrational motion about their equilibrium positions.Thus, these particles perform elastic oscillations.Elastic waves are focused on particles that move in unison to produce a mechanical wave.In solids, elastic waves can propagate in four principle modes that are based on the way the particles oscillate.These waves can propagate as longitudinal, shear, and surface waves and in the thin materials as plate waves.
In longitudinal waves, the oscillations occur in the longitudinal direction or the direction of wave propagation.Since compressional and dilatational forces are active in these waves, they are also called pressure or compressional waves.They are also sometimes called density waves because their particle density fluctuates as they move.Compression waves can be generated in liquids, as well as solids, because the energy travels through the atomic structure by a series of compressions and expansion movements [34].In the transverse or shear wave, the particles oscillate at a right angle or perpendicular to the direction of propagation.Shear waves require a solid material for effective propagation and therefore are not effectively propagated in materials such as liquids or gasses [34].Waves in an isotropic elastic solid are governed by the vector Navier-Lame equations.

Generalization of Kelvin-Voigt Model.
In Kelvin-Voigt model for viscoelastic materials, the viscosity term is augmented to (14), leading to the generalized equation (15).But this model is a generalization of one-dimensional constitutive equation (14).In order to obtain a model for viscoelastic materials in most general form, the constitutive equation ( 2) should be extended.In this paper this equation is rewritten as in which   is called the first constant of viscosity.Also by symmetric of ( 16) this equation can be written as in which   is the second constant of viscosity.In the simple case ( 13), ( 16) is reduced as in which  is a positive constant.By comparing ( 18) and ( 15), it can be easily inferred that In the next sections, it has been shown that ( 16) and ( 17) can justify the physics.
Mathematical Problems in Engineering 2.3.Unifying Navier-Lame and Navier-Stokes Equations.At this point, by comparing ( 21) and ( 12), it is apparent that ( 21) can be written in two parts.The first and second parts are concerned with the solid and fluid properties of the body, respectively.In fact, (21) governs the displacements of the particles in both solids and fluids and in bodies with both solid and fluid properties, but not in materials involving two or more different phases.Equation ( 21) is a unification of Navier-Lame and Navier-Stokes equations; therefore, ( 16) and ( 17) can be used as constitutive equation to all viscoelastic materials.

Longitudinal Waves in Materials
Equation ( 21) governs displacement of integrated individual particles of any continuum body.To study the propagation of longitudinal stress waves in a body, we note that ( 21) can be written in one dimension as follows: in which where the only initial condition is (0, 0) = 0.This equation can be solved by separation of variables method, so (, ) = ()().Using this in ( 22) leads to in which  is constant.Now, the following five cases are considered.
Case 1 ( is positive).In this case one can assume  = 1/ 2 ,  ̸ = 0; hence, (24) becomes in which The solution, therefore, can be written as By applying the initial condition, one arrives at If  2 = − 1 , then the velocity of the particle, placed at the origin, is zero, since Therefore, in this case, there is no motion in the medium.
On the other hand, since the medium is assumed to be elastic, the movement of the particle, located at the origin, stops after a while.Hence, if  2 = − 1 , at time , then (29) can be written as Since  1 is positive and  2 is negative, the second part of ( 30) is not zero and so  2 = − 1 .This means that there is no motion in the medium.Therefore, in this case, the equation does not have a nontrivial solution.
Case 2. If  = 0, then the solution is trivial.
Case 3. If  is negative and  2  4 − 4 2 is positive, then one can assume  = −1/ 2 ,  ̸ = 0, so in which Thus the solution is as follows: Also, the velocity is Now, from the initial condition, one obtains If  1 = 0, then   (0, ) = 0, which means that there is no motion and the solution is trivial.Otherwise,  2 = − 1 , which yields Since after a time  the particle that is at the origin stops, then Consequently, the solution can be written as On the other hand, it can be written as follows: Case 4 ( is negative and  2  4 − 4 2 = 0).In this case the solution is trivial.
Case 5 ( is negative and  2  4 − 4 2 < 0).In this case one can assume  = −1/ 2 ,  ̸ = 0, so that in which Hence, the solution is Also, the velocity is given as From the initial condition one obtains If  1 = 0, then   (0, ) = 0, and hence there is no motion in the medium and the solution is trivial.Therefore,  1 ̸ = 0, and  1 must be zero, so that Since after a time  the particle which is at the origin stops, then Consequently the solution can be written as in which Therefore, the solution can be written as (, ) =  −( 2  2 )/2 cos( − ) sin(). (49)

4.1.
Longitudinal Waves in a Perfectly Elastic Solid.In a metal like aluminum, whose characteristics are listed in Table 1, the wave that propagates in the medium is of the form (49).
If one considers  = 2,  = 0.003, and  = 1, the depiction of the longitudinal wave that propagates from the origin is shown in Figure 1.It is clear that it is a wave with no attenuation.

Longitudinal Waves in Viscoelastic
Solids.If the viscoelastic properties are characterized as in Table 2, the depiction of the longitudinal stress wave is shown in Figure 2. It is also clear that this wave is attenuated.

Longitudinal Waves in a Viscoelastic Newtonian Fluid.
In a liquid with parameters shown in Table 3, the wave propagating in the medium is of the form (39).So then the wave characteristics are discussed in the following for material with properties as mentioned in Table 3.
If one assumes  = 2,  = 0.01, and  = 0.003, then the displacement graph of particles in this medium is as shown in Figure 3.In Figure 3(a) the displacement of particles at the origin without returning to their initial position and in Figure 3(b) the displacement of a particle in each position  without returning to its initial position are seen.4, the graph of longitudinal stress wave in this medium is as shown in Figure 4.In Figure 4(a), the distance between the particle originally located in the origin and the same particle moved away from the origin in time  is seen.This distance starts to decrease with the particles returning to the origin.In Figure 4(b), the distance between the particle at the time  and the same particle approaching its equilibrium position at the origin is shown; in Figure 4(c) the distance between the particles located in their respective equilibrium position and the same particles moved away from the equilibrium position in time  is shown.This distance then starts to decrease with the particles returning to their respective equilibrium position.

Conclusions
When a material is subjected to an impulsive force, there occur infinitesimal displacements of the integrated particles of that body.In solids and in fluids the Navier-Lame and the Navier-Stokes equations govern these displacements, respectively.In the present work it was shown that through the use of a generalized form of Kelvin-Voigt model of viscoelasticity the unification of these two equations is possible and one obtains an equation that can represent solid materials, fluids, and soft materials.A new concept of viscoelasticity was defined in which every material has some degree of both solid and fluid properties depending on the particular parameters of viscosity.This is of particular importance in relation to the soft materials.Using this generalized equation, propagation of stress disturbance, pulse, and waves in solids as well as fluids can be studied in one single frame.

Figure 1 :
Figure 1: The vibration of particles.

Figure 2 :
Figure 2: Vibration of a particle with attenuation.

Figure 3 :
Figure 3: The displacement of particles. u

Figure 4 :
Figure 4: The distance between the particles at the time .

Table 1 :
Parameters of aluminum that is considered to be perfectly elastic.

Table 2 :
Parameters of aluminum with fluid properties pertaining to viscoelasticity.

Table 3 :
Parameters of a Newtonian liquid like water which is perfectly fluid with no solid parameters.

Table 4 :
Parameters of a liquid like water which is not perfectly fluid and has solid parameters.