# Shock Wave Solutions for Some Nonlinear Flow Models Arising in the Study of a Non-Newtonian Third Grade Fluid

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

---

## Abstract

This study is based upon constructing a new class of
closed-form shock wave solutions for some nonlinear problems arising in the
study of a third grade fluid model. The Lie symmetry reduction technique has
been employed to reduce the governing nonlinear partial differential
equations into nonlinear ordinary differential equations. The reduced
equations are then solved analytically, and the shock wave solutions are
constructed. The conditions on the physical parameters of the flow problems
also fall out naturally in the process of the derivation of the solutions.


---

## Body

## 1. Introduction

A shock wave is a disturbance that propagates through a media. Shock is yielded when a disturbance is made to move through a fluid faster than the speed of sound (the celerity) of the medium. This can occur when a solid object is forced through a fluid. It represents a sharp discontinuity of the parameters that delineate the media. Unlike solutions where 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 this is an extremely scarce phenomena. Some interesting communications dealing with the shock wave solutions are found in [1–3].During the past several decades, the study of the nonlinear problems dealing with the flow of non-Newtonian fluids has attracted considerable attention. This interest is due to several important applications in engineering and industry such as reactive polymer flows in heterogeneous porous media, electrochemical generation of elemental bromine in porous electrode systems, manufacture of intumescent paints for fire safety applications, extraction of crude oil from the petroleum products, synthetic fibers, and paper production [4]. Due to the diverse physical nature of non-Newtonian fluids, there is no single constitutive expression which describe the physical behavior of all non-Newtonian fluid models. Because of this important issue, several models of non-Newtonian fluids have been proposed in the literature. Together with this factor, the mathematical modelling of non-Newtonian incompressible fluid flows gives rise to nonlinear and complicated differential equations. As a consequence of this nonlinearity factor, the exact (closed-form) solutions of these sort of problems are scarce in general. Several techniques and methods have been developed in the recent few years to construct the solutions of these non-Newtonian fluid flow problems. Some of these useful methods are variational iteration method, Adomian decomposition method, homotopy perturbation method, homotopy analysis method, and semi-inverse variational method. Literature survey witnesses that, despite all these methods, the exact (closed-form) solutions of the non-Newtonian fluid flow problems are still rare in the literature, and it is not surprising that new exact (closed-form) solutions are most welcome, provided they correspond to physically realistic situations. Some interesting and useful communications in this area are made in the studies [5–14].One of the widely accepted models amongst non-Newtonian fluids is the class of Rivlin-Ericksen fluids of differential type [15]. Rivlin-Ericksen fluids of differential type have acquired the special status in order to describe the several nonstandard features such as normal stress effects, rod climbing, shear thinning, and shear thickening. In the literature much attention has been given to the flow of a second grade fluid [16–19]. A second grade fluid model is the simplest subclass of non-Newtonian fluids for which one can reasonably hope to obtain an analytical solution. In most of the flow situations, the governing equations for a second grade fluid are linear. Although a second grade fluid model is able to predict the normal stress differences, it does not take into account the shear thinning, and shear thickening phenomena that many fluids reveal. Therefore some experiments may be well described by third grade fluid [20–24]. The mathematical model of a third grade fluid represents a more realistic description of the behavior of non-Newtonian fluids. A third grade fluid model represents a further attempt towards the study of the flow properties of non-Newtonian fluids. Therefore, a third grade fluid model has been considered in this study. This model is known to capture the non-Newtonian effects such as shear thinning or shear thickening as well as normal stress.In this particular study, we have constructed the shock wave solutions of some nonlinear PDEs arising in the study of third grade fluid flow in porous medium. We know that the flow of non-Newtonian fluids in a porous medium has wide range of engineering applications. These include enhanced oil recovery, ceramic processing, and geothermal engineering. Motivated by these facts, we have developed some shock wave solutions of three nonlinear problems dealing with the unsteady flow of third grade fluid in a porous half space.


## 2. Geometry of the Models

Consider a Cartesian coordinate frame of referenceOXYZ with x-axis along the direction of the flow and y-axis pointing in the vertically upward direction. The third grade fluid occupies the porous space y>0 and is in contact with an infinite plate at y=0. Since the plate is infinite in the XZ-plane, therefore all the physical quantities except the pressure depend on y only. We have taken three different problems on the same flat plate geometry. The geometry of the problems is shown in Figure 1.Figure 1
Geometry of the physical models and coordinate system.


