On generalized Bödewadt flow of TiO 2 /water nanofluid over a permeable surface with temperature jump

Abstract

Titanium dioxide (TiO₂) nanomaterial has numerous applications in the fields of cosmetics, medicines, coatings, inks, plastics, food, and textiles. Therefore, the problem of heat and mass transport on Bödewadt flow of TiO₂/water nanofluid over a rotating disk subjected to wall suction is studied. The impact of chemical reaction with partial slip and temperature jump conditions are also considered. For the numerical solution to the problem, the similarity variables are added to transform the three-dimensional flow equations into a favorable set of ordinary differential equations. The impacts of shear stresses, rates of heat and mass transport, and cooling efficiency of nanofluid on the flow are investigated by employing a bvp4c routine in Matlab software. Additionally, the plots for two-dimensional streamlines are presented to visualize the impact of slip velocity and rotation. Through asymptotic analysis, it is found that the presence of similarity solutions for nanofluid over the disk can occur only if the disk is driven to a significant amount of suction. The skin friction factor grows by enhancing the nanoparticle volume fraction $ϕ$ with a slight reduction in heat and mass transport rates. The fluid temperature is reduced by augmenting $ϕ$ providing the cooling efficiency of TiO₂. The fluid concentration falls significantly when a chemical reaction occurs at a faster rate.

Keywords

Bödewadt flow nanofluid wall suction chemical reaction rotating disk

Introduction

A great deal of attention has been paid in the past decade to the traditional Von Karman flow and heat transfer, for which the motion is driven by a rotating disk. Its cousin problem constitutes one of the classical problems of fluid mechanics such that the motion is superimposed owing to the fluid rotating with a uniform angular velocity at a larger distance from a stationary disk. However, although this case also possesses theoretical as well as practical importance, for example in understanding the dynamics of tornadoes and hurricanes, and rotor–stator systems in turbines, less research has been conducted to understand its physical insight. Therefore, the current work is devoted to the three-dimensional revolving flow motion and heat-induced by a stationary disk, but unlike the classical case as considered in the past studies, a disk with a uniform stretching in the radial direction is taken into account here for the first time in the literature. Bödewadt boundary layer flow occurs due to a rotating flow over a stationary disk and it represents a full analytical solution of the Navier–Stokes equations since it was first theoretically investigated by Bödewadt.¹ For the example of the conventional rotating disk problem,² the fluid forced outwards by the centrifugal force is replaced by a fluid stream in the axial direction. On the other hand, a reverse effect is observed for the revolving flow over a stationary disk, so that the fluid drawn to the axis of rotation is swept upwards, a phenomenon as a result of the radial pressure gradient being balanced by the centrifugal force. In the book by Schlichting³ it was cited that “The secondary flow which accompanies rotation near a solid wall can be clearly observed in a teacup: after the rotation has been generated by vigorous stirring and again after the flow has been left to itself for a short while, the radial inward flow field near the bottom will be formed. Its existence can be inferred from the fact that tea leaves settle in a little heap near the center at the bottom.” Bödewadt flow and heat transfer over a stretching stationary disk were reported by Turkyilmazoglu.⁴ Thermally developed generalized Bödewadt flow containing nanoparticles over a rotating surface with slip condition was studied by Abbas et al.⁵ Mahyuddin et al.⁶ discussed the Bödewadt flow and heat transfer in nanofluid over a permeable and radially stretching disk. A widespread inspection of this category of flows with its several uses can be found through^7–13 and references therein.

Nanofluids are relatively new generation heat transfer fluids that exhibit much higher thermal conductivity in comparison to conventional coolants even at low particle concentrations. Due to this reason, they have importance in a number of industrial sectors including transportation, power generation, micro-manufacturing, thermal therapy for cancer treatment, chemical and metallurgical sectors, as well as heating, cooling, ventilation, and air-conditioning. To be more specific magnetic nanofluids also known as ferrofluids can serve to deliver drugs and radiation to cancer patients without destroying the health tissues which is a significant side effect of conventional cancer treatment.¹⁴ Nanofluids are pretty useful in tackling cooling problems in thermal systems. Their high thermal conductivity can also be utilized in engine oils, automatic transmission fluids, coolants, and lubricants. The concept of NF for the enhancement of thermal conductivity of heat transfer fluids was first proposed by Choi and Eastman.¹⁵ In the study of nanofluids, the models of Tiwari and Das¹⁶ and Buongiorno et al.¹⁷ were frequently applied by many researchers. Khan et al.¹⁸ formulated the problem of heat transfer of nanofluids past a stretching disk using the Tiwari and Das model and the Buongiorno model. They found that nanofluids improve the rate of skin friction and heat transfer. Mustafa¹⁹ investigated the impacts of MHD, partial slip, and temperature jump on the nanofluid flow with the transfer of heat and mass above a rotating disk. Further, he observed the influences of the Brownian motion and thermophoresis on the fluid flow. The peristaltic motion of a couple stress nanomaterials in a tapered channel was investigated by Rafiq and Abbas.²⁰ In this study, authors adopted lubrication approximation theory to simplify the normalized equations. Alkuhayli²¹ analyzed the heat transfer analysis of a hybrid nanofluid flow on a rotating disk considering thermal radiation effects. Some recent studies in this direction were reported in Refs.^22–25

