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

Authors: Ahmad Barzkar; Hojatollah Adibi
Journal: Mathematical Problems in Engineering
 (2015)

---

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


---

## Body

## 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 T
=
(
σ
i
j
) and infinitesimal strain tensor e
=
(
e
i
j
). The general form of Hook’s law in components is(1)
σ
i
j
=
C
i
j
l
k
e
l
k
,where C
i
j
l
k is the fourth-order stiffness tensor [3]. If the body is isotropic and homogenous, then this law is simply written as [3](2)
T
=
λ
θ
I
+
2
μ
ein which θ is the trace of e, I 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 u
=
(
u
1
,
u
2
,
u
3
), is related to the strain tensor by the relation [4](3)
e
=
1
2
grad
⁡
u
+
grad
⁡
u
T
.Newton’s second law of motion is based on the conservation of momentum and expressed as [5](4)
div
⁡
T
=
ρ
∂
2
u
∂
t
2in which ρ is the mass density of the body. Now, in an isotropic, homogenous elastic solid, by combining (2), (3), and (4) we have(5)
u
t
t
=
λ
+
μ
ρ
grad
⁡
div
⁡
u
+
μ
ρ
∇
2
uwhich 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(6)
∂
ρ
∂
t
+
div
⁡
ρ
V
=
0
,where the vector field V 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](7)
div
⁡
T
=
ρ
∂
V
∂
t
.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](8)
T
=
-
p
I
+
λ
tr
⁡
D
I
+
2
μ
D
,where the pressure p is induced by tension, that is, the difference between the dynamical and thermodynamical pressures, λ and μ are independent parameters characterizing viscosity, I is the unit tensor, and D is the rate of strain which is a symmetric tensor of order two.Finally the kinematical relation of fluids describing the relation between tensorD and velocity vector field V is given as(9)
D
=
1
2
grad
⁡
V
+
grad
⁡
V
T
.Moreover, the pressure can be written as [8](10)
p
=
-
λ
+
2
3
μ
tr
⁡
⁡
D
.It can be deduced from (7), (8), (9), and (10) that(11)
∂
V
∂
t
=
λ
+
μ
ρ
grad
⁡
div
⁡
V
+
μ
ρ
∇
2
V
.However, in terms of displacement vector field u,(12)
u
t
t
=
λ
+
μ
ρ
grad
⁡
div
⁡
u
˙
+
μ
ρ
∇
2
u
˙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 form(13)
T
=
σ
0
0
0
0
0
0
0
0
,
e
=
e
1
0
0
0
e
2
0
0
0
e
3which, according to the generalized Hook’s law, are related by(14)
σ
=
E
e
1
,where E 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(15)
σ
=
E
e
1
+
η
∂
e
1
∂
t
,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.


## 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 as(16)
T
=
λ
θ
I
+
2
μ
e
+
λ
f
θ
˙
Iin which λ
f is called the first constant of viscosity. Also by symmetric of (16) this equation can be written as(17)
T
=
λ
θ
I
+
2
μ
e
+
λ
f
θ
˙
I
+
2
μ
f
e
˙in which μ
f is the second constant of viscosity. In the simple case (13), (16) is reduced as(18)
σ
=
E
e
1
+
μ
λ
f
λ
+
μ
e
1
˙
-
μ
α
λ
+
μ
exp
⁡
⁡
-
λ
+
μ
λ
f
tin which α is a positive constant. By comparing (18) and (15), it can be easily inferred that(19)
η
=
μ
λ
f
λ
+
μ
.In 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:(20)
u
t
t
=
λ
+
μ
ρ
grad
⁡
div
⁡
u
+
μ
ρ
∇
2
u
+
λ
f
ρ
grad
⁡
div
⁡
u
˙
,where u
˙
=
∂
u
/
∂
t. 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: (21)
u
t
t
=
λ
+
μ
ρ
grad
⁡
div
⁡
u
+
μ
ρ
∇
2
u
+
λ
f
+
μ
f
ρ
grad
⁡
div
⁡
u
˙
+
μ
f
ρ
∇
2
u
˙
,where λ
f and μ
f 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.