## 3. Problems to Be Investigated

### 3.1. Unsteady Flow of a Third Grade Fluid over a Flat Rigid Plate with Porous Medium

Following the methodology of [22, 23], the unsteady incompressible flow of a third grade fluid over the rigid plate with porous medium is governed by
(1)(ρ+α1φκ)∂u∂t=μ∂2u∂y2+α1∂3u∂y2∂t+6β3(∂u∂y)2∂2u∂y2-2β3φκ(∂u∂y)2u-φκμu,
where u is the velocity component in x-direction, t is time, y is spatial variable, ρ is the density, μ is the coefficient of viscosity, α1 and β3 are the material constants (for details on these material constants and the conditions that are satisfied by these constants, the reader is referred to [24]), φ is the porosity and κ is the permeability of the porous medium.In order to solve (1), the relevant boundary and initial conditions are
(2)u(0,t)=u0g(t),t>0,(3)u(∞,t)=0,t>0,(4)u(y,0)=f(y),y>0,
where u0 is the reference velocity and g(t) and f(y) are yet functions to be determined. The first boundary condition (2) is the no-slip condition, and the second boundary condition (3) says that the main stream velocity is zero. This is not a restrictive assumption since we can always measure velocity relative to the main stream. The initial condition (4) indicates initially that the fluid is moving with some nonuniform velocity f(y).


#### 3.1.1. Reduction of the Governing Equation

We know that from the principal of Lie symmetry methods [25, 26] 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 (1) admits Lie point symmetry generators, ∂/∂t (translation in t) and ∂/∂y (translation in y). Let X1 and X2 be time-translation and space-translation symmetry generators, respectively. Then the solution corresponding to the generator
(5)X=X1+mX2=∂∂t+m∂∂y,(m>0)
would represent travelling wave solution with constant wave speed m. Travelling wave solutions are characterized by the fact that the profiles of these solutions at different time instants are obtained from one another by appropriate shifts (translations) along the y-axis. Consequently, a Cartesian coordinate system moving with the constant speed can be introduced in which the profile of the derived quantity is stationary. m=0 represents the stationary or steady-state solutions.The characteristic system corresponding to (5) is
(6)dym=dt1=du0.
Solving (6), invariant solutions are given by
(7)u(y,t)=F(η)withη=y-mt,
where F(η) is an arbitrary function of the characteristic variable η=y-mt. Making use of (7) into (1) results in a third-order ordinary differential for F(η) as follows:
(8)(ρ+α1φκ)mdFdη+μd2Fdη2-mα1d3Fdη3+6β3(dFdη)2d2Fdη2-2β3φκF(dFdη)2-μφκF=0.
Thus the original third-order nonlinear PDE (1) is reduced to a third-order ODE (5) along certain curves in the y-t plane. These curves are called characteristic curves or just the characteristic.


#### 3.1.2. Shock Wave Solution

In this section, we show that the travelling wave solutions of (1) approach a shock wave solution. Now we construct the shock wave solution of the reduced equation (5). The starting hypothesis for shock wave solution is given by
(9)F(η)=aexp(bη),
where a and b are the free parameters to be determined. Substituting (9) into (8), we obtain
(10)[m(ρ+α1φκ)b+μb2-mα1b3-μφκ]+e2bη[6β3a2b4-2β3φκa2b2]=0.
Separating (10) in the powers of e0 and e2bη, we find
(11)e0:m(ρ+α1φκ)b+μb2-mα1b3-μφκ=0,(12)e2Bη:β3a2[3b4-φκb2]=0,withβ3a2≠0.
From (12), we deduce
(13)b=φ3κ.
Using the value of b in (11), we obtain
(14)m(ρ+α1φκ)φ3κ+μ(φ3κ)-mα1(φ3κ)φ3κ-μφκ=0.
Finally, the solution for F(η) (provided the condition (14) holds) is written as
(15)F(η)=aexp[φ3κη].
So the solution u(y,t) which satisfies the condition (14) is written as
(16)u(y,t)=aexp[φ3κ(y-mt)]withm>0.Remark 1.
Note that the solution (16) is the shock wave solution to the governing PDE (1). The previous solution is valid under the particular condition on the physical parameters of the flow given in (14) (which is some kind of dispersion relation in m). 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
(17)∂u∂y→∞asy→∞.Remark 2.
Note that the solution (16) also satisfies the particular initial and the boundary condition; that is,
(18)u(0,t)=g(t)=exp[-φ3κmt],t>0,u(y,0)=f(y)=exp[φ3κy],y>0,
with
(19)g(0)=f(0)=a=1.
The functions g(t) and f(y) depend on the physical parameters of the flow.Remark 3.
We also observe that the physical significance of the imposing condition (14) is that it gives the speed of travelling shock wave. From (14), we deduce
(20)m=2μ(φ/3κ)1/2(ρ+(2α1φ/3κ))>0.


