Next Article in Journal
Estimation and Analysis of Higher-Order Harmonics in Advanced Integrated Circuits to Implement Noise-Free Future-Generation Micro- and Nanoelectromechanical Systems
Next Article in Special Issue
Entropy Generation Analysis and Radiated Heat Transfer in MHD (Al2O3-Cu/Water) Hybrid Nanofluid Flow
Previous Article in Journal
Additive Manufacturing of Bulk Nanocrystalline FeNdB Based Permanent Magnets
Previous Article in Special Issue
Electroosmotic Mixing of Non-Newtonian Fluid in a Microchannel with Obstacles and Zeta Potential Heterogeneity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Analysis of Thermal Radiative Maxwell Nanofluid Flow Over-Stretching Porous Rotating Disk

by
Shuang-Shuang Zhou
1,
Muhammad Bilal
2,*,
Muhammad Altaf Khan
3 and
Taseer Muhammad
4,5
1
School of Science, Hunan City University, Yiyang 413000, China
2
Department of Mathematics, City University of Science and Information Technology, Peshawar 25000, Pakistan
3
Institute for Groundwater Studies, Faculty of Natural and Agricultural Sciences, University of Free State, Bloemfontein 9300, South Africa
4
Department of Mathematics, College of Sciences, King Khalid University, Abha 61413, Saudi Arabia
5
Mathematical Modelling and Applied Computation Research Group (MMAC), Department of Mathematics, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
*
Author to whom correspondence should be addressed.
Micromachines 2021, 12(5), 540; https://doi.org/10.3390/mi12050540
Submission received: 21 March 2021 / Revised: 18 April 2021 / Accepted: 3 May 2021 / Published: 10 May 2021
(This article belongs to the Special Issue Non-Newtonian Microfluidics)

Abstract

:
The fluid flow over a rotating disk is critically important due to its application in a broad spectrum of industries and engineering and scientific fields. In this article, the traditional swirling flow of Von Karman is optimized for Maxwell fluid over a porous spinning disc with a consistent suction/injection effect. Buongiorno’s model, which incorporates the effect of both thermophoresis and Brownian motion, describes the Maxwell nanofluid nature. The dimensionless system of ordinary differential equations (ODEs) has been diminished from the system of modeled equations through a proper transformation framework. Which is numerically computed with the bvp4c method and for validity purposes, the results are compared with the RK4 technique. The effect of mathematical abstractions on velocity, energy, concentration, and magnetic power is sketched and debated. It is perceived that the mass transmission significantly rises with the thermophoresis parameter, while the velocities in angular and radial directions are reducing with enlarging of the viscosity parameter. Further, the influences of thermal radiation Rd and Brownian motion parameters are particularly more valuable to enhance fluid temperature. The fluid velocity is reduced by the action of suction effects. The suction effect grips the fluid particles towards the pores of the disk, which causes the momentum boundary layer reduction.

1. Introduction

The researchers have been interested in Maxwell nanofluid flow over a porous spinning disc because of its many uses in engineering and innovation. Non-Newtonian fluids are important in a variety of manufactured liquids, including plastics, polymers, pulps, toothpaste and fossil fluids. To simulate the analysis of these liquids, a variety of models have been suggested. Shear stress and shear rate are linked in non-Newtonian liquids because of their nonlinear existence. The momentum equation in these fluids involves dynamic nonlinear terms, making it difficult to solve. A variety of mathematical models exist in the literature to simulate the performance of these fluids.
In the present era, the role of nanotechnology to fulfill the increasing demand for energy and face energy challenges is remarkable. The usage of nanoparticles in ordinary base fluid (water, kerosene oil, etc.) effectively enhances the heat transfer and improves their thermal properties. The applications of nanoparticles in certain fields of engineering and industry are in the cooling systems of electronic devices and in cancer therapy, heat exchangers, transformer cooling, nuclear reactors, and space cooling systems. Due to the ability of oil wetting and dispersing, they are also used for cleaning purposes, in power generation, microfabrication, hyperthermia, and metallurgical purposes.

2. Literature Review

