Novel Analysis of Two Kinds Hybrid Models in Ferro Martial Inserting Variable Lorentz Force Past a Heated Disk:An Implementation of Finite Element Method

Enran Hou,Umar Nazir,Samaira Naz,Muhammad Sohail,4,?,Muhammad Nadeem,Jung Rye Lee,Choonkil Parkand Ahmed M.Galal

1College of Mathematics,Huaibei Normal University,Huaibei,235000,China

2Department of Applied Mathematics and Statistics,Institute of Space Technology,P.O.Box 2750,Islamabad,44000,Pakistan

3Department of Mathematics,Government College University Faisalabad,Faisalabad,38000,Pakistan

4Department of Mathematics,Khwaja Fareed University of Engineering&Information Technology,Rahim Yar Khan,64200,Pakistan

5School of Mathematics and Statistics,Qujing Normal Uniersity,Qujing,655011,China

6Department of Data Science,Daejin University,Kyunggi,11159,Korea

7Research Institute for Natural Sciences,Hanyang University,Seoul,04763,Korea

8Mechanical Engineering Department,College of Engineering,Prince Sattam Bin Abdulaziz University,Wadi Addawaser,11991,Saudi Arabia

9Production Engineering and Mechanical Design Department,Faculty of Engineering,Mansoura University,Mansoura,P.O 35516,Egypt

ABSTRACT In this article,the rheology of Ferro-fluid over an axisymmetric heated disc with a variable magnetic field by considering the dispersion of hybrid nanoparticles is considered.The flow is assumed to be produced by the stretching of a rotating heated disc.The contribution of variable thermophysical properties is taken to explore the momentum,mass and thermal transportation.The concept of boundary layer mechanism is engaged to reduce the complex problem into a simpler one in the form of coupled partial differential equations system.The complex coupled PDEs are converted into highly nonlinear coupled ordinary differential equations system(ODEs)and the resulting nonlinear flow problem is handled numerically.The solution is obtained via finite element procedure(FEP) and convergence is established by conducting the grid-independent survey.The solution of converted dimensionless problem containing fluid velocity,temperature and concentration field is plotted against numerous involved emerging parameters and their impact is noted.From the obtained solution,it is monitored that higher values of magnetic parameter retard the fluid flow and escalating values of Eckert number results in to enhance temperature profile.Ferro-fluid flow and heat energy for the case of the Yamada Ota hybrid model are higher than for the case of the Hamilton Crosser hybrid model.Developing a model is applicable to the printing process,electronic devices,temperature measurements,engineering process and food-making process.The amount of mass species is reduced vs.incline impacts of chemical reaction and Schmidt parameter.

KEYWORDS Ferro-fluid; non-constant magnetic field; heated disc; thermal properties of nanoparticles; hybrid correlations among nanoparticles

Nomenclature:

Symbols/notation Used for

w,uVelocity components

ufNon-uniform axisymmetric velocity

σElectrical conductivity

νKinematic viscosity

kThermal conductivity

QHeat absorption/generation

T∞Ambient temperature

CConcentration

DMass diffusion

φ1,φ2Volume fractions?1,?2Variable thermal conductivity and concentration numbers

∞ Infinity

PrPrandtl number

HsHeat source parameter

θDimensionless temperature

ReReynolds number

NUNusselt number

FEM Finite element method

TiO2Titanium oxide

AStretching parameter

c,aConstants

x,rSpace coordinates(m)

B0Strength of magnetic field

ρFluid density

CpHeat capacity

dConsistent constant

TFluid temperature(K)

K0Chemical reaction number

C∞Ambient concentration

Tw,CwWall temperature and concentration(K)

LLength(m)

ηIndependent variable

MMagnetic field

EcEckert number

ScSchmidt number

?Dimensionless concentration

CfSkin friction coefficient

SHSherwood number

C2H6O2Ethylene glycol

SiO2Silicon dioxide

PDEs Partials differential equations

UwWall velocity(ms?1)

1 Introduction

