# Impact of Ferromagnetic Nanoparticles Submerged in Chemically Reactive Viscoelastic Fluid Transport Influenced by Double Magnetic Dipole

Authors: N. Kousar; Taher A. Nofal; W. Tahir; S. M. Bilal; Ndolane Sene
Journal: Journal of Nanomaterials
 (2022)

---

## Abstract

The present article summarized the effects of double magnetic dipole for chemically reactive viscoelastic fluid in the presence of two different ferromagnetic particles, namely, nickel zinc ferriteNiZnFe2O4 and magnetite ferrite (Fe3O4). Due to double magnetic dipole, an external magnetic field is applied normal to the flow. Blood is used as base fluid due to its viscoelastic fluid properties. The Cattaneo-Christov heat flux model is used for heat transport phenomena. The physical model is formulated in the form of partial differential equations which are then converted into ordinary differential equations using the suitable transformations. The system is solved numerically using shooting method along with Runge-Kutta-Fehlberg method. The characteristics of different parameters like the strength of homogeneous-heterogeneous reactions (k1andk2), ferrohydrodynamic interaction (β1), Schmidt number (Sc), Deborah number (α1a), and thermal relaxation time (α1c) on velocity, temperature, and concentration profiles are analyzed through graphs and in tabular form. It has been observed that as magnetic dipole creates a force which attracts the ferrite particles, hence, it slows down the velocity profile. Concentration field depresses due to the presence of strength of heterogeneous reaction parameter k2.It is also noted that by expanding values of thermal relaxation time (α1c), the temperature profile shows a reverse behavior.


---

## Body

## 1. Introduction

Ferrofluids lie in the category of smart materials those consist of micron-sized colloidal magnetic nanoparticles that are saturated in a nonmagnetic base fluid. The most fascinating feature of these fluids is its highly magnetizing ability when an external magnetic field is applied. Pappell [1] initially highlighted the characteristic of ferrofluid in 1963. He utilized ferrofluid in a weightless atmosphere in terms of liquid rocket fuel, which is pinched near a pump inlet by applying an external magnetic field. Ferrofluids have many appealing applications in electrical instruments such as in hard disk, rotating X-ray tubes, shafts, and rods. The role of ferrofluids in biomedical equipment’s is no doubt incredible, which is helpful in the process of wound treatment, asthma treatment, removal of cancer with hyperthermia, and many more. For these kinds of procedures, the homogeneous and heterogeneous have great importance. The homogenous and heterogeneous reactions take places in the bulk and occur on the catalyst surface, respectively. Whether the process is homogeneous or heterogeneous depends upon the situation that they exist in the absolute majority of the fluid or a part of catalytic surfaces. In essence, the homogeneous process is continuously intact within the given phase, whereas heterogeneous reactions have restricted boundaries. There are also a variety of chemically reacting structures that involve both h-h reactions termed as such as catalysis, burning, and biochemical reacting systems. In this flow, due to an external flow applied on the outer region of boundary layer, there exists a reaction in that region as explained by [2, 3]. Initially, Andersson and Valnes [4] studied the viscous ferrofluid with magnetic dipole effects on the stretching surface. They examined the influence of magnetothermomechanical coupling on the fluid. Numerical studies related to ferrofluid inside a channel with cold walls focusing on a line source dipole delineated by Ganguly et al. [5]. In a porous medium, Sharma et al. [6] specified the attributes of convection by focusing dust particles on a ferromagnetic fluid. The results corresponding to characteristics of heat transfer in a Darcy number ferrofluid flow with porous wall is analyzed by Strek [7]. The dipole effects in a ferrofluid flow origin below the channel with isothermal wall were studied by Strek and Jopek [8]. They focused on time dependent heat transfer and analyzed the spatial alteration creating magnetization because of gradient effects, depending temperature magnetic susceptibility for ferrofluids. Sadiq et al. [9] studied the Casson fluid model to study the thermal performance under the Brownian and thermoheretic effects of hybrid nanoparticles. The dipole’s magnetization in the series flow of ferrofluid taking elongated surface highlighting thermal radiation is explored by Makinde and Aziz [10]. Aminfar et al. [11] implement a non-acting magnetic field inside a vertical tube during the numerical studies of mixed convection of ferrofluid. During his studies, he used both + and - types of magnetic gradients; the (-) gradient of magnetic field reacts the same as the buoyancy forces and boosts the Nusselt number, whereas the magnetic field with a (+) gradient decelerates it. A stagnation point flow along elongated sheet in the presence of heterogeneous-homogeneous reactions that is analyzed by Bachok et al. [12] analyzed a chemically reactive stagnation point flow over the elongated sheet. In his paper, he discussed the case of fluid having less kinetic viscosity in which he observed that when the extending velocity is low as compared to free steam velocity, a boundary layer is achieved, whereas the inverted boundary layer attains in the case when stretching the velocity exceeds the free steam velocity. Affixing 2-phase mixture model and effecting control volume method, Aminfar et al. [13] executed the transversal non acting magnetic field accomplishing electric current during the inspection of flowing ferrofluid inside a duct. They showed the flourishing behavior of average heat transfer coefficient. The shrinking surface is taken by Kameswaran et al. [14] during the investigation of homogeneous-heterogeneous reaction in a porous medium. By taking into account the radiation effect in an extended sheet, Titus and Abraham [15] used ferrofluid. It also concludes that due to magnetic field local vortex vary the advection energy transport also boosting heat transfer ability. Sheikholeslami et al. [16] conferred the aspects of non-uniform magneto hydrodynamic flow of ferroliquid utilizing convective heat transport. Hayat et al. [17] also detailed the effects of MHD nanofluid with heterogeneous-homogeneous reaction and velocity slip condition. It is concluded that in the case of concentration, both homogeneous and heterogeneous parameters show a reverse behavior. Imtiaz et al. [18] reflected the impact of homogeneous and heterogeneous reactions in the examination of MHD flow in a curved stretchable surface. He depicted the increasing behavior between curvature parameter and fluid velocity. Due to many applications, the boundary layer flow of non-Newtonian fluids with different effects has added an enormous attraction in the recent years [19–33].The ultimate goal of this research is to scrutinize the effects of magnetic double dipole in the examination of Maxwell ferrofluid which is composed of base fluid blood, whereas ferrite particles were used asNiZnFe2O4andFe3O4 under the high lightening impact of homogenous and heterogeneous reaction. Also, the concentration levels of these particles are about 20% into the base fluid. The basic theme of including these ferrite particles is to enhance the strength of heat transferring phenomena which is the inspiration of this research. Brief literature survey is summarized in Section 1. Section 2 is focused on the mathematical formulation of problem. In Section 3, the computational procedure is discussed in details. The possible outcomes of the current study are deliberated in Section 4. Section 5 presents the key features of this article.


