Numerical exploration of radiative heat transfers in a magneto bioconvection Maxwell fluid passing through a spinning disk

Abstract

When applied to a rotating disk, Maxwell nanofluids have a wide range of applications in various fields. They provide increased heat transfer efficiency in cooling systems, which is essential for preserving the ideal operating temperatures of spinning gears like disk brakes and turbines. Furthermore, Maxwell nanofluids in microfluidic devices provide improved fluid manipulation and control, improving performance in applications such as microscale pumps and lab-on-a-chip systems. This article has numerically examined the magneto-bio-convection flow of Maxwell, which includes nanofluid, past a disk surface based on these applications. We present the model equations in PDE format, then shift them to ODEs using the appropriate variables. We utilize the bvp4c approach to obtain numerical solutions to the modeled equations. We also take into account the effects of thermal radiation, heat sources, Brownian motion, thermophoresis, chemical reactivity, and activation energy. We find that a higher stretching/shrinking variable has enhanced the radial velocity profile, while simultaneously diminishing the axial and angular velocity field. While the velocity profiles in the axial direction and redial direction have reduced with larger magnitude of Maxwell fluid variable. The thermal distributions have risen due to higher magnetic, Brownian motion, thermophoresis, heat sources, and thermal radiation components. Thermophoresis, magnetic and thermal radiation, heat sources, and Brownian motion all contribute to an increase in heat transfer rates.

Keywords

Maxwell fluid MHD heat source Brownian motion thermal radiation thermophoresis chemical reaction activation energy

Introduction

A Maxwell fluid is a type of viscoelastic fluid, which means it displays both viscous as well as elastic features. As Maxwell fluid flow, characterized by its behavior governed by the Maxwell model, manifests in a dynamic relationship between elasticity and viscosity, profoundly influencing velocity profiles and heat transfer characteristics.¹ Abdelsalam et al.² studied comparatively the rheological features of convective Maxwell flow fluid on a porous sheet. Yasin et al.³ inspected Maxwell fluid flow on a sheet with slip constraints subject to Hall effects and Darcy’s law. Contrasting to Newtonian fluids, Maxwell fluid exhibits time-dependent responses to applied stresses, transitioning from elastic to viscous behavior over a characteristic relaxation time was explored by Muhammad et al.⁴ This transition from elastic to viscous behavior results in unique velocity distributions compared to purely viscous fluids. In regions of high stress or strain rates where deformation occurs rapidly, the elastic component dominates, causing the fluid to resist deformation and store energy. Consequently, the velocity distribution deviates from the parabolic profile typically observed in Newtonian fluids, with higher velocities observed at the center of the flow channel and reduced velocities near solid boundaries. This alteration in velocity distribution arises from the elastic properties of the Maxwell fluid, which resist deformation near boundaries, leading to velocity gradients that differ from those predicted by classical fluid dynamics.⁵ Additionally, the time-dependent nature of Maxwell fluid flow introduces complexities in velocity evolution, with the fluid transitioning between elastic and viscous responses.⁶ This dynamic behavior can influence flow stability, turbulence characteristics, and the development of velocity profiles over time. Faizan et al.⁷ conducted research on the Maxwell nanofluid as it traveled across a stretching/shrinking disc in the direction of a stagnation point. Reddy et al.⁸ investigated the flow characteristics of three distinct nanofluids on a stretching surface in the presence of a non-uniform heat source, as well as the influence of thermal radiation.

Magnetohydrodynamics (MHD) is a multidisciplinary field that examines electrically conducted fluids, with magnetic fields effects. In MHD, the fluid’s motion is affected not only by conventional fluid dynamics but also by electromagnetic forces arising from interactions between the fluid’s motion and the magnetic field.⁹ These electromagnetic forces significantly impact velocity and temperature distributions within the fluid. For instance, when a conducted fluid flows in a magnetic field, it experiences Lorentz forces, which can alter the fluid’s velocity profile by inducing drag forces and modifying flow patterns.¹⁰ Furthermore, MHD effects can influence temperature distributions through joule heating, where electrical currents induced in the fluid by the magnetic field dissipate energy as heat.¹¹ This phenomenon can lead to non-uniform temperature distributions within the fluid, with regions experiencing stronger magnetic fields exhibiting higher temperatures due to increased Joule heating as observed by Bilal et al.¹² The impressions of MHD on velocity and thermal distributions are essential in several fields, like astrophysics,¹³ plasma physics,¹⁴ and engineering, where accurate characterization of fluid behavior under the influence of magnetic fields is critical for augmenting performance and attaining preferred consequences.¹⁵ Khan et al.¹⁶ conducted a comprehensive study on the production of entropy and thermal flow analysis for MHD fluid flow on a gyrating cylinder, considering the impacts of Joule heating. Their research investigated the multifaceted association among fluid dynamics, electromagnetic forces, and thermal effects in the presence of a magnetic field. Through their analysis, they demonstrated that the velocity profiles exhibited a notable decline, attributed to the influence of electromagnetic forces induced by magnetic field. Li et al.¹⁷ conducted an investigation into the numerical solution for heat transfer near a stagnation point of Maxwell nanofluid flow over a porous rotating disk. Reddy et al.¹⁸ investigated the MHD stagnation point flow of a Williamson hybrid nanofluid over a porous extended sheet.

Heat sources for fluid flow on a surface can arise from various mechanisms, including external heating sources like electrical heaters or solar radiation, as well as internal sources like chemical reactions or frictional heating. These heat sources play a significant role in influencing the thermal transference manner within the fluid.¹⁹ In regions where the heat source is present, the fluid temperature increases due to absorption of thermal energy, leading to the formation of thermal layers at boundaries characterized by higher temperatures near the surface.²⁰ This temperature gradient drives convective heat transfer, causing the heated fluid to transfer thermal energy to cooler regions through advection and diffusion. Additionally, the existence of a heat source induces changes in flow patterns and velocity distributions, further impacting the heat transfer process.^21,22 The influences of heat sources on fluid and thermal flows are crucial for various engineering applications.^23,24 Pasha et al.²⁵ studied the effects of nonlinear radiative magnetized fluid flow with impacts of heat sink, source, and chemical reactivity and observed that thermal distribution have augmented with growth in the strength of source factor. Sajid et al.²⁶ debated on trihybrid nanofluid flow on a porous sheet using radiation and thermal source factor. Kumar and Sharma²⁷ studied optimization of entropy for fluid flow on a gyrating disk with impacts of magnetic field and non-uniform thermal source and observed that thermal distribution has augmented with development in thermal source and magnetic factors.