Nanoparticles dispersed in the base fluid can improve the thermos-physical properties of nanofluids.Nanofluids outperform other working fluids in thermal conductivity and convective heat transfer.As a result,there has been an increase in the study of nanofluids in thermal systems in recent years.The shape of nanoparticle influences nanofluid-based thermal systems and hydraulic performance.The influence of nanoparticle shape on the performance of single-particle and hybrid nanofluid flows in various thermal environments has received a lot of attention.The thermophysical properties of nanofluids containing various-shaped nanoparticles have received much attention.The introduction summarizes various studies on the heat transfer characteristics of single-particle nanofluids with various particle shapes.Thermal energy conservation is gaining popularity to reduce heat loss while increasing heat transfer.Researchers are increasingly turning to nanofluids because conventional working fluids have limited improvement potential in heat transfer systems.The heat transfer of nanofluids in thermal systems such as tubes has received much attention.Algehyne et al.[1]studied the thermal transport of a Maxwell hybrid nanofluid on an infinitely stretchable vertical porous sheet.Hybrid nanofluids can be used to track thermal transport.A non-Fourier heat flux model can be used to exploit thermal failure.External heat sources must also be taken into account.The mathematical formulation of the model is based on a coupled nonlinear partial differential equations system.The partial differential equations model was converted to an ordinary differential equations model using transformations.Finite elements are then used to solve the nonlinear ordinary differential equations system.A grid-independent survey and a comparative analysis confirm the approach and solution’s dependability.The physical properties of the solution are investigated and recorded concerning several relevant variables.The higher the number in Forchheimer and Darcy’s porous numbers,the slower the motion.Rather than nanofluids,hybrid nanostructures accelerate thermal growth.Hybrid nanoparticles have a much thicker momentum layer than nanoparticles.A pseudoplastic model with tri-hybrid nanoparticles was used by Hou et al.[2]to model the thermal energy and solute transport problems to improve heat transfer.Thermal and solute transport was governed by the well-known Fourier and Fick’s laws.The coupled partial differential equations system was simplified using the boundary layer theory.The modeled partial differential equations were transformed into ordinary differential equations using a similarity transformation and then numerically solved using a Finite Element Procedure.The result was compared to many other factors.There is also a survey without grids.The fluid was discovered to be hotter than the tri-hybrid nanoparticles.The temperature and concentration of the fluid increased as the Dufour number increased.The fluid temperature increased as the heat generation parameter and Eckert numbers increased.The buoyancy force and porosity parameters were used to polarize the fluid velocity.Wang et al.[3]used a porous model of a stretched surface to study heat and mass transfer.The relationship between Darcy Forchheimer and the Forchheimer family is mentioned.Cattaneo and Christov’s double diffusion theories are used to derive the thermal and mass transport expressions.The thermal conductivity of tri-hybrid nanoparticles embedded in a pseudo-plastic material is also investigated.The phenomenon of boundary layers is the basis for this complicated model.Correlation can be used to convert tri-hybrid nanoparticle model partial differential equations to ordinary differential equations.Ordinary differential equations that have been transformed cannot be solved precisely.The numerical problem is solved using the finite element method.The researchers compare their findings to those previously published to ensure their approach is valid.To summarize,this research shows how finite element methods can solve complex engineering problems.Nanoparticles are also required to improve thermal performance and control mass diffusion.Hou et al.[4] investigated heat and mass transport in a pseudo-plastic model with Soret and Dufour effects over a porous surface with entropy analysis.The thermal energy and solute transport problems were modeled using heat generation/absorption and chemical reactions,including tri-hybrid nanoparticles to improve heat transfer.Thermal and solute transport was governed by the well-known Fourier and Fick’s laws.Using Cartesian coordinates and boundary layer theory,a set of coupled partial differential equations was simplified in the physical model.The modeled PDEs were transformed into corresponding ODEs using a similarity transformation and then numerically solved using the Finite Element Procedure.The result was compared to several other factors.There is also a survey without grids.The fluid was discovered to be hotter than the trihybrid nanoparticles.The temperature and concentration of the fluid increased as the Dufour number increased.The fluid temperature increased as the heat generation parameter and Eckert numbers increased.The buoyancy force and porosity parameters were used to polarize the fluid velocity.Hou et al.[5] studied a porous stretched surface with a hydromagnetic Carreau–Yasuda liquid,the effects of Soret and Dufour.Heat generation/absorption,viscous dissipation,and joule heating have significant effects.A chemical reaction can be seen during mass transportation.The conductivity of the Carreau–Yasuda liquid was investigated.PDEs were the governing laws in their early stages.The ODEs were then solved using the FEM (Finite Element Method).Grid-independent analysis was used to assess the methodology chosen.When parameter values were increased,significant physical effects were discovered.The rates of heat and mass transfer were calculated and discussed in detail.Nazir et al.[6] investigated the thermal efficiency of MoS2SiO2Williamson fluid with nano-sized hybrid particles.The Williamson fluid is thought to be thermally relaxed.These properties reduce the fluid’s temperature field.The fluid-particle resistive force controls the boundary layer.The resistive forces of Darcy and Forchheimer can also be compared.A hybrid nanofluid may be used to improve heat transfer.Hybrid nanofluid produces less heat than mono nanofluid.Relaxation impairs one’s ability to focus.Just as heat-generating chemical reactions affect temperature,they also affect mass mobility.It transports mass and heat better than a mono nanofluid.Using hybrid nanofluid instead of nanofluid as a heat-conducting material can significantly improve cooling system efficiency.Nazir et al.[7]have considered the Carreau-Yasuda model.Carreau and Yasuda’s shear rate-dependent model is used numerically to predict fluid flow dynamics.The use of hybrid nanostructures improves two-dimensional heat transfer.The pure Carreau-Yasuda fluid and the fluid with nanostructures are compared.The numerical simulations are carried out by the FEM software,which saves the results as graphs and numerical data.In both types of fluids,the wall shear stress is proportional to the magnetic field intensity.As the temperature of fluid rises,the heat transfer rate decreases due to Ohmic dissipation.The heat transfer rate of a hybrid nanofluid decreases faster than that of a mono nanofluid.The concentration profile clearly shows that generative reactions behave differently from destructive reactions.The Weissenberg number increases the fluid motion proclivity.In a hybrid nanofluid,the diffusion of wall movement has a more significant effect than the diffusion of solid boundary movement in a mono nanofluid.Rayleigh convection slows atom and ion movement in mono nano and hybrid nanofluids.Shear stress increases as the magnetic field strength increases.MHD fluids have higher shear stress at the wall than non-MHD fluids.Changes in magnetic field intensity have an impact on temperature profiles.The magnetic force causes electrons to slow down.As a result,boosting the magnetic field’s strength reduces convective heat energy transport.The melting parabolic stretched surface to focus on Powell–Eyring fluid flow was discussed by Chu et al.[8].Boundary layer theory contains flow rheology expressions.In this study,similarity transformations were used to convert PDEs to ODEs.These transformed ODEs were solved using the finite element method,including the symbolic computational.The grid-independent analysis demonstrates the utility and effectiveness of FEM.The surface drag and heat transfer coefficients show that the FEM is accurate and dependable.The nonlinear problem method is highly effective and can be applied to various nonlinear problems.The magnetic parameter is not inversely proportional to the thermal profile,but the Prandtl number behaves inversely.Waqas et al.[9] bio-convective flow of modified viscoelastic second grade and handled the complex equations numerically.They found that higher values of the Peckelt number reduce the density profile of microorganism.Nazir et al.[10]investigated the thermal performance of tri-hybrid nanomaterial in Sisko martial involving variable thermal properties.They implemented a numerical approach to obtain numerical results.Algehyne et al.[11]investigated the effects of thermal transport on a pseudo-plastic material by heating and stretching a porous sheet.To stimulate energy conservation,heat generation and absorption are used.Prandtl’s boundary-layer approach was used to model the phenomenon in Cartesian coordinates to make the problem easier to understand.The rheology is modeled using nonlinear coupled PDEs.Similarity transformation is used to convert these derived PDEs to ODEs.The converted ODEs with their new parameters were numerically approximated using the finite element method.In this study’s limiting case,it was discovered that the obtained solution was consistent with previously published findings.The finite element method is an excellent method for solving nonlinear problems in mathematical modeling,among other things.Hafeez et al.[12] created it to use a system of partial differential equations to simulate the thermal and solute properties of the Casson liquid in a Darcy porous medium.Convective boundary conditions determine how models behave in these situations.A 108-tolerance finite element method is used to solve these problems numerically.The subject of discussion is a mathematical algorithm(also known as a finite element approach).There is apparent convergence of 300 elements.Computer simulations study flow dynamics,heat transfer,and mass transport.The effects of temperature and concentration gradients on mass and heat energy transport are studied using simulations.It is fascinating to watch the numbers for Dufour and Soret change over time.Abdelmalek et al.[13]performed both the numerical and parametric analysis,as evidenced by their simulations.Suspending hybrid nanostructures significantly improved thermal conductivity and mass transport.When thermal radiation is increased,the temperature drops significantly.The Weissenberg parameter affects the velocity of a Carreau liquid in a magnetic field.Specie diffusion from the wedge into the fluid increases when the chemical reaction is generative;the opposite trend occurs when the chemical reaction is destructive.Imran et al.’s[14]research focused on the unsteady two-dimensional flow of a second-grade hybrid nanofluid over a stretchable moving horizontal porous plate.The entropy of the system is taken into account.The magnetic field and Joule heating are also considered.This experiment uses silicon carbide and titanium oxide nanoparticles dispersed in oil.Small particles(such as spheres,bricks,tetrahedrons,and platelets) are also studied for their shape factors.PDEs are represented mathematically by accounting for Joule heating and viscous dissipation effects.The current PDEs are converted into ODEs using MATLAB’s bvp4c (shooting) scheme.The effects of various physical constraints on momentum and temperature fields are graphed and numerically investigated.Manzoor et al.’s [15] work on improving heat transfer in industrial and engineering applications has piqued the interest of researchers.For example,in heat exchangers and electronic devices,the rate of thermal transport is critical.Oil,water,and ethylene glycol are examples of low thermal conductivity base fluids.Hybrid convective heat transfer has a wide range of manufacturing applications.Researchers investigated heat transfer in hybrid nanoparticles with magnetic dipoles using blood as a base fluid.Using similarity transformations,the flow problem of PDEs can be transformed into the structure of nonlinear ODEs.Several recent contribution covering the modeling of physical aspects are covered in [16–21] and references therein.Mahbood et al.[22] developed a system for performing his calculations.The MHD couple stress fluid was studied in greater depth using the fourth–fifth-order Runge–Kutta–Fehlberg equations,including nonlinear thermal radiation,velocity slip,and entropy generation.Mahbood et al.[23] investigated the heat and mass transfer properties of a Brinkmann-type nanofluid flow composed of sodium alginate and Fe3O4in a vertical rotating frame.The Runge-Kutta Fehlberg method(RKF-45)solves dimensionless nonlinear partial differential equations.Ayub et al.[24]investigated the second-order slip phenomenon and the flow of cross nanofluids that comes to a halt at a certain point using a spectral relaxation method.Brownian motion,joules,viscous heating,thermophoresis,and mixed convection are all considered.The way MHD transferred heat was studied by Shah et al.[25].Thermal conductivity allows a nanofluid to pass through using thermophoresis,Brownian motion,and a binary chemical reaction in cylindrical coordinates.The velocity and activation energy of nanofluid are determined using perpendicular magnetic dipole field measurements.Wang et al.[26] estimated numerical consequences of traveling waves in terms of chemical kinetics using equations of time-fractional fishers.Some recently published contributions are mentioned in[27–32].