## 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 as(16)
T
=
λ
θ
I
+
2
μ
e
+
λ
f
θ
˙
Iin which λ
f is called the first constant of viscosity. Also by symmetric of (16) this equation can be written as(17)
T
=
λ
θ
I
+
2
μ
e
+
λ
f
θ
˙
I
+
2
μ
f
e
˙in which μ
f is the second constant of viscosity. In the simple case (13), (16) is reduced as(18)
σ
=
E
e
1
+
μ
λ
f
λ
+
μ
e
1
˙
-
μ
α
λ
+
μ
exp
⁡
⁡
-
λ
+
μ
λ
f
tin which α is a positive constant. By comparing (18) and (15), it can be easily inferred that(19)
η
=
μ
λ
f
λ
+
μ
.In 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:(20)
u
t
t
=
λ
+
μ
ρ
grad
⁡
div
⁡
u
+
μ
ρ
∇
2
u
+
λ
f
ρ
grad
⁡
div
⁡
u
˙
,where u
˙
=
∂
u
/
∂
t. 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: (21)
u
t
t
=
λ
+
μ
ρ
grad
⁡
div
⁡
u
+
μ
ρ
∇
2
u
+
λ
f
+
μ
f
ρ
grad
⁡
div
⁡
u
˙
+
μ
f
ρ
∇
2
u
˙
,where λ
f and μ
f 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: (22)
u
t
t
=
c
2
u
x
x
+
d
2
u
t
x
x
,
∞
<
x
<
+
∞
,
t
≥
0in which(23)
c
2
=
λ
s
+
2
μ
s
ρ
,
d
2
=
λ
f
+
2
μ
f
ρ
,where the only initial condition is u
0,0
=
0. This equation can be solved by separation of variables method, so u
x
,
t
=
X
x
T
(
t
). Using this in (22) leads to(24)
X
X
′′
=
c
2
T
T
′′
+
d
2
T
′
T
′′
=
γin which γ is constant. Now, the following five cases are considered.Case 1 (γ is positive).
In this case one can assumeγ
=
1
/
k
2, k
≠
0; hence, (24) becomes(25)
X
′′
-
k
2
X
=
0
T
′′
-
k
2
d
2
T
′
-
k
2
c
2
T
=
0
→
yields
X
x
=
c
1
e
k
x
+
c
2
e
-
k
x
T
t
=
d
1
e
r
1
t
+
d
2
e
r
2
tin which(26)
r
1
=
1
2
k
2
d
2
+
k
k
2
d
4
+
4
c
2
,
r
2
=
1
2
k
2
d
2
-
k
k
2
d
4
+
4
c
2
.The solution, therefore, can be written as(27)
u
x
,
t
=
c
1
e
k
x
+
c
2
e
-
k
x
d
1
e
r
1
t
+
d
2
e
r
2
t
.By applying the initial condition, one arrives at(28)
u
0,0
=
c
1
+
c
2
d
1
+
d
2
=
0
⟶
c
2
=
-
c
1
or
d
2
=
-
d
1
.If c
2
=
-
c
1, then the velocity of the particle, placed at the origin, is zero, since(29)
u
t
0
,
t
=
c
1
+
c
2
d
1
r
1
e
r
1
t
+
d
2
r
2
e
r
2
t
.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, ifd
2
=
-
d
1, at time T, then (29) can be written as(30)
u
t
0
,
T
=
d
1
c
1
+
c
2
r
1
e
r
1
T
-
r
2
e
r
2
T
=
0
.Since r
1 is positive and r
2 is negative, the second part of (30) is not zero and so c
2
=
-
c
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 k
2
d
4
-
4
c
2 is positive, then one can assume γ
=
-
1
/
k
2, k
≠
0, so(31)
X
′′
+
k
2
X
=
0
T
′′
+
k
2
d
2
T
′
+
k
2
c
2
T
=
0
→
yields
X
x
=
c
1
cos
⁡
⁡
k
x
+
c
2
sin
⁡
k
x
T
t
=
d
1
e
r
1
t
+
d
2
e
r
2
tin which(32)
r
1
=
1
2
-
k
2
d
2
+
k
k
2
d
4
-
4
c
2
,
r
2
=
1
2
-
k
2
d
2
-
k
k
2
d
4
-
4
c
2
.Thus the solution is as follows:(33)
u
x
,
t
=
c
1
cos
⁡
⁡
k
x
+
c
2
sin
⁡
k
x
d
1
e
r
1
t
+
d
2
e
r
2
t
.Also, the velocity is(34)
u
t
x
,
t
=
c
1
cos
⁡
⁡
k
x
+
c
2
sin
⁡
k
x
r
1
d
1
e
r
1
t
+
r
2
d
2
e
r
2
t
.Now, from the initial condition, one obtains(35)
u
0,0
=
c
1
d
1
+
d
2
=
0
→
yields
c
1
=
0
,
or
d
2
=
-
d
1
.If c
1
=
0, then u
t
0
,
t
=
0, which means that there is no motion and the solution is trivial. Otherwise, d
2
=
-
d
1, which yields(36)
u
t
0
,
t
=
c
1
d
1
r
1
e
r
1
t
-
r
2
e
r
2
t
.Since after a time T the particle that is at the origin stops, then(37)
r
1
e
r
1
T
-
r
2
e
r
2
T
=
0
→
yields
T
=
1
r
1
-
r
2
ln
⁡
⁡
r
2
r
1
.Consequently, the solution can be written as(38)
u
x
,
t
=
d
1
c
1
cos
⁡
⁡
k
x
+
c
2
sin
⁡
k
x
e
r
1
t
-
e
r
2
t
.On the other hand, it can be written as follows:(39)
u
x
,
t
=
α
1
cos
⁡
⁡
k
x
+
α
2
sin
⁡
k
x
e
r
1
t
-
e
r
2
t
.Case 4 (γ is negative and k
2
d
4
-
4
c
2
=
0).
In this case the solution is trivial.Case 5 (γ is negative and k
2
d
4
-
4
c
2
<
0).
In this case one can assumeγ
=
-
1
/
k
2, k
≠
0, so that(40)
X
′′
+
k
2
X
=
0
T
′′
+
k
2
d
2
T
′
+
k
2
c
2
T
=
0
→
yields
X
x
=
c
1
cos
⁡
⁡
k
x
+
c
2
sin
⁡
k
x
T
t
=
e
-
(
k
2
d
2
t
)
/
2
d
1
cos
⁡
⁡
r
t
+
d
2
sin
⁡
r
tin which(41)
r
=
1
2
k
4
c
2
-
k
2
d
4
.Hence, the solution is(42)
u
x
,
t
=
e
-
(
k
2
d
2
t
)
/
2
c
1
cos
⁡
⁡
k
x
+
c
2
sin
⁡
k
x
×
d
1
cos
⁡
⁡
r
t
+
d
2
sin
⁡
r
t
.Also, the velocity is given as(43)
u
t
x
,
t
=
e
-
(
k
2
d
2
t
)
/
2
c
1
cos
⁡
⁡
k
x
+
c
2
sin
⁡
k
x
×
r
d
2
-
1
2
k
2
d
2
d
1
cos
⁡
⁡
r
t
-
r
d
1
+
1
2
k
2
d
2
d
2
sin
⁡
r
t
.From the initial condition one obtains(44)
u
0,0
=
c
1
d
1
=
0
→
yields
c
1
=
0
,
or
d
1
=
0
.If c
1
=
0, then u
t
0
,
t
=
0, and hence there is no motion in the medium and the solution is trivial. Therefore, c
1
≠
0, and d
1 must be zero, so that(45)
u
t
0
,
t
=
e
-
k
2
d
2
t
/
2
c
1
d
2
r
cos
⁡
⁡
r
t
-
1
2
k
2
d
2
sin
⁡
r
t
.Since after a time T the particle which is at the origin stops, then(46)
r
cos
⁡
⁡
r
T
-
1
2
k
2
d
2
sin
⁡
r
T
=
0
⟶
T
n
=
1
r
tan
-
1
⁡
2
r
k
2
d
2
+
n
π
.Consequently the solution can be written as(47)
u
x
,
t
=
2
c
1
d
2
c
1
2
+
c
2
2
×
e
-
k
2
d
2
t
/
2
cos
⁡
⁡
k
x
-
φ
sin
⁡
r
tin which(48)
φ
=
cos
⁡
-
1
⁡
c
1
c
1
2
+
c
2
2
.Therefore, the solution can be written as(49)
u
x
,
t
=
β
e
-
(
k
2
d
2
t
)
/
2
cos
⁡
⁡
k
x
-
φ
sin
⁡
r
t
.


