Next Article in Journal
Core/Shell Glycine-Polyvinyl Alcohol/Polycaprolactone Nanofibrous Membrane Intended for Guided Bone Regeneration: Development and Characterization
Next Article in Special Issue
Evaluation of the Durability of Slippery, Liquid-Infused Porous Surfaces in Different Aggressive Environments: Influence of the Chemical-Physical Properties of Lubricants
Previous Article in Journal
Yttria-Coated Tungsten Fibers for Use in Tungsten Fiber-Reinforced Composites: A Comparative Study on PVD vs. CVD Routes
Previous Article in Special Issue
Evaluation of Coating Film Formation Process Using the Fluorescence Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Nanomechanical Concepts in Magnetically Guided Systems to Investigate the Magnetic Dipole Effect on Ferromagnetic Flow Past a Vertical Cone Surface

1
Department of Mathematics, Faculty of Science, King Mongkut’s University of Technology Thonburi, 126 Pracha Uthit Road, Bang Mod, Thung Khru, Bangkok 10140, Thailand
2
Department of Mathematical Sciences, Faculty of Physical Sciences, Bayero University, Kano 700241, Nigeria
3
Department of Mathematical Sciences, University of Lakki Marwat, Lakki Marwat 28420, Pakistan
4
KMUTTFixed Point Research Laboratory, Room SCL 802 Fixed Point Laboratory, Department of Mathematics, Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT), Bangkok 10140, Thailand
5
Center of Excellence in Theoretical and Computational Science (TaCS-CoE), Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT), 126 Pracha-Uthit Road, Bang Mod, Thung Khru, Bangkok 10140, Thailand
6
Department of Mathematics and Statistics, Hazara University Mansehra, Khyber Pakhtunkhwa 21120, Pakistan
*
Author to whom correspondence should be addressed.
Coatings 2021, 11(9), 1129; https://doi.org/10.3390/coatings11091129
Submission received: 16 August 2021 / Revised: 10 September 2021 / Accepted: 12 September 2021 / Published: 16 September 2021

Abstract

:
Because of the floating magnetic nanomaterial, ferrofluids have magneto-viscous properties, enabling controllable temperature changes as well as nano-structured fluid characteristics. The study’s purpose is to evolve and solve a theoretical model of bioconvection nanofluid flow with a magnetic dipole effect in the presence of Curie temperature and using the Forchheimer-extended Darcy law subjected to a vertical cone surface. The model also includes the nonlinear thermal radiation, heat suction/injection, viscous dissipation, and chemical reaction effects. The developed model problem is transformed into nonlinear ordinary differentials, which have been solved using the homotopy analysis technique. In this problem, the behavior of function profiles are graphically depicted and explained for a variety of key parameters. For a given set of parameters, tables representthe expected numerical values and behaviors of physical quantities. The nanofluid velocity decreases as the ferrohydrodynamic, local inertia, and porosity parameters increase and decrease when the bioconvection Rayleigh number increases. Many key parameters improved the thermal boundary layer and temperature. The concentration is low when the chemical reaction parameter and Schmidt number rises. Furthermore, as the bioconvection constant, Peclet and Lewis numbers rise, so does the density of motile microorganisms.

1. Introduction