Available studies ensure that no report is presented on heat and mass transport for axisymmetric flow produced due to the stretching of a heated rotating disk with Yamada Ota hybrid and Hamilton Crosser models containing the variable properties via finite element scheme.This work fills this gap.Section 1 contains the comprehensive literature review,modeling is covered in Section 2 and the scheme is explained in Section 3 and implemented,the obtained solution is discussed in detail in Section 4 and important results are highlighted in Section 5.

2 Formulation of the Developed Model

The model is prepared using the following assumptions:

■A rheology of Ferro-fluid is considered;

■A axisymmetric heated disc is taken out;

■Variable magnetic field and variable wall temperature is addressed;

■Yamada Ota hybrid and Hamilton Crosser models are used;

■Viscous dissipation and Joule heating influences are observed;

■Motion of fluid particles is induced using the movement of wall;

■Thermal properties regarding nanoparticles are inserted;

■An efficient approach based on the finite element method is used;

■Chemical reaction is occurred into mass species;

■Non constant thermal conductivity and non-constant mass diffusion are addressed;

■Fig.1 represents the flow behavior of the current model.

Figure 1:Geometry of model

A variable magnetic field [0,0,B0dr?1/3] is imposed at surface of wall.It is noticed that fluidic motion is generated by stretching surface in the radial direction along with velocityUw=ar1/3whereas ambient fluid is generated with velocityuf=cr1/3.The formulated PDEs which present the continuity equation,momentum,heat and mass transport are derived via boundary layer approximations and these obtained PDEs[29,31,32]are

