Home Heat transfer characteristics of cobalt ferrite nanoparticles scattered in sodium alginate-based non-Newtonian nanofluid over a stretching/shrinking horizontal plane surface
Article Open Access

Heat transfer characteristics of cobalt ferrite nanoparticles scattered in sodium alginate-based non-Newtonian nanofluid over a stretching/shrinking horizontal plane surface

  • Samia Elattar , Umair Khan EMAIL logo , Aurang Zaib , Anuar Ishak , Norah Alwadai and Ahmed M. Abed
Published/Copyright: January 30, 2024

Abstract

Magnetite and cobalt ferrite (CoFe2O4) nanoparticles are frequently utilized in several applications, including magnetic drug delivery, hyperthermia, magnetic resonance imaging, etc. In the current investigation, the magnetohydrodynamic three-dimensional heat transfer (HT) flow induced by a non-Newtonian Eyring–Powell fluid is incorporated by a carrier sodium alginate (NAC6H7O6)-based CoFe2O4 nanoparticles over a deformable (stretching/shrinking) horizontal plane surface with orthogonal shear stress and power-law velocity. The HT analysis along with the substantial effect of irregular heat source/sink as well as entropy generation is also performed. The similarity variables altered the posited leading equations into ordinary differential (similarity) equations. The function bvp4c in Matlab is then used to solve these equations numerically for various parameter values. Results indicate that, in general, there are two alternative solutions for the phenomenon of suction and deformable parameters. In addition, the essential thermal evaluation is enhanced owing to the significance of CoFe2O4 nanoparticles, magnetic parameter, and irregular heat source/sink.

1 Introduction

Nanofluids (NFs) are being used in industrial and scientific processes. They make up a distinct group of nanomaterials. Compared to metals, liquids used for heat transfer (HT), such as water, engine oil, and ethylene glycol, have low thermal conductivity. A considerable increase in heat conduction can therefore be achieved by distributing solid metal particles in heat diffusion [1]. The effects of slip on the 2D flow of an NF past an axisymmetric stretchable sheet with suction were scrutinized by Gorder et al. [2]. Das [3] discovered the performance of convective heat transport (CHT) incorporated NF via a porous stretched sheet with impacts of slip and heat absorption/generation. Water-based copper and alumina nanoparticles are considered to simulate the behavior of HT fluid flow. Noghrehabadi et al. [4] studied a computational simulation of the NF flow across a porous stretchable cylinder and utilized a two-phase model subject to a magnetic field (MF). The magnetohydrodynamic (MHD) time-dependent free convective flow of water-based dissimilar nanoparticles toward a vertical permeable stretchable surface was examined by Freidoonimehr et al. [5]. Naramgari and Sulochana [6] investigated the MHD NF flow across an exponential stretching sheet in a porous media. The impacts of heat generation/absorption and thermal radiation were included in this scrutiny. The investigation concluded that the exponential parameter enhances the effects of mass and HT. Makinde et al. [7] calculated the impression of MHD on the NF flow across an irregular stretchy sheet with a chemical reaction. The impact of heat sinks/sources as well as radiation was presented to talk about the flow behavior. Acharya et al. [8] highlighted the characteristics of Hall current on the MHD hybrid NFs flow through a rotatable disc due to thermal radiation. From the analysis, it is exposed that the HT uplifts for the radiative parameter and the temperature decreases for the Hall parameter. Zaib et al. [9] inspected the buoyancy flow and heat transport progression incorporated micropolar NF across a wedge. They found that dual solutions exist in the opposing flow case only. Recently, Khan et al. [10] surveyed the influence of radiation on a buoyancy flow caused by an NF via a stretchable porous bent sheet. The slip migration, binary reaction, and activation energy are all taken into consideration to explore the specifics of water-based Al2O3 nanoparticles. Recently, Divya et al. [11] considered the influence of suction/injection on the radiative flow and HT of NF past a Riga plate with heat source/sink and Cattaneo–Christov model. They found that the velocity and the temperature declined owing to suction.

A subclass of NFs known as magnetic NFs (MNFs) is made up of colloidal suspensions of nanoscale magnetic particles that exhibit both magnetic and fluid properties. These particles are typically 5–10 nm in size. Magnetic nanoparticles are typically synthesized from metals and their oxides, such as Mn–Zn ferrite, cobalt ferrite (CoFe2O4), and iron oxide, in either polar or non-polar fluid transporters, including oil, water, and ethylene glycol in a range of shapes and sizes. The capacity of MNFs to quickly change viscosity is its main advantage [12]. A variety of applications for magnetic nanoparticles have a great deal of potential due to their special material characteristics. Magnetic nanoparticles are particularly appealing for biological applications due to their physiochemical features, including a large surface area-to-volume ratio and the possibility of surface functioning [13]. They are intriguing substances with numerous uses, such as magnetic separation, biosensors, magnetic resonance imaging, targeted drug delivery, HT/dissipation, magneto-optic sensors, etc. Khan et al. [14] observed the characteristics of the HT and SPF of CoFe2O4 nanoparticles with viscous dissipation toward a stretchable sheet. They professed that the profile of temperature reduces due to magnetic and Eckert numbers. The enhancement of HT by utilizing a non-Newtonian fluid through sodium alginate (NAC6H7O6)-based non-magnetic and magnetic particles was explored by Hussanan et al. [15]. Joshi et al. [16] assessed the influence of the high-order chemical reaction on the magneto flow of water-based hybrid NFs across a permeable bidirectional stretched sheet with volumetric heat generation. The results obtained from the current problem show that the thermal profile grows as the volume fraction increases, while the velocity decreases. The characteristic of magnetic flow induced by hybrid NFs past a stretchable sheet with the mechanism of slip effect was inspected by Pandey et al. [17]. They observed that the velocity augments due to the free convective parameter. Tadesse et al. [18] investigated the 2D Darcy–Forchheimer flow (DFF) toward a stagnation point past a shrinkable sheet induced by magnetic ferrofluids with convective heating and dissipation effects. The findings suggest that there are two possible solutions for the shrinking sheet. Tshivhi and Makinde [19] explored the influence of MF on the flow and coolant HT made of NF across a heated movable sheet. It is discovered that there are dual solutions within a particular range of shrinking sheets. Shafiq et al. [20] discussed the impression of erratic radiation on the DFF and HT incorporated by CoFe2O4 Casson fluid past a rotatable stretched sheet. Usman et al. [21] examined the water-based CoFe2O4 and copper nanoparticles by considering the theory of Darcy–Forchheimer through an exponentially stretched sheet with radiation effect. The nonlinear heat source/sink effects on an unsteady flow through a deformable sheet induced by water-based CoFe2O4 nanoparticles with porous media were scrutinized by Waini et al. [22]. Dual solutions are found for a particular strength of the shrinking parameter along with the unsteady parameter. Tshivhi and Makinde [23] examined the combined impact of Joule heating and MF on the dissipative flow of NF across a slippery and shrinking surface and obtained dual solutions. The stability analysis is also performed to check the stable solution and found that the upper solution is stable, whereas the lower solution is unstable. Recently, Murtaza et al. [24] explored the significance of CoFe2O4 nanoparticles through irregular viscosity and second-order slip past an extendable sheet in a porous media.

For a more extensive understanding of engineering and manufacturing technology issues, it is crucial to take into account research on boundary layer phenomena that focus on non-Newtonian heat-transfer substances. These substances advance widely during many industrial processes, including missile guidance, fluidization, nuclear reactor cooling, rain erosion, aerosol, and paint spraying. Although the existence of such substances is extremely problematic, numerous fundamental models/problems have been created and investigated to determine the proper flow behavior. One of the subclasses of the non-Newtonian fluid models is the Eyring–Powell fluid (EPF) model. As opposed to other non-Newtonian fluids, the model of EPF has a substantial advantage because, in place of empirical relations, it is derived from the key theory known as the “kinetic theory of gases.” More importantly, it advances at high shear rates as a viscous fluid [25]. The subject of EPF across a constantly movable porous surface with a parallel free stream has been addressed by Jalil et al. [26]. Hayat et al. [27] explored the significance of convective boundary conditions (BCs) on the flow and HT induced by EPF across a movable surface with a stable free stream. They presented an analytic solution. The features of fluid flow and thermal transport across a shrinkable surface were investigated by Roşca and Pop [28]. They performed the stability analysis to find a stable solution. Rahimi et al. [29] studied the viscous flow of a non-Newtonian EPF across a stretchable sheet and exploited the collocation technique to acquire the posited outcome. It is observed that raising the fluid material parameter causes the velocity to drop while expanding the EPF material parameter causes it to increase. Reddy et al. [30] reviewed the impression of chemical reaction on the MHD flow of EPF past a stretching slandering sheet with solar energy, solutal, thermal, and velocity slips. It has been observed that larger values of radiative HT have a greater impact on the temperature. Abbas and Megahed [31] explored the influence of radiation on the dissipative flow of EPF past a stratified surface in a porous media and obtained a numerical solution by the Chebyshev spectral technique. It is demonstrated that the temperature profile is dramatically impacted by the existence of viscous dissipation and thermal radiation in the fluid flow. Recently, Aljabali et al. [32] investigated the relationship between the motions of dusty scale particles under the impact of buoyancy convective flow of EPF. The finding suggests that the particles of dusty have an influence on fluid motion, which caused a slowdown in fluid transmission.

The current study uses a model developed by Tiwari and Das [33] to investigate the EPF flow of NAC6H7O6 (sodium alginate)-based CoFe2O4 NF via a deformable (stretchable/shrinkable) horizontal plane surface (DHPS). To fully understand the nuances of NF flow, irregular heat source/sink and magnetic effects are also accomplished. The modified system of ordinary differential equations is obtained using the similarity parameters, and it is then numerically determined using the bvp4c. In this examination, bifurcation values for the substantial constraints are also presented. The attained bifurcation values are renowned as the value at which laminar boundary-layer flow (BLF) separates from turbulent BLF. Since, the non-unique solutions or dual solutions frequently occurred due to the case of the shrinking/stretching sheet, the mass transpiration velocity and the posited unconfined vorticity inside the boundary layer. As a result, when the key outcome rate is reached, operations for some items can be designed to produce the required results, and productivity will be optimized. Therefore, it is advantageous to investigate the sodium alginate-based CoFe2O4 nanomaterials with EPF, including the controlling factors that can slow the process of flow separation and speed up the rate of HT. This endeavor is important because it will serve as a future resource for engineers, fluid mechanics experts, scientists, and practitioners. It will also serve as a foundational analysis for actual applications.