### 3.2. Unsteady Magnetohydrodynamic (MHD) Flow of Third Grade Fluid in a Porous Medium

By employing the same geometry as we have explained in Section2, in this problem we extend the previous model by considering the fluid to be electrically conducting under the influence of a uniform magnetic field applied transversely to the flow. We provide the closed-form solution of the problem by reducing the governing nonlinear PDE into an ODE with the help of Lie reduction technique.The time-dependent magnetohydrodynamic flow of a third grade fluid in a porous half space in the absence of the modified pressure gradient takes the form(21)(ρ+α1φκ)∂v∂t=μ∂2v∂y2+α1∂3v∂y2∂t+6β3(∂v∂y)2∂2v∂y2-2β3φκ(∂v∂y)2u-φκμv-σB02v,
where v is the velocity component in x-direction, σ is the electrical conductivity, and B0 is the uniform applied magnetic field. In order to solve (21), the relevant time and space dependent velocity boundary conditions are
(22)u(0,t)=u0g(t),t>0,u(∞,t)=0,t>0,u(y,0)=f(y),y>0.
As it can be seen (21) also admits Lie point symmetry generators, ∂/∂t (translation in t) and ∂/∂y (translation in y). Let X1 and X2 be time-translation and space-translation symmetry generators, respectively. The invariant solution corresponding to the generator X=X1+mX2 is given by
(23)v(y,t)=G(ξ)withξ=y-mt.
Using (23) into (21) yields
(24)(ρ+α1φκ)mdGdξ+μd2Gdξ2-mα1d3Gdξ3+6β3(dGdξ)2d2Gdξ2-2β3φκG(dGdξ)2-μφκG-σB02G=0.
Following the same methodology adopted to solve the previous problem, the reduced ODE (24) admits an exact solution of the form
(25)G(ξ)=aexp[φ3κξ],
provided that
(26)m(ρ+α1φκ)φ3κ+μ(φ3κ)-mα1(φ3κ)φ3κ-μφκ-σB02=0.
Thus the solution of the PDE (21) which satisfies the condition (26) is written as
(27)v(y,t)=exp[φ3κ(y-mt)]withm>0.Remark 4.
Note that the previous solution (27) also satisfies the boundary and initial conditions given in (22). The imposing physical condition (26) gives the speed of travelling shock wave. From (26), we find
(28)m=2μ(φ/3κ)+σB02(φ/3κ)1/2[ρ+(2α1φ/3κ)]>0.
If we set B0=0 (no magnetic field), we recover the condition given in (14).


### 3.3. Unsteady Magnetohydrodynamic (MHD) Flow of Third Grade Fluid in a Porous Medium with Plate Suction/Injection

