Mathematical Problems in Engineering

Volume 2015 (2015), Article ID 789238, 9 pages

http://dx.doi.org/10.1155/2015/789238

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

Department of Applied Mathematics, Faculty of Mathematics and Computer Sciences, Amirkabir University of Technology, 424 Hafez Avenue, P.O. Box 15875-4413, Tehran, Iran

Received 2 September 2014; Accepted 14 December 2014

Academic Editor: Florin Pop

Copyright © 2015 Ahmad Barzkar and Hojatollah Adibi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

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.

#### 1. 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 iswhere is the fourth-order stiffness tensor [3]. If the body is isotropic and homogenous, then this law is simply written as [3]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 , is related to the strain tensor by the relation [4]Newton’s second law of motion is based on the conservation of momentum and expressed as [5]in which is the mass density of the body. Now, in an isotropic, homogenous elastic solid, by combining (2), (3), and (4) we havewhich 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 aswhere 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]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 asMoreover, the pressure can be written as [8]It can be deduced from (7), (8), (9), and (10) thatHowever, 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).(iii)Acoustic waves experience attenuation.

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–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 formwhich, according to the generalized Hook’s law, are related bywhere 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 equationwhere 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.

#### 2. New Mathematical Model

##### 2.1. 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 asin which is called the first constant of viscosity. Also by symmetric of (16) this equation can be written asin which is the second constant of viscosity. In the simple case (13), (16) is reduced asin which is a positive constant. By comparing (18) and (15), it can be easily inferred thatIn the next sections, it has been shown that (16) and (17) can justify the physics.

##### 2.2. Generalization of Navier-Lame Equation

In the first section, the Navier-Lame equation was deduced by combining (2), (3), and (4). Now (2) is replaced by (16); therefore, by combining (16), (3), and (4), a generalized Navier-Lame equation is obtained as follows:where . Also if (2) is replaced by (17), then, by combining (17), (3), and (4), another generalization of Navier-Lame equation can be deduced as follows: where and resemble the form of equations for and , respectively.

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

#### 3. 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 whichwhere the only initial condition is . This equation can be solved by separation of variables method, so . Using this in (22) leads toin which is constant. Now, the following five cases are considered.

*Case 1 ( is positive). *In this case one can assume , ; hence, (24) becomesin whichThe solution, therefore, can be written asBy applying the initial condition, one arrives atIf , then the velocity of the particle, placed at the origin, is zero, sinceTherefore, 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 , at time , then (29) can be written asSince is positive and is negative, the second part of (30) is not zero and so . 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 , then the solution is trivial.

*Case 3. *If is negative and is positive, then one can assume , , soin whichThus the solution is as follows:Also, the velocity isNow, from the initial condition, one obtainsIf , then , which means that there is no motion and the solution is trivial. Otherwise, , which yieldsSince after a time the particle that is at the origin stops, thenConsequently, the solution can be written asOn the other hand, it can be written as follows:

*Case 4 ( is negative and ). *In this case the solution is trivial.

*Case 5 ( is negative and ). *In this case one can assume , , so thatin whichHence, the solution isAlso, the velocity is given asFrom the initial condition one obtainsIf , then , and hence there is no motion in the medium and the solution is trivial. Therefore, , and must be zero, so thatSince after a time the particle which is at the origin stops, thenConsequently the solution can be written asin whichTherefore, the solution can be written as

#### 4. Results

##### 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).