## 2. Mathematical Formulation

Consider a steady 2D flow of an incompressible non-Newtonian Maxwell electrically conducting a ferromagnetic fluid running on a flat surface inx direction as shown in Figure 1. The location of dipole is set in theydirection as the displacement between the dipole and surface is taken asd, while the magnetic fieldHgenerated with the presence of magnetic dipole directions asx-axis. The scalar potential of permanent magnetic dipole which influences ferrofluid defined in [2] is
(1)δa=−I02πTan−1y+dx+Tan−1y+dx,where I0 is the representation of dipole moment per unit length. The relation between gradient of magnetic scalar potential δa and applied magnetic field is held as H1=−∇δa. The components of H1 are
(2)Hx=‐∂δa∂x=‐I02πy+dx2+y+d2+y−dx2+y−d2,Hy=‐∂δa∂y=‐I02πxx2+y−d2+xx2+y+d2.Figure 1
Geometry of flow model.The force field strength has a resultant magnitudeH, with their element forms that are termed as
(3)H=∂δa∂x2+∂δa∂y212,(4)∇Hx=∂δa/∂x∂2δa/∂x2+∂δa/∂y∂2δa/∂x∂y∂δa/∂x2+∂δa/∂y21/2,(5)∇Hy=∂δa/∂x∂2δa/∂x∂y+∂δa/∂y∂2δa/∂y2∂δa/∂x2+∂δa/∂y21/2.Experiencing Equation (1) in Equations (4)-(5), we get the following form by expanding up to order x2:
(6)∇Hy=0.The assuming effects of surface’s wall is∂δa/∂xy=0=∂2δa/∂y2y=0=0, with including assumption that x≫d We get the transformed form of equation (4)
(7)∇Hx=‐I0π1x2.The change in magnetizationM can be expressed as
(8)M=K1Tc‐T.