This particular model is an extension of previous two problems with combined effects of plate suction/injection and MHD nature of the fluid. Thus, for flow under consideration, we seek a velocity of the form(29)V=[w(y,t),-W,0],
where w denotes the velocity of the fluid in x-direction and W>0 indicates suction velocity and W<0 blowing or injection velocity.The unsteady MHD flow of a third grade fluid in a porous half space with plate suction/injection is governed by(30)(ρ+α1φκ)∂w∂t=(ρ+α1φκ)W∂w∂y+μ∂2w∂y2+α1∂3w∂y2∂t+6β3(∂w∂y)2∂2w∂y2-α1W∂3w∂y3-2β3φκ(∂w∂y)2u-φκμw-σB02w.The PDE (30) is solved subject to the same boundary and initial conditions specified for the previous models. The invariant solution of the previous nonlinear PDE under the linear combination of time-translation and space-translation symmetry generators is given by
(31)w(y,t)=H(θ)withθ=y-mt.
Inserting (31) in (30), we get a third-order nonlinear ODE in H(θ), namely,
(32)(ρ+α1φκ)mdHdθ+(ρ+α1φκ)WdHdθ+μd2Hdθ2-mα1d3Hdθ3+6β3(dHdθ)2d2Hdθ2-α1W∂3H∂θ3-2β3φκH(dHdθ)2-(μφκ+σB02)H(θ)=0.
Following the same procedure used to tackle the first problem, (32) admits the exact solution of the form
(33)H(θ)=aexp[φ3κθ],
provided
(34)0=m(ρ+α1φκ)φ3κ+(ρ+α1φκ)Wφ3κ+μ(φ3κ)-(m+W)α1(φ3κ)φ3κ-(μφκ+σB02).
The solution of w(y,t) is written as
(35)w(y,t)=exp[φ3κ(y-mt)]withm>0.Remark 5.
The previous solution is valid only under the particular condition on the physical parameters given in (34). The condition (34) also gives the speed of shock wave. Thus from (34), we obtain
(36)m=2μ(φ/3κ)+σB02-(ρ+(α1φ/3κ))φ/3κWφ/3κ[ρ+(2α1φ/3κ)-(α1φ/3κ)W]>0.
Note that, if we set W=B0=0 (with no porosity and magnetic field), we recover the previous two conditions given in (14) and (26).Remark 6.
We note that the shock wave solutions (16), (27), and (35) are the same, but the imposing conditions on the physical parameters of the flow given in (14), (26), and (34) under which these solutions are valid are different. This means that in each case the speed of the travelling shock wave is different. Therefore, the graphical behavior of these solutions is the same which shows the shock wave behavior of the flow problems. However, the imposing conditions contain the magnetic field, suction/blowing, porosity, and second grade and the third grade parameters. Thus these closed-form shock wave solutions are valid for the particular values of these parameters.


## 3.1. Unsteady Flow of a Third Grade Fluid over a Flat Rigid Plate with Porous Medium

Following the methodology of [22, 23], the unsteady incompressible flow of a third grade fluid over the rigid plate with porous medium is governed by
(1)(ρ+α1φκ)∂u∂t=μ∂2u∂y2+α1∂3u∂y2∂t+6β3(∂u∂y)2∂2u∂y2-2β3φκ(∂u∂y)2u-φκμu,
where u is the velocity component in x-direction, t is time, y is spatial variable, ρ is the density, μ is the coefficient of viscosity, α1 and β3 are the material constants (for details on these material constants and the conditions that are satisfied by these constants, the reader is referred to [24]), φ is the porosity and κ is the permeability of the porous medium.In order to solve (1), the relevant boundary and initial conditions are
(2)u(0,t)=u0g(t),t>0,(3)u(∞,t)=0,t>0,(4)u(y,0)=f(y),y>0,
where u0 is the reference velocity and g(t) and f(y) are yet functions to be determined. The first boundary condition (2) is the no-slip condition, and the second boundary condition (3) says that the main stream velocity is zero. This is not a restrictive assumption since we can always measure velocity relative to the main stream. The initial condition (4) indicates initially that the fluid is moving with some nonuniform velocity f(y).


### 3.1.1. Reduction of the Governing Equation

We know that from the principal of Lie symmetry methods [25, 26] 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 (1) admits Lie point symmetry generators, ∂/∂t (translation in t) and ∂/∂y (translation in y). Let X1 and X2 be time-translation and space-translation symmetry generators, respectively. Then the solution corresponding to the generator
(5)X=X1+mX2=∂∂t+m∂∂y,(m>0)
would represent travelling wave solution with constant wave speed m. Travelling wave solutions are characterized by the fact that the profiles of these solutions at different time instants are obtained from one another by appropriate shifts (translations) along the y-axis. Consequently, a Cartesian coordinate system moving with the constant speed can be introduced in which the profile of the derived quantity is stationary. m=0 represents the stationary or steady-state solutions.The characteristic system corresponding to (5) is
(6)dym=dt1=du0.
Solving (6), invariant solutions are given by
(7)u(y,t)=F(η)withη=y-mt,
where F(η) is an arbitrary function of the characteristic variable η=y-mt. Making use of (7) into (1) results in a third-order ordinary differential for F(η) as follows:
(8)(ρ+α1φκ)mdFdη+μd2Fdη2-mα1d3Fdη3+6β3(dFdη)2d2Fdη2-2β3φκF(dFdη)2-μφκF=0.
Thus the original third-order nonlinear PDE (1) is reduced to a third-order ODE (5) along certain curves in the y-t plane. These curves are called characteristic curves or just the characteristic.