Implication of desired boundary conditions[29]is considered as

Hybrid correlations among nanoparticles[30]are delivered as in Table 1.

Table 1: Physical properties density,heat capacity and thermal conductivity[2,3]

Similarity variables[29]are

Thermal conductivity and mass diffusion are functions of temperature[31]which are addressed as

Using similarity variables in system of PDEs and obtained ODEs[29]are

Dimensionless boundary conditions[29]are

Dimensionless parameters are derived as

Drag force coefficient for Ferro-fluid[29,32]is

Temperature gradient and rate of mass diffusion under hybrid nanofluid and variable properties are modeled[29,32]as

3 Numerical Procedure

Model problem along with boundary conditions is simulated by an efficient technique called finte element method.Numberious computational fluid dynamics problems are solved by FEM.It is noticed that code regarding finite element method is designed on MAPLE 18.Moreover,Maple is known as numerical and symobolic computing enviroment.It deals with numberious areas,such as numerical analysis,data processing,symobilc mathematics and visulizations.Basicscally,six steps of FEM are discussed here and these steps are skecched by Fig.2.

Figure 2:Illustration of FEM’s steps

Step I:The subdomains are obtained by dividing the domain into elements which are connected with neighboring elements at nodes.The discretization of the elements gives approximation solution over each element rather than over whole domain and approximation solution is considered as a linear polynomial.

