# An Exact Solution of MHD Boundary Layer Flow of Dusty Fluid over a Stretching Surface

Authors: Mudassar Jalil; Saleem Asghar; Shagufta Yasmeen
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

This paper deals with the boundary layer flow of electrically conducting dusty fluid over a stretching surface in the presence of applied magnetic field. The governing partial differential equations of the problem are transformed to nonlinear nondimensional coupled ordinary differential equations using suitable similarity transformations. The problem is now fully specified in terms of characterizing parameters known as fluid particle interaction parameter, magnetic field parameter, and mass concentration of dust particles. An exact analytical solution of the resulting boundary value problem is presented that works for all values of the characterizing parameters. The effects of these parameters on the velocity field and the skin friction coefficient are presented graphically and in the tabular form, respectively. We emphasize that an approximate numerical solution of this problem was available in the literature but no analytical solution was presented before this study.


---

## Body

## 1. Introduction

During the last few decades several exact analytical solutions are developed for the one-dimensional flow of Newtonian and non-Newtonian fluids; however, only limited success has been achieved for two- and three-dimensional flows. This is because of the nonlinear nature of the Navier Stokes equations for the viscous fluid and more complex constitutive equations for non-Newtonian fluids. Further difficulty arises when the governing equations are coupled nonlinear equations. Although numerous numerical and approximate techniques are available to compute solution of the fluid dynamics problems, the advantage of the exact analytical solution lies in mathematical ingenuity and the strength of the solution to explain physics of fluid flow. The exact analytical solution can also be considered as bench mark for further numerical and approximate investigations. Some of the exact solutions for the laminar flow of Newtonian and non-Newtonian fluids are given in [1–8].Dusty fluid model flows have been a subject of special interest in recent studies due to their two-phase nature. This phenomenon occurs in fluid (liquid or gas) flows containing a distribution of solid particles. For example, motion of the dusty air in fluidization problems and the chemical process in which raindrops are formed by coalescence of small dust particles. Cosmic dust, which is formed due to the mixing of dust particles and gas, is primary precursor for planetary systems. The production of tails of comet 238 is due to emission of ionized gas and the dust particles from the comet body. The application of the dusty fluid can also be visualized in the processes such as nuclear reactor cooling, atmospheric fallout, powder technology, dust collection, acoustics, paint spray, rain erosion, sedimentation, performance of solid fuel rock nozzles, and guided missiles. These facts have expedited the consideration of modeling, solving, and analyzing the flow of dusty fluids.Keeping interest in two-phase flows, many researchers worked the dusty fluid model for various flow configuration and the boundary conditions. However, realizing the difficulty of nonlinear coupled equations, no attempt has been made in working out any analytical solution. Therefore, the solutions given by them are through numerical and approximate schemes. A brief background history of dusty fluid is presented now. The study on laminar flow of dusty fluid is initiated by Saffman [9]. Chakrabarti [10] investigated the flow of dusty gas in the boundary layer region. The flow of dusty fluid over semi-infinite plate was discussed by Datta and Mishra [11]. Vajravelu and Nayfeh [12] analyzed the hydromagnetic flow of dusty fluid over stretching surface in the presence of suction velocity. Gireesha et al. [13] obtained the numerical solution for the two-dimensional boundary layer flow and heat transfer of dusty fluid over stretching surface. Numerical solution of MHD flow and heat transfer of dusty fluid over a linearly stretching sheet was given by Gireesha et al. [14]. Ramesh et al. [15] considered MHD boundary layer flow of dusty fluid over inclined stretching sheet. Later on, Gireesha et al. [16] examined thermal radiation effects of MHD flow of dusty fluid over exponentially stretching sheet. We once again mention that all of these studies involve numerical or approximate solutions, whereas no exact analytical solution is provided so far.The study of MHD flow is significant due to its many engineering applications such as the cooling of reactors, electrostatic precipitation, power generators, MHD pumps, accelerators, petroleum industry, and the design of heat exchangers. Furthermore, MHD flow plays important role in petroleum industries, agriculture, geophysics, astrophysics, solar physics, metrology, and the motion of earth’s core. Keeping in view the importance of MHD flow and the importance of exact solution, we present exact analytical solution for the MHD boundary layer flow of dusty fluid over a stretching surface. Another distinctive aspect of this study is the presentation of mathematical result for the nonlinear coupled partial differential equations which are usually handled numerically (as has been in the case of dusty fluid). The mathematical methodology consists of converting the modeled equations into self-similar form characterized by dusty and magnetic parameters, using the known similarity transformations. The coupled equations are then solved analytically based on the solution given by Crane [1]. The exact solution is calculated for the flow field and the skin friction in terms of the characterizing parameters such as mass concentration of dust particles, the fluid particle interaction parameter, and the magnetic parameter.


## 2. Description of the Problem