### 3.1.2. Shock Wave Solution

In this section, we show that the travelling wave solutions of (1) approach a shock wave solution. Now we construct the shock wave solution of the reduced equation (5). The starting hypothesis for shock wave solution is given by
(9)F(η)=aexp(bη),
where a and b are the free parameters to be determined. Substituting (9) into (8), we obtain
(10)[m(ρ+α1φκ)b+μb2-mα1b3-μφκ]+e2bη[6β3a2b4-2β3φκa2b2]=0.
Separating (10) in the powers of e0 and e2bη, we find
(11)e0:m(ρ+α1φκ)b+μb2-mα1b3-μφκ=0,(12)e2Bη:β3a2[3b4-φκb2]=0,withβ3a2≠0.
From (12), we deduce
(13)b=φ3κ.
Using the value of b in (11), we obtain
(14)m(ρ+α1φκ)φ3κ+μ(φ3κ)-mα1(φ3κ)φ3κ-μφκ=0.
Finally, the solution for F(η) (provided the condition (14) holds) is written as
(15)F(η)=aexp[φ3κη].
So the solution u(y,t) which satisfies the condition (14) is written as
(16)u(y,t)=aexp[φ3κ(y-mt)]withm>0.Remark 1.
Note that the solution (16) is the shock wave solution to the governing PDE (1). The previous solution is valid under the particular condition on the physical parameters of the flow given in (14) (which is some kind of dispersion relation in m). 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
(17)∂u∂y→∞asy→∞.Remark 2.
Note that the solution (16) also satisfies the particular initial and the boundary condition; that is,
(18)u(0,t)=g(t)=exp[-φ3κmt],t>0,u(y,0)=f(y)=exp[φ3κy],y>0,
with
(19)g(0)=f(0)=a=1.
The functions g(t) and f(y) depend on the physical parameters of the flow.Remark 3.
We also observe that the physical significance of the imposing condition (14) is that it gives the speed of travelling shock wave. From (14), we deduce
(20)m=2μ(φ/3κ)1/2(ρ+(2α1φ/3κ))>0.


## 3.1.1. Reduction of the Governing Equation

We know that from the principal of Lie symmetry methods [25, 26] 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 (1) admits Lie point symmetry generators, ∂/∂t (translation in t) and ∂/∂y (translation in y). Let X1 and X2 be time-translation and space-translation symmetry generators, respectively. Then the solution corresponding to the generator
(5)X=X1+mX2=∂∂t+m∂∂y,(m>0)
would represent travelling wave solution with constant wave speed m. Travelling wave solutions are characterized by the fact that the profiles of these solutions at different time instants are obtained from one another by appropriate shifts (translations) along the y-axis. Consequently, a Cartesian coordinate system moving with the constant speed can be introduced in which the profile of the derived quantity is stationary. m=0 represents the stationary or steady-state solutions.The characteristic system corresponding to (5) is
(6)dym=dt1=du0.
Solving (6), invariant solutions are given by
(7)u(y,t)=F(η)withη=y-mt,
where F(η) is an arbitrary function of the characteristic variable η=y-mt. Making use of (7) into (1) results in a third-order ordinary differential for F(η) as follows:
(8)(ρ+α1φκ)mdFdη+μd2Fdη2-mα1d3Fdη3+6β3(dFdη)2d2Fdη2-2β3φκF(dFdη)2-μφκF=0.
Thus the original third-order nonlinear PDE (1) is reduced to a third-order ODE (5) along certain curves in the y-t plane. These curves are called characteristic curves or just the characteristic.


## 3.1.2. Shock Wave Solution