Fluid flow with Brownian motion and thermophoresis encompasses a complex relationship that significantly influence particle dynamics and transport in various systems. Brownian motion, driven by haphazard collisions between suspended particles and surrounding fluid molecules, leads to erratic particle movement that affects dispersion and mixing. At higher temperatures, the kinetic energy of particles increases, intensifying Brownian motion and enhancing particle diffusion rates.²⁸ Consequently, temperature gradient induces changes in particle distribution and concentration profiles within the fluid.²⁹ Thermophoresis, meanwhile, involves the directed motion of particles in response to temperature gradients. This phenomenon is influenced by factors like particle mass, temperature gradient magnitude, and fluid properties.³⁰ In regions of temperature variation, thermophoretic forces act on particles, causing them to migrate toward zones of greater or lesser temperatures depending on their features.³¹ Temperature gradients also influence fluid properties like thickness and density, which in turn affect the thermophoretic behavior of particles. Furthermore, concentration gradients play a vital role in both Brownian motion and thermophoresis by influencing the spatial distribution and transport of particles within the fluid.³² Variations in particles’ concentration impact collision frequencies and thus the extent of Brownian motion, while concentration gradients also affect the thermophoretic mobility of particles. Moreover, interactions between Brownian motion, thermophoresis, and fluid flow leads to complex phenomena like particles’ accumulation, deposition, and dispersion, with implications for various industrial and natural processes ranging from nanoparticle transport in biological systems to aerosol dynamics in atmospheric science.³³ The coupled effects of Brownian motion, thermophoresis, temperature, and concentration gradients is crucial for precisely modeling and calculating particle behavior in complex fluid environments, with implications for fields including environmental engineering, nanotechnology, and biophysics.³⁴ Shahzad et al.³⁵ studied impacts of thermophoresis as well as Brownian motion on Darcy-Forchheimer flow of fluid through two disks with Cattaneo-Christov model. Many such study can be seen in References ^36–40.

Fluid flow with chemical reaction is a dynamic process where the movement of the fluid medium interacts intimately with chemical kinetics, resulting in distinct concentration profiles within the flow system of fluid.⁴¹ As the fluid flows, it transports reactants and products throughout the domain, influencing the distribution of fluid’s concentrations as noticed by Nasir et al.⁴² Convective transport, driven by fluid velocity and turbulence, accelerates the mixing of reactants chemically, facilitating collisions and promoting reaction rates in regions of high flow as observed by Li et al.⁴³ Chemical reactions occurring within the fluid alter the concentrations of fluid, leading to spatial variations in species concentrations over time.⁴⁴ Khan et al.⁴⁵ studied fluid flow on a movable needle with impacts of chemical reactivity Hall current and magnetic field. Many applications of fluid flow with impressions of chemical reactivity involve, cooling towers, catalysis and fog dispersion, etc. as discussed by He et al.⁴⁶ Factors like reaction rate constants, reactant concentrations, and temperature profoundly impact the extent and rate of chemical reactions, further shaping the concentration profiles. The interaction between fluid flow and chemical reaction dynamics results in augmentation of concentration profiles as observed by Razaq et al.⁴⁷ Mathematical modeling combined with physical approaches in case of fluid flow with impacts of chemical reactions, are essential tools to design of efficient and sustainable processes in fields like chemical engineering, environmental science, and materials science. Mishra et al.⁴⁸ simulated numerically nanofluid flow with significance of chemical reactions.

Nanofluids flows over a revolving disk have variety of uses. Nanofluids have the potential to upsurge the heat transfer rates in the field of thermal management, which makes them useful for cooling systems in spinning machinery such as industrial disk brakes or computer hard drives. Furthermore, nanofluids can enhance the lubricating qualities of lubrication systems, lowering wear and friction on rotating parts like gears and bearings. Furthermore, nanofluids in aerospace applications can help reduce drag and enhance the aerodynamic performance of propellers or spinning disks in aircraft engines. All things considered, a variety of industries have used nanofluids over spinning disks to improve heat transmission, lubrication, and aerodynamics. The current analysis investigates the laminar, time-independent and incompressible flow of Maxwell nanofluid over a gyratory disk surface. With the help of a rotating and stretching disk, we will study the thermal and solutal energy transport in the Von Kármán swirling flow of Maxwell nanofluids in this effort with various flow effects.

Problem formulation

Assume that the laminar and time-independent flow of Maxwell nanofluid over a spinning disk surface. A cylindrical coordinate system is assumed to analyze the flow analysis by assuming the u, v, w along r-, φ-, z-directions. Here, r, φ, and z are the radial, azimuthal, and axial directions. The stretching velocity of the disk is denoted by $u_{w} = cr$ where $c$ is the stretching constant. The disk swirls with an angular velocity $v = Ω r$ where $Ω$ , is the swirling rate (Figure 1). A magnetic field of strength $B_{0}$ is applied normal (z-direction) to the fluid flow. Furthermore, the effects of thermal radiation, heat source, thermophoresis, Brownian motion, chemical reaction, and Activation energy are taken into consideration. The main equation take the form as appended below^49–51:

\frac{\partial u}{\partial r} + \frac{u}{r} + \frac{\partial w}{\partial z} = 0,

(1)

\begin{matrix} u \frac{\partial u}{\partial r} + w \frac{\partial u}{\partial z} - \frac{1}{r} v^{2} = ν_{f} \frac{\partial^{2} u}{\partial z^{2}} - λ \\ (\begin{matrix} u^{2} \frac{\partial^{2} u}{\partial r^{2}} + 2 uw \frac{\partial^{2} u}{\partial r \partial z} + w^{2} \frac{\partial^{2} u}{\partial z^{2}} \\ - 2 \frac{1}{r} uv \frac{\partial v}{\partial r} - 2 \frac{1}{r} vw \frac{\partial v}{\partial z} + \frac{1}{r} v^{2} \frac{\partial u}{\partial r} + \frac{1}{r^{2}} u v^{2} \end{matrix}) \\ - \frac{σ_{f} B_{0}^{2}}{ρ_{f}} (u + λ w \frac{\partial u}{\partial z}), \end{matrix}

