Thermal transport exploration of ternary hybrid nanofluid flow in a non-Newtonian model with homogeneous-heterogeneous chemical reactions induced by vertical cylinder

Studying the combination of convection and chemical processes in blood flow can have significant applications like understanding physiological processes, drug delivery, biomedical devices, and cardiovascular diseases, and implications for various fields can lead to developing new treatments, devices, and models. This research paper investigates the combined effect of convection, heterogeneous-homogeneous chemical processes, and shear rate on the flow behavior of a ternary hybrid Carreau bio-nanofluid passing through a stenosed artery. The ternary hybrid Carreau bio-nanofluid consists of three different types of nanoparticles dispersed in a Carreau fluid model, miming the non-Newtonian behavior of blood. This assumed study generates a system of PDEs that are processed with similarity transformation and converted into ODEs. Furthermore, these ODEs are solved with bvp4c. The results show that the convection, heterogeneous-homogeneous chemical processes, and shear rate significantly impact the bio-nano fluid’s flow behavior and the stenosed artery’s heat transfer characteristics.


Introduction
A chemical mechanism including both homogeneous and heterogeneous chemical reactions leads to the development of atherosclerotic plaques in arteries.Homogeneous chemical processes occur within a single phase, such as the bloodstream oxidizing low-density lipoprotein (LDL).These reactions involve the adhesion, migration, and transformation of monocytes into macrophages, the uptake and oxidation of LDL, and the release of pro-inflammatory cytokines.So, it's essential to know how different types of chemical reactions work together to make suitable atherosclerosis treatments.The articles [1][2][3] are helpful to understand these concepts in this regard.Vaidya et al. 4 presented the theory of how blood flows through small arteries using a mixture called MHD Bingham fluid.It looks at how the fluid moves due to peristalsis, a natural squeezing motion.The article also examines how several chemical reactions occur in this process.Tanveer et al. 5 made a theoretical exploration of peristaltic activity within a magnetohydrodynamics (MHD) based blood flow system comprising non-Newtonian fluid behavior.Later on, some practical, theoretical investigation of peristaltic mechanism is the Sutterby fluid model concerning thermal transport involving chemical reaction was discussed by Imran et al. 6 Keeping in focus subjected to cubic autocatalysis, Shah et al. 7 elucidated the intricacies of temporal variations in cross-nanofluid behavior on a melting surface.Khan et al. 8 investigated activation energy and thermal performance in the presence of Nano lubricants.In this study, a chemical reaction is added to check the mass transport of fluid.Zafar et al. 9 revealed the impact of an irreversible chemical reaction with Prandtl nanofluid attached to the stretched surface.Some of the latest studies [10][11][12] are related to the chemical process in the stagnation point flow of a water-based graphene-oxide, the significance of thermophoresis in particle deposition, and cross-liquid flow past a moving wedge.
The property of infinite shear rate viscosity is a crucial and fundamental rheological attribute that is define the fluidic behavior.The high shear rate that approaches infinity signifies the apparent viscosity of a substance, as it operates at a pace significantly greater than the ordinary range of shear rates encountered within practical applications.Infinite shear rate viscosity is significant in understanding the flow properties of complex fluids such as polymers, colloids, and biological fluids, which can exhibit non-Newtonian behavior at high shear rates.Multiple studies in this regard were conducted over time which are listed in References 13-16.The infinite shear rate viscosity and heat transport properties of magnetized Carreau nanofluid are examined by Ayub et al. 17 The Carreau nanofluid model, characterized by infinite shear rate viscosity, has several noteworthy energy transportation features, including its inclination and magnetization properties.These features are particularly studied by Shah et al. 18 He also listed many advantages of the proposed model when applied to a wedge geometry.Later on, Khan and Sardar 19 talked about how a certain type of fluid flows over a wedge shape.The fluid has ''infinite shear rate viscosity,'' affecting how it moves over the wedge.He concluded that it was steady and happening in two dimensions.Wahab et al. 20 researched the infinite shear rate viscosity of Carreau fluid with dual aspects.This research aimed at the homogeneous or heterogeneous and inclined magnetic characteristics of an infinite shear rate viscosity model of Carreau fluid in conjunction with nanoscale heat transport.
Ternary hybrid nanofluids has demonstrated significant potential in enhancing fundamental fluids' thermal transport characteristics in numerous heat transfer contexts.In Ternary hybrid nanofluids, mixing three types of nanoparticles with liquid could improve blood flow in terms of thickness and force.Ternary hybrid nanofluids can potentially increase the base fluid's thermal conductivity and heat transfer rate, which can help regulate blood flow in areas of the body affected by ischemia.Adding nanoparticles made of carbon to three different fluids could help blood carry more oxygen.This may help people who have anemia feel better.However, the synthesis and characterization of Ternary hybrid nanofluids for biomedical applications are still in the early stages of research.The following research articles [21][22][23][24] tell give more knowledge and show great significance.Ahmed et al. 25 made exploration about the heat transfer enhancement of a novel mixed metal oxide ZnO-Al 2 O 3 -TiO 2 nanofluid in a square flow conduit.Thirumalaisamy et al. 26 conducted an analysis of fluid flow and heat transfer in an inclined rectangular porous cavity using a ternary aqueous Fe 3 O 4 -MWCNT-Cu nanofluid.Ramesh et al. 27 made study related to the effects of heat source/sink and porous medium in a stretchable convergent/divergent channel with a ternary nanofluid.Alharbi et al. 28 investigated the flow pattern and heat transportation of a Darcy ternary-hybrid nanofluid over an extending cylinder.In this study they used induction effects and numerical solution is fetched using computational methods.Mumtaz et al. 29 used unique approach to examine the radiative ternary nanofluid flow on curved geometry with second-order velocity slip constraints and cross-diffusion.Additionally, many there are latest studies [30][31][32] which are related to thermosolutal performance, entropy generation, effects of temperature and particle volume concentration and magnetic flux density and heat source/sink of a ternary hybrid nanofluid in a different geometries.

