Next Article in Journal
Facile Preparation of Bilayer Titanium Silicate (TS-1) Zeolite Membranes by Periodical Secondary Growth
Next Article in Special Issue
Impact of Velocity Second Slip and Inclined Magnetic Field on Peristaltic Flow Coating with Jeffrey Fluid in Tapered Channel
Previous Article in Journal
Effect of Trajectory Curvature on the Microstructure and Properties of Surfacing Wall Formed with the Process of Wire Arc Additive Manufacturing
Previous Article in Special Issue
Convective Heat Transfer and Magnetohydrodynamics across a Peristaltic Channel Coated with Nonlinear Nanofluid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Second-Order Slip and Double Stratification Coatings on 3D MHD Williamson Nanofluid Flow with Cattaneo–Christov Heat Flux

1
Department of Computer Science, Bahria University, Islamabad 44000, Pakistan
2
Department of Mechanical Engineering, Sejong University, Seoul 143-747, Korea
3
Department of Mathematics and Computer Science, Faculty of Science, Beirut Arab University, Beirut 11072809, Lebanon
4
Department of ICT Convergence Rehabilitation Engineering, Soonchunhyang University, Asan 31538, Korea
5
Department of Computer Science and Engineering, Soonchunhyang University, Asan 31538, Korea
6
Faculty of Science, Jiangsu University, Zhenjiang 212013, China
*
Author to whom correspondence should be addressed.
Coatings 2019, 9(12), 849; https://doi.org/10.3390/coatings9120849
Submission received: 1 November 2019 / Revised: 5 December 2019 / Accepted: 9 December 2019 / Published: 11 December 2019
(This article belongs to the Special Issue Recent Trends in Coatings and Thin Film–Modeling and Application)

Abstract

:
The present research examines the impact of second-order slip with thermal and solutal stratification coatings on three-dimensional (3D) Williamson nanofluid flow past a bidirectional stretched surface and envisages it analytically. The novelty of the analysis is strengthened by Cattaneo–Christov (CC) heat flux accompanying varying thermal conductivity. The appropriate set of transformations is implemented to get a differential equation system with high nonlinearity. The structure is addressed via the homotopy analysis technique. The authenticity of the presented model is verified by creating a comparison with the limited published results and finding harmony between the two. The impacts of miscellaneous arising parameters are deliberated through graphical structures. Some useful tabulated values of arising parameters versus physical quantities are also discussed here. It is observed that velocity components exhibit an opposite trend with respect to the stretching ratio parameter. Moreover, the Brownian motion parameter shows the opposite behavior versus temperature and concentration distributions.

1. Introduction

The heat transfer phenomenon plays a vital role when the temperature varies between different bodies or parts of the same body. Heat can be transferred by three methods: convection of fluids, conduction in solids, and radiation. To examine body heat transfer, the principle of heat transfer can be applied to the human body. Here, we can quote one example: the metal pan is used to transfer heat from the stove to food. Some applications of heat transfer are cooking food over metal pots, boiling milk in metal pots, and thermal treatment of pain by a hot water bag. Straughan [1] considered the Cattaneo–Christov (CC) model for heat flux and thermal convection over a Newtonian fluid. Khan et al. [2] analyzed it numerically by engaging the bvp4c MATLAB-based function on a Sisko fluid flow accompanied by generalized Fick’s and Fourier’s laws over a nonlinear stretched surface. Hayat et al. [3] examined analytically the magnetohydrodynamics (MHD) flow of Jeffrey fluid past a variable thick surface via the impacts of the chemical reaction and the CC model in a stratified medium. Waqas et al. [4] investigated the CC heat flux model for energy equation formulation rather than Fourier’s law of heat conduction. It was observed that a variable thermal conductivity remains inversely proportional to a temperature profile. The Soret–Dufour effects on walls with the second-grade fluid flow between inclined parallel plain walls were inspected by Khan et al. [5]. Heat transfer characteristics of the incompressible flow of the second-grade fluid flow produced by a stretching sheet were analyzed by Ghadikolaei et al. [6]. Khan et al. [7] investigated the behavior of homogenous–heterogeneous reactions against heat transfer flow due to a stretching sheet. They noticed that homogenous–heterogeneous reactions reduce fluid concentration.
In modern engineering processes, especially in metallurgical engineering and metalworking practices, the role of MHD is fundamental for electrically conducting fluids. The magnetic field function is crucial in cooling the hot plasma inside a nuclear reactor vessel. Similarly, the magnetic field is employed for the mixing of metals inside an electrical furnace [8]. Chamkha et. al. [9] analyzed the magnetic field effect on the mixed convection unsteady flow in an ambient fluid past a cone rotating with an unsteady angular velocity. Pullepu et al. [10] analyzed the free convection flow with variable surface temperature over a nonisothermal vertical cone. Akbar et al. [11] examined the two dimensional (2D) electrically conducting flow of the hyperbolic tangent fluid past a stretching surface. They observed that an increment in the Hartmann number decelerates the fluid velocity in the domain of the stretching sheet. Seini et al. [12] highlighted the magnetic field impact over a stretching surface accompanied by the appearance of slip velocity near the stagnation point flow. They perceived that the impact of the magnetic field is more significant on the velocity profile. Ravindran et al. [13] considered the impact of a transverse magnetic field and heat generation and absorption on time-dependent mixed convection flow over a porous cone with a chemical reaction. Boland et al. [14] simulated MHD flow of viscous fluid over a circular cylinder covered with a permeable layer. They adopted the Darcy–Brinkman–Forchheimer model to study the flow inside a porous medium. Ellahi et al. [15] investigated the influence of a Hall current on MHD Jeffrey fluid flow over a nonuniform duct. Mishra et al. [16] explored heat transfer and mass in the appearance of a magnetic field of viscoelastic fluid flow. They determined that the behavior of the magnetic field against the velocity profile is opposite to temperature distribution and concentration profiles. Hussain et al. [17] analyzed numerically the influence of the applied magnetic field on a non-Newtonian fluid flow past a stretching surface.
Nanofluids are vital in many engineering applications, such as in biomedical engineering and many chemical processes. Nanofluid is composed of nanometer-sized particles with a diameter of less than 100 nm and some conventional fluid. The basic aim of using nanofluids is to upgrade the heat transfer and thermal conductivity to attain better cooling. Khan et al. [18] highlighted nanofluid flow over a stretching sheet. Makinde [19] extended the work of Khan et al. [18] to convective boundary condition in nanofluid flow. Nadeem et al. [20] analyzed the second-grade nanofluid (nonorthogonal stagnation point) flow in the direction of the stretching surface. The impact of the variable magnetic field on the nanofluid flow between two disks was explored by Hayami et al. [21]. They found the analytical solution via the homotopy perturbation method and observed that the temperature of the boundary layer thickness decreases with the increase of the Brownian motion parameter and thermophoretic factor. Nanofluid flow in a permeable medium over a convectively heated permeable shrinking sheet was examined by Hayat et al. [22]. Sheikholeslami et al. [23,24] considered the behavior of the magnetic field on the free and forced convection flow of nanofluids respectively by making use of the two-phase model. Hassan et al. [25] elaborated on convective transport of heat transfer in a nanofluid through a porous medium. They concluded that convective heat transfer is improved by nanoparticle concentration, and the magnetic field impacts second-order slip flow. Nayak et al. [26] focused on the numerical solution of the three-dimensional (3D) nanofluid flow with nonlinear thermal radiation with convective conditions and slip. Hosseini et al. [27] observed the nanofluid MHD flow in a microchannel heat sink via the KKL (Koo–Kleinsteuer–Li) model. They noticed that the interaction between nanoparticles and the solid phase enhances the Nusselt number. In recent years, several scientists have used nanofluid heat transfer in their studies [28,29,30,31,32,33,34].
Fluids are basically divided into two groups: Newtonian [35] and non-Newtonian fluids [36]. Fluids that abide by Newton’s law of viscosity are termed Newtonian fluids. However, a contradiction to Newton’s law of viscosity is observed in non-Newtonian fluids. Applications of non-Newtonian fluids may be found in many industrial and engineering areas, glass fiber, hot rolling, casting, and paper production. Amongst many non-Newtonian fluids, Williamson fluid possesses shear thinning property (i.e., viscosity tends to decrease when shear stress increases). Abundant articles may be found in the literature that highlight the importance of Williamson fluid in numerous scenarios. Ramzan et al. [37] examined Williamson nanofluid flow over a Riga plate. They found that with an increase in Williamson fluid parameter velocity, distribution decreases. Ramzan et al. [38] also analyzed the numerical solution of the 2D MHD stagnation point of Williamson fluid flow under the effect of homogeneous–heterogeneous reactions over a linearly stretched surface and found opposite behavior of temperature and velocity distribution against the Williamson fluid parameter. Nadeem et al. [39,40] analyzed the 2D Williamson fluid flow over a stretching sheet considering the influence of nanosized particles, also characterized as Williamson nanofluid. They studied the Williamson nanofluid peristaltic flow in a curved channel, including compliant walls.
Homotopy analysis method (HAM) was suggested by Liao [41] in 1992 to solve highly nonlinear differential equations. This technique has an edge over the rest of the contemporary techniques. HAM is one of the best and simplest technique for obtaining the convergent series solution for weakly, as well as, highly nonlinear differential equations. This technique includes the concept of homotopy from topology. HAM is used for finding a convergent series solution with high nonlinearity. Homotopy discriminates itself from other methods in the following ways:
  • Freedom to choose large or small parameters;
  • Guaranteed series solution convergence; and
  • Freedom to choose linear operators and base function.