(2)

\begin{matrix} u \frac{\partial v}{\partial r} + w \frac{\partial v}{\partial z} + \frac{1}{r} uv = ν_{f} \frac{\partial^{2} v}{\partial z^{2}} - λ \\ (\begin{matrix} u^{2} \frac{\partial^{2} v}{\partial r^{2}} + 2 uw \frac{\partial^{2} v}{\partial r \partial z} + w^{2} \frac{\partial^{2} v}{\partial z^{2}} + 2 \frac{1}{r} uv \frac{\partial u}{\partial r} \\ + 2 \frac{1}{r} vw \frac{\partial u}{\partial z} - \frac{2}{r^{2}} u^{2} v - \frac{1}{r^{2}} v^{3} + \frac{1}{r^{2}} v^{2} \frac{\partial v}{\partial r} \end{matrix}) \\ - \frac{σ_{f} B_{0}^{2}}{ρ_{f}} (v + λ w \frac{\partial v}{\partial z}), \end{matrix}

(3)

\begin{matrix} u \frac{\partial T}{\partial r} + w \frac{\partial T}{\partial z} = \frac{k_{f}}{{(ρ C_{p})}_{f}} \frac{\partial^{2} T}{\partial z^{2}} + τ (D_{B} \frac{\partial T}{\partial z} \frac{\partial C}{\partial z} + \frac{D_{T}}{T_{\infty}} {(\frac{\partial T}{\partial z})}^{2}) \\ - \frac{1}{{(ρ C_{p})}_{f}} \frac{\partial q_{r}}{\partial z} + \frac{Q_{0}}{{(ρ C_{p})}_{f}} (T - T_{\infty}), \end{matrix}

(4)

\begin{matrix} u \frac{\partial C}{\partial r} + w \frac{\partial C}{\partial z} = D_{B} \frac{\partial^{2} C}{\partial z^{2}} + \frac{D_{T}}{T_{\infty}} \frac{\partial^{2} T}{\partial z^{2}} \\ - k_{r}^{2} (C - C_{\infty}) {(\frac{T}{T_{\infty}})}^{n} \exp (\frac{- E_{a}}{K^{*} T}), \end{matrix}

(5)

with boundary conditions:

\begin{matrix} u_{w} = cr, w = 0, v = Ω r, T = T_{w}, C = C_{w} at z = 0, \\ u \to 0, v \to 0, C \to C_{\infty}, T \to T_{\infty} as z \to \infty . \end{matrix}

(6)

Figure 1.

Graphical view of the flow phenomenon.

Above $q_{r}$ is described as:

q_{r} = - \frac{16 σ^{*} T_{\infty}^{3}}{3 k^{*}} \frac{\partial T}{\partial z} .

(7)

The similarity transformations are defined as⁵¹:

\begin{matrix} u = F (ξ) r Ω, v = G (ξ) r Ω, w = H (ξ) \sqrt{Ω ν_{f}}, \\ θ (ξ) = \frac{T - T_{\infty}}{T_{w} - T_{\infty}}, ϕ (ξ) = \frac{C - C_{\infty}}{C_{w} - C_{\infty}}, ξ = \sqrt{\frac{Ω}{ν_{f}}} z . \end{matrix}

(8)

By implementing equation (8) we have from above

H' (ξ) + 2 F (ξ) = 0,

(9)

\begin{matrix} F ″ (ξ) + G^{2} (ξ) - F^{2} (ξ) - H (ξ) F' (ξ) - M (F (ξ) + α F' (ξ) H (ξ)) \\ - α (F ″ (ξ) H^{2} (ξ) + 2 F (ξ) F' (ξ) H (ξ) - 2 G (ξ) G' (ξ) H (ξ)) = 0, \end{matrix}

(10)

\begin{matrix} G ″ (ξ) - 2 F (ξ) G (ξ) - G' (ξ) H (ξ) - M (G (ξ) + α G' (ξ) H (ξ)) \\ - α (H^{2} (ξ) G ″ (ξ) + 2 F (ξ) H (ξ) G' (ξ) + 2 G (ξ) F' (ξ) H (ξ)) = 0, \end{matrix}

(11)

\begin{matrix} (1 + \frac{4}{3} Rd) θ ″ (ξ) - \Pr H (ξ) θ' (ξ) + \Pr Nb θ' (ξ) ϕ' (ξ) \\ + \Pr Nt θ'^{2} (ξ) + Q \Pr θ (ξ) = 0, \end{matrix}

(12)

\begin{matrix} ϕ ″ (ξ) - ScH (ξ) ϕ' (ξ) + \frac{Nt}{Nb} θ ″ (ξ) - Sc K_{r} {(1 + δ θ (ξ))}^{n} \\ \exp (- \frac{E}{(1 + δ θ (ξ))}) ϕ (ξ) = 0, \end{matrix}

(13)

with constraints at boundary:

\begin{matrix} F (ξ = 0) = γ, G (ξ = 0) = 1, H (ξ = 0) = 0, \\ θ (ξ = 0) = 1, ϕ (ξ = 0) = 1, \\ F (ξ \to \infty) \to 0, G (ξ \to \infty) \to 0, θ (ξ \to \infty) \to 0, \\ ϕ (ξ \to \infty) \to 0 . \end{matrix}

(14)

In above equations, $M (= \frac{σ_{f} B_{0}^{2}}{ρ_{f} Ω})$ is magnetic factor, $Q (= \frac{Q_{0}}{{(ρ C_{p})}_{f} Ω})$ is heat source factor, $Rd (= \frac{4 σ^{*} T_{\infty}^{3}}{k_{f} k^{*}})$ thermal radiation factor, $α (= λ Ω)$ is Maxwell parameter, $Nt (= \frac{τ D_{T} (T_{w} - T_{\infty})}{ν_{f} T_{\infty}})$ is thermophoresis factor, $Nb (= \frac{τ D_{B} (C_{w} - C_{\infty})}{ν_{f}})$ is Brownian motion factor, $Sc (= \frac{ν_{f}}{D_{B}})$ is Schmidt number, $K_{r} (= \frac{k_{r}^{2}}{c})$ is chemical reaction factor, $\Pr (= \frac{ν_{f} {(ρ C_{p})}_{f}}{k_{f}})$ is Prandtl number, $E (= \frac{E_{a}}{K^{*} T_{\infty}})$ is Activation energy factor, and $δ (= \frac{T_{w} - T_{\infty}}{T_{\infty}})$ is temperature difference factor.