Motivation
Heart diseases are a significant reason why many people die or get sick all around the world.Even doctors and scientists still have difficulty finding ways to deal with these health problems.Atherosclerosis, the buildup of plaque in the arteries, is a common cardiovascular disease that can cause stenosis, the narrowing of the arteries, and reduce blood flow.This research paper aims to study how different things affect the way a special fluid flows through a narrow part of the body called an artery.We're interested in the fluid's movement, the chemical reactions that happen inside it, and how fast the fluid is moving.

Novelty
This research presents an exceptional aspect of blood flow because it looks at how flow rate and chemical processes influence the posited blood velocity in arteries.This is a unique or novel avenue for enhancing blood flow in stenotic arteries and, which could help make it easier for blood to flow through clogged-up arteries.Furthermore, analysis of the effects of various factors and especially the inclusion of the heterogeneous-homogeneous chemical processes within the blood flow behavior and heat transfer characteristics of the bio-nanofluid provide a comprehensive understanding of the complex fluid dynamics in stenosed arteries.

Mathematical formulation of the flow problem
It is assumed that blood is moving in a straight, flat path inside an artery having length L o 2 units long.Blood moves through the x-axis and r-axis.The blood is flowing with velocity u w = u 0 x L o and the artery is affected by a magnetic field at an angle v to the base axis.To examine the temperature of blood, various factors like heat absorption/release, and thermal radiation must be considered.Ternary nanofluid has been utilized to examine the effects of three different types of nanoparticles on blood flow within an artery.The considered artery has the maximum length l and the radius 33 , where its mathematical form represents the area Moreover, assuming that E and F act as autocatalysts, isothermal reactions will exhibit autocatalysis.
The reaction taking place on the catalytic surface follows a single-order isothermal process, expressed as.
It's important to remember that away from the field, the rate = k s G a G b 2 will be zero and affect the artery.k 1 and k s are numbers that tell us how fast reactions happen.These reactions can either happen in one place or many places at once.G a and G b are numbers that are used to show concentration.More simply E and F tell us how much of two different things are in a chemical reaction (Figure 1).
Therefore, equation of continuity, momentum and heat is given [34][35][36][37][38][39] as under: The appropriate boundary conditions are (see References 40-42): In above equations, the correlations of the thermophysical characteristics for ternary hybrid nanofluid are given by [43][44][45][46] : Electrical conductivity Moreover, the numerical data of the thermophysical characteristics of the three distinct nanoparticles such as TiO 2 , Au, and A galong with based fluid (blood) are given in Table 1.
According to Roseland approximation, the thermal radiative heat flux is mathematically expressed as: where s ÃÃ and k ÃÃ are named as Stefan Boltzmann constant and the mean absorption coefficient, respectively.Moreover, the term T 4 is further exercised by using the Taylor series at the point T ' and ignoring the higher order power terms which becomes T 4 ffi 4TT 3 ' À 3T 4 ' (Figure 2).

Solution procedure of the similarity equations
For conversion of PDEs into ODEs, the following similarities 41,48,49 are introduced as: By giving inputs of above equation into equations ( 4)-( 7), following results appeared in the form of ODEs are as follows: The associated transformed BCs are: In addition, the above equations comprised distinct influential parameters which are mathematically defined as: The important component of the engineering physical quantities are the skin friction coefficient and local Nusselt number which are defined as follows 50 : where the wall shear stress is defined as: The wall heat flux is given by Substituting the similarity variables and the equations ( 23) and (24) in equation ( 22), the reduce form of the skin friction and local Nusselt number can take place the form as: In which, A 1 , A 2 , A 3 , A 4 and A 5 are given by:

Numerical scheme
Several numerical techniques [51][52][53][54][55][56][57][58][59][60][61] exist in the literature which are useful for solving ODEs and fractional ODEs.The Runge-Kutta fifth order technique, commonly abbreviated as RK-5, is a computational approach/numerical approach.RK-5 (Runge-Kutta 5th order) is commonly employed to solve differential equations.The method in question constitutes a fourthorder Runge-Kutta enhancement at the higher-order level, thereby rendering it capable of greater precision than its counterpart for a given step size.This denotes an advancement in accuracy concerning numerical methods in mathematical computing.The RK-5 methodology is employed to evaluate the solution of a differential equation by performing an iterative process of intermediary value computations.The intermediate values are derived from the slope of the solution at multiple points within the domain and are subsequently integrated to yield the ensuing approximation of the solution.The RK-5 technique is extensively utilized to resolve differential equations, particularly in situations requiring heightened precision.
To use the RK-5 method in MATLAB for solving a differential equation, we need to do the following steps:

Validation of the numerical scheme
In this portion of the work, we have shown the authentication of the given scheme for the limiting cases.Therefore, the present study outcomes of the heat transfer rate with available prior work of Sarwar and Hussain 62 for the various values of g due to the absence of the other parameters such as Weissenberg number, thermal radiation parameter, heat generation parameter, and temperature ratio parameter is taken to be unity.The outcomes are shown in Table 2, where the existing work are excellent match with prior research work can give confidence that the scheme is accurate and can be applicable for the unavailable work.

Analysis of the results
The Runge-Kutta fifth order technique, commonly abbreviated as RK-5, is a computational approach/ numerical approach.RK-5 (Runge-Kutta 5th order) is widely employed to solve differential equations.The method in question constitutes a fourth-order Runge-Kutta enhancement at the higher-order level, thereby rendering it capable of greater precision compared to its counterpart for a given step size.This denotes an advancement in accuracy with respect to numerical methods in mathematical computing.The RK-5 methodology is employed to evaluate the solution of a differential equation by performing an iterative process of intermediary value computations.The intermediate values are derived from the slope of the solution at multiple points within the domain and are subsequently integrated to yield the ensuing approximation of the solution.The RK-5 technique is extensively utilized to resolve differential equations, particularly in situations requiring heightened precision.
The Weissenberg coefficient is elevated, the viscoelastic characteristics of blood can cause it to display non-Newtonian behavior.The flow conditions of the blood may bring about varying viscosity, thus resulting in intricate flow patterns comprising of turbulence and the forming of vortices.Intricate flow patterns can create heightened resistance to blood flow, ultimately leading to a decrease in the flow rate.This information is shown in Figure 3(a) and (b).Regarding the velocity profile, the elicitation of v and M is depicted in Figure 4(a) and (b).This figures show that the blood flow goes down when the value of both numbers/parameters gets higher.Strong magnets/inclined magnetic effect causes a force called Lorentz force.This force can slow blood flow and change how much blood moves through your body.This can also affect your blood pressure.The impact g in relation to blood velocity is depicted in Figure 5(a) and (b) while accounting for the presence of the assisting/opposing force parameter.The larger the numerical value entered into g, the greater the resulting velocity of the fluid.The occurrence of stenosis, leading to the constriction of an artery, results in alterations in the curvature of the vessel.This phenomenon can induce modifications in the blood flow and its associated patterns, consequently having noteworthy physiological implications.Figure 6(a) and (b) shows that how the Q parameter affect the blood flow temperature, with assisting/opposing force parameter.The temperature of blood flow within a constricted artery can be impacted by the heat generation parameter.When an artery becomes stenosed, the flow of blood can be hindered because the vessel has narrowed, resulting in elevated levels of shear stress and turbulence.The generation of heat within the blood can influence the temperature of the circulating blood.Higher heat generation in stenosed arteries increases blood flow temperature with physiological effects.Heat Step 1: Define the differential equation as a MATLAB function.
function dydt = odefun (t, y) % Define the differential equation dydt =.% expression for the derivative of y at time t Step 2: Set up the initial conditions for the problem.t 0 =.% initial time y 0 =.% initial value of y Step 3: Choose the step size.h =.% step size Step 4: Set up a loop to iterate over the time steps and compute the solution using the RK-5 method.t = t 0 ; % initialize time.y = y 0 ; % initialize solution.while t \ tf % tf is the final time.A statistical representation of the heat transportation rate during the analysis of blood flow with varying R d is presented in Figure 9. Figure 10 shows how quickly heat moves in the bloodstream, depending on the position of u w .which clearly that fluctuation of heat transmission during the analysis of blood flow as u w changes.Moreover, Table 3 shows physical characteristics with symbols and ranges given as intervals.Table 4 shows how a particular thing affected different profiles and explains why they reacted that way.

Final remarks
In this research paper, we have delved into the intricate mechanisms of heat transport and synergistic integration of convective chemical processes within a cutting-edge Ternary Hybrid Carreau Bio-Nanofluid,     focusing on its application in blood.This research attempt gave us valuable insights into the unique thermal behavior and chemical interactions occurring in this innovative fluid system.This study contributes to understanding heat transport in complex fluids and paves the way for future advancements in bio-nanofluid dynamics.The wise conclusion is listed as below; 1. Lorentz force reduces the flow rate.2. The curved path of the blood reduces rubbing against the vessel walls.3. High We values make blood more elastic, which increases resistance and decreases blood flow.4. As the strength of a homogeneous reaction parameter increases, the reaction rate increases.5.A stronger heterogeneous reaction can speed up chemical reactions, potentially impacting blood flow.6.The strength of a heterogeneous reaction parameter rises, and the rate of the chemical reaction increases, which can have both positive and negative effects on blood flow.

Limitation of study
This study may encounter several limitations that could affect the findings' scope, reliability, and generalizability.Some potential limitations are included; 1.The model may involve simplifying assumptions, such as assuming steady-state conditions or neglecting certain physical phenomena, which might limit the realism of the study.2. The numerical methods employed to the governing equations may have limitations, such as convergence issues, which can impact the accuracy and stability of the simulations.3. Assumptions about the Thermophysical properties of the bio-nanofluid may not accurately represent real-world conditions where these properties can vary with temperature and concentration.4. The study may make assumptions about the kinetics of chemical reactions, and the chosen reaction mechanisms may not fully capture the complexity of real chemical processes.5. Modeling homogeneous and heterogeneous chemical reactions simultaneously could be challenging, and the assumptions about the interplay between these reactions may introduce uncertainties.6.The lack of experimental data for validation could limit the verification of the model's accuracy and applicability to real-world situations.7. The study may assume specific boundary conditions that might not be entirely realistic for certain practical scenarios, leading to potential discrepancies between the model and actual experimental setups.

Future work
Here are some potential future work suggestions for the current research.
1. Exploration may be made to find the numerical solution of a multi-phase system.
2. The influence of a magnetic field can be judged with Cross fluid model.3. The optimal parameters for enhancing heat transfer and chemical reactions can be identified.4. Validation with the theoretical findings through experimental studies may be taken place in future studies.5. Investigating the impact of different nanoparticle shapes on the heat transfer and chemical reactions within the hybrid Carreau bionanofluid may also be considered.6. Artificial neural network technique may be utilized to solve the mathematical model.

Figure 1 .
Figure 1.Geometry of blood flow via arteries.

Figure 2 .
Figure 2. Overview of the study mechanism.
) and (b) and 8(a) and (b) illustrate the regulation of u(h) within blood circulation k 2 , L b using with assisting/opposing force parameter.Both parameters result in elevated concentrations.Chemical reactions may transpire within the blood flow, including the binding of oxygen to hemoglobin in erythrocytes or the enzymatic degradation of glucose to yield energy.The reactions transpire endogenously within the human organism and are subject to regulation by enzymatic and other biochemical mechanisms.Simply put, we may say that chemical changes can happen in our blood.Some examples are when oxygen attaches to red blood cells or when glucose gets broken down for energy.The body has natural reactions that are controlled by enzymes and other chemicals.

Figure 3 .
Figure 3. (a and b) Velocity profiles f (h) with several values of We.

Figure 4 .
Figure 4. (a and b) Velocity profiles f (h) with several values of vand M.

Figure 6 .
Figure 6.(a and b) Temperature profiles u(h) with several values of Q.

Figure 5 .
Figure 5. (a and b) Velocity profiles f (h) with several values of g.

Figure 7 .
Figure 7. (a and b) Results k 2 on u(h).

Figure 8 .
Figure 8. (a and b) Results of L b on u(h).

Figure 10 .
Figure 10.Rate of heat transport in blood flow analysis with variation of u w .

Table 2 .
Comparison of the heat transfer rate for the several values of g due to the absence of the several influential parameters when u w = 1.

Table 3 .
Physical parameter, symbol, and range in intervals.

Table 4 .
Tabular results for parameter impact on profiles with explanations for their behavior.

Table 5 .
Numerical assessment of attached physical parameters on Local Nusselt Number and Local Skin Friction of Ternary hybrid/ Ternary Hybrid Carreau nanofluid (blood).