There have been many attempts in the literature to discuss the varied fluid problems utilizing the homotopy analysis method [42,43,44,45,46,47,48,49].
A literature review discloses that copious literature may be quoted in the case of 2D non-Newtonian flows. Less work is available on 3D geometries, and this group becomes narrower if we talk about 3D Williamson nanofluid flows. The subject matter of 3D MHD flow of Williamson nanofluid over a bidirectional stretched surface with second-order slip and double stratification is even more rarely discussed. The structure of this paper is as follows: In Section 2, we present the mathematical model. In Section 3, Section 4 and Section 5, we discuss the homotopic scheme in detail with zeroth and mth order solutions. In Section 6, we address the convergence analysis. In Section 7, we present the results and discuss their physical importance, and finally, we provide concluding remarks.

2. Mathematical Modeling

Here, we consider the steady 3D Williamson nano liquid flow with velocities of U w = a x in the x-direction and V w = b y in the y-direction, respectively, over a bidirectional extended sheet. While a and b are constants (Figure 1), concentration buoyancy force and thermal are used by the fluid with double stratification phenomena to study heat and mass transfers.
The following are the governing boundary layer equations:
u x + v y + w z = 0
u u x + v u y + w u z = ν u z z + 2 Γ ν u z u z z σ B 0 2 ρ u + [ g [ α 1 ( T T ) + α 2 ( T T ) 2 ] + g [ α 3 ( C C ) + α 4 ( C C ) 2 ] ]
u v x + v v y + w v z = ν v z z + 2 Γ ν v z v z z σ B 0 2 ρ v
q + λ E ( q t + V . q q . V + ( . V ) q ) = ( k T )
J + λ C ( J t + V . J J . V + ( . V ) J ) = D B C
q + λ E ( V . q q . V ) = ( k T )
J + λ C ( V . J J . V ) = D B C
under the supervision of above-mentioned consideration and the impression of thermophoresis and Brownian-motion, Equations (6) and (7) takes the form:
u T x + v T y + w T z + λ E ϕ E = 1 ρ c p z ( k T z ) + τ [ D B C z T z + D T T ( T z ) 2
u C x + v C y + w C z + λ C ϕ C = D B C z z + D T T T z z
where
ϕ E = u 2 T x x + v 2 T y y + w T z z + 2 u v T x y + 2 u w T x z + 2 v w T y z + [ ( u u x + v u y + w u z ) T x + ( u v x + v v y + w v z ) T y + ( u w x + v w y + w w z ) T z ]
ϕ C = u 2 C x x + v 2 C y y + w C z z + 2 u v C x y + 2 u w C x z + 2 v w C y z + [ ( u u x + v u y + w u z ) C x + ( u v x + v v y + w v z ) C y + ( u w x + v w y + w w z ) C z ]
following boundary conditions supports the above-mentioned system of equations:
u = U w + U s l i p      v = V w + V s l i p          w = 0 T = T w = T 0 + d 1 x         C = C w = C 0 + d 2 x                               at       z = 0 u 0     v 0      T T = T 0 + e 1 x      C C = C 0 + e 2 x            as       z
here
U s l i p = 2 3 ( 3 α l 3 α 3 2 1 l 2 K n ) Λ u z 1 4 ( l 4 + 2 K n 2 ( 1 l 2 ) ) Λ 2 u z z = A u z + B u z z
V s l i p = 2 3 ( 3 α l 3 α 3 2 1 l 2 K n ) Λ v z 1 4 ( l 4 + 2 K n 2 ( 1 l 2 ) ) Λ 2 v z z = C v z + D v z z ,
where A = 2 3 ( 3 α l 3 α 3 2 1 l 2 K n ) Λ , B = 1 4 ( l 4 + 2 K n 2 ( 1 l 2 ) ) Λ 2 , C = 2 3 ( 3 α l 3 α 3 2 1 l 2 K n ) Λ , D = 1 4 ( l 4 + 2 K n 2 ( 1 l 2 ) ) Λ 2 , l = min [ 1 K n , 1 ] , α describes momentum accommodation coefficient and varies from 0 α 1 K n denotes Knudsen number and Λ denotes molecular mean free path. On the basis of definition of l , we found for any particular estimates of K n we own 1 l 0 . The molecular mean free path is always positive. Therefore, we know that B, D ˂ 0 and C and A are positive.
To solve Equations (1),(3) and (8),(9), following similarity transformations are introduced:
u = a x f ( η )             v = a y g ( η )                w = a ν ( f ( η ) + g ( η ) ) θ ( η ) = T T T w T           ϕ ( η ) = C C C w C               η = a ν z
here, f , g , θ and ϕ are the non-dimensional form for both velocities, temperature and the concentration.
Condition for incompressibility is self-satisfied and Equations (2),(3) and (8),(9) reduce to:
f f 2 + ( f + g ) f + W e f f + λ ( 1 + β 2 θ ) θ + λ N r ( 1 + β 3 ϕ ) ϕ H a f = 0
g g 2 + ( f + g ) g + W e g g H a g = 0
( 1 + ε θ ) θ + ε θ 2 + Pr N b θ ϕ + Pr N t θ 2 Pr f ( S 1 + θ ) + Pr ( f + g ) θ δ t Pr ( ( f + g ) 2 θ 2 f θ ( f + g ) + ( f 2 f ( f + g ) ) ( S 1 + θ ) + ( f + g ) ( f + g ) θ ) = 0
ϕ + N t N b θ Pr L e f ( S 2 + ϕ ) + Pr L e ( f + g ) ϕ Pr L e δ c ( ( f + g ) 2 ϕ 2 f ( f + g ) ϕ + ( f 2 f ( f + g ) ) ( S 2 + ϕ ) + ( f + g ) ( f + g ) ϕ )
and boundary conditions hold the form
f ( 0 ) = 0    f ( 0 ) = 1 + γ 1 f ( 0 ) + γ 2 f ( 0 )     g ( 0 ) = 0 g ( 0 ) = 1 + γ 3 g ( 0 ) + γ 4 g ( 0 )     θ ( 0 ) = 1 S 1     ϕ ( 0 ) = 1 S 2 f ( ) 0     g ( ) 0    θ ( ) = 0    ϕ ( ) = 0                                          a s       z         
where the parameters given above are defined as follows:
λ = G r x Re x 2       G r x = g β T ( T w T ) x 3 ν 2    Re x = u w x ν N r = α 3 ( C w C 0 ) α 1 ( T w T 0 )    S 1 = e 1 d 1 γ 1 = A a ν      γ 2 = B a ν     γ 3 = C a ν      γ 4 = D a ν      β 2 = α 2 α 1 ( T w T 0 )    β = b a β 3 = α 4 α 3 ( C w C 0 )          α 1 = β T     W e = U w Γ 2 c ν      Pr = μ c p k    H a = σ B 0 2 ρ a S 2 = e 2 d 2 L e = α D B     N t = τ D T ( T w T 0 ) T ν    N b = τ D B ( C w C 0 ) T ν    δ c = λ C a     δ t = λ E a
C f x is the coefficients of Skin friction in x- and C f y in the y-direction are represented as follows:
C f x = τ w x ρ U w 2         C f y = τ w y ρ U w 2       
where   τ w x | z = 0 = u z + Γ 2 ( u z ) 2      and   τ w y | z = 0 = v z + Γ 2 ( v z ) 2
Coefficients of Skin friction in dimensionless forms are:
C f x Re 1 2 = [ f + W e 2 ( f ) 2 ] η = 0 C f y Re 1 2 = [ g + W e 2 ( g ) 2 ] η = 0

3. Homotopic Solutions

For the considered problem, ( L f , L g , L θ , L ϕ ) are the linear operators and ( f 0 , g 0 , θ 0 , ϕ 0 ) are the initial guesses expressed in the following form:
f 0 ( η ) = A 1 + γ 1 γ 2 ( 1 exp ( η ) )        g 0 ( η ) = β 1 + γ 3 γ 4 ( 1 exp ( η ) ) where        1 + γ 1 γ 2 0    1 + γ 3 γ 4 0 θ 0 ( η ) = ( 1 S 1 ) ( 1 exp ( η ) )                ϕ 0 ( η ) = ( 1 S 2 ) ( 1 exp ( η ) )
L f ( f ) = d 3 f d η 3 d f d η        L g ( g ) = d 3 g d η 3 d g d η L θ ( θ ) = d 2 θ d η 2 θ                L ϕ ( ϕ ) = d 2 ϕ d η 2 ϕ
these operators satisfy the following condition:
L f [ C 1 + C 2 exp ( η ) + C 3 exp ( η ) ] = 0       L g [ C 4 + C 5 exp ( η ) + C 6 exp ( η ) ] = 0 L θ [ C 9 exp ( η ) + C 10 exp ( η ) ] = 0

4. Zeroth Order Deformation

The zeroth order deformation problem is defined as follows:
( 1 p ) L f [ f ( η ; p ) f 0 ( η ) ] = p f N f [ f ( η ; p ) , g ( η ; p ) ]
( 1 p ) L g [ g ( η ; p ) g 0 ( η ) ] = p g N g [ f ( η ; p ) , g ( η ; p ) ]
( 1 p ) L θ [ θ ( η ; p ) θ 0 ( η ) ] = p θ N θ [ f ( η ; p ) , g ( η ; p ) , θ ( η ; p ) , ϕ ( η ; p ) ]
( 1 p ) L ϕ [ ϕ ( η ; p ) ϕ 0 ( η ) ] = p ϕ N ϕ [ f ( η ; p ) , g ( η ; p ) , θ ( η ; p ) , ϕ ( η ; p ) ]
f ( 0 ; p ) = 0       f ( 0 ; p ) = 1 + γ 1 f ( 0 ; p ) + γ 2 f ( 0 ; p )      f ( ; p ) = 0 g ( 0 ; p ) = 0       g ( 0 ; p ) = β + γ 3 g ( 0 ; p ) + γ 2 g ( 0 ; p )         g ( ; p ) = 0 θ ( 0 ; p ) = 1 S 1       θ ( ; p ) = 0        ϕ ( 0 ; p ) = 1 S 2              ϕ ( ; p ) = 0
N f [ f ( η ; p ) , g ( η ; p ) ] = 3 f ( η ; p ) η 3 ( f ( η ; p ) η ) 2 + ( f + g ) 2 f ( η ; p ) η 2 + W e 2 f ( η ; p ) η 2 3 f ( η ; p ) η 3 + λ ( 1 + β 2 θ ) θ + λ N r ( 1 + β 3 ϕ ) ϕ M f ( η ; p ) η          N g [ f ( η ; p ) , g ( η ; p ) ] = 3 g ( η ; p ) η 3 ( g ( η ; p ) η ) 2 + ( f + g ) 2 g ( η ; p ) η 2 + W e 2 g ( η ; p ) η 2 3 g ( η ; p ) η 3 M g ( η ; p ) η
N θ [ f ( η ; p ) , g ( η ; p ) , θ ( η ; p ) , ϕ ( η ; p ) ] = [ ( 1 + ε θ ) 2 θ ( η ; p ) η 2 + ε ( θ ( η ; p ) η ) 2 + Pr N b θ ( η ; p ) η ϕ ( η ; p ) η + Pr N t ( θ ( η ; p ) η ) 2 + Pr ( f + g ) θ ( η ; p ) η δ t Pr ( ( f + g ) 2 2 θ ( η ; p ) η 2 2 f ( η ; p ) η θ ( η ; p ) η ( f + g ) + ( ( c ) 2 2 f ( η ; p ) η 2 ( f + g ) ) ( S 1 + θ ) + ( f + g ) ( f ( η ; p ) η g ( η ; p ) η ) θ ( η ; p ) η )
N ϕ [ f ( η ; p ) , g ( η ; p ) , θ ( η ; p ) , ϕ ( η ; p ) ] = 2 ϕ ( η ; p ) η 2 + N t N b 2 θ ( η ; p ) η 2 PrLe f ( η ; p ) η ( S 2 + θ ) + PrLe ( f + g ) ϕ ( η ; p ) η PrLe δ c ( ( f + g ) 2 2 ϕ ( η ; p ) η 2 2 f ( η ; p ) η ( f + g ) ϕ ( η ; p ) η + ( ( f ( η ; p ) η ) 2 2 f ( η ; p ) η 2 ( f + g ) ) ( S 2 + ϕ ) + ( f + g ) ( f ( η ; p ) η g ( η ; p ) η ) ϕ ( η ; p ) η )
Here p [ 0 , 1 ] is embedding parameter and f , g , θ   and   ϕ are the non-zero auxiliary parameters.

5. m th -Order Deformation Problems

Here, we have
L f [ f m ( η ) χ m f m 1 ( η ) ] = f R f n ( η )
L g [ g m ( η ) χ m g m 1 ( η ) ] = g R g n ( η )
L θ [ θ m ( η ) χ m θ m 1 ( η ) ] = θ R θ n ( η )
L ϕ [ ϕ m ( η ) χ m ϕ m 1 ( η ) ] = ϕ R ϕ n ( η )
R f m ( η ) = f m 1 f m 1 2 + k = 0 m 1 ( f m 1 k + g m 1 k ) f k + W e k = 0 m 1 f m 1 k f k + λ θ m 1 + β 2 λ k = 0 m 1 θ m 1 k θ k + λ N r ϕ m 1 + λ N r β 3 k = 0 m 1 ϕ m 1 k ϕ k H a f m 1
R g m ( η ) = g m 1 g m 1 2 + k = 0 m 1 ( f m 1 k + g m 1 k ) g k + W e k = 0 m 1 g m 1 k g k H a g m 1
R θ m ( η ) = θ m 1 + ε k = 0 m 1 θ m 1 k θ k + ε k = 0 m 1 θ m 1 k θ k + Pr N b k = 0 m 1 θ m 1 k ϕ k + Pr N t k = 0 m 1 θ m 1 k θ k Pr ( S 1 + θ ) f m 1 + Pr k = 0 m 1 ( f m 1 k θ k +    g m 1 k θ k ) δ t Pr k = 0 m 1 f m 1 k l = 0 k f k 1 θ 1 δ t Pr k = 0 m 1 g m 1 k l = 0 k g k 1 θ 1 2 δ t Pr k = 0 m 1 f m 1 k l = 0 k g k 1 θ 1 2 δ t Pr k = 0 m 1 f m 1 k l = 0 k f k 1 θ 1 2 δ t Pr k = 0 m 1 g m 1 k l = 0 k f k 1 θ 1 S 1 δ t Pr k = 0 m 1 f m k 1 f k 1 + S 1 δ t Pr k = 0 m 1 f m 1 k f k + S 1 δ t Pr k = 0 m 1 g m 1 k f k δ t Pr k = 0 m 1 f m k 1 l = 0 k f k 1 θ 1 + δ t Pr k = 0 m 1 f m 1 k l = 0 k f k 1 θ 1 + δ t Pr k = 0 m 1 g m 1 k l = 0 k f k 1 θ 1 + δ t Pr k = 0 m 1 f m 1 k l = 0 k f k 1 θ 1 + δ t Pr k = 0 m 1 f m 1 k l = 0 k g k 1 θ 1 + δ t Pr k = 0 m 1 g m 1 k l = 0 k f k 1 θ 1 + δ t Pr k = 0 m 1 g m 1 k l = 0 k g k 1 θ 1
R ϕ m ( η ) = ϕ m 1 + N t N b k = 0 m 1 θ m 1 k θ k PrLeS 2 f m 1 PrLe k = 0 m 1 f m 1 k θ k f m 1 + PrLe k = 0 m 1 f m 1 k ϕ k + PrLe k = 0 m 1 g m 1 k ϕ k PrLe δ c k = 0 m 1 f m 1 k l = 0 k f k 1 ϕ 1 PrLe δ c k = 0 m 1 g m 1 k l = 0 k g k 1 ϕ 1 2 δ c PrLe k = 0 m 1 f m 1 k l = 0 k g k 1 ϕ 1 2 δ c PrLe k = 0 m 1 f m 1 k l = 0 k f k 1 ϕ 1 2 δ c PrLe k = 0 m 1 g m 1 k l = 0 k f k 1 ϕ 1 + S 2 δ c PrLe k = 0 m 1 f m k 1 f k 1 + S 2 δ c PrLe k = 0 m 1 f m 1 k f k + S 2 δ c PrLe k = 0 m 1 g m 1 k f k PrLe δ c k = 0 m 1 f m k 1 l = 0 k f k 1 ϕ 1 + δ c PrLe k = 0 m 1 f m 1 k l = 0 k f k 1 ϕ 1 + δ c PrLe k = 0 m 1 g m 1 k l = 0 k f k 1 ϕ 1 + δ c PrLe k = 0 m 1 f m 1 k l = 0 k f k 1 ϕ 1 + δ c PrLe k = 0 m 1 f m 1 k l = 0 k g k 1 ϕ 1 + δ c PrLe k = 0 m 1 g m 1 k l = 0 k f k 1 ϕ 1 + δ c PrLe k = 0 m 1 g m 1 k l = 0 k g k 1 ϕ 1
χ m = { 0 ,       m 1 1 ,        m > 1
the final solutions can be transcribed in the subsequent forms:
f m ( η ) = f m ( η ) + D 1 + D 2 e η + D 3 e η g m ( η ) = g m ( η ) + D 4 + D 5 e η + D 6 e η θ m ( η ) = θ m ( η ) + D 7 e η + D 8 e η ϕ m ( η ) = ϕ m ( η ) + D 9 e η + D 10 e η
where f m , g m , θ m ,   and   ϕ m symbolize the special solutions.

6. Convergence Analysis

HAM is used to obtain the solution of higher order nonlinear problems or those in series form. It gives several choices to control and modify the convergence region for the series solutions. Figure 2 represents the curves behavior of all distributions. Characteristic parameters f , g , θ and ϕ have permissible ranges 1.6 f 0.4 , 2.15 g 0.2 , 2.75 θ 0.8 and 2.6 ϕ 0.6 when γ = 0.2, Le = 1, N t = 0.2, N b = 0.3, Pr = 1.0, ε = 0.3, λ = 0.002, β1= β3 = 0.2, β = 0.1 and M = 0.2 Table 1 represents the numerical results obtained for series solutions depicting the convergence of approximations up to the 25th order of approximations, that is, enough for series solution convergence. It can be verified that the graphical depiction in Figure 2 and the tabular results in Table 1 are in total consensus.
Table 2 was developed to validate the results obtained in the current model for skin friction in both directions by comparison with Malik et al. [50], and excellent harmony in both outcomes is achieved.

7. Results and Discussion

In this section, we analyze the impact of appearing factors on particular distributions in Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11, Figure 12, Figure 13, Figure 14, Figure 15, Figure 16, Figure 17, Figure 18, Figure 19, Figure 20 and Figure 21. The behavior on velocity profiles of β (ratio parameter) is described in Figure 3 and Figure 4. It is noticed that contradictory behavior shown by both velocities ( f , g ) for an increasing rate of β . As β = b a , a was smaller for higher values of β , which specified a decreasing velocity rate along the x-direction, or b with higher values specified an increasing rate along the y-direction. In Figure 5 and Figure 6,     δ t and    δ c illustrate the influence of thermal relaxation and the concentration relaxation factor on temperature concentration and distributions. We found that both concentration and temperature fields associated with the thickness of the boundary layers were the functions of decreasing δ c and δ t , respectively. Furthermore, δ c = 0 and δ t = 0 existing model will transform into classical laws of Fick’s and Fourier’s respectively. The influence of thermal conductivity ε on the temperature distribution is described in Figure 7. For higher values of ε , an increasing rate for the thermal boundary layer is found, which in result increases the temperature distribution. In Figure 8 the impact of Lewis number Le on concentration field is described. The strength of Lewis number depends on smaller estimations of mass diffusivity than the thermal diffusivity, which shows that exhausted Brownian motion coefficient decreases nanoparticle concentration profile. Figure 9 highlighted the influence of mixed convective factor λ on the velocity field ( g ). Higher estimations of λ produce stronger buoyancy force, which indicates an increasing rate in the velocity field ( g ). The behavior of Prandtl number Pr against temperature distribution is presented in Figure 10. It is inspected that heat diffusion is very slow from the heated surface for higher estimates of Pr than smaller estimations of Pr. Therefore, temperature decreases with increasing values of Pr. Figure 11 and Figure 12 show the influence of Ha (Hartmann number) on both velocity profiles ( f , g ). Retardation in the fluid motion is seen due to resistance effered by strong Lorentz force. This act finally points out the decreasing rate on both velocity distributions. Figure 13 and Figure 14 demonstrate the impact of Brownian motion factor N b on concentration and temperature distribution. For the larger estimates of N b , fluid temperature increases and rapidly reduces the deposition of particles far away from the fluid on the stretched surface. Due to which it increases and decreases concentration. The influence of the N t on the concentration distribution is described in Figure 15. When the estimates of N t are high and they are directly proportional to the temperature distribution, as a result it enhances the concentration distribution and its concentration associated with thickness of boundary layer. Figure 16 illustrates the influence of N t on the temperature profile. When the values of N t are higher, nanoparticles move from ambient fluid with higher temperature to the ambient fluid with lower temperature, and as a results temperature is higher in the boundary layer region. Finally, we identified the thickness of augmented thermal boundary layer. Figure 17 and Figure 18 are drawn to display the influence of Williamson fluid parameter We on both velocity profiles. Increasing values of We decrease both velocities profiles. By increasing Williamson factor, relaxation time enhances. It causes to increase liquid viscosity, which results in decrease in the velocity profile. In Figure 19 the impact of stretching ratio factor β and mixed convective parameter λ on coefficient of Skin friction along x -direction is displayed. Which highlighted that the coefficient of Skin friction shpws increasing behavior versus β and λ . Similar tendency can be observed in Hartmann number H a and mixed convective parameter λ, against coefficient of Skin friction in x -direction, as displayed in Figure 20. Analysis of the impact of λ and Nr on Skin friction is described in Figure 21. It is noted that a thinner boundary layer is associated with larger λ , which result in higher velocity gradient near the wall. That’s why Skin friction reduces against λ .

8. Concluding Remarks

Three-dimensional Williamson nanofluid flow was investigated considering the Cattaneo–Christov heat flux model. The originality of the envisioned mathematical model was boosted by considering the influence of double stratification and second-order slip. HAM was applied to obtain the problem solution in series form. The salient outcomes of the problem are as follows:
  • The stretching ratio parameter had an opposite impact on both velocities.
  • Increasing values of concentrations and temperature distributions decreased the thermal and concentration relaxation parameters, respectively.
  • The higher temperature was in direct proportion with the thermal conductivity parameter.
  • Velocity increased for values of the mixed convection parameter.
  • For the large values of the Prandtl number, fluid temperature decreased.
  • Both velocity components were decreasing functions of the Hartmann number.
  • Skin friction coefficients against the x- and y- directions displayed an accelerating tendency with respect to values of the stretching ratio and mixed convection parameters.

Author Contributions

Conceptualization, M.R. and A.L.; methodology, M.R.; software, S.K.; validation, S.Y., S.N., D.L. and Y.N; funding acquisition, Y.N.

Funding

This research was supported by Basic Science Research Program through the National Research Foundation of Korea(NRF) funded by the Ministry of Education(NRF-2017R1D1A3B03028309) and also supported by the Soonchunhyang University Research Fund.

Conflicts of Interest

The authors declare no conflicts of interest regarding this publication.

Nomenclature

k0Elastic parameter
RexLocal Reynold parameter
S1, S2Thermal and concentration stratification parameter
kThermal conductivity
CfySkin friction coefficients in the y-direction
HaHartmann number
LeLewis number
NbBrownian motion parameter
NtThermophoresis parameter
WeWilliamsons fluid parameter
NrRatio of concentration to buoyancy forces
GrxGrashof number
CfxCoefficients of skin friction in the x-direction
u, v, wVelocity components
α3,α4Linear and nonlinear coefficients of concentration expansions
JMass flux
TAmbient temperature
DBBrownian diffusion coefficient
θTemperature parameter
A,B,C,DConstants
UwVelocity along x-axis
σElectrical conductivity
k0Elastic parameter
α1Normal stress moduli
γ1, γ3First-order slip parameter
γ2, γ4Second-order slip parameter
νKinematic viscosity
βStretching ratio parameter
δt, δcThermal and concentration relaxation parameters
λMixed convection parameter
ρDensity of fluid
λcRelaxation time of mass flux
λERelaxation time of heat flux
β23Nonlinear temperature’s and concentration’s convection parameter λ
α1, α2Linear and nonlinear coefficients of thermal expansions
qNormal heat flux
CAmbient concentration
DTThermophoretic diffusion coefficient
f, gNondimensional velocity parameters
φConcentration parameter
ρDensity of the fluid

References

  1. Straughan, B. Thermal convection with the Cattaneo–Christov model. Int. J. Heat Mass Transf. 2010, 53, 95–98. [Google Scholar] [CrossRef]
  2. Khan, W.A.; Khan, M.; Alshomrani, A.S.; Ahmad, L. Numerical investigation of generalized Fourier’s and Fick’s laws for Sisko fluid flow. J. Mol. Liq. 2016, 224, 1016–1021. [Google Scholar] [CrossRef]
  3. Hayat, T.; Khan, M.I.; Farooq, M.; Alsaedi, A.; Waqas, M.; Yasmeen, T. Impact of Cattaneo–Christov heat flux model inflow of variable thermal conductivity fluid over a variable thicked surface. Int. J. Heat Mass Transf. 2016, 99, 702–710. [Google Scholar] [CrossRef]
  4. Waqas, M.; Hayat, T.; Farooq, M.; Shehzad, S.A.; Alsaedi, A. Cattaneo-Christov heat flux model for the flow of variable thermal conductivity generalized Burgers fluid. Int. J. Heat Mass Transf. 2016, 220, 642–648. [Google Scholar] [CrossRef]
  5. Khan, A.U.; Ahmed, N.; Mohyud-Din, S.T. Thermo-diffusion and diffusion-thermo effects on the flow of second-grade fluid between two inclined plane walls. J. Mol. Liq. 2016, 224, 1074–1082. [Google Scholar]
  6. Ghadikolaei, S.S.; Hosseinzadeh, K.; Yassari, M.; Sadeghi, H.; Ganji, D.D. Analytical and numerical solution of non-Newtonian second-grade fluid flow on a stretching sheet. Therm. Sci. Eng. Prog. 2018, 5, 309–316. [Google Scholar] [CrossRef]
  7. Khan, I.; Malik, M.Y.; Salahuddin, T.; Khan, M.; Rehman, K.U. Homogenous–heterogeneous reactions in MHD flow of Powell–Eyring fluid over a stretching sheet with Newtonian heating. Neural Comput. Appl. 2018, 30, 3581–3588. [Google Scholar] [CrossRef] [Green Version]
  8. Ibrahim, W.; Shankar, B.; Nandeppanavar, M.M. MHD stagnation point flow and heat transfer due to nanofluid towards a stretching sheet. Int. J. Heat Mass Transf. 2013, 56, 1–9. [Google Scholar] [CrossRef]
  9. Chamkha, A.J.; Al-Mudhaf, A. Unsteady heat and mass transfer from a rotating vertical cone with a magnetic field and heat generation or absorption effects. Int. J. Therm. Sci. 2005, 44, 267–276. [Google Scholar] [CrossRef]
  10. Pullepu, B.; Chamkha, A.J.; Pop, I. Unsteady laminar free convection flows past a non-isothermal vertical cone in the presence of a magnetic field. Chem. Eng. Commun. 2012, 199, 354–367. [Google Scholar] [CrossRef]
  11. Akbar, N.S.; Nadeem, S.; Haq, R.U.; Khan, Z.H. Numerical solutions of Magnetohydrodynamic boundary layer flow of tangent hyperbolic fluid towards a stretching sheet. Indian J. Phys. 2013, 87, 1121–1124. [Google Scholar] [CrossRef]
  12. Seini, I.Y.; Makinde, O.D. Boundary layer flow near stagnation-points on a vertical surface with a slip in the presence of the transverse magnetic field. Int. J. Numer. Methods Heat Fluid Flow 2014, 24, 643–653. [Google Scholar] [CrossRef]
  13. Ravindran, R.; Ganapathirao, M.; Pop, I. Effects of chemical reaction and heat generation/absorption on unsteady mixed convection MHD flow over a vertical cone with non-uniform slot mass transfer. Int. J. Heat Mass Transf. 2014, 73, 743–751. [Google Scholar] [CrossRef]
  14. Bovand, M.; Rashidi, S.; Esfahani, J.A.; Saha, S.C.; Gu, Y.T.; Dehesht, M. Control of flow around a circular cylinder wrapped with a porous layer by magnetohydrodynamic. J. Magn. Magn. Mater. 2016, 401, 1078–1087. [Google Scholar] [CrossRef]
  15. Ellahi, R.; Bhatti, M.M.; Pop, I. Effects of the hall and ion slip on MHD peristaltic flow of Jeffrey fluid in a non-uniform rectangular duct. Int. J. Numer. Methods Heat Fluid Flow 2016, 26, 1802–1820. [Google Scholar] [CrossRef]
  16. Mishra, S.R.; Pattnaik, P.K.; Bhatti, M.M.; Abbas, T. Analysis of heat and mass transfer with MHD and chemical reaction effects on viscoelastic fluid over a stretching sheet. Indian J. Phys. 2017, 91, 1219–1227. [Google Scholar] [CrossRef]
  17. Hussain, A.; Malik, M.Y.; Awais, M.; Salahuddin, T.; Bilal, S. Computational and physical aspects of MHD Prandtl-Eyring fluid flow analysis over a stretching sheet. Neural Comput. Appl. 2019, 31, 425–433. [Google Scholar] [CrossRef]
  18. Khan, W.A.; Pop, I. Boundary-layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 2010, 53, 2477–2483. [Google Scholar] [CrossRef]
  19. Makinde, O.D.; Aziz, A. Boundary layer flow of a nanofluid past a stretching sheet with a convective boundary condition. Int. J. Therm. Sci. 2011, 50, 1326–1332. [Google Scholar] [CrossRef]
  20. Nadeem, S.; Mehmood, R.; Akbar, N.S. Non-orthogonal stagnation point flow of a nano non-Newtonian fluid towards a stretching surface with heat transfer. Int. J. Heat Mass Transf. 2013, 57, 679–689. [Google Scholar] [CrossRef]
  21. Hatami, M.; Jing, D.; Song, D.; Sheikholeslami, M.; Ganji, D.D. Heat transfer and flow analysis of nanofluid flow between parallel plates in the presence of a variable magnetic field using HPM. J. Magn. Magn. Mater. 2015, 396, 275–282. [Google Scholar] [CrossRef]
  22. Hayat, T.; Imtiaz, M.; Alsaedi, A. MHD 3D flow of a nanofluid in the presence of convective conditions. J. Mol. Liq. 2015, 212, 203–208. [Google Scholar] [CrossRef]
  23. Sheikholeslami, M.; Ganji, D.D.; Javed, M.Y.; Ellahi, R. Effect of thermal radiation on magnetohydrodynamics nanofluid flow and heat transfer by means of a two-phase model. J. Magn. Magn. Mater. 2015, 374, 36–43. [Google Scholar] [CrossRef]
  24. Sheikholeslami, M.; Rokni, H.B. Nanofluid two-phase model analysis in the existence of induced magnetic field. Int. J. Heat Mass Trans. 2017, 107, 288–299. [Google Scholar] [CrossRef]
  25. Hassan, M.; Marin, M.; Alsharif, A.; Ellahi, R. Convective heat transfer flow of a nanofluid in a porous medium over wavy surface. Phys. Lett. A 2018, 382, 2749–2753. [Google Scholar] [CrossRef]
  26. Nayak, M.K.; Shaw, S.; Pandey, V.S.; Chamkha, A.J. Combined effects of slip and convective boundary condition on MHD 3D stretched flow of nanofluid through porous media inspired by non-linear thermal radiation. Indian J. Phys. 2018, 92, 1–12. [Google Scholar] [CrossRef]
  27. Hosseini, S.R.; Sheikholeslami, M.; Ghasemian, M.; Ganji, D.D. Nanofluid heat transfer analysis in a microchannel heat sink (MCHS) under the effect of a magnetic field by means of the KKL model. Powder Technol. 2018, 324, 36–47. [Google Scholar] [CrossRef]
  28. Lu, D.; Ramzan, M.; Mohammad, M.; Howari, F.; Chung, J.D. A Thin Film Flow of Nanofluid Comprising Carbon Nanotubes Influenced by Cattaneo-Christov Heat Flux and Entropy Generation. Coatings 2019, 9, 296. [Google Scholar] [CrossRef] [Green Version]
  29. Li, Z.; Shafee, A.; Ramzan, M.; Rokni, H.B.; Al-Mdallal, Q.M. Simulation of natural convection of Fe 3 O 4-water ferrofluid in a circular porous cavity in the presence of a magnetic field. Eur. Phys. J. Plus 2019, 134, 77. [Google Scholar] [CrossRef]
  30. Suleman, M.; Ramzan, M.; Ahmad, S.; Lu, D.; Muhammad, T.; Chung, J.D. A Numerical Simulation of Silver–Water Nanofluid Flow with Impacts of Newtonian Heating and Homogeneous–Heterogeneous Reactions Past a Nonlinear Stretched Cylinder. Symmetry 2019, 11, 295. [Google Scholar] [CrossRef] [Green Version]
  31. Lu, D.; Li, Z.; Ramzan, M.; Shafee, A.; Chung, J.D. Unsteady squeezing carbon nanotubes based nano-liquid flow with Cattaneo–Christov heat flux and homogeneous–heterogeneous reactions. Appl. Nanosci. 2019, 9, 169–178. [Google Scholar] [CrossRef]
  32. Ramzan, M.; Sheikholeslami, M.; Saeed, M.; Chung, J.D. On the convective heat and zero nanoparticle mass flux conditions in the flow of 3D MHD Couple Stress nanofluid over an exponentially stretched surface. Sci. Rep. 2019, 9, 562. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Li, Z.; Sheikholeslami, M.; Shafee, A.; Ramzan, M.; Kandasamy, R.; Al-Mdallal, Q.M. Influence of adding nanoparticles on solidification in a heat storage system considering radiation effect. J. Mol. Liq. 2019, 273, 589–605. [Google Scholar] [CrossRef]
  34. Farooq, U.; Lu, D.C.; Munir, S.; Suleman, M.; Ramzan, M. Flow of Rheological Nanofluid over a Static Wedge. J. Nanofluids 2019, 8, 1362–1366. [Google Scholar] [CrossRef]
  35. Qayyum, S.; Khan, M.I.; Hayat, T.; Alsaedi, A. Comparative investigation of five nanoparticles in flow of viscous fluid with Joule heating and slip due to rotating disk. Phys. B Condens. Matter 2018, 534, 173–183. [Google Scholar] [CrossRef]
  36. Nguyen, T.; van der Meer, D.; van den Berg, A.; Eijkel, J.C. Investigation of the effects of time periodic pressure and potential gradients on viscoelastic fluid flow in circular narrow confinements. Microfluid. Nanofluid. 2017, 21, 37. [Google Scholar] [CrossRef] [Green Version]
  37. Ramzan, M.; Bilal, M.; Chung, J.D. Radiative Williamson nanofluid flow over a convectively heated Riga plate with chemical reaction-A numerical approach. Chin. J. Phys. 2017, 55, 1663–1673. [Google Scholar] [CrossRef]
  38. Ramzan, M.; Bilal, M.; Chung, J.D. MHD stagnation point Cattaneo–Christov heat flux in Williamson fluid flow with homogeneous–heterogeneous reactions and convective boundary condition—A numerical approach. J. Mol. Liq. 2017, 225, 856–862. [Google Scholar] [CrossRef]
  39. Nadeem, S.; Hussain, S.T. Flow and heat transfer analysis of Williamson nanofluid. Appl. Nanosci. 2014, 4, 1005–1012. [Google Scholar] [CrossRef] [Green Version]
  40. Nadeem, S.; Maraj, E.N.; Akbar, N.S. Investigation of peristaltic flow of Williamson nanofluid in a curved channel with compliant walls. Appl. Nanosci. 2014, 4, 511–521. [Google Scholar] [CrossRef] [Green Version]
  41. Liao, S.J. Beyond Perturbation; Chapman & Hall/CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  42. Jafarimoghaddam, A. On the Homotopy Analysis Method (HAM) and Homotopy Perturbation Method (HPM) for a nonlinearly stretching sheet flow of Eyring-Powell fluids. Eng. Sci. Technol. Int. J. 2019, 22, 439–451. [Google Scholar] [CrossRef]
  43. Freidoonimehr, N.; Rostami, B.; Rashidi, M.M. Predictor homotopy analysis method for nanofluid flow through expanding or contracting gaps with permeable walls. Int. J. Biomath. 2015, 8, 1550050. [Google Scholar] [CrossRef]
  44. Ray, A.K.; Vasu, B.; Bég, O.A.; Gorla, R.S.; Murthy, P.V.S.N. Homotopy semi-numerical modeling of non-Newtonian nanofluid transport external to multiple geometries using a revised Buongiorno Model. Inventions 2019, 4, 54. [Google Scholar] [CrossRef] [Green Version]
  45. Shukla, N.; Rana, P.; Bég, O.A.; Singh, B.; Kadir, A. Homotopy study of magnetohydrodynamic mixed convection nanofluid multiple slip flow and heat transfer from a vertical cylinder with entropy generation. Propuls. Power Res. 2019, 8, 147–162. [Google Scholar] [CrossRef]
  46. Shqair, M. Solution of different geometries reflected reactors neutron diffusion equation using the homotopy perturbation method. Results Phys. 2019, 12, 61–66. [Google Scholar] [CrossRef]
  47. Ramzan, M.; Farooq, M.; Alsaedi, A.; Hayat, T. MHD three-dimensional flow of couple stress fluid with Newtonian heating. Eur. Phys. J. Plus 2013, 128, 49. [Google Scholar] [CrossRef]
  48. Hussain, T.; Shehzad, S.A.; Hayat, T.; Alsaedi, A.; Al-Solamy, F.; Ramzan, M. Radiative hydromagnetic flow of Jeffrey nanofluid by an exponentially stretching sheet. PLoS ONE 2014, 9, e103719. [Google Scholar] [CrossRef]
  49. Ramzan, M.; Yousaf, F. Boundary layer flow of three-dimensional viscoelastic nanofluid past a bi-directional stretching sheet with Newtonian heating. AIP Adv. 2015, 5, 057132. [Google Scholar] [CrossRef]
  50. Malik, M.Y.; Bilal, S.; Salahuddin, T.; Rehman, K.U. Three-dimensional Williamson fluid flow over a linearly stretching surface. Math. Sci. Lett. 2017, 6, 53–61. [Google Scholar] [CrossRef]
Figure 1. Schematic flow diagram.
Figure 1. Schematic flow diagram.
Coatings 09 00849 g001
Figure 2. curves for f , g , θ    and ϕ .
Figure 2. curves for f , g , θ    and ϕ .
Coatings 09 00849 g002
Figure 3. Impact of β against f’.
Figure 3. Impact of β against f’.
Coatings 09 00849 g003
Figure 4. Impact of β against g’.
Figure 4. Impact of β against g’.
Coatings 09 00849 g004
Figure 5. Impact of δc against ϕ.
Figure 5. Impact of δc against ϕ.
Coatings 09 00849 g005
Figure 6. Impact of δt against θ.
Figure 6. Impact of δt against θ.
Coatings 09 00849 g006
Figure 7. Impact of ∈ against θ.
Figure 7. Impact of ∈ against θ.
Coatings 09 00849 g007
Figure 8. Impact of Le against ϕ.
Figure 8. Impact of Le against ϕ.
Coatings 09 00849 g008
Figure 9. Impact of λ against f’.
Figure 9. Impact of λ against f’.
Coatings 09 00849 g009
Figure 10. Impact of Pr against θ.
Figure 10. Impact of Pr against θ.
Coatings 09 00849 g010
Figure 11. Impact of H a against f’.
Figure 11. Impact of H a against f’.
Coatings 09 00849 g011
Figure 12. Impact of H a against g’.
Figure 12. Impact of H a against g’.
Coatings 09 00849 g012
Figure 13. Impact of Nb against ϕ.
Figure 13. Impact of Nb against ϕ.
Coatings 09 00849 g013
Figure 14. Impact of Nb against θ.
Figure 14. Impact of Nb against θ.
Coatings 09 00849 g014
Figure 15. Impact of Nt against ϕ.
Figure 15. Impact of Nt against ϕ.
Coatings 09 00849 g015
Figure 16. Impact of Nt against θ.
Figure 16. Impact of Nt against θ.
Coatings 09 00849 g016
Figure 17. Impact of We against f’.
Figure 17. Impact of We against f’.
Coatings 09 00849 g017
Figure 18. Impact of We against g’.
Figure 18. Impact of We against g’.
Coatings 09 00849 g018
Figure 19. Impact of λ and β against Skin friction.
Figure 19. Impact of λ and β against Skin friction.
Coatings 09 00849 g019
Figure 20. Impact of λ and H a against Skin friction.
Figure 20. Impact of λ and H a against Skin friction.
Coatings 09 00849 g020
Figure 21. Impact of Nr and λ against Skin friction.
Figure 21. Impact of Nr and λ against Skin friction.
Coatings 09 00849 g021
Table 1. Convergence analysis of series solution using different order approximations when M = 0.2, γ1 =0.2, γ2 = 0.2, N b = 0.3, N t = 0.2, λ = 0.002, Le = 1.0, Pr = 1.0 δt = 0.2, δc = 0.2, ε =0.1, β = 0.1.
Table 1. Convergence analysis of series solution using different order approximations when M = 0.2, γ1 =0.2, γ2 = 0.2, N b = 0.3, N t = 0.2, λ = 0.002, Le = 1.0, Pr = 1.0 δt = 0.2, δc = 0.2, ε =0.1, β = 0.1.
Order of Approximations f ( 0 ) g ( 0 ) θ ( 0 ) ϕ ( 0 )
11.135090.072260.311470.53000
51.270390.073290.356900.46424
101.335920.073440.360060.44888
201.373890.073600.360990.44350
251.387310.073680.361580.44215
301.387310.073680.361580.44215
Table 2. Comparative estimates of We with Malik et al. [50] for skin friction along both directions in the limiting case, i.e., by considering the Hartmann number, second-order slip, temperature, and concentration profiles to zero.
Table 2. Comparative estimates of We with Malik et al. [50] for skin friction along both directions in the limiting case, i.e., by considering the Hartmann number, second-order slip, temperature, and concentration profiles to zero.
We [ f + W e 2 ( f ) 2 ] η = 0 [ g + W e 2 ( g ) 2 ] η = 0
[40]Present Outcomes[50]Present Outcomes
0.11.09341.09330.46610.4660
0.21.26951.26950.48410.4841
0.31.33401.33410.50250.5024
0.41.49151.49160.52200.5221

Share and Cite

MDPI and ACS Style

Ramzan, M.; Liaquet, A.; Kadry, S.; Yu, S.; Nam, Y.; Lu, D. Impact of Second-Order Slip and Double Stratification Coatings on 3D MHD Williamson Nanofluid Flow with Cattaneo–Christov Heat Flux. Coatings 2019, 9, 849. https://doi.org/10.3390/coatings9120849

AMA Style

Ramzan M, Liaquet A, Kadry S, Yu S, Nam Y, Lu D. Impact of Second-Order Slip and Double Stratification Coatings on 3D MHD Williamson Nanofluid Flow with Cattaneo–Christov Heat Flux. Coatings. 2019; 9(12):849. https://doi.org/10.3390/coatings9120849

Chicago/Turabian Style

Ramzan, Muhammad, Asma Liaquet, Seifedine Kadry, Sungil Yu, Yunyoung Nam, and Dianchen Lu. 2019. "Impact of Second-Order Slip and Double Stratification Coatings on 3D MHD Williamson Nanofluid Flow with Cattaneo–Christov Heat Flux" Coatings 9, no. 12: 849. https://doi.org/10.3390/coatings9120849

APA Style

Ramzan, M., Liaquet, A., Kadry, S., Yu, S., Nam, Y., & Lu, D. (2019). Impact of Second-Order Slip and Double Stratification Coatings on 3D MHD Williamson Nanofluid Flow with Cattaneo–Christov Heat Flux. Coatings, 9(12), 849. https://doi.org/10.3390/coatings9120849

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