The Nusselt and Sherwood number are defined as:

\sqrt{{Re}_{r}} C_{f} = \frac{\sqrt{τ_{rz}^{2} + τ_{θ z}^{2}}}{ρ_{f}}, N u_{r} = - (1 + \frac{16 σ^{*} T_{\infty}^{3}}{3 k_{f} k^{*}}) \frac{r {\frac{\partial T}{\partial z} |}_{z = 0}}{(T_{w} - T_{\infty})} .

(15)

S h_{r} = - \frac{r {\frac{\partial C}{\partial z} |}_{z = 0}}{(C_{w} - C_{\infty})} .

(16)

Using equation (8), (15) and (16) are reduced as:

\sqrt{{Re}_{r}} Cf = \sqrt{f'^{2} (0) + g'^{2} (0)}, \frac{N u_{r}}{\sqrt{{Re}_{r}}} = - (1 + \frac{4}{3} Rd) θ' (ξ = 0),

(17)

\frac{S h_{r}}{\sqrt{{Re}_{r}}} = - ϕ' (ξ = 0) .

(18)

Here ${Re}_{r} (= \frac{r^{2} Ω}{ν_{f}})$ is the local Reynolds number.

Numerical solution

To solve the modeled equations by means of numerical method, we have chosen bvp4c MATLAB built-in function. This method is applicable to solve boundary value problems arising in physical and engineering fields. To ensure the accuracy of the bvp4c, error tolerance of 10⁻⁴ is chosen. To apply this method, let us assume that:

\begin{matrix} H (ξ) & = ϒ_{1}, H' (ξ) = ϒ_{1}', \\ F (ξ) & = ϒ_{2}, F' (ξ) = ϒ_{3}, F ″ (ξ) = ϒ_{3}', \\ G (ξ) & = ϒ_{4}, G' (ξ) = ϒ_{5}, G ″ (ξ) = ϒ_{5}', \\ θ (ξ) & = ϒ_{6}, θ' (ξ) = ϒ_{7}, θ ″ (ξ) = ϒ_{7}', \\ ϕ (ξ) & = ϒ_{8}, ϕ' (ξ) = ϒ_{9}, ϕ ″ (ξ) = ϒ_{9}' . \end{matrix}

(19)

Then

ϒ_{1}' = - 2 ϒ_{1},

(20)

ϒ_{3}' = - \frac{{\begin{matrix} - {(ϒ_{2})}^{2} + {(ϒ_{4})}^{2} - M (ϒ_{4} + α ϒ_{1} ϒ_{3}) - \\ ϒ_{1} ϒ_{3} - 2 α ϒ_{1} ϒ_{2} ϒ_{3} + 2 α ϒ_{4} ϒ_{5} ϒ_{2} \end{matrix}}}{(1 - α {(ϒ_{1})}^{2})},

(21)

ϒ_{5}' = - \frac{{\begin{matrix} - 2 ϒ_{1} ϒ_{2} - M (ϒ_{4} + α ϒ_{1} ϒ_{5}) - ϒ_{1} ϒ_{5} \\ - 2 α ϒ_{1} ϒ_{2} ϒ_{5} - 2 α ϒ_{1} ϒ_{3} ϒ_{4} \end{matrix}}}{(1 - α {(ϒ_{1})}^{2})},

(22)

ϒ_{7}' = - \frac{{- \Pr ϒ_{1} ϒ_{7} + \Pr Nb ϒ_{9} ϒ_{7} + \Pr Nt {(ϒ_{7})}^{2} + Q \Pr ϒ_{6}}}{(1 + 4 / 3 Rd)},

(23)

\begin{matrix} ϒ_{9}' = \\ - {- Sc ϒ_{1} ϒ_{9} + \frac{Nt}{Nb} ϒ_{7}' - Sc K_{r} {(1 + δ ϒ_{6})}^{n} \exp (- \frac{E}{(1 + δ ϒ_{6})}) ϒ_{8}}, \end{matrix}

(24)

with boundary conditions:

\begin{matrix} ϒ_{1} - 0, ϒ_{2} - γ, ϒ_{4} - 1, ϒ_{6} - 1, ϒ_{8} - 1, \\ ϒ_{2} - 0, ϒ_{4} - 0, ϒ_{6} - 0, ϒ_{8} - 0 . \end{matrix}

(25)

The applied numerical technique has the following advantages:

• A variety of BVPs, including stiff and highly nonlinear problems, may be solved using the robust solver bvp4c. It makes use of collocation and adaptive mesh refinement techniques to guarantee precise answers for a variety of issue kinds.

• In comparison to previous approaches, the bvp4c solver requires comparatively few function evaluations to solve BVPs, indicating its computational efficiency. Because of its efficiency, it may be used to solve problems involving a lot of equations or intricate geometries.

• By employing fourth-order collocation techniques, bvp4c offers very accurate solutions to BVPs. It uses error control techniques and automatically modifies the mesh density to make sure the calculated solution satisfies the required accuracy standards.

• bvp4c is a versatile and adaptable tool that can handle mixed boundary conditions, discontinuities, and singularities in systems of nonlinear differential equations.

Validation

To insure the accuracy, the present results are compared with Sparrow and Gregg⁵² for varying $\Pr$ when $M = 0$ , $α = 0$ , $Q = 0$ , $Rd = 0$ , $Nt = 0$ , $Nb = 0$ , $Sc = 0$ , $K_{r} = 0$ , and $E = 0$ . The present results are displayed in Table 1. We can see from this table that the present and past results are in great agreement which ensure the accuracy of the applied numerical method.

Table 1.

Validation of present results with published results by Sparrow and Gregg⁵² for varying $\Pr$ when $M = 0$ , $α = 0$ , $Q = 0$ , $Rd = 0$ , $Nt = 0$ , $Nb = 0$ , $Sc = 0$ , $K_{r} = 0$ , and $E = 0$ .

$\Pr$	Sparrow and Gregg⁵²	Present
1.0	0.3925	0.3926097547
10.0	1.1341	1.1341001232
100.0	2.6871	2.6871110345

Discussion of results