The rotating disk phenomena are widely used in centrifugal filtration, turbomachines, the braking system of vehicles, jet motors, sewing machines, turbine systems, heat exchangers and computer disk drives, etc. Von Karman [1] for the first time introduced similarity transformation. To solve Navier-Stokes equations, he studied fluid flow over an infinite rotating disk. Cochran [2] employed Von Karman’s similarity transformation to incompressible fluid over a rotating disk and examined the asymptotic solution. Wagner [3] investigated the mechanism of heat transfer over the rotating disk, by considering Von Karman’s velocity distribution, and analyzed convection in the non-turbulent flow. Turkyimazogl [4] looked at fluid movement over a spinning disc that was stretching under the influence of a static electric field. Liang et al. [5] reported a comparative study between semi-analytical model and experimental data to yield the best settlement. Millsaps and Pohlhausen [6] used Von Karman’s similarity approach and analyzed heat distribution with the consequences of entropy generation over revolving disk. The three-dimensional (3D) magnetohydrodynamics (MHD) stagnation flow of ferrofluid, the numerical solution was revealed by Mustafa et al. [7]. Mustafa et al. [8], by taking MHD nanofluid over rotating surface with the effects of partial slips, observed that the boundary layer thickness and momentum transport are reduced due to slip effects. Rashidi et al. [9] have used a spinning disc to perform a viscous dissipation review for MHD nanofluid.
The people of the modern world are facing many challenges due to the increasing demand for energy by the latest technologies. Firstly Choi [10] presented the nanofluids terminology. The Brownian motion and thermophoresis mechanisms bring about a significant role in improving the thermal properties of base fluid presented by Buongiorno [11]. Turkyilmazolglu [12] analytically studied the energy and momentum equations of nanofluid flow, to deduce heat and flow transport. Pourmehran [13] considered Cu and Al2O3 nanoparticles to study heat and flow transfer in the microchannel. Hatami et al. [14] reported the heat transfer in nanofluid with the phenomena of natural convection. The Oldroyd-B fluid with nanoparticles over stretching sheet surface was reported by Nadeem et al. [15]. Aziz and Afify [16] have used the technique of the Lie group, to study non-Newtonian nanofluids. Yang et al. [17] studied the convective heat with Buongiorno Model’s for nanofluid in the concentric annulus.
The study of a Newtonian fluid, due to its wide applicability in different fields of science and engineering, attracted the attention of scientists and researchers during the last few decades. Its major role is in geophysics, polymer solution, paper production, cosmetic processes, exotic lubricants, paints, suspensions, colloidal solutions, nuclear and chemical industries, pharmaceuticals, oil reservoirs, bioengineering, etc. [18]. Xiao et al. [19] attempted a fractal model for the capillary flow through a torturous capillary with the non-smooth surface in porous media. Attia [20] evaluated the Reiner-Rivlin numerical simulations for thermal convection over a porous spinning disc qualitatively. Griffiths [21] tested the Newtonian fluid Carreau viscosity model and high shear stresses on spinning discs. The numerical analysis of Reiner-Rivlin fluid flow for heat transfer and slip flow over a spinning disc is treated by Mustafa and Tabassum [22]. The micropolar fluid for thermophoretic diffusion generated by the rotation of the disk was examined by Doh and Muthtamilselvan [23].
Darcy’s law is a mathematical equation that explains how fluid flows through a porous medium. Henry Darcy developed the law based on the effects of studies on the flow of water across sand beds, laying the groundwork for hydrogeology, a branch of earth sciences [24]. Fourier’s law in heat conduction, Ohm’s law in electrical networks, and Fick’s law in diffusion theory are all examples of this law. Morris Muskat [25] improved Darcy’s equation for a single-phase flow by incorporating viscosity into Darcy’s single (fluid) phase equation. It is easy to see that viscous fluids have a harder time passing through a porous medium than less viscous fluids. Rasool et al. [26,27,28] numerically simulated the Darcy-Forchheimer effect on MHD nanoliquid flow between stretching non-linear sheets. Rasool et al. [29] scrutinized the consequences of thermal radiation, chemical reaction and Dufour-Soret on incompressible steady Darcy-Forchheimer flow of nanoliquid. Shafiq et al. [30] studied nanofluid flow under the influence of convective boundary conditions and thermal slip over a spinning frame. They found that the axial and transverse velocity fields all drop significantly due to the Forchheimer number’s strong retardation. Skin friction is intensified by the Forchheimer number and porosity ratios, while skin friction is diminished by all slip parameters.
Viscoelastic fluids are a subclass of Newtonian fluid having memory effects. The intensity of energy discharged by these fluids is mainly accountable for recovery after the stress is removed. The Maxwell flow regime is the most basic viscoelastic fluid model, expressing memory effects by fluid relaxation time [31]. The attitude of the current model is very close to that of other geomaterials and polymers models. The aim of the present work is to provide a mathematical model for unsteady boundary layer flow of non-Newtonian Maxwell nanofluid with the heat transmission over a porous spinning disc. The present work has many industrial and engineering applications, which increases its worth. Using a resemblance method, the system of ODEs is limited to a structure of PDEs. A boundary value solver (bvp4c) technique is used to draw a numerical solution to the problem while RK4 method has been applied for validity.

3. Formulation of the Problem

Consider an unsteady hybrid nanoliquid flow over a stretching porous spinning disc. The magnetic force B 0 is introduced to the disc vertically. The disc rotates and stretches at different speeds ( u , v ) = ( c r , c Ω ) , where c and w are the spinning and extending rates, respectively. The disk temperature is represented by τ w . The formulation of the problems is conducted in ( r , φ , z ) cylindrical coordinates, where u , v , w is velocity component increasing in ( r , φ , z ) direction. At z the axis, the motion of the disk is assumed to be axisymmetric. The thermal radiation is significant in modeling the energy equation. The viscosity of a fluid is taken to be temperature-dependent μ ( τ ) = μ 0 e ζ ( τ τ 0 ) . The concentration and temperature are represented by ( C w , τ w ) and ( C , τ ) represent the concentration and temperature above the disk surface.

3.1. Governing Equations

Under the presuppositions stated above, the flow equations are as observes [31]:
. V = 0 ,
ρ f ( V . ) V = p + . S + j × B ,
( V . ) τ = α 2 τ + τ * ( D B C . τ + D τ τ τ . τ ) + . q r a d ,
( V . ) C = D B 2 C + D τ τ 2 τ ,
. B = 0 ,
ρ B t = ρ + ( V × B ) + ρ σ μ 2 2 V .
where D B , D T ,   ρ f , V and α are the coefficients of Brownian motion, thermophoretic diffusion, fluid density, velocity, and thermal and diffusivity respectively. Equations (1)–(6) are simplified because of the boundary layer approximation concept:
u r + u r + w z = 0 ,
u t + u u r + w u z v 2 r = [ 1 ρ ( μ ( τ ) u z ) λ 1 ( u 2 2 u r 2 + w 2 u z 2 + 2 u w 2 u r z 2 u v r v r 2 v w r v z + u v 2 r 2 + v 2 r 2 u r ) σ B 0 2 ρ ( u + w λ 1 u z ) ] ,
v t + u v r + w v z u v r = [ 1 ρ ( μ ( τ ) v z ) λ 1 ( u 2 2 v r 2 + w 2 2 v z 2 + 2 u w 2 v r z + 2 u v r u r + 2 v w r u z + 2 u 2 v r 2 v 2 r 2 v r ) σ B 0 2 ρ ( v + w λ 1 v z ) ] ,
τ t + u τ r + w τ z = k ρ c p ( 2 τ z 2 ) + τ * ( D B τ z C z + D τ τ ( τ z ) 2 ) 1 ρ c p q r z ,
C t + u C r + w C z = D B ( 2 C z 2 ) + D τ τ ( 2 τ z 2 ) ,
B r t = [ w B r z B r w z + u B z z + B z u z + 1 σ μ 2 ( 2 B r r 2 + 2 B r z 2 + 1 r B r r B r r 2 ) ] ,
B z t = [ w B r r + B r w r + 1 r w B r u B z r B z u r 1 r u B z + 1 σ μ 2 ( 2 B z r 2 + 2 B z z 2 + 1 r B z r ) ] .
The temperature difference within a flow is assumed to be small, therefore higher-order in Taylor series and ignored at τ , by using Rosseland approximation, the simple form of radiation heat flux is as follow [31]:
q r = 4 σ * τ 4 3 k * z = 16 σ * τ 3 3 k * τ z ,
using Equation (14) in (10), we get
τ t + u τ r + w τ z = k ρ c p ( 2 τ z 2 ) + τ * ( D B τ z C z + D τ τ ( τ z ) 2 ) 16 σ * τ 3 3 k * τ z .
The boundary conditions are:
u = c r , v = Ω r ,   w = W ,   C = C w ,   T = T w ,   B r = 1 , B z = 1         a t z = 0 u 0 ,   v 0 ,   w 0 ,   C C ,   T T ,   B r 0 , B z 0     a s         Z .

3.2. Similarity Transformation

Considering the following transformation, to reduce the system of PDEs to the system of ODEs:
u = c r 1 α t F ( η ) , ν = Ω r 1 α t G ( η ) ,   w = c v 1 α t H ( η ) , η = c v z ,   B r = c r M 0 1 α t M ( η ) , B z = M 0 ( 2 ν f c ) 1 2 1 α t N ( η ) ,     T = ( T ) + Θ ( η ) ( T W T ) C = ( C ) + ϕ ( η ) ( C W C ) . }
The following system of ODEs is obtained by using Equation (17) in Equations (7)–(13) and (15) and (16):
F = F Θ + e δ Θ δ { S ( F η 2 + F ) + F 2 + H F G 2 } + β 1 e δ Θ δ { H F + 2 F F H 2 H G G } M e δ Θ δ ( F + β 1 H F ) ,
G = δ G Θ + e δ Θ { S ( G η 2 + G ) + 2 F G + H G } + β 1 e δ Θ { 2 F H G + 2 F H G } + M e δ Θ ( G + β 1 H G ) 1 + β 1 H 2 e δ Θ ,
Θ = 4 R d Θ 2 Θ 2 ( Θ 1 ) 3 6 Θ 2 Θ ( Q w 1 ) 2 3 Θ 2 ( Q w 1 ) Pr ( S 2 ( Θ η ) H Θ + N b Θ ϕ + N t Θ 3 ) ( 1 + 4 3 ) R d + 4 3 R d Θ 3 ( Q w 1 ) 3 + 3 Θ 2 ( Q w 1 ) 2 + 3 Θ ( Θ 1 ) ,
ϕ = S c ( A ϕ η + H ϕ ) N t N b Θ ,
M = B t [ H M + M H + F N + N F + S ( M η 2 + M ) ] ,
M = B t [ 2 H M + 2 N F S 2 ( N η + N ) ] .
the transforms conditions are:
F ( 0 ) = 1 ,   G ( 0 ) = ω ,   H = W s ,   Θ ( 0 ) = 1 ,   f ( 0 ) = 1 ,   M ( 0 ) = 0 ,   N ( 0 ) = 1 , F ( ) = 0 ,   G ( ) = 0 ,   H ( 0 ) = 0 ,   Θ ( ) = 0 ,   f ( ) = 0 ,   M ( ) = 0 ,   N ( ) = 0 .
where P r = ν / α is the Prandtl number, S c = v / D B is the Schmidt number, β 1 = λ 1 c the Deborah number, while suction/injection parameter, Brownian motion, variable viscosity, thermophoresis parameter, magnetic parameter, temperature ratio, and thermal radiation are defined as [31]:
W s = w c ν , N b = τ * D B ( C w - C ) v , δ = ζ ( τ w - τ ) , N t = τ * D τ ( τ w - τ ) τ v , M = σ B 0 2 c p , R d = 4 σ * T 3 k k * .
The skin friction C f x , local Sherwood S h r and Nusselt number N u r are mathematically can be written as [32,33]:
C f = τ z r 2 + τ z φ 2 ρ ( Ω r ) 2 ,
S h r = r ( C w C ) ( τ C ) | z = 0 ,
N u r = r ( τ w τ ) [ 1 + 16 σ * τ 3 3 k k * ] ( τ z ) | z = 0 .
The skin friction, Sherwood and Nusselt numbers have a non-dimensional structure as:
C f Re r 1 / 2 = f 2 ( 0 ) + g 2 ( 0 ) ,
R e 1 2 S h r = ϕ ( 0 ) ,
R e 1 2 N u r = ( 1 + 4 3 R d { 1 + ( Θ w 1 ) Θ ( 0 ) } 3 ) Θ ( 0 ) .
Here R e = r u ν is the local Reynold number.

4. Solution Procedures

The higher-order model equation is brought down to first order by choosing variables:
χ 1 = H ,   χ 2 = F ,   χ 3 = F ,   χ 4 = G ,   χ 5 = G ,   χ 6 = θ ,   χ 7 = θ ,   χ 8 = ϕ ,   χ 9 = ϕ ,   χ 10 = M ,   χ 11 = M ,   χ 13 = N ,   χ 14 = N . }
χ 1 = 2 χ 2 , χ 1 = χ 3 ,
χ 3 = χ 7 χ 3 + e δ θ δ { S ( η χ 3 2 + χ 2 ) + χ 2 2 + χ 1 χ 2 2 χ 4 2 } + β 1 e δ θ δ { χ 1 χ 3 + 2 χ 1 χ 2 χ 3 2 χ 1 χ 4 χ 5 } M 0 e δ θ δ ( χ 2 + β 1 χ 1 χ 3 ) , χ 4 = χ 5 ,
χ 5 = δ χ 5 χ 7 + e δ θ { S ( χ 5 η 2 ) + 2 χ 2 χ 4 + χ 1 χ 5 } + β 1 { 2 χ 1 χ 2 χ 5 2 χ 1 χ 4 χ 2 } + M 0 e δ θ ( χ 4 + β 1 χ 1 χ 5 ) 1 + β 1 e δ θ χ 1 2 ,
χ 6 = χ 7 , χ 7 = 4 3 R d χ 7 2 ( θ w 1 ) { 3 χ 6 2 ( θ w 1 ) 2 + 6 χ 6 ( θ w 1 ) 3 } Pr { N b χ 7 χ 9 + N t χ 7 2 A χ 7 χ 10 χ 1 χ 7 } 1 4 3 R d 4 3 R d χ 6 ( θ w 1 ) { ( θ w 1 ) χ 6 2 3 χ 6 ( θ w 1 ) + 3 } , χ 8 = χ 9 ,
χ 9 = N t N b ( 4 3 R d χ 7 2 ( θ w 1 ) { 3 χ 6 2 ( θ w 1 ) 2 + 6 χ 6 ( θ w 1 ) 3 } Pr { N b χ 7 χ 9 + N t χ 7 2 A χ 7 χ 10 χ 1 χ 7 } 1 4 3 R d 4 3 R d χ 6 ( θ w 1 ) { ( θ w 1 ) χ 6 2 3 χ 6 ( θ w 1 ) + 3 } ) S c ( A χ 9 χ 10 + χ 91 χ 9 ) ,
χ 10 = χ 11 , χ 11 = B t [ 2 χ 11 χ 2 χ 1 χ 11 + χ 2 χ 14 + χ 13 χ 2 + S ( η χ 12 2 + χ 11 ) ] , χ 12 = χ 13 , χ 13 = B t [ 2 χ 11 χ 1 + 2 χ 13 χ 2 + S 2 ( η χ 14 + χ 13 ) ]
The boundary conditions are:
χ 1 ( 0 ) = 1 , χ 2 ( 0 ) = 0 , χ 4 ( 0 ) = 1 , χ 6 ( 0 ) = 1 , χ 8 ( 0 ) = 1 , χ 11 ( 0 ) = 0 , χ 13 ( 0 ) = 1 , χ 2 ( ) = 0 , χ 4 ( ) = 0 , χ 6 ( ) = 0 , χ 8 ( ) = 0 , χ 11 ( ) = 1 , χ 13 ( ) = 0 .

5. Results and Discussion

The discussion section is devoted to understanding better the graphical and physical description. The system of non-linear Equations (18)–(23) along with their boundary conditions, Equation (24), are solved through numerical method bvp4c. The configuration of the problem is described in Figure 1. The velocities, energy profile, concentration distribution ϕ ( η ) , magnetic strength in the radial direction M ( η ) , and azimuthal magnetic strength N ( η ) are explored graphically through different physical constraints Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9 and Figure 10. While keeping P r = 6.7 , ω = 1.0 , δ = 0.5 , β = 0.1 , M = 1.1 , N b = 0.5 , N t = 0.7 , Q w = 1.2 , S c = 2.0 , and R d = 0.5 .
Figure 2a–c depicts the behavior velocities profiles against the variation of Deborah number β . All three velocity shows decreasing behavior for incremented of β . Deborah number is the measure of the content evaluation period to content recreational time, so having optimum stress relaxation or eliminating observation time increases the value of β . It reflects the fluid’s solid-like reaction. The hydrodynamic boundary layer thins out, and the velocity experiences more resistance.
The implications of the suction factor ( W s < 0 ) on the velocities profile and heat spectrum on the disc surface are depicted in Figure 3a–c. It seems to be that as the suction velocity rises, the velocity decreases. Because the suction velocity draws the fluid particles towards the pores in the disk, which causes the momentum boundary layer reduction. The enhancement of suction velocity also decreases the fluid temperature.
Figure 4a–d are drawn to depict the effects of injection velocity parameter ( W s > 0 ) on the radial, azimuthal, tangential velocity profile and temperature distribution. It can be seen that the enhancing of injection velocity increases the velocity and fluid temperature.
The radial, azimuthal, tangential velocity profile and temperature distribution variation are illustrated in Figure 5a–d. Figure 5a depicts the dominant behavior of radial velocity with ω . Physically, fluid particles are moved in a radial direction owing to centrifugal force with the enhancement of parameters ω . The increasing value of ω shows that the rotation parameter becomes greater than extending. The radial velocity exceeds the disc stretching velocity too close to the disc stretching surface. It does, however, gradually fade away from the disc. It is concluded that the impact of centrifugal force is limited and dominant in the vicinity of the disk’s surface. Besides, angular velocity G ( η ) near the disk flourished shown in Figure 5b. The axial velocity increasing against the strengthening of ω is illustrated in Figure 5c. Figure 5d demonstrates the decreases of temperature θ ( η ) with ω . Physically, a faster spinning disc reduces the width of the thermal boundary, which play a significant role in the cooling of the system.
Figure 6a–d are rough sketches that show the effects of the viscosity factor on the velocity and temperature spectrum on the disk’s surface. The viscosity parameter causes radial velocity, tangential velocity, and temperature profile to increase, while azimuthal velocity reduces.
Figure 7a,b indicates how the heat transfer changes when the Prandtl number Pr and the thermal radiation factor change. The thermal dispersion of a strong Prandtl fluid is low, while the thermal diffusivity of a low Prandtl fluid is high. Figure 7b is drawn to explore the influence of radiation parameters R d on the thermal mechanism of fluid. With a higher value of the radiation parameter, an increasing trend in heat is analyzed. Because, with enhancing of thermal radiation the fluid absorbs more heat, and as a result increment occurs in boundary layer thickness and fluid temperature. Figure 7c is sketched to illustrate temperature θ ( η ) variation versus temperature ratio Q w . While Figure 7d indicates the consequence of molecular diffusion on the thermal performance of a Maxwell nanofluid. The Brownian motion produces random movement between fluid particles, which generate more heat, and as a result, the fluid temperature increases.
Figure 8a,b accordingly represents the action of the concentration field as a function of the thermophoresis term and Brownian. It can be shown that with Nt, the nanofluid concentration field raises while with Nb, the concentration profile drops. The suction parameter ( W s ) < 0 and Schmidt number Sc effect are illustrated in Figure 8c and Figure 8d respectively. The concentration field reduces with an increase of both suction velocity and Schmidt number. The consequences of different parameters ( B t , W s ) on the axial and radial magnetic strength profiles are illustrated through Figure 9a–d respectively. It can observe that both parameters Batchlor number Bt and injection parameter Ws positively effect the magnetic strength profile along axial and radial direction.
Figure 10a,b express the nature of heat θ ( 0 ) and mass transfer ϕ ( 0 ) against radiation parameter Rd and Brownian motion parameter Nb, respectively. Because of the improving effect of radiation, the fluid temperature also rises, which enhances the heat transmission rate. Table 1 displays the numerical outcomes for skin friction and compared it with the published literature. Table 2 illustrates the comparison of bvp4c and RK4 techniques for the numerical outcomes. The numerical outputs for Sherwood number S h r and Nusselt number N u r are plotted in Table 3.

6. Conclusions

In the present mathematical model, the Maxwell nanoliquid flow over a porous spinning disk with suction/injection effects has been examined. The flow is studied in the context of magnetization and radiation. The numerical results are found through bvp4c, while for comparison purposes, the computation is carried out via the RK4 technique. From the above studies the following conclusions have been drawn:
  • The centrifugal force is more effective and dominant near the disk surface, so causes the maximum velocity in the neighborhood of the surface of the disk.
  • The fluid velocity reduces by the action of suction effects. Because the suction effect gripping the fluid particles towards the pores of the disk, which causes the momentum boundary layer reduction.
  • The axial, radial and tangential velocities are boosts with the increases of injection parameter.
  • The enhancement in rotation parameter also boosts the azimuthal and radial flows.
  • The thermal energy profile enhances by the consequence of Brownian motion parameter Nb. The Brownian motion produces random movement between fluid particles, which generate more heat, as a result, the fluid temperature increases.
  • The fluid temperature decline with Prandtl number P r , while incline with thermal radiation parameter, R d .
  • The nanofluid concentration field enhances with N t , while N b causes the reduction of concentration profile.

Author Contributions

M.B. wrote the original manuscript and performed the numerical simulations. M.A.K. reviewed the mathematical results and S.-S.Z., T.M. and M.A.K. restructured the manuscript and revised the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Deanship of Scientific Research at King Khalid University, Abha, Saudi Arabia for funding this work through research groups program under grant number R.G.P-1/142/42 and the APC was funded by S.S.Z.

Acknowledgments

The authors extend their appreciation to the Deanship of Scientific Research at King Khalid University, Abha, Saudi Arabia for funding this work through research groups program under grant number R.G.P-1/142/42.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kármán, T.V. Über laminare und turbulente Reibung. ZAMM J. Appl. Math. Mech. Z. Für Angew. Math. Mech. 1921, 1, 233–252. [Google Scholar] [CrossRef] [Green Version]
  2. Cochran, W. The flow due to a rotating disc. Math. Proc. Camb. Philos. Soc. 1934, 30, 365–375. [Google Scholar] [CrossRef]
  3. Wagner, C. Heat transfer from a rotating disk to ambient air. J. Appl. Phys. 1948, 19, 837–839. [Google Scholar] [CrossRef]
  4. Turkyilmazoglu, M. MHD fluid flow and heat transfer due to a stretching rotating disk. Int. J. Therm. Sci. 2012, 51, 195–201. [Google Scholar] [CrossRef]
  5. Liang, M.; Fu, C.; Xiao, B.; Luo, L.; Wang, Z. A fractal study for the effective electrolyte diffusion through charged porous media. Int. J. Heat Mass Transf. 2019, 137, 365–371. [Google Scholar] [CrossRef]
  6. Millsaps, K.; Pohlhausen, K. Heat transfer by laminar flow from a rotating plate. J. Aeronaut. Sci. 1952, 19, 120–126. [Google Scholar] [CrossRef]
  7. Mustafa, I.; Javed, T.; Ghaffari, A. Heat transfer in MHD stagnation point flow of a ferrofluid over a stretchable rotating disk. J. Mol. Liq. 2016, 219, 526–532. [Google Scholar] [CrossRef]
  8. Mustafa, M. MHD nanofluid flow over a rotating disk with partial slip effects: Buongiorno model. Int. J. Heat Mass Transf. 2017, 108, 1910–1916. [Google Scholar] [CrossRef]
  9. Rashidi, M.; Abelman, S.; Mehr, N.F. Entropy generation in steady MHD flow due to a rotating porous disk in a nanofluid. Int. J. Heat Mass Transf. 2013, 62, 515–525. [Google Scholar] [CrossRef]
  10. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Argonne National Lab.: Lemont, IL, USA, 1995. [Google Scholar]
  11. Buongiorno, J. Convective transport in nanofluids. J. Heat Transf. 2006, 128, 240–250. [Google Scholar] [CrossRef]
  12. Turkyilmazoglu, M. Analytical solutions of single and multi-phase models for the condensation of nanofluid film flow and heat transfer. Eur. J. Mech. B Fluids 2015, 53, 272–277. [Google Scholar] [CrossRef]
  13. Pourmehran, O.; Rahimi-Gorji, M.; Hatami, M.; Sahebi, S.; Domairry, G. Numerical optimization of microchannel heat sink (MCHS) performance cooled by KKL based nanofluids in saturated porous medium. J. Taiwan Inst. Chem. Eng. 2015, 55, 49–68. [Google Scholar] [CrossRef]
  14. Hatami, M.; Song, D.; Jing, D. Optimization of a circular-wavy cavity filled by nanofluid under the natural convection heat transfer condition. Int. J. Heat Mass Transf. 2016, 98, 758–767. [Google Scholar] [CrossRef]
  15. Nadeem, S.; Haq, R.U.; Akbar, N.S.; Lee, C.; Khan, Z.H. Numerical study of boundary layer flow and heat transfer of Oldroyd-B nanofluid towards a stretching sheet. PLoS ONE 2013, 8, e69811. [Google Scholar] [CrossRef]
  16. Afify, A.A.; Abd El-Aziz, M. Lie group analysis of flow and heat transfer of non-Newtonian nanofluid over a stretching surface with convective boundary condition. Pramana 2017, 88, 31. [Google Scholar] [CrossRef]
  17. Yang, C.; Li, W.; Nakayama, A. Convective heat transfer of nanofluids in a concentric annulus. Int. J. Therm. Sci. 2013, 71, 249–257. [Google Scholar] [CrossRef]
  18. Hayat, T.; Qayyum, S.; Alsaedi, A.; Asghar, S. Radiation effects on the mixed convection flow induced by an inclined stretching cylinder with non-uniform heat source/sink. PLoS ONE 2017, 12, e0175584. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Xiao, B.; Huang, Q.; Chen, H.; Chen, X.; Long, G. A fractal model for capillary flow through a single tortuous capillary with roughened surfaces in fibrous porous media. Fractals 2021, 29, 2150017. [Google Scholar]
  20. Attia, H.A. Numerical study of flow and heat transfer of a non-Newtonian fluid on a rotating porous disk. Appl. Math. Comput. 2005, 163, 327–342. [Google Scholar] [CrossRef]
  21. Griffiths, P.T. Flow of a generalised Newtonian fluid due to a rotating disk. J. Non-Newton. Fluid Mech. 2015, 221, 9–17. [Google Scholar] [CrossRef] [Green Version]
  22. Tabassum, M.; Mustafa, M. A numerical treatment for partial slip flow and heat transfer of non-Newtonian Reiner-Rivlin fluid due to rotating disk. Int. J. Heat Mass Transf. 2018, 123, 979–987. [Google Scholar] [CrossRef]
  23. Doh, D.; Muthtamilselvan, M. Thermophoretic particle deposition on magnetohydrodynamic flow of micropolar fluid due to a rotating disk. Int. J. Mech. Sci. 2017, 130, 350–359. [Google Scholar] [CrossRef]
  24. Whitaker, S. Flow in porous media I: A theoretical derivation of Darcy’s law. Transp. Porous Media 1986, 1, 3–25. [Google Scholar] [CrossRef]
  25. Muskat, M. The flow of homogeneous fluids through porous media. Soil Sci. 1938, 46, 169. [Google Scholar] [CrossRef]
  26. Rasool, G.; Shafiq, A.; Alqarni, M.S.; Wakif, A.; Khan, I.; Bhutta, M.S. Numerical Scrutinization of Darcy-Forchheimer Relation in Convective Magnetohydrodynamic Nanofluid Flow Bounded by Nonlinear Stretching Surface in the Perspective of Heat and Mass Transfer. Micromachines 2021, 12, 374. [Google Scholar] [CrossRef]
  27. Rasool, G.; Shafiq, A. Numerical exploration of the features of thermally enhanced chemically reactive radiative Powell–Eyring nanofluid flow via Darcy medium over non-linearly stretching surface affected by a transverse magnetic field and convective boundary conditions. Appl. Nanosci. 2020, 1–18. [Google Scholar] [CrossRef]
  28. Rasool, G.; Zhang, T. Darcy-Forchheimer nanofluidic flow manifested with Cattaneo-Christov theory of heat and mass flux over non-linearly stretching surface. PLoS ONE 2019, 14, e0221302. [Google Scholar] [CrossRef]
  29. Rasool, G.; Shafiq, A.; Baleanu, D. Consequences of Soret–Dufour effects, thermal radiation, and binary chemical reaction on Darcy Forchheimer flow of nanofluids. Symmetry 2020, 12, 1421. [Google Scholar] [CrossRef]
  30. Shafiq, A.; Rasool, G.; Khalique, C.M. Significance of thermal slip and convective boundary conditions in three-dimensional rotating Darcy-Forchheimer nanofluid flow. Symmetry 2020, 12, 741. [Google Scholar] [CrossRef]
  31. Ahmed, J.; Khan, M.; Ahmad, L. Stagnation point flow of Maxwell nanofluid over a permeable rotating disk with heat source/sink. J. Mol. Liq. 2019, 287, 110853. [Google Scholar] [CrossRef]
  32. Ferdows, M.; Shamshuddin, M.; Zaimi, K. Dissipative-Radiative Micropolar Fluid Transport in a NonDarcy Porous Medium with Cross-Diffusion Effects. CFD Lett. 2020, 12, 70–89. [Google Scholar] [CrossRef]
  33. Devi, S.S.U.; Mabood, F. Entropy anatomization on Marangoni Maxwell fluid over a rotating disk with nonlinear radiative flux and Arrhenius activation energy. Int. Commun. Heat Mass Transf. 2020, 118, 104857. [Google Scholar] [CrossRef]
Figure 1. Stretchable porous rotating disk.
Figure 1. Stretchable porous rotating disk.
Micromachines 12 00540 g001
Figure 2. The out-turn of Deborah number β on the axial, radial and azimuthal velocity profiles, respectively.
Figure 2. The out-turn of Deborah number β on the axial, radial and azimuthal velocity profiles, respectively.
Micromachines 12 00540 g002
Figure 3. The out-turn of the suction parameter ( W s < 0 ) on the axial velocity, radial velocity and temperature distribution profiles, respectively.
Figure 3. The out-turn of the suction parameter ( W s < 0 ) on the axial velocity, radial velocity and temperature distribution profiles, respectively.
Micromachines 12 00540 g003
Figure 4. The out-turn of injection parameter ( W s > 0 ) on the axial, radial and azimuthal velocity profiles and temperature distribution, respectively.
Figure 4. The out-turn of injection parameter ( W s > 0 ) on the axial, radial and azimuthal velocity profiles and temperature distribution, respectively.
Micromachines 12 00540 g004
Figure 5. The out-turn of rotation parameter ω on the axial, radial and azimuthal velocity and temperature profiles, respectively.
Figure 5. The out-turn of rotation parameter ω on the axial, radial and azimuthal velocity and temperature profiles, respectively.
Micromachines 12 00540 g005
Figure 6. The out-turn of viscosity parameter δ on the axial, radial and azimuthal velocity and temperature profiles, respectively.
Figure 6. The out-turn of viscosity parameter δ on the axial, radial and azimuthal velocity and temperature profiles, respectively.
Micromachines 12 00540 g006
Figure 7. The out-turn of different parameters ( P r , R d , Θ , N b ) on the temperature distribution profile.
Figure 7. The out-turn of different parameters ( P r , R d , Θ , N b ) on the temperature distribution profile.
Micromachines 12 00540 g007
Figure 8. The out-turn of different parameters ( N t , N b , W s , S c ) on the concentration profile.
Figure 8. The out-turn of different parameters ( N t , N b , W s , S c ) on the concentration profile.
Micromachines 12 00540 g008
Figure 9. The out-turn of different parameters ( B t , W s ) on the axial and radial magnetic strength profiles, respectively.
Figure 9. The out-turn of different parameters ( B t , W s ) on the axial and radial magnetic strength profiles, respectively.
Micromachines 12 00540 g009
Figure 10. The out-turn of radiation parameter and Brownian coefficient versus the heat and mass transmission profiles, respectively.
Figure 10. The out-turn of radiation parameter and Brownian coefficient versus the heat and mass transmission profiles, respectively.
Micromachines 12 00540 g010
Table 1. The numerical outcomes for skin friction F ( 0 ) .
Table 1. The numerical outcomes for skin friction F ( 0 ) .
ω Mustafa et al. [7]Ahmed et al. [31]Present Paper
0−1.1737−1.1379−1.1380
1−0.9483−0.9485−0.9487
2−0.3262−0.3264−0.3266
53.19373.119373.11939
1012.720912.720912.7811
2040.905740.905740.9058
Table 2. Comparison between Runge Kutta order four and bvp4c method.
Table 2. Comparison between Runge Kutta order four and bvp4c method.
η RK4bvp4cAbsolute Error
1.01.0000001.0000008.146310 × 10 13
1.21.1998311.1998313.387821 × 10 9
1.40.9884590.9884592.845561 × 10 9
1.60.8791890.8791892.813281 × 10 9
1.80.5393930.5393933.287961 × 10 9
Table 3. The comparison of R K 4 and Bvp4c for Sherwood S h r and Nusselt number N u r , while keeping ω = 1.0 , δ = 0.5 , β = 0.1 , M = 1.1 , N b = 0.5 , N t = 0.7 , Q w = 1.2 and R d = 0.5 .
Table 3. The comparison of R K 4 and Bvp4c for Sherwood S h r and Nusselt number N u r , while keeping ω = 1.0 , δ = 0.5 , β = 0.1 , M = 1.1 , N b = 0.5 , N t = 0.7 , Q w = 1.2 and R d = 0.5 .
S h r N u r
Pr Bvp4c R K 4 S c Bvp4c R K 4
3.00.77182840.77182851.01.4579861.457995
4.00.69998850.69998851.51.5312111.531220
5.00.62900890.62900882.01.5969491.596949
6.00.56323720.56323702.51.6856391.685639
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhou, S.-S.; Bilal, M.; Khan, M.A.; Muhammad, T. Numerical Analysis of Thermal Radiative Maxwell Nanofluid Flow Over-Stretching Porous Rotating Disk. Micromachines 2021, 12, 540. https://doi.org/10.3390/mi12050540

AMA Style

Zhou S-S, Bilal M, Khan MA, Muhammad T. Numerical Analysis of Thermal Radiative Maxwell Nanofluid Flow Over-Stretching Porous Rotating Disk. Micromachines. 2021; 12(5):540. https://doi.org/10.3390/mi12050540

Chicago/Turabian Style

Zhou, Shuang-Shuang, Muhammad Bilal, Muhammad Altaf Khan, and Taseer Muhammad. 2021. "Numerical Analysis of Thermal Radiative Maxwell Nanofluid Flow Over-Stretching Porous Rotating Disk" Micromachines 12, no. 5: 540. https://doi.org/10.3390/mi12050540

APA Style

Zhou, S. -S., Bilal, M., Khan, M. A., & Muhammad, T. (2021). Numerical Analysis of Thermal Radiative Maxwell Nanofluid Flow Over-Stretching Porous Rotating Disk. Micromachines, 12(5), 540. https://doi.org/10.3390/mi12050540

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