Numerical analysis of Casson nanofluid three-dimensional flow over a rotating frame exposed to a prescribed heat flux with viscous heating.

This study investigates heat transfer characteristics and three-dimensional flow of non-Newtonian Casson nanofluid over a linearly stretching flat surface in the rotating frame of a reference. The current model includes the Buongiorno nanofluid model comprises nanoparticles' haphazard motion and thermo-migration. It also considered mechanisms for viscous heating and constant heat flux at the boundary. The nonlinear partial differential system modeling includes the non-Newtonian Casson fluid model and the boundary layer approximation. The system governing equations were nondimensionalized and numerically solved. A parametric study was conducted to analyze the significance of dimensionless parameters on velocities, the concentration, temperatures, Nusselt number, friction factors, and Sherwood number. The study reveals that the Casson nanoliquid temperature enhanced significantly due to the mechanisms of haphazard motion and thermo-migration. The momentum layer thickness of nano Casson fluid reduced due to the rotation phenomenon while the thermal layer structure amended notably. In the absence of rotation, there is no transverse velocity. The thermal layer structure is enhanced owing to the viscous heating process. The intense haphazard motion and thermo-migration mechanisms lead to maximum heat transfer rate at the plate. In addition, results show that the Coriolis force strength elevation shows similar axial and transverse velocities behavior. In addition, the nanoparticle concentration is observed higher due to the rotation aspect and Casson fluid parameter. Furthermore, the Casson fluid factor decreases with velocities, but the trend is the opposite for the high Casson fluid factor. The thermal and solute layer thickness growth is due to the nanoparticles' thermo-diffusion. In conclusion, the larger rotation factor increases the friction factors. The maximum plate heat transfer rate is when higher Nb and Nt are higher.

Introduction

Nowadays, the flow of non-Newtonian fluids plays a crucial role in several industrial applications such as the polymer industry, 3D printers, the production of plastics, paper manufacturing, etc., and numerous technological processes. The study of dynamics of non-Newtonian material is quite challenging to physicists, engineers, and mathematicians owing to the intricacy of these fluids. Numerous non-Newtonian material models are proposed, as a single constitutive equation does not display all possessions of non-Newtonian materials. The power-law material and two or three-grade fluid models are extensively studied in the literature (Andersson and Dandapat[1]; Hassanien[2]; Sadeghy and Sharifi[3]; Serdar and Salih Dokuz[4]; Sajid et al.[5,6]). The power-law model forecasts shear thickening and shear thinning performance. However, it is insufficient in articulating normal stress performance. A two-grade material model exhibit the normal stress effects. However, it is inadequate to express the performance of shear-thinning/thickening (see for more details Aksoy et al.[7]). The Maxwell model is a rate-type non-Newtonian material model that can envisage stress relaxation, but this model is incapable to study shear-dependent viscosity (see Hayat et al.[8]). The Casson material model is another type of non-Newtonian material model which is capable to predict the yield stress and Casson material is a shear-thinning material (see Fung[9]; Dash et al.[10]). In recent years, the Casson fluid model has been used extensively to study the dynamics of non-Newtonian materials owing to its wide range of industrial applications. Pramanik[11] studied the boundary layer transmission of a non-Newtonian Casson material and heat transfer over an exponentially elongating surface with mass suction or injection. He reported that enlarging values of the Casson material factor suppress the velocity layer thickness. Shaw et al.[12] examined the nonlinear convection transmission of Casson material toward a horizontal flat plate with Newton boundary conditions. They remarked that the friction factor at the plate was enhanced by the Casson material factor and buoyancy ratio factor. The Casson material model was utilized by Tamoor et al.[13] to investigate the heat transfer characteristics along with dissipation, magnetism, and Joule heating over a cylinder. They observed that the curvature factor has an analogous stimulus on the velocity, skin friction, temperature, and Nusselt number. The importance of adjustable viscosity and thermal conductivity is presented by Animasaun[14] in the study of the non-Darcy flow of dissipating Casson material with thermophoresis. Raju and Sandeep[15] considered the time-related three-dimensional dynamics of Casson–Carreau materials past an elongated surface. The influence of porous matrix on the dynamics of Casson material is analyzed by Nadeem et al.[16] by considering three-dimensional flow and magnetism aspects. Shehzad et al.[17] extended the problem of Nadeem et al.[16] by considering the internal heat production mechanism. Analysis of cross-diffusion and Rosseland thermal radiation aspects on the 3D dynamics of Casson material over a heated surface is performed by Zia et al.[18]. They remarked that the velocity layer structures were diminished by the Casson material factor. Prasad et al.[19] investigated the 3-directional dynamics of Casson material over a porous slender flat plate subjected to heat production or captivation. Very recently, Salahuddin et al.[20] examined the Casson material flow subjected to the activation energy, internal energy change, and thermophysical properties.