Step II:The shape functions are orthogonal nodal basis which possess the property Kronecker delta.The approximation solution is obtained from the product of shape functions.The linear shape functions are used here.

Step III:The weak formulation of the transformed problem is covered in this step.

Step IV:The elements stiffness matrix is generated due to using of Galerkin approach in weak form whereas concept of assembly process is helped to generate elements stiffness matrices.

Step V:Finally,a global stiffness matrix is obtained via assembly process.Picard linearization method is used to linearize the system of non-linear equations.

Step VI:Finally,algebraic equations are achieved using the Picard linearization method.Code is designed in MAPLE 18.The table is tabulated to check the convergence of the developed analaysis.It is observed that outcomes regarding concentration,velocity and temperature profiles at the mid of each element are simulated in the Table 2.Moreover,the number of elements is increased against velocity,temperature and concentration profiles and obtained outcomes are noticed.It is found that the results of velocity,temperature and concentration profiles are not changed after 300 elements.Therefore,the results of the developed model are converged at 300 elements and all numerical and graphical simulations are obtained for 300 elements.

Table 2:Numerical simulations of concentration,heat energy and velocity at mid of each 300 elements when M=0.001,Pr=206,A=0.3,?1=0.02,?2=0.08,Ec=1.4,Hs=?1.4,Sc=0.7,Kc=0.2

4 Results and Discussion

Mass diffusion and heat energy in Ferro-fluid including variable properties are analyzed in the heated disk.Yamada Ota hybrid and Hamilton Crosser models are inserted.Heat source and Joule heating are observed along with chemical reaction.The feature based on the stagnation point is studied.A complex model is developed while the complex model is simulated with the help of a numerical approach called the finite element method.The ranges of dimensionless parameters 0.001≤M≤0.4,≤,0≤A≤0.3,206≤Pr≤208,?1.8≤Hs≤0.3,0.0≤Sc≤0.8,0.0≤Kc≤0.8,0.0≤Ec≤2.0,0.0≤?1≤1.5,0.0≤?1≤1.7 are used in simulations.Variation in various parameters is accumulated on flow,mass diffusion and heat energy which are explained below.

4.1 Consequences Regarding Flow and Thermal Energy

In this subsection,flow analysis and heat energy analysis are studiedvs.various parameters.It is mentioned that two kinds of models regarding the hybrid approach are studied.Moreover,solid curves are plotted for representation of the Yamada Ota hybrid approach while Hamilton Crosser is captured to address the consequence of the Hamilton Crosser model.Fig.3 is addressed to predict the behavior of the magnetic number on velocity curves.It is investigated that flow becomes slow downvs.higher impact of the magnetic parameter.It is because Lorentz force is applied along the vertical direction of the heated disk.Therefore,flow is slowed down due to negative Lorentz force.Physically,the magnetic parameter is formulated using the concept of Lorentz force in momentum equations.It is noticed that the magnetic parameter is implemented along perpendicular direction on the wall.Therefore,the fluid becomes significantly viscous and the thickness associated with momentum layers decreases.The role of the Eckert number is implemented on the temperature field presented in Fig.4.Maximum amount of heat energy is gainedvs.higher impact of Eckert number.It is studied that appearance of Eckert number is produced because of viscous dissipation.More heat dissipates when Eckert number is increased.Therefore,direct proportional relation is estimated among heat energy and Eckert number.Moreover,maximum production regarding thermal energy is generated for the Yamada Ota hybrid model rather than heat energy for the case Hamilton Crosser hybrid model.Eckert number is a dimensionless number formulated due to in viscous dissipation in the energy equation.An increment in Eckert number results in viscous dissipation which is enhanced.Hence,heat energy is significantly dissipated when the Eckert number is increased.Fluid particles gain more heat energyvs.the impact of Eckert number.Fig.5 addresses an impact of heat source parameter on thermal curves.Thermal curves are plotted for positive numerical values as well as negative values of heat source parameter.It is investigated that negative numerical values for heat source parameter are predicted for heat absorption while positive numerical values for heat source parameter are observed for heat generation phenomenon.It is observed that the heat generation parameter is implemented at the wall of the surface.Hence,thermal energy can be adjusted against the variation in heat source parameter.Moreover,thickness related to thermal layers is enhanced against higher numerical values of heat source number.The observation of magnetic parameter on heat energy curves is analyzed in Fig.6.It is observed that heat energy is enhanced into fluidic particles when magnetic parameter is increased.The occurrence of magnetic parameter on heat energy is generated due to applying the impact of Joule heating.Physically,the Joule heating process is occurred by interaction among moving fluidic particles that make current.Thermal layers thickness is reduced when the magnetic parameter is increased.Moreover,maximum production regarding thermal energy is generated for Yamada Ota hybrid model rather than heat energy for the case Hamilton Crosser hybrid model.

Figure 3:An influence of M on velocity curves when Pr=208,A=0.0,?1=0.03,?2=0.08,Ec=1.4,Hs=?1.8,Sc=0.7,Kc=0.3

Figure 4:An influence of Ec on thermal curves when M=0.001,Pr=206,A=0.0,?1=0.04,?2=0.08,Hs=?1.4,Sc=0.7,Kc=0.4

Figure 5:An influence of Hs on thermal curves when M=0.03,Pr=206,A=0.0,?1=0.02,?2=0.08,Ec=1.4,Sc=0.7,Kc=0.2

Figure 6:An influence of M on thermal curves when Pr=206,A=0.0,?1=0.02,?2=0.08,Ec=1.4,Hs=?1.4,Sc=0.7,Kc=0.2

4.2 Consequences Regarding Mass Diffusion

Figs.7–9 are plotted to determine influences ofSc,Kcand?2on concentration field including models related to Hamilton Crosser and Yamada Ota hybrid models.Fig.7 is addressed the effect of Schmidt number on concentration field.It is estimated that Curves associated with concentration field are declinedvs.an impact of Schmidt number.This is because Schmidt number is ratio among momentum and mass diffusivities.Therefore,an inversely proportional relation is investigatedvs.an impact of Schmidt number.Schmidt number is known as dimensionless number.Mathematically,is formulated mass diffusion equation.Hence,ration amongνfandDfis called Schmidt number.Moreover,maximum production regarding concentration is generated for the Yamada Ota hybrid model rather than concentration for the case Hamilton Crosser hybrid model.Fig.8 is prepared to predict the amount regarding concentrationvs.chemical reaction number.Diffusion into fluid particles becomes slowed down when the chemical reaction parameter is increased.The chemical reaction number is known as coefficient of temperature based on the concentration field.It can be visualized by variation inKc.It is observed thatKc <0 is the case when solute particles are generated during the chemical reaction.Such kind of chemical reaction is known as a generative chemical reaction.Further,the chemical reaction is called destructive chemical reaction when solute is used in chemical reaction.The role of?2on mass diffusion is investigated by Fig.9.It is noticed that process of mass diffusion is enhanced when?2is increased.Mathematically,variable mass diffusion number has direct relation against temperature difference.Further,?2is the dimensionless parameter,which is formulated variable mass diffusion in mass diffusion equation.Mass diffusion for the case of?2=0 is higher than mass diffusion for the case of?2/=0.Thickness associated with mass diffusion is enhancedvs.higher numerical values of?2.

Figure 7: An influence of Sc on concentration curves when M=0.04,Pr=206,A=0.0,?1=0.02,?2=0.08,Ec=1.4,Hs=?1.4,Kc=0.8

Figure 8: An influence of Kc on concentration curves when M=0.001,Pr=206,A=0.0,?1=0.02,?2=0.08,Ec=1.9,Hs=?1.4,Sc=0.9

Figure 9: An influence of ?2 on concentration curves when M=0.3,Pr=209,A=0.0,?1=0.02,Ec=1.4,Hs=?1.8,Sc=0.7,Kc=0.2

4.3 Investigation of Nusselt Number,Flow Rate and Sherwood Number

Table 3 investigates numerical aspects of Nusselt number,flow rate and Sherwood number against the variation in magnetic number,heat source parameter,Dufour number and Schmidt number.Flow rate,mass diffusion rate and heat energy rate are declined when the heat source parameter and magnetic parameter are increased.But heat transfer rate and mass diffusion rate are boosted against argument numerical values of Dufour number.

Table 3:Numerical aspects of skin friction coefficient,Sherwood number and Nusselt number against M,Hs,Sc and Kc when Pr=206,A=0.3,?1=0.02,?2=0.08,Ec=1.4

Table 3 (continued)Variation in parameters?Re1 2 Cf?Re?1 2 NU?Re?1 2 SH Kc0.60.3174218242 1.3466823431.389752571 1.20.3389184670 1.3506041611.391472897

5 Conclusions

Dispersion of hybrid nanoparticles is carried out in this research with variable properties for the fluid flow through a heated disc.The arising problem is tackled numerically.The remarkable findings are:

■Ferro-fluid flow and heat energy for the case of the Yamada Ota hybrid model are higher than for the case of the Hamilton Crosser hybrid model;

■Convergence of the present problem is ensured by 300 elements;

■Developing a model is applicable to printing process,electronic devices,temperature measurements,engineering process and food making;

■Heat energy is boosted against variation in Eckert number and heat source number,but production regarding thermal energy is declinedvs.variation in Joule heating process;

■The amount of mass species is reducedvs.incline impacts of chemical reaction and Schmidt parameter;

■Mass diffusion rate and thermal energy rate are increased against variation in magnetic parameter,but mass diffusion rate and thermal energy rate are declinedvs.heat source number.

Availability of Data and Material:The data used to support this study are included in the Manuscript.

Funding Statement:The authors received no specific funding for this study.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.