In this section, we show that the travelling wave solutions of (1) approach a shock wave solution. Now we construct the shock wave solution of the reduced equation (5). The starting hypothesis for shock wave solution is given by
(9)F(η)=aexp(bη),
where a and b are the free parameters to be determined. Substituting (9) into (8), we obtain
(10)[m(ρ+α1φκ)b+μb2-mα1b3-μφκ]+e2bη[6β3a2b4-2β3φκa2b2]=0.
Separating (10) in the powers of e0 and e2bη, we find
(11)e0:m(ρ+α1φκ)b+μb2-mα1b3-μφκ=0,(12)e2Bη:β3a2[3b4-φκb2]=0,withβ3a2≠0.
From (12), we deduce
(13)b=φ3κ.
Using the value of b in (11), we obtain
(14)m(ρ+α1φκ)φ3κ+μ(φ3κ)-mα1(φ3κ)φ3κ-μφκ=0.
Finally, the solution for F(η) (provided the condition (14) holds) is written as
(15)F(η)=aexp[φ3κη].
So the solution u(y,t) which satisfies the condition (14) is written as
(16)u(y,t)=aexp[φ3κ(y-mt)]withm>0.Remark 1.
Note that the solution (16) is the shock wave solution to the governing PDE (1). The previous solution is valid under the particular condition on the physical parameters of the flow given in (14) (which is some kind of dispersion relation in m). 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
(17)∂u∂y→∞asy→∞.Remark 2.
Note that the solution (16) also satisfies the particular initial and the boundary condition; that is,
(18)u(0,t)=g(t)=exp[-φ3κmt],t>0,u(y,0)=f(y)=exp[φ3κy],y>0,
with
(19)g(0)=f(0)=a=1.
The functions g(t) and f(y) depend on the physical parameters of the flow.Remark 3.
We also observe that the physical significance of the imposing condition (14) is that it gives the speed of travelling shock wave. From (14), we deduce
(20)m=2μ(φ/3κ)1/2(ρ+(2α1φ/3κ))>0.


## 3.2. Unsteady Magnetohydrodynamic (MHD) Flow of Third Grade Fluid in a Porous Medium

By employing the same geometry as we have explained in Section2, in this problem we extend the previous model by considering the fluid to be electrically conducting under the influence of a uniform magnetic field applied transversely to the flow. We provide the closed-form solution of the problem by reducing the governing nonlinear PDE into an ODE with the help of Lie reduction technique.The time-dependent magnetohydrodynamic flow of a third grade fluid in a porous half space in the absence of the modified pressure gradient takes the form(21)(ρ+α1φκ)∂v∂t=μ∂2v∂y2+α1∂3v∂y2∂t+6β3(∂v∂y)2∂2v∂y2-2β3φκ(∂v∂y)2u-φκμv-σB02v,
where v is the velocity component in x-direction, σ is the electrical conductivity, and B0 is the uniform applied magnetic field. In order to solve (21), the relevant time and space dependent velocity boundary conditions are
(22)u(0,t)=u0g(t),t>0,u(∞,t)=0,t>0,u(y,0)=f(y),y>0.
As it can be seen (21) also admits Lie point symmetry generators, ∂/∂t (translation in t) and ∂/∂y (translation in y). Let X1 and X2 be time-translation and space-translation symmetry generators, respectively. The invariant solution corresponding to the generator X=X1+mX2 is given by
(23)v(y,t)=G(ξ)withξ=y-mt.
Using (23) into (21) yields
(24)(ρ+α1φκ)mdGdξ+μd2Gdξ2-mα1d3Gdξ3+6β3(dGdξ)2d2Gdξ2-2β3φκG(dGdξ)2-μφκG-σB02G=0.
Following the same methodology adopted to solve the previous problem, the reduced ODE (24) admits an exact solution of the form
(25)G(ξ)=aexp[φ3κξ],
provided that
(26)m(ρ+α1φκ)φ3κ+μ(φ3κ)-mα1(φ3κ)φ3κ-μφκ-σB02=0.
Thus the solution of the PDE (21) which satisfies the condition (26) is written as
(27)v(y,t)=exp[φ3κ(y-mt)]withm>0.Remark 4.
Note that the previous solution (27) also satisfies the boundary and initial conditions given in (22). The imposing physical condition (26) gives the speed of travelling shock wave. From (26), we find
(28)m=2μ(φ/3κ)+σB02(φ/3κ)1/2[ρ+(2α1φ/3κ)]>0.
If we set B0=0 (no magnetic field), we recover the condition given in (14).


## 3.3. Unsteady Magnetohydrodynamic (MHD) Flow of Third Grade Fluid in a Porous Medium with Plate Suction/Injection