This section shows the obtained results from the present analysis and their physical explanation. The obtained results are presented in Figures 2 to 17 and Tables 2 and 3. The default values of factor used in the study are selected as $M = 0.5$ , $γ = 0.1$ , $α = 0.1$ , $Q = 0.5$ , $Rd = 0.2$ , $Nt = 0.1$ , $Nb = 0.1$ , $\Pr = 0.72$ , $Sc = 1.0$ , $K_{r} = 0.1$ , $E = 1.0$ , and $δ = 1.0$ . Figures 2 and 3 show the effect of stretching velocity factor ( $γ$ ) on radial velocity ( $F (ξ)$ ) and Azimuthal velocity ( $G (ξ)$ ) profiles, respectively. From these Figures, we can see that the greater $γ$ enhances ( $F (ξ)$ ) and decreases the velocity profiles in angular direction. The reason is that the stretching rate heightens in radial direction via higher ratio factor and the swirling rate reduces in angular direction due to higher ratio factor. Figures 4 and 5 demonstrate the influence of Maxwell parameter ( $α$ ) on radial velocity ( $F (ξ)$ ) and Azimuthal velocity ( $G (ξ)$ ) profiles, respectively. From these Figures, we can see that the greater $α$ reduces the velocity panels in radial and angular directions. These results have scientific justification since the fluid grows increasingly solid-like with larger values of $α$ due to the significant relaxation time value in the fluid. As a result, there is less movement that is fluid in every direction. Figures 6 to 8 show the effect ( $M$ ) on ( $F (ξ)$ ), ( $G (ξ)$ ), and $θ (ξ)$ profiles. We can see from these Figures that the greater $M$ lessens the velocity panels in radial and angular directions. Further, from Figure 8, we can see that the greater values of $M$ increases the thermal distribution. The Lorentz force, which works perpendicular to both the velocity and magnetic field vectors, is responsible for the drop in velocity profiles in both radial and tangential directions caused by the magnetic factor. The fluid flow is resisted by this force, which causes the velocity components in those directions to drop. Similarly, the higher Lorentz force enhances the heat transmission rate at the sheet surface which results higher temperature profile. The coefficient of heat production, denoted by $Q$ has a significant impact on the dimensionless temperature profile of the fluid, denoted by $θ (ξ)$ . According to Figure 9, it is obvious that as $Q$ increases, so does the temperature $θ (ξ)$ . The thermal layer width at boarder line will grow with upsurge in the fluid temperature $θ (ξ)$ . The external heat source has a major influence on thermal gradient of fluid, which ultimately causes an upsurge in $θ (ξ)$ as well as the thermal state of the fluid. The width of thermal layer at borderline grows to a greater extent when a substantial quantity of heat energy is created by the fluid particles themselves. Figure 10 illustrates the temperature fluctuations caused by radiation $Rd$ . It is observed that the $θ (ξ)$ of the Maxwell nanofluid climb as the amount of $Rd$ heightens. As $Rd$ shows that $Rd \propto \frac{1}{k^{*}}$ so increases the amount of $Rd$ means the reduction in $k^{*}$ causes more heat is transported to working fluid. In addition to this, temperature and thermal boundary layer thickness both upsurge as a result of larger values of the thermal radiation parameter. When the thermophoresis constraint $Nt$ growths, the $θ (ξ)$ also augment as shown in Figure 11. The reason behind this fact is that the fluid’s thermal resistance tends to increase due to the thermophoretic distribution, which causes the non-dimensional fluid temperature to rise. Specifically, the temperature difference is the root cause of this phenomenon. As the thermophoresis constraint $Nt$ upsurges, particles in the fluid migrate from the colder area to the warmer one which causes the temperature profile to climb. Here, the thermophoresis constraint $Nt$ is directly related to the temperature distribution. The nanoparticle concentration profile $ϕ (ξ)$ is a function of $Nt$ is shown in Figure 12. When $Nt$ is raised, then the concentration field is upswing. Physically $Nt$ is a key factor in the process of mass transmission, therefore $Nt$ is of vital importance. Nanoparticles are more likely to be displaced from hotter zone to colder one when $Nt$ is improved, leading to a rise in boundary layer width. Figure 13 portrays the effects of $Nb$ on $θ (ξ)$ . From the Figure it is seen that the $θ (ξ)$ profiles get more pronounced as the values of $Nb$ amplified. As $Nb$ enhances, the solid nanoparticles’ irregular mobility also intensifies. The dynamic process strengthens the thermal layer by converting the interior energy of nanoparticles into heat and kinetic energy. Thus, as seen in Figure 16, the thermal properties rise as the $Nb$ factor upturns. Figure 14 shows the effect of $Nb$ on $ϕ (ξ)$ . As the value of $Nb$ amplified, the profile $ϕ (ξ)$ declines. A positive value of $Nb$ boost the nanoparticles’ random mobility, which ultimately enhances fluid mixing as a result the solutal concentration, becomes weakens and reduces mass diffusion. Consequently heightened $Nb$ leads to increased random motion among nanoparticles. This enhanced mixing potentially yields a more uniform dispersion of nanoparticles and reduced concentration gradients, thereby thinning the concentration boundary layer and reducing its concentration. The behavior of $(K_{r} > 0)$ on concentration profiles $ϕ (ξ)$ is illustrated in Figure 15. As the chemical reaction $K_{r}$ grows, the reactants are utilized leading to a reduction in their concentration within the boundary. As shown in Figure 16, the activation energy $E$ effect on $ϕ (ξ)$ may be imagined. The graph clearly shows that when $E$ is made larger, the concentration profile rises. The Arrhenius function degrades as $E$ increases exponentially, leading to an uptick in the concentration field $ϕ (ξ)$ and the encouragement of the generative chemical process. Physically when a chemical reaction takes place at low temperatures and high activation energies, the reaction rate decreases, making the process go more slowly. The effect of $Sc$ on the concentration profile $ϕ (ξ)$ is revealed in Figure 17. This displays that upsurge in $Sc$ declines $ϕ (ξ)$ . The reason behind this is that the Schmidt number $(Sc = \frac{ν_{f}}{D_{B}})$ is a measure of the ratio between momentum diffusivity and mass diffusivity. There is a negative relationship between mass diffusivity and momentum diffusivity, which results in reduced scalar diffusivity, as the Schmidt number $Sc$ enhanced.

Figure 2.

Impact of $γ$ on $F (ξ)$ .

Figure 3.

Impact of $γ$ on $G (ξ)$ .