Fluids that are often magnetized by the existence of an exterior magnetic field are known as ferrofluids, which is an abbreviation for fluid and ferromagnetic particles. These fluids are made up of colloidal fluids formed of nanosized ferromagnetic or ferrimagnetic particles that have been stopped inside the fluid transporter. Brownian motion causes particle suspension and must not start moving under normal conditions. Besides that, to avoid clogging, each ferromagnetic particle is encased in a solvent, and the nano-scaled ferromagnetic particles have a weak magnetic attraction whenever the surfactant’s Van der Waals force adequately stopped aggregation or clustering. Numerous applications of ferromagnetic fluids have emerged in a variety of fields. Heat transfer agents, angular momentum changers, friction reducers, and so on are used in electronic equipment, analytical techniques, and medical science; some examples can be found in the references [1,2,3]. Because of these numerous applications, many researchers and scientists have been focused on this subject. Andersson and Vanes [4] first investigated the influence caused by magnetic dipoles on ferrofluids. Zeeshan et al. [5] investigated the convective heat transfer flow of ferromagnetic fluids with partial slip effects using a stretching sheet. Hayat et al. [6] reported on radiation and magnetic dipole effects of Williamson ferromagnetic fluid flow across a stretched surface.
A nanofluid is a nanometer-sized particle suspended in a fluid. Choi [7] established the basic extension of “nanofluid”, and scientific results verified that heat transfer can be significantly enhanced through the mixture of tiny metallic nanomaterials with the working fluids. A few studies in particular on nanofluids have been conducted. Ellahi [8] performed an analytical study and concluded that the temperature variable and viscosity affects MHD non-Newtonian nanofluid flow in a pipe. Ellahi et al. [9] presented peristaltic nanofluid flow with entropy generation via a medium of porosity. Hayat et al. [10] investigated the flow of third-grade nanofluids caused by a rotating stretchable disk containing a heat source and a chemical reaction. Awais et al. [11] explored the effects of magnetohydrodynamics on peristaltic ciliary-induced flow coatings to rheological hybrid nanofluids. Reddy et al. [12] studied the boundary layer naturally convective MHD nanofluid flow along a vertical cone under the influence of chemical reaction and heat suction/injection.
Due to disorganized frameworks and destabilization, low-density microorganisms remain on the surface of a fluid, causing bioconvection. Because nanoparticles move differently than motile microorganisms, the cumulative importance of nanomaterials and bioconvection is such that they play a vital role in microfluidic devices. Bioconvection is a novel manufacturing and fluid mechanic with a biological phenomenon involving gyrotactic microorganisms. As a result, it becomes an interesting field of research to which many researchers continue to pay attention. Alsaedi et al. [13] investigated stratified magnetohydrodynamic nanofluid flow, causing bioconvection in gyrotactic microorganisms. Hayat et al. [14] researched the magnetohydrodynamic (MHD) nonlinear radiative nanofluid flow with gyrotactic microorganisms. Nadeem et al. [15] reported on the Rosseland assessment for ferromagnetic fluid with involvement of magnetic dipoles and gyrotactic microorganisms. Bhatti and Michaelides [16] researched thermo-bioconvection nanofluid flows across a Riga plate as a function of Arrhenius activation energy. Waqas et al. [17] have also numerically simulated the magnetized non-Newtonian bioconvection nanofluid flow along stretching cylinders/plates.
Combining mass and heat fluxes in liquid saturated porous media is crucial among a wide range of engineering procedures such as heating systems, oil and gas reservoirs, and chemical catalytic reactor designs [18,19]. The dragging force, the Darcy–Forchheimer technique, is a widely popular method for simulating fluid passed through a porous medium with high velocity. In the literature, flow through a cone in Darcy–Forchheimer porous media has already been analyzed by many researchers. Kumar et al. [20] investigated the non-Darcy MHD viscoelastic fluid flow through a flat plate and a vertical cone. Chamkha et al. [21] explored the non-Newtonian natural convective nanofluid flow over a saturated cone in a non-Darcy porous medium with uniform volume fraction and heat fluxes. Mallikarjuna et al. [22] researched the impacts of radiation, thermophoresis and transpiration on convective non-Darcy flow via a rotating cone. Durairaj et al. [23] investigated the chemically reacting Casson fluid of a non-Darcy porous medium flow through a flat plate and a vertical cone saturated with heat generating/absorbing. Patrulescu et al. [24] investigated a convection flow due to a vertical plate embedded in a bi-disperse non-Darcy porous medium.
According to a recent literature review, despite important applications in extrusion systems, geothermics, organic compounds, geophysics, improved manufacturing techniques, material processing, and improved energy generation, research on viscous ferrofluid flows via a linear vertical cone with consideration of Darcy–Forchheimer porous media has been studied by very few researchers in the past. The aim of this research is to use the Forchheimer-extended Darcy law to explore the effect of magnetic dipole and porosity relations in the boundary layer of a ferromagnetic nanofluid flow via a vertical cone surface. The study of the effect of the magnetic dipole on ferromagnetic nanofluid flow via the vertical cone surface makes this work different from the existing literature. A nonlinear ordinary differential equation replaces the governing equations and is solved using HAM techniques. Initially, Liao [25,26,27] presented a homotopy analysis method with HAM. The method has fast convergent solutions with many advantages over some existing methods. Various researchers have been drawn to it as a result of its rapid convergence [28,29,30,31]. The results collected for all associated parameters on all profiles are shown graphically. The validation of the results by comparing them to previously published material in the literature is an important feature of the presented model. In this regard, illustrious coherence has been attained.

2. Materials and Methods

An incompressible electrically conducting viscous nanofluid flow via a vertical cone with bioconvection is explored in two dimensions as an axisymmetric, steady, natural convective ferrofluid flow. Furthermore, it is postulated that temperature and concentration are non-uniform at the surface due to the influence of heat generation/absorption, chemical reaction and viscous dissipation. The flow is electrically magnetized by a magnetic dipole, and a Darcy–Forchheimer porous medium model is also used. Thermal radiation exists as a unidirectional flux in the transverse to the cone surface (s-direction). In comparison to the s-direction, the radiation heat flux in the x-direction is considered neglected. The x-axis of the chosen coordinate system corresponds to the direction of flow over the cone surface. T w is taken to be the temperature at the cone’s surface ( s = 0 ), and the concentration is governed by the condition D B C s + D T T T s = 0 at the cone’s surface, where T is the temperature and C is the concentration and N density of microorganisms in the ambient nanofluid.
The boundary layer equation [12,14,21] based on the assumptions stated above are the equations of continuity and momentum as well as energy, concentration, and microorganisms:
( r u ) x + { ( r w ) } s = 0 ,
w u s + u u x = μ f ρ f 2 u s 2 μ f k o * u ρ f C b k o * u 2 + λ o M H x + g [ β T ( T T ) + β C ( C C ) + β N ( N N ) ] c o s ( α ) ,
w T s + u T x = α f 2 T s 2 + μ f ( ρ c p ) f u s 2 + u H x + w H s λ 0 ( ρ c p ) f T M T τ D B C s T s + D T T T s 2 1 ( ρ c p ) f q r s ,
w C s + u C x = D B 2 C s 2 + D T T 2 T s 2 K r ( C C ) ,
w N s + u N x + b W c C w C C s N s + N 2 C s 2 = D n 2 N s 2 ,
with initial boundary conditions
u = 0 , w = W w , T = T w , D B C s = D T T T s , N = N w   a t   s = 0 ,
u 0 , w 0 , T T , C C , N N ,   a s   s ,
where (u,w) are the velocity components in the x-direction (radial) and s-direction (transverse), respectively; T, C, N are the temperature, concentration, and gyrotactic microorganism, respectively; the diffusion coefficients named Brownian, thermophoresis, and microorganism correspond to D B , D T , and D n , respectively; while τ is the ratio of heat capacitance, fluid density is ρ f , thermal conductivity of fluid is k f , electrical conductivity of fluid is σ , the dynamic viscosity is μ f , thermal diffusivity of base fluid is α f , magnetic permeability is λ o , heat capacitance of fluid is ( ρ c p ) f , first order chemical reaction parameter is K r , speed of gyrotactic cell is W c , and b is chemotaxis.