## 4. Results

### 4.1. Longitudinal Waves in a Perfectly Elastic Solid

In a metal like aluminum, whose characteristics are listed in Table1, the wave that propagates in the medium is of the form (49).Table 1
Parameters of aluminum that is considered to be perfectly elastic.
λ
s
λ
f
μ
s
μ
f
ρ
5
0.6
×
10
6 KN/m2
0
27.6
×
10
6 KN/m2
0
56000 Kg/m3If one considersk
=
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.Figure 1
The vibration of particles.
(a)
Particle at the origin
(b)
Particle at each positionx


### 4.2. Longitudinal Waves in Viscoelastic Solids

If the viscoelastic properties are characterized as in Table2, the depiction of the longitudinal stress wave is shown in Figure 2. It is also clear that this wave is attenuated.Table 2
Parameters of aluminum with fluid properties pertaining to viscoelasticity.
λ
s
λ
f
μ
s
μ
f
ρ
50.6 × 106 KN/m2
3016
27.6 × 106 KN/m2
2130
56000 Kg/m3Figure 2
Vibration of a particle with attenuation.
(a)
Particle at the origin
(b)
Particle at each positionx


### 4.3. Longitudinal Waves in a Viscoelastic Newtonian Fluid

In a liquid with parameters shown in Table3, 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.Table 3
Parameters of a Newtonian liquid like water which is perfectly fluid with no solid parameters.
λ
s
λ
f
μ
s
μ
f
ρ
0
30.12 × 102
0
21.3 × 102
996.3 Kg/m3If one assumesk
=
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 x without returning to its initial position are seen.Figure 3
The displacement of particles.
(a)
Particle at the origin
(b)
Particle at each positionx


### 4.4. Longitudinal Waves in an Ideal Newtonian Fluid

If a liquid is viscoelastic with viscosity coefficients shown in Table4, 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 T 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 T 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 T is shown. This distance then starts to decrease with the particles returning to their respective equilibrium position.Table 4
Parameters of a liquid like water which is not perfectly fluid and has solid parameters.
λ
s
λ
f
μ
s
μ
f
ρ
50.6
30.12 × 102
27.6
21.3 × 102
996.3 Kg/m3Figure 4
The distance between the particles at the timeT.
(a)
The distance between the particle originally located in the origin and the same particle moved away from the origin in timeT
(b)
The distance between the particle at the timeT and the same particle approaching its equilibrium position at the origin
(c)
The distance between the particles located in their respective equilibrium position and the same particles moved away from the equilibrium position in timeT


## 4.1. Longitudinal Waves in a Perfectly Elastic Solid

In a metal like aluminum, whose characteristics are listed in Table1, the wave that propagates in the medium is of the form (49).Table 1
Parameters of aluminum that is considered to be perfectly elastic.
λ
s
λ
f
μ
s
μ
f
ρ
5
0.6
×
10
6 KN/m2
0
27.6
×
10
6 KN/m2
0
56000 Kg/m3If one considersk
=
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.Figure 1
The vibration of particles.
(a)
Particle at the origin
(b)
Particle at each positionx


## 4.2. Longitudinal Waves in Viscoelastic Solids

If the viscoelastic properties are characterized as in Table2, the depiction of the longitudinal stress wave is shown in Figure 2. It is also clear that this wave is attenuated.Table 2
Parameters of aluminum with fluid properties pertaining to viscoelasticity.
λ
s
λ
f
μ
s
μ
f
ρ
50.6 × 106 KN/m2
3016
27.6 × 106 KN/m2
2130
56000 Kg/m3Figure 2
Vibration of a particle with attenuation.
(a)
Particle at the origin
(b)
Particle at each positionx


## 4.3. Longitudinal Waves in a Viscoelastic Newtonian Fluid

In a liquid with parameters shown in Table3, 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.Table 3
Parameters of a Newtonian liquid like water which is perfectly fluid with no solid parameters.
λ
s
λ
f
μ
s
μ
f
ρ
0
30.12 × 102
0
21.3 × 102
996.3 Kg/m3If one assumesk
=
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 x without returning to its initial position are seen.Figure 3
The displacement of particles.
(a)
Particle at the origin
(b)
Particle at each positionx


## 4.4. Longitudinal Waves in an Ideal Newtonian Fluid

If a liquid is viscoelastic with viscosity coefficients shown in Table4, 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 T 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 T 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 T is shown. This distance then starts to decrease with the particles returning to their respective equilibrium position.Table 4
Parameters of a liquid like water which is not perfectly fluid and has solid parameters.
λ
s
λ
f
μ
s
μ
f
ρ
50.6
30.12 × 102
27.6
21.3 × 102
996.3 Kg/m3Figure 4
The distance between the particles at the timeT.
(a)
The distance between the particle originally located in the origin and the same particle moved away from the origin in timeT
(b)
The distance between the particle at the timeT and the same particle approaching its equilibrium position at the origin
(c)
The distance between the particles located in their respective equilibrium position and the same particles moved away from the equilibrium position in timeT


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



---

*Source: 789238-2015-04-27.xml*