Next Article in Journal
Construction of N-CDs and Calcein-Based Ratiometric Fluorescent Sensor for Rapid Detection of Arginine and Acetaminophen
Next Article in Special Issue
Fe-Cr-Nb-B Ferrofluid for Biomedical Applications
Previous Article in Journal
Analysis of the Effect of Electrode Materials on the Sensitivity of Quartz Crystal Microbalance
Previous Article in Special Issue
Evaluation of Multiple Semi-Twisted Tape Inserts in a Heat Exchanger Pipe Using Al2O3 Nanofluid
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nanoparticle Aggregation and Thermophoretic Particle Deposition Process in the Flow of Micropolar Nanofluid over a Stretching Sheet

1
School of Mathematics and Statistics, Xuzhou University of Technology, Xuzhou 221018, China
2
Department of Mathematics, Nanchang Institute of Technology, Nanchang 330044, China
3
Department of Mathematics, Davangere University, Davangere 577002, India
4
Department of Mathematical Sciences, Faculty of Science and Technology, University Kebangsaan Malaysia, UKM, Bangi 43600, Malaysia
5
Department of Mathematics and Social Sciences, Sukkur IBA University, Sukkur 65200, Pakistan
6
Department of Mathematical Sciences, Federal Urdu University of Arts, Science & Technology, Gulshan-e-Iqbal Karachi 75300, Pakistan
7
Department of Physics, College of Sciences, University of Bisha, P.O. Box 344, Bisha 61922, Saudi Arabia
8
Physics Department, Faculty of Science, Al-Azhar University, Assiut 71524, Egypt
9
Laboratory of Nano-Smart Materials for Science and Technology (LNSMST), Department of Physics, Faculty of Science, King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
10
Research Center for Advanced Materials Science (RCAMS), King Khalid University, P.O. Box 9004, Abha 61413, Saudi Arabia
11
Nanoscience Laboratory for Environmental and Biomedical Applications (NLEBA), Metallurgical Laboratory 1, Department of Physics, Faculty of Education, Ain Shams University, Roxy, Cairo 11757, Egypt
12
Radiological Sciences Department, College of Applied Medical Sciences, King Khalid University, P.O. Box 9004, Abha 61421, Saudi Arabia
13
BioImaging Unit, Space Research Centre, Michael Atiyah Building, University of Leicester, Leicester LE1 7 RH, UK
14
Mechanical Engineering Department, College of Engineering, Prince Sattam Bin Abdulaziz University, Wadiaddawaser 11991, Saudi Arabia
15
Production Engineering and Mechanical Design Department, Faculty of Engineering, Mansoura University, Mansoura 35516, Egypt
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(6), 977; https://doi.org/10.3390/nano12060977
Submission received: 12 February 2022 / Revised: 3 March 2022 / Accepted: 14 March 2022 / Published: 16 March 2022
(This article belongs to the Special Issue Applications of Nanofluids – II)

Abstract

:
The purpose of this research is to investigate the consequence of thermophoretic particle deposition (TPD) on the movement of a TiO2/water-based micropolar nanoliquid surface in the existence of a porous medium, a heat source/sink, and bioconvection. Movement, temperature, and mass transfer measurements are also performed in the attendance and nonappearance of nanoparticle aggregation. The nonlinear partial differential equations are transformed into a system of ordinary differential equations using appropriate similarity factors, and numerical research is carried out using the Runge-Kutta-Felhberg 4th/5th order and shooting technique. The obtained results show that improved values of the porous constraint will decline the velocity profile. Improvement in heat source/sink parameter directly affects the temperature profile. Thermophoretic parameter, bioconvection Peclet number, and Lewis number decrease the concentration and bioconvection profiles. Increases in the heat source/sink constraint and solid volume fraction will advance the rate of thermal dispersion. Nanoparticle with aggregation exhibits less impact in case of velocity profile, but shows a greater impact on temperature, concentration, and bioconvection profiles.

1. Introduction