2 Mathematical formulation of the problem

Let us consider the three-dimensional incompressible MHD EPF (non-Newtonian) model with HT characteristics conveying CoFe2O4 nanoparticles in a regular NAC6H7O6-based fluid over a DHPS subjected to irregular heat source/sink and suction effects. The following equations guide the EPF model:

(1) div V = 0 ,

(2) ρ nf D V D t = p + div S ,

where V signifies the velocity vector field, ρ nf indicates the NF density, D D t represents the operator for the material time derivative, p represents the NF pressure, and S stands for the tensor of Cauchy stress. It is indicated that the constitutive equation for the NF model of EPF is as follows [34]:

(3) S = μ nf V + 1 β b Sinh 1 1 δ b V .

The final term on the right-hand side of Eq. (3) can be expressed as [34], which is the second-order estimation of the function:

(4) sinh 1 1 δ b V 1 δ b V 1 6 1 δ b V 3 , 1 δ b V 1 ,

where μ nf is the viscosity of the NF, and β b and δ b signify the material constants of the EPF model. The physical configuration of the non-Newtonian three-dimensional flow problem is exposed in Figure 1, where the geometry is designed via the Cartesian coordinate system x b , y b , and z b , which is used with corresponding velocities u b , v b , and w b . Here, x b – axis and y b – axis coordinates are reserved literally for the respective SMW (streamwise) and SNW (spanwise) directions, while the z b – axis is perpendicular to the DHPS. Assume that the velocities of both SMW and SNW at the HPS (horizontal plane surface) are represented by u wb ( x b ) = A x b m and v wb ( x b ) = B x b m , in which, A , B resemble the stretching rates [35] in the corresponding x b – and y b – directions. Moreover, w b ( x b ) = w 0 corresponds the suction/injection velocity with w 0 < 0 , and w 0 > 0 represent the cases of injection and suction, respectively, whereas the impermeable HPS is represented by w 0 = 0 . Also, it is assumed that the constant temperature located in the HPS and ambient temperature are referred to the corresponding notations such as T wb and T . A constant MF B 0 is applied normal to the HPS. The magnetic Reynolds number is treated as a very small quantity to ignore the induced MF. For thermal improvement, the CoFe2O4 nanoparticle is considered into the regular-based (sodium alginate) fluid to form the requisite NF.

Figure 1 
               Physical model and Cartesian coordinates.
Figure 1

Physical model and Cartesian coordinates.

With the help of the above-stated assumptions and the boundary layer approximations, the main simplified equations in terms of PDEs of the NF can be written as follows [32,34]:

(5) u b x b + v b y b + w b z b = 0 ,

(6) ρ nf u b u b x b + v b u b y b + w b u b z b = μ nf + 1 β b δ b 2 u b z b 2 1 2 β b δ b 3 u b z b 2 2 u b z b 2 σ nf B 0 2 u b ,

(7) ρ nf u b v b x b + v b v b y b + w b v b z b = μ nf + 1 β b δ b 2 v b z b 2 1 2 β b δ b 3 v b z b 2 2 v b z b 2 σ nf B 0 2 v b ,

(8) ( ρ c p ) nf u b T b x b + v b T b y b + w b T b z b = k nf 2 T b z b 2 + k nf u wb ( x b ) x b v nf [ A c ( T wb T ) e ξ + B c ( T b T ) ] ,

subject to the following physical BCs are

(9) u b = γ c u wb ( x b ) , v b = γ c v wb ( x b ) , w b = w 0 , T b = T wb at z b = 0 , u b 0 , v b 0 , T b T as z b .

Here, δ b and β b refer to the EPF or material constraints, T b is the temperature of the NF, A c denotes the “exponentially decaying space coefficients,” and B c is called the “time-dependent heat absorption/generation.” Additionally, the negative and positive simultaneous values of A c and B c formed the posited phenomena of heat sink and heat source, respectively. Therefore, the DHPS is denoted by γ c with γ c > 0 for “stretching case,” and γ c < 0 for “shrinking case” while the case γ c = 0 represents the static HPS.

In addition, the lasting mathematical notations used in the aforesaid equations for the thermophysical properties (TPPs) of the NF are σ nf the “electrical conductivity,” ρ nf the “density,” k nf the “TCN,” μ nf the “viscosity,” and ( ρ c p ) nf the “specific heat.” Hence, Table 1 shows the relationships of these TPPs of the NF.

Table 1

Requisite posited correlations of NAC6H7O6-based CoFe2O4 nanoparticles [18,19,36]

Properties NF correlations
Absolute viscosity μ nf = μ f ( 1 φ s c ) 2.5
Density ρ nf = ρ f φ s c ρ s c ρ f + ( 1 φ s c )
Heat capacity ( ρ c p ) nf = ( ρ c p ) f φ s c ( ρ c p ) s c ( ρ c p ) p + ( 1 φ s c )
Thermal conductivity k nf = k f ( k s c + 2 k f ) 2 φ s c ( k f k s c ) ( k s c + 2 k f ) + φ s c ( k f k s c )
Electrical conductivity σ nf = σ f 1 + 3 σ s σ f 1 σ s σ f + 2 σ s σ f 1 φ s c

Furthermore, the symbol φ s c parades the solid nanoparticle volume fraction. However, the specific condition φ s c = 0 reduces the equations of the NF model into a carrier-based NAC6H7O6 (sodium alginate) fluid. Likewise, the subscript notations like nf , s c , and f , respectively, called the NF, nanoparticles, and the carrier-based NAC6H7O6 fluid. However, c p refers to the requisite posited “heat capacitance at uniform pressure.” The experimentation data set of the regular base NAC6H7O6 fluid and the CoFe2O4 are quantitatively mentioned in Table 2.

Table 2

TPPs of the NAC6H7O6-based fluid and the CoFe2O4 nanoparticles [37,38]

Properties ρ ( kg / m 3 ) c p (J/kg K) k (W/m K) σ ( S / m ) Pr
NAC6H7O6 989 4,175 0.6376 2.60 × 10−4 6.45
CoFe2O4 4,907 700 3.7 1.1 × 107

To comfort the scrutinization of the model, therefore, we posit a similarity ansatz satisfying the above “continuity equation” (5) as follows:

(10) u b = A x b m F ( ξ ) , v b = B x b m G ( ξ ) , ξ = ( A / υ f ) 1 / 2 C ( x b ) z b , w b = A υ f C m x b m 1 F ( ξ ) + x b m C x b C { ξ F ( ξ ) F ( ξ ) } ,

where ξ refers to the “pseudo-similarity variable” and a prime indicates the differentiation by variable ξ . Also, C = C ( x b ) and C x b = d C ( x b ) d x b . Using Eq. (10), the SMW and SNW Eqs. (6) and (7) for momentum are changed into the following reduced form as follows:

(11) ρ f ρ nf μ nf μ f + 1 β b δ b μ f C 2 F ρ f ρ nf A 3 x b 2 m C 4 2 β b δ b 3 μ f υ f F 2 F σ nf / σ f ρ nf / ρ f σ f B 0 2 ρ f A F + m x b m 1 C x b C x b m F F m x b m 1 F 2 = 0 ,

(12) ρ f ρ nf μ nf μ f + 1 β b δ b μ f C 2 G ρ f ρ nf A B 2 x b 2 m C 4 2 β b δ b 3 μ f υ f G 2 G σ nf / σ f ρ nf / ρ f σ f B 0 2 ρ f A G + m x b m 1 C x b C x b m F G m x b m 1 F G = 0 .

To obtain the similarity equations (see Weidman [35] and Khan et al. [39]), the following term was adopted as follows:

(13) C x b C = K c m x b .

Now, using integration by parts, we obtain

(14) C = D c x b K c m ,

where D c and K c are requisite arbitrary constants. In addition, we put β b = β c x b 3 m 1 2 , and δ b = δ c x b 3 m 1 2 function with constants EPF parameters β c and δ c , while σ f = σ f x b m 1 2 is the base fluid electrical conductivity which is assumed as a function of x b that makes the momentum equations are dimensionless. Substituting Eqs. (13), (14), and the above functional parameters into both (SMW and SNW) Eqs. (11) and (12) yields

(15) ρ f ρ nf μ nf μ f + 1 β c δ c μ f F ρ f ρ nf 1 2 β c δ c μ f A 3 x b 2 m D c 2 x b 2 K c m δ c 2 x b 3 m 1 υ f F 2 F σ nf / σ f ρ nf / ρ f σ f B 0 2 x b m 1 ρ f A D c 2 x b 2 K c m F + m x b m 1 K c m x b m 1 D c 2 x b 2 K c m F F m x b m 1 D c 2 x b 2 K c m F 2 = 0 ,

(16) ρ f ρ nf μ nf μ f + 1 β c δ c μ f G ρ f ρ nf 1 2 β c δ c μ f A B 2 x b 2 m D c 2 x b 2 K c m δ c 2 x b 3 m 1 υ f G 2 G σ nf / σ f ρ nf / ρ f σ f B 0 2 x b m 1 ρ f A D c 2 x b 2 K c m G + m x b m 1 K c m x b m 1 D c 2 x b 2 K c m F G m x b m 1 D c 2 x b 2 K c m G F = 0 .

Thus, the coefficient of F F and the coefficient of F G are equally set to unity, so that,

(17) D c = m + 1 2 , and K c = m 1 2 m ,

and Eqs. (15) and (16) become as follows:

(18) 1 ρ nf / ρ f μ nf μ f + Σ a F m + 1 4 Σ a Σ b ρ nf / ρ f F 2 F + F F 2 m m + 1 F 2 2 m + 1 σ nf / σ f ρ nf / ρ f M c F = 0 ,

(19) 1 ρ nf / ρ f μ nf μ f + Σ a G m + 1 4 Σ a Σ b ρ nf / ρ f E c 2 G 2 G + F G 2 m m + 1 F G 2 m + 1 σ nf / σ f ρ nf / ρ f M c G = 0 .

Moreover, the distinguished parameters in the dimensionless form are defined as follows: Σ a = 1 β c δ c μ f and Σ b = A 3 δ c 2 υ f the EPF parameter, M c = σ f B 0 2 ρ f A the magnetic parameter, and E c = B A the stretching rate ratio parameter.

In addition, using Eq. (10) along with S ( ξ ) = T b T T wb T , the energy Eq. (8) is changed to the following obligatory form:

(20) C 2 Pr k nf / k f ( ρ c p ) nf / ( ρ c p ) f S + m x b m 1 C x b C x b m F S + 1 Pr ( k nf / k f ) ( ρ nf / ρ f ) x b m 1 ( ( ρ c p ) nf / ( ρ c p ) f ) ( μ nf / μ f ) [ A c e ξ + B c S ] = 0 .

Using Eqs. (13, (14), and (17), the same practice is followed to renovate Eq. (20), so that we obtain

(21) k nf k f S + Pr ( ρ c p ) nf ( ρ c p ) f F S + 2 m + 1 ρ nf ρ f k nf k f μ nf μ f [ A c e ξ + B c S ] = 0 .

Furthermore, the suction/injection factor for the HPS is given by

(22) w 0 = A υ f C m + 1 2 x b m 1 f wc ,

and the changed BCs are

(23) G ( 0 ) = 0 , F ( 0 ) = f wc , G ( 0 ) = F ( 0 ) = γ c , S ( 0 ) = 1 , G ( ξ ) 0 , F ( ξ ) 0 , S ( ξ ) 0 as ξ ,

where f wc refers to the suction/injection factor with f wc > 0 refers to the case of suction and f wc < 0 refers to the case of injection, whereas f wc = 0 refers to the impermeable HFS. Also, the Prandtl number is represented by Pr = υ f / α f .

In this examination, it is pointed ready that a given model is compacted to the “third extension of the Bank’s problem” (see Weidman [35]); if f wc = 0 (impermeable surface), M c = 0 (without MF), γ c = 1 (stretching surface), Σ a = 0 , Σ b = 0 (Newtonian model), and φ s c = 0 (carrier-based fluid) are examined. However, Weidman [35] did not include the HT characteristics because the equation of energy was left out or omitted. In addition, due to the case of viscous fluid, with no effects of the MF, permeable flat plate and Newtonian model are considered, Eqs. (18) and (23) can also be simplified to Miklavčič and Wang [40] (linear shrinking) and Fang [41] (power-law shrinking). These published works also did not explore the requisite energy equation.

The gradients of engineering physical interests are the skin friction coefficient in the SMW ( C f ) , the skin friction coefficient in the SNW ( C g ) , and the local HT rate (HTR) Nu x b , which is defined as follows:

(24) C f = 1 ρ f u wb 2 μ nf + 1 β b δ b u b z b 1 6 β b δ b 3 u b z b 3 z b = 0 , C g = 1 ρ f v wb 2 μ nf + 1 β b δ b v b z b 1 6 β b δ b 3 v b z b 3 z b = 0 , and Nu x b = x b k f ( T w b T ) k nf T b z b z b = 0 .

Using (10), (13), (14), and (17) with S ( ξ ) = T b T T wb T into the aforementioned Eq. (24), the following obtained dimensionless form of the gradients takes place:

(25) Re x b x b 1 / 2 C f = m + 1 2 μ nf μ f + Σ a F ( 0 ) m + 1 12 Σ a Σ b F ( 0 ) 3 , Re y b x b 1 / 2 C g = m + 1 2 1 E c μ nf μ f + Σ a G ( 0 ) m + 1 12 E c 2 Σ a Σ b G ( 0 ) 3 , and Re x b x b 1 / 2 Nu x b = m + 1 2 k nf k f S ( 0 ) .

Hence, Re x b x b = u wb ( x b ) x b υ f and Re y b x b = v wb ( x b ) x b υ f are the local Reynolds number.

3 Entropy analysis

The expression of volumetric entropy generation induced by Eyring–Powell NF is given by

(26) EG = k nf T 2 T b z b 2 + 1 T μ nf + 1 β b δ b u b z b + v b z b 2 1 2 β b δ b 3 u b z b + v b z b 4 + σ nf B 0 2 T ( u b 2 + v b 2 ) 2 .

The characteristics of entropy generation can be expressed as

(27) EG 0 = k f ( T wb T ) 2 x b 2 T 2 .

Similarity transformation (10) can be used to write the dimensionless form of the entropy generation as

(28) NG = EG EG 0 = k nf k f m + 1 2 Re x b x b S 2 + μ nf μ f + Σ a m + 1 2 Br b Re x b x b Ω b ( F + E c G ) 2 Σ a Σ b 2 m + 1 2 2 Br b Re x b x b Ω b ( F + E c G ) 4 + σ nf σ f M c Br b Re x b x b Ω b ( F ' 2 + E c 2 G 2 ) ,

where E c = B A , Br b = μ f u wb 2 k f ( T wb T ) , Ω b = ( T wb T ) T stand for stretching rate ratio parameter, Brinkman number, and temperature difference parameter, respectively.

4 Temporal stability analysis

We have examined the stability of dual (upper and lower) solutions as time progresses. Initially, Merkin [42] introduced the procedure of the stability analysis, and later on, Weidman et al. [43] followed. First, we begin with the newly established variables for the assessed process as follows:

(29) u b = A x b m F ( ξ , τ b ) ξ , v b = B x b m G ( ξ , τ b ) ξ , ξ = ( A / υ f ) 1 / 2 C ( x b ) z b , S ( ξ , τ b ) = T b T T wb T w b = A υ f C m x b m 1 F ( ξ , τ b ) + x b m C x b C ξ F ( ξ , τ b ) ξ F ( ξ , τ b ) ,

where τ b = A x b m 1 t a signifies the new-fangled non-dimensional time. In addition, the following can occur in the state by applying the aforementioned Eq. (29) to the time-dependent or unstable form of the necessary leading Eqs. (6)–(8) as well as Eqs. (13), (14), and (17). Thus,

(30) 1 ρ nf / ρ f μ nf μ f + Σ a 3 F ξ 3 m + 1 4 Σ a Σ b ρ nf / ρ f 2 F ξ 2 2 3 F ξ 3 + F 2 F ξ 2 2 m m + 1 F ξ 2 2 m + 1 σ nf / σ f ρ nf / ρ f M c F ξ 2 m + 1 2 F ξ τ b = 0 ,

(31) 1 ρ nf / ρ f μ nf μ f + Σ a 3 G ξ 3 m + 1 4 Σ a Σ b ρ nf / ρ f E c 2 2 G ξ 2 2 3 G ξ 3 + F 2 G ξ 2 2 m m + 1 F ξ G ξ 2 m + 1 σ nf / σ f ρ nf / ρ f M c G ξ 2 m + 1 2 G ξ τ b = 0 ,

(32) k nf / k f Pr ( ρ c p ) nf / ( ρ c p ) f 2 S ξ 2 + F S ξ + 2 m + 1 ρ nf ρ f k nf k f Pr ( ρ c p ) nf ( ρ c p ) f μ nf μ f [ A c e ξ + B c S ] 2 m + 1 S τ b ,

according to the relevant BCs:

(33) G ( 0 , τ b ) = 0 , F ( 0 , τ b ) = f wc , F ( 0 , τ b ) ξ = G ( 0 , τ b ) ξ 2 = γ c , S ( 0 , τ b ) = 1 , G ( ξ , τ b ) ξ 0 , F ( ξ , τ b ) ξ 0 , S ( ξ , τ b ) 0 as ξ .

Following that the subsequent relationships are regarded as follows (see Weidman et al. [43]):

(34) F ( ξ , τ b ) = F 0 ( ξ ) + e Γ b τ b f ( ξ ) , G ( ξ , τ b ) = G 0 ( ξ ) + e Γ b τ b g ( ξ ) , S ( ξ , τ b ) = S 0 ( ξ ) + e Γ b τ b s ( ξ ) ,

where the stable solution of F = F 0 , G = G 0 , and S = S 0 can effortlessly perturb PDEs (6)–(8). Furthermore, the functions f ( ξ ) , g ( ξ ) , and s ( ξ ) in Eq. (34) are moderately minor as associated with F 0 ( ξ ) , G 0 ( ξ ) , and S 0 ( ξ ) . The stability of the results is concluded by the sign (positive or negative) of the eigenvalue Γ b . Through setting, τ b = 0 , we have f ( ξ ) = f 0 ( ξ ) , g ( ξ ) = g 0 ( ξ ) , and s ( ξ ) = s 0 ( ξ ) . Using Eq. (34) in Eqs. (30)–(32) yields the linearized eigenvalue problem

(35) 1 ρ nf / ρ f μ nf μ f + Σ a f 0 m + 1 4 Σ a Σ b ρ nf / ρ f ( 2 F 0 F 0 f 0 + F 0 2 f 0 ) + F 0 f 0 + f 0 F 0 4 m m + 1 F 0 f 0 2 m + 1 σ nf / σ f ρ nf / ρ f M c f 0 + 2 m + 1 Γ b f 0 = 0 ,

(36) 1 ρ nf / ρ f μ nf μ f + Σ a g 0 m + 1 4 Σ a Σ b ρ nf / ρ f E c 2 ( 2 G 0 G 0 g 0 + G 0 2 g 0 ) + F 0 g 0 + f 0 G 0 2 m m + 1 ( F 0 g 0 + f 0 G 0 ) 2 m + 1 σ nf / σ f ρ nf / ρ f M c g 0 + 2 m + 1 Γ b g 0 = 0 ,

(37) k nf / k f Pr ( ρ c p ) nf / ( ρ c p ) f s 0 + F 0 s 0 + f 0 S 0 + 2 m + 1 ρ nf ρ f k nf k f Pr ( ρ c p ) nf ( ρ c p ) f μ nf μ f B c s 0 + 2 m + 1 Γ b s 0 ,

subject to BCs

(38) g 0 ( 0 ) = 0 , f 0 ( 0 ) = 0 , f 0 ( 0 ) = g 0 ( 0 ) = 0 , s 0 ( 0 ) = 0 , g 0 ( ξ ) 0 , f 0 ( ξ ) 0 , s 0 ( ξ ) 0 as ξ .

Adjusting the value of f 0 ( 0 ) , g 0 ( 0 ) or s 0 ( 0 ) yields the equivalent positive or negative values of Γ b in the aforementioned Eqs. (35), (36), and (37). Maintaining generality without losing, we set f 0 ( 0 ) = 1 and discover the eigenvalues Γ b by solving the set of related Eqs. (35)–(38) (see Khan et al. [44] and Duguma et al. [45]).

5 Numerical method

The dimensionless group of extremely difficult nonlinear Eqs. (18), (19), and (21) along with BCs (23) cannot be solved analytically or exactly. Consequently, the collocation formula has been used to numerically solve these dimensionless form equations through the fourth-order boundary value problem (bvp4c). This code, which is included in MATLAB, provides the solution for the well-known 3-stage Lobatto IIIA formula, which is likewise based on the finite difference method. The new variables must first be introduced for the bvp4c package to function. Then, we can use them to obtain the set of first-order equations from the governing dimensionless equations. Several reasonable initial guesses were needed for the software to get an improved result. Furthermore, the BCs when ξ approaches to infinity are replaced by those at ξ = ξ . Asymptotically convergence criteria and precision up to 10−6 level error tolerance are met in every instance during the iteration process. More about this procedure can be seen in the article Shampine et al. [46].

6 Analysis of the results

This section of the study shows the full physical clarifications of our entire results for all of the different contained embedded parameters. Table 3 displays a comparison of the HTR for the different values of Pr based on the accessible works when φ s c = f wc = Σ a = Σ b = M c = A c = B c = 0 , m = 1.0 , and γ c = 1.0 . The available and prior findings are nicely aligned and produce an excellent sound based on the results of the above table. The smallest numerical eigenvalues Γ b for a diverse value of γ c are displayed in Table 4. According to this table, the upper branch (UB) solution is physically and temporally feasible since the value of the least eigenvalue is positive. It is evident from the negative value of γ c that the lower branch (LB) solution is unstable and not consistent with physical reality. This perfect match can reassure us that the current executed scheme is well enough to apply and locate the unobtainable results. Figures 215 are prepared to examine the influence of the sundry physical constraints on the shear stress in the SMW, SNW, and HTR. Also, the entropy profiles are portrayed in Figures 1618. In addition, the smallest eigenvalues Γ b are also shown by the graph (Figure 19). The considered model proposes two solutions, an UB solution and an LB solution for the variance of the sundry parameters. The UB solution is denoted by solid blue lines, and the LB solution is denoted by dashed blue lines. On the other hand, the intersection of the UB and LB is commonly referred to as a bifurcation point.

Table 3

Numerical comparison of HT values for the selected sundry choices of Pr when f wc = Σ a = Σ b = M c = A c = B c = φ s c = 0 , γ c = 1.0 , and m = 1.0

Pr Devi and Devi [47] Wang [48] Khan and Pop [49] Present results
2.0 0.91135 0.9114 0.9113 0.9113543
6.13 1.75968 1.7596836
7.0 1.89540 1.8954 1.8954 1.8954001
20.0 3.35390 3.3539 3.3539 3.3539025
Table 4

Smallest eigenvalues Γ b for a diverse value of γ c when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0

γ c Smallest eigenvalues γ c
Upper solution Lower solution
‒0.3 1.9435 −1.6435
−0.5 1.8317 −1.5345
−0.8 1.6673 −1.3672
−1.0 1.3534 −1.1427
−1.05 1.2234 −1.0584
−1.1 1.0234 −0.9145
−1.12 0.9367 −0.8246
−1.14 0.8466 −0.7657
−1.16 0.7534 −0.6567
−1.18 0.6393 −0.5389
−1.2 0.4534 −0.4358
−1.21 0.3638 −0.2167
−1.214 0.0145 −0.0250
Figure 2 
               Shear stress in the SMW direction versus 
                  
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                        
                        {f}_{\text{wc}}
                     
                   for distinct values of 
                     
                        
                        
                           m
                        
                        m
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , and 
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                           =
                           ‒
                           1.0
                        
                        {\gamma }_{\text{c}}=&#x2012;1.0
                     
                  .
Figure 2

Shear stress in the SMW direction versus f wc for distinct values of m when Σ a = 0.1 , Σ b = 0.5 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , φ s c = 0.025 , and γ c = 1.0 .

Figure 3 
               Shear stress in the SNW direction versus 
                  
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                        
                        {f}_{\text{wc}}
                     
                   for distinct values of 
                     
                        
                        
                           m
                        
                        m
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , and 
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                           =
                           ‒
                           1.0
                        
                        {\gamma }_{\text{c}}=&#x2012;1.0
                     
                  .
Figure 3

Shear stress in the SNW direction versus f wc for distinct values of m when Σ a = 0.1 , Σ b = 0.5 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , φ s c = 0.025 , and γ c = 1.0 .

Figure 4 
               Variation of HTR versus 
                  
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                        
                        {f}_{\text{wc}}
                     
                   for distinct values of 
                     
                        
                        
                           m
                        
                        m
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , and 
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                           =
                           ‒
                           1.0
                        
                        {\gamma }_{\text{c}}=&#x2012;1.0
                     
                  .
Figure 4

Variation of HTR versus f wc for distinct values of m when Σ a = 0.1 , Σ b = 0.5 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , φ s c = 0.025 , and γ c = 1.0 .

Figure 5 
               Shear stress in the SMW direction versus 
                  
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                        
                        {f}_{\text{wc}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                        
                        {\Sigma }_{\text{b}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , and 
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                           =
                           ‒
                           1.0
                        
                        {\gamma }_{\text{c}}=&#x2012;1.0
                     
                  .
Figure 5

Shear stress in the SMW direction versus f wc for distinct values of Σ b when Σ a = 0.1 , m = 2.0 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , φ s c = 0.025 , and γ c = 1.0 .

Figure 6 
               Shear stress in the SNW direction versus 
                  
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                        
                        {f}_{\text{wc}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                        
                        {\Sigma }_{\text{b}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , and 
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                           =
                           ‒
                           1.0
                        
                        {\gamma }_{\text{c}}=&#x2012;1.0
                     
                  .
Figure 6

Shear stress in the SNW direction versus f wc for distinct values of Σ b when Σ a = 0.1 , m = 2.0 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , φ s c = 0.025 , and γ c = 1.0 .

Figure 7 
               Variation of HTR versus 
                  
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                        
                        {f}_{\text{wc}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                        
                        {\Sigma }_{\text{b}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , and 
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                           =
                           ‒
                           1.0
                        
                        {\gamma }_{\text{c}}=&#x2012;1.0
                     
                  .
Figure 7

Variation of HTR versus f wc for distinct values of Σ b when Σ a = 0.1 , m = 2.0 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , φ s c = 0.025 , and γ c = 1.0 .

Figure 8 
               Shear stress in the SMW direction versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                        
                        {\varphi }_{{s}_{c}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 8

Shear stress in the SMW direction versus γ c for distinct values of φ s c when Σ a = 0.1 , m = 2.0 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 9 
               Shear stress in the SNW direction versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                        
                        {\varphi }_{{s}_{c}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 9

Shear stress in the SNW direction versus γ c for distinct values of φ s c when Σ a = 0.1 , m = 2.0 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 10 
               Variation of HTR versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                        
                        {\varphi }_{{s}_{c}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 10

Variation of HTR versus γ c for distinct values of φ s c when Σ a = 0.1 , m = 2.0 , M c = 0.1 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 11 
               Shear stress in the SMW direction versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                        
                        {M}_{\text{c}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 11

Shear stress in the SMW direction versus γ c for distinct values of M c when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 12 
               Shear stress in the SNW direction versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                        
                        {M}_{\text{c}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 12

Shear stress in the SNW direction versus γ c for distinct values of M c when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 13 
               Variation of HTR versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                        
                        {M}_{\text{c}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 13

Variation of HTR versus γ c for distinct values of M c when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 14 
               Variation of HTR versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           ,
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           >
                           0
                        
                        {A}_{\text{c}}^{\ast },{B}_{\text{c}}^{\ast }\gt 0
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  
                  
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 14

Variation of HTR versus γ c for distinct values of A c , B c > 0 when Σ a = 0.1 m = 2.0 , φ s c = 0.025 , E c = 0.5 , M c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 15 
               Variation of HTR versus 
                  
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                   for distinct values of 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           ,
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           <
                           0
                        
                        {A}_{\text{c}}^{\ast },{B}_{\text{c}}^{\ast }\lt 0
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                   
                  
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 15

Variation of HTR versus γ c for distinct values of A c , B c < 0 when Σ a = 0.1 m = 2.0 , φ s c = 0.025 , E c = 0.5 , M c = 0.1 , Σ b = 0.5 , and f wc = 2.0 .

Figure 16 
               Variation of entropy for distinct values of 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                        
                        {M}_{\text{c}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 Br
                              
                              
                                 b
                              
                           
                           =
                           1.0
                        
                        {\text{Br}}_{\text{b}}=1.0
                     
                  , 
                     
                        
                        
                           
                              
                                 Ω
                              
                              
                                 b
                              
                           
                           =
                           0.01
                        
                        {\Omega }_{\text{b}}=0.01
                     
                  , 
                     
                        
                        
                           
                              
                                 Re
                              
                              
                                 
                                    
                                       x
                                    
                                    
                                       b
                                    
                                 
                                 
                                    
                                       x
                                    
                                    
                                       b
                                    
                                 
                              
                           
                           =
                           10
                        
                        {\mathrm{Re}}_{{x}_{\text{b}}{x}_{\text{b}}}=10
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 16

Variation of entropy for distinct values of M c when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , Br b = 1.0 , Ω b = 0.01 , Re x b x b = 10 , and f wc = 2.0 .

Figure 17 
               Variation of entropy for distinct values of 
                     
                        
                        
                           
                              
                                 Br
                              
                              
                                 b
                              
                           
                        
                        {\text{Br}}_{\text{b}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Ω
                              
                              
                                 b
                              
                           
                           =
                           0.01
                        
                        {\Omega }_{\text{b}}=0.01
                     
                  , 
                     
                        
                        
                           
                              
                                 Re
                              
                              
                                 
                                    
                                       x
                                    
                                    
                                       b
                                    
                                 
                                 
                                    
                                       x
                                    
                                    
                                       b
                                    
                                 
                              
                           
                           =
                           10
                        
                        {\mathrm{Re}}_{{x}_{\text{b}}{x}_{\text{b}}}=10
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 17

Variation of entropy for distinct values of Br b when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , M c = 0.1 , Ω b = 0.01 , Re x b x b = 10 , and f wc = 2.0 .

Figure 18 
               Variation of entropy for distinct values of 
                     
                        
                        
                           
                              
                                 Ω
                              
                              
                                 b
                              
                           
                        
                        {\Omega }_{\text{b}}
                     
                   when 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 a
                              
                           
                           =
                           0.1
                        
                        {\Sigma }_{\text{a}}=0.1
                     
                  , 
                     
                        
                        
                           m
                           =
                           2.0
                        
                        m=2.0
                     
                  , 
                     
                        
                        
                           
                              
                                 φ
                              
                              
                                 
                                    
                                       s
                                    
                                    
                                       c
                                    
                                 
                              
                           
                           =
                           0.025
                        
                        {\varphi }_{{\text{s}}_{\text{c}}}=0.025
                     
                  , 
                     
                        
                        
                           
                              
                                 E
                              
                              
                                 c
                              
                           
                           =
                           0.5
                        
                        {E}_{\text{c}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 A
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           
                              
                                 B
                              
                              
                                 c
                              
                              
                                 ⁎
                              
                           
                           =
                           0.1
                        
                        {A}_{\text{c}}^{\ast }={B}_{\text{c}}^{\ast }=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Σ
                              
                              
                                 b
                              
                           
                           =
                           0.5
                        
                        {\Sigma }_{\text{b}}=0.5
                     
                  , 
                     
                        
                        
                           
                              
                                 Br
                              
                              
                                 b
                              
                           
                           =
                           1.0
                        
                        {\text{Br}}_{\text{b}}=1.0
                     
                  , 
                     
                        
                        
                           
                              
                                 M
                              
                              
                                 c
                              
                           
                           =
                           0.1
                        
                        {M}_{\text{c}}=0.1
                     
                  , 
                     
                        
                        
                           
                              
                                 Re
                              
                              
                                 
                                    
                                       x
                                    
                                    
                                       b
                                    
                                 
                                 
                                    
                                       x
                                    
                                    
                                       b
                                    
                                 
                              
                           
                           =
                           10
                        
                        {\mathrm{Re}}_{{x}_{\text{b}}{x}_{\text{b}}}=10
                     
                  , and 
                     
                        
                        
                           
                              
                                 f
                              
                              
                                 wc
                              
                           
                           =
                           2.0
                        
                        {f}_{\text{wc}}=2.0
                     
                  .
Figure 18

Variation of entropy for distinct values of Ω b when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , Br b = 1.0 , M c = 0.1 , Re x b x b = 10 , and f wc = 2.0 .

Figure 19 
               Variation of smallest eigenvalues 
                     
                        
                        
                           
                              
                                 Γ
                              
                              
                                 b
                              
                           
                        
                        {\Gamma }_{\text{b}}
                     
                   for a diverse value of 
                     
                        
                        
                           
                              
                                 γ
                              
                              
                                 c
                              
                           
                        
                        {\gamma }_{\text{c}}
                     
                  .
Figure 19

Variation of smallest eigenvalues Γ b for a diverse value of γ c .

Figures 24 illustrate the influence of m on the shear stress coefficient in the SMW direction, shear stress coefficient in the SNW direction, and HTR of the NAC6H7O6-based CoFe2O4 nanoparticles, respectively. The UB and LB outcomes are created in these graphs for the circumstance of the suction parameter. The products of friction factor in both directions of the UBS for NAC6H7O6-based CoFe2O4 nanoparticles improve with the remarkable augmentation of m versus mass suction parameter f wc . However, opposite monotonic behaviors are established for the LBS with superior consequences of m . Besides, the condition is noted that m is the power exponent constraint for the momentum distribution, whereas m = 1 and m > 1 refer to the instance of linear and nonlinear (power-law) stretching/shrinking velocity, respectively. As proved in Figure 4, the higher power-law index parameter value upsurges the HTR in both (UB and LB) outcomes. Additionally, the rate of HT is higher for the power-law index parameter values m = 2 , 3 as compared to the linear case m = 1 . Also, the solid gray ball in all three above-stated graphs denotes the point of bifurcation where the UB and LB solutions are merged. At this particular point, the outcomes of the problem are unique and mathematically it is denoted by f wc = f wc C , where f wc C is the bifurcation or critical point. This bifurcation point will be changed via changing the values of m , see the located graphs. Meanwhile, the outcomes are not unique for the case of f wc > f wc C , and no solutions exist for the case of f wc < f wc C when we do the variations in the power-law index parameter f wc = f wc C . With higher values of m , the absolute or magnitude of the bifurcation values are augmented. This design highlights that the BL separation decline from the HPS with greater impacts of m .

The effect of Σ b on the friction factor coefficients in the directions of SMW and SNW of the NAC6H7O6-based CoFe2O4 nanoparticles for the UBS and LBS are presented in Figures 5 and 6. From these graphical structures, it is recognized that the friction factor coefficients in both directions decline for the UBS with growing the material or non-Newtonian (EPF) parameter Σ b but a notable augmentation is observed in the LBS. On the other hand, the HTR with rising values of Σ b for the outcomes of UB as well as LB is shown in Figure 7. Superior Σ b central to a considerable improvement in the values of HTR for the LBS, whereas retardation is noticed for the branch of the US. Generally, the absolute viscosity μ f and the material parameter Σ b follow the inverse relation. Therefore, a significant increase in the value of Σ b produces a reduction in the value of μ f , and as a consequence, the velocity profile is boosted. Thus, the friction factor coefficients in both routes of the NAC6H7O6-based CoFe2O4 nanoparticles should, therefore, decrease as the velocity increases. In addition, the bifurcation value developments in magnitude or absolute due to stronger impacts of the material parameter Σ b . Thus, a delay in the separation of the BLF from the HPS is perceived because of the higher consequence in the values of Σ b .

Figures 8 and 9 illustrate the outcome of the friction factor coefficients in both directions of the regular NAC6H7O6-based CoFe2O4 nanoparticles for the UBS and LBS with growing values of φ s c . The outcome in both graphs shows that the shear stress for the UBS declines in the situation of stretching and escalates in the situation of shrinking due to the superior effects of φ s c . However, continuous reductions are seen for the LBS in the case of deformable with increasing φ s c . Figure 10 exhibits the influences of φ s c on the HTR of the NAC6H7O6-based CoFe2O4 nanoparticles for the UB and LB solutions. It is professed that the outcome of HTR escalates owing to the advanced values of φ s c . Generally, the EPF flow TCN develops with higher impacts of φ s c , as a fallout, the HT upsurges. Additionally, two separate (UB and LB) branches meet at a single location known as a critical point. At this position, the results are exceptional (or unique) and scientifically it is expressed as γ c = γ c C , where γ c C corresponds the critical point. For the specific domain areas such as γ c C < γ c < , and < γ c < γ c C , the dual outcomes and no outcomes exist, respectively (see Elatter et al. [50]). Owing to greater impacts of φ s c , the bifurcation or critical values in magnitude declines. This behavior is triggered due to the stronger effects of φ s c , and hence the rise of the flow separation from the HPS.

The consequence of the magnetic parameter M c on the friction factor coefficients in both directions and HTR of the NAC6H7O6-based CoFe2O4 nanoparticles for the UB and LB solutions are revealed in Figures 1113, respectively. For rising M c , the values of friction factor coefficients in both paths upsurge for the UBS, while the trends are reversed for the LBS. More especially, the outcomes for the UBS exist if we go forward from shrinking to stretching plane surface, but for the case of stretching, the LBS diminished due to higher value of M c . However, the HTR is also enriching for the UB and as well as for the LB results owing to the superior impacts of M c . Physically, the Lorentz force operates along the fluid flow due to the implementation of an MF in the laminar flow condition. The Lorentz force generates internal resistance within the boundary layer, resulting in a dragging force. Therefore, the fluid particles undergo a reduction in momentum due to the consequence of the Lorentz force, as a result, the velocity profiles decline. Also, the friction force and the velocity follow the basic law of physics called a direct relation. According to this relation, the friction factor coefficients intensify in both routes. Besides, the magnitude of the bifurcation values uplifts with higher magnetic parameter and hence delays the BL flow separation from the HPS.

Figures 14 and 15 portray the impact of the internal heat source parameter, A c , B c > 0 , and heat sink constraint, A c , B c < 0 , on the HTR of the NAC6H7O6)-based CoFe2O4 nanoparticles for the dual branches, respectively. The greater impact of A c , B c > 0 leads to a substantial decline in the UB as well as the LB outcomes. However, a very opposite pattern or tendency is realized in Figure 15 when we increase the values of A c , B c < 0 . Owing to alteration the simultaneous positive and negative values of A c and B c , as a result, the critical values are the same (Figures 14 and 15). In light of this, the difference between the UB and LB solution curves is preferable for the situation of a heat source as opposed to the effect of a heat sink parameter.

Figures 1618 show the behavior of the entropy generation for both solution branches due to the different values of M c , Br b , and Ω b . These profiles suggest that the entropy accelerates by either improving magnetic or Brinkman numbers. Meanwhile, it is decelerated due to the higher impact of the temperature difference parameter Ω b . Physically, entropy is generated by all irreversibility mechanisms, and as a result, the Lorentz forces cause entropy to increase. Increasing the Brinkman number leads to greater amounts of entropy generation created by the irreversibility of the fluid friction.

Figure 19 shows the smallest eigenvalues Γ b against the shrinking/stretching parameter γ c for the particular selected magnetic parameter value such as M c = 0.10 when Σ a = 0.1 , m = 2.0 , φ s c = 0.025 , E c = 0.5 , A c = B c = 0.1 , Σ b = 0.5 , and f wc = 2.0 . From the graph, it is seen that for positive values of Γ b the term e Γ b τ b 0 as Γ b and it explodes for higher negative values. In other words, the sign of the smallest eigenvalue Γ b for the first solution is positive, which signifies the stability characteristics of the solutions. Meanwhile, the values of smallest eigenvalue Γ b for the second solution are negative, which indicates the second solution long-term instability. As a result, the upper solution is physically more trustworthy than the lower solution since its stability reflects its physical reliability.

7 Final remarks

This article aims to investigate the NAC6H7O6-based CoFe2O4 nanoparticles on magneto flow and HT toward a nonlinear DHPS with irregular heat source/sink and power-law velocity. The dual or binary (UB and LB) solutions were obtained using the bvp4c solver. The reliability and validity of results were tested through stability analysis in a few instances where it was determined that the current results were equivalent to the prior ones. Numerous portraits have been plotted to examine the effects of the various relevant constraints on the shear stress in the respective SMW and SNW directions, rate of HT and entropy generation profiles. We can draw the following important conclusions from this investigation as follows:

  • Multiple solutions are provoked only if the strength of suction is applied and in the region of the shrinking sheet.

  • The shear stress in both directions as well as the HTR augment due to the power-law index m and the suction parameter.

  • The shear stress in both directions and HTR decline due to the EPF parameter.

  • The sodium alginate-based CoFe2O4 nanoparticles boost the HTR in both resolutions while augmenting the shear stress in SMW and SNW directions in the UB solution and decline in the LB solution.

  • The presence of a magnetic parameter ultimately augments the shear stress in the UB solution and declines in the LB solution, whereas the HTR uplifts in both solutions.

  • The HTR declines in the presence of source parameter and uplift due to sink parameter in both UB and LB solutions.

  • The entropy generation accelerates due to magnetic and Brinkman numbers, while the opposite trend is observed due to the temperature difference parameter.

  • From the temporal stability analysis, it is recognized that the UB solutions are stable and physically acceptable.



Acknowledgments

This work was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University, through the Research Groups Program Grant no. (RGP-1444-0060).

  1. Funding information: This work was funded by the Deanship of Scientific Research at Princess Nourah bint Abdulrahman University, through the Research Groups Program Grant no. (RGP-1444-0060).

  2. Author contributions: Conceptualization, S.E., and U.K.; methodology, S.E., A.I., and U.K.; software, A.I., A.Z., S.E., and U.K.; validation, N.A., A.M.A., U.K., A.Z., and A.I.; formal analysis, A.I., A.M.A., N.A., and A.Z.; investigation, A.I., and A.Z.; resources, A.I.; data curation, A.I., A.Z., U.K., and S.E., writing – original draft preparation, N.A., U.K., and A.M.A.; writing – review and editing, N.A., A.Z., and A.M.A.; visualization, A.Z.; supervision, A.I.; project administration, N.A., A.M.A. and S.E.; funding acquisition, N.A., A.M.A. and S.E. All authors have accepted responsibility for the entire content of this manuscript and approved its submission.

  3. Conflict of interest: The authors state no conflict of interest.

  4. Data availability statement: The datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

References

[1] Nima NI, Salawu SO, Ferdows M, Shamshuddin MD, Alsenafi A, Nakayama A. Melting effect on non‐Newtonian fluid flow in gyrotactic microorganism saturated non‐Darcy porous media with variable fluid properties. Appl Nanosci. 2020;10:3911–24.10.1007/s13204-020-01491-ySearch in Google Scholar

[2] Gorder RAV, Sweet E, Vajravelu K. Nano boundary layers over stretching surfaces. Commun Nonlinear Sci Numer Simul. 2010;15(6):1494–500.10.1016/j.cnsns.2009.06.004Search in Google Scholar

[3] Das K. Slip flow and convective heat transfer of nanofluids over a permeable stretching surface. Comput Fluids. 2012;64:34–42.10.1016/j.compfluid.2012.04.026Search in Google Scholar

[4] Noghrehabadi A, Ghalambaz M, Izadpanahi E, Pourrajab R. Effect of magnetic field on the boundary layer flow, heat, and mass transfer of nanofluids over a stretching cylinder. J Heat Mass Transf Res. 2014;1(1):9–16.Search in Google Scholar

[5] Freidoonimehr N, Rashidi MM, Mahmud S. Unsteady MHD free convective flow past a permeable stretching vertical surface in a nano-fluid. Int J Therm Sci. 2015;87:136–45.10.1016/j.ijthermalsci.2014.08.009Search in Google Scholar

[6] Naramgari S, Sulochana C. Dual solutions of radiative MHD nanofluid flow over an exponentially stretching sheet with heat generation/absorption. Appl Nanosci. 2016;6(1):131–9.10.1007/s13204-015-0420-zSearch in Google Scholar

[7] Makinde OD, Mabood F, Ibrahim SM. Chemically reacting on MHD boundary layer flow of nanofluids over a nonlinear stretching sheet with heat source/sink and thermal radiation. Therm Sci. 2018;22:495–506.10.2298/TSCI151003284MSearch in Google Scholar

[8] Acharya N, Bag R, Kundu PK. Influence of Hall current on radiative nanofluid flow over a spinning disk: A hybrid approach. Phys E: Low-dimens Syst Nanostruct. 2019;111:103–12.10.1016/j.physe.2019.03.006Search in Google Scholar

[9] Zaib A, Haq RU, Sheikholeslami M, Khan U. Numerical analysis of effective Prandtl model on mixed convection flow of γAl2O3–H2O nanoliquids with micropolar liquid driven through wedge. Phys Scr. 2020;95:035005.10.1088/1402-4896/ab5558Search in Google Scholar

[10] Khan U, Zaib A, Madhukesh JK, Elattar S, Eldin SM, Ishak A, et al. Features of radiative mixed convective heat transfer on the slip flow of nanofluid past a stretching bended sheet with activation energy and binary reaction. Energies. 2022;15:7613.10.3390/en15207613Search in Google Scholar

[11] Divya S, Eswaramoorthi S, Loganathan K. Numerical computation of Ag/Al2O3 nanofluid over a Riga Plate with heat sink/source and non-Fourier heat flux model. Math Comput Appl. 2023;28(1):20.10.3390/mca28010020Search in Google Scholar

[12] Genc S, Derin B. Synthesis and rheology of ferrofluids: A review. Curr Opin Chem Eng. 2014;3:118–24.10.1016/j.coche.2013.12.006Search in Google Scholar

[13] Issa B, Obaidat IM, Albiss BA, Haik Y. Magnetic nanoparticles: Surface effects and properties related to biomedicine applications. Int J Mol Sci. 2013;14(11):21266–305.10.3390/ijms141121266Search in Google Scholar PubMed PubMed Central

[14] Khan ZH, Khan WA, Qasim M, Shah IA. MHD stagnation point ferrofluid flow and heat transfer toward a stretching sheet. IEEE Trans Nanotechnol. 2014;13(1):35–40.10.1109/TNANO.2013.2286991Search in Google Scholar

[15] Hussanan A, Qasim M, Chen Z-M. Heat transfer enhancement in sodium alginate based magnetic and non-magnetic nanoparticles mixture hybrid nanofluid. Phys A. 2020;550:123957.10.1016/j.physa.2019.123957Search in Google Scholar

[16] Joshi N, Upreti H, Pandey AK, Kumar M. Heat and mass transfer assessment of magnetic hybrid nanofluid flow via bidirectional porous surface with volumetric heat generation. Int J Appl Comput Math. 2021;7:64.10.1007/s40819-021-00999-3Search in Google Scholar

[17] Pandey AK, Upreti H, Joshi N, Uddin Z. Effect of natural convection on 3D MHD flow of MoS2–GO/H2O via porous surface due to multiple slip mechanisms. J Taibah Univ Sci. 2022;16(1):749–62.10.1080/16583655.2022.2113729Search in Google Scholar

[18] Tadesse FB, Makinde OD, Enyadene LG. Hydromagnetic stagnation point flow of a magnetite ferrofluid past a convectively heated permeable stretching/shrinking sheet in a Darcy–Forchheimer porous medium. Sādhanā. 2021;46:115.10.1007/s12046-021-01643-ySearch in Google Scholar

[19] Tshivhi KS, Makinde OD. Magneto-nanofluid coolants past heated shrinking/stretching surfaces: Dual solutions and stability analysis. Results Eng. 2021;10:100229.10.1016/j.rineng.2021.100229Search in Google Scholar

[20] Shafiq A, Lone SA, Sindhu TN, Nonlaopon K. Statistical modelling for the Darcy–Forchheimer flow of Casson cobalt ferrite-water/ethylene glycol nanofluid under nonlinear radiation. Symmetry. 2022;14:1717.10.3390/sym14081717Search in Google Scholar

[21] Usman M, Amin S, Saeed A. Magnetohydrodynamic hybrid nanofluid flow with the effect of Darcy–Forchheimer theory and slip conditions over an exponential stretchable sheet. Adv Mech Eng. 2022;14(8):1–15.10.1177/16878132221116479Search in Google Scholar

[22] Waini I, Khan U, Zaib A, Ishak A, Pop I, Akkurt N. Time-dependent flow of water-based CoFe2O4-Mn-ZnFe2O4 nanoparticles over a shrinking sheet with mass transfer effect in porous media. Nanomaterials. 2022;12:4102.10.3390/nano12224102Search in Google Scholar PubMed PubMed Central

[23] Tshivhi KS, Makinde OD. Dual solutions and stability analysis for Buongiorno model of magnetohydrodynamics nanofluid flow past a heated shrinking slippery surface. J Nanofluids. 2023;12(6):1485–94.10.1166/jon.2023.2032Search in Google Scholar

[24] Murtaza S, Kumam P, Bilal M, Sutthibutpong T, Rujisamphan N, Ahmad Z. Parametric simulation of hybrid nanofluid flow consisting of cobalt ferrite nanoparticles with second-order slip and variable viscosity over an extending surface. Nanotechnol Rev. 2023;12:20220533.10.1515/ntrev-2022-0533Search in Google Scholar

[25] Powell RE, Eyring H. Mechanism for relaxation theory of viscosity. Nat Lond. 1944;154:427–8.10.1038/154427a0Search in Google Scholar

[26] Jalil M, Asghar S, Imran SM. Self similar solutions for the flow and heat transfer of Powell–Eyring fluid over a moving surface in a parallel free stream. Int J Heat Mass Transf. 2013;65:73–9.10.1016/j.ijheatmasstransfer.2013.05.049Search in Google Scholar

[27] Hayat T, Iqbal Z, Qasim M, Obaidat S. Steady flow of an Eyring Powell fluid over a moving surface with convective boundary conditions. Int J Heat Mass Transf. 2012;55:1817–22.10.1016/j.ijheatmasstransfer.2011.10.046Search in Google Scholar

[28] Roşca AV, Pop I. Flow and heat transfer of Powell–Eyring fluid over a shrinking surface in a parallel free stream. Int J Heat Mass Transf. 2014;71:321–7.10.1016/j.ijheatmasstransfer.2013.12.020Search in Google Scholar

[29] Rahimi J, Ganji DD, Khaki M, Hosseinzadeh Kh. Solution of the boundary layer flow of an Eyring-Powell non-Newtonian fluid over a linear stretching sheet by collocation method. Alex Eng J. 2017;56(4):621–7.10.1016/j.aej.2016.11.006Search in Google Scholar

[30] Reddy SRR, Reddy PBA, Bhattacharyya K. Effect of nonlinear thermal radiation on 3D magneto slip flow of Eyring-Powell nanofluid flow over a slendering sheet with binary chemical reaction and Arrhenius activation energy. Adv Powder Technol. 2019;30:3203–13.10.1016/j.apt.2019.09.029Search in Google Scholar

[31] Abbas W, Megahed AM. Powell-Eyring fluid flow over a stratified sheet through porous medium with thermal radiation and viscous dissipation. AIMS Math. 2021;6(12):13464–79.10.3934/math.2021780Search in Google Scholar

[32] Aljabali A, Kasim ARM, Arifin NS, Ariffin NAN, Ching DLC, Waini I, et al. Two-phase flow of Eyring–Powell fluid with temperature dependent viscosity over a vertical stretching sheet. Mathematics. 2022;10:3111.10.3390/math10173111Search in Google Scholar

[33] Tiwari RK, Das MK. Heat transfer augmentation in a two-sided lid-driven differentially heated square cavity utilizing nanofluids. Int J Heat Mass Transf. 2007;50:2002–18.10.1016/j.ijheatmasstransfer.2006.09.034Search in Google Scholar

[34] Patel M, Timol MG. Numerical treatment of Powell-Eyring fluid flow using method of satisfaction of asymptotic boundary conditions (MSABC). Appl Numer Math. 2009;59:2584–92.10.1016/j.apnum.2009.04.010Search in Google Scholar

[35] Weidman P. Flows induced by power-law stretching surface motion modulated by transverse or orthogonal surface shear. C R Mѐc. 2017;345(2):169–76.10.1016/j.crme.2016.10.016Search in Google Scholar

[36] Upreti H, Pandey AK, Kumar M, Makinde OD. Ohmic heating and non-uniform heat source/sink roles on 3D Darcy–Forchheimer flow of CNTs nanofluids over a stretching surface. Arab J Sci Eng. 2020;45:7705–17.10.1007/s13369-020-04826-7Search in Google Scholar

[37] Alghamdi M, Wakif A, Thumma T, Khan U, Baleanu D, Rasool G. Significance of variability in magnetic field strength and heat source on the radiative-convective motion of sodium alginate-based nanofluid within a Darcy-Brinkman porous structure bounded vertically by an irregular slender surface. Case Stud Therm Eng. 2021;28:101428.10.1016/j.csite.2021.101428Search in Google Scholar

[38] Sadighi S, Afshar H, Jabbari M, Ashtiani HAD. Heat and mass transfer for MHD nanofluid flow on a porous stretching sheet with prescribed boundary conditions. Case Stud Therm Eng. 2023;49:103345.10.1016/j.csite.2023.103345Search in Google Scholar

[39] Khan U, Zaib A, Pop I, Waini I, Ishak A. MHD flow of a nanofluid due to a nonlinear stretching/shrinking sheet with a convective boundary condition: Tiwari–Das nanofluid model. Int J Numer Methods Heat Fluid Flow. 2022;32(10):3233–58.10.1108/HFF-11-2021-0730Search in Google Scholar

[40] Miklavčič M, Wang C. Viscous flow due to a shrinking sheet. Q Appl Math. 2006;64(2):283–90.10.1090/S0033-569X-06-01002-5Search in Google Scholar

[41] Fang T. Boundary layer flow over a shrinking sheet with power-law velocity. Int J Heat Mass Transf. 2008;51(25–26):5838–43.10.1016/j.ijheatmasstransfer.2008.04.067Search in Google Scholar

[42] Merkin JH. Mixed convection boundary layer flow on a vertical surface in a saturated porous medium. J Eng Math. 1980;14(4):301–13.10.1007/BF00052913Search in Google Scholar

[43] Weidman PD, Kubitschek DG, Davis AMJ. The effect of transpiration on self-similar boundary layer flow over moving surfaces. Int J Eng Sci. 2006;44(11–12):730–7.10.1016/j.ijengsci.2006.04.005Search in Google Scholar

[44] Khan U, Zaib A, Ishak A, Alotaibi AM, Eldin SM, Akkurt N, et al. Stability analysis of buoyancy magneto flow of hybrid nanofluid through a stretchable/shrinkable vertical sheet induced by a micropolar fluid subject to nonlinear heat sink/source. Magnetochemistry. 2022;8(12):188.10.3390/magnetochemistry8120188Search in Google Scholar

[45] Duguma KA, Makinde OD, Enyadene LG. Stability analysis of dual solutions of convective flow of Casson nanofluid past a shrinking/stretching slippery sheet with thermophoresis and Brownian motion in porous media. J Math. 2023;2023:5954860.10.1155/2023/5954860Search in Google Scholar

[46] Shampine LF, Gladwell I, Thompson S. Solving ODEs with MATLAB. Cambridge: Cambridge University Press; 2003.10.1017/CBO9780511615542Search in Google Scholar

[47] Devi SSU, Devi SPA. Heat transfer enhancement of Cu-Al2O3/water hybrid nanofluid flow over a stretching sheet. J Niger Math Soc. 2017;36:419–33.Search in Google Scholar

[48] Wang CY. Free convection on a vertical stretching surface. J Appl Math Mech (ZAMM). 1989;69:418–20.10.1002/zamm.19890691115Search in Google Scholar

[49] Khan W, Pop I. Boundary-layer flow of a nanofluid past a stretching sheet. Int J Heat Mass Transf. 2010;53:2477–83.10.1016/j.ijheatmasstransfer.2010.01.032Search in Google Scholar

[50] Elatter S, Khan U, Zaib A, Ishak A, Saleh W, Abed AM. Scrutinization of waste discharge concentrations in Eyring-Powell nanofluid past a deformable horizontal plane surface. Water. 2023;15(19):3419.10.3390/w15193419Search in Google Scholar

Received: 2023-11-08
Revised: 2024-01-02
Accepted: 2024-01-04
Published Online: 2024-01-30

© 2024 the author(s), published by De Gruyter

This work is licensed under the Creative Commons Attribution 4.0 International License.

Articles in the same Issue

  1. Regular Articles
  2. Numerical study of flow and heat transfer in the channel of panel-type radiator with semi-detached inclined trapezoidal wing vortex generators
  3. Homogeneous–heterogeneous reactions in the colloidal investigation of Casson fluid
  4. High-speed mid-infrared Mach–Zehnder electro-optical modulators in lithium niobate thin film on sapphire
  5. Numerical analysis of dengue transmission model using Caputo–Fabrizio fractional derivative
  6. Mononuclear nanofluids undergoing convective heating across a stretching sheet and undergoing MHD flow in three dimensions: Potential industrial applications
  7. Heat transfer characteristics of cobalt ferrite nanoparticles scattered in sodium alginate-based non-Newtonian nanofluid over a stretching/shrinking horizontal plane surface
  8. The electrically conducting water-based nanofluid flow containing titanium and aluminum alloys over a rotating disk surface with nonlinear thermal radiation: A numerical analysis
  9. Growth, characterization, and anti-bacterial activity of l-methionine supplemented with sulphamic acid single crystals
  10. A numerical analysis of the blood-based Casson hybrid nanofluid flow past a convectively heated surface embedded in a porous medium
  11. Optoelectronic–thermomagnetic effect of a microelongated non-local rotating semiconductor heated by pulsed laser with varying thermal conductivity
  12. Thermal proficiency of magnetized and radiative cross-ternary hybrid nanofluid flow induced by a vertical cylinder
  13. Enhanced heat transfer and fluid motion in 3D nanofluid with anisotropic slip and magnetic field
  14. Numerical analysis of thermophoretic particle deposition on 3D Casson nanofluid: Artificial neural networks-based Levenberg–Marquardt algorithm
  15. Analyzing fuzzy fractional Degasperis–Procesi and Camassa–Holm equations with the Atangana–Baleanu operator
  16. Bayesian estimation of equipment reliability with normal-type life distribution based on multiple batch tests
  17. Chaotic control problem of BEC system based on Hartree–Fock mean field theory
  18. Optimized framework numerical solution for swirling hybrid nanofluid flow with silver/gold nanoparticles on a stretching cylinder with heat source/sink and reactive agents
  19. Stability analysis and numerical results for some schemes discretising 2D nonconstant coefficient advection–diffusion equations
  20. Convective flow of a magnetohydrodynamic second-grade fluid past a stretching surface with Cattaneo–Christov heat and mass flux model
  21. Analysis of the heat transfer enhancement in water-based micropolar hybrid nanofluid flow over a vertical flat surface
  22. Microscopic seepage simulation of gas and water in shale pores and slits based on VOF
  23. Model of conversion of flow from confined to unconfined aquifers with stochastic approach
  24. Study of fractional variable-order lymphatic filariasis infection model
  25. Soliton, quasi-soliton, and their interaction solutions of a nonlinear (2 + 1)-dimensional ZK–mZK–BBM equation for gravity waves
  26. Application of conserved quantities using the formal Lagrangian of a nonlinear integro partial differential equation through optimal system of one-dimensional subalgebras in physics and engineering
  27. Nonlinear fractional-order differential equations: New closed-form traveling-wave solutions
  28. Sixth-kind Chebyshev polynomials technique to numerically treat the dissipative viscoelastic fluid flow in the rheology of Cattaneo–Christov model
  29. Some transforms, Riemann–Liouville fractional operators, and applications of newly extended M–L (p, s, k) function
  30. Magnetohydrodynamic water-based hybrid nanofluid flow comprising diamond and copper nanoparticles on a stretching sheet with slips constraints
  31. Super-resolution reconstruction method of the optical synthetic aperture image using generative adversarial network
  32. A two-stage framework for predicting the remaining useful life of bearings
  33. Influence of variable fluid properties on mixed convective Darcy–Forchheimer flow relation over a surface with Soret and Dufour spectacle
  34. Inclined surface mixed convection flow of viscous fluid with porous medium and Soret effects
  35. Exact solutions to vorticity of the fractional nonuniform Poiseuille flows
  36. In silico modified UV spectrophotometric approaches to resolve overlapped spectra for quality control of rosuvastatin and teneligliptin formulation
  37. Numerical simulations for fractional Hirota–Satsuma coupled Korteweg–de Vries systems
  38. Substituent effect on the electronic and optical properties of newly designed pyrrole derivatives using density functional theory
  39. A comparative analysis of shielding effectiveness in glass and concrete containers
  40. Numerical analysis of the MHD Williamson nanofluid flow over a nonlinear stretching sheet through a Darcy porous medium: Modeling and simulation
  41. Analytical and numerical investigation for viscoelastic fluid with heat transfer analysis during rollover-web coating phenomena
  42. Influence of variable viscosity on existing sheet thickness in the calendering of non-isothermal viscoelastic materials
  43. Analysis of nonlinear fractional-order Fisher equation using two reliable techniques
  44. Comparison of plan quality and robustness using VMAT and IMRT for breast cancer
  45. Radiative nanofluid flow over a slender stretching Riga plate under the impact of exponential heat source/sink
  46. Numerical investigation of acoustic streaming vortices in cylindrical tube arrays
  47. Numerical study of blood-based MHD tangent hyperbolic hybrid nanofluid flow over a permeable stretching sheet with variable thermal conductivity and cross-diffusion
  48. Fractional view analytical analysis of generalized regularized long wave equation
  49. Dynamic simulation of non-Newtonian boundary layer flow: An enhanced exponential time integrator approach with spatially and temporally variable heat sources
  50. Inclined magnetized infinite shear rate viscosity of non-Newtonian tetra hybrid nanofluid in stenosed artery with non-uniform heat sink/source
  51. Estimation of monotone α-quantile of past lifetime function with application
  52. Numerical simulation for the slip impacts on the radiative nanofluid flow over a stretched surface with nonuniform heat generation and viscous dissipation
  53. Study of fractional telegraph equation via Shehu homotopy perturbation method
  54. An investigation into the impact of thermal radiation and chemical reactions on the flow through porous media of a Casson hybrid nanofluid including unstable mixed convection with stretched sheet in the presence of thermophoresis and Brownian motion
  55. Establishing breather and N-soliton solutions for conformable Klein–Gordon equation
  56. An electro-optic half subtractor from a silicon-based hybrid surface plasmon polariton waveguide
  57. CFD analysis of particle shape and Reynolds number on heat transfer characteristics of nanofluid in heated tube
  58. Abundant exact traveling wave solutions and modulation instability analysis to the generalized Hirota–Satsuma–Ito equation
  59. A short report on a probability-based interpretation of quantum mechanics
  60. Study on cavitation and pulsation characteristics of a novel rotor-radial groove hydrodynamic cavitation reactor
  61. Optimizing heat transport in a permeable cavity with an isothermal solid block: Influence of nanoparticles volume fraction and wall velocity ratio
  62. Linear instability of the vertical throughflow in a porous layer saturated by a power-law fluid with variable gravity effect
  63. Thermal analysis of generalized Cattaneo–Christov theories in Burgers nanofluid in the presence of thermo-diffusion effects and variable thermal conductivity
  64. A new benchmark for camouflaged object detection: RGB-D camouflaged object detection dataset
  65. Effect of electron temperature and concentration on production of hydroxyl radical and nitric oxide in atmospheric pressure low-temperature helium plasma jet: Swarm analysis and global model investigation
  66. Double diffusion convection of Maxwell–Cattaneo fluids in a vertical slot
  67. Thermal analysis of extended surfaces using deep neural networks
  68. Steady-state thermodynamic process in multilayered heterogeneous cylinder
  69. Multiresponse optimisation and process capability analysis of chemical vapour jet machining for the acrylonitrile butadiene styrene polymer: Unveiling the morphology
  70. Modeling monkeypox virus transmission: Stability analysis and comparison of analytical techniques
  71. Fourier spectral method for the fractional-in-space coupled Whitham–Broer–Kaup equations on unbounded domain
  72. The chaotic behavior and traveling wave solutions of the conformable extended Korteweg–de-Vries model
  73. Research on optimization of combustor liner structure based on arc-shaped slot hole
  74. Construction of M-shaped solitons for a modified regularized long-wave equation via Hirota's bilinear method
  75. Effectiveness of microwave ablation using two simultaneous antennas for liver malignancy treatment
  76. Discussion on optical solitons, sensitivity and qualitative analysis to a fractional model of ion sound and Langmuir waves with Atangana Baleanu derivatives
  77. Reliability of two-dimensional steady magnetized Jeffery fluid over shrinking sheet with chemical effect
  78. Generalized model of thermoelasticity associated with fractional time-derivative operators and its applications to non-simple elastic materials
  79. Migration of two rigid spheres translating within an infinite couple stress fluid under the impact of magnetic field
  80. A comparative investigation of neutron and gamma radiation interaction properties of zircaloy-2 and zircaloy-4 with consideration of mechanical properties
  81. New optical stochastic solutions for the Schrödinger equation with multiplicative Wiener process/random variable coefficients using two different methods
  82. Physical aspects of quantile residual lifetime sequence
  83. Synthesis, structure, IV characteristics, and optical properties of chromium oxide thin films for optoelectronic applications
  84. Smart mathematically filtered UV spectroscopic methods for quality assurance of rosuvastatin and valsartan from formulation
  85. A novel investigation into time-fractional multi-dimensional Navier–Stokes equations within Aboodh transform
  86. Homotopic dynamic solution of hydrodynamic nonlinear natural convection containing superhydrophobicity and isothermally heated parallel plate with hybrid nanoparticles
  87. A novel tetra hybrid bio-nanofluid model with stenosed artery
  88. Propagation of traveling wave solution of the strain wave equation in microcrystalline materials
  89. Innovative analysis to the time-fractional q-deformed tanh-Gordon equation via modified double Laplace transform method
  90. A new investigation of the extended Sakovich equation for abundant soliton solution in industrial engineering via two efficient techniques
  91. New soliton solutions of the conformable time fractional Drinfel'd–Sokolov–Wilson equation based on the complete discriminant system method
  92. Irradiation of hydrophilic acrylic intraocular lenses by a 365 nm UV lamp
  93. Inflation and the principle of equivalence
  94. The use of a supercontinuum light source for the characterization of passive fiber optic components
  95. Optical solitons to the fractional Kundu–Mukherjee–Naskar equation with time-dependent coefficients
  96. A promising photocathode for green hydrogen generation from sanitation water without external sacrificing agent: silver-silver oxide/poly(1H-pyrrole) dendritic nanocomposite seeded on poly-1H pyrrole film
  97. Photon balance in the fiber laser model
  98. Propagation of optical spatial solitons in nematic liquid crystals with quadruple power law of nonlinearity appears in fluid mechanics
  99. Theoretical investigation and sensitivity analysis of non-Newtonian fluid during roll coating process by response surface methodology
  100. Utilizing slip conditions on transport phenomena of heat energy with dust and tiny nanoparticles over a wedge
  101. Bismuthyl chloride/poly(m-toluidine) nanocomposite seeded on poly-1H pyrrole: Photocathode for green hydrogen generation
  102. Infrared thermography based fault diagnosis of diesel engines using convolutional neural network and image enhancement
  103. On some solitary wave solutions of the Estevez--Mansfield--Clarkson equation with conformable fractional derivatives in time
  104. Impact of permeability and fluid parameters in couple stress media on rotating eccentric spheres
  105. Review Article
  106. Transformer-based intelligent fault diagnosis methods of mechanical equipment: A survey
  107. Special Issue on Predicting pattern alterations in nature - Part II
  108. A comparative study of Bagley–Torvik equation under nonsingular kernel derivatives using Weeks method
  109. On the existence and numerical simulation of Cholera epidemic model
  110. Numerical solutions of generalized Atangana–Baleanu time-fractional FitzHugh–Nagumo equation using cubic B-spline functions
  111. Dynamic properties of the multimalware attacks in wireless sensor networks: Fractional derivative analysis of wireless sensor networks
  112. Prediction of COVID-19 spread with models in different patterns: A case study of Russia
  113. Study of chronic myeloid leukemia with T-cell under fractal-fractional order model
  114. Accumulation process in the environment for a generalized mass transport system
  115. Analysis of a generalized proportional fractional stochastic differential equation incorporating Carathéodory's approximation and applications
  116. Special Issue on Nanomaterial utilization and structural optimization - Part II
  117. Numerical study on flow and heat transfer performance of a spiral-wound heat exchanger for natural gas
  118. Study of ultrasonic influence on heat transfer and resistance performance of round tube with twisted belt
  119. Numerical study on bionic airfoil fins used in printed circuit plate heat exchanger
  120. Improving heat transfer efficiency via optimization and sensitivity assessment in hybrid nanofluid flow with variable magnetism using the Yamada–Ota model
  121. Special Issue on Nanofluids: Synthesis, Characterization, and Applications
  122. Exact solutions of a class of generalized nanofluidic models
  123. Stability enhancement of Al2O3, ZnO, and TiO2 binary nanofluids for heat transfer applications
  124. Thermal transport energy performance on tangent hyperbolic hybrid nanofluids and their implementation in concentrated solar aircraft wings
  125. Studying nonlinear vibration analysis of nanoelectro-mechanical resonators via analytical computational method
  126. Numerical analysis of non-linear radiative Casson fluids containing CNTs having length and radius over permeable moving plate
  127. Two-phase numerical simulation of thermal and solutal transport exploration of a non-Newtonian nanomaterial flow past a stretching surface with chemical reaction
  128. Natural convection and flow patterns of Cu–water nanofluids in hexagonal cavity: A novel thermal case study
  129. Solitonic solutions and study of nonlinear wave dynamics in a Murnaghan hyperelastic circular pipe
  130. Comparative study of couple stress fluid flow using OHAM and NIM
  131. Utilization of OHAM to investigate entropy generation with a temperature-dependent thermal conductivity model in hybrid nanofluid using the radiation phenomenon
  132. Slip effects on magnetized radiatively hybridized ferrofluid flow with acute magnetic force over shrinking/stretching surface
  133. Significance of 3D rectangular closed domain filled with charged particles and nanoparticles engaging finite element methodology
  134. Robustness and dynamical features of fractional difference spacecraft model with Mittag–Leffler stability
  135. Characterizing magnetohydrodynamic effects on developed nanofluid flow in an obstructed vertical duct under constant pressure gradient
  136. Study on dynamic and static tensile and puncture-resistant mechanical properties of impregnated STF multi-dimensional structure Kevlar fiber reinforced composites
  137. Thermosolutal Marangoni convective flow of MHD tangent hyperbolic hybrid nanofluids with elastic deformation and heat source
  138. Investigation of convective heat transport in a Carreau hybrid nanofluid between two stretchable rotatory disks
  139. Single-channel cooling system design by using perforated porous insert and modeling with POD for double conductive panel
  140. Special Issue on Fundamental Physics from Atoms to Cosmos - Part I
  141. Pulsed excitation of a quantum oscillator: A model accounting for damping
  142. Review of recent analytical advances in the spectroscopy of hydrogenic lines in plasmas
  143. Heavy mesons mass spectroscopy under a spin-dependent Cornell potential within the framework of the spinless Salpeter equation
  144. Coherent manipulation of bright and dark solitons of reflection and transmission pulses through sodium atomic medium
  145. Effect of the gravitational field strength on the rate of chemical reactions
  146. The kinetic relativity theory – hiding in plain sight
  147. Special Issue on Advanced Energy Materials - Part III
  148. Eco-friendly graphitic carbon nitride–poly(1H pyrrole) nanocomposite: A photocathode for green hydrogen production, paving the way for commercial applications
Downloaded on 10.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/phys-2023-0182/html
Scroll to top button