Magnetic Dipole

The magnetic field features impacted the ferrofluid flow with magnetic dipole effects detected mostly by magnetic scalar potential Φ 1 , as given in Equation (8):
Φ 1 = γ 2 π x x 2 + ( s + c ) 2 ,
Considering H x and H s to be the components of magnetic field, with γ as the magnetic field strength at the source, see Equations (9) and (10):
H x = Φ 1 x = γ 2 π x 2 ( s + d ) 2 [ x 2 + ( s + d ) 2 ] 2
H s = Φ 1 s = γ 2 π 2 x ( s + d ) [ x 2 + ( s + d ) 2 ] 2 .
As the strength of a magnetic body is normally approximately equal to the H x and H s gradients, it is therefore given as in (11):
H = H x 2 + H s 2 .
Equation (12) displays the approximate linearized relation of the magnetization M as function of temperature T,
M = K 1 ( T T ) ,
with K 1 identified as the ferromagnetic coefficient. Figure 1 depicts the physical configuration of the heated ferrofluid.
Considering the following transformations, given the stream function as Φ ( x , s ) , such that
u = 1 r Φ s w = 1 r Φ x
with Φ ( x , s ) = ν f r R a x f ( ζ ) and R a x is the Rayleigh number given by R a x = ρ f β T g ( T w T ) x 3 c o s ( α ) ν f 2 , therefore, the following are given:
u = ν f R a x 1 2 x f ( ζ ) , w = ν f R a x 1 4 x ζ f ( ζ ) f ( ζ ) , ζ = s x R a x 1 4 , ( T w T ) θ ( ζ ) = ( T T ) , ( C w C ) ϕ ( ζ ) = ( C C ) , ( N w N ) χ ( ζ ) = ( N N ) .
Taking r to be approximately the cone local radius, for the thermal boundary layer becoming thin, it will be along the x coordinate with r = x s i n ( α ) .
By using the above transformations, Equation (1) will be satisfactory, and Equations (2)–(5) will be
f P 1 f + f f F r f 2 + 2 β ( ζ + α 1 ) 4 ( 1 + θ ) + N c ϕ + R a b χ = 0 ,
1 + R d θ + P r f θ + N b ϕ θ + N t ( θ ) 2 + 2 P r β λ ( θ ϵ ) ( f ζ f ) ( ζ + d α 1 ) 3 + P r β λ ( θ ϵ ) 2 f ( ζ + α 1 ) 4 + 4 ( ζ f f ) ( ζ + α 1 ) 5 + P r E c ( f ) 2 = 0 ,
ϕ + N t N b θ + S c f ϕ δ S c ϕ = 0 ,
χ P e ϕ χ + ϕ χ + δ n ϕ + L b f χ = 0 .
Moreover, with the new boundary conditions:
f = 1 , f = S , θ = 1 , N b ϕ + N t θ = 0 , χ = 1 a t ζ = 0 , f 0 , θ 0 , ϕ 0 , χ 0 , a s ζ ,
where α 1 is dimensionless distance, N c is the ratio due to buoyancy force, R a b is the bioconvection Rayleigh number, β is the ferrohydrodynamic interaction parameter, ε the Curie temperature, λ is the heat dissipation parameter, S is the heat generation/absorption parameter ( S > 0 for suction and S < 0 for injection), N b and N t are the Brownian motion and thermophoresis parameters, the Prandtl number is P r , the Eckert number is E c , the radiation parameter is R d , the local inertia parameter is F r , the chemical reaction parameter is δ , the porosity parameter is P 1 , the Schmidt number is S c , the Lewis number of bioconvection is L b , the Peclet number P e , δ n is the bioconvection constant, and quantities are defined by
S = W w x ν f R a 1 2 , N c = g β C ( C w C ) x 3 c o s ( α 1 ) ν f 2 R a x , R a b = g β N ( N w N ) x 3 c o s ( α ) ν f 2 R a x , P 1 = μ f k o * , L i = C b k o * , β = γ λ o K ρ f ( T w T ) 2 π μ f 2 , P r = μ f c p k f , E c = W w 2 c p ( T w T ) , λ = μ f 2 ρ f ( T w T ) R a x 3 4 , δ = K 1 x 2 ν f R a x 1 2 , R d = 16 σ * T 3 3 k * k f , S c = ν f D B , P e = b W c D n , L b = ν f D n , ϵ = T T T w , d = a R a x x , N t = τ D T ( T w T ) T α f , N b = τ D B ( C w C ) α f , δ n = N N w N .
The local Nusselt, Sherwood, and local Density expressions, as well as the coefficient of skin friction, will be computed by
C f = 2 τ w ρ f W w 2 , N u x = q h x k f ( T w T ) , S h x = q m x D B ( C w C ) , S n x = q n x D n ( N w N ) ,
τ w = μ f u s | s = 0 , q h = [ k f T s + q r ] | s = 0 , q m = D B C s | s = 0 , q n = D n N s | s = 0 ,
R a x 1 4 C f = 2 f ( 0 ) , N u = R a x 1 4 ( 1 + R d ) θ ( 0 ) , S h = R a x 1 4 ϕ ( 0 ) , S n = R a x 1 4 χ ( 0 ) .

