Next Article in Journal
The Influence of Environmental Change (Crops and Water) on Population Redistribution in Mexico and Ethiopia
Previous Article in Journal
A Four-Probe Method Using Different Probe Spacings for Measurement and Exact Reconstruction of Parallel Profiles
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Finite Element Analysis of Thermo-Diffusion and Multi-Slip Effects on MHD Unsteady Flow of Casson Nano-Fluid over a Shrinking/Stretching Sheet with Radiation and Heat Source

1
School of Energy and Power, Xi’an Jiaotong University, No.28West Xianning Road, Xi’an 710049, China
2
Department of Applied Mathematics, Northwestern Polytechnical University, Dongxiang Road, Beilin District, Xi’an 710129, China
3
School of Management, Xi’an Jiaotong University, No.28West Xianning Road, Xi’an 710049, China
4
School of Mathematics, Northwest University, No.229 North Taibai Avenue, Xi’an 710069, China
*
Author to whom correspondence should be addressed.
Submission received: 22 September 2019 / Revised: 24 November 2019 / Accepted: 26 November 2019 / Published: 30 November 2019
(This article belongs to the Section Nanotechnology and Applied Nanosciences)

Abstract

:
In this article, we probe the multiple-slip effects on magnetohydrodynamic unsteady Casson nano-fluid flow over a penetrable stretching sheet, sheet entrenched in a porous medium with thermo-diffusion effect, and injection/suction in the presence of heat source. The flow is engendered due to the unsteady time-dependent stretching sheet retained inside the porous medium. The leading non-linear partial differential equations are transmuted in the system of coupled nonlinear ordinary differential equations by using appropriate transformations, then the transformed equations are solved by using the variational finite element method numerically. The velocity, temperature, solutal concentration, and nano-particles concentration, as well as the rate of heat transfer, the skin friction coefficient, and Sherwood number for solutal concentration, are presented for several physical parameters. Next, the effects of these various physical parameters are conferred with graphs and tables. The exact values of flow velocity, skin friction, and Nusselt number are compared with a numerical solution acquired with the finite element method (FEM), and also with numerical results accessible in literature. In the end, we rationalize the convergence of the finite element numerical solution, and the calculations are carried out by reducing the mesh size.

1. Introduction

The basic idea of the non-slip condition is concerned with the Navier Stokes theory. To study the temperature and velocity, numerous authors have found both the analytical and numerical solution by the implementation of non-slip boundary conditions. The utmost significance of slip-conditions in nano-channels has stimulated, to a great extent, interest in the study of vibrating values [1]. Now it is not unknown that if fluid contains concentrated suspensions, in that case the slip could be stirring. Soltani and Yilmazer [2] have performed by parallel disk rheometer with prominence on the wall slip phenomena on the rheological characterization of extremely filled suspensions consisting of Newtonian matrix and diverse with two different power sizes of aluminum and two different sizes of bed glass. In the articulations of fluid-like as suspension, emulsion foams, and polymer solution, partial velocity slip can occur on the stretching boundary. In an assortment of industry procedures, the influences of slip can crop-up at the boundary of pipes, surfaces, and walls. The Navier Stoke velocity slip conditions are customary to approach in the study of slip phenomenon.
There are hundreds of applications in industries and engineering, and magneto-hydrodynamic fluid flows through stretching sheet have achieved much importance in recent days, according to Mabbod et al. [3]. Applications as well as the liquid coating on the photographic films, boundary layers throughout the liquid film in the concentration procedure, and the aerodynamic excrescence of plastic sheets exist. With this, the extensive range of applications of magnetohydrodynamic flow that can be found in copious fields like in electronic cooling process, in boilers, heat lagging, and metal extrusion, geothermal system, nuclear process, micro-magnetohydrodynamic pumps, underground water system, in energy storage units, biological conveyance, and in the thermal energy procedure has played a very important role. Unsteady flow due to the stretching sheet has been scrutinized by Pop and Na [4]. Afterward, Sheridan et al. [5] investigated the significance of variable wall temperature and variable heat flux in the boundary layer flow over the unsteady stretching surface with similarity transformation. An unsteady stagnation point flow of the viscous fluid caused by a stretching sheet under the influence of slip condition has been reported by Bhattacharyya et al. [6]. The impact of the magnetic field on the two-dimensional flow of nano-fluid with and without slip condition was conferred by Khan et al. [7] and Mohyud-Din et al. [8,9], correspondingly. Khan et al. [10] investigated two-dimensional electrically conducting a flow of nano-fluid due to the stretching sheet under the impact of convective boundary condition.
The wide-ranging applications of non-Newtonian liquids in numerous industries such as molten plastics, nuclear waste disposal, artificial fibers, foodstuffs, transpiration cooling, and petroleum reservoirs have gained fascination to analyze their characteristic comportment. Casson fluid is one sort of fluid model for non-Newtonian fluid. Casson fluid has an inestimable viscosity at zero rates of shear and have zero viscosity at an infinite rate of shear and the yield stress under which there is no flow occurs. Consequently, Casson-fluid could be defined such as shear-thinning liquid or pseudo-plastic fluid. The flow characteristic of blood can precisely be defined by Casson fluid at the small shear rates. Sharada and Shankar [11] scrutinized numerically Soret and Dufour effects, thermal radiation, and chemical reaction on fluid flow, and heat and mass transfers of Casson fluid over an exponentially stretching surface. Ullah et al. [12] observed that unsteady MHD diverse convection slip flow of Casson fluid towards non-linearly stretching sheet soaked in a porous medium in existence of slip and convective boundary conditions occurs. Furthermore, Ullah et al. [13] acquired numerical results for hydro-magnetic Falkner Skan flow of the Casson fluid past a moving wedge, along with heat transfer. Pushpalatha et al. [14] deliberated thermal diffusion effects and radiation on Casson-fluid flow with convective boundary conditions. Pushpalatha et al. [15] deliberated the thermos-diffusion consequences on the Casson fluid over an unsteady stretching surface in the existence of thermal radiations and the magnetic field. They found that the Casson parameters and unsteadiness parameter have the propensity to denigrate the velocity dissemination. Uddin et al. [16] studied the free convective boundary layer flow of nano-fluid over a penetrable upward horizontal plate in a porous medium with thermal convective boundary conditions.
Recently, Rahman et al. [17] explored the hydro-magnetic slip flow of water-based nano-fluid over a wedge with a convective surface in the existence of heat generation or absorption. Ibrahim and Shankar [18] analyzed magneto-hydrodynamic boundary layer flow and heat transfer of nano-fluid past over penetrable stretching sheet using velocity, thermal, and solute slip boundary conditions. Das [19] presented an imperative numerical inquiry on the convective heat transfer performance of nanofluids over a penetrable stretching surface in the presence of partial slip, thermal buoyancy, and internal heat generation/absorption. The transfer of heat, investigated through the thermal radiation, has enormous applications in different technical processes, satellites, nuclear power-plants, astronomical vehicles, gas turbines, comprising missiles, and in several propulsion devices for jets. Linear radiations are not systematically in effect for the enormous temperature alteration. Abbas et al. [20] deliberated the consequences of radiation in the existence of a uniform magnetic field of nano-fluid on the curved stretching surface by integrating the slip effect. In recent times, Makinde et al. [21] deliberated a numerical study of the effects of radiation on chemically reacting MHD nano-fluid prejudiced by the heat source/sink, and the collective heat and mass transfer analysis aimed at mixed convection flow over the vertical surface, with radiation and chemical reaction explained by Ibrahim et al. [22]. Prasannakumara et al. [23] deliberated the velocity slip effects, temperature jump, solutal slip, and thermal radiation, on steady flow, and the transfer of heat and mass of incompressible Jeffrey nano-fluid over the horizontal stretching surface. Imtiaz et al. [24] scrutinized the unsteady MHD flow of the curved stretching surface.
The literature remains silent whether we have to investigate the impact of multiple slips on MHD unsteady Casson nano-fluid flow, heat, and mass transfer in the presence of heat source with thermo-diffusion effect over a stretching/shrinking sheet. The objective of this article was to prolong the recent work of Fazle Mabood and Standford Shateyi [25]. Appropriate similarities have been used for transformation, and the governing non-linear partial differential equations are rendered in the non-dimensional and non-linear system of ODEs. The resulting system of non-linear ODEs has been solved numerically with an efficient and validated variational finite element method (FEM) with the boundary conditions. We also use a special case for the existing model to compare our outcomes along with previous studies. A parametric study has been performed to inquire into the mass and heat transfer characteristics and also the impact of different parameters of the flow. After that, we have provided a numerical comparison of our results and discussed them with graphs. In future studies, transient flow with slip effects in the presence of mixed convection and chemical reaction at the sheet can be examined.