Figure 4.

Impact of $α$ on $F (ξ)$ .

Figure 5.

Impact of $α$ on $G (ξ)$ .

Figure 6.

Impact of $M$ on $F (ξ)$ .

Figure 7.

Impact of $M$ on $G (ξ)$ .

Figure 8.

Impact of $M$ on $θ (ξ)$ .

Figure 9.

Impact of $Q$ on $θ (ξ)$ .

Figure 10.

Impact of $Rd$ on $θ (ξ)$ .

Figure 11.

Impact of $Nt$ on $θ (ξ)$ .

Figure 12.

Impact of $Nt$ on $ϕ (ξ)$ .

Figure 13.

Impact of $Nb$ on $θ (ξ)$ .

Figure 14.

Impact of $Nb$ on $ϕ (ξ)$ .

Figure 15.

Impact of $K_{r}$ on $ϕ (ξ)$ .

Figure 16.

Impact of $E$ on $ϕ (ξ)$ .

Figure 17.

Impact of $Sc$ on $ϕ (ξ)$ .

Table 2.

Effects of $α, γ$ , and $M$ on $\sqrt{{Re}_{r}} Cf$ .

$α$	$γ$	$M$	$\sqrt{{Re}_{r}} Cf$
0.4			1.1031
0.9			1.2343
1.3			1.4019
1.9			1.5068
	0.2		1.5982
	0.7		1.6423
	1.1		1.7421
	1.5		1.8975
		0.8	1.9889
		1.4	2.1431
		1.7	2.1982
		2.1	2.3364

Table 3.

Effects of $Rd$ , $Q$ , $M$ , $Nb$ , and $Nt$ on $\frac{N u_{r}}{\sqrt{{Re}_{r}}}$ .

$Rd$	$Q$	$M$	$Nb$	$Nt$	$\frac{N u_{r}}{\sqrt{{Re}_{r}}}$
0.1					1.5322
0.2					1.5436
0.3					1.5525
0.4					1.5697
	0.1				1.2864
	0.2				1.3075
	0.3				1.3254
	0.4				1.3443
		0.1			1.6974
		0.2			1.7032
		0.3			1.7164
		0.4			1.7256
			0.1		1.6356
			0.2		1.6132
			0.3		1.5953
			0.4		1.5785
				0.1	1.1853
				0.2	1.2542
				0.3	1.3246
				0.4	1.3800

Figures 18 to 20 show the surface plot of $\sqrt{{Re}_{r}} C f_{r}$ , $\frac{N u_{r}}{\sqrt{{Re}_{r}}}$ , and $\frac{S h_{r}}{\sqrt{{Re}_{r}}}$ . It is revealed that the enhancement against the rising values of $α$ and $γ$ has been elaborated in Figure 18. Figure 19 shows that when $Q$ and $Rd$ increases results enhanced with $\frac{N u_{r}}{\sqrt{{Re}_{r}}}$ . $\frac{S h_{r}}{\sqrt{{Re}_{r}}}$ grows rapidly with the larger values of $K_{r}$ and $δ$ as seen in Figure 20.

Figure 18.

Impact of $α, γ$ on ${Re}_{r}^{1 / 2} C f_{r}$ .

Figure 19.

Impact of $Q, Rd$ on ${Re}_{r}^{- 1 / 2} N u_{r}$ .

Figure 20.

Impact of $δ, K_{r}$ on ${Re}_{r}^{- 1 / 2} S h_{r}$ .

Table 2 reveals the role of $γ$ , $α$ , and M on $\sqrt{{Re}_{r}} C f_{r}$ . It is noted that the greater values of $γ$ , $α$ , and M enhances the surface drag. Table 3 shows the effects of $Rd$ , $Q$ , $M$ , $Nb$ , and $Nt$ on $\frac{N u_{r}}{\sqrt{{Re}_{r}}}$ . From this Table, we have observed that greater values of $Rd$ , $Q$ , $M$ , and $Nt$ enhances the rate of heat transfer at the sheet surface while the greater values of $Nb$ reduces it. Table 4 shows the effects of $Nb$ and $Nt$ on $\frac{S h_{r}}{\sqrt{{Re}_{r}}}$ . From this Table, we observed that the higher values of $Nb$ increases the mass transfer rate and the higher values of $Nt$ reduce the mass transfer rate.

Table 4.

Effects of $Nb$ and $Nt$ on $\frac{S h_{r}}{\sqrt{{Re}_{r}}}$ .

$Nb$	$Nt$	$δ$	$K_{r}$	$E$	$\frac{S h_{r}}{\sqrt{{Re}_{r}}}$
0.1					1.6424
0.2					1.6633
0.3					1.6864
0.4					1.7035
	0.1				1.8675
	0.2				1.8589
	0.3				1.8064
	0.4				1.7532
		0.1			2.6543
		0.4			2.7456
		0.7			2.8012
		1.0			2.8489
			1.0		2.9235
			1.2		2.9476
			1.4		2.9645
			1.6		2.9865
				1.0	2.1043
				1.3	2.1224
				1.6	2.1654
				1.9	2.1834

Conclusions

In this effort, we have examined the laminar, time-independent, and incompressible flow of Maxwell nanofluid over a rotating disk surface. The problem is modeled in a cylindrical coordinate system. The model equations are formed as PDEs which then transformed into ODEs by using suitable similarity variables. The modeled equations are solved numerically by using bvp4c MATLAB built-in function. Furthermore, the effects of thermal radiation, heat source, thermophoresis, Brownian motion, chemical reaction, and Activation energy are taken into consideration. The final remarks are listed as:

• It is found that the greater ratio factor has enhanced the velocity profile in radial direction while decreased the velocity profiles in angular direction.

• The greater Maxwell factor has reduced the velocity profile in radial and angular directions.

• The greater magnetic factor has reduced the velocity profile in radial and angular directions.

• The greater magnetic, Brownian motion, thermophoresis, heat source, and thermal radiation factors have increased the thermal distributions.

• The greater Schmidt number, chemical reaction, and Brownian motion factors have reduced the concentration profiles while the greater thermophoresis factor and activation energy factor have increased the concentration profiles.

• The heat transfer rate is increased by the thermal radiation, heat source, magnetic, Brownian motion, and thermophoresis factors.

• The higher Brownian motion factor has increased the mass transfer rate and the higher thermophoresis factor has reduced the mass transfer rate.

Footnotes

Appendix

Acknowledgements

The authors extend their appreciation to University of Tabuk, Saudi Arabia.

Handling Editor: Aarthy Esakkiappan

Declaration of conflicting interests

The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Funding

The author(s) received no financial support for the research, authorship, and/or publication of this article.

ORCID iDs

Anwar Saeed

Arshad Khan

Data availability

All data used in this manuscript have been presented within the article.

References

Salahuddin

Mahmood

Khan

, et al. A permeable squeezed flow analysis of Maxwell fluid near a sensor surface with radiation and chemical reaction. Chem Phys 2022; 562: 111627.

Abdelsalam

Abbas

Megahed

, et al. A comparative study on the rheological properties of upper convected Maxwell fluid along a permeable stretched sheet. Heliyon 2023; 9: e22740.

Yasin

Hina

Naz

A modern study on peristaltically induced flow of Maxwell fluid considering modified Darcy’s law and Hall effect with slip condition. Alex Eng J 2023; 76: 835–850.

Muhammad

Ahmad

Farooq

, et al. Computational investigation of magnetohydrodynamics boundary of Maxwell fluid across nanoparticle-filled sheet. Al-Salam J Eng Technol 2023; 2: 88–97.

Khan

Shah

Khan

, et al. Numerical investigation of MHD Cattaneo–Christov thermal flux frame work for Maxwell fluid flow over a steady extending surface with thermal generation in a porous medium. Int J Thermofluids 2023; 20: 100418.

Abbas

Shatanawi

Shatnawi

TAM

. Transportation of nanomaterial Maxwell fluid flow with thermal slip under the effect of Soret–Dufour and second-order slips: nonlinear stretching. Sci Rep 2023; 13: 2182.

Faizan

Zaib

Khan

, et al. Modeling of thermal and solute transport within a Maxwell fluid in contact with a porous rotating disc. Eur Phys J Spec Top 2024; 233: 1519–1537.

Reddy

Meenakumari

Sucharitha

, et al. Heat and mass transfer analysis of conducting non-Newtonian nanofluid flows over an elongating sheet with a non-uniform heat source. Mod Phys Lett B 2024; 2450349.

Jalili

Sharif Mousavi

Jalili

, et al. Thermal evaluation of MHD Jeffrey fluid flow in the presence of a heat source and chemical reaction. Int J Mod Phys B 2024; 38: 2450113.

10.

Lone

Khan

Raiza

, et al. A semi-analytical solution of the magnetohydrodynamic blood-based ternary hybrid nanofluid flow over a convectively heated bidirectional stretching surface under velocity slip conditions. AIP Adv 2024; 14: 045013.

11.

Reddy

Shankar Goud

Nisar

, et al. Heat absorption/generation effect on MHD heat transfer fluid flow along a stretching cylinder with a porous medium. Alex Eng J 2023; 64: 659–666.

12.

Bilal

Tariq

Urva

, et al. A novel nonlinear diffusion model of magneto-micropolar fluid comprising Joule heating and velocity slip effects. Waves Random Complex Media. Epub ahead of print 30 May 2022. DOI: 10.1080/17455030.2022.2079761.

13.

Dang

Zhang

Yan

, et al. A new tool for understanding the solar wind–venus interaction: three-dimensional multifluid MHD model. Astrophys J 2023; 945: 91.

14.

Yang

Wang

, et al. Performance analysis of thermoacoustic plasma MHD generation. Energy 2023; 263: 125647.

15.

Shoukat

Zubair

Farman

, et al. Impacts of joule heating with Cattaeno-Christove heat flux model in a MHD flow of Erying-Powell fluid on a Riga plate. Alex Eng 2023; 64: 741–748.

16.

Khan

Shah

Alzahrani

, et al. Entropy generation and thermal analysis for rotary motion of hydromagnetic Casson nanofluid past a rotating cylinder with Joule heating effect. Int Commun Heat Mass Transf 2020; 119: 104979.

17.

Faizan

Ali

, et al. Modelling and analysis of heat transfer in MHD stagnation point flow of Maxwell nanofluid over a porous rotating disk. Alex Eng J 2024; 91: 237–248.

18.

Reddy

Ajithkumar

Zafar

SSB

, et al. Magnetohydrodynamic stagnation point flow of Williamson hybrid nanofluid via stretching sheet in a porous medium with heat source and chemical reaction. Proc IMechE, Part E: J Process Mechanical Engineering. Epub ahead of print 24 March 2024. DOI: 10.1177/09544089241239583.

19.

Swain

Mohammed Ibrahim

Dharmaiah

, et al. Numerical study of nanoparticles aggregation on radiative 3D flow of maxwell fluid over a permeable stretching surface with thermal radiation and heat source/sink. Results Eng 2023; 19: 101208.

20.

Manvi

Tawade

Biradar

, et al. The effects of MHD radiating and non-uniform heat source/sink with heating on the momentum and heat transfer of Eyring-Powell fluid over a stretching. Results Eng 2022; 14: 100435.

21.

Gautam

Rajput

Bhattacharyya

, et al. Comparative study of two non-Newtonian fluids with bioconvective induced MHD flow in presence of multiple slips, heat source/sink and nonlinear thermal radiation. Chem Eng J Adv 2022; 12: 100365.

22.

Ahmed

Zaib

Ali

, et al. Cattaneo–Christov Double Diffusion (CCDD) on Sutterby nanofluid with irreversibility analysis and motile microbes due to a RIGA plate. Micromachines 2022; 13: 1497.

23.

Gülmüş

Muratçobanoğlu

Mandev

, et al. Experimental and numerical investigation of flow and thermal characteristics of aluminum block exchanger using surface-modified and recycled nanofluids. Int J Numer Methods Heat Fluid Flow 2023; 33: 2685–2709.

24.

Sajid

Jamshed

Eid

, et al. Magnetized Cross tetra hybrid nanofluid passed a stenosed artery with nonuniform heat source (sink) and thermal radiation: novel tetra hybrid Tiwari and Das nanofluid model. J Magn Magn Mater 2023; 569: 170443.

25.

Pasha

Hussain

Alam

, et al. Impact of magnetized non-linear radiative flow on 3D chemically reactive sutterby nanofluid capturing heat sink/source aspects. Case Stud Therm Eng 2023; 41: 102610.

26.