3. HAM Solutions Methodology

The homotopy analysis method (HAM) was applied to solve Equations (15)–(18). Shijun Liao developed this technique in 1992. It is often valid, regardless of whether there are a limited number of parameters or otherwise. It can be used to solve both weakly and strongly nonlinear problems. It offers a wide range of options for selecting the base functions of solutions, as well as discretion in choosing the linear operators. However, it provides a convenient method for ensuring the convergence of series solutions. Therefore, this method differs from other techniques, with examples like Adomain decomposition and the delta expansion methods. In the introduction section, some studies on the approach were presented.
Taking the initial guesses of the f ( ζ ) , θ ( ζ ) , ϕ ( ζ ) , and χ ( ζ ) with the auxiliary linear operators respectively as
f 0 ( ζ ) = 1 e ζ , θ 0 ( ζ ) = B i 1 + B i e ζ , ϕ 0 ( ζ ) = N t N b e ζ , χ 0 ( ζ ) = e ζ .
and
L f = f f , L θ = θ θ , L ϕ = ϕ ϕ , L χ = χ χ ,
the properties are satisfied as given below
L f ( Λ 1 + Λ 2 e ζ + Λ 3 e ζ ) = 0 , L θ ( Λ 4 e ζ + Λ 5 e ζ ) = 0 , L ϕ ( Λ 6 e ζ + Λ 7 e ζ ) = 0 , L χ ( Λ 8 e ζ + Λ 9 e ζ ) = 0
with arbitrary constants Λ i , i [ 1 , 9 ] . The Zeroth order form of the problem is given by
( 1 p ) L f [ f ( ζ ; p ) f 0 ( ζ ) ] = p h f N f [ f ( ζ , p ) , θ ( ζ , p ) , ϕ ( ζ , p ) , χ ( ζ , p ) ] ,
( 1 p ) L θ [ θ ( ζ ; p ) θ 0 ( ζ ) ] = p h θ N θ [ θ ( ζ , p ) , f ( ζ , p ) , ϕ ( ζ , p ) ] ,
( 1 p ) L ϕ [ ϕ ( ζ , p ) ϕ 0 ( ζ ) ] = p h ϕ N ϕ [ ϕ ( ζ , p ) , θ ( ζ , p ) , f ( ζ , p ) ] ,
( 1 p ) L χ [ χ ( ζ , p ) χ 0 ( ζ ) ] = p h χ N χ [ χ ( ζ , p p ) , ϕ ( ζ , p ) , f ( ζ , p ) ] ,
with p [ 0 , 1 ] as the embedded parameter, and nonlinear operators N f , N θ , N ϕ , and N χ obtained by using Equations (15)–(18).
The problems’ equivalent m order of the deformation are
L f [ f m ( ζ , p ) η m f m 1 ( ζ ) ] = h f R f , m ( ζ ) ,
L θ [ θ m ( ζ , p ) η m θ m 1 ( ζ ) ] = h θ R θ , m ( ζ ) ,
L ϕ [ ϕ m ( ζ , p ) η m ϕ m 1 ( ζ ) ] = h ϕ R ϕ , m ( ζ ) ,
L χ [ χ m ( ζ , p ) η m χ m 1 ( ζ ) ] = h χ R χ , m ( ζ ) ,
f m = S , f m = 0 , θ m B i θ m = 0 , N b ϕ m + N t θ m = 0 , χ m = 0 , a t   ζ = 0 f m = 0 , θ m = 0 , ϕ m = 0 , χ m = 0 a s ζ .
η m = 0 , if m 1 1 , if m > 1 ,
where R f m ( ζ ) , R θ m ( ζ ) , R ϕ m ( ζ ) , R χ m ( ζ ) can be obtained using Equations (15)–(18).
The general solutions are given by
f m ( ζ ) = f m s ( ζ ) + Λ 1 + Λ 2 e ζ + Λ 3 e ζ ,
θ m ( ζ ) = θ m s ( ζ ) + Λ 4 e ζ + Λ 5 e ζ ,
ϕ m ( ζ ) = ϕ m s ( ζ ) + Λ 6 e ζ + Λ 7 e ζ ,
χ m ( ζ ) = χ m s ( ζ ) + Λ 8 e ζ + Λ 9 e ζ ,
where ( f m s ( ζ ) , θ m s ( ζ ) , ϕ m s ( ζ ) , χ m s ( ζ ) ) are special solutions.

4. Analysis of Convergence of the Solutions

The convergence Table 1 is organized for each profile up towards the 35th order of approximation. Table 2 compares the current work to the published work and reveals that there is very close agreement.

5. Results and Discussion