This particular model is an extension of previous two problems with combined effects of plate suction/injection and MHD nature of the fluid. Thus, for flow under consideration, we seek a velocity of the form(29)V=[w(y,t),-W,0],
where w denotes the velocity of the fluid in x-direction and W>0 indicates suction velocity and W<0 blowing or injection velocity.The unsteady MHD flow of a third grade fluid in a porous half space with plate suction/injection is governed by(30)(ρ+α1φκ)∂w∂t=(ρ+α1φκ)W∂w∂y+μ∂2w∂y2+α1∂3w∂y2∂t+6β3(∂w∂y)2∂2w∂y2-α1W∂3w∂y3-2β3φκ(∂w∂y)2u-φκμw-σB02w.The PDE (30) is solved subject to the same boundary and initial conditions specified for the previous models. The invariant solution of the previous nonlinear PDE under the linear combination of time-translation and space-translation symmetry generators is given by
(31)w(y,t)=H(θ)withθ=y-mt.
Inserting (31) in (30), we get a third-order nonlinear ODE in H(θ), namely,
(32)(ρ+α1φκ)mdHdθ+(ρ+α1φκ)WdHdθ+μd2Hdθ2-mα1d3Hdθ3+6β3(dHdθ)2d2Hdθ2-α1W∂3H∂θ3-2β3φκH(dHdθ)2-(μφκ+σB02)H(θ)=0.
Following the same procedure used to tackle the first problem, (32) admits the exact solution of the form
(33)H(θ)=aexp[φ3κθ],
provided
(34)0=m(ρ+α1φκ)φ3κ+(ρ+α1φκ)Wφ3κ+μ(φ3κ)-(m+W)α1(φ3κ)φ3κ-(μφκ+σB02).
The solution of w(y,t) is written as
(35)w(y,t)=exp[φ3κ(y-mt)]withm>0.Remark 5.
The previous solution is valid only under the particular condition on the physical parameters given in (34). The condition (34) also gives the speed of shock wave. Thus from (34), we obtain
(36)m=2μ(φ/3κ)+σB02-(ρ+(α1φ/3κ))φ/3κWφ/3κ[ρ+(2α1φ/3κ)-(α1φ/3κ)W]>0.
Note that, if we set W=B0=0 (with no porosity and magnetic field), we recover the previous two conditions given in (14) and (26).Remark 6.
We note that the shock wave solutions (16), (27), and (35) are the same, but the imposing conditions on the physical parameters of the flow given in (14), (26), and (34) under which these solutions are valid are different. This means that in each case the speed of the travelling shock wave is different. Therefore, the graphical behavior of these solutions is the same which shows the shock wave behavior of the flow problems. However, the imposing conditions contain the magnetic field, suction/blowing, porosity, and second grade and the third grade parameters. Thus these closed-form shock wave solutions are valid for the particular values of these parameters.


## 4. Shock Wave Behavior of the Solutions

Figures2 and 3 show the shock wave behavior of the solutions (16), (27), and (35) in 2D and 3D, respectively. From the graphs it is quite clear that the slope of the velocity profile is approaching to infinity representing the sudden sharp discontinuity in the velocity field.Figure 2
Profile of the shock wave solutions (16), (27), and (35) with φ=4, κ=0.2, m=1, and t=π/2.Figure 3
3D profile of the shock wave solutions (16), (27), and (35) with φ=4, κ=0.2, m=1, and t=π/2.


## 5. Final Comments

In this paper, we have presented closed-form shock wave solutions for some nonlinear problems which describe the phenomena of third grade fluids. In each case the governing nonlinear PDEs are reduced to nonlinear ODEs by using the Lie point symmetry (which is translation) in thet and y directions. The reduced ODEs are then solved analytically. We observe that the shock wave solutions (16), (27), and (35) are the same, but the imposing conditions on the physical parameters of the flow models given in (14), (26), and (34) under which these solutions are valid are different. These solutions do not directly contain the parameter which is responsible for showing the behavior of third grade fluid parameter on the flow. However, the imposing conditions under which these solutions are valid do contain the third grade parameter. To emphasize, we can say that these solutions are valid for the particular values of third grade fluid parameter. The results obtained describe the mathematical structure of the shock wave behavior of the flow problems. The models considered in this study are prototype, but the obtained solutions are going to be very helpful in carrying out further analysis of the shock wave characteristic 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: 602902-2013-06-23.xml*