2. Mathematical Modeling

The unsteady two-dimensional MHD flow of an incompressible fluid and Casson nano-fluid viscous flow over an electrically conducted stretching/shrinking sheet in the presence of heat source was considered. We choose an xy-ordinate system as the measurement of the sheet will be taken with x-, y-axis, and the y-axis along the vertical direction of the sheet, as given in Figure 1. The sheet is moving with non-uniform velocity U ( x , t ) = a x / ( 1 λ t ) since a is the stretching rate along the x-axis but λ t is positively constant, as by property λ t < 1 . The transverse magnetic field is supposed be the function of distance from the origin and is defined as B ( x ) = B 0 x with B 0 0 , where B 0 is the strength of the magnetic field. The induced magnetic field is negligible as compared to the applied magnetic field. The free stream temperature is supposed to be T , the free mass concentration is C , and the nano-particle concentration is E . The governing equations for flow can be put into the form as in ([25,26])
u x + v y = 0 ,
u t + u u x + v u y = ν ( 1 + 1 β ) 2 u y 2 σ B 2 ( x ) u ρ + g β T ( T T ) + g β C ( C C ) + g β E ( E E ) ν κ u
T t + u T x + v T y = α ( 1 + 16 T 3 σ * 3 k * κ ) 2 T y 2 + τ [ D B E y T y + D T T ( T y ) 2 ] + Q ( T T ) ρ C p + D T C 2 C y 2
C t + u C x + v C y = D m 2 C y 2 + D T 2 T y 2
E t + u E x + v E y = D B 2 E y 2 + D T T 2 T y 2
where x and y are the co-ordinates in the x-axis and y-axis; u and v are the velocity components along the x-axis and y-axis respectively; α , ν , σ , ρ are the thermal diffusivity, kinetic viscosity, electrical conductivity, and density of fluid, respectively; g is the acceleration due to gravity; β T is the thermal expansion coefficient; β C is the solutal concentration expansion coefficient; β E is the nano-particle concentration expansion coefficient; T is the temperature; C is the solutal concentration; E is the nano-particle concentration; D M is the molecular diffusivity; D T is the thermal diffusivity; D B is the Brownian diffusivity; σ * is the Stefan-Boltzmann constant; k * is the mean absorption coefficient; T is the stream temperature; Q is the uniform volumetric heat generation/absorption; C p is the specific heat of fluid; C s is the concentration of susceptibility; and K T is the thermal diffusion ratio. The boundary conditions for the above mathematical model are (see [25]):
u = U ( x , t ) + U s l i p , v = v w , T = T w ( x , t ) + T s l i p , C = C w ( x , t ) + C s l i p , E = E w ( x , t ) + E s l i p , a s y = 0
u 0 , T T , C C , E E , a s y ,
where x and y are the co-ordinates; u and v are the velocity components along the x-axis and y-axis, respectively; α , ν , σ , ρ are the thermal diffusivity, kinetic viscosity, electrical conductivity and density of fluid, respectively; g is the acceleration due to gravity; β T is thermal expansion coefficient; β C is the solutal concentration expansion coefficient; β E is the nano-particle concentration expansion coefficient; T is the temperature; C is the solutal concentration; E is the nano-particle concentration; D M is the molecular diffusivity; D T is the thermal diffusivity; D B is the Brownian diffusivity; σ * is the Stefan-Boltzmann constant; k * is the mean absorption coefficient; T is the stream temperature; Q is the uniform volumetric heat generation/absorption; C p is the specific heat of fluid; C s is the concentration of susceptibility; and K T is the thermal diffusion ratio.
Where the injection/suction velocity as v w = v 0 / x , T w ( x , t ) is the temperature of sheet and C w ( x , t ) , E w ( x , t ) are concentrations at surface of the below form (see [25]):
  • T w ( x , t ) = T + T 0 a x 2 ν ( 1 λ t ) 2
  • C w ( x , t ) = C + C 0 a x 2 ν ( 1 λ t ) 2
  • E w ( x , t ) = E + E 0 a x 2 ν ( 1 λ t ) 2