The role of β on velocity profile is detected in Figure 2. The ferromagnetic parameter emphasizes the effect of the magnetic dipole’s external magnetic field on fluid dynamics. As the magnetic field acts as a deforming force, the axial velocity decreases. Figure 3 depicts the important properties of the F r on the dimensionless velocity. Obviously, increasing F r causes the velocity of fluid layers to decrease. From a physical standpoint, the local inertia parameter generates resistance forces against the motion of fluid particles. Moreover, as the local inertia parameter increases, so does the velocity. Figure 4 depicts the impact of increasing P 1 on the velocity. The porosity parameter is defined as the kinematic viscosity to permeability strength of porous space ratio. As the porosity parameter increases, the velocity curves definitely decrease. Figure 5 depicts the effect of the R a b on the velocity profile. It was discovered that when the values of R a b rise, the velocity accelerates. As a result of the buoyancy forces caused by bioconvection, the fluid velocity increases by increasing the bioconvection Raleigh number. Figure 6 shows that the temperature on the boundary layers enhances as a result of increasing the values of β . This is due to an interaction of the fluid’s movement and the interference of the ferromagnetic particles. The interplay reduces the velocity while frictional heating increases between the fluid layers, resulting in an increase in the thermal boundary layer thickness. The effect of λ on the temperature profile is depicted in Figure 7. In this case, temperature is displayed as an increasing function of λ . Usually, as the values of λ increase, so does the thermal conductivity, and thus the temperature. Figure 8 displays the role of F r on temperature profile. A rise in F r results in a rise in temperature and the thermal boundary layer thickness. Figure 9 reveals that as B i increases, so does the thickness of the thermal boundary layer, and the temperature also enhances. A higher Biot number contributes to more convection, which leads to the enhancement of the temperature and thermal boundary layer thickness. Figure 10 shows the effect of S < 0 on the dimensionless temperature. This figure reveals that when increasing S < 0 , the temperature and the thickness of the thermal boundary layers both decrease. As an outcome, suction is removed from the warm fluid in the boundary layer region to a large extent. Moreover, the opposite trend is observed in the temperature profile, as shown in Figure 11 with S > 0 . This is attributable to the fact that the temperature of the fluid is raised by injecting warm fluid into the boundary layer region. Figure 12 depicts the effects of R d on the temperature profile for various R d values. An increase in the R d causes a rise in the temperature, and the effect of thermal radiation improves the medium’s thermal diffusive. Besides that, for higher N t values, temperatures rise in the boundary layer region (Figure 13). The thermophoretic force developed in the boundary layer regime is an outcome of the temperature gradient; such forces entail the diffusion of nanoparticles out of the higher temperature area to a lower temperature area, leading to a thermal boundary layer thickness enhancement. Figure 14 depicts the N b characteristics on a temperature profile. Usually, a rise in N b improves the motion of fluid particles randomly, resulting in more heat generation. As a result, the temperature rises. Figure 15 depicts the temperature distribution by raising E c values. The E c defines the connection among the flow of kinetic energy and heat enthalpy variation. As a consequence, raising E c also raises the kinetic energy. Moreover, temperature is well understood to be defined as average kinetic energy. As a result, the fluid’s temperature rises. This graph shows that as E c increases, so does the temperature. Figure 16 shows the role of P r on the temperature profile. It has been discovered that raising the P r lowers the temperature of the fluid flow. Large P r values clearly result in the thinning of thermal boundary layers. As δ increases, so does the concentration profile. The concentration becomes less effective as the delta values increase, as shown in Figure 17. The physical effect of N b on a concentration profile is depicted in Figure 18. Brownian motion does play a role in determining the efficiency of heat transfer during nanofluid flow. The nanoparticles collide with one another and transfer energy due to the random motion of a nanofluid. As a result, as N b levels rise, the concentration profile falls. The effect of N t on nanoparticle concentration is depicted in Figure 19. The concentration field rises in this case due to an increase in N t . Larger N t causes an increase in thermophoresis forces, which further frequently carries nanomaterials from higher to lower temperature regions. As a result, the concentration decreases. The effect of S c on concentration is depicted in Figure 20. Sc denotes the momentum-to-mass diffusivity ratio, which measures the relative efficacy of momentum and mass transport through diffusion within concentration boundary layers. Figure 21 and Figure 22 show the effects of P e and L b on the microorganism profile. According to these figures, the microorganism field decreases with increase of both numbers. According to Table 3, as P r estimates increase, so does the Nusselt number. Table 4 shows that as δ and S c increase, so does the Sherwood number. The results of motile microorganism density are increased by increasing L b and R a b , as shown in Table 5.

6. Conclusions

A vertical cone has been used to study the flow of ferromagnetic nanofluid with bioconvection and magnetic dipole elements. There is also evidence of the Darcy–Forchheimer flow model. Viscosity dissipation, Brownian motion, chemical reaction, and thermophoresis are involved. Using appropriate transformations, nonlinear PDEs were reduced to a set of nonlinear ODEs. The important outcomes are listed below:
  • The velocity decreases as the ferromagnetic interaction parameter, porosity parameter, and local inertia parameter increase.
  • The temperature rises as the ferromagnetic interaction, heat dissipation, injection, thermal radiation parameters and Eckert number are raised, and reduces when the Prandtl number and sunction parameter are raised.
  • When the Brownian motion, chemical reaction parameters and Schmidt number increase, the concentration decreases, while it increases when the thermophoresis parameter increases.
  • The motile microorganism density decreases through raising the Peclet number and bioconvection Lewis number.

Author Contributions

A.H.U. and W.K. formulated and solved the problem. A.H.U. and Z.S. wrote the draft and the revised manuscript. A.H.U. and Z.S. sketched the graphs. U.W.H. and P.K. performed the investigations, supervision, and funding. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Availability exists for entirety of the data.