Non-Newtonian materials are frequently encountered in a plethora of technological applications, including the production of crude soft elastic materials, the behavior of lubricants, sludges, pigments, packaged foods, and the mobility of organic fluids. Rheological fluids include blood, sauces, solvent cosmetics, colors, some lubricants, and rupturing mud. The relationship between stress and strain within those fluids is nonlinear. These liquids are substantially more complex to study than Newtonian fluids. As a result, various nonlinear frameworks for non-Newtonian liquids have been proposed. Even Navier–Stokes formulas are complex, containing a lot of restrictions, and the findings of resulting equalities are more difficult to acquire [1]. Eringen [2] shows that the Navier–Stokes theory fails to precisely describe the properties of different liquids indicated by microscopic characteristics arising from the micromotions and local structure of liquid constituents. These fluids are used to investigate the properties of contaminated lubricants, oils, polymeric mixtures, nanocrystals, animal blood containing rigid cells, and a variety of other biological solutions. Based on the wide range of applications, many researchers are working on this fluid model recently are addressed here [3,4,5,6,7,8].
Nanofluids have a wide range of roles in industry and engineering. The use of nanoparticles in fluid flow analysis has shown several uses such as power generation, microscale electronics, chemical processes, etc. Depending on the type of nanoparticles used in the composition procedure, nanofluid is also known as nanomaterials. Improved heat conductivity, which is the core and fundamental property of nanomaterials, is crucial. The most often used base fluids in this regard are water, toluene, ethylene, and kerosene oil, and several related research have recently been published to actually apply the concept of nanofluids. Nanofluids have distinct properties that make them particularly useful in a variety of thermal control activities. They have better thermal conductivity as well as a better convective heat transfer coefficient than the base fluid. Bahiraei et al. [9,10] examined the heat transport and entropy analysis in the presence of nanofluid. Khan et al. [11] explored the unstable stagnation point motion of hybrid nanoparticles over a spinning disk. The impacts of a magnetic dipole in the attendance of the KKL model to stimulate radiative nano liquid flow across a stretched sheet via Kumar et al. [12] Alhadhrami et al. [13] studied the effect of TPD on nanoliquid Glauert wall jet slip flow. In the existence of chemically reactive activation energy, Madhukesh et al. [14] examined the Marangoni-Bio convection movement of Casson nano liquid across a porous media. Muhammad and Nadeem [15] discussed convective heat transfer in the presence of magnetic dipole by using nickelzinc ferrite, manganese zinc ferrite, and magnetite ferrite with ethylene glycol. Guo et al. [16] examined the heat transfer characteristics in a nanofluid containing γ-Fe2O3 nanoparticles. Xie et al. [17] studied thermal conductivity improvement in aluminum nanoparticles suspended in a base fluid. They found that thermal conductivity is improved with the suspension of nanoparticles. The researchers have always found it difficult to model the thermophysical parameters in such a way that they are consistent with the experimental findings. Using standard models, the gain in the thermal conductivity of nanoliquids was grossly overestimated. The relative viscosity of aggregated nanoparticles was modeled using the fractal approach [18]. Recently, Mackolil and Mahanthesh [19] investigated the effects of nanoparticle aggregation and temperature-dependent surface tension on Marangoni convection in TiO2–EG nanoliquid. Here are some of the most important papers on nanoparticle aggregation [20,21,22].
Thermal transportation and dynamics of fluid flow research have scientific and practical applications in cable design, plastic sheet engineering, fiberglass structure, and thermal spouting. The heat source/sink (HS/S) is another major factor that determines heat transmission. When a heat source or heat sink is used, the heat dispersal throughout the entire field changes. In this view, many researchers do significant works on HS/S. Recently, Saleh et al. [23] probed the magnetic dipole and HS/S effects on Maxwell hybrid nanofluid flow across a stretched sheet. Madhukesh et al. [24] inspected the behavior of water transporting SWCNT nanoparticles and swimming microorganisms across a Riga plate with an HS/S. Micropolar liquid moving through a convectively hot surface, containing an nth order chemical process and an HS/S, was studied by Sajid et al. [25]. Khan et al. [26] investigated the mixed convection movement of pair stress nanofluid with HS/S effects across an oscillating stretchable surface. Chu et al. [27] looked at mathematical simulations of time-dependent micro-rotation blood circulation induced by a curved moving surface via gold nanoparticle conduction with non-uniform HS/S.
When there is a temperature gradient, the suspended particle will tend to travel from high to low-temperature areas. The force that causes this phenomenon is known as the thermophoretic force. The phenomenon of thermophoresis is important in the mass transfer mechanism of numerous devices that use microscopic micron-sized particles and huge temperature differences in the fields. Many researchers have shown interest in this concept over various geometries. Shehzad et al. [28] examined the forced convective Maxwell fluid flow across a spinning disc while thermophoretic particles motion. Kumar et al. [29] investigated the effects of thermophoretic particle’s deposition (TPD) on heat and mass transfer in Casson fluid flow through a moving thin needle. Some of the works on the notion of TPD can be found in [30,31,32].
Microbes are monocellular organisms that are significantly more productive than plants at lessening the greenhouse gases and absorption of CO2 from the environment. Monocellular organisms, on the other hand, can only be employed in a mixture of nanofluids. The fertility of the soil is improved by using microorganisms. Microbes, such as shielded cells, are considered as gyrotactic microorganisms for moving to a liquid for increasing density stratification based on the gradient thickness. Gyrotactic microorganisms are added to nanofluids to promote nanoparticle preservation and mass propagation. Ali et al. [33] investigated the influence of combined convection and magnetohydrodynamic (MHD) flow on the dynamics of Casson nano liquid in the stagnation point of a revolving sphere with finite element modelling. Alqarni et al. [34] evaluated the impact of bioconvection in 3D viscoelastic nanoliquid movement due to exponentially stretched surface with nonlinear radiative heat transmission and variable thermal conductivity. Revised Fourier’s and Fick’s Laws for duplicating mixed bioconvective motion of radiative-reactive Walters-B liquid transmission of minute particles subject to Lorentz force on radiative-reactive Walters-B fluids were investigated by Wakif et al. [35]. Farooq et al. [36] investigated the thermally radioactive bioconvection current of Carreau nano liquid using modified Cattaneo–Christov expressions and an exponential space-based heat source.
According to the previous section affordable literature study, no study on micropolar nanoparticles flow across a stretched sheet in the existence of the porous media, HS/S, TPD, and bioconvection has been done. The present paper is presented to study the influence of nanoparticle aggregation and without aggregation on various parameters over respective profiles. The governing equations are converted into a system of ODEs and solved numerically. The influence of various parameters is studied in detail. The current examination is conducted to obtain the answers to the following questions.
  • What is the influence of nanoparticles with and without aggregation on porosity parameter over a velocity profile?
  • What is the impact of nanoparticles with and without aggregation over a thermal distribution?
  • What is the effect of HS/S parameter, thermophoretic parameter, bioconvection Peclet number and bioconvection Lewis number on respective profiles?

2. Mathematical Formulation