Recently, a vast number of industries merged nanomaterials with their broad applications for their superior properties and variety of implementation methodologies. Examples of such applications include thermal energy storage systems, electronic cooling, advanced nuclear thermal systems, hybrid-powered machines, solar liquid heating/cooling, and optoelectronics, etc. Nanofluids are designed by partying tiny solid particles (Al2O3, TiO2, ZnO2, Cu, Al, Au, etc.) and base fluid (such as ethylene glycol, kerosene, oil, and water). Conveying nanoparticles into the functioning fluids improves the thermal performance of functioning fluids, this fact was first remarked by Choi and Eastman[21]. More industrial applications of using nanofluids are presented Refs.[22-34]. Buongiorno[35] investigated nanoliquids convective heat transfer utilizing a mathematical model that involves two leading mechanisms Brownian motion and thermodiffusion. By utilizing the model proposed by Buongiorno, the boundary layer dynamics of nanofluid toward an upright flat plate were investigated by Kuznetsov and Nield[36]. They reported that the Brownian motion and thermodiffusion mechanisms are accountable for augmenting the temperature boundary layer structure. Minkowycz problem was revisited by Nield and Kuznetsov[37] by employing the Buongiorno model and considering porous matrix effects. Boundary layer transmission of nanomaterial over a horizontal plate is investigated by Khan and Pop[38]. They employed the Keller box method to treat the nonlinear problem and found that the mechanism of the Brownian movement of nanoparticles reduces the concentration of the nanoparticle’s layer structure.

The three-directional transmission of magneto-nano Casson material over a plate with Rosseland radiation and Newton boundary condition is investigated by Gupta and Sharma[39]. They utilized the Buongiorno nanofluid model and found that the temperature of Casson fluid enhances by conveying nanoparticles. Nadeem et al.[40] examined the three-dimensional dynamics of magneto-nanofluid over a plate that is elongated with linear velocity. They also scrutinized the influence of the Newton boundary condition. Saeed et al.[41] studied the thin-film dynamics of Casson nanofluid past an inclined spinning disk subjected to Rosseland heat and internal heat production. Mahanthesh et al.[42] used the Buongiorno model to investigate the dynamics of two-phase particulate nanofluids over a vertical plate. Naga Santoshi et al.[43] presented the numerical simulations of the 3-directional flow of Casson–Carreau nanomaterial. However, the published studies related to 3-directional flow Casson nanofluid flow using the Buongiorno model over a rotating stretchable surface are quite limited. In addition, a review paper that deals with the heating/cooling processes of suspended nanomaterials in lubricants and refrigerants are presented by Yang et al.[44].

With regards to the thermo-migration of nanoparticles in the motion of various fluids, all relevant published papers in the literature have reported the reason behind this is that the concentration decreases with thermophoresis after the formation of a relatively particle-free layer near the surface due to the migration of tiny particles through thermophoresis. The impact of thermophoresis shall be ignored when considering either the energy flux results from the composition gradient or the mass flux created by a temperature gradient. The presence of different responses to the force caused by a temperature gradient is adequate to augment the temperature distribution due to an increase in thermophoresis For Newtonian liquids, the impact of thermophoresis on the nanoparticle’s concentration decreases. In contrast, the impact increases for nonNewtonian liquids. The current paper utilizes a similar model as Owhaib et al.[45], using the bvp5c algorithm and Buongiorno nanoliquids model. They mathematically studied the boundary value problem of an elongated rotating plate with constant heat flux with thermal radiation. Among the physical observation, they conclude that the nanofluid single-phase model is inadequate to study their problem.

The main objective of the current investigation is to study the 3D rotating flow of Casson liquid over an elongated flat frame exposed to a constant heat flux condition with viscous heating. The nanoparticle's haphazard movement and thermo-migration are incorporated in the model utilizing the Buongiorno model. The model is solved using the MATLAB bvp5c algorithm. The obtained results are analyzed and presented in graphs to display the impact of various parameters such as Nb, Nt, Ec, Ca, Pr, Le, and Ro on velocity, temperature, and concentration profiles. In addition, the heat transfer rate and the friction coefficients are to be scrutinized and presented as well.

Mathematical formulation