Acknowledgments

This research project is supported by Thailand Science Research and Innovation (TSRI) Basic Research Fund: Fiscal year 2022. The authors are thankful to KMUTT for providing the technical and financial support. The first author is supported by the Petchra Pra Jom Klao Doctoral Scholarship Academic for Ph.D. program of King Mongkut’s University of Technology Thonburi (KMUTT) [Grant No. 13/2562].

Conflicts of Interest

The authors have no conflict of interests.

References

  1. Raj, K.; Moskowitz, R. Commercial applications of ferrofluids. J. Magn. Magn. Mater. 1990, 85, 233–245. [Google Scholar] [CrossRef]
  2. Goldman, A. Handbook of Modern Ferromagnetic Materials; Springer Science and Business Media: Berlin, Germany, 2012. [Google Scholar]
  3. Hathway, D.B. dB-sound. Eng. Mag. 1979, 13, 42–44. [Google Scholar]
  4. Andersson, H.I.; Valnes, O.A. Flow of a heated ferrofluid over a stretching sheet in the presence of a magnetic dipole. Acta Mech. 1998, 128, 39–47. [Google Scholar] [CrossRef]
  5. Zeeshan, A.; Majeed, A.; Ellahi, R.; Zia, Q.M.Z. Mixed convection flow and heat transfer in ferromagnetic fluid over a stretching sheet with partial slip effects. Therm. Sci. 2018, 22, 2515–2526. [Google Scholar] [CrossRef]
  6. Hayat, T.; Ahmad, S.; Khan, M.I.; Alsaedi, A. Exploring magnetic dipole contribution on radiative flow of ferromagnetic Williamson fluid. Results Phys. 2018, 8, 545–551. [Google Scholar] [CrossRef]
  7. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Argonne National Lab: Lemont, IL, USA, 1995.
  8. Ellahi, R. The effects of MHD and temperature dependent viscosity on the flow of non-Newtonian nanofluid in a pipe: Analytical solutions. Appl. Math. Model. 2013, 37, 1451–1467. [Google Scholar] [CrossRef]
  9. Ellahi, R.; Raza, M.; Akbar, N.S. Study of peristaltic flow of nanofluid with entropy generation in a porous medium. J. Porous Media 2017, 20, 461–478. [Google Scholar] [CrossRef]
  10. Hayat, T.; Ahmad, S.; Khan, M.I.; Alsaedi, A. Modeling and analyzing flow of third grade nanofluid due to rotating stretchable disk with chemical reaction and heat source. Phys. B Condens. Matter. 2018, 537, 116–126. [Google Scholar] [CrossRef]
  11. Awais, M.; Shah, Z.; Perveen, N.; Ali, A.; Kumam, P.; Thounthong, P. MHD effects on ciliary-induced peristaltic flow coatings with rheological hybrid nanofluid. Coatings 2020, 10, 186. [Google Scholar] [CrossRef] [Green Version]
  12. Reddy, P.S.; Sreedevi, P.; Chamkha, A.J. MHD natural convection boundary layer flow of nanofluid over a vertical cone with chemical reaction and suction/injection. Comput. Therm. Sci. Int. J. 2017, 9, 165–182. [Google Scholar] [CrossRef]
  13. Alsaedi, A.; Khan, M.I.; Farooq, M.; Gull, N.; Hayat, T. Magnetohydrodynamic (MHD) stratified bioconvective flow of nanofluid due to gyrotactic microorganisms. Adv. Powder Technol. 2017, 28, 288–298. [Google Scholar] [CrossRef]
  14. Hayat, T.; Bashir, Z.; Qayyum, S.; Alsaedi, A. Nonlinear radiative flow of nanofluid in presence of gyrotactic microorganisms and magnetohydrodynamic. Int. J. Numer. Methods Heat Fluid Flow 2019, 29, 3039–3055. [Google Scholar] [CrossRef]
  15. Nadeem, S.; Ijaz, M.; El-Kott, A.; Ayub, M. Rosseland analysis for ferromagnetic fluid in presence of gyrotactic microorganisms and magnetic dipole. Ain Shams Eng. J. 2020, 11, 1295–1308. [Google Scholar] [CrossRef]
  16. Bhatti, M.M.; Michaelides, E.E. Study of Arrhenius activation energy on the thermo-bioconvection nanofluid flow over a Riga plate. J. Therm. Anal. Calorim. 2021, 143, 2029–2038. [Google Scholar] [CrossRef]
  17. Waqas, H.; Khan, S.A.; Alghamdi, M.; Alqarni, M.S.; Muhammad, T. Numerical simulation for bio-convection flow of magnetized non-Newtonian nanofluid due to stretching cylinder/plate with swimming motile microorganisms. Eur. Phys. J. Spec. Top. 2021, 230, 1239–1256. [Google Scholar] [CrossRef]
  18. Ingham, D.B.; Pop, I. (Eds.) Transport Phenomena in Porous Media III; Elsevier: Oxford, UK, 2005; Volume 3. [Google Scholar]
  19. Zeng, Z.; Grigg, R. A criterion for non-Darcy flow in porous media. Transp. Porous Media 2006, 63, 57–69. [Google Scholar] [CrossRef]
  20. Kumar, B.R.; Sivaraj, R. MHD viscoelastic fluid non-Darcy flow over a vertical cone and a flat plate. Int. Commun. Heat Mass Transf. 2013, 40, 1–6. [Google Scholar] [CrossRef]
  21. Chamkha, A.; Abbasbandy, S.; Rashad, A.M. Non-Darcy natural convection flow for non-Newtonian nanofluid over cone saturated in porous medium with uniform heat and volume fraction fluxes. Int. J. Numer. Method Heat Fluid Flow 2015, 25, 422–437. [Google Scholar] [CrossRef]
  22. Mallikarjuna, B.; Rashad, A.M.; Hussein, A.K.; Raju, S.H. Transpiration and thermophoresis effects on Non-Darcy convective flow past a rotating cone with thermal radiation. Arab. J. Sci. Eng. 2016, 41, 4691–4700. [Google Scholar] [CrossRef]
  23. Durairaj, M.; Ramachandran, S.; Mehdi, R.M. Heat generating/absorbing and chemically reacting Casson fluid flow over a vertical cone and flat plate saturated with non-Darcy porous medium. Int. J. Numer. Method Heat Fluid Flow 2017, 27, 156–173. [Google Scholar] [CrossRef]
  24. Pătrulescu, F.O.; Groşan, T.; Pop, I. Natural convection from a vertical plate embedded in a non-Darcy bidisperse porous medium. J. Heat Transf. 2020, 142, 012504. [Google Scholar] [CrossRef]
  25. Liao, S.J. The Proposed Homotopy Analysis Method for the Solution of Nonlinear Problems. Ph.D. Thesis, Shanghai Jiao Tong University, Shanghai, China, 1992. [Google Scholar]
  26. Liao, S.J. An explicit, totally analytic approximate solution for Blasius viscous flow problems. Int. J. Non-Linear Mech. 1999, 34, 759–778. [Google Scholar] [CrossRef]
  27. Liao, S. Beyond Perturbation: Introduction to the Homotopy Analysis Method; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  28. Irfan, M.; Khan, M.; Khan, W.A.; Sajid, M. Thermal and solutal stratifications in flow of Oldroyd-B nanofluid with variable conductivity. Appl. Phys. A 2018, 124, 674. [Google Scholar] [CrossRef]
  29. Tlili, I.; Waqas, H.; Almaneea, A.; Khan, S.U.; Imran, M. Activation energy and second order slip in bioconvection of Oldroyd-B nanofluid over a stretching cylinder: A proposed mathematical model. Processes 2019, 7, 914. [Google Scholar] [CrossRef] [Green Version]
  30. Usman, A.H.; Khan, N.S.; Humphries, U.W.; Ullah, Z.; Shah, Q.; Kumam, P.; Thounthong, P.; Khan, W.; Kaewkhao, A.; Bhaumik, A. Computational optimization for the deposition of bioconvection thin Oldroyd-B nanofluid with entropy generation. Sci. Rep. 2021, 11, 11641. [Google Scholar] [CrossRef] [PubMed]
  31. Hayat, T.; Waqas, M.; Shehzad, S.A.; Alsaedi, A. Mixed convection flow of viscoelastic nanofluid by a cylinder with variable thermal conductivity and heat source/sink. Int. J. Numer. Method Heat Fluid Flow 2016, 26, 214–234. [Google Scholar] [CrossRef]
Figure 1. A picture scheme of the problem.
Figure 1. A picture scheme of the problem.
Coatings 11 01129 g001
Figure 2. Influence of β on f ( ζ ) .
Figure 2. Influence of β on f ( ζ ) .
Coatings 11 01129 g002
Figure 3. Influence of F r on f ( ζ ) .
Figure 3. Influence of F r on f ( ζ ) .
Coatings 11 01129 g003
Figure 4. Influence of P 1 on f ( ζ ) .
Figure 4. Influence of P 1 on f ( ζ ) .
Coatings 11 01129 g004
Figure 5. Influence of R a b on f ( ζ ) .
Figure 5. Influence of R a b on f ( ζ ) .
Coatings 11 01129 g005
Figure 6. Influence of β on θ ( ζ ) .
Figure 6. Influence of β on θ ( ζ ) .
Coatings 11 01129 g006
Figure 7. Influence of λ on θ ( ζ ) .
Figure 7. Influence of λ on θ ( ζ ) .
Coatings 11 01129 g007
Figure 8. Influence of F r on θ ( ζ ) .
Figure 8. Influence of F r on θ ( ζ ) .
Coatings 11 01129 g008
Figure 9. Influence of B i on θ ( ζ ) .
Figure 9. Influence of B i on θ ( ζ ) .
Coatings 11 01129 g009
Figure 10. Influence of S < 0 on θ ( ζ ) .
Figure 10. Influence of S < 0 on θ ( ζ ) .
Coatings 11 01129 g010
Figure 11. Influence of S > 0 on θ ( ζ ) .
Figure 11. Influence of S > 0 on θ ( ζ ) .
Coatings 11 01129 g011
Figure 12. Influence of R d on θ ( ζ ) .
Figure 12. Influence of R d on θ ( ζ ) .
Coatings 11 01129 g012
Figure 13. Influence of N t on θ ( ζ ) .
Figure 13. Influence of N t on θ ( ζ ) .
Coatings 11 01129 g013
Figure 14. Influence of N b on θ ( ζ ) .
Figure 14. Influence of N b on θ ( ζ ) .
Coatings 11 01129 g014
Figure 15. Influence of E c on θ ( ζ ) .
Figure 15. Influence of E c on θ ( ζ ) .
Coatings 11 01129 g015
Figure 16. Influence of P r on θ ( ζ ) .
Figure 16. Influence of P r on θ ( ζ ) .
Coatings 11 01129 g016
Figure 17. Influence of δ on ϕ ( ζ ) .
Figure 17. Influence of δ on ϕ ( ζ ) .
Coatings 11 01129 g017
Figure 18. Influence of N b on ϕ ( ζ ) .
Figure 18. Influence of N b on ϕ ( ζ ) .
Coatings 11 01129 g018
Figure 19. Influence of N t on ϕ ( ζ ) .
Figure 19. Influence of N t on ϕ ( ζ ) .
Coatings 11 01129 g019
Figure 20. Influence of S c on ϕ ( ζ ) .
Figure 20. Influence of S c on ϕ ( ζ ) .
Coatings 11 01129 g020
Figure 21. Influence of P e on χ ( ζ ) .
Figure 21. Influence of P e on χ ( ζ ) .
Coatings 11 01129 g021
Figure 22. Influence of L b on χ ( ζ ) .
Figure 22. Influence of L b on χ ( ζ ) .
Coatings 11 01129 g022
Table 1. Convergence of HAM solutions with order approximations.
Table 1. Convergence of HAM solutions with order approximations.
Order of Approximations f ( 0 ) θ ( 0 ) ϕ ( 0 ) χ ( 0 )
11.876330.63331.122220.69596
51.935210.767481.123270.69583
101.933850.768041.123490.69575
151.931350.76581.123540.69575
251.930260.76581.123540.69575
301.930260.76581.123540.69575
351.930260.76581.123540.69575
Table 2. Comparison of θ ( 0 ) and ϕ ( 0 ) by alternating buoyancy ratio parameter N c with published work.
Table 2. Comparison of θ ( 0 ) and ϕ ( 0 ) by alternating buoyancy ratio parameter N c with published work.
N c θ ( 0 ) Present ϕ ( 0 ) Present
Reddy et al. [12]Reddy et al. [12]
0.10.325980.326011.483941.48381
0.20.324050.324111.467891.46801
0.30.322290.322311.452141.45215
0.40.321250.321291.435981.43594
0.50.318680.318671.419381.41940
Table 3. Influence of β , ε , λ , S , R d , F r , E c , and P r on θ ( 0 ) .
Table 3. Influence of β , ε , λ , S , R d , F r , E c , and P r on θ ( 0 ) .
β ε λ S Rd Fr Ec Pr θ ( 0 )
0.20.30.20.40.30.20.10.70.44701
0.3 0.43167
0.5 0.40752
0.4 0.39766
0.5 0.38914
0.6 0.37245
0.5 0.45806
0.8 0.45739
1.1 0.45623
0.8 1.03123
1.2 0.87433
1.6 0.64876
0.5 0.45342
0.7 0.44534
0.9 0.43998
0.6 0.64554
1.0 0.63854
1.4 0.61291
0.3 0.38453
0.5 0.38123
0.7 0.37941
6.70.75651
7.70.80612
8.70.86432
Table 4. Influence of δ , S c , N b on ϕ ( 0 ) .
Table 4. Influence of δ , S c , N b on ϕ ( 0 ) .
δ Sc Nb ϕ ( 0 )
0.20.30.51.28733
0.4 1.26931
0.6 1.24887
0.5 1.32742
0.7 1.31022
0.9 1.30271
1.01.33075
1.51.41186
2.01.47572
Table 5. Influence of N δ , P e , L b , R a b on χ ( 0 ) .
Table 5. Influence of N δ , P e , L b , R a b on χ ( 0 ) .
δ n PeLb Ra b χ ( 0 )
0.10.30.30.20.57238
0.2 0.56103
0.3 0.55271
0.5 0.42714
1 0.45102
1.5 0.49327
0.5 0.53185
1 0.62386
1.5 0.66671
0.30.75408
0.40.76965
0.50.77121
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Usman, A.H.; Shah, Z.; Kumam, P.; Khan, W.; Humphries, U.W. Nanomechanical Concepts in Magnetically Guided Systems to Investigate the Magnetic Dipole Effect on Ferromagnetic Flow Past a Vertical Cone Surface. Coatings 2021, 11, 1129. https://doi.org/10.3390/coatings11091129

AMA Style

Usman AH, Shah Z, Kumam P, Khan W, Humphries UW. Nanomechanical Concepts in Magnetically Guided Systems to Investigate the Magnetic Dipole Effect on Ferromagnetic Flow Past a Vertical Cone Surface. Coatings. 2021; 11(9):1129. https://doi.org/10.3390/coatings11091129

Chicago/Turabian Style

Usman, Auwalu Hamisu, Zahir Shah, Poom Kumam, Waris Khan, and Usa Wannasingha Humphries. 2021. "Nanomechanical Concepts in Magnetically Guided Systems to Investigate the Magnetic Dipole Effect on Ferromagnetic Flow Past a Vertical Cone Surface" Coatings 11, no. 9: 1129. https://doi.org/10.3390/coatings11091129

APA Style

Usman, A. H., Shah, Z., Kumam, P., Khan, W., & Humphries, U. W. (2021). Nanomechanical Concepts in Magnetically Guided Systems to Investigate the Magnetic Dipole Effect on Ferromagnetic Flow Past a Vertical Cone Surface. Coatings, 11(9), 1129. https://doi.org/10.3390/coatings11091129

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