### 2.1. Analysis of Flow

We have considered a viscoelastic fluid flow over a flat surface with double magnetic dipoles, which are placing at a space (d) from the wall and perpendicular to the surface. Consider l as the length of plate l>>d with the wall. The temperature changes linearly with the plate length that is defined as Tw=Tc1−x/l. In the boundary layer flow, the connection between homogeneous and heterogeneous reactions adding 2-species chemically, named, AandB, are taken following Chaudhary and Merkin [2] as
(9)A+2B⟶3B,rate=kcab2.A⟶B,rate=ksaConcentrationsAandB are represented by aandb, and kii=c,s are termed as rate constants. The process is considered to be isothermal for both reactions.The corresponding flow equations in the presence of magnetic dipole are(10)∂u∂x+∂v∂y=0,u∂u∂x+v∂u∂y+λau2∂2u∂x2+v2∂2u∂y2+2uv=μ0ρnfM∂H∂x+μnfρnf∂2u∂y2,u∂T∂x+v∂T∂y+λcu∂u∂x∂T∂x+v∂v∂y∂T∂y+u∂v∂x∂T∂y+v∂u∂y∂T∂x+2uv∂2T∂x∂y+u2∂2T∂x2+v2∂T2∂y2=knfρcpnf∂2T∂y2,u∂a∂x+v∂a∂y=DA∂2a∂y2‐kcab2,u∂b∂x+v∂b∂y=DB∂2a∂y2+kcab2.The admissible flow conditions are(11)uy=0=0,vy=0=0,Ty=0=Tc1−xl,DA∂a∂yy⟶0=ksa0,DB∂a∂yy⟶0=−ksa0,uy⟶∞=u0,Ty⟶∞=Tc,ay⟶∞⟶a0,by⟶∞⟶0.where (u,v) is the velocity components alongxandydirections, (μ0) is the magnetic permeability, (μnf) is the dynamic viscosity, (vnf) is the kinematic viscosity of nanofluid, (ρnf) is the nanofluid density,ρcpnfis the specific heat, (knf) represented as thermal conductivity of the nanofluid, (λa) is the relaxation time of fluid, and (λc) identifies as the relaxation time of heat flux.


### 2.2. Thermophysical Properties

Using similarity transform as defined by [23],
(12)η=u0vx12y,T=Tc1−xlθη,u=u0f′η,v=vu0x12ηf′η2−fη2,a=a0gη,b=a0hη.By implementing similarity variables express in Eq. (12), our modified equations take the subsequent form:
(13)11−ϕ2.51−ϕ+ϕρs/ρf−α1af24f‴+ff″2−α1a4f′2f″η+2ff′f″−β1θ1−ϕ+ϕρs/ρf=0,(14)knf/kf1−ϕ+ϕρCps/ρCpf−α1cPrf24θ1″+Prfθ′2−f′θ+α1c42ff″+2ff′θ′=0,(15)1Scg″+fg′−k1gh2=0,(16)δSch″+fh′+k1gh2=0,with nondimensional boundary conditions:
(17)fη=f′η=0,θη=1,g′η=k2gη,δh′η=−k2gηatη=0,f′η⟶1,θη⟶0,gη=1,hη=0asη⟶∞.In Eqs. (13)-(16), β1 is the ferrohydrodynamic interaction, Pr defines the Prandtl number, α1a mean the Deborah number, α1c is represented here as a nondimensional thermal relaxation time, Sc the shows Schmidt number, and k1andk2 describe the strength of homogeneous and heterogeneous reaction parameters, which are defined here as
(18)β1=I0μ0KTcπρlμ02,Pr=ρcpνk,α1a=λ1au0x−1,α1c=λ1cu0xSc=νfDA,k1=kca02S,k2=ksDAνfS.DAandDB are equivalent, i.e., δ=1; then, we can write as
(19)gη+hη=1.Utilizing Eqs. (15)-(16), we get the equation as
(20)1Scg″+fg′−k1g1−g2=0.The corresponding boundary equation of the concentration field yields the following form:(21)g′η=k2gη,η⟶0,gη=1,η⟶∞.The conversions of the wall drag and convective heat transfer into a nondimensional form are scaled as(22)Cf=−2τwρnfu02,Nu=xknfkfTc−Tw∂T∂yy=0.The wall shear stress is termed as(23)τw=μnf∂u∂yy=0,qw=−knf∂T∂yy=0.The dimensionless expressions of Eq. (22) are attained:
(24)12Re12xCf=11−ϕ−2.5f′′0,Re−12xNux=−knfkfθ1′0,where,
(25)Rex=ρnfu0xμnf.