This study considers the laminar, steady-state, three-directional flow of an elongated non-Newtonian Casson nanofluid surface over a rotating frame. The nanoliquids surface is exposed to a constant surface heat flux condition at the boundary. The Cartesian coordinate framework is aligned with -plane and the fluid region is considered at . The incompressible Casson nanoliquid rotates consistently about -axis with an unvarying rate . Figure 1 shows a schematic diagram of the problem under investigation.

Figure 1

Schematic diagram of the problem under consideration.

Deploying Prandtl’s boundary layer approach, the current problem under investigation governing equations are defined as below (note Refs.[35,46]):where -dimensionless Casson fluid factor, and are velocities along and ‐directions, is the kinematic viscosity, is the dynamic viscosity, is the density, is the temperature, is the nanoparticle concentration, is the thermal diffusivity, is the thermal conductivity, is the specific heat of the Casson fluid, is the specific heat of the nanoparticles, is the coefficient of Brownian diffusion, is the coefficient of thermo-migration diffusion and is the ambient temperature.

The pertinent boundary conditions arewhere a is stretching rate and is constant heat flux.

Upon considering

Equation (1) is satisfied and the other Eqs. (2)–(6) yieldwhere —similarity variable, and are dimensionless axial velocity, transverse velocity, temperature, and nanoparticle concentration correspondingly, (Prandtl number), (rotation parameter), (Lewis number), (thermophoresis parameter), (Brownian motion parameter), and (Eckert number).

The dimensionless expressions of friction factors, Nusselt number, and Sherwood number are given below;where is Reynolds number.

Numerical technique

The nonlinear boundary value problem defined in Eqs. (8)–(12) is solved numerically utilizing the bvp5c routine by letting , , and to obtain the following single-order differential system:with

The system first-order governing equations are solved using the MATLAB—bvp5c routine (see Ref.[47]). This routine integrates a system of differential equations of the form y′ = f(x,y), subject to the boundary conditions. It is using the finite difference method with achievable accuracy of 10–8. Condition at infinity is rescaled to 5. For method validation, the obtained numeric solution − θ'(0) is compared with the published studies at and . Table 1 data reveals that the current study results are in good agreement with the data in the literature. In the next section, a parametric analysis is performed.

Table 1

Comparison of − θ'(0) values with of Khan and Pop[38] and Gorla and Sidawi[46] when and .

Pr_l	Khan and Pop[38]	Gorla and Sidawi[46]	Present (bvp5c)
0.07	0.0663	0.0656	0.06562
0.2	0.1691	0.1691	0.16909
0.7	0.4539	0.5349	0.45392
2	0.9113	0.9114	0.91136
7	1.8954	1.8905	1.89542

Interpretation of the outcomes

The Casson nanoliquid 3D rotating flow over a linearly stretched surface in the presence of viscous heating, Brownian movement, thermo-migration, and heat flux condition is studied. The study numeric solution presented in Figs. 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26 and 27 show the parametric analysis of the prime parameters such as rotation parameter (), Casson fluid factor (), thermo-migration parameter (), Eckert number (), Lewis number (), Brownian movement factor () and Prandtl factor () on velocities , temperature , nanoparticle concentration , friction coefficients at the plate , Sherwood number () and () Nusselt number distributions. Figures 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26 and 27 are plotted for , and .

Figure 2

Effect of on .

Figure 3

Effect of on .

Figure 4

Effect of on .

Figure 5

Effect of on .

Figure 6

Effect of on .

Figure 7

Effect of on .

Figure 8

Effect of on .

Figure 9

Effect of on .

Figure 10

Effect of on .

Figure 11

Effect of on .

Figure 12

Effect of on .

Figure 13

Effect of on .

Figure 14

Effect of on .

Figure 15

Effect of on .

Figure 16

Effect of on .

Figure 17

Effect of on .

Figure 18

Effect of on .

Figure 19

Effect of on .

Figure 20

Effect of and on .

Figure 21

Effect of and on .

Figure 22

Effect of and on .

Figure 23

Effect of and on .

Figure 24

Effect of and on .

Figure 25

Effect of and on .

Figure 26

Effect of and on .

Figure 27

Effect of and on .