where T 0 , C 0 , and E 0 are the reference temperature, reference solutal concentration, and reference nano-particle concentration respectively, such that 0 T 0 T w , 0 C 0 C w and 0 E 0 E w , with these above expressions being valid if ( 1 λ t ) > 0 .
Usually, the stream function Ψ is defined as u = Ψ y and v = Ψ x that the Equation (1) is satisfied. We introduce similarity transformations to solve the above equations (see [25,26]):
η = a ν ( 1 λ t ) y , ψ = a ν ( 1 λ t ) x f ( η ) , θ ( η ) = T T T w T , ϕ ( η ) = C C C w C , ξ ( η ) = E T E w E
In view of the above similarity transformations of Equation (7), PDEs from (1)–(5) transform into the following system of nonlinear ODEs:
( 1 + 1 β ) f + f f f 2 σ ( η 2 f + f ) M f + λ 1 θ + λ 2 ϕ + λ 3 ξ k p f = 0 ,
( 1 + R ) 1 P r θ f θ + f θ σ η 2 θ + 2 θ + N b ϕ θ + N t θ 2 + Q θ + N d ϕ = 0
1 S c ϕ f ϕ + f ϕ σ η 2 ϕ + 2 ϕ + S r θ = 0
ξ L e f ξ f ξ + σ η 2 ξ + 2 ξ + N t N b θ = 0
for the above problem the transformed boundary conditions are:
f ( 0 ) = f w ,
f ( 0 ) = 1 + S f f ( 0 ) ,
θ ( 0 ) = 1 + S θ θ ( 0 ) ,
ϕ ( 0 ) = 1 + S ϕ ϕ ( 0 ) ,
ξ ( 0 ) = 1 + S ξ ξ ( 0 ) ,
f ( ) 0 , θ ( ) 0 , ϕ ( ) 0 , ξ ( ) 0 ,
where unsteadiness parameter σ = λ / a , M is the magnetic field parameter, P r is the Prandtl number, N b is the Brownian motion parameter, N t is the thermo-phoresis parameter, β is the Deborah number, S c is the Schmidt number, R is the thermal radiation parameter, λ 1 , λ 2 , and λ 3 are the buoyancy parameters, k p is the permeability parameter, N d is the Dufour parameter, S r is the Soret parameter, and f w is the Suction/injection parameter.
The primes show the differentiation with respect to η . The parameters used in Equations (8)–(11) are explained as:
M = [ σ ( 1 λ t ) ρ a B 0 ] 1 / 2 , P r = ν / α , N b = τ D B ( C w C ) / ν , N t = ( τ D T ( T w T ) / ν T , β = k 0 a / ( 1 λ t ) , S c = ν D m , L e = σ D B , R = ( 16 σ * T 3 ) / ( 3 k * K ) , σ = λ / a , λ 1 = ( g β T T 0 ) / ( 2 a ν ) , λ 2 = ( g β C C 0 ) / ( 2 a ν ) , λ 3 = ( g β E E 0 ) ( 2 a ν ) , S r = ( D T T 0 ) / ν C 0 , N d = ( D m K T C 0 ) / ( C s C p T 0 ) , k p = ν ( 1 λ t ) / a , f w = v w ( 1 γ t ) 1 / 2 / ( ν a ) 1 / 2 ,
Also, the local skin friction co-efficient, local Nusselt number, and local Sherwood number are explained as below
C f = μ ρ ( U w ) 2 ( u y ) y = 0
N u = x k ( T T ) [ k ( T y ) y = 0 4 σ 3 k ( T 4 y ) y = 0 ]
S h = x ( C C ) ( C y ) y = 0
when Equation (7) is substitited into Equations (17)–(19), the final dimensionless form is obtained; C f r = ( R e x ) 1 / 2 = f ( 0 ) , N u r = N u / ( R e x ) 1 / 2 = ( 1 + R ) θ ( 0 ) , S h r = S h / ( R e x ) 1 / 2 = ϕ ( 0 ) where the local Reynolds number is R e x , reduced skin friction is C f r , reduced Nusselt number is N u r , and the reduced Sherwood number is S h r .

3. Finite Element Method Solutions

We solve the system of non-linear boundary value problem that is given in Equations (8)–(11) numerically by applying the finite element method (FEM) subjected with the boundary conditions (12)–(16). The FEM is extremely effectual and has been applied to study different problems in fluid mechanics, CFD, mass transfer, heat transfer, solid mechanics, and also in many other fields. The general detail of the finite element method (FEM) can be found in Raddy [26,27,28]. Reddy [29] gives a general detail of the variational finite element method, which also found that the finite element method (FEM) is employed in commercial software like ADINA, ANSYS, MATLAB, and ABAQUS. Swapna et al. and Rana et al. [30,31] explain that the variational finite element method solves the boundary value problem very efficiently and accurately. To solve the non-linear boundary value problem (8)–(11) by using finite element method (FEM), along with boundary conditions (12)–(16), to apply FEM, first we have to consider:
f = g .
The Equations (8)–(11) take the form
( 1 + 1 β ) g + f g g 2 σ ( η 2 g + g ) M g + λ 1 θ + λ 2 ϕ + λ 3 ξ k p g = 0 ,
( 1 + R ) 1 P r θ g θ + f θ σ η 2 θ + 2 θ + N b ϕ θ + N t θ 2 + Q θ + N d ϕ = 0
1 S c ϕ g ϕ + f ϕ σ η 2 ϕ + 2 ϕ + S r θ = 0
ξ L e g ξ f ξ + σ η 2 ξ + 2 ξ + N t N b θ = 0 ,
The corresponding boundary conditions now reduce to the following form:
f ( 0 ) = f w , g ( 0 ) = 1 + S f g ( 0 ) , θ ( 0 ) = 1 + S θ θ ( 0 ) , ϕ ( 0 ) = 1 + S ϕ ϕ ( 0 ) , ξ ( 0 ) = 1 + S ξ ξ ( 0 ) ,
h ( ) 0 , θ ( ) 0 , ϕ ( ) 0 , ξ ( ) 0 ,

3.1. Variational Formulations

The variational formulation is associated with Equations (20)–(24) over a linear element ( η c , η c + 1 ) is given by
η c η c + 1 w 1 { f g } d η = 0 ,
η c η c + 1 w 2 { ( 1 + 1 β ) g + f g g 2 σ ( η 2 g + g ) M g + λ 1 θ + λ 2 ϕ + λ 3 ξ k p g } d η = 0 ,
η c η c + 1 w 3 { ( 1 + R ) 1 P r θ g θ + f θ σ η 2 θ + 2 θ + N b ϕ θ + N t θ 2 + Q θ + N d ϕ } d η = 0 ,
η c η c + 1 w 4 { 1 S c ϕ g ϕ + f ϕ σ η 2 ϕ + 2 ϕ + S r θ } d η = 0 ,
η c η c + 1 w 5 { ξ L e g ξ f ξ + σ η 2 ξ + 2 ξ + N t N b θ } d η = 0 ,
where w 1 , w 2 , w 3 , w 4 , and w 5 are weight functions and maybe viewed as variations in f, g, θ , ϕ , ξ resp.

3.2. Finite Element Formulations

The equation of finite element model that is gained from the Equations (27)–(31) by exchanging the FEM approximation in the form:
f ¯ = j = 1 3 f ¯ j e ψ j e , g ¯ = j = 1 3 g ¯ j e ψ j e , θ ¯ = j = 1 3 θ ¯ j e ψ j e , ϕ ¯ = j = 1 3 ϕ ¯ j e ψ j e
with w 1 = w 2 = w 3 = w 4 = w 5 = ψ i e ( i = 1 , 2 , 3 ) , where the shape function ψ i e are the shape functions for the element ( η c , η c + 1 ) and are taken as
ψ 1 e = ( η c + 1 η c 2 η ) ( η c + 1 η ) ( η c + 1 η c ) 2 , ψ 2 e = 4 ( η η c ) ( η c + 1 η ) ( η c + 1 η c ) 2 , ψ 3 e = ( η c + 1 η c 2 η ) ( η η c ) ( η c + 1 η c ) 2 , η c η η c + 1 .
Therefore, the finite element model equations are given by
[ W 11 ] [ W 12 ] [ W 13 ] [ W 14 ] [ W 15 ] [ W 21 ] [ W 22 ] [ W 23 ] [ W 24 ] [ W 25 ] [ W 31 ] [ W 32 ] [ W 33 ] [ W 34 ] [ W 35 ] [ W 41 ] [ W 42 ] [ W 43 ] [ W 44 ] [ W 45 ] [ W 51 ] [ W 52 ] [ W 53 ] [ W 54 ] [ W 55 ] f g θ ϕ ξ = b 1 b 2 b 3 b 4 b 5
where [ W m n ] and b m (m,n = 1,2,3,4,5) are the matrices and are given as:
[ W i j 11 ] = η c η c + 1 ψ i e d ψ j e d η d η , W i j 12 = η c η c + 1 ψ i e ψ j e d η , [ W i j 13 ] = [ W i j 14 ] = 0 , [ W i j 15 ] = [ W i j 21 ] = 0 , [ W i j 22 ] = ( 1 + 1 β ) η c η c + 1 d ψ i e d η d ψ j e d η d η + η c η c + 1 f ¯ ψ i e d ψ j e d η d η η c η c + 1 g ¯ ψ i e ψ j e d η A η 2 η c η c + 1 ψ i e d ψ j e d η d η η c η c + 1 g ¯ ψ i e ψ j e d η A η c η c + 1 ψ i e ψ j e d η M η c η c + 1 ψ i e ψ j e d η k p η c η c + 1 ψ i e ψ j e d η [ W i j 23 ] = λ 1 η c η c + 1 ψ i e ψ j e d η , W i j 24 = λ 2 η c η c + 1 ψ i e ψ j e d η , W i j 25 = λ 3 η c η c + 1 ψ i e ψ j e d η ,
[ W i j 31 ] = 0 , [ W i j 32 ] = 0 , [ W i j 33 ] = ( 1 + R ) 1 P r η c η c + 1 d ψ i e d η d ψ j e d η d η η c η c + 1 g ¯ ψ i e ψ j e d η + η c η c + 1 f ¯ ψ i e d ψ j e d η d η A η 2 η c η c + 1 ψ i e d ψ j e d η d η 2 A η c η c + 1 ψ i e ψ j e d η + N b η c η c + 1 ϕ ¯ ψ i e d ψ j e d η d η + N t η c η c + 1 θ ¯ ψ i e d ψ j e d η d η + Q η c η c + 1 ψ i e ψ j e d η , [ W i j 34 ] = N d η c η c + 1 d ψ i e d η d ψ j e d η d η , [ W i j 35 ] = 0 ,
[ W i j 41 ] = [ W i j 42 ] = 0 , [ W i j 43 ] = S r η c η c + 1 d ψ i e d η d ψ j e d η d η , [ W i j 44 ] = 1 S c η c η c + 1 d ψ i e d η d ψ j e d η d η + η c η c + 1 f ¯ ψ i e d ψ j e d η d η η c η c + 1 g ¯ ψ i e ψ j e d η d η A η 2 η c η c + 1 ψ i e d ψ j e d η d η 2 A η c η c + 1 ψ i e ψ j e d η , [ W i j 45 ] = [ W i j 51 ] = [ W i j 52 ] = [ W i j 54 ] = 0 , [ W i j 53 ] = N t N b η c η c + 1 d ψ i e d η d ψ j e d η d η
[ W i j 55 ] = η c η c + 1 d ψ i e d η d ψ j e d η d η + L e η c η c + 1 f ¯ ψ i e d ψ j e d η d η L e η c η c + 1 g ¯ ψ i e ψ j e d η d η L e A η 2 η c η c + 1 ψ i e d ψ j e d η d η 2 L e A η c η c + 1 ψ i e ψ j e d η ,
and
b i 1 = 0 , b i 2 = ( ψ d g d η ) η c η c + 1 , b i 3 = ( 1 + R ) 1 P r ( ψ d θ d η ) η c η c + 1 N d ( ψ d ϕ d η ) η c η c + 1 , b i 4 = 1 S c ( ψ d ϕ d η ) η c η c + 1 S r ( ψ d θ d η ) η c η c + 1 , b i 5 = ( ψ d ξ d η ) η c η c + 1 N t N b ( ψ d θ d η ) η c η c + 1 ,
where f ¯ = j = 1 3 f ¯ j ψ j e , g ¯ = j = 1 3 g ¯ j ψ j e , θ ¯ = j = 1 3 θ ¯ j ψ j e , and ϕ ¯ = j = 1 3 ϕ ¯ j ψ j e are assumed to be known.
After the assemblage of element equations, a consequential system of non-linear equations is obtained, after that, it requires an iterative scheme to solve it for an efficient solution. The calculation of f, g, θ , and ϕ are then conceded out for a higher level, and then proceeding until the required 0.00005 is not obtained. The results in Table 1 show the convergency, as we computed, for the increasing number of elements, n = 300 , 500 , 1000 , 1300 , 1500 , 1800 , 1900 , 2100 . From the results, it is clear that as the number of elements increases further, no significant difference in the values of f, g, θ , ϕ , and ξ can be seen as the number of elements increases beyond 1800, so the outcomes at n = 1800 elements are reported.

4. Results and Discussion

The numerical calculations have been performed for the velocity, temperature, solutal, and Casson nano-fluid volume fraction functions for a different assessment of physical parameters, as magnetic parametre M, Prandtl number P r , Unsteadiness σ , permeability k p , Brownian motion N b , thermophoresis N t , thermal radiation R, Dufour N d , Schmidt S c , buoyancy λ 1 , λ 2 , λ 3 , Suction/Injection f w , Soret S r , Lewis number L e , hydro-dynamic slip S f , thermal slip S θ , solutal slip S ϕ , and nano-particles concentration slip S ξ .
The detail of the present results and the appraisal of flow velocity is made with the exact solution that is given by Crane [32] as f ( η ) = 1 e ( η ) under the special case ( M = 0 , σ = 0 , β , λ 1 = λ 2 = λ 3 = 0 , S f = 0 , f w = 0 , k p = 0 ). The finite element method’s outcomes have decent concurrence with the exact solution, which approves the validity of the finite element method. It can be seen clearly in Table 2, and in Table 3, that the skin friction coefficient attained by the finite element method is equivalent to the numerical outcomes of Gireesha et al. [33] and the exact solution of Mudassar et al. [34] under special case σ = 0 , β = 0 , λ 1 = λ 2 = λ 3 = 0 , S f = 0 , f w = 0 , k p = 0 . To confirm the accuracy of the presented numerical results, the results obtained by the finite element method for skin friction co-efficient for steady and unsteady flow have been compared with the numerical results that have previously been reported in studies and shown in Table 4. Regarding our results, there is admirable conformity among our outcomes and previously available research articles, which approves the cogency and the accuracy of the current results that are obtained by the finite element method (FEM). Table 5 describes the results of the heat transfer rate that are acquired by the finite element method, which are compared with the results of earlier studies and with the accurate solution of Ishak et al. [35] in a special case (Nt = Nb = 0). We observe that our results are in complete agreement and the grid invariance test has performed to sustain accuracy for 4 decimal points. In Table 6, local skin friction coefficient f ( 0 ) , the rates of heat transfer θ ( 0 ) , and mass transfer ϕ ( 0 ) that are acquired by FEM, which are also compared with the published research work, show excellent correlation.
Figure 2 shows the effect of M on the velocity function with no hydro-dynamic slip and with hydro-dynamic slip. The results show that in both cases the velocity component decreases as the value of the magnetic parameter M increases. Physically, the magnetic parameter M produced Lorentz force like a drag force, which causes the motion of the fluid to decelerate. Figure 2 also shows that in both cases the suction f w decelerates the momentum boundary layer thickness. The boundary layer thickness can be controlled with the help of suction f w .
Figure 3 and Figure 4 exemplify that the velocity boundary layer thickness increases with enhancing the thermal buoyancy parameter λ 1 and the nano-particles concentration buoyancy parameter λ 3 in the non-existence of hydro-dynamic slip and the existence of hydro-dynamic slip. In these two cases, we have to perceive that the momentum boundary layer is prolonged with the increasing values of buoyancy λ 1 and λ 3 . Also in Figure 3 and Figure 4, the velocity profile is increasing in the case of steady and unsteady flow and alike behavior of solutal buoyancy has been observed in Figure 5. It is observed that increasing the value of radiation parameter R is the reason to improve the velocity profile and incentive of instability in parameter k p , and the velocity profile is explained in Figure 6.
Also, we perceive in both cases that as the permeability of k p increases, this causes a decline in the velocity profile with no hydro-dynamic slip and with hydro-dynamic slip. Furthermore, we also perceive in Figure 6 that the suction decreases the momentum boundary layer thickness. Similar behavior of the Deborah number β has been perceived on the velocity profile in Figure 7. Figure 8 indicates the impact of M on the temperature profile with the non-existence of thermal slip and with existence thermal slip. In the figure, it is clear that with an increasing value of M the temperature profile also increases in both cases. From Figure 8 it is also clear that suction f w decreases the thermal boundary layer thickness in both cases and the thickness of the thermal boundary layer can be controlled with the help of suction f w . Figure 9 depicts the similar behavior of M. It is perceived that the temperature increases with the discrepancy in parameter R.
Figure 10 demonstrates the influence of σ and P r on the temperature profile with and without the existence of thermal slip, with increasing P r , σ , and S θ causing the thermal boundary layer to decrease. The description of the effect of N b , f w , and S θ with the temperature profile is explained in Figure 11. It is clear that with the increase in N b , the thermal boundary layer of fluid flow increases in both cases. It can also be perceived that f w decreases the thermal boundary layer thickness in both cases. In Figure 12, there is the influence of M and S c on the solutal profile with the non-existence of solutal slip and with the existence of solutal slip. From the figure, it is clear that with the increasing value of M, the solutal profile also increases in both cases. Furthermore, the effect of S c over the solutal profile declines the profile. The impact of the existence and non-existence of solutal slip and increasing values of S r is depicted in Figure 13. It is seen that increasing the value of the Soret parameter increases the solutal profile in both the steady and unsteady cases. It is also seen that increasing the value of the slip results in decreasing the thickness of the solutal profile boundary.
Figure 14 indicates the impact of N t on the nano-particle volume fraction profile without the existence of nano-particles concentration slip and with the existence of nano-particles concentration slip. It is clear from the figure that in both the steady and unsteady cases, increasing values of N t cause a decline in the nano-particle volume fraction profile. The impacts of σ and S γ over the nano-particle volume fraction profile have been illustrated in Figure 14. It is noticed that the nano-particle volume fraction profile decreases in increments in the values of σ and S γ . Figure 15 and Figure 16 show the impact of Lewis number L e , unsteadiness parameter σ , and suction parameter f w on the nano-particle volume fraction profile, and interpretations about the nano-particle volume fraction profile are made. It is seen that increasing values of L e and f w in the boundary layer of nano-particles causes the volume fraction profile to decline. In addition, it is clear from the figures that the nano-particle volume fraction profile decreases both with and without the existence of nano-particles concentration slip.
Figure 17 shows the impact of M, σ , and S f on the skin friction co-efficient and also shows that the skin friction co-efficient decreases with the increase of slip parameter, magnetic and unsteadiness parameter. In Figure 18, the skin friction coefficient increases with the increment in the thermal buoyancy, suction parameters, and solutal buoyancy. Figure 19 represents the increment in the value of Nusselt number with increasing values of magnetic, radiation and thermal slip parameters. As perceive in Figure 20 that when the values of solutal buoyancy and suction intensity increases with the existence and the non-existence absence of thermal buoyancy. Figure 21 demonstrates that the impact of S c , S r , and S ϕ on the reduced Sherwood number. Also from the Figure 22, it clear that the reduced Sherwood number decreases with increment in the values of the unsteadiness parameter σ , magnetic parameter M and the suction parameters f w .
In Table 7 we analyze the variation of physical parameters M, λ 1 , λ 2 , λ 3 , σ , P r , L e on skin friction co-efficient f ( 0 ) , Nusselt number θ ( 0 ) , and Sherwood number ξ ( 0 ) .
The following results are concluded from Table 7:
(i) The skin-friction coefficient is increasing while reducing the local Nusselt and Sherwood numbers through improvement in the Magnetic parameter; (ii) The increment in thermal buoyancy parameters λ 1 , λ 2 , λ 3 causes the skin-friction coefficient to decrease while increasing the local Nusselt number and Sherwood number; (iii) With the increasing unsteadiness parameter σ , the skin-friction coefficient, local Nusselt, and Sherwood numbers are also increasing; (iv) The skin-friction coefficient is increasing with the increment in Prandtl number and the local Nusselt number and Sherwood number also increase; (v) The skin-friction coefficient is increasing with the increasing Lewis number and also increment in the local Nusselt number and Sherwood number.

5. Conclusions

In this paper, we have developed a mathematical model to investigate the unsteady two-dimensional magnetohydrodynamic (MHD) flow and heat transfer of an incompressible electrically conducted fluid. This study was conducted to analyze the multiple slip effects on magneto-hydrodynamic unsteady Casson nano-fluid flow over a non-linear porous shrinking/stretching sheet in the presence of a heat source with Soret effect. Appropriate similarities have been used for transformations and the governing partial differential equations (PDEs) are rendered into a system of ordinary differential equations (ODEs). The resulting system of ordinary differential equations was solved numerically with an efficient and validated finite element method (FEM). We also use a special case of the present model to compare our results with previous studies. A parametric study has been performed to explore the mass and heat transfer characteristics, and also the impact of different parameters, of the flow. The following are the results that can be concluded from the present study.
  • The cause of reducing the fluid velocity near the region of the boundary layer is the increment in values of slip, suction, magnetic field, and unsteady parameters.
  • The increment in thermal buoyancy parameters λ 1 , λ 2 , λ 3 and slip parameters cause decreasing the skin-friction coefficient while increasing the local Nusselt number and Sherwood number.
  • Decreasing the value of the Sherwood number causes an increasing value of the solutal slip parameter, Schmidt number, and Soret number, but the effect is the opposite with increasing values of the unsteadiness, magnetic, and suction parameters.
  • With the increasing value of the magnetic parameter M, the slip velocity parameter and suction parameter are found to be reduced in the velocity profile.
  • The skin-friction coefficient decreases with the increasing value of slip, magnetic, and unsteadiness parameters, but the effect is the opposite for increasing values of thermal buoyancy, suction parameter, and solutal buoyancy.
  • The increment of the values of buoyancy parameter, injection parameter, and heat source parameter cause an increase in the velocity profile.

Author Contributions

L.A. modeled the problem and wrote the manuscript. X.L. thoroughly checked the mathematical modeling and English corrections. B.A. helped in MATLAB coding. S.M., and S.A. writing—review and editing. X.L. contributed to the results and discussions. All authors finalized the manuscript after its internal evaluation.

Funding

This work was supported by the National Natural Science Foundation of China (No. 51676152), Equipment Advance Research Field Foundation (No. 61402070302), and Fundamental Research Funds for the Central Universities (No. zrzd2017012).

Acknowledgments

The first author acknowledge with thanks the Liu Xiaomin for his continuous guidance throughout the work at Xian Jiaotong University China, National Natural Science Foundation of China, Equipment Advance Research Field Foundation, Fundamental Research Funds for the Central Universities for technical and financial support for this research project. Authors would also like to ac- knowledge and express their gratitude to the Chinese Government Scholarship Council (CSC) for the scholarship award.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Tretheway, D.C.; Meinhart, C.D. Apparent fluid slip at hydrophobic microchannel walls. Phys. Fluids 2002, 14, L9–L12. [Google Scholar] [CrossRef]
  2. Soltani, F.; Yilmazer, L. Slip velocity and slip layer thickness in flow of concentrated suspensions. J. Appl. Polym. Sci. 2010, 70, 515–522. [Google Scholar] [CrossRef]
  3. Mabood, F.; Das, K. Melting heat transfer on hydromagnetic flow of a nanofluid over a stretching sheet with radiation and second-order slip. Eur. Phys. J. Plus 2016, 131, 1–12. [Google Scholar] [CrossRef]
  4. Pop, I.; Na, T.Y. Unsteady flow past a stretching sheet. Mech. Res. Commun. 1996, 23, 413–422. [Google Scholar] [CrossRef]
  5. Sharidan, S.; Mahmood, M.; Pop, I. Similarity solutions for the unsteady boundary layer flow and heat transfer due to a stretching sheet. Appl. Mech. Eng. 2006, 11, 647. [Google Scholar]
  6. Bhattacharyya, K.; Mukhopadhyay, S.; Layek, G. Slip effects on an unsteady boundary layer stagnation-point flow and heat transfer towards a stretching sheet. Chin. Phys. Lett. 2011, 28, 094702. [Google Scholar] [CrossRef]
  7. Khan, U.; Ahmed, N.; Mohyud-Din, S.T. Heat transfer effects on carbon nanotubes suspended nanofluid flow in a channel with non-parallel walls under the effect of velocity slip boundary condition: A numerical study. Neural Comput. Appl. 2017, 28, 37–46. [Google Scholar] [CrossRef]
  8. Mohyud-Din, S.T.; Zaidi, Z.A.; Khan, U.; Ahmed, N. On heat and mass transfer analysis for the flow of a nanofluid between rotating parallel plates. Aerosp. Sci. Technol. 2015, 46, 514–522. [Google Scholar] [CrossRef]
  9. Mohyud-Din, S.; Khan, U.; Ahmed, N.; Hassan, S. Magnetohydrodynamic flow and heat transfer of nanofluids in stretchable convergent/divergent channels. Appl. Sci. 2015, 5, 1639–1664. [Google Scholar] [CrossRef]
  10. Makinde, O.; Mishra, S. On stagnation point flow of variable viscosity nanofluids past a stretching surface with radiative heat. Int. J. Appl. Comput. Math. 2017, 3, 561–578. [Google Scholar] [CrossRef]
  11. Sharada, K.; Shankar, B. MHD mixed convection flow of a Casson fluid over an exponentially stretching surface with the effects of Soret, Dufour, thermal radiation and chemical reaction. World J. Mech. 2015, 5, 165. [Google Scholar] [CrossRef]
  12. Ullah, I.; Bhattacharyya, K.; Shafie, S.; Khan, I. Unsteady MHD mixed convection slip flow of Casson fluid over nonlinearly stretching sheet embedded in a porous medium with chemical reaction, thermal radiation, heat generation/absorption and convective boundary conditions. PLoS ONE 2016, 11, e0165348. [Google Scholar] [CrossRef] [PubMed]
  13. Ullah, I.; Khan, I.; Shafie, S. Hydromagnetic Falkner-Skan flow of Casson fluid past a moving wedge with heat transfer. Alex. Eng. J. 2016, 55, 2139–2148. [Google Scholar] [CrossRef]
  14. Pushpalatha, K.; Sugunamma, V.; Reddy, J.R.; Sandeep, N. Heat and mass transfer in unsteady MHD Casson fluid flow with convective boundary conditions. Int. J. Adv. Sci. Technol. 2016, 91, 19–38. [Google Scholar] [CrossRef]
  15. Pushpalatha, K.; Reddy, J.R.; Sugunamma, V.; Sandeep, N. Numerical study of chemically reacting unsteady Casson fluid flow past a stretching surface with cross diffusion and thermal radiation. Open Eng. 2017, 7, 69–76. [Google Scholar] [CrossRef]
  16. Uddin, M.J.; Khan, W.; Ismail, A.M. Free convection boundary layer flow from a heated upward facing horizontal flat plate embedded in a porous medium filled by a nanofluid with convective boundary condition. Transp. Porous Media 2012, 92, 867–881. [Google Scholar] [CrossRef]
  17. Rahman, M.; Al-Lawatia, M.; Eltayeb, I.; Al-Salti, N. Hydromagnetic slip flow of water based nanofluids past a wedge with convective surface in the presence of heat generation (or) absorption. Int. J. Therm. Sci. 2012, 57, 172–182. [Google Scholar] [CrossRef]
  18. Das, K. Slip flow and convective heat transfer of nanofluids over a permeable stretching surface. Comput. Fluids 2012, 64, 34–42. [Google Scholar] [CrossRef]
  19. Ibrahim, W.; Shankar, B. MHD boundary layer flow and heat transfer of a nanofluid past a permeable stretching sheet with velocity, thermal and solutal slip boundary conditions. Comput. Fluids 2013, 75, 1–10. [Google Scholar] [CrossRef]
  20. Abbas, Z.; Wang, Y.; Hayat, T.; Oberlack, M. Slip effects and heat transfer analysis in a viscous fluid over an oscillatory stretching surface. Int. J. Numer. Methods Fluids 2009, 59, 443–458. [Google Scholar] [CrossRef]
  21. Makinde, O.D.; Mabood, F.; Ibrahim, M.S. Chemically reacting on MHD boundary-layer flow of nanofluids over a non-linear stretching sheet with heat source/sink and thermal radiation. Therm. Sci. 2018, 22, 495–506. [Google Scholar] [CrossRef]
  22. Ibrahim, S.; Kumar, P.; Lorenzini, G.; Lorenzini, E.; Mabood, F. Numerical study of the onset of chemical reaction and heat source on dissipative MHD stagnation point flow of Casson nanofluid over a nonlinear stretching sheet with velocity slip and convective boundary conditions. J. Eng. Thermophys. 2017, 26, 256–271. [Google Scholar] [CrossRef]
  23. Prasannakumara, B.; Krishnamurthy, M.; Gireesha, B.; Gorla, R.S. Effect of multiple slips and thermal radiation on MHD flow of Jeffery nanofluid with heat transfer. J. Nanofluids 2016, 5, 82–93. [Google Scholar] [CrossRef]
  24. Imtiaz, M.; Hayat, T.; Alsaedi, A.; Hobiny, A. Homogeneous-heterogeneous reactions in MHD flow due to an unsteady curved stretching surface. J. Mol. Liq. 2016, 221, 245–253. [Google Scholar] [CrossRef]
  25. Mabood, F.; Shateyi, S. Multiple Slip Effects on MHD Unsteady Flow Heat and Mass Transfer Impinging on Permeable Stretching Sheet with Radiation. Model. Simul. Eng. 2019, 2019, 1–11. [Google Scholar] [CrossRef]
  26. Malarselvi, A.; Bhuvaneswari, M.; Sivasankaran, S.; Ganga, B.; Hakeem, A.A. Effect of Slip and Convective Heating on Unsteady MHD Chemically Reacting Flow Over a Porous Surface with Suction. In Applied Mathematics and Scientific Computing; Springer: Cham, Switzerland, 2019; pp. 357–365. [Google Scholar]
  27. Singh, K.; Kumar, M. Effect of viscous dissipation on double stratified MHD free convection in micropolar fluid flow in porous media with chemical reaction, heat generation and ohmic heating. Chem. Process Eng. Res. 2015, 31, 75–80. [Google Scholar]
  28. Rana, P.; Bhargava, R.; Bég, O.A. Finite element simulation of unsteady magneto-hydrodynamic transport phenomena on a stretching sheet in a rotating nanofluid. Proc. Inst. Mech. Eng. Part N J. Nanoeng. Nanosyst. 2013, 227, 77–99. [Google Scholar] [CrossRef]
  29. Reddy, J.N. An Introduction to the Finite Element Method; McGraw-Hill Education: New York, NY, USA, 1993. [Google Scholar]
  30. Swapna, G.; Kumar, L.; Rana, P.; Singh, B. Finite element modeling of a double-diffusive mixed convection flow of a chemically-reacting magneto-micropolar fluid with convective boundary condition. J. Taiwan Inst. Chem. Eng. 2015, 47, 18–27. [Google Scholar] [CrossRef]
  31. Gupta, D.; Kumar, L.; Beg, O.A.; Singh, B. Finite-element analysis of transient heat and mass transfer in microstructural boundary layer flow from a porous stretching sheet. Comput. Therm. Sci. Int. J. 2014, 6. [Google Scholar] [CrossRef]
  32. Crane, L.J. Flow past a stretching plate. Z. Für Angew. Math. Und Phys. ZAMP 1970, 21, 645–647. [Google Scholar] [CrossRef]
  33. Gireesha, B.; Ramesh, G.; Bagewadi, C. Heat transfer in MHD flow of a dusty fluid over a stretching sheet with viscous dissipation. J. Appl. Sci. Res. 2012, 3, 2392–2401. [Google Scholar]
  34. Jalil, M.; Asghar, S.; Yasmeen, S. An exact solution of MHD boundary layer flow of dusty fluid over a stretching surface. Math. Probl. Eng. 2017, 2017. [Google Scholar] [CrossRef]
  35. Ishak, A.; Nazar, R.; Pop, I. Boundary layer flow and heat transfer over an unsteady stretching vertical surface. Meccanica 2009, 44, 369–375. [Google Scholar] [CrossRef]
  36. Uddin, M.; Rana, P.; Bég, O.A.; Ismail, A.M. Finite element simulation of magnetohydrodynamic convective nanofluid slip flow in porous media with nonlinear radiation. Alex. Eng. J. 2016, 55, 1305–1319. [Google Scholar] [CrossRef]
  37. Chamkha, A.; Aly, A.; Mansour, M. Similarity solution for unsteady heat and mass transfer from a stretching surface embedded in a porous medium with suction/injection and chemical reaction effects. Chem. Eng. Commun. 2010, 197, 846–858. [Google Scholar] [CrossRef]
  38. Ali, M.E. Heat transfer characteristics of a continuous stretching surface. Wärme-Und Stoffübertragung 1994, 29, 227–234. [Google Scholar] [CrossRef]
  39. Pal, D. Combined effects of non-uniform heat source/sink and thermal radiation on heat transfer over an unsteady stretching permeable surface. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 1890–1904. [Google Scholar] [CrossRef]
  40. Haile, E.; Shankar, B. Heat and mass transfer in the boundary layer of unsteady viscous nanofluid along a vertical stretching sheet. J. Comput. Eng. 2014, 2014. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Physical sketch with xy-coordinates.
Figure 1. Physical sketch with xy-coordinates.
Applsci 09 05217 g001
Figure 2. Effects of M, f w , and S f on f .
Figure 2. Effects of M, f w , and S f on f .
Applsci 09 05217 g002
Figure 3. Effects of λ 1 , σ , and S f = 0 on f .
Figure 3. Effects of λ 1 , σ , and S f = 0 on f .
Applsci 09 05217 g003
Figure 4. Effects of λ 3 and σ on f . (a) S f = 0 (b) S f = 0.1 .
Figure 4. Effects of λ 3 and σ on f . (a) S f = 0 (b) S f = 0.1 .
Applsci 09 05217 g004
Figure 5. Effects of λ 2 and R on f . (a) S f = 0 (b) S f = 1 .
Figure 5. Effects of λ 2 and R on f . (a) S f = 0 (b) S f = 1 .
Applsci 09 05217 g005
Figure 6. Effects of k p , f w , and S f on f .
Figure 6. Effects of k p , f w , and S f on f .
Applsci 09 05217 g006
Figure 7. Effects of β , f w , and S f on f .
Figure 7. Effects of β , f w , and S f on f .
Applsci 09 05217 g007
Figure 8. Effects of M, f w , and S θ on θ .
Figure 8. Effects of M, f w , and S θ on θ .
Applsci 09 05217 g008
Figure 9. Effects of M, R, and S θ on θ .
Figure 9. Effects of M, R, and S θ on θ .
Applsci 09 05217 g009
Figure 10. Effects of σ , P r , and S θ on θ .
Figure 10. Effects of σ , P r , and S θ on θ .
Applsci 09 05217 g010
Figure 11. Effects of N b , f w , and S θ on θ .
Figure 11. Effects of N b , f w , and S θ on θ .
Applsci 09 05217 g011
Figure 12. Effects of M, S c , and S ϕ on ϕ .
Figure 12. Effects of M, S c , and S ϕ on ϕ .
Applsci 09 05217 g012
Figure 13. Effects of S r , σ , and S ϕ on ϕ .
Figure 13. Effects of S r , σ , and S ϕ on ϕ .
Applsci 09 05217 g013
Figure 14. Effects of N t , σ , and S ξ on ξ .
Figure 14. Effects of N t , σ , and S ξ on ξ .
Applsci 09 05217 g014
Figure 15. Effects of L e , f w , and S ξ on ξ .
Figure 15. Effects of L e , f w , and S ξ on ξ .
Applsci 09 05217 g015
Figure 16. Effects of L e , σ , and S ξ on ξ .
Figure 16. Effects of L e , σ , and S ξ on ξ .
Applsci 09 05217 g016
Figure 17. Effects of M, σ , and S f on the skin friction coefficient.
Figure 17. Effects of M, σ , and S f on the skin friction coefficient.
Applsci 09 05217 g017
Figure 18. Effects of f w , λ 1 , and λ 2 on the skin friction coefficient.
Figure 18. Effects of f w , λ 1 , and λ 2 on the skin friction coefficient.
Applsci 09 05217 g018
Figure 19. Effects of M, R, and S θ on the reduced nusselt number.
Figure 19. Effects of M, R, and S θ on the reduced nusselt number.
Applsci 09 05217 g019
Figure 20. Effects of f w , λ 1 , and λ 2 on the reduced nusselt number.
Figure 20. Effects of f w , λ 1 , and λ 2 on the reduced nusselt number.
Applsci 09 05217 g020
Figure 21. Effects of S c , S r , and S ϕ on the reduced Sherwood number.
Figure 21. Effects of S c , S r , and S ϕ on the reduced Sherwood number.
Applsci 09 05217 g021
Figure 22. Effects of M, σ , and f w on the reduced Sherwood number.
Figure 22. Effects of M, σ , and f w on the reduced Sherwood number.
Applsci 09 05217 g022
Table 1. The convergence of the results of finite element method (FEM) when P r = 1 , M = σ = 0.2 , β = 0.1 , λ 1 = λ 2 = λ 3 = 0.2 , S f = S θ = S ϕ = S ξ = S t = S r = R = S g = S p = 0.5 , k P = N u = N t = N b = 0.1 , S c = 10 , L e = 5 , f w = 0.3 .
Table 1. The convergence of the results of finite element method (FEM) when P r = 1 , M = σ = 0.2 , β = 0.1 , λ 1 = λ 2 = λ 3 = 0.2 , S f = S θ = S ϕ = S ξ = S t = S r = R = S g = S p = 0.5 , k P = N u = N t = N b = 0.1 , S c = 10 , L e = 5 , f w = 0.3 .
Number of Elements f ( 3 ) g ( 3 ) θ ( 3 ) ϕ ( 3 ) ξ ( 3 )
3000.824390.162410.167330.068760.03046
5000.824550.162500.197290.068750.03046
10000.824660.162570.197260.068740.03045
13000.824690.162590.197260.068730.03045
15000.824700.162600.197250.068730.03045
18000.824710.162600.197250.068730.03045
19000.824720.162610.197250.068740.03045
21000.824720.162610.197250.068740.03045
Table 2. The comparison of the exact solution of Crane and FEM for the flow velocity f ( η ) .
Table 2. The comparison of the exact solution of Crane and FEM for the flow velocity f ( η ) .
η Crane [32]FEM η Crane [32]FEM
01170.00090.0009
10.13790.137980.00030.0003
20.13530.135390.00010.0001
30.04980.0498100.00000.0000
40.01830.0018110.00000.0000
50.00670.0067120.00000.0000
60.00250.0025130.00000.0000
Table 3. Assessment of skin friction coefficient for different values of M.
Table 3. Assessment of skin friction coefficient for different values of M.
MGireesha et al. [33]Mudassar et al. [34]FEM (Our Results)Error in %
β = 0 Exact Solution (a)(b) | ( b a a ) | × 100
0.01.0001.0000001.00000780.00078
0.21.0951.0954451.09544620.00010
0.51.2241.2247451.22474520.00001
1.01.4141.4142141.41421420.00001
1.21.4831.4832401.48323850.00010
1.51.5811.5811391.58113920.00001
2.01.7321.7320511.73205150.00002
Table 4. Assessment of skin friction coefficient for various values of M and σ when f w = λ 1 = λ 2 = λ 3 = k p = S f = 0 .
Table 4. Assessment of skin friction coefficient for various values of M and σ when f w = λ 1 = λ 2 = λ 3 = k p = S f = 0 .
MMabood and DasFazle and StanfordFEM σ Chamkha et al.Fazle and StanfordFEM
[36] [25] (Our Results) [37] [25] (Our Results)
0−1.000008−1.0000084−1.00000820.2--1.068027
11.41421351.414213561.414213530.4--1.134687
52.44948972.449489742.449489630.6--1.199118
103.31662473.316624793.316624630.81.2615121.2610421.261042
507.14142847.141428437.141428391.0--1.320522
10010.04987510.049875610.04987511.21.3780521.3777241.377724
50022.38302922.383029322.38302831.4--1.432836
100031.63858431.638584031.63858331.6--1.486039
Table 5. Comparison of θ ( 0 ) for various values of P r when M = f w = S θ = σ = λ 1 = λ 2 = λ 3 = R = N d = k p = S f = 0 .
Table 5. Comparison of θ ( 0 ) for various values of P r when M = f w = S θ = σ = λ 1 = λ 2 = λ 3 = R = N d = k p = S f = 0 .
Pr AliFazle and StanfordIshak et al.Dulal Pal.Haile et al.Ishak et al. [35]FEMError in %
[38] [25] [35] [39][40] (a) Exact Solution(b) Our Results | ( b a a ) | × 100
0.720.80580.8088---0.80863134980.80863392990.0003
1.000.96911.00001.00001.00001.00041.0000000001.00000802130.0008
3.001.91441.92371.92371.92361.92341.9236825941.92367772210.0003
10.03.70063.72073.72073.72073.72053.7206739013.72066816830.0002
100--12.294112.294012.296212.29408326012.2940516590.0003
Table 6. Comparison of f ( 0 ) , θ ( 0 ) , and ϕ ( 0 ) for different values of Pr, and Sc when M = σ = λ 3 = β = R = N b = N t = N d = S r = 0 , f w = 0.5 , λ 1 = k p = 1 , λ 2 = 2 , S f = S θ = S ϕ = S ξ = 0 .
Table 6. Comparison of f ( 0 ) , θ ( 0 ) , and ϕ ( 0 ) for different values of Pr, and Sc when M = σ = λ 3 = β = R = N b = N t = N d = S r = 0 , f w = 0.5 , λ 1 = k p = 1 , λ 2 = 2 , S f = S θ = S ϕ = S ξ = 0 .
Pr Sc Chamkha et al. [37]FEM(Our Results)
f ( 0 ) , θ ( 0 ) , ϕ ( 0 ) f ( 0 ) θ ( 0 ) ϕ ( 0 )
0.710.220.273771.1583930.7401400.273711.1583890.740160
0.600.496771.1010871.3122480.496691.1010821.312192
0.940.599411.0769311.7134430.599351.0768991.713358
0.30.620.386390.6219801.3550450.386550.6218921.355036
0.71 0.504381.0990321.3375710.504291.0990181.337267
1.0 0.552441.3848751.3310270.552371.3847531.331856
3.0 0.690422.9669991.3167100.690362.9668941.316382
Table 7. Various mathematical values of physical constraints M, λ 1 , λ 2 , λ 3 , σ , P r , L e when R = f w = S r = S t = S f = S p = S g = 0.5 , k p = 0.1 , N u = N b = N t = Q = 0.1 , β = 0.2 , S c = 10 .
Table 7. Various mathematical values of physical constraints M, λ 1 , λ 2 , λ 3 , σ , P r , L e when R = f w = S r = S t = S f = S p = S g = 0.5 , k p = 0.1 , N u = N b = N t = Q = 0.1 , β = 0.2 , S c = 10 .
M λ 1 λ 2 λ 3 σ Pr Le f ( 0 ) θ ( 0 ) ξ ( 0 )
0.50.20.20.20.5150.435640.708631.33066
1.0 0.476950.701921.32746
1.5 0.512390.696121.32468
0.50.20.20.20.5150.435640.708631.33066
0.5 0.422230.710521.33163
0.8 0.409010.712371.33258
0.50.20.20.20.5150.435640.708631.33066
0.5 0.430950.709281.33099
0.8 0.426280.709931.33133
0.5o.20.20.20.5150.435640.708631.33066
0.5 0.432170.709061.33089
0.8 0.428700.709491.33113
0.50.20.20.20.5150.435640.708631.33066
1.0 0.469260.795561.37699
1.5 0.499030.860941.41349
0.50.20.20.20.5150.435640.708631.33066
3 0.438190.888131.28962
5 0.439251.002881.25652
0.50.20.20.20.5150.435640.708631.33066
70.436350.709191.44132
100.436890.709931.54699

Share and Cite

MDPI and ACS Style

Ali, L.; Liu, X.; Ali, B.; Mujeed, S.; Abdal, S. Finite Element Analysis of Thermo-Diffusion and Multi-Slip Effects on MHD Unsteady Flow of Casson Nano-Fluid over a Shrinking/Stretching Sheet with Radiation and Heat Source. Appl. Sci. 2019, 9, 5217. https://0-doi-org.brum.beds.ac.uk/10.3390/app9235217

AMA Style

Ali L, Liu X, Ali B, Mujeed S, Abdal S. Finite Element Analysis of Thermo-Diffusion and Multi-Slip Effects on MHD Unsteady Flow of Casson Nano-Fluid over a Shrinking/Stretching Sheet with Radiation and Heat Source. Applied Sciences. 2019; 9(23):5217. https://0-doi-org.brum.beds.ac.uk/10.3390/app9235217

Chicago/Turabian Style

Ali, Liaqat, Xiaomin Liu, Bagh Ali, Saima Mujeed, and Sohaib Abdal. 2019. "Finite Element Analysis of Thermo-Diffusion and Multi-Slip Effects on MHD Unsteady Flow of Casson Nano-Fluid over a Shrinking/Stretching Sheet with Radiation and Heat Source" Applied Sciences 9, no. 23: 5217. https://0-doi-org.brum.beds.ac.uk/10.3390/app9235217

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop