Next Article in Journal
Drug Delivery Systems and Flavonoids: Current Knowledge in Melanoma Treatment and Future Perspectives
Next Article in Special Issue
Replenishment of the Gas in a Hydrophobically-Structured Surface by Mass Transfer at the Liquid-Gas Interface for Improving the Stability of Entrapped Gas
Previous Article in Journal
Double-Sided Sapphire Optrodes with Conductive Shielding Layers to Reduce Optogenetic Stimulation Artifacts
Previous Article in Special Issue
The Effect of Fin Shape on the Heat Transfer and the Solution Time of a Microchannel Evaporator in a CO2 Air Conditioning System—A Numerical Investigation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nonlinear Thermal Diffusion and Radiative Stagnation Point Flow of Nanofluid with Viscous Dissipation and Slip Constrains: Keller Box Framework Applications to Micromachines

1
Department of Industrial Engineering, Jazan University, Jazan 82822, Saudi Arabia
2
Department of Mathematics, University of Wah, Wah Cantt 47040, Pakistan
3
Department of Mathematics, COMSATS University Islamabad, Sahiwal 57000, Pakistan
4
Department of Mathematics and Statistics, Riphah International University I-14, Islamabad 44000, Pakistan
5
Department of Mechanical Engineering, Lebanese American University, Beirut 1102, Lebanon
6
Department of Industrial & Systems Engineering, College of Engineering, Princess Nourah bint Abdulrahman University, P.O. Box 84428, Riyadh 11671, Saudi Arabia
7
Center of Research, Faculty of Engineering, Future University in Egypt, New Cairo 11835, Egypt
8
Mechanical Engineering Department, College of Engineering and Islamic Architecture, Umm Al-Qura University, Makkah 21955, Saudi Arabia
9
Research Unity: Materials, Energy and Renewable Energies, Faculty of Science of Gafsa, University of Gafsa, Gafsa 2100, Tunisia
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(11), 1839; https://doi.org/10.3390/mi13111839
Submission received: 31 August 2022 / Revised: 25 October 2022 / Accepted: 25 October 2022 / Published: 27 October 2022
(This article belongs to the Special Issue Heat and Mass Transfer in Microchannels)

Abstract

:
The radiated flow of magnetized viscous fluid subject to the viscous dissipation phenomenon is numerically studied. The radiative phenomenon is addressed with nonlinear relations. Further, analysis is performed by using the slip effects and convective thermal flow constraints. The transformed problem is numerically evaluated using the Keller Box method. The physical parameter effects, such as the magnetic parameter for the velocity profile, Prandtl number, Brownian motion parameter and Biot number for the energy profile and Lewis number, and the thermophoresis parameter for the concentration profile are discussed. The obtained results suggest applications in enhancing the heat transfer phenomenon, thermal system, energy generation, heat transmission devices, power generation, chemical reactions, etc.

1. Introduction

The interaction of magnetic forces has significance in engineering applications, such as machine cooling, the lubrication industry, oil pumping, etc. The concept is associated with the impact of electromagnetic body forces on fluid flows. The primary role of magnetic force is in the suction phenomenon, wall motion and blowing. The flow of electrically conducted fluids is usually considered to make the process more effective. The impact of an externally applied magnetic field of a certain strength causes electrical conductivity in the flow of the considered fluids. This is the traditional type of MHD flow mechanism. In sea water, an externally applied magnetic field has much lower thermal conductivity, so an external electrical field is applied to overcome this problem to control the flow. An intersected electric and magnetic field causes the production of a wall-parallel Lorentz force with the objective of varying the procedure of the pressure gradient driven in the boundary layer. Bhatti et al. [1] evaluated the magnetized flow of nanofluids with the enhancement of the solar energy framework. Sabu et al. [2] discussed the particles’ shapes for MHD flow in a rotating flow with the passive control phenomenon. The Williamson nanofluid flow with MHD effects was examined by Song et al. [3]. Jie et al. [4] reported sinusoidal wave flows due to oscillations under a magnetic force impact. Wang et al. [5] performed a bioconvective investigation with the slip phenomenon created by Maxwell liquids.
Nanotechnology is extensively used in industry because nanometer-sized particles provide a number of physical and chemical resources. The term nanofluid was first introduced by Choi [6] in 1995. “Nanofluids” are beneficial for reducing abrasion and reducing in operating components, including pumps and compressors, and result in gasoline savings of more than 6%. It is quite possible that great financial savings can ultimately be made. The appearance of nanofluids is the focal point of research on the sliding of nanofluids in the presence of nanoparticles. Nanofluids consist of particles that are nanometer-sized that slide inside the base fluids. Nanofluids are artificially created colloids containing a base fluid and nanoparticles. Nanoparticles are widely used in improving heat transfer systems, thermal energy, new fuels, hybrid engines, prescription drugs, electricity production and fuel reduction and most cancer remedies. Convection heat transfer fluids, as well as mixtures of oil, ethylene glycol and water, have thermal conductivity, which is important for the constant switching between heat transfer and surface heat transfer. This latest form of liquid was first discovered by increases in thermal conductivity in conventional liquids. The heat transfer mechanism in nanofluid is imprecise, but predictable, as it clearly relies on various additional structures, such as surfactant effects, particle size and particle dispersion of dispersed debris. In addition to progress, nanofluids are regularly used in biomedicine for the identification of tumor cells using nano-scale preparation methods and delivery systems, and breaking blockages in the circulatory system in supply pathways in thallium research. Rasheed et al. [7] modeled the curved surface flow problem with Prandtl nanofluid via convective constraints. Zeeshan [8] reported the entropy-generation significance for viscous nanomaterials with a nonlinear radiation impact. Rasheed et al. [9] evaluated Jeffrey nanofluid due to a vertically moving cylinder with dissipation impact. Zeeshan et al. [10] predicted a rotating surface with a thin film of nanoparticles. Hybrid nanofluid analysis with ethylene glycol as the base liquid was explored by Zeeshan et al. [11].
The radiative phenomenon is an interesting mode of heat transfer consisting of the transfer of energy as electromagnetic waves. At high temperatures, the radiation effects are important and have different applications in thermal systems, heat transfer devices, chemical processes, fusion and fission reactions, atomic blasts, nuclear systems, etc. In solar power, gas turbines, aircrafts, space vehicles and missile technology, the application of the radiative phenomenon is also observed. Zhao et al. [12] addressed the features of thermal radiation associated with street canyons. Ramesh et al. [13] presented the decomposition of the thermophoretic phenomenon for carbon nanotubes with radiative analysis. Famakinwa et al. [14] used nonlinear thermal radiation for examining the parallel plate flow of a hybrid nanomaterial with variable features. Ketchate1 et al. [15] observed the onset of nonlinear thermal radiation for aluminum oxide and water-based nanoparticles in porous channels. Some modern developments regarding the fluid flow and nanoparticle concepts are listed in Refs. [16,17,18,19,20,21].
The aim of the current work is to report the stagnation point flow of nanofluid due to a moving surface with magnetic force and porous medium effects. The heat transfer phenomenon is observed by considering the viscous dissipation and radiative phenomenon. The nonlinear thermal radiation relationships are followed, instead of the traditional linear radiation phenomenon. The flow is subject to slip effects. Further, thermal analysis is performed using the convective boundary conditions. For the formulated flow system, the Keller Box numerical scheme is implemented to compute the numerical simulations. It is noted that the Keller Box scheme is widely used in the flow problem; however, numerical simulations for the nanofluid problem with these flow features have not yet been reported. The heat and mass transfer rate and wall shear stress are assessed using a useful resource of plotting the Sherwood Number and Nusselt number relative to several physical elements.

2. Mathematical Modeling

A two-dimensional, steady boundary layer incompressible flow of nanofluid is considered over a stretching sheet in the presence of solar radiation (thermal radiation). Flow occurs due to the linear velocity of a surface, which is mathematically written as u w ( x ) = a x and is considered through a porous medium. The y -axis of the coordinate system should be perpendicular to the stretching sheet denoted as the x-axis. Assume that   B 0 is the uniform magnetic field strength, which is applied perpendicularly to the stretching sheet. T   represents the fluid temperature and C represents the concentration. T w is the fluid’s reference temperature, while T is the fluid’s ambient temperature. The effects of viscous dissipation are also incorporated during flow. Under these assumptions, the governing equations after using conservation law are modeled [8,15].
u x + v y = 0 ,
  u u x + v v y = u u x + v f 2 u y 2 σ e B 0 2 ρ f ( u u ) ν k p ( u u ) ,  
u T x + v T y = α 2 T y 2 + v f C f ( u y ) 2 1 ( ρ C ) f ( q r y ) + σ e B 0 2 ( ρ C ) f ( u u ) 2 + τ [ D B T y C y + D T T ( T y ) 2 ]  
u C x + v C y = D B 2 C y 2 + D T T 2 T y 2 ,
With boundary conditions given as:
y = 0   : u = u w ( x ) + L u y = a x + L u y ,   v = 0 , k T y = h ( T T w ) ,   C = C w ,
  y   : u u ( x ) = b x ,   T T ,   C C .  
The mathematical expression for radiative heat flux is given as follows [13,14]:
q r = 4 σ * 3 k * T 4 y = 16 σ * T 3 3 k * T y ,  
Therefore, Equation (3) becomes:
u T x + v T y = y [ ( α + 16 σ * T 3 3 ( ρ C ) f k * ) T y ] + v f C f ( u y ) 2 + σ e B 0 2 ( ρ C ) f ( u u ) 2 + τ [ D B T y C y + D T T ( T y ) 2 ]
In Equation (7) on the R.H.S., the first term is expressed as [14,15]:
  α y [ ( 1 + R d ( 1 + ( θ w 1 ) θ ) 3 ) T y ]
where R d = 16 σ * T 3 3 k k * is the radiation parameter. Equation (7) then becomes:
u T x + v T y = α y [ ( 1 + R d ( 1 + ( θ w 1 ) θ ) 3 ) T y ] + v f C f ( u y ) 2 + σ e B 0 2 ( ρ C ) f ( u u ) 2 + τ [ D B T y C y + D T T ( T y ) 2 ]
By using the similarity transformation
η = a υ y ,   u = a x f ( η ) ,   v = a v f f ,   θ ( η ) = T T T w T ,   φ ( η ) = C C C w C   ,
the continuity equation will be identically satisfied and Equations (2)–(4) can be written as:
f + f f { f } 2 + λ 2 + M ( λ f ) K f = 0  
1 P r [ θ ( 1 + R d ( 1 + ( θ w 1 ) θ ) 3 ) + 3 θ 2 ( θ w 1 ) R d ( 1 + ( θ w 1 ) θ ) 2 ] + f θ + E C * f 2 + M E C * ( λ f ) 2 + N b θ ϕ + N t θ 2 = 0  
ϕ + L e f φ + N t N b θ = 0
With the thermal slip boundary conditions:
f ( η ) = 0 ,   f ( η ) = 1 + β f ( η ) ,     θ ( η ) = B i ( 1 θ ( η ) ) ,     φ ( η ) = 1           as   η = 0
  η : f ( η ) λ , θ ( η ) 0 , φ ( η ) 0  
The involved physical parameters are defined as:
where M = σ B 0 ρ f a is the magnetic parameter, K = ν k p is the porosity parameter, β = L a v f represents the slip parameter, λ = b a represents the stretching parameter, B i = h k v f a is the Biot number, N b = τ D B ( C w C ) v f is the Brownian motion parameter, N t = τ D T ( T w T ) T v f is the thermophoresis parameter, E c * = U W 2 C p ( T w T ) is the local Eckert number, L e = v D B is the Lewis number and P r = v f a is the Prandtl number.

3. Keller Box Scheme

The nonlinear coupled system justified via Equations (10)–(12) is numerically solved with the Keller Box scheme. This scheme was used due to its high accuracy. First, the governing problem is reduced into the first-order equations
f = u
u = v
θ = p
φ = g
Equations (1)–(3) become:
v + f v u 2 + λ 2 + M ( λ u ) K u = 0
1 P r [ p ( 1 + R d ( 1 + ( θ w 1 ) θ ) 3 ) + 3 p 2 ( θ w 1 ) R d ( 1 + ( θ w 1 ) θ ) 2 ] + f p + N b p g + N t p 2 + E c * v 2 + M E c * ( λ u ) 2 = 0
g + L e f g + N t N b p = 0
Therefore, BC becomes:
f ( 0 ) = 0 ,       u ( 0 ) = 1 + β v ,       p ( 0 ) = B i [ 1 θ ( 0 ) ] ,       φ ( 0 ) = 1           at   η = 0
u ( ) λ ,       θ ( ) 0 ,       ϕ ( ) 0 at   η
By applying the Keller Box method:
f j f j 1 h j = u j + u j 1 2
f j f j 1 h j 2 ( u j + u j 1 ) = 0
u j u j 1 h j = v j + v j 1 2
u j u j 1 h j 2 ( v j + v j 1 ) = 0
θ j θ j 1 h j = p j + p j 1 2
θ j θ j 1 h j 2 ( p j + p j 1 ) = 0
φ j φ j 1 h j = g j + g j 1 2
φ j φ j 1 h j 2 ( g j + g j 1 ) = 0
Using Equations (10)–(13) in Equations (5)–(7), we get:
v j v j 1 + h j [ 1 4 ( f j + f j 1 ) ( v j + v j 1 ) ( u j + u j 1 2 ) 2 + λ 2 + M { λ ( u j + u j 1 2 ) } K 2 { u j + u j 1 } ] = 0
1 P r [ ( p j p j 1 h j ) ( 1 + R d ( 1 + ( θ w 1 ) ( θ j + θ j 1 2 ) ) 3 ) + 3 ( p j + p j 1 2 ) 2 ( θ w 1 ) R d ( 1 + ( θ w 1 ) ( θ j + θ j 1 2 ) ) 2 ] + ( f j + f j 1 2 ) ( p j + p j 1 2 ) + N b ( p j + p j 1 2 ) ( g j + g j 1 2 ) + N t ( p j + p j 1 2 ) 2 + E c * ( v j + v j 1 2 ) 2 + M E c * [ ( λ ( u j + u j 1 2 ) ) 2 ] = 0
p j p j 1 + h j ( 1 + R d ( 1 + ( θ w 1 ) ( θ j + θ j 1 2 ) ) 3 ) [ Pr 4 { ( f j + f j 1 ) ( p j + p j 1 ) + N b ( p j + p j 1 ) ( g j + g j 1 ) + N t ( p j + p j 1 ) 2 + E c * ( v j + v j 1 ) 2 + M E c * [ ( 2 λ ( u j + u j 1 ) ) 2 ] } + 3 ( p j + p j 1 2 ) 2 ( θ w 1 ) R d ( 1 + ( θ w 1 ) ( θ j + θ j 1 2 ) ) 2 ] = 0
g j g j 1 + h j 4 L e ( f j + f j 1 ) ( g j + g j 1 ) + N t N b ( p j p j 1 ) = 0

4. Newton Method for Linearization

Linearize the system of nonlinear equations:
f j ( i + 1 ) = f j ( i ) + δ f j ( i )
u j ( i + 1 ) = u j ( i ) + δ u j ( i )
v j ( i + 1 ) = v j ( i ) + δ v j ( i )
p j ( i + 1 ) = p j ( i ) + δ p j ( i )
θ j ( i + 1 ) = θ j ( i ) + δ θ j ( i )
φ j ( i + 1 ) = φ j ( i ) + δ φ j ( i )
Therefore, Equation (18) becomes:
δ f j δ f j 1 h j 2 [ δ u j + δ u j 1 ] = ( r 1 ) j
δ u j δ u j 1 h j 2 [ δ v j + δ v j 1 ] = ( r 2 ) j
δ θ j δ θ j 1 h j 2 [ δ p j + δ p j 1 ] = ( r 3 ) j
δ φ j δ φ j 1 h j 2 [ δ g j + δ g j 1 ] = ( r 4 ) j
where
( r 1 ) j = f j 1 f j + h j u j 1 2
( r 2 ) j = u j 1 u j + h j v j 1 2
( r 3 ) j = θ j 1 θ j + h j p j 1 2
( r 4 ) j = p j 1 p j + h j g j 1 2
Equation (14) becomes:
δ v j [ A 1 ] + δ v j 1 [ A 2 ] + δ f j [ A 3 ] + δ f j 1 [ A 4 ] + δ u j [ A 5 ] + δ u j 1 [ A 6 ] = r 5
Equation (15) becomes:
δ p j [ B 1 ] + δ p j 1 [ B 2 ] + δ f j [ B 3 ] + δ f j 1 [ B 4 ] + δ g j [ B 5 ] + δ g j 1 [ B 6 ] + δ v j [ B 7 ] + δ v j 1 [ B 8 ] = r 6
Equation (17) becomes:
δ g j [ C 1 ] + δ g j 1 [ C 2 ] + δ f j [ C 3 ] + δ f j 1 [ C 4 ] + δ p j [ C 5 ] + δ p j 1 [ C 6 ] = r 7
where
A 1 = 1 + h j 2 f j 1 2 ,   A 2 = 1 + h j 2 f j 1 2 , A 3 = A 4 = h j 2 v j 1 2 ,   A 5 = A 6 = 1 2 + M 2 λ K 2
B 1 = 1 + 1 2 f j 1 2 + N b 2 g j 1 2 + N t p j 1 2 ,   B 2 = 1 + 1 2 f j 1 2 + N b 2 g j 1 2 + N t p j 1 2
B 3 = B 4 = 1 2 p j 1 2 + M E c * ( f j 1 2 + λ ) ,   B 5 = B 6 = N b 2 p j 1 2 ,   B 7 = B 8 = E c * p v 1 2
C 1 = 1 + h j 2 L e f j 1 2 ,   C 2 = 1 h j 2 L e f j 1 2 ,   C 3 = C 4 = h j 2 L e g j 1 2 ,   C 5 = C 6 = N t N b
and
r 5 = v j 1 v j h j [ ( f j 1 2 ) ( v j 1 2 ) ( u j 1 2 ) 2 λ 2 M { λ 2 + ( u j 1 2 ) 2 λ ( u j + u j 1 ) } + K 2 ( u j + u j 1 )
r 6 = p j 1 p j +   P r 1 + R d ( 1 + ( θ w 1 ) ( θ j + θ j 1 2 ) ) 3 [ ( p j 1 2 ) ( f j 1 2 ) + E c * ( v j 1 2 ) 2 + M E c * { λ 2 + ( f j 1 2 ) 2 + λ ( f j + f j 1 ) } + N b 2 ( f j 1 2 ) ( g j 1 2 ) + N t ( p j 1 2 ) 2 ] + 3 ( θ w 1 ) R d ( 1 + ( θ w 1 ) ( θ j + θ j 1 2 ) ) 2 ( p j 1 2 ) 2
r 7 = g j 1 g j h j 2 L e ( g j 1 2 ) ( f j 1 2 ) N t N b ( p j p j 1 )
f ( 0 ) = 0 ,       u ( 0 ) = 1 + β v ,       p ( 0 ) = B i [ 1 θ ( 0 ) ] ,       φ ( 0 ) = 1               at   η = 0
u ( ) λ ,       θ ( ) 0 ,       ϕ ( ) 0         at   η

5. Computational Method

In this section, the arising PDEs that are nonlinear along with the BCs are solved by using the Keller Box method. For this, non-linear equations are first converted into ordinary differential equations of the first order. The conversion process is carried out by using the following transformation:
( y 1 , y 2 , y 3 , y 4 , y 5 , y 6 , y 7 , y 8 , y 9 , y 10 ) = ( f , f , f , f θ , θ , θ , , , )
( y 1 y 2 y 3 y 4 y 5 y 6 y 7 y 8 ) = ( 1 y 3 y 4 [ y 3 2 y 2 y 4 λ 2 M ( λ y 3 ) K y 4 ] y 6 1 1 + R d ( 1 + z y 4 ) 3 [ Pr { y 1 y 5 + Nby 5 y 6 + Nt ( y 5 ) 2 + E c * y 3 2 + M E c ( λ y 2 ) 2 } + 3 y 5 2 z R d ( 1 + z y 4 ) 2 ] y 8 ( L e y 2 y 8 + N t N b y 6 ) )
( y 1 ( 0 ) y 2 ( 0 ) y 3 ( 0 ) y 4 ( 0 ) y 5 ( 0 ) y 6 ( 0 ) y 7 ( 0 ) ) = ( u 1 1 + β f λ u 2 B i ( 1 y 4 ) 1 u 3 ) .
The above system of first-order ODEs and corresponding ICs is resolved using the Keller Box method. The corresponding values of the mysterious ICs u 1 , u 2   a n d   u 3 are considered by Newton’s approach to BCs f ( η ) λ , θ ( η ) 0 , ϕ ( η ) 0   a s   n and the boundary conditions are satisfied. The calculations were conducted using the mathematical software MATLAB.

6. Results and Discussion

A two-dimensional, steady boundary layer flow of MHD nanofluid through an extending sheet with the occurrence of thermal radiation was analyzed. We examined the graphical results of the problem from an actual perspective. The variations in speed, temperature and nanoparticles concentration with governing parameters, for example N b , the parameter for Brownian motion; the Prandtl number P r ,   N t , the parameter for thermophoresis; the Biot number Bi and the magnetic parameter ( M ), were determined using the Keller Box technique and the mathematical calculations were processed utilizing the MATLAB program. The outcomes are presented in Tables and the diagrams displayed in Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12.
Figure 1 shows that the velocity in the boundary area was reduced by the slip parameter. Figure 2 shows that, by increasing the values of K , the velocity profile declined, i.e., with an increase in the value of K , the boundary layer thickness increased, so the velocity decreased with the increase in the value of K because the effect of the porous medium, which would act against the flow, also increased and enhanced the deceleration of the flow.
The results of magnetic field parameter M over the dimensionless velocity and temperature are displayed in Figure 3 and Figure 4, respectively. Figure 4 shows the effect of M on the velocity profile; it caused the velocity profile to decrease with a decrease in force, which resulted in a reduction in speed. In addition, higher values of M decreased the boundary layer thickness. In the Lorentz force, expansion was caused by the magnetic field; subsequently, the highly restricting power made it difficult for the liquid to stream without any hindrance. The magnetic field acted against the vehicle cycle. Figure 4 shows the impact of the magnetic field on the temperature profiles, which increased with an increase in M . The Lorentz force tended to increase the temperature in nanofluid movement. Thus, the thermal boundary layer thickness increased for a more grounded magnetic field. Figure 5 illustrates the impact of stretching boundary λ on the speed profile. It is apparent that expanding the upper side resulted in an increase in the speed field. However, expanding the upper side decreased the thickness of the force and limited the layers for both the non-Newtonian and Newtonian liquid cases.
Figure 6 and Figure 7 illustrate the impact of Brownian motion parameter N b on the temperature and nanoparticle volume concentration profile, respectively. Figure 6 presents the effects of different values of N b on the thermal boundary layer. At this phase, Brownian motion occurred due to the size of the nanoparticles in the nanofluids. In the transfer of heat, the motion of the particles in the fluid played a significant role. In the flow, the nanoparticles’ motion was affected by an increase in N b . The temperature of the nanofluid increased when the kinetic energy of the nanoparticles was increased by the strength of this chaotic motion. Figure 7 illustrates the variations in the nanoparticle concentration with N b . For stronger Brownian motion, a slight decrease in φ occurred accounted away from the sheet. The impact of thermophoresis parameter N t on the temperature profile and the concentration profile is presented in Figure 8 and Figure 9. Figure 8 illustrates the impact of the thermophoresis parameter of nanoparticles and thermal radiation on the thermal boundary layer. Additionally, it shows that the expansion in the thermophoresis boundary resulted in an increase in the thermophoretic force that transported nanoparticles from warm to cool spaces; however, the temperature profile increased due to this. Figure 9 shows the impact of N t on the nanoparticle concentration profile. With an increase in N t , there was an increase in the wideness of the concentration boundary and overshoot formed near the wall. Figure 10 shows the impact of the Prandtl number on the nanofluid’s thermal boundaries. Practically, if P r increased, there would be a decrease in the thermal diffusivity, and this principle would diminish the ability to transfer energy, decreasing the thermal boundary layer. Figure 11 depicts the impact of B i on the temperature profile. As the Biot number increased, the thermal boundary increased. Figure 12 illustrates the impact of Le on the nanoparticle concentration profile. At higher values of Le, the concentration of the nanoparticles greatly decreased. As L e increased, the thermal boundary layer thickness also increased; however, the thickness of the concentration boundary layer decreased.
Table 1 indicates that, for a weaker Brownian movement, N b = 0.1 , the Nusselt number decreased with an increase in Nt. As the Brownian motion increased (for example, N b changed from 0.1 to 0.4), the decrease in the Nusselt number became critical with an expansion in N t . This is reasonable as the strengthened Brownian movement was comparable to the instant movement of nanoparticles that stretched the wall to the dormant fluid. The heat movement rate at the wall θ’(0) was an increasing function of the radiation boundary. It was additionally seen that, for a given N b value, the decrease in Nusselt number with N t was greater in the case of R d = 0 when compared with R d = 1 . These observations were obtained from the results of Zeeshan [8]. Table 2 provides information for the decreased Sherwood number. The Sherwood number decreased when N t increased. This was because of the way that stronger thermophoresis impacts would generally move the nanoparticle volume fraction away from the surface, as the hot surface repulsed the sub-micron-sized particles, subsequently forming a moderately molecule-free layer close to the surface. For weaker Brownian movements and stronger Brownian movements, this decrease was constant.

7. Conclusions

This research is about the numerical examination of the impacts of thermal radiation on the MHD nanofluid flow when an extending sheet that is not linear is saturated in a permeable medium. The Keller Box method was used to obtain a numerical solution for governing non-linear equations after transformation by using similarity transformation. The mass and heat transfer rates, shear stress and other factors, such as the nanoparticle concentration profile, temperature and velocity, were computed numerically by using MATLAB, and their graphical representation is also provided. The MHD nanofluid flow was used to investigate the effects of slip with velocity when thermal radiation is applied. Some of the observations of interest in this study are
The increase in the stretching sheet parameter suggests an increase in the velocity profile;
The increase in the parameter of velocity slip decreased the velocity profile;
The increase in the porosity parameter reduced the velocity profile;
When the Brownian motion parameter increased, the nanoparticle concentration decreased, whereas the temperature increased;
The nanoparticle concentration and temperature both increased with an increase in the thermophoresis parameter;
By increasing the Prandtl number, the temperature profile decreased;
The increase in the Biot number indicated an increase in the temperature profile;
As L e increased, the concentration decreased;
By increasing the Brownian motion parameter, the Sherwood number increased in contrast to the increase in the thermophoresis parameter;
These results can be further extended for different non-Newtonian fluid models, such as Jeffrey fluid, viscoelastic fluid, couple stress fluid, Maxwell fluid, Oldroyd-B fluid, etc. Moreover, some other features can be discussed for current flow problems, such as entropy generation, heat source, activation energy, bioconvection, etc.

Author Contributions

All authors are equally contributed in the mathematical modeling, data analysis and final drafting. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: 22UQU4340474DSR15. Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R163), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data are fully available in the manuscript.

Acknowledgments

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: 22UQU4340474DSR15. Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2022R163), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a Stretching rate T Fluid temperature
C Concentration T w Surface temperature
T Free stream temperature u Free stream velocity
( u , v ) Velocity components ( x , y ) Cartesian coordinates
σ e Electrically conductivity B 0 Magnetic force strength
T Fluid temperature C Fluid concentration
D B Brownian diffusion q r Radiative flux
v f Fluid viscosity α Thermal conductivity
L Slip factor M Magnetic parameter
R d Radiation parameter β Slip parameter
B i Biot number N b Brownian motion parameter
K Porosity parameter λ Stretching parameter
N t Thermophoresis parameter L e Lewis number and
E c * Local Eckert number P r Prandtl number

References

  1. Bhatti, M.M.; Bég, O.A.; Abdelsalam, S.I. Computational framework of magnetized MgO-Ni/water-based stagnation nanoflow past an elastic stretching surface: Application in solar energy coatings. Nanomaterials 2022, 12, 1049. [Google Scholar] [CrossRef] [PubMed]
  2. Sabu, A.S.; Wakif, A.; Areekara, S.; Mathew, A.; Shah, N.A. Significance of nanoparticles’ shape and thermo-hydrodynamic slip constraints on MHD alumina-water nanoliquid flows over a rotating heated disk: The passive control approach. Int. Commun. Heat Mass Transf. 2021, 129, 105711. [Google Scholar] [CrossRef]
  3. Song, Y.-Q.; Hamid, A.; Sun, T.-C.; Khan, M.I.; Qayyum, S.; Kumar, R.N.; Prasannakumara, B.C.; Khan, S.U.; Chinram, R. Unsteady mixed convection flow of magneto-Williamson nanofluid due to stretched cylinder with significant non-uniform heat source/sink features. Alexandria Eng. J. 2022, 61, 195–206. [Google Scholar] [CrossRef]
  4. Jie, Z.; Khan, M.I.; Al-Khaled, K.; El-Zahar, E.R.; Acharya, N.; Raza, A.; Khan, S.U.; Xia, W.-F.; Tao, N.-X. Thermal transport model for Brinkman type nanofluid containing carbon nanotubes with sinusoidal oscillations conditions: A fractional derivative concept. Waves Random Complex Media 2022, 32, 1–20. [Google Scholar] [CrossRef]
  5. Wang, F.; Ahmad, S.; Al Mdallal, Q.; Alammari, M.; Khan, M.N.; Rehman, A. Natural bio-convective flow of Maxwell nanofluid over an exponentially stretching surface with slip effect and convective boundary condition. Sci. Rep. 2022, 12, 2220. [Google Scholar] [CrossRef]
  6. Choi, S.U.S. Nanofluids: A new field of scientific research and innovative applications. Heat Transf. Eng. 2008, 29, 429–431. [Google Scholar] [CrossRef]
  7. Rasheed, H.U.; Islam, S.; Zeeshan; Khan, J.; Abbas, T.; Mohmand, M.I. Numerical solution of chemically reactive and thermally radiative MHD Prandtl nanofluid over a curved surface with convective boundary conditions. Z. Angew. Math. Mech. 2021, e202100125. [Google Scholar] [CrossRef]
  8. Zeeshan. Second law and entropy generation analysis of magnetized viscous fluid flow over a permeable expandable sheet with nonlinear thermal radiation: Brownian and thermophoresis effect. Adv. Mech. Eng. 2022, 14, 1–9. [Google Scholar] [CrossRef]
  9. Rasheed, H.U.; Al-Zubaidi, A.; Islam, S.; Saleem, S.; Khan, Z.; Khan, W. Effects of Joule heating and viscous dissipation on magnetohydrodynamic boundary layer flow of jeffrey nanofluid over a vertically stretching cylinder. Coatings 2021, 11, 353. [Google Scholar] [CrossRef]
  10. Zeeshan; Rasheed, H.U.; Khan, W.; Khan, I.; Alshammari, N.; Hamadneh, N. Numerical computation of 3D Brownian motion of thin film nanofluid flow of convective heat transfer over a stretchable rotating surface. Sci. Rep. 2022, 12, 2708. [Google Scholar] [CrossRef]
  11. Zeeshan; Khan, I.; Weera, W.; Mohamed, A. Heat transfer analysis of Cu and Al2O3 dispersed in ethylene glycol as a base fluid over a stretchable permeable sheet of MHD thin-film flow. Sci. Rep. 2022, 12, 8878. [Google Scholar] [CrossRef] [PubMed]
  12. Zhao, H.; Xu, G.; Shi, Y.; Li, J.; Zhang, Y. The characteristics of dynamic and non-uniform thermal radiation experienced by pedestrians in a street canyon. Build. Environ. 2022, 222, 109361. [Google Scholar] [CrossRef]
  13. Ramesh, G.; Madhukesh, J.; Shah, N.A.; Yook, S.-J. Flow of hybrid CNTs past a rotating sphere subjected to thermal radiation and thermophoretic particle deposition. Alexandria Eng. J. 2022, in press. [CrossRef]
  14. Famakinwa, O.A.; Koriko, O.K.; Adegbie, K.S. Effects of viscous dissipation and thermal radiation on time dependent incompressible squeezing flow of CuO-Al2O3water hybrid nanofluid between two parallel plates with variable viscosity. J. Comput. Math. Data Sci. 2022, 5, 100062. [Google Scholar] [CrossRef]
  15. Ketchate, C.G.N.; Kapen, P.T.; Fokwa, D.; Tchuen, G. Stability analysis of mixed convection in a porous horizontal channel filled with a Newtonian Al2O3/Water nanofluid in presence of magnetic field and thermal radiation. Chin. J. Phys. 2022, 79, 514–530. [Google Scholar] [CrossRef]
  16. Abbasi, A.; Farooq, W.; Tag-ElDin, E.S.M.; Khan, S.U.; Khan, M.I.; Guedri, K.; Elattar, S.; Waqas, M.; Galal, A.M. Heat transport exploration for hybrid nanoparticle (Cu, Fe3O4)-based blood flow via tapered complex wavy curved channel with slip features. Micromachines 2022, 13, 1415. [Google Scholar] [CrossRef]
  17. Waqas, H.; Oreijah, M.; Guedri, K.; Khan, S.U.; Yang, S.; Yasmin, S.; Khan, M.I.; Bafakeeh, O.T.; Tag-ElDin, E.S.M.; Galal, A.M. Gyrotactic motile microorganisms impact on pseudoplastic nanofluid flow over a moving Riga surface with exponential heat flux. Crystals 2022, 12, 1308. [Google Scholar] [CrossRef]
  18. Shahid, M.; Javed, H.M.A.; Ahmad, M.I.; Qureshi, A.A.; Khan, M.I.; Alnuwaiser, M.A.; Ahmed, A.; Khan, M.A.; Tag-ElDin, E.S.M.; Shahid, A.; et al. A brief assessment on recent developments in efficient electrocatalytic Nitrogen reduction with 2D non-metallic nanomaterials. Nanomaterials 2022, 12, 3413. [Google Scholar] [CrossRef]
  19. Manzoor, N.; Qasim, I.; Khan, M.I.; Ahmed, M.W.; Guedri, K.; Bafakeeh, O.T.; Tag-Eldin, E.S.M.; Galal, A.M. Antibacterial applications of low pressure plasma on degradation of multidrug resistant V. Cholera. Appl. Sci. 2022, 12, 9737. [Google Scholar] [CrossRef]
  20. Mamatha, S.U.; Devi, R.L.V.R.; Ahammad, N.A.; Shah, N.A.; Rao, B.M.; Raju, C.S.K.; Khan, M.I.; Guedri, K. Multi-linear regression of triple diffusive convectively heated boundary layer flow with suction and injection: Lie group transformations. Int. J. Mod. Phys. B 2022, in press. [Google Scholar] [CrossRef]
  21. Kiranakumar, H.V.; Thejas, R.; Naveen, C.S.; Khan, M.I.; Prasanna, G.D.; Reddy, S.; Oreijah, M.; Guedri, K.; Bafakeeh, O.T.; Jameel, M. A review on electrical and gas-sensing properties of reduced graphene oxide-metal oxide nanocomposites. Biomass Convers. Biorefinery 2022, in press. [CrossRef]
Figure 1. Impact of β on the velocity profile.
Figure 1. Impact of β on the velocity profile.
Micromachines 13 01839 g001
Figure 2. Impact of K on the velocity profile.
Figure 2. Impact of K on the velocity profile.
Micromachines 13 01839 g002
Figure 3. Impact of M on the velocity profile.
Figure 3. Impact of M on the velocity profile.
Micromachines 13 01839 g003
Figure 4. Impact of M on the temperature profile.
Figure 4. Impact of M on the temperature profile.
Micromachines 13 01839 g004
Figure 5. Impact of λ on the velocity profile.
Figure 5. Impact of λ on the velocity profile.
Micromachines 13 01839 g005
Figure 6. Impact of N b on the temperature profile.
Figure 6. Impact of N b on the temperature profile.
Micromachines 13 01839 g006
Figure 7. Impact of N b on the concentration profile.
Figure 7. Impact of N b on the concentration profile.
Micromachines 13 01839 g007
Figure 8. Impact of N t on the temperature profile.
Figure 8. Impact of N t on the temperature profile.
Micromachines 13 01839 g008
Figure 9. Impact of N t on the concentration profile.
Figure 9. Impact of N t on the concentration profile.
Micromachines 13 01839 g009
Figure 10. Impact of P r on the temperature profile.
Figure 10. Impact of P r on the temperature profile.
Micromachines 13 01839 g010
Figure 11. Impact of B i on the temperature profile.
Figure 11. Impact of B i on the temperature profile.
Micromachines 13 01839 g011
Figure 12. Impact of L e on the concentration profile.
Figure 12. Impact of L e on the concentration profile.
Micromachines 13 01839 g012
Table 1. Values of the decreased Nusselt number θ ( 0 ) when P r = 7 ,   M = 0.5 ,   λ = 0.5 ,   E c = 0.2 ,   θ w = 1.5 ,   L e = 10   and B i = 0.1 .
Table 1. Values of the decreased Nusselt number θ ( 0 ) when P r = 7 ,   M = 0.5 ,   λ = 0.5 ,   E c = 0.2 ,   θ w = 1.5 ,   L e = 10   and B i = 0.1 .
N b R d N t
0.10.20.30.4
0.10
1
0.07543
0.3597
0.07418
0.3584
0.07271
0.3570
0.07095
0.3556
0.20
1
0.06623
0.3486
0.06365
0.3469
0.06042
0.3452
0.05623
0.3434
0.30
1
0.05249
0.3362
0.04712
0.3341
0.04013
0.3319
0.03023
0.3296
0.40
1
0.03307
0.3225
0.02328
0.3200
0.01018
0.3173
−0.006327
0.3145
Table 2. Values of the decreased Sherwood number φ ( 0 ) when P r = 7 ,   M = 0.5 ,   λ = 0.5 ,   E c = 0.2 ,   θ w = 1.5 ,   L e = 10   and B i = 0.1 .
Table 2. Values of the decreased Sherwood number φ ( 0 ) when P r = 7 ,   M = 0.5 ,   λ = 0.5 ,   E c = 0.2 ,   θ w = 1.5 ,   L e = 10   and B i = 0.1 .
N b R d N t
0.10.20.30.4
0.10
1
2.43
2.37
2.535
2.396
2.65
2.424
2.783
2.456
0.20
1
2.417
2.364
2.5
2.384
2.597
2.406
2.714
2.429
0.30
1
2.418
2.363
2.504
2.381
2.609
2.401
2.737
2.422
0.40
1
2.424
2.362
2.518
2.380
2.630
2.399
2.756
2.419
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bafakeeh, O.T.; Ahmad, B.; Noor, S.; Abbas, T.; Khan, S.U.; Khan, M.I.; Elattar, S.; Eldin, S.M.; Oreijah, M.; Guedri, K. Nonlinear Thermal Diffusion and Radiative Stagnation Point Flow of Nanofluid with Viscous Dissipation and Slip Constrains: Keller Box Framework Applications to Micromachines. Micromachines 2022, 13, 1839. https://doi.org/10.3390/mi13111839

AMA Style

Bafakeeh OT, Ahmad B, Noor S, Abbas T, Khan SU, Khan MI, Elattar S, Eldin SM, Oreijah M, Guedri K. Nonlinear Thermal Diffusion and Radiative Stagnation Point Flow of Nanofluid with Viscous Dissipation and Slip Constrains: Keller Box Framework Applications to Micromachines. Micromachines. 2022; 13(11):1839. https://doi.org/10.3390/mi13111839

Chicago/Turabian Style

Bafakeeh, Omar T., Bilal Ahmad, Skeena Noor, Tasawar Abbas, Sami Ullah Khan, Muhammad Ijaz Khan, Samia Elattar, Sayed M. Eldin, Mowffaq Oreijah, and Kamel Guedri. 2022. "Nonlinear Thermal Diffusion and Radiative Stagnation Point Flow of Nanofluid with Viscous Dissipation and Slip Constrains: Keller Box Framework Applications to Micromachines" Micromachines 13, no. 11: 1839. https://doi.org/10.3390/mi13111839

APA Style

Bafakeeh, O. T., Ahmad, B., Noor, S., Abbas, T., Khan, S. U., Khan, M. I., Elattar, S., Eldin, S. M., Oreijah, M., & Guedri, K. (2022). Nonlinear Thermal Diffusion and Radiative Stagnation Point Flow of Nanofluid with Viscous Dissipation and Slip Constrains: Keller Box Framework Applications to Micromachines. Micromachines, 13(11), 1839. https://doi.org/10.3390/mi13111839

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