In general, it is assumed that the flow velocities stick to the wall of the body. However, in practical applications of science and engineering, it is sometimes impossible to satisfy the no-slip condition exactly on the disk, particularly when the surface of the disk is filled with impurities or the disk itself is rough. On the other hand, if the amplitude of the impurities or the characteristic scale of the protuberances is small as compared with the boundary layer thickness on the disk, the no-slip condition may be approximated by a partial slip condition applied to the envelope of the roughnesses, in accordance with the partial slip condition for rough surfaces proposed by Navier.²⁶ In micro-scale fluid dynamics, because of the micro-scale dimensions, the fluid flow behavior belongs to the slip flow regime and greatly differs from the traditional flow.²⁷ For example, in the situation when the fluid is particulate such as emulsions, suspensions, foams, and polymer solutions,²⁸ the no-slip condition is inadequate. In such cases, the suitable boundary condition is the partial slip. In general, the slip coefficients may be different in different directions. Whenever a uniform roughness on the surface is considered the slip coefficients may be equal. Sparrow et al.²⁹ considered the flow of a Newtonian fluid due to the rotation of a porous-surfaced disk and for that purpose replaced the conventional no-slip boundary conditions at the disk surface with a set of linear slip-flow conditions. A substantial reduction in torque then occurred as a result of the surface slip. This problem was recently reconsidered by Miklavcic and Wang³⁰ who pointed out that the same slip flow boundary conditions as those used by Sparrow et al.²⁹ also could be used for slightly rarefied gases or flow over grooved surfaces. The importance of slip in stretching sheet problems was also investigated by Wang.³¹ The presence of velocity slip on the wall may cause temperature rise or temperature jump, which must be taken into account in real-life applications.³² Turkyilmazoglu³³ recently demonstrated for the stretching sheet problem that an increase in the magnitude of temperature slip greatly affects the temperature profiles, resulting in variations in the rate of heat transfer.

The combination of heat and mass transfer has an essential role in the process of many applications in different fields. It controls temperature and humidity in the agricultural lands. Also, it prevents crops from damage by freezing, helps in evaporation of wet surfaces, etc. It also helps to purify blood in the liver and kidneys. This combination also separates chemicals in the process of distillation. Because of its many applications in different areas, many scientists and mathematicians paid attention to a model of stretchable disks. A numerical analysis of the time-dependent flow of viscous fluid due to a stretchable rotating disk with heat and mass transfer was reported by Ibrahim.³⁴ Agarwal³⁵ investigated the heat and mass transfer in electrically conducting micropolar fluid flow between two stretchable disks. Numerical simulation of heat and mass transfer in magnetic nanofluid flow by a rotating disk with variable fluid properties was reported by Sharma et al.²⁴ Few recent significant research regarding heat and mass transfer demonstrated in Refs.^36–44

Due to the numerous uses of titanium dioxide (TiO₂) in cosmetics, medicines, coatings, inks, plastics, and textiles. The Bödewadt flow of nanofluid over a rotating disk with heat and mass transport phenomena is examined in this paper. The impacts of wall suction, chemical reaction, velocity slip, and thermal slip are also investigated. This type of study has never been conducted and is unavailable in the literature. The novelty of the problem is illustrated in Table 1. The investigation is made by applying the Tiwari and Das nanofluid model. In section 2, a physical model of the nanofluid is proposed. In section 3, the similarity transformation is applied to convert the governing equations into a set of coupled ordinary differential equations, which are then evaluated numerically with the help of a Matlab solver bvp4c. In section 4, the impacts of physical parameters on the fluid velocities, temperature, concentration, skin friction factor, rate of heat transfer, and rate of mass transfer are discussed theoretically and presented in tabular and graphical forms. The results of the study could help to better understand heat transfer in similar conditions and have applications in fields such as energy engineering and thermal management.

Table 1.

The novelty of the current work.

Ref. no.	Venkatachala and Nath⁹	Mustafa¹⁹	Rafiq et al.²⁰	Sharma et al.²⁴	Present
Rotating disk	Yes	Yes	Yes	Yes	Yes
Nanofluid	Yes	No	No	Yes	Yes
Suction	Yes	Yes	Yes	No	Yes
Thermal slip	No	No	Yes	Yes	Yes
Chemical reaction	No	No	No	No	Yes

Problem formulation

Consider the three-dimensional steady, incompressible, and axisymmetric flow above the permeable disk with slip boundary and temperature jump conditions. Titanium dioxide is chosen as the nanomaterial while drinking water is selected as the base fluid. A constant angular velocity is assumed for the rotation of fluid around the vertical axis. The flow system is described using $(r, θ, z)$ coordinates. The axial symmetry will cause the derivatives along $θ$ to be dropped. A physical example of a rotating disk is displayed in Figure 1. The concentration and uniform temperature at the surface of the disk are $C_{w}$ and $T_{w}$ , respectively while $T_{\infty}$ and $C_{\infty}$ are the temperature and concentration of the far field, respectively.

Figure 1.

Flow geometry.

The basic equations that govern the flow are given as³⁶:

Mass conservation equation:

\frac{\partial}{\partial r} (ru) + \frac{\partial}{\partial z} (rw) = 0 .

(1)

Momentum equation:

ρ_{nf} (\frac{\partial u}{\partial r} u - \frac{v^{2}}{r} + \frac{\partial u}{\partial z} w) + \frac{\partial p}{\partial r} = μ_{nf} (\frac{\partial^{2} u}{\partial z^{2}} + \frac{\partial u}{\partial r} \frac{1}{r} - \frac{u}{r^{2}} + \frac{\partial^{2} u}{\partial r^{2}}),

(2)

ρ_{nf} (\frac{\partial v}{\partial r} u + \frac{vu}{r} + \frac{\partial v}{\partial z} w) = μ_{nf} (\frac{\partial^{2} v}{\partial z^{2}} + \frac{\partial v}{\partial r} \frac{1}{r} - \frac{v}{r^{2}} + \frac{\partial^{2} v}{\partial r^{2}}),

(3)

ρ_{nf} (\frac{\partial w}{\partial r} u + \frac{\partial w}{\partial z} w) + \frac{\partial p}{\partial z} = μ_{nf} (\frac{\partial^{2} w}{\partial z^{2}} + \frac{\partial w}{\partial r} \frac{1}{r} + \frac{\partial^{2} w}{\partial r^{2}}) .

(4)

Energy equation:

\frac{\partial T}{\partial r} u + \frac{\partial T}{\partial z} w = α_{nf} (\frac{\partial^{2} T}{\partial z^{2}} + \frac{\partial T}{\partial r} \frac{1}{r} + \frac{\partial^{2} T}{\partial r^{2}}) .

(5)

Concentration equation:

\frac{\partial C}{\partial r} u + \frac{\partial C}{\partial z} w = D_{B} (\frac{\partial^{2} C}{\partial z^{2}} + \frac{\partial C}{\partial r} \frac{1}{r} + \frac{\partial^{2} C}{\partial r^{2}}) - k_{r} (C_{w} - C_{\infty}) .

(6)

The boundary conditions (BCs) with the partial slip and temperature jump at the wall are expressed as³³:

\begin{matrix} For z = 0 u = τ_{rz} β_{1}, v = τ_{θ z} β_{1}, w = - w_{0}, \\ T = \frac{\partial T}{\partial z} β_{2} + T_{w}, C = C_{w}, \\ For z \to \infty u \to 0, v \to r Ω, T \to T_{\infty}, C \to C_{\infty} . \end{matrix}

(7)

Where $v_{nf}$ and $ρ_{nf}$ represent the kinematic viscosity and density of the NF, respectively, $D_{B}$ denotes the mass diffusivity, the thermal diffusivity of NF is denoted by $α_{nf}$ , $β_{1}$ , and $β_{2}$ signify the coefficients of velocity slip and thermal slip, respectively, and the physical names of all other terms are displayed in the nomenclature. The mathematical expressions used for NF are expressed in Table 2 whereas Table 3 displays some physical characteristics of base fluid and nanomaterials.

Table 2.

Correlations of nanofluid.^22,24

Fluid property	Correlations
Density	$\frac{ρ_{nf}}{ρ_{f}} = \frac{ρ_{s}}{ρ_{f}} φ + 1 - φ$
Heat capacity	$\frac{{(ρ C_{p})}_{nf}}{{(ρ C_{p})}_{f}} = \frac{{(ρ C_{p})}_{s}}{{(ρ C_{p})}_{f}} φ + 1 - φ$
Dynamic viscosity	$\frac{μ_{nf}}{μ_{f}} = \frac{1}{{[- (φ - 1)]}^{2.5}}$
Kinematic viscosity	$v_{nf} = \frac{μ_{nf}}{ρ_{nf}}$
Thermal conductivity	$\frac{k_{nf}}{k_{f}} = \frac{((k_{s} + 2 k_{f}) / (k_{f} - k_{s})) - 2 φ}{((k_{s} + 2 k_{f}) / (k_{f} - k_{s})) + 2 φ}$
Thermal diffusivity	$α_{nf} = \frac{k_{nf}}{{(ρ C_{p})}_{nf}}$

Table 3.

Thermal and physical properties of NPs and base fluid.²¹

Properties	TiO2	Water
ρ (kg/m³)	783	997.12
C_p (J/kg K)	2090	4179
k (W/mK)	0.145	0.6129

Note that the centrifugal force is balanced by the radial pressure gradient outside the boundary layer and hence we have $\partial p / \partial r = ρ_{nf} r Ω^{2}$ .

Similarity transformations

Now the Von Karman transformations² are introduced with the dimensionless velocities $[H (η), F (η), G (η)]$ , temperature $θ (η)$ , and concentration $ψ (η)$ to transform the governing equations into a favorable set of ordinary differential equations as:

\begin{matrix} \begin{matrix} [u (r, z), v (r, z), w (r, z)] = [r Ω F (η), r Ω G (η), \sqrt{v_{f} Ω} H (η)], \\ T (r, z) = T_{\infty} + (T_{w} - T_{\infty}) θ (η), C (r, z) = C_{\infty} + (C_{w} - C_{\infty}) ψ (η), \\ with η = z \sqrt{\frac{Ω}{v_{f}}} . \end{matrix} \end{matrix}

(8)

By utilizing the above transformations, equations (1)–(7) now become

H' = - 2 F,

(9)

F ″ = \frac{α_{1}}{α_{2}} (1 + H F' + F^{2} - G^{2}),

(10)

G ″ = \frac{α_{1}}{α_{2}} (H G' + 2 F G),

(11)

θ ″ = \frac{α_{4}}{α_{3}} P_{r} H θ',

(12)

ψ ″ = S_{c} H ψ' + k_{c} S_{c} ψ,

(13)

with BCs:

\begin{matrix} For η = 0 : H = - A, F - λ F' = 0, G - λ G' = 0, \\ θ - ξ θ' = 1, ψ = 1 . \\ For η \to \infty : F \to 0, G \to 1, θ \to 0, ψ \to 0 . \end{matrix}

(14)

The dimensionless terms arising in equations (9)–(14) are defined in Table 4, while the constant ratios $α_{1}, α_{2}, α_{3}$ , and $α_{4}$ are given as

\begin{matrix} α_{1} = \frac{1}{{(1 - ϕ)}^{2.5}}, α_{2} = ϕ \frac{ρ_{s}}{ρ_{f}} + 1 - ϕ, \\ α_{3} = \frac{{(ρ C_{p})}_{s}}{{(ρ C_{p})}_{f}} ϕ + 1 - ϕ, α_{4} = \frac{(k_{s} + 2 k_{f}) / (k_{f} - k_{s}) - 2 φ}{(k_{s} + 2 k_{f}) / (k_{f} - k_{s}) + φ} . \end{matrix}

(15)

Table 4.

The dimensionless quantities with notations and definitions.

Dimensionless numbers/parameters	Notations	Definitions
Wall suction parameter	$A$	$\frac{w_{0}}{\sqrt{v_{f} Ω}}$
Velocity slip parameter	$λ$	$β_{1} ρ_{f} \sqrt{v_{f} Ω}$
Thermal slip parameter	$ξ$	$β_{2} \sqrt{\frac{Ω}{v_{f}}}$
Prandtl number	$P_{r}$	$\frac{ν_{f} {(ρ C_{p})}_{f}}{k_{f}}$
Schmidt number	$S_{c}$	$\frac{ν_{f}}{D_{B}}$
Chemical reaction parameter	$k_{c}$	$\frac{k_{r}}{Ω}$

The practical terms of significance for the considered problem are the skin friction factor $C_{f}$ , the local Nusselt number $Nu$ , and the local Sherwood number $Sh$ , defined in the following sub-sections.

Skin friction factor

The skin friction factor for this problem is defined by

C_{f} = \frac{\sqrt{τ_{r}^{2} + τ_{θ}^{2}}}{ρ_{f} {(r Ω)}^{2}},

(16)

where $τ_{r}$ and $τ_{θ}$ represent the radial and tangential wall stresses, respectively, and both of these are expressed as

τ_{r} = μ_{nf} {(\frac{\partial u}{\partial z})}_{z = 0}, τ_{θ} = μ_{nf} {(\frac{\partial v}{\partial z})}_{z = 0} .

(17)

Using transformations (8) and equation (18), and equation (17) yields

{Re}_{r}^{1 / 2} C_{f} = \frac{\sqrt{G' {(0)}^{2} + F' {(0)}^{2}}}{{(1 - ϕ)}^{2.5}},

(18)

here $R e_{r} = r^{2} Ω / v_{f}$ signifies the local Reynolds number.

Nusselt number

Nusselt number for the present problem is defined as

Nu = \frac{r q_{w}}{(T_{w} - T_{\infty}) k_{f}},

(19)

here $q_{w}$ signifies the heat flux from the disk and is expressed as

q_{w} = - k_{nf} {(\frac{\partial T}{\partial z})}_{z = 0} .

(20)

Utilizing transformation (8) with equation (20), equation (19) gives

{Re}_{r}^{- 1 / 2} Nu = - \frac{k_{nf}}{k_{f}} θ' (0) .

(21)

Sherwood number

Sherwood number for this study is defined below

Sh = \frac{r q_{m}}{D_{B} (C_{w} - C_{\infty})},

(22)

where $q_{m}$ indicates the mass flux which is given by

q_{m} = - D_{B} {(\frac{\partial C}{\partial z})}_{z = 0} .

(23)

With the use of transformation (8) and equation (23) in equation (22), one finds

{Re}_{r}^{- 1 / 2} Sh = - ψ' (0) .

(24)

Numerical technique

The profiles of velocities, temperature, and concentration of the nanofluid are evaluated by employing the bvp4c routine in Matlab software. This built-in package is very easy to implement and reliable. It comprises a finite difference code that uses the three-stage Lobatto-IIIa formula to produce a continuous solution. The solution derived from this technique is accurate up to the fourth order in the entire domain. For the evaluation of the solution, the thickness of the boundary layer $η_{\infty}$ and initial guess are selected depending on the pertinent parameters. Next, the bvp4c solver is implemented by converting equations (9)–(14) into a system of first-order equations by substituting:

\begin{matrix} y (1) = H, y (2) = F, y (3) = F', y (4) = G, y (5) = G', \\ y (6) = θ, y (7) = θ', y (8) = ψ, y (9) = ψ' . \end{matrix}

(25)

Meanwhile, the first-order equations are attained as:

y (1)' = - 2 y (2),

(26)

y (2)' = y (3),

(27)

y (3)' = \frac{α_{1}}{α_{2}} (1 + y (1) y (3) + y (2) y (2) - y (4) y (4)),

(28)

y (4)' = y (5),

(29)

y (5)' = \frac{α_{1}}{α_{2}} (y (1) y (5) + 2 y (2) y (4)),

(30)

y (6)' = y (7),

(31)

y (7)' = \frac{α_{4}}{α_{3}} P_{r} y (1) y (7),

(32)

y (8)' = y (9),

(33)

y (9)' = S_{c} y (1) y (9) + k_{c} S_{c} y (8),

(34)

with BCs:

\begin{matrix} ya (1) + A = 0, ya (2) - λ ya (3) = 0, ya (4) - λ ya (5) = 0, \\ ya (6) - ξ ya (7) - 1 = 0, ya (8) - 1 = 0, \\ yb (2) \to 0, yb (4) - 1 \to 0, yb (6) \to 0, yb (8) \to 0 . \end{matrix}

(35)

The system is given in equations (26)–(35) is inserted in the bvp4c solver to obtain the numerical solution. The computations are made with a tolerance of $10^{- 4}$ and illustrated in plots.

Numerical results and discussion

In this section, the influences of the physical parameters are discussed theoretically and illustrated graphically. The nanofluid flow is governed by the momentum, energy, and concentration equations with the velocity slip, thermal slip, and chemical reaction. The calculations are executed numerically by adopting a Matlab routine bvp4c. The code is validated by making a comparison with the published work of Rahman and Andersson.³⁶ Table 5 demonstrates that the present results of $- θ' (0)$ are in perfect agreement with those of Rahman and Andersson.³⁶ In this study, $P_{r} = 6.2$ is taken as unity in the entire results. The numerical values of ${Re}_{r}^{1 / 2} C_{f}$ , ${Re}_{r}^{- 1 / 2} Nu$ , and ${Re}_{r}^{- 1 / 2} Sh$ for different values of nanoparticle volume fraction without suction $(A = 0)$ and with suction $(A > 0)$ are depicted in Table 6. It is noticed that increasing $ϕ$ values result in growth in the skin friction factor and produce a slight reduction in the heat as well as mass transport rates. Furthermore, the heat transfer rate remains uniform when there is no suction at the surface.

Table 5.

Comparison of $- θ' (0)$ with the results of Rahman and Andersson³⁶ by reducing $ϕ = S_{c} = k_{c} = λ = ξ = 0$ .

$P_{r}$	$A = 1$		$A = 2$		$A = 3$
$P_{r}$	Rahman and Andersson³⁶	Present	Rahman and Andersson³⁶	Present	Rahman and Andersson³⁶	Present
0.5	0.1934	0.2021	0.9239	0.9239	1.4736	1.4736
0.7	0.2984	0.3034	1.3154	1.3154	2.0713	2.0713
1	0.4943	0.4960	1.9114	1.9114	2.9705	2.9705
2	1.4013	1.4013	3.9215	3.9215	5.9743	5.9743
7	6.7857	6.7857	13.9633	13.9633	20.9879	20.9879

Table 6.

The numerical values of ${Re}_{r}^{1 / 2} C_{f}$ , ${Re}_{r}^{- 1 / 2} Nu$ , and ${Re}_{r}^{- 1 / 2} Sh$ at $λ = 0.5$ , $ξ = 0.1$ , $S_{c} = 0.5$ , and $k_{c} = 1$ for ascending values of nanoparticle volume fraction with and without suction.

$A$	$ϕ$	${Re}_{r}^{1 / 2} C_{f}$	${Re}_{r}^{- 1 / 2} Nu$	${Re}_{r}^{- 1 / 2} Sh$
0	0	0.8755	0	0.6078
	0.01	0.8898	0	0.6074
	0.02	0.9044	0	0.6070
	0.03	0.9193	0	0.6066
	0.04	0.9347	0	0.6062
	0.05	0.9502	0	0.6059
1	0	0.9557	3.6639	0.9297
	0.01	0.9710	3.6310	0.9284
	0.02	0.9867	3.5983	0.9271
	0.03	1.0027	3.5658	0.9259
	0.04	1.0190	3.5335	0.9246
	0.05	1.0358	3.5014	0.9234
2	0	1.0851	5.5173	1.3316
	0.01	1.1015	5.4607	1.3305
	0.02	1.1183	5.4045	1.3293
	0.03	1.1355	5.3487	1.3281
	0.04	1.1530	5.2934	1.3269
	0.05	1.1709	5.2384	1.3257

Figure 2(a) to (d) represents the impact of $ϕ$ on the fluid axial velocity $H (η)$ , radial velocity $F (η)$ , tangential velocity $G (η)$ , temperature $θ (η)$ , and concentration $ψ (η)$ , respectively. The fluid is expected to flow in a vertical direction and, therefore, $H (η)$ is negative. Figure 2(a) yields that the axial velocity of the fluid initially grows up to a maximum as the distance from the disk grows and when $z \to \infty$ , it asymptotically reaches a constant value. Furthermore, the velocity of the fluid in the axial direction upsurges by augmenting the values of $ϕ$ . Figure 2(b) depicts that the flow is radially outward close to the disk and appears to be radially inward at a specific distance from the disk. Moreover, the fluid radial velocity shows a downfall when the values of $ϕ$ are enhanced. It is observed in Figure 2(c) that the fluid tangential velocity is inversely proportional to the volume fraction of nanoparticles. The fluid radial and tangential velocity declined due to the collision of NPs, which reduces the velocity boundary layer (BL) thickness. In Figure 2(d), it is noticed that the temperature of the fluid is reduced on the ascending values of $ϕ$ . This is due to the reason that the thermal conductivity of TiO₂ nanomaterials is very low which causes the fluid temperature to decline.

Figure 2.

Impact of $ϕ$ on (a) $H (η)$ , (b) $F (η)$ , (c) $G (η)$ , and (d) $θ (η)$ .

Figure 3(a) to (e) displays the distribution of fluid velocities, temperature, and concentration with and without the effects of wall suction. The viscous drag causes the axial velocity of the fluid to decline close to the disk while the axial velocity remains constant at a distance far from the disk. For the conservation principle to be valid, the fluid particles close to the disk are forced to travel in the direction of the rotational axis due to the axially independent radial pressure gradient. As a result, in the Bödewadt flow, the axial flow is directed upward. The components of velocity exhibit oscillatory behavior and when suction is applied, these oscillations are damped and have significantly smaller amplitude. With a moderate wall suction velocity, such oscillatory behavior in the flow fields is suppressed. The asymptotic value of the velocity profiles is reached at a closer distance above the disk as a result of this impact. Therefore, it can be claimed that the wall suction’s damping effect is likely to dominate any absolute or convective type instabilities that are known from published studies. On higher values of the suction parameter, the radial and circumferential velocity components yield an increasing trend while the axial component of velocity reduces significantly with negative values revealing the vertically downward motion of the fluid. Figure 3(a) reveals that using higher $A$ values causes a greater volume of the fluid to be dragged toward the disk. Rafiq et al.¹² and Rafiq and Hashmi¹³ reported that the solution for temperature and concentration is possible only when $H (\infty)$ is negative. In Figure 3(d) it is noticed that $θ (η)$ has no meaningful solution for $A = 0$ and $A = 0.5$ , when $H (η) > 0$ (Figure 3(a)) while a physically acceptable solution for $θ (η)$ is achieved for $A \geq 1$ , when $H (η) < 0$ (Figure 3(a)). Moreover, the increased wall suction velocity suppresses the temperature and concentration of the fluid with the improvement in the heat and mass transmission from the disk as illustrated in Figure 3(d) and (e).

Figure 3.

Impact of $A$ on (a) $H (η)$ , (b) $F (η)$ , (c) $G (η)$ , (d) $θ (η)$ , and (e) $ψ (η)$ .

Figure 4(a) to (c) represents the fluctuations of velocity profiles with the impacts of $λ$ . In the presence of slip, it is possible to see the oscillatory nature of the three velocity components that Bödewadt¹ observed. As the slip effect becomes stronger, these oscillations are damped and their amplitude is further suppressed. Figure 4(a) illustrates that the radial velocity near the disk is reduced by augmenting the values of $λ$ but enhanced at a distance far from the disk. Moreover, ascending values of $λ$ cause the fluid tangential velocity to grow up and the fluid axial velocity to reduce. Also, note that the axial flow is significantly affected by the slip at a distance far from the disk.

Figure 4.

Impact of $λ$ on (a) $H (η)$ , (b) $F (η)$ , and (c) $G (η)$ .

Figure 5(a) portrays a declining trend in thermal boundary layer thickness as well as surface temperature by enhancing the values of the thermal slip parameter. It is because as the thermal slip grows then there is less heat transmission from the disk to the surrounding fluid layers, due to which the fluid temperature is dropped. Figure 5(b) depicts that the profiles of concentration exhibit a decreasing trend to large values of $k_{c}$ and $S_{c}$ . Since $S_{c}$ is inversely proportional to $D_{B}$ , therefore, greater $S_{c}$ values cause a reduction in fluid mass diffusivity and as a result, the concentration of the fluid drops. Moreover, higher values of the chemical reaction parameter produce a more destructive chemical reaction, which diffuses the liquid molecules, and hence the fluid concentration is reduced.

Figure 5.

Impact of (a) $ξ on θ (η)$ and (b) $S_{c}$ and $k_{c}$ on $ψ (η)$ .

The fluctuations of the volumetric flow rate $- H (\infty)$ , radial wall stress $F' (0)$ , tangential wall stress $G' (0)$ , and heat flux $- θ' (0)$ with changing values of $ϕ$ and $A$ are demonstrated in Figure 6(a) to (d). It is noticed that all of these profiles are augmented by enhancing the values of $A$ . Moreover, the ascending values of $ϕ$ result in a slight reduction in the volumetric flow rate as well as the radial and tangential wall stress but the radial wall stress slightly grows in the absence of wall suction. Note that, when $ϕ$ grows, the heat flux is slightly enhanced but it remains constant in the absence of wall suction as demonstrated in Figure 6(d).

Figure 6.

Impact of $A$ and $ϕ$ on (a) $- H (\infty)$ , (b) $F' (0)$ , (c) $G' (0)$ , and (d) $- θ' (0)$ .

The two-dimensional streamlines in the xz-plane are depicted in Figure 7(a) and (b) to visualize the effects of no-slip velocity and slip velocity. The streamlines with the no-slip condition are portrayed in Figure 7(a) whereas the streamlines with the impacts of velocity slip are displayed in Figure 7(b). The existence of slip can be seen in the figure which retards the flow of fluid. Figure 8 portrays the two-dimensional streamlines in the xy-plane. In this figure, the rotation effects are evident from the streamlines.

Figure 7.

Two-dimensional streamlines in xz-plane for $y = 0$ m at (a) $λ = 0$ and (b) $λ = 1.0$ .

Figure 8.

Two-dimensional streamlines in xy-plane for $z = 0$ m.

Conclusions

The focus of this study is on analyzing the heat and mass transfer phenomena that occur during the revolving Bödewadt flow of a nanofluid generated via a rotating disk. The novelty of the work lies in the investigation of heat and mass transfer in a nanofluid that flows due to a rotating disk, taking into account the effects of chemical reaction, partial slip, and temperature jump conditions. The study uses similarity transformations and the bvp4c method to solve the governing equations and presents graphical illustrations to aid in physical analysis. The motivation for this study arises from the widespread utility of nanofluids and the need for precise thermophysical and chemical properties to meet certain requirements necessary for several applications. The study aims to fill the knowledge gap by investigating the heat and mass transfer phenomena occurring during the flow of a nanofluid generated via a rotating disk, which has not been extensively analyzed in the literature. The validation of results and graphical depictions are essentially provided to test the accuracy of the solution and to facilitate physical analysis. These findings could be useful for understanding heat transfer in similar conditions in fields like energy engineering and thermal management. The key findings of this study are summarized below:

The skin friction factor increases with higher $ϕ$ values, which also causes a slight decrease in heat and mass transmission rates and when there is no suction at the surface, the rate of heat transfer remains constant.

The skin friction factor, the local Nusselt number, and the local Sherwood number all rise as the suction grows.

The amplitude of oscillatory profiles is inversely related to the axial distance, wall suction, and velocity slip.

The radial and circumferential velocity is directly related to the wall suction while the axial velocity is inversely related to it.

When the values of velocity slip are augmented, then the fluid tangential velocity upsurges, the axial velocity diminishes, and a significant impact is seen on the axial flow.

There is no boundary layer solution for heat transfer problems in pure Bodewadt flow with moderate or zero suction. The flow becomes physically possible when strong suction is applied.

The fluid temperature and the thermal BL thickness are both reduced by a significant rise in $ϕ$ . TiO₂ nanomaterials are also used as coolants in industrial equipment nowadays to enhance heat performance.

The fluid concentration is reduced by augmenting the Schmidt number and it is also declined by enhancing the chemical reaction parameter.

Footnotes

Appendix

Notation

$D_{B}$	Mass diffusivity $(m^{2} s^{- 1})$	$A$	Wall suction parameter
$F, G, H$	Dimensionless velocity functions	Sh	Sherwood number
$C$	Fluid concentration $(mol L^{- 1})$	Greek symbols
$T$	Fluid temperature $(K)$	$Ω$	Angular velocity $(s^{- 1})$
$q_{w}, q_{m}$	Heat and mass flux $(W m^{- 2})$	$α_{1 - 4}$	Constant ratios
$w_{0}$	Mass flux coefficient	$ρ$	Density $(kg m^{- 3})$
$r, θ, z$	Cylindrical polar coordinates $(m)$	$θ$	Dimensionless temperature function
$C_{f}$	Skin friction coefficient	$μ$	Dynamic viscosity $(kg m^{- 1} s^{- 1})$
$C_{p}$	Specific heat capacity $(J k g^{- 1} K^{- 1})$	$ν$	Kinematic viscosity $(m^{2} s^{- 1})$
$C_{w}, C_{\infty}$	Surface and ambient concentration $(mol L^{- 1})$	$τ_{r}$	Radial stress $(kg m^{- 1} s^{- 2})$
$T_{w}, T_{\infty}$	Surface and ambient temperature $(K)$	$τ_{θ}$	Tangential stress $(kg m^{- 1} s^{- 2})$
$k$	Thermal conductivity $(W m^{- 1} K^{- 1})$	$ϕ$	The volume fraction of $Ti O_{2}$ nanoparticles
$k_{r}$	Chemical reaction rate $(s^{- 1})$	$α$	Thermal diffusivity $(m^{2} s^{- 1})$
$u, v, w$	Velocity components $(m s^{- 1})$	$β_{2}$	Thermal slip coefficient
Dimensionless numbers/parameters		$ξ$	Thermal slip parameter
$k_{c}$	Chemical reaction parameter	$β_{1}$	Velocity slip coefficient
Nu	Nusselt number	$λ$	Velocity slip parameter
$P_{r}$	Prandtl number	Subscripts
$R e_{r}$	Reynolds number	$f$	Fluid
$S_{c}$	Schmidt number	$nf$	Nanofluid

Acknowledgements

We are thankful to the reviewers for their encouraging comments and constructive suggestions to improve the quality of the manuscript.

Handling Editor: Chenhui Liang

Declaration of conflicting interests

The author(s) received no financial support for 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

Zaheer Abbas

Shahana Siddique

Muhammad Yousuf Rafiq

Aqeel U Rehman

References

Bödewadt

. Die drehströmung über festem grunde. ZAMM-J Appl Math Mech/Zeitschrift für Angewandte Mathematik und Mechanik 1940; 20: 241–253.

Von Kármán

. Uber laminare und turbulente Reibung. Z Angew Math Mech 1921; 1: 233–252.

Schlichting

Gersten

. Boundary-layer theory. Berlin, Heidelberg, Germany: Springer, 2016.

Turkyilmazoglu

. Bödewadt flow and heat transfer over a stretching stationary disk. Int J Mech Sci 2015; 90: 246–250.

Abbas

Rafiq

Hasnain

, et al. Thermally developed generalized Bödewadt flow containing nanoparticles over a rotating surface with slip condition. Int Com-mun Heat Mass Transf 2021; 122: 105143.

Mahyuddin

Lok

Ahmad

. Bödewadt flow and heat transfer in nanofluid over a permeable and radially stretching disk. Sains Malays 2022; 51: 619–632.

King

Lewellen

. Boundary-layer similarity solutions for rotating flows with and without magnetic interaction. Phys Fluids 1964; 7: 1674–1680.

Kuo

. Axisymmetric flows in the boundary layer of a maintained vortex. J Atmos Sci 1971; 28: 20–41.

Venkatachala

Nath

. Mass transfer effects on the generalised vortex flow over a stationary surface with or without magnetic field. Proc Indian Acad Sci Sec C: Eng Sci 1980; 3: 211–223.

10.

Mustafa

Khan

Hayat

, et al. On Bödewadt flow and heat transfer of nanofluids over a stretching stationary disk. J Mol Liq 2015; 211: 119–125.

11.

Rahman

Andersson

. Heat transfer in generalized vortex flow over a permeable surface. Int J Heat Mass Transf 2018; 118: 1090–1097.

12.

Rafiq

Mustafa

Farooq

. Numerical assessment of Bödewadt flow and heat transfer over a permeable disk with variable fluid properties. Physica A: Stat Mech Appl 2019; 534: 122138.

13.

Rafiq

Hashmi

. Bödewadt flow over a permeable disk with homogeneous-heterogeneous reactions: a numerical study. Appl Sci 2019; 9: 4046.

14.

Wong

De Leon

. Applications of nanofluids: current and future. Adv Mech Eng 2010; 2: 519659.

15.

Choi

Eastman

. Enhancing thermal conductivity of fluids with nanoparticles (No. ANL/MSD/CP-84938; CONF-951135-29). Argonne, IL United States: Argonne National Lab, 1995.

16.

Tiwari

Das

. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int J Heat Mass Transf 2007; 50: 2002–2018.

17.

Buongiorno

Venerus

Prabhat

, et al. A benchmark study on the thermal conductivity of nanofluids. J Appl Phys 2009; 106: 094312.

18.

Khan

Mustafa

Hayat

, et al. Numerical study for Bödewadt flow of water based nanofluid over a deformable disk: buongiorno model. Indian J Phys 2017; 91: 527–533.

19.

Mustafa

. MHD nanofluid flow over a rotating disk with partial slip effects: Buongiorno model. Int J Heat Mass Transf 2017; 108: 1910–1916.

20.

Rafiq

Abbas

Ullah

. Peristaltic mechanism of couple stress nanomaterial in a tapered channel. Ain Shams Eng J 2022; 13: 101779.

21.

Alkuhayli

NAM

. Heat transfer analysis of a hybrid nanofluid flow on a rotating disk considering thermal radiation effects. Case Stud Therm Eng 2023; 49: 103131.

22.

Abbas

Hussain

Rafiq

, et al. Oscillatory slip flow of Fe3O4 and Al2O3 nanoparticles in a vertical porous channel using Darcy’s law with thermal radiation. Heat Transfer 2020; 49: 3228–3245.

23.

Ibrahim

Gamachu

. Entropy generation in radiative magneto-hydrodynamic mixed convective flow of viscoelastic hybrid nanofluid over a spinning disk. Heliyon 2022; 8: e11854.

24.

Sharma

Vijay

Mabood

, et al. Numerical simulation of heat and mass transfer in magnetic nanofluid flow by a rotating disk with variable fluid properties. Int Com-mun Heat Mass Transfer 2022; 133: 105977.

25.

Sharma

Mishra

, et al. Darcy-Forchheimer hybrid nanofluid flow over the rotating Riga disk in the presence of chemical reaction: artificial neural network approach. Alex Eng J 2023; 76: 101–130.

26.

Navier

CLMH

. Memoir on the laws of the movement of fluids. Mem Acad R Sci Inst Fr 1822; 6: 389–440.

27.

Gad-el-Hak

. The fluid mechanics of microdevices—the Freeman scholar lecture. J Fluids Eng 1999; 121: 5–33.

28.

Yoshimura

Prud’homme

. Wall slip corrections for Couette and parallel disk viscometers. J Rheology 1988; 32: 53–67.

29.

Sparrow

Beavers

Hung

. Flow about a porous-surfaced rotating disk. Int J Heat Mass Transf 1971; 14: 993–996.

30.

Miklavčič

Wang

. The flow due to a rough rotating disk. Z Angew Math Phys 2004; 55: 235–246.

31.

Wang

. Analysis of viscous flow due to a stretching sheet with surface slip and suction. Nonlinear Anal Real World Appl 2009; 10: 375–380.

32.

Hong

Asako

. Some considerations on thermal boundary condition of slip flow. Int J Heat Mass Transf 2010; 53: 3075–3079.

33.

Turkyilmazoglu

. Analytic heat and mass transfer of the mixed hydrodynamic/thermal slip MHD viscous flow over a stretching sheet. Int J Mech Sci 2011; 53: 886–896.

34.

Ibrahim

. Numerical analysis of time-dependent flow of viscous fluid due to a stretchable rotating disk with heat and mass transfer. Results in Phys 2020; 18: 103242.

35.

Agarwal

. Heat and mass transfer in electrically conducting micropolar fluid flow between two stretchable disks. Materials Today: Proc 2021; 46: 10227–10238.

36.

Rahman

Andersson

. On heat transfer in Bö dewadt flow. Int J Heat Mass Transf 2017; 112: 1057–1061.

37.

Arqub

. Numerical solutions for the Robin time-fractional partial differential equations of heat and fluid flows based on the reproducing kernel algorithm. Int J Numer Meth Heat Fluid Flow 2018; 28: 828–856.

38.

Arqub

Shawagfeh

. Application of reproducing kernel algorithm for solving Dirichlet time-fractional diffusion-Gordon types equations in porous media. J Porous Media 2019; 22: 411–434.

39.

Abu Arqub

. Numerical simulation of time-fractional partial differential equations arising in fluid flows via reproducing Kernel method. Int J Numer Meth Heat Fluid Flow 2020; 30: 4711–4733.

40.

Arqub

Al-Smadi

. Numerical solutions of Riesz fractional diffusion and advection-dispersion equations in porous media using iterative reproducing kernel algorithm. J Porous Media 2020; 23: 783–804.

41.

Ahmed

Khan

Ahmad

. Radiative heat flux effect in flow of Maxwell nanofluid over a spiraling disk with chemically reaction. Phys A: Stat Mech Appl 2020; 551: 123948.

42.

Sadiq

Hayat

. Entropy optimized flow of Reiner-Rivlin nanofluid with chemical reaction subject to stretchable rotating disk. Alex Eng J 2022; 61: 3501–3510.

43.

Abbas

Rafiq

Khaliq

, et al. Dynamics of the thermally radiative and chemically reactive flow of Sisko fluid in a tapered channel. Adv Mech Eng 2022; 14: 16878132221129735.

44.

Abbas

Shakeel

Rafiq

, et al. Rheology of peristaltic flow in couple stress fluid in an inclined tube: heat and mass transfer analysis. Adv Mech Eng 2022; 14: 16878132221139984.