## 2.1. Analysis of Flow

We have considered a viscoelastic fluid flow over a flat surface with double magnetic dipoles, which are placing at a space (d) from the wall and perpendicular to the surface. Consider l as the length of plate l>>d with the wall. The temperature changes linearly with the plate length that is defined as Tw=Tc1−x/l. In the boundary layer flow, the connection between homogeneous and heterogeneous reactions adding 2-species chemically, named, AandB, are taken following Chaudhary and Merkin [2] as
(9)A+2B⟶3B,rate=kcab2.A⟶B,rate=ksaConcentrationsAandB are represented by aandb, and kii=c,s are termed as rate constants. The process is considered to be isothermal for both reactions.The corresponding flow equations in the presence of magnetic dipole are(10)∂u∂x+∂v∂y=0,u∂u∂x+v∂u∂y+λau2∂2u∂x2+v2∂2u∂y2+2uv=μ0ρnfM∂H∂x+μnfρnf∂2u∂y2,u∂T∂x+v∂T∂y+λcu∂u∂x∂T∂x+v∂v∂y∂T∂y+u∂v∂x∂T∂y+v∂u∂y∂T∂x+2uv∂2T∂x∂y+u2∂2T∂x2+v2∂T2∂y2=knfρcpnf∂2T∂y2,u∂a∂x+v∂a∂y=DA∂2a∂y2‐kcab2,u∂b∂x+v∂b∂y=DB∂2a∂y2+kcab2.The admissible flow conditions are(11)uy=0=0,vy=0=0,Ty=0=Tc1−xl,DA∂a∂yy⟶0=ksa0,DB∂a∂yy⟶0=−ksa0,uy⟶∞=u0,Ty⟶∞=Tc,ay⟶∞⟶a0,by⟶∞⟶0.where (u,v) is the velocity components alongxandydirections, (μ0) is the magnetic permeability, (μnf) is the dynamic viscosity, (vnf) is the kinematic viscosity of nanofluid, (ρnf) is the nanofluid density,ρcpnfis the specific heat, (knf) represented as thermal conductivity of the nanofluid, (λa) is the relaxation time of fluid, and (λc) identifies as the relaxation time of heat flux.


## 2.2. Thermophysical Properties

Using similarity transform as defined by [23],
(12)η=u0vx12y,T=Tc1−xlθη,u=u0f′η,v=vu0x12ηf′η2−fη2,a=a0gη,b=a0hη.By implementing similarity variables express in Eq. (12), our modified equations take the subsequent form:
(13)11−ϕ2.51−ϕ+ϕρs/ρf−α1af24f‴+ff″2−α1a4f′2f″η+2ff′f″−β1θ1−ϕ+ϕρs/ρf=0,(14)knf/kf1−ϕ+ϕρCps/ρCpf−α1cPrf24θ1″+Prfθ′2−f′θ+α1c42ff″+2ff′θ′=0,(15)1Scg″+fg′−k1gh2=0,(16)δSch″+fh′+k1gh2=0,with nondimensional boundary conditions:
(17)fη=f′η=0,θη=1,g′η=k2gη,δh′η=−k2gηatη=0,f′η⟶1,θη⟶0,gη=1,hη=0asη⟶∞.In Eqs. (13)-(16), β1 is the ferrohydrodynamic interaction, Pr defines the Prandtl number, α1a mean the Deborah number, α1c is represented here as a nondimensional thermal relaxation time, Sc the shows Schmidt number, and k1andk2 describe the strength of homogeneous and heterogeneous reaction parameters, which are defined here as
(18)β1=I0μ0KTcπρlμ02,Pr=ρcpνk,α1a=λ1au0x−1,α1c=λ1cu0xSc=νfDA,k1=kca02S,k2=ksDAνfS.DAandDB are equivalent, i.e., δ=1; then, we can write as
(19)gη+hη=1.Utilizing Eqs. (15)-(16), we get the equation as
(20)1Scg″+fg′−k1g1−g2=0.The corresponding boundary equation of the concentration field yields the following form:(21)g′η=k2gη,η⟶0,gη=1,η⟶∞.The conversions of the wall drag and convective heat transfer into a nondimensional form are scaled as(22)Cf=−2τwρnfu02,Nu=xknfkfTc−Tw∂T∂yy=0.The wall shear stress is termed as(23)τw=μnf∂u∂yy=0,qw=−knf∂T∂yy=0.The dimensionless expressions of Eq. (22) are attained:
(24)12Re12xCf=11−ϕ−2.5f′′0,Re−12xNux=−knfkfθ1′0,where,
(25)Rex=ρnfu0xμnf.