The influence of rotation aspect () on velocities , temperature , and nanoparticle concentration is accessible in Figs. 2, 3, 4 and 5 correspondingly. The ratio of angular velocity rate and elongating velocity rate is termed as a rotation factor. Higher values of rotation factor designate to lesser elongating velocity rate, as a result, the magnitude of velocities are reduced by increasing values of . Figures 4 and 5 show that temperature and nanoparticle concentration increased with increasing rotation factor (). In addition, there is no transverse velocity in the absence of rotation factor () and the axial velocity layer thickness was found to be higher. In contrast, the thermal and nanoparticle concentration layer structure is suppressed in the absence of rotation factor (). Figures 6, 7, 8 and 9 designate how the Casson fluid factor () disturbs the velocities , temperature , and nanoparticle concentration fields. A decreasing tendency is detected for the axial velocity layer structure when values are increased (see Fig. 6), while this trend is opposite for the thermal and nanoparticles concentration layer structure (see Figs. 8, 9). In Fig. 7, the magnitude of the transverse velocity () decreases near the elongated plate, and then it increases for larger values of .

The importance of Eckert number () on temperature and nanoparticle concentration distributions are outlined in Figs. 10 and 11 correspondingly. Improving the tendency of temperature and nanoparticle concentration layer structure is found for progressive values of . The kinetic energy is directly related to the Eckert number (). As increases, the kinetic energy strength increases, consequently, the magnitude of temperature and nanoparticle concentration enhances for higher values of . The consequence of on temperature and nanoparticle concentration distributions are delineated in Figs. 12 and 13. By definition, higher corresponds to stronger thermal diffusivity and lower mass diffusivity. Stronger thermal diffusivity is accountable for thicker thermal layer structure (see Fig. 12), similarly, a lower mass diffusivity is accountable for thinner nanoparticle concentration layer structure (see Fig. 13).

To analyze the effects of nanoparticles on the dynamics of Casson fluid through Brownian movement () and thermo-diffusion of nanoparticles (), Figs. 14, 15, 16 and 17 are presented. The arbitrary movement of nanoparticles generates supplementary internal heat in the Casson fluid due to collision, and this is accountable for improvement in the thermal field (see Fig. 14). In Fig. 15, the nanoparticle concentration field diminished significantly for enlarging values of . On the other hand. The thermophoresis mechanism is related to thermal gradient and that produces additional heat in the Casson fluid system and thereby the magnitude of temperature and nanoparticles concentration fields increased (see Figs. 16, 17). Thus one can conclude that the thermal layer structure of Casson fluid improved notably by conveying nanoparticles. Figures 18 and 19 display the influence of Prandtl number () on temperature and nanoparticle concentration . A significant decrement occurred in via higher . Higher corresponds to the lower thermal conductivity of the Casson fluid, which causes a decrement in the magnitude of the temperature field (see Fig. 18). In Fig. 19, the nanoparticle concentration increases near the plate and decreases far away from the plate by increasing numeric values of .

Figures 20 and 21 visualize the impact of rotation factor () and Casson fluid factor () on the plate friction factors along and directions ( respectively when and . Here both and are upsurged for increasing values of . This is because the velocities layer structure is thinner at the plate for larger values of . While the and are diminished for cumulative values of (see Fig. 21). Figures 22 and 23 are drawn to examine the stimulus of and on and correspondingly. Here the heat transfer rate at the plate is found maximum when both and values are set higher (see Fig. 22). Whereas, the enhancing and declining tendency of is seen for larger values of and respectively (see Fig. 23). Figures 24 and 25 depict the influences of and on and correspondingly. Figure 24 exhibits that increases for larger and . In Fig. 25, the mass transfer rate at the plate is found maximum when is kept lower value and the value of set to higher. Figures 26 and 27 elucidate the performance of and for distinct values of and . Figure 26 exhibits that the increases for larger and . However, the is an increasing property of and decreasing property of .

Final remarks

The major findings of the current research work are given below:

An elevation in the strength of Coriolis force exhibits similar behavior for both axial and transverse velocities.

Stronger Coriolis force exerted by rotation improves the thermal layer structure.

The nanoparticle concentration is found higher due to the rotation aspect and Casson fluid parameter.

The viscous dissipation effect leads to an enhancement of nanoparticle concentration and temperature profiles.

There is no transverse velocity in the absence of rotation factor while the axial velocity layer thickness was found to be higher in the absence of rotation factor.

The Casson fluid factor exhibits decreasing behavior of velocities, while the opposite trend is observed for the larger Casson fluid factor.

The thermo-diffusion due to nanoparticles is accountable for thermal and solute layer thickness growth.

Temperature layer thickness enlarges owing to the nanoparticles’ random movement.

Temperature and nanoparticle’s concentration layer thickness improves for larger Eckert number.

The friction factor increased for the larger rotation factor.

The heat transfer rate at the plate is found maximum when both and values are set higher.