Consider a steady of two-dimensional laminar boundary layer flow of an electrically conducting viscous incompressible dusty fluid over a semi-infinite surface. The surface is stretching with the velocityuw(x)=cx, where the positive constant c is the stretching rate. Cartesian coordinate system is located in such a way that x-axis and y-axis are taken along (and normal to) the surface, respectively, while the origin of the system is located at the leading edge. The shape of dust particles is assumed to be spherical with uniform size and constant number density. Considering these physical assumptions, along with the boundary layer approximations, the governing equations for the flow of dusty fluid together with the boundary conditions are given as (see [14])(1)∂u∂x+∂v∂y=0,u∂u∂x+v∂u∂y=ν∂2u∂y2+γτup-u-σB02ρu,∂up∂x+∂vp∂y=0,up∂up∂x+vp∂up∂y=1τu-up,y=0;u=uwx,v=0,y⟶∞;u=0,up=0,vp=v,where u and up (similarly v and vp) are the x and y components of fluid and the dust particles velocity, respectively. Parameter γ=mN/ρ denotes the mass concentration of dust particles; τ=m/K is the relaxation time of particle phase and ρ, ν, σ, and B0 represent the density, kinematic viscosity, electrical conductivity of the fluid, and induced magnetic field, respectively. The terms appearing in γ and τ are the Stokes resistance (K), the number density (N), and the mass of the dust particle (m).We define the following similarity transformations [14]:(2)η=cνy,u=cxf′η,v=-cνfη,up=cxg′η,vp=-cνgη,where prime (′) represents differentiation with respect to η. Introducing transformations (2) in (1), we get the following self-similar boundary value problem:(3)f′′′+ff′′-f′2+βγg′-f′-Mf′=0,gg′′-g′2+βf′-g′=0,f0=0,f′0=1,f′∞=0,g′η=0,gη=fηasη⟶∞,where β and M are the fluid particle interaction parameter and magnetic parameter, respectively, and are defined by(4)β=1cτ=Kcm,M=σB02cρ.The skin friction coefficient Cf and the shear stress τw are defined as(5)Cf=-τwρuw2,τw=μ∂u∂yy=0.Using transformation (2), the dimensionless skin friction coefficient is defined by Rex-1/2Cf=-f′′(0), where Rex=cx2/υ is the local Reynolds number.In the next section, we will present the exact analytical solution of these equations.


## 3. Exact Analytical Solution

Based upon the solution given by Crane [1], we propose that the nonlinear coupled equations (3) admit exponentially decaying solution of the form(6)f=α1+α2e-α3η,g=α4+α5e-α6η,where αi,s (for i=1,…,6) are arbitrary constants. The values of these constants can be determined by putting the solution (6) in (3). This yields(7)α1=-α2=α4=1ζ,α3=α6=-ζ,α5=-βζ1+β,where the term ζ is given as(8)ζ=1+M+βγ1+β.Substituting the values of αi,s (for i=1,…,6) from (7) into (6), we get closed form exact solution:(9)f=1ζ1-e-ζη,g=1ζ1-β1+βe-ζη,It is worth noting that the solution given in (9) works for all values of β, γ, and M. The physical quantities like the velocities of fluid and the dust particles and the skin friction coefficient are now expressed as(10)f′η=e-ζη,g′η=β1+βe-ζη,-f′′0=ζ=1+M+βγ1+β.The numerical solution of (3) exists in the literature, but the exact analytical closed from solution is presented here for the first time. Now using these results one can easily interpret the effects of dusty fluid parameters on the physical quantities without using any numerical scheme. It will not be out of place to make a comparison with the numerical results presented by Gireesha et al. [17] using Runge Kutta Fehlberg fourth-fifth-order method (RKF45 Method). An excellent match is found between the two results as shown in Table 1.Table 1
Comparison of skin friction coefficient for various values ofβ and M with Gireesha et al. [17].
M
Gireesha et al. [17]β = 0
Exact solution
Gireesha et al. [17]β = 0.5
Exact solution
0.2
1.000
1.000000
1.034
1.033505
0.2
1.095
1.095445
1.126
1.126114
0.5
1.224
1.224745
1.252
1.252251
1.0
1.414
1.414214
1.438
1.438101
1.2
1.483
1.483240
1.506
1.506032
1.5
1.581
1.581139
1.602
1.602540
2.0
1.732
1.732051
1.751
1.751609We also present the graphs for velocity profiles (Figures1(a) and 1(b)) for different values of parameters M and β and compare these results with the numerical results given in [14]. This comparison will further validate the authenticity of the two solutions. From Figure 1, we observe that both the fluid velocity and the velocity of dust particles decrease with increasing M. This is because increasing magnetic field the opposing Lorentz force increases resulting in the decrease of the fluid velocity. Effects of fluid particle interaction β, mass concentration of dust particles γ, and the magnetic parameter M on the skin friction coefficient are given in Table 2. It is observed that skin friction increases with increasing values of β,γ, and M. This happens because increase in these parameters produces resistance to the flow and the skin friction increases.Table 2
Skin friction coefficient for various values ofβ, γ, and M.
γ
β
M
-
f
′
′
(
0
)
0.2
1.437591
0.5
0.5
1.0
1.471960
1.0
1.527525
0.2
1.420094
0.1
0.5
1.0
1.425950
1.0
1.431782
2.0
1.741647
0.1
0.5
5.0
2.456284
10.0
3.321646Figure 1
Effects ofM (a) and β (b) on the velocity profiles.
(a)
(b)


## 4. Conclusion

An exact analytical solution for the MHD boundary layer flow of dusty fluid over a linearly stretching surface is presented. This result is uniformly valid for all values of physical characterizing parameters arising in the flow of dusty fluid. Expressions for the velocity and the skin friction are presented for the flow configuration. The results are compared with the existing study giving good agreement. This work provides a base for further research in undertaking analytical solutions of dusty fluid flow over stretching surface. Exact solution for MHD flow for viscoelastic non-Newtonian fluid can be obtained using the analysis of this paper that will be reported in a subsequent study.



---

*Source: 2307469-2017-03-06.xml*