## 3. Numerical Methodology

The saturation of two ferrite particles into Maxwell fluid under the action of a double magnetic dipole is examined with homogeneous and heterogeneous reactions. The transformed nonlinear systems of ODEs are solved numerically using the shooting method along with the RK-45 algorithm. For the implementation of the shooting method, one should convert the boundary value problem into an initial value problem. The reduce initial value problem is further converted into a system of first order differential equations and then solved by choosing the missing conditions as an initial guess. Therefore, the suitable transformations are used to obtain initial value problem. The new set of variables for first order system of equations is defined as(26)y1,y2,y3,y4,y5,y6,y7=f,f′,f″,θ,θ′,g,g′.This substitution yields(27)y1′=y2,y1′=y2,(28)y2′=y3,(29)y3′=1/11−ϕ2.51−ϕ+ϕρs/ρf−α1ay124∗−y1y32+α1a4y22y3η+2y1y2y3+βy41−ϕ+ϕρs/ρf,(30)y4′=y5(31)y5′=1/knf/kf1−ϕ+ϕ∗ρCps/ρCpf−α1cPrf24∗Pr∗−y1y52+y2y4−α1c42y1y3+2y3y2y5,(32)y6′=y7(33)y6′=‐Scy1y7−k1y61−y272.which are subject to the following conditions:
(34)y10=0,y20=0,y30=ω1unknowninitialcondition,y40=1,y50=ω2unknowninitialcondition,y70=k2y60,y80=ω3unknowninitialcondition.To solve the above system of Eqs. (27)-(33), the values of ω1, ω2, and ω3 are unknown, so by taking a suitable initial guess, the convergent numerical solution is obtained. It is important to note that if the boundary residuals are fewer than the tolerance error 10-6, the calculated solution converges. If the computed results do not satisfy this requirement, the initial estimates are changed by using Newton’s technique, and the procedure is repeated until the solution fulfills the specified convergence threshold. The thermophysical properties of the base fluid and ferrite particles for numerical procedure are defined in Table 1.Table 1
Thermophysical properties of blood, nickel zinc ferrite, and magnetite ferrite.
ρkg/m3CpJ/kgKkW/mKBlood1060.037700.52Nickel zinc ferrite48007106.3Magnetite ferrite51806709.7


## 4. Results and Discussion

The mathematical model of chemically reactive viscoelastics fluid over the stretching surface is solved numerically. The numerical solutions are checked by applying the grid independence test. The analysis is performed at different tolerance levels for Nusselt number and skin friction coefficient. The numerical values of Nusselt number and skin friction are presented in Table2.Table 2
Grid independence for Nusselt number and skin friction forβ=0.1,λ=0.01,γ=1.0,Pr=6.0..
ζNusselt numberSkin friction10‐23.94752.552310‐43.94712.551410‐63.94712.5514Table3 represents the performance of quantity of engineering interest (Nusselt number) against ferrohydrodynamic interaction (β) and Prandtl number (Pr) in both cases nickel zinc ferrite NiZnFe2O4 and magnetite ferrite Fe2O4. It is depicted that the increment occurs in the Nusselt number due to the increasing values ofPr, with better results of magnetite ferriteFe2O4as compared to nickel zinc ferriteNiZnFe2O4. Table 4 displays the response of mass flux coefficient against the strength of homogeneous and heterogeneous parameters and Schmidt number. It is perceived that the mass flux coefficient shows a positive trend towards the strength of homogeneous (k1) and heterogeneous parameters (k2) while an opposite behavior seen in the case of the Schmidt number (Sc). The higher magnitude was observed in the case of nickel zinc ferrite.Table 3
Obtained values of‐θ′0 by altering values of β and Pr.
βPr‐θ′0NiZnFe2O40.03.05541.053.01802.03.007443.00001.053.018063.1079Fe3O40.03.09041.053.04892.03.009343.00601.053.048963.1479Table 4
Obtained values ofg′0 by altering values of k1,k2 and Sc.
k1k2Scg′0NiZnFe2O40.11.37260.30.11.39310.51.41450.11.37260.10.21.11210.30.89770.11.37260.30.12.08050.52.50900.11.37260.10.21.11210.30.8977Fe3O40.11.28580.30.11.30510.51.32530.11.28580.10.21.05100.30.85550.11.28580.30.11.94810.52.34860.11.28580.10.21.05100.30.8555The analyzation of ferrohydrodynamic interaction (β1) on velocity field is seen in Figure 2. It is depicted that the velocity decreases due to the increasing values of ferrohydrodynamic interaction parameter (β1). This opposite behavior is due to the action of Lorentz forces which resist the flow and provide more resistance to transportation phenomena. From the graph, it is cleared that the magnetite ferrite offers more resistance due to which the velocity profile decreases in the case of magnetite ferrite Fe3O4. Figure 3 results an inclination in temperature with altering values of β1. It is detected that the temperature profile boosts up as soon as the ferrohydrodynamic interaction parameter (β1) enhances. Reasoned is that Lorentz forces which produces under the action of magnetic field have the potential to produce resistance which in term produce heating among the layers of fluids hence thermal boundary layer thickness rises. It is also perceived that the magnetite ferrite Fe3O4 offers more resistance due to which the rapid temperature enhancement occurred in the case of the magnetite ferrite Fe3O4. The sketch of velocity with augmenting values of Deborah number (α1a) is displayed in Figure 4. It is detected that the increasing Deborah number results in a decreasing velocity profile. As (α1a) precisely relate to relaxation time of fluid, immediately (α1a) increasing fluid relaxation time shot up and provide extra blocking to the fluid motion which results in thinning momentum boundary layer thickness. The maximum velocity is depicted in the case of nickel zinc ferrite NiZnFe2O4. The graph of Deborah number (α1a) against the temperature field is seen in Figure 5. It is perceived that the augmented values of the Deborah number (α1a) amplifies the temperature field. This phenomena can be depicted as the enlarging Deborah number (α1a) tends to have a larger relaxation time of fluid which directly means to offer resistance and generates heat which enhances the temperature profile. The effects of the thermal relaxation parameter (α1c) upon the velocity profile is observed in Figure 6. It is clear from the figure that the increasing thermal relaxation parameter (α1c) is responsible for the increasing velocity profile. This fact is explained as by up surging thermal relaxation parameter (α1c) momentum of fluid particles also enhances which in term increase kinetic energy and velocity profile. Figure 7 enlightened the conduct of thermal relaxation time (α1c) on the temperature field. By expanding values of (α1c) temperature profile show reverse behavior, because enlarging thermal relaxation time implies materials particles demand extra pace to exchanging energy to their connecting particles, thus cause reduction in temperature profile. The impression of Prandtl number on velocity is realized in Figure 8. The increasing behavior of velocity is observed by the altering values of the Prandtl number, which is due to straight forward relation between the Prandtl number and momentum diffusivity. As soon as the Prandtl increases, the momentum diffusivity also increases; thus, in the result, fluid motion and momentum boundary layer thickness also boost up. Association between Prandtl number and temperature field clearly delineated in Figure 9, the contrary response of temperature towards Prandtl number, actually Prandtl number exhibit inverse relation to thermal diffusivity of fluid, So by enhancing Prandtl number directly means to lessen diffused heat betwixt fluid layers hence reduces temperature profile. The response of Schmidt number (Sc) towards the concentration field is characterized in Figure 10. The increase in Schmidt number (Sc) consequently increases the concentration field because the concentration field relates with a viscous diffusion rate to molecular diffusion rate, so when the Schmidt number (Sc) magnifies the viscous diffusion rate, it also magnifies which in turn enhances the concentration field. Figure 11 shows the connection between the strength of homogeneous reaction parameter k1 and concentration field. It is observed by enhancing the strength of homogeneous parameterk1concentration field that is going to shrink. This conduct is explained as reactants are consumed during the process of homogeneous reaction which depresses concentration profile. Also, the minimum concentration was observed in case of magnetite ferrite Fe3O4. The graph of concentration field against the altering values of strength of heterogeneous reaction parameter k2 is displayed in Figure 12. An increase in k2 causes a reduction in concentration field, while higher concentration is found in the case of nickel zinc ferrite NiZnFe2O4. Figure 13 shows the impact of Pr number upon a heat transfer rate. As we increase the Prandtl number, the Nusselt number also increases because by increasing Prandtl number rate of momentum, the diffusivity also increases which gives rise to the kinetic energy and also the heat transfer through convection. The impact of strength of homogeneous parameter k1 and Schmidt number (Sc) in the mass flux coefficient g′0 is carried out in Figure 14. It is seen that (Sc) shows a negative response, while k1 displays an augmented amplitude against (g′0). Figure 15 depicts the trend of (g′0) that relates with the strength of the heterogeneous parameter and Schmidt number (Sc). As it is remarked because of the amplification in (k2), the mass flux coefficient (g′0) magnifies due to the accumulation of particles by the generation of the heterogeneous mixture. Also in the case of (Sc), we observed that as soon as (Sc) increases, the Sherwood number (g′0) reduces by reducing the mass diffusivity.Figure 2
Velocity profile with fluctuation inβ1.Figure 3
Temperature profile with fluctuation inβ1.Figure 4
Velocity profile with fluctuation inα1a.Figure 5
Temperature profile with fluctuation inα1a.Figure 6
Velocity profile with fluctuation inα1c.Figure 7
Temperature profile with fluctuation inα1c.Figure 8
Velocity profile with fluctuation inPr.Figure 9
Temperature profile with fluctuation inPr.Figure 10
Concentration profile with fluctuation inSc.Figure 11
Concentration profile with fluctuation ink1.Figure 12
Concentration profile with fluctuation ink2.Figure 13
Nusselt number with mutation ofβ and Pr..Figure 14
Mass flux coefficient with mutation ofk1andSc..Figure 15
Mass flux coefficient with mutation ofk2andSc..


## 5. Concluding Remarks

In this article, we have examined the response of double magnetic dipole in a chemically reactive viscoelastic fluid over a flat sheet by considering the two nanomagnetic ferrite particles. For heat transfer rate, we used the Cattaneo-Christov heat flux model which is the generalization of Fourier law by including thermal relaxation term. The thermal relaxation time converts the energy transport in form of thermal waves with finite speed. The modeled PDEs are converted into the system of ODEs and then solved numerically by shooting method. The results are depicted graphically with the impact of important parameters. The major outcomes of this study are as follows:(i)
Due to the force generated by magnetic dipoles, the velocity of the fluid reduces, and the temperature increases by increasing the ferrohydrodynamic interaction parameterβ1(ii)
For the large value of Deborah number (α1a), the velocity profile decreases(iii)
The thermal relaxation parameter (α1c) shows a positive response towards velocity while the negative response towards temperature because particles take a surplus time to shifting energy(iv)
The strength of homogeneous reaction parameterk1 undermines concentration(v)
The concentration field depresses due to the presence of strength of heterogeneous reaction parameterk2(vi)
It has been observed that heat transfer rate increases in the presence of magnetite ferriteNiZnFe2O4 as compared to nickel zinc ferrite NiZnFe2O4 by increasing Pr(vii)
From physical point of view, Pr is the proportion of momentum diffusivity to thermal diffusivity. The contrary response of temperature towards the Prandtl number is observed; actually, the Prandtl number exhibits an inverse relation to thermal diffusivity of the fluid, so by enhancing Prandtl number directly means to lessen the diffused heat between fluid layers and hence reduce the temperature profile



---

*Source: 2558419-2022-04-01.xml*