Sajid

Pasha

Jamshed

, et al. Radiative and porosity effects of trihybrid Casson nanofluids with Bödewadt flow and inconstant heat source by Yamada-Ota and Xue models. Alex Eng J 2023; 66: 457–473.

27.

Kumar

Sharma

Entropy optimization analysis of Marangoni convective flow over a rotating disk moving vertically with an inclined magnetic field and nonuniform heat source. Heat Transf 2023; 52: 1778–1805.

28.

Sharma

Khanduri

Mishra

, et al. Combined effect of thermophoresis and Brownian motion on MHD mixed convective flow over an inclined stretching surface with radiation and chemical reaction. Int J Mod Phys B 2023; 37: 2350095.

29.

Waqas

Khan

Ali

, et al. Numerical computation of Brownian motion and thermophoresis effects on rotational micropolar nanomaterials with activation energy. Propuls Power Res 2023; 12: 397–409.

30.

Sudarmozhi

Iranian

Alessa

Investigation of melting heat effect on fluid flow with brownian motion/thermophoresis effects in the occurrence of energy on a stretching sheet. Alex Eng J 2024; 94: 366–376.

31.

Rashed

Mahmoud

Mabrouk

SM.

Enhanced flow and temperature profiles in ternary hybrid nanofluid with gyrotactic microorganisms: a study on magnetic field, Brownian motion, and thermophoresis phenomena. J Appl Comput Mech 2024; 10: 597–609.

32.

Thabet

Khan

Abd-Alla

, et al. Thermal enhancement, thermophoretic diffusion, and Brownian motion impacts on MHD micropolar nanofluid over an inclined surface: numerical simulation. Numer Heat Transf Part A Appl. Epub ahead of print 3 November 2023. DOI: 10.1080/10407782.2023.2276319.

33.

Almeida

Gireesha

Venkatesh

Magnetohydrodynamic flow of a micropolar nanofluid in association with Brownian motion and thermophoresis: irreversibility analysis. Heat Transf 2023; 52: 2032–2055.

34.

Tawade

Guled

Noeiaghdam

, et al. Effects of thermophoresis and Brownian motion for thermal and chemically reacting Casson nanofluid flow over a linearly stretching sheet. Results Eng 2022; 15: 100448.

35.

Shahzad

Imran

Tahir

, et al. Brownian motion and thermophoretic diffusion impact on Darcy-Forchheimer flow of bioconvective micropolar nanofluid between double disks with Cattaneo-Christov heat flux. Alex Eng J 2023; 62: 1–15.

36.

Zafar

Zaib

Lone

, et al. Numerical analysis of Williamson nanofluid over lubricated surface due to microorganism with thermal radiation. J Therm Anal Calorim. Epub ahead of print 19 February 2024. DOI: 10.1007/s10973-023-12874-y.

37.

Faizan

Ali

Loganathan

, et al. Entropy analysis of sutterby nanofluid flow over a riga sheet with gyrotactic microorganisms and Cattaneo–Christov double diffusion. Mathematics 2022; 10: 3157.

38.

Asaigeethan

Vaithiyalingam

Loganathan

, et al. Numerical study of Carreau fuzzy nanofluid across a stretching cylinder using a modified version of Buongiorno’s nanofluid model. Alex Eng J 2024; 101: 318–329.

39.

Zafar

Zaib

Lone

, et al. Radiative heat in a Williamson fluid flow through a lubricated surface containing swimming microorganism. J Therm Anal Calorim. Epub ahead of print 29 January 2024. DOI: 10.1007/s10973-023-12861-3..

40.

Faizan

Ajithkumar

Reddy

, et al. A theoretical analysis of the ternary hybrid nano-fluid with Williamson fluid model. Ain Shams Eng J 2024; 15: 102839.

41.

Kumar

Dharmendar Reddy

Computational modelling of radiative Maxwell fluid flow over a stretching sheet containing nanoparticles with chemical reaction. J Indian Chem Soc 2023; 100: 100877.

42.

Nasir

Berrouk

Aamir

, et al. Significance of chemical reactions and entropy on Darcy-forchheimer flow of H₂O and C₂H₆O₂ convening magnetized nanoparticles. Int J Thermofluids 2023; 17: 100265.

43.

Raghunath

Alfaleh

, et al. Effects of activation energy and chemical reaction on unsteady MHD dissipative Darcy–Forchheimer squeezed flow of Casson fluid over horizontal channel. Sci Rep 2023; 13: 2666.

44.

Shah

SAGA

Hassan

Karamti

, et al. Effect of thermal radiation on convective heat transfer in MHD boundary layer Carreau fluid with chemical reaction. Sci Rep 2023; 13: 4117.

45.

Khan

Kumam

Khan

, et al. Chemically reactive nanofluid flow past a thin moving needle with viscous dissipation, magnetic effects and hall current. PLoS One 2021; 16: e0249264.

46.

Arain

Khan

, et al. Theoretical exploration of heat transport in a stagnant power-law fluid flow over a stretching spinning porous disk filled with homogeneous-heterogeneous chemical reactions. Case Stud Therm Eng 2023; 50: 103406.

47.

Razaq

Hayat

Khan

, et al. ATSS model based upon applications of Cattaneo-Christov thermal analysis for entropy optimized ternary nanomaterial flow with homogeneous-heterogeneous chemical reactions. Alex Eng J 2023; 79: 390–401.

48.

Mishra

Anwar

Kumam

, et al. Numerical investigation of chemically reacting jet flow of hybrid nanofluid under the significances of bio-active mixers and chemical reaction. Heliyon 2023; 9: e17678.

49.

Ahmed

Khan

Ahmad

Joule heating effects in thermally radiative swirling flow of Maxwell fluid over a porous rotating disk. Int J Thermophysics 2019; 40: 1–18.

50.

Hafeez

Khan

Ahmed

, et al. Flow of Oldroyd-B fluid over a rotating disk through porous medium with Soret–Dufour effects. Arab J Sci Eng 2020; 45: 5949–5957.

51.

Khan

Hafeez

Lone

, et al. Darcy flow of convective and radiative Maxwell nanofluid over a porous disk with the influence of activation energy. Heliyon 2023; 9: e18003.

52.

Sparrow

Gregg

JL.

Heat transfer from a rotating disk to fluids of any Prandtl number. J Heat Transf 1959; 81: 249–251.