The steady, laminar, incompressible boundary layer flow of a micropolar nanofluid movement over a continuous stretching sheet in the presence of the porous medium, HS/S, TPD and bioconvection is considered. The flow configuration of the model is schematically revealed in Figure 1, where x and y represents the respective Cartesian coordinate axes. The x axis is taken along the stretching surface and y axis is measured normal to it. Therefore, the wall uniform stretching velocity of the sheet is specified by u w = a x , where a is a positive constant. In addition, the constant temperature of the fluid at the wall and ambient scenario are represented via T w and T , respectively. The constant concentration and motile microorganism at the wall surface of the sheet is dignified via C w and N w , respectively, while C and N denoted the respective uniform ambient concentration and the motile microorganism. Under these requisite posited assumptions, the governing equations for continuity, momentum, temperature, concentration and bioconvection in terms of PDEs are stated as (see, e.g., in [37,38,39]).
u x + v y = 0 ,
u u x + v u y = K 1 ρ n f H y υ n f K * u + 2 u y 2 K 1 + μ n f 1 ρ n f ,
u H x + v H y = 2 H y 2 Ω 1 ρ n f j 2 H + u y K 1 ρ n f j ,
u T x + v T y = k n f ρ C p n f 2 T y 2 + Q 1 ρ C p n f T T ,
u C x + v C y = D B 2 C y 2 y V T C C ,
u N x + v N y = D m 2 N y 2 y b W c C w C N C y ,
along with subject to the boundary conditions
u = u w , v = 0 , N = N w , H = n 1 u y , T = T w , C = C w   at   y = 0 ,
u 0 , T T , N N , H 0 , C C   as   y .
Furthermore, the spin gradient viscosity is given by Ω 1 = μ f μ n f μ f + α 1 2 j μ n f + K 1 2 j (see, e.g., in [31]), the micropolar parameter is given as α 1 = K 1 μ f and j = υ f a signifies the reference length. The microrotation parameter is denoted by n 1 , where n 1 0 , 1 . In the absence of microrotation parameter, the microelements are highly concentrated and unable to spin close to the wall surface.
According to the work in [32], the thermophoretic velocity is demarcated as
V T = υ f K 2 T r T y
where K 2 and T r are the thermophoretic constant and reference temperature, respectively.
According to experimental and empirical investigations (see, e.g., in [19,40,41]), the nanoparticle aggregation factor is important in the dynamics and heat transmission of nanofluid flows. When the aggregation component was taken into consideration, the measurement results of the nanomaterial agreed precisely. With the aggregation kinetic factor, the nanoparticle volume fraction becomes
ϕ a 1 = ϕ r a 1 r p 1 3 D 1 .
The thermal conductivity is computed by combining the Bruggeman and modified Maxwell models, as indicated in Table 1 (see [42]) and the appropriate expressions are provided below (see, e.g., in [43,44,45]):
k a 1 k f = 0.25 k p 1 k f 3 ϕ i 1 1 + 3 1 ϕ i 1 1 + k p 1 k f 3 ϕ i 1 1 + 3 1 ϕ i 1 1 2 + 8 k p 1 k f 0.5 ,
ϕ i 1 = r a 1 r p 1 D 1 3 ,
ρ a 1 = 1 ϕ i 1 ρ f + ϕ i 1 ρ s ,
ρ C p a 1 = 1 ϕ i 1 ρ C p f + ϕ i 1 ρ C p s .
From Equations (12)–(14), the subscript a 1 denotes aggregates and p 1 signifies nanoparticles. Fractional index D 1 = 1.8 , r a 1 r p 1 = 3.34 , ϕ m 1 = 0.605 reflect fast-moving flows and about the mono disperse system are the commonly accepted values.
Furthermore, the similarity transformations for the considered requisite problem are defined as
Ψ = x f η a υ f , u = Ψ y , v = Ψ x ,   H = a x a υ f g η , η = a υ f y ,   θ η = T T w T w T , φ η = N N w N w N , χ η = C C w C w C .
Implementing the Equation (15) into leading governing equations, we get the following reduced form of similarity equations:
f 1 A 1 + α 1 + A 2 f f f 2 + α 1 g β 1 A 1 f = 0 ,
g 1 A 1 + α 1 2 + A 2 f g f g α 1 2 g + f = 0 ,
k h n f k f θ + A 3 Pr f θ + Pr H s θ = 0 ,
χ S c + f χ τ θ χ + θ χ = 0 ,
φ + L b f φ P e χ φ + χ ( σ + φ ) = 0 ,
with subject transformed BCs:
f 0 = 0 , g 0 = n 1 f 0 , θ 0 = χ 0 = φ 0 = f 0 = 1 ,   at   η = 0 ,
f η 0 , θ η 0 , φ η 0 , g η 0 , χ η 0   as   η .
Similarly, Equations (16)–(22) comprised distinct controlling parameters which are namely and symbolically given as follows: β = υ f / a K * signifies porous constraint, H s = Q 1 / a ρ C p f signifies heat source/sink constraint, Pr = υ f / α f signifies Prandtl number, S c = υ f / D B signifies Schmidt number, τ = K 2 T w T / T r signifies thermophoretic constraint, L b = υ f / D m signifies bioconvection Lewis number, P e = b W c / D m signifies bioconvection Peclet number, σ = N / N w N signifies concentration difference parameter and the other thermophysical constants expressions are demarcated as
A 1 = 1 ϕ a 1 ϕ m 1 2.5 ϕ m 1 , A 2 = 1 ϕ a 1 + ϕ a 1 ρ a 1 ρ f and A 3 = 1 + ρ C p a 1 ρ C p f ϕ a 1 ϕ a 1 .
The important engineering factors like Skin friction, wall couple stress factor, Nusselt number, Sherwood number and density number of motile microorganisms is given as
C f = τ w ρ f u w 2 , C s = M w μ f u w , N u x = x q w k f T w T S h x = x q m D B C w C , N n x = x q n D m N w N
The terms τ w , M w , q w , q m ,   and   q n defined as
τ w = μ n f + K 1 u y + K 1 H y = 0
M w = H y y = 0 × μ n f + K 1 2 j
q w = T y y = 0 × k n f
q m = C y y = 0 × D B
q n = D m N y y = 0
Using Equations (15) and (21) the above expressions reduced as follows:
C f = 1 A 1 + α 1 1 n 1 f 0 Re , C s = 1 A 1 + α 1 2 g 0 , N u x = k n f k f θ 0 Re , S h x = Re χ 0 , N n x = Re φ 0 .

3. Numerical Process and Authentication of Code

The RKF-45 scheme and the procedure of shooting are exercised to tackle the set of Equations (16)–(20) and BCs (21) and (22). Because the resultant equations are of higher order and have two points. We must first convert this to an initial value concern in order to solve it. Let us consider the new variables:
f = h 1 1 , f = h 1 2 , f = h 1 3 , f = h 1 3 , g = h 1 4 , g = h 1 5 , g = h 1 5 , θ = h 1 6 , θ = h 1 7 , θ = h 1 7 , χ = h 1 8 , χ = h 1 9 , χ = h 1 9 , φ = h 1 10 , φ = h 1 11 , φ = h 1 11 .
h 3 = A 2 h 1 1 h 1 3 h 1 2 2 + α 1 h 1 5 β A 1 h 1 2 / 1 A 1 + α 1
h 5 = A 2 h 1 1 h 1 5 h 1 2 h 1 4 α 1 2 h 1 4 + h 1 3 / 1 A 1 + α 1 2
h 7 = A 3 h 1 1 h 1 7 + H s h 1 6 / k n f k f Pr
h 9 = S c h 1 1 h 1 9 τ h 1 7 h 1 8 + h 1 7 h 1 9
h 11 = L b h 1 1 h 1 11 P e h 1 9 h 1 11 + h 1 9 ( σ + h 1 10 ) .
with
h 1 1 0 = 0 , h 1 2 0 = 1 , h 1 3 0 = h 1 4 0 / n 1 , h 1 4 0 = n 1 h 1 3 0 , h 1 5 0 = γ 1 , h 1 6 0 = 1 , h 1 7 0 = γ 2 , h 1 8 0 = 1 , h 1 9 0 = γ 3 , h 1 10 0 = 1 , h 1 11 0 = γ 4 .
The transformed initial value problem Equations (30)–(35)are numerically solved by using the thermo-physical characteristics of nanofluid provided in Table 1 (see [42]) and thermophysical properties of nanomaterial and base fluid given in Table 2 (see [38]) by varying the parameter values, α 1 = 0.1 , β = 0.1 , H s = 0.5 , S c = 0.8 , τ = 0.1 , L b = P e = 1 , σ = 0.1 , and n 1 = 1 quantitatively solved by estimating the unknowns using the shooting approach. With a step size of 0.1 and an acceptance inaccuracy of approximately 10 8 . The numerical results of the present study have been compared with results of previous literature [46,47,48] and found to be the better match (see Table 3).

4. Analysis of Results

This section outlines that how the numerous regulating factors determine velocity and temperature profiles, concentration, motile microorganism profiles, and engineering coefficients are examined in greater details. The analysis is conducted in two modes namely, (1) with aggregation ϕ i 1 1 , and (2) without aggregation over a porous, micro- polar constraint, heat sink/source constraint, thermophoretic constraint, bio-convection Lewis number and bio-convection Peclet number.
Figure 2 illustrates the effect of porous parameter β on a dimensionless velocity profile f . Inclination in the values of β will declines f . This is because that the drag force faces the profiles of velocity movement in the existence of porous space, causing the speed to decline. Moreover, the diagram illustrates that in the absence of nanoparticle aggregation, the fluid velocity is lower than in the presence of nanoparticle aggregation.
Figure 3 demonstrations the influence of micropolar constraint α 1 over a velocity f and micro rotation velocity fields g . Further, noted from the figure that when the micropolar constraint α 1 is small, the viscous impact is limited to a very thin layer close to the wall, and the thickness of the boundary layer is equivalent to α 1 and inversely relative to micro rotation velocity g . Upgrading in the values of α 1 will declines the velocity profile but escalates the micro rotation velocity form. It is also observed that velocity is more in the presence of aggregation than absence of aggregation in f but, reverse trend is seen in g .
The discrepancy of HS/S constraint over a temperature profile is illustrated in the Figure 4. Where heat dispersion is a major critical factor in engineering, this profile plays a vital role. Up gradation in the HS/S parameter will enhance the thermal dispersal in the given system. The heat sink will behave like a functional as a confine changer, transferring the heat created by the body into the nanoliquid. Consequently, in a given circumstance of a HS, the thermal dispersion is diffident, but in the occasion of a heat source, the heat is created by the body surface. In the occurrence of a heat source, higher thermal act is seen than in the absence of a heat sink. The graphic clearly shows that heat dispersion in nanoparticles with aggregation is greater than in nanoparticles without aggregation.
Figure 5 shows the discrepancy of thermophoretic constraint over a concentration profile. Escalating values of will diminishes the concentration. The motion of the nanoparticles enhances as the temperature gradient rises, ensuing in a lower concentration. The figure also shows that nanoparticles with aggregation have a lower concentration than nanoparticles without aggregation.
Figure 6 illustrates the variation of bioconvection Peclet number P e and bioconvection Lewis number L b over bioconvection profile φ . Rise in the values in P e and L b will decreases the φ . The percentage of thermal diffusivity with mass diffusivity is known as the Lewis number. It is sometimes stated as a ratio between the Schmidt number and the Prandtl number. P e and L b results in decrease in microorganism dispersion, while density and BLT for motile microorganisms are decreasing. The figure also shows that nanoparticles with aggregation have a lower bioconvection profile than nanoparticles without aggregation.
Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show the effect of several dimensionless constraints on important primary engineering coefficients. Figure 7 depicts the effect of micropolar constraint on surface drag force for various porosity constraint values. It is detected from the figure that surface drag force is less in case of nanoparticles without aggregation than with aggregation. This is because the existence of a porous parameter causes drag on the fluid flow, and an increase in the micropolar constraint increases the boundary layer thickness. As a result, surface drag force declines. Figure 8 exhibit the influence of micro-rotation parameter on wall couple stress factor for various values of micropolar constraint. It is seen from the figure that the rise in the values of micro-rotation parameter and micropolar constraint will enhances the wall couple stress in the system. This is due to enhancement in the vortex viscosity factor.
Figure 9 displays the consequence of solid volume fraction ϕ over a temperature distribution rate for numerous values of heat source/sink constraint H s . Escalating in the values of H s and ϕ will upraises rate of thermal circulation. Augmentation in the precise values of ϕ will improves the BLT and increased H s will improves the heat propagation. As a result, rate thermal distribution improves.
Figure 10 depicts the nature of ϕ over a mass transfer rate for numerous values of thermophoretic constraint τ . An increase in ϕ and τ will decrease the rate of mass transmission. As the values of thermophoretic constraint increases, temperature difference gradient also increases, as a result, the particle motion in the system increases. Oppositely, enhancement in the solid volume fraction will enhances the BLT. Due to presence of thermophoretic constraint, the rate of mass transfer diminishes. Figure 11 illustrates the consequence bioconvection Lewis number L b on density number of motile microorganisms for various values of bioconvection Peclet number P e . It is observed that, increased values of P e will enhance the density number of motile microorganisms.

5. Concluding Remarks

In the present investigation micropolar nanoliquid flow past a continuous stretching sheet subject to porous medium, HS/S, TPD, and bioconvection is examined. Here, the analysis is made in the presence and absence of nanoparticle aggregation. The system of equations that represents governing equations are transformed into ODEs with utilizing suitable similarity variables and those equations are solved numerically with the aid of computing software. The influence of various dimensionless parameters is studied with respective profiles. The major outcomes of the present investigation are as follows:
  • Enhancement in the porous parameter will diminishes the velocity profile due to presence of porous medium which drags the fluid motion.
  • For escalating values of α 1 velocity is more in the presence of aggregation than absence of aggregation in f but, reverse trend is seen in g .
  • Improvement in the values of H s will improve the profiles of temperature and heat transfer rate. Heat transfer gradually increases from heat sink to source.
  • Nanoparticles with aggregation have a lower concentration than nanoparticles without aggregation in the presence of thermophoretic parameter.
  • Density number of motile microorganisms will be decreased by improved values of P e .
  • Nanoparticle with aggregation is lesser impact in velocity profile but shows more impact in temperature, concentration and bioconvection profiles.

Author Contributions

Conceptualization, J.K.M. and U.K.; methodology, J.K.M.; software, J.K.M. and U.K.; validation, J.K.M., Y.Y., U.K. and F.W.; formal analysis, A.-H.A.-A., A.Z., A.M.G., I.S.Y. and M.S.A.; investigation, A.M.G., Y.Y., A.-H.A.-A. and I.S.Y.; resources, A.M.G.; data curation, A.Z.; writing—original draft preparation, U.K., J.K.M. and F.W.; writing—review and editing, Y.Y., M.S.A., A.Z., A.M.G. and M.S.A.; visualization, Y.Y., A.Z., M.S.A. and F.W.; supervision, A.Z.; project administration, A.-H.A.-A.; funding acquisition, A.-H.A.-A., I.S.Y. and M.S.A. All authors have read and agreed to the published version of the manuscript.

Funding

The authors express their appreciation to The Research Center for Advanced Materials Science (RCAMS) at King Khalid University, Saudi Arabia, for funding this work under the grant number RCAMS/KKU/018-20. Also, this work is partially supported by the National Natural Science Foundation of China (No. 12071408) and the Natural Science Foundation of Jiangsu 419 Province (No. BK20201149). Furthermore, this work is also partially supported by the Natural Science Foundation of Anhui Province (No. 2008085MA11).

Data Availability Statement

Not Applicable.

Acknowledgments

The Research Center for Advanced Materials Science (RCAMS) at King Khalid University, Saudi Arabia, for funding this work under the grant number RCAMS/KKU/018-20. Also, this work is partially supported by the National Natural Science Foundation of China (No. 12071408) and the Natural Science Foundation of Jiangsu 419 Province (No. BK20201149). Furthermore, this work is also partially supported by the Natural Science Foundation of Anhui Province (No. 2008085MA11). We would like to thank the reviewers for their thoughtful comments and efforts towards improving our paper.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

a Stretching Constant n 1 Microrotation Parameter
C Concentration N w The density of motile microorganisms at the wall
C w & C Wall and ambient concentration N The ambient density of motile microorganisms
C p Specific heat N n x Density number of motile microorganisms
C f Skin friction N u x Nusselt number
C s Wall couple stress factor P e Bioconvection Peclet number
D Diffusivity Pr Prandtl number
D 1 Fractal index r a 1 Radii of aggregates
f η Dimensionless velocity profile r p 1 Radii of nanoparticle
g η Dimensionless microrotation velocity profile Re local Reynolds number
H s Heat source/sink constraint S c Schmidt number
j Microinertia density S h x Sherwood number
k Thermal conductivity T Temperature
k a 1 Thermal aggregation conductivity T r Reference temperature
K 1 Coefficient of vortex viscosity T w & T Wall and ambient temperature
K 2 Thermophoretic constant u w Uniform velocity
K * Permeability of porous media V T Thermophoretic velocity
L b Bioconvection Lewis number u & v Velocity components
N The density of motile microorganism x & y Coordinate axis
Greek symbols
Ω 1 Spin gradient viscosity α 1 Micropolar parameter
β Porous parameter τ Thermophoretic parameter
ν Kinematic viscosity σ Concentration difference parameter
μ Dynamic viscosity ρ Density
ϕ m 1 Extreme volume fraction ϕ a 1 Effective volume fraction of aggregates
η Similarity variable Ψ Stream function
θ ( η ) Dimensionless temperature profile φ η Dimensionless bioconvection profile
χ ( η ) Dimensionless concentration profile ϕ Solid volume fraction
Subscripts
a 1 aggregates p 1 Nanoparticle
f Fluid n f Nanofluid

References

  1. Mehmood, R.; Rana, S.; Akbar, N.; Nadeem, S. Non-aligned stagnation point flow of radiating Casson fluid over a stretching surface. Alex. Eng. J. 2018, 57, 939–946. [Google Scholar] [CrossRef]
  2. Eringen, A.C. Theory of Micropolar Fluids. J. Math. Mech. 1966, 16, 1–18. [Google Scholar] [CrossRef]
  3. Habib, D.; Abdal, S.; Ali, R.; Baleanu, D.; Siddique, I. On bioconvection and mass transpiration of micropolar nanofluid dynamics due to an extending surface in existence of thermal radiations. Case Stud. Therm. Eng. 2021, 27, 101239. [Google Scholar] [CrossRef]
  4. Lund, L.A.; Omar, Z.; Khan, U.; Khan, I.; Baleanu, D.; Nisar, K.S. Stability Analysis and Dual Solutions of Micropolar Nanofluid over the Inclined Stretching/Shrinking Surface with Convective Boundary Condition. Symmetry 2020, 12, 74. [Google Scholar] [CrossRef] [Green Version]
  5. Madhukesh, J.K.; Prasannakumara, B.C.; Kumar, R.S.V.; Rauf, A.; Shehzad, S.A. Flow of hydromagnetic micropolar-casson nanofluid over porous disks influenced by cattaneo-christov theory and slip effects. J. Porous Media 2022, 25, 35–49. [Google Scholar] [CrossRef]
  6. Mabood, F.; Shamshuddin; Mishra, S. Characteristics of thermophoresis and Brownian motion on radiative reactive micropolar fluid flow towards continuously moving flat plate: HAM solution. Math. Comput. Simul. 2021, 191, 187–202. [Google Scholar] [CrossRef]
  7. Zaib, A.; Khan, U.; Shah, Z.; Kumam, P.; Thounthong, P. Optimization of entropy generation in flow of micropolar mixed convective magnetite (Fe3O4) ferroparticle over a vertical plate. Alex. Eng. J. 2019, 58, 1461–1470. [Google Scholar] [CrossRef]
  8. Khan, M.I.; Waqas, H.; Khan, S.U.; Imran, M.; Chu, Y.-M.; Abbasi, A.; Kadry, S. Slip flow of micropolar nanofluid over a porous rotating disk with motile microorganisms, nonlinear thermal radiation and activation energy. Int. Commun. Heat Mass Transf. 2021, 122, 105161. [Google Scholar] [CrossRef]
  9. Bahiraei, M.; Naseri, M.; Monavari, A. A second law analysis on flow of a nanofluid in a shell-and-tube heat exchanger equipped with new unilateral ladder type helical baffles. Powder Technol. 2021, 394, 234–249. [Google Scholar] [CrossRef]
  10. Bahiraei, M.; Naseri, M.; Monavari, A. Thermal-hydraulic performance of a nanofluid in a shell-and-tube heat exchanger equipped with new trapezoidal inclined baffles: Nanoparticle shape effect. Powder Technol. 2021, 395, 348–359. [Google Scholar] [CrossRef]
  11. Khan, U.; Zaib, A.; Abu Bakar, S.; Ishak, A. Unsteady stagnation-point flow of a hybrid nanofluid over a spinning disk: Analysis of dual solutions. Neural Comput. Appl. 2022, 89, 1–18. [Google Scholar] [CrossRef]
  12. Kumar, R.N.; Suresha, S.; Gowda, R.J.P.; Megalamani, S.B.; Prasannakumara, B.C. Exploring the impact of magnetic dipole on the radiative nanofluid flow over a stretching sheet by means of KKL model. Pramana-J. Phys. 2021, 95, 180. [Google Scholar] [CrossRef]
  13. Alhadhrami, A.; Alzahrani, H.A.H.; Kumar, R.N.; Gowda, R.J.P.; Sarada, K.; Prasanna, B.M.; Madhukesh, J.K.; Madhukeshwara, N. Impact of thermophoretic particle deposition on Glauert wall jet slip flow of nanofluid. Case Stud. Therm. Eng. 2021, 28, 101404. [Google Scholar] [CrossRef]
  14. Madhukesh, J.K.; Ramesh, G.K.; Prasannakumara, B.C.; Shehzad, S.A.; Abbasi, F.M. Bio-Marangoni convection flow of Casson nanoliquid through a porous medium in the presence of chemically reactive activation energy. Appl. Math. Mech.-Engl. Ed. 2021, 42, 1191–1204. [Google Scholar] [CrossRef]
  15. Muhammad, N.; Nadeem, S. Ferrite nanoparticles Ni-ZnFe2O4, Mn-ZnFe2O4 and Fe2O4 in the flow of ferromagnetic nanofluid. Eur. Phys. J. Plus 2017, 132, 377. [Google Scholar] [CrossRef]
  16. Guo, S.-Z.; Li, Y.; Jiang, J.-S.; Xie, H.-Q. Nanofluids Containing γ-Fe2O3 Nanoparticles and Their Heat Transfer Enhancements. Nanoscale Res. Lett. 2010, 5, 1222–1227. [Google Scholar] [CrossRef] [Green Version]
  17. Xie, H.; Wang, J.; Xi, T.; Liu, Y.; Ai, F.; Wu, Q. Thermal conductivity enhancement of suspensions containing nanosized alumina particles. J. Appl. Phys. 2002, 91, 4568–4572. [Google Scholar] [CrossRef]
  18. Wolthers, W.; Duits, M.H.G.; ven den Ende, D.; Mellema, J. Shear history dependence of the viscosity of aggregated colloidal dispersions. J. Rheol. 1996, 40, 799–811. [Google Scholar] [CrossRef]
  19. Mackolil, J.; Mahanthesh, B. Sensitivity analysis of Marangoni convection in TiO2–EG nanoliquid with nanoparticle aggregation and temperature-dependent surface tension. J. Therm. Anal. 2021, 143, 2085–2098. [Google Scholar] [CrossRef]
  20. Mahanthesh, B.; Thriveni, K.; Rana, P.; Muhammad, T. Radiative heat transfer of nanomaterial on a convectively heated circular tube with activation energy and nanoparticle aggregation kinematic effects. Int. Commun. Heat Mass Transf. 2021, 127, 105568. [Google Scholar] [CrossRef]
  21. Kolsi, L.; Oztop, H.F.; Ghachem, K.; Almeshaal, M.A.; Mohammed, H.A.; Babazadeh, H.; Abu-Hamdeh, N. Numerical Study of Periodic Magnetic Field Effect on 3D Natural Convection of MWCNT-Water/Nanofluid with Consideration of Aggregation. Processes 2019, 7, 957. [Google Scholar] [CrossRef] [Green Version]
  22. Mackolil, J.; Mahanthesh, B. Inclined magnetic field and nanoparticle aggregation effects on thermal Marangoni convection in nanoliquid: A sensitivity analysis. Chin. J. Phys. 2021, 69, 24–37. [Google Scholar] [CrossRef]
  23. Saleh, B.; Madhukesh, J.K.; Kumar, R.S.V.; Afzal, A.; Abdelrhman, Y.; Aly, A.A.; Gowda, R.J.P. Aspects of magnetic dipole and heat source/sink on the Maxwell hybrid nanofluid flow over a stretching sheet. Proc. Proc. Inst. Mech. Eng. Part E J. Process Mech. Eng. 2022, 09544089211056243. [Google Scholar] [CrossRef]
  24. Madhukesh, J.; Ramesh, G.; Aly, E.H.; Chamkha, A.J. Dynamics of water conveying SWCNT nanoparticles and swimming microorganisms over a Riga plate subject to heat source/sink. Alex. Eng. J. 2021, 61, 2418–2429. [Google Scholar] [CrossRef]
  25. Sajid, T.; Jamshed, W.; Shahzad, F.; Eid, M.R.; Alshehri, H.M.; Goodarzi, M.; Akgül, E.K.; Nisar, K.S. Micropolar fluid past a convectively heated surface embedded with nth order chemical reaction and heat source/sink. Phys. Scr. 2021, 96, 104010. [Google Scholar] [CrossRef]
  26. Khan, S.U.; Shehzad, S.A.; Rauf, A.; Ali, N. Mixed convection flow of couple stress nanofluid over oscillatory stretching sheet with heat absorption/generation effects. Results Phys. 2018, 8, 1223–1231. [Google Scholar] [CrossRef]
  27. Chu, Y.-M.; Khan, U.; Shafiq, A.; Zaib, A. Numerical Simulations of Time-Dependent Micro-Rotation Blood Flow Induced by a Curved Moving Surface Through Conduction of Gold Particles with Non-uniform Heat Sink/Source. Arab. J. Sci. Eng. 2021, 46, 2413–2427. [Google Scholar] [CrossRef]
  28. Shehzad, S.; Mabood, F.; Rauf, A.; Tlili, I. Forced convective Maxwell fluid flow through rotating disk under the thermophoretic particles motion. Int. Commun. Heat Mass Transf. 2020, 116, 104693. [Google Scholar] [CrossRef]
  29. Kumar, R.N.; Gowda, R.J.P.; Madhukesh, J.K.; Prasannakumara, B.C.; Ramesh, G.K. Impact of thermophoretic particle deposition on heat and mass transfer across the dynamics of Casson fluid flow over a moving thin needle. Phys. Scr. 2021, 96, 075210. [Google Scholar] [CrossRef]
  30. Shankaralingappa, B.M.; Madhukesh, J.K.; Sarris, I.E.; Gireesha, B.J.; Prasannakumara, B.C. Influence of Thermophoretic Particle Deposition on the 3D Flow of Sodium Alginate-Based Casson Nanofluid over a Stretching Sheet. Micromachines 2021, 12, 1474. [Google Scholar] [CrossRef]
  31. Hafeez, A.; Khan, M.; Ahmed, J. Oldroyd-B fluid flow over a rotating disk subject to Soret–Dufour effects and thermophoresis particle deposition. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2021, 235, 2408–2415. [Google Scholar] [CrossRef]
  32. Chen, S.-B.; Shahmir, N.; Ramzan, M.; Sun, Y.-L.; Aly, A.A.; Malik, M.Y. Thermophoretic particle deposition in the flow of dual stratified Casson fluid with magnetic dipole and generalized Fourier’s and Fick’s laws. Case Stud. Therm. Eng. 2021, 26, 101186. [Google Scholar] [CrossRef]
  33. Ali, B.; Shafiq, A.; Manan, A.; Wakif, A.; Hussain, S. Bioconvection: Significance of mixed convection and mhd on dynamics of Casson nanofluid in the stagnation point of rotating sphere via finite element simulation. Math. Comput. Simul. 2021, 194, 254–268. [Google Scholar] [CrossRef]
  34. Alqarni, M.S.; Waqas, H.; Alghamdi, M.; Muhammad, T. Importance of bioconvection in 3D viscoelastic nanofluid flow due to exponentially stretching surface with nonlinear radiative heat transfer and variable thermal conductivity. J. Therm. Anal. 2021, 147, 4805–4819. [Google Scholar] [CrossRef]
  35. Wakif, A.; Animasaun, I.L.; Khan, U.; Alshehri, A.M. Insights into the Generalized Fourier’s and Fick’s Laws for Simulating Mixed Bioconvective Flows of Radiative-Reactive Walters-B Fluids Conveying Tiny Particles subject to Lorentz Force. Preprint 2021. [Google Scholar] [CrossRef]
  36. Farooq, U.; Waqas, H.; Khan, M.I.; Khan, S.U.; Chu, Y.-M.; Kadry, S. Thermally radioactive bioconvection flow of Carreau nanofluid with modified Cattaneo-Christov expressions and exponential space-based heat source. Alex. Eng. J. 2021, 60, 3073–3086. [Google Scholar] [CrossRef]
  37. Naganthran, K.; Basir, F.M.; Thumma, T.; Ige, E.O.; Nazar, R.; Tlili, I. Scaling group analysis of bioconvective micropolar fluid flow and heat transfer in a porous medium. J. Therm. Anal. 2021, 143, 1943–1955. [Google Scholar] [CrossRef]
  38. Kumar, T.; Kumar, B.R.; Makinde, O.D.; Kumar, A.V. Magneto-Convective Heat Transfer in Micropolar Nanofluid over a Stretching Sheet with Non-Uniform Heat Source/Sink. Defect Diffus. Forum 2018, 387, 78–90. [Google Scholar] [CrossRef]
  39. Alam Khan, N.; Sultan, F.; Alam Khan, N. Heat and Mass Transfer of Thermophoretic MHD Flow of Powell–Eyring Fluid over a Vertical Stretching Sheet in the Presence of Chemical Reaction and Joule Heating. Int. J. Chem. React. Eng. 2015, 13, 37–49. [Google Scholar] [CrossRef]
  40. Mahanthesh, B. Flow and heat transport of nanomaterial with quadratic radiative heat flux and aggregation kinematics of nanoparticles. Int. Commun. Heat Mass Transf. 2021, 127, 105521. [Google Scholar] [CrossRef]
  41. Chen, J.; Zhao, C.; Wang, B. Effect of nanoparticle aggregation on the thermal radiation properties of nanofluids: An experimental and theoretical study. Int. J. Heat Mass Transf. 2020, 154, 119690. [Google Scholar] [CrossRef]
  42. Mahanthesh, B.; Thriveni, K. Nanoparticle aggregation effects on radiative heat transport of nanoliquid over a vertical cylinder with sensitivity analysis. Appl. Math. Mech.-Engl. Ed. 2021, 42, 331–346. [Google Scholar] [CrossRef]
  43. Ellahi, R.; Hassan, M.; Zeeshan, A. Aggregation effects on water base Al2O3-nanofluid over permeable wedge in mixed convection: Aggregation effects on water base Al2O3. Asia-Pacific J. Chem. Eng. 2016, 11, 179–186. [Google Scholar] [CrossRef]
  44. Chen, H.; Witharana, S.; Jin, Y.; Kim, C.; Ding, Y. Predicting thermal conductivity of liquid suspensions of nanoparticles (nanofluids) based on rheology. Particuology 2009, 7, 151–157. [Google Scholar] [CrossRef]
  45. Chen, H.; Ding, Y.; He, Y.; Tan, C. Rheological behaviour of ethylene glycol based titania nanofluids. Chem. Phys. Lett. 2007, 444, 333–337. [Google Scholar] [CrossRef]
  46. Pal, D.; Mandal, G. Thermal radiation and MHD effects on boundary layer flow of micropolar nanofluid past a stretching sheet with non-uniform heat source/sink. Int. J. Mech. Sci. 2017, 126, 308–318. [Google Scholar] [CrossRef]
  47. Fauzi, E.L.A.; Ahmad, S.; Pop, I. Flow over a permeable stretching sheet in micropolar nanofluids with suction. In Proceedings of the AIP Conference Proceedings, American Institute of Physics, Penang, Malaysia, 10 July 2014; Volume 1605, pp. 428–433. [Google Scholar] [CrossRef]
  48. Nazar, R.; Amin, N.; Filip, D.; Pop, I. Stagnation point flow of a micropolar fluid towards a stretching sheet. Int. J. Non-linear Mech. 2004, 39, 1227–1235. [Google Scholar] [CrossRef]
Figure 1. Physical model of the problem.
Figure 1. Physical model of the problem.
Nanomaterials 12 00977 g001
Figure 2. Sway of f for different values of β .
Figure 2. Sway of f for different values of β .
Nanomaterials 12 00977 g002
Figure 3. Sway of f & g for different values of α 1 .
Figure 3. Sway of f & g for different values of α 1 .
Nanomaterials 12 00977 g003
Figure 4. Sway of θ for different values of H s .
Figure 4. Sway of θ for different values of H s .
Nanomaterials 12 00977 g004
Figure 5. Sway of χ for different values of τ .
Figure 5. Sway of χ for different values of τ .
Nanomaterials 12 00977 g005
Figure 6. Sway of φ for different values of P e & L b .
Figure 6. Sway of φ for different values of P e & L b .
Nanomaterials 12 00977 g006
Figure 7. Nature of Skin friction C f over Micropolar parameter α 1 for numerous values of porous parameter β .
Figure 7. Nature of Skin friction C f over Micropolar parameter α 1 for numerous values of porous parameter β .
Nanomaterials 12 00977 g007
Figure 8. Nature of wall couple stress factor C s over Micropolar parameter α 1 for numerous values ofmicrorotation parameter n 1 .
Figure 8. Nature of wall couple stress factor C s over Micropolar parameter α 1 for numerous values ofmicrorotation parameter n 1 .
Nanomaterials 12 00977 g008
Figure 9. Nature of Nusselt number N u x over solid volume fraction ϕ for numerous values of heat source/sink parameter H s .
Figure 9. Nature of Nusselt number N u x over solid volume fraction ϕ for numerous values of heat source/sink parameter H s .
Nanomaterials 12 00977 g009
Figure 10. Nature of Sherwood number S h x over solid volume fraction ϕ for numerous values of thermophoretic parameter τ .
Figure 10. Nature of Sherwood number S h x over solid volume fraction ϕ for numerous values of thermophoretic parameter τ .
Nanomaterials 12 00977 g010
Figure 11. Nature of density number of motile microorganisms N n x over bioconvection Peclet number P e for numerous values of bioconvection Lewis number L b .
Figure 11. Nature of density number of motile microorganisms N n x over bioconvection Peclet number P e for numerous values of bioconvection Lewis number L b .
Nanomaterials 12 00977 g011
Table 1. Thermophysical properties of nanomaterials.
Table 1. Thermophysical properties of nanomaterials.
S. NoPropertiesExpressions for Nanofluid
01Viscosity μ n f μ f = 1 ϕ a 1 ϕ m 1 2.5 ϕ m 1
02Density ρ n f ρ f = 1 ϕ a 1 + ϕ a 1 ρ a 1 ρ f
03Heat capacity ρ c p n f ρ c p f = ϕ a 1 ρ c p a 1 ρ c p f + 1 ϕ a 1
04Thermal conductivity k n f k f = 2 k f + k a 1 + 2 ϕ a 1 k a 1 k f 2 k f + k a 1 ϕ a 1 k a 1 k f
Table 2. Thermo-physical properties of nanomaterial and base fluid.
Table 2. Thermo-physical properties of nanomaterial and base fluid.
PropertiesTitanium DioxideWater
ρ k g / m 3 4250997.1
C p J   k g K 686.24179
k W /   m K 8.95380.613
Pr-6.2
Table 3. Assessment of Re C f with published results for sundry values of α 1 and n 1 in the neglecting of A 1 , A 2 and β .
Table 3. Assessment of Re C f with published results for sundry values of α 1 and n 1 in the neglecting of A 1 , A 2 and β .
n 1 α 1 Present Study
0.00.0−1.0000
-1.0−1.3678
-2.0−1.6211
-4.0−2.0040
0.50.0−1.0000
-1.0−1.2246
-2.0−1.4140
-4.0−1.7319
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yu, Y.; Madhukesh, J.K.; Khan, U.; Zaib, A.; Abdel-Aty, A.-H.; Yahia, I.S.; Alqahtani, M.S.; Wang, F.; Galal, A.M. Nanoparticle Aggregation and Thermophoretic Particle Deposition Process in the Flow of Micropolar Nanofluid over a Stretching Sheet. Nanomaterials 2022, 12, 977. https://doi.org/10.3390/nano12060977

AMA Style

Yu Y, Madhukesh JK, Khan U, Zaib A, Abdel-Aty A-H, Yahia IS, Alqahtani MS, Wang F, Galal AM. Nanoparticle Aggregation and Thermophoretic Particle Deposition Process in the Flow of Micropolar Nanofluid over a Stretching Sheet. Nanomaterials. 2022; 12(6):977. https://doi.org/10.3390/nano12060977

Chicago/Turabian Style

Yu, Yangyang, Javali K. Madhukesh, Umair Khan, Aurang Zaib, Abdel-Haleem Abdel-Aty, Ibrahim S. Yahia, Mohammed S. Alqahtani, Fuzhang Wang, and Ahmed M. Galal. 2022. "Nanoparticle Aggregation and Thermophoretic Particle Deposition Process in the Flow of Micropolar Nanofluid over a Stretching Sheet" Nanomaterials 12, no. 6: 977. https://doi.org/10.3390/nano12060977

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop