# Shock Wave Solution for a Nonlinear Partial Differential Equation Arising in the Study of a Non-Newtonian Fourth Grade Fluid Model

Authors: Taha Aziz; A. Fatima; F. M. Mahomed
Journal: Mathematical Problems in Engineering
 (2013)

---

## Abstract

This study focuses on obtaining a new class of closed-form shock wave solution
also known as soliton solution for a nonlinear partial differential equation which governs
the unsteady magnetohydrodynamics (MHD) flow of an incompressible fourth grade fluid model. The travelling wave symmetry formulation of the model leads to a shock wave
solution of the problem. The restriction on the physical parameters of the flow problem also
falls out naturally in the course of derivation of the solution.


---

## Body

## 1. Introduction

A shock wave (also called shock front or simply “shock”) is a type of propagating disturbance through a media. It actually expresses a sharp discontinuity of the parameters that delineate the media. Unlike solitons, the energy is a conserved quantity and thus remains constant during its propagation. Shock wave dissipates energy relatively quickly with distance. One source of a shock wave is when the supersonic jets fly at a speed that is greater than the speed of sound. This results in the drag force on aircraft with shocks. These waves also appear in various interesting phenomena in real-life situations. For example, solitons appear in the propagation of pulses through optical fibers. Another example is where cnoidal waves appear in shallow water waves although an extremely scarce phenomena.The dynamics and mechanics of non-Newtonain fluid flow problems have been an important area of interest in the recent few years. The flow phenomena of non-Newtonian fluids occur in a variety of industrial and technological applications. Because of the diverse physical structure and behavior of non-Newtonian fluids, there is no single mathematical expression which describes all the characteristics of non-Newtonian fluids. Due to this fact, several models of non-Newtonian fluids have been proposed. Apart from this fact, the model equations of the problem dealing with the flow of non-Newtonian fluids are higher-order nonlinear and complex in nature. Several methods have been developed in the recent years to obtain the solutions of these sort of flow problems. Some of these techniques are variational iteration method, the Adomian decomposition method, homotopy perturbation method, homotopy analysis method, and semi-inverse variational method. But all these techniques fail to develop the exact (closed-form) solutions of the non-Newtonian fluid flow problems.One of the simplest classes of non-Newtonian fluid models is the second grade fluid [1]. Although the second grade model is found to predict the normal stress differences, it does not take into account the shear thinning and shear thickening phenomena due to its constant apparent shear viscosity. For this reason, some experiments may be well described through fluids of grade three or four. Very little attention has been given to date to the flows of fourth grade fluid [2]. This model is known as the most generalized model amongst the differential-type non-Newtonian fluid models [3]. The fourth grade fluid model describes most of the non-Newtonian flow properties at one time. This model is known to capture the interesting non-Newtonian flow properties such as shear thinning and shear thickening that many other non-Newtonian models do not show. This model is also capable of predicting the normal stress effects that lead to phenomena like “die-swell” and “rod-climbing” [4]. With these facts in mind, we have considered a fourth grade fluid model in this study. In general, the model equations of the problem dealing with the flow of fourth grade fluids are higher-order nonlinear equations. The literature survey witnesses that very limited studies are reported in the literature, up to now, dealing with the flow problems of fourth grade fluid and these investigations that  further narrow are down when we discussed the closed-form solutions of these problems. However, some useful and interesting communications in this direction are made in the studies [5–11].In this study, we have used an interesting method to construct the solution of nonlinear problem arising in the study of non-Newtonian fluid. We have explored the shock wave behavior of the problem which deals with the unsteady flow of fourth grade fluid. We have also taken into account the magnetohydrodynamic nature of the fluid by applying uniform magnetic field as an external body force. This concept is introduced so that our solution can be easily reduced to the problem that deals with the effects of body forces.


## 2. Mathematical Structure of the Model

The unsteady MHD flow of an incompressible fluid is governed by law of conservation of mass and momentum; namely,(1)divV=0,(2)ρdVdt=divT-σB02V.
In the above equations, V is the velocity vector, ρ the density of the fluid, d/dt the total time derivative, and T the Cauchy stress tensor. We have considered a uniform magnetic field of strength B0, which is applied in the transverse direction of the flow as an external body force by assuming that the induced magnetic field and the external field are negligible.For fourth grade fluid model, the Cauchy stress tensor satisfies the constitutive equations [3]:
(3)T=-pI+∑j=1n‍Sjwithn=4,
where p is the pressure, I the identity tensor, and Sj the extra stress tensor as(4)S1=μA1,S2=α1A2+α2A12,S3=β1A3+β2(A1A2+A2A1)+β3(trA12)A1,S4=γ1A4+γ2(A3A1+A1A3)+γ3A22+γ4(A2A12+A12A2)+γ5(trA2)A2+γ6(trA2)A12+[γ7trA3+γ8tr(A2A1)]A1.
Here, μ is the dynamic viscosity; αi(i=1,2),βi(i=1,2,3), and γi(i=1,2,...,8) are material constants. The Rivlin-Ericksen tensors A1 to A4 are defined by
(5)A1=(gradV)+(gradV)T,(6)An=dAn-1dt+An-1(gradV)+(gradV)TAn-1(n>1),
in which grad is the gradient operator.Consider the unsteady MHD flow of an incompressible fourth grade fluid which occupies the half-spacey>0 over an infinite rigid plate which lies in the xz-plane. The x-axis and y-axis are chosen parallel and perpendicular to the plate. By taking the velocity field (u(y,t),0,0), the conservation of mass equation is identically satisfied. To obtain the governing PDE in u, substituting (3)–(6) into (2) and rearranging, we obtain the following model equation in the absence of the modified pressure gradient:
(7)ρ∂u∂t=μ∂2u∂y2+α1∂3u∂y2∂t+β1∂4u∂y2∂t2+6(β2+β3)(∂u∂y)2∂2u∂y2+γ1∂5u∂y2∂t3+2(3γ2+γ3+γ4+γ5+3γ7+γ8)∂∂y×[(∂u∂y)2∂2u∂y∂t]-σB02u.


## 3. Reduction of the Model Equation

We know that from the principal of the Lie symmetry methods that if a differential equation is explicitly independent of any dependent or independent variable, then this particular differential equation remains invariant under the translation symmetry corresponding to that particular variable. We noticed that (7) admits the Lie point symmetry generators, ∂/∂t (time-translation) and ∂/∂y (space-translation in y). For a detailed analysis, the readers are referred to [12, 13].LetX1 and X2 be time-translation and space-translation symmetry generators, respectively. Then, the solution corresponding to the generator
(8)X=X1+mX2=∂∂t+m∂∂y(m>0)
would represent the travelling wave solution with constant wave speed m. The Langrangian system corresponding to (8) is
(9)dym=dt1=du0.
Solving (9), invariant solutions are given by
(10)u(y,t)=f(η)withη=y-mt,
where f(η) is an arbitrary function of the characteristic variable η=y-mt. Making use of  (10) into (7) results in a fifth-order ordinary differential for f(η)(11)-mρdfdη=μd2fdη2-α1md3fdη3+β1m2d4fdη4+6(β2+β3)(dfdη)2d2fdη2-γ1m3d5fdη5-2m(3γ2+γ3+γ4+γ5+3γ7+γ8)ddη×[(dfdη)2d2fdη2]-σB02f.
Thus, the original fifth-order nonlinear PDE (7) reduced to a fifth-order ODE (11) along certain curves in the y-t plane. These curves are called characteristic curves or just the characteristic.


## 4. Shock Wave Solution

Now, we obtain shock wave solution of the reduced equation (11). The starting hypothesis for shock wave solution is given by [14–17]
(12)f(η)=Aexp(Bη),
where A and B are the free parameters to be determined. Inserting (12) in (11), we obtain
(13)0=[mρB+μB2-mα1B3+β1m2B4-γ1m3B5-σB02]+e2Bη[6(β2+β3)A2B4-6m(3γ2+γ3+γ4+γ5+3γ7+γ8)A2B5].
Separating (13) in the powers of e0 and e2Bη, we find the following:
(14)e0:mρB+μB2-mα1B3+β1m2B4-γ1m3B5-σB02=0,(15)e2Bη:6(β2+β3)A2B4-6m(3γ2+γ3+γ4+γ5+3γ7+γ8)A2B5=0.
From (15), we deduce
(16)B=(β2+β3)m(3γ2+γ3+γ4+γ5+3γ7+γ8).
Using the value of B in (14), we obtain
(17)0=ρ(β2+β3)(3γ2+γ3+γ4+γ5+3γ7+γ8)+μ(β2+β3)2m2(3γ2+γ3+γ4+γ5+3γ7+γ8)2-α1(β2+β3)3m2(3γ2+γ3+γ4+γ5+3γ7+γ8)3+β1(β2+β3)4m2(3γ2+γ3+γ4+γ5+3γ7+γ8)4-γ1(β2+β3)5m2(3γ2+γ3+γ4+γ5+3γ7+γ8)5-σB02.
Thus, the solution for f(η) (provided the condition (17) holds) can be written as
(18)f(η)=Aexp[(β2+β3)m(3γ2+γ3+γ4+γ5+3γ7+γ8)η].
So, the solution u(y,t) which satisfies the condition (17) is written as
(19)u(y,t)=Aexp[(β2+β3)(y-mt)m(3γ2+γ3+γ4+γ5+3γ7+γ8)]u(y,t)=Aexpu(y,t)=Aexpu(y,t)withm>0.
We observe that the solution (19) does satisfy the physically relevant boundary and initial conditions(20a)u(0,t)=E(t),(20b)u(y,0)=F(y),(20c)∂u(y,0)∂t=G(y),(20d)∂2u(y,0)∂t2=H(y),where(21a)E(t)=Aexp(-(β2+β3)mtm(3γ2+γ3+γ4+γ5+3γ7+γ8)),(21b)F(y)=Aexp((β2+β3)ym(3γ2+γ3+γ4+γ5+3γ7+γ8)),(21c)G(y)=-(β2+β3)(3γ2+γ3+γ4+γ5+3γ7+γ8)F(y),(21d)H(y)=((β2+β3)(3γ2+γ3+γ4+γ5+3γ7+γ8))2F(y),with
(22)A=E(0)=F(0).
The functions E(t), F(y), G(y), and H(y) depend on the physical parameters of the flow model. The boundary condition (20a) is the no-slip condition at y=0. The initial velocity E(0) of the rigid plate can be prescribed, but its velocity for t>0, E(t), cannot be arbitrary and is given by (21a). Similarly, the initial velocity profile u(y,0) cannot be arbitrary and is given by (21b).Note that the solution (19) is the soliton solution or the shock wave solution to the governing PDE (7). The above solution is valid under the particular condition on the physical parameters of the flow model given in (17). This solution does show the hidden shock wave behavior of the flow problem with slope of the velocity field or the velocity gradient approaches to infinity such that
(23)∂u∂y⟶∞asy>0.
If we denote
(24)ϵ=(β2+β3)(3γ2+γ3+γ4+γ5+3γ7+γ8),
thus, the imposing condition (17) can be written as
(25)m2=[-γ1ϵ5+β1ϵ4-α1ϵ3+μϵ2](σB02-ρϵ),withσB02≠ρϵ.
We observe that the condition (25) gives the speed m of the travelling shock wave. The range of the values of ϵ for which m is real depends not only on the zeros of the cubic polynomial in ϵ on the numerator of (25) but also on the sign of the denominator. The shock wave behavior of solution (19) is observed from Figure 1.Figure 1
Shock wave behaviour of solution (19), where β2=2,β3=1,m=1,t=π/2, and γ2=γ3=γ4=γ5=γ7=γ8=0.2.


## 5. Concluding Remarks

In this study, we have obtained the mathematical structure of the closed-form shock wave solution of higher-order nonlinear PDE arising in the study of a fourth grade non-Newtonian fluid model. Translational symmetry generators in variablest and y have been utilized to perform reduction of governing nonlinear partial differential equation into ordinary differential equation, and, thereafter, the closed-form shock wave solution has been constructed. We have considered a prototype model of the flow problem, but the solution is going to be very helpful in carrying out further analysis of the shock wave behavior associated with the non-Newtonian fluid flow models. The method that we have adopted is also prosperous for tackling wide range of nonlinear problems in non-Newtonian fluid mechanics.



---

*Source: 573170-2013-05-09.xml*