The nanofluid has a large number of nanoparticles and nanometer-sized molecules. Tiny particles are disseminated (added) in fluids to increase their heat transfer capabilities, even though the structure and content of nano-molecules depend on carbides, metals, and carbon nanotubes. The combinations that make up the nanofluid are typically combined with particles of average nanoscale size. Nanotubes, nanofibers, nanowires, nanoparticles, nanosheets, nanorods, and other materials are routinely made with nanofluids. The most effective and practical methods have been developed and put into practice for the modeling of nanofluid flow models, while the addition of specific materials of a particular kind has increased the thermal conductivity of various fluids. Several engineering procedures use nanoparticles in fluids to improve heat transfer. Equipment like heat exchangers in heavy machinery, automobiles, and industries is highly dependent on effective energy transfer. With these applications in mind, Choi and Eastman (Choi and Eastman, 1995) developed nanofluids to enhance the heat transfer capabilities of ordinary fluids. According to Eastman et al. (Eastman et al., 2001), adding copper nanoparticles to ethylene glycol at a concentration of 0.3 volumes results in a 40% increase in thermal conductivity. Khan and Pop (Khan and Pop, 2010) have numerically addressed the issue of laminar fluid flow that develops when a flat surface is stretched in a nanofluid. It was discovered that the reduced heat transfer is a diminishing function of every single dimensionless number. The numerical solution of the boundary layer flow caused by a linearly extending sheet in a nanofluid was inspected by Makinde and Aziz (Makinde and Aziz, 2011). They observed that the convective heating, thermophoresis, and Brownian motion all get more intense as the local temperature rises, which causes the thermal boundary layer to thicken. The thermal conductivity performance of carbon nanotubes in fluid flow over a stretching sheet was discovered by Haq et al. (Haq et al., 2015). Sheikholeslami et al. (Sheikholeslami et al., 2016) discovered the impacts of Lorentz forces on free convective flow in the presence of nanofluid with thermal radiation. A two-dimensional time-independent flow conveying nanofluid toward a thin needle was deliberated by Soid et al. (Soid et al., 2017) where the existence of multiple solutions was reported. Bakar et al. (Bakar et al., 2018) looked into the stability analysis of mass suction impacts through a shrinking/stretching cylinder considering the nanoparticles. Kamal et al. (Kamal et al., 2019) investigated the flow of nanomaterial through a stretchable/shrinkable sheet with a chemical reaction effect. A two-dimensional magneto mixed convection flow induced by a shrinking/stretching plate in a nanofluid was investigated by Jumana et al. (Jumana et al., 2020). Johan and Mansur (Johan and Mansur, 2021) examined the features of dusty nanomaterial flow and thermal transport analysis past a stretchable sheet with a slip boundary. They used three kinds of particles, namely, alumina, copper, and titania. Shahzad et al. (Shahzad et al., 2022) scrutinized the influence of the slip effect on the flow and heat transport incorporated by copper nanofluid with different shape factors through a heated stretched sheet and found that the platelet shape factor particles have a greater rate of heat transfer compared to other shape factors. Shamshuddin et al. (Shamshuddin et al., 2024a) examined the ohmic and ferromagnetic effects on the flow of nanofluid through a porous rotating disk. They discovered that in normal fluid and nanofluid, the increased suction velocity suppresses the fluid velocity produced by the surface’s porosity. Recently, Shamshuddin et al. (Shamshuddin et al., 2024b) explored the impact of electromagnetic on the radiative flow of tri-hybrid nanofluid through a bidirectional stretchy surface in penetrable media. The factors of the surface frictional force of the ternary hybrid nanofluids are found to decrease with increasing the Lorentz force and number of surface pores.

The problem of sheets stretching or shrinking in a viscous fluid has caught a lot of concentration as it has numerous applications in physics, engineering, and other scientific disciplines. It frequently occurs in practical issues that have attracted a lot of research attention because of their wide range of significance in fields like the production of glass fibre, glass blowing, metal extrusion, transportation, microfluidics, paper production, hot rolling, space, and acoustics (see Fisher (Fisher, 1976)). The BLF over an ongoing solid kind surface flowing at uniform motion was initially studied by Sakiadis (Sakiadis, 1961) in light of these applications. Numerous authors (Chen and Char, 1988; Ishak et al., 2009; Mi, 2015; Guo and Fu, 2019; Wang, 2019; Zhao et al., 2019; Zi and Wang, 2019) have thought about different elements of this problem and found similarity solutions since the groundbreaking research conducted by Crane (Crane, 1970), who provided an exact solution for the 2D steady flow caused by a stretchable surface in a quiescent fluid. Instead of focusing on the scenario of a stretched sheet, researchers instead looked at the scenario of a shrinking sheet. According to Goldstein (Goldstein, 1965), this new kind of flow of shrinking sheet is fundamentally a backward flow. The steady flow across a shrinkable sheet was studied by Miklavčič and Wang (Miklavčič and Wang, 2006). They discovered that mass suction is necessary to continue the flow across a shrinkable sheet. Waini et al. (Waini et al., 2019) examined the dependency of time-varying flow along with thermal transport across a shrinking/stretching sheet incorporated with hybrid nanofluids and provided multiple solutions.

The examination of cross-flow began following the early studies by Prandtl (Prandtl, 1946a) and Blasius (Blasius, 1908) that included the flow over a smooth surface caused by thin viscosity. Prandtl (Prandtl, 1946b) is believed to be the initial researcher to publish the findings for uniform pressure gradients flowing through a finite yawed cylinder. Weidman et al. (Weidman, 2017) examined the boundary layer via cross-flow generated by transverse plate motions. The work of Weidman was recently expanded by (Roşca et al., 2021) by taking into account rotational stagnation-point flow that transports hybrid nanofluids along a permeable shrinking or stretching surface. It was discovered that both stretching and shrinking surfaces can have multiple solutions to the fundamental similarity equations.

The phenomena of an irregular heat sink or source have applications in both engineering and medicine, involving the recovery of crude oil, the construction of thrust bearings, and the cooling of metallic sheets, etc. In the presence of an irregular heat sink/source, Tawade et al. (Tawade et al., 2016) addressed the motion of the MHD unsteady thin film and heat transfer past a stretchable sheet. It was determined that irregular heat parameters are crucial to the effectiveness of heat transfer. Thumma et al. (Thumma et al., 2017) revealed that the stretching of a sheet caused the MHD convective motion of nanofluid to have a changeable heat sink or source. To obtain the solution, a well-known Keller-Box numerical approach is utilized. Kumar et al. (Kumar et al., 2020) looked at the movement of hybrid ferrofluids film and heat transfer in the inclusion of radiation and erratic heat source/sink (EHSE/EHSK). The rate of heat transfer is believed to be greater in hybrid ferrofluids than in ferrofluids. In addition, the velocity of the fluid and temperature tend to decline as film thickness parameters increase. Areekara et al. (Areekara et al., 2021) investigated the impact of an irregular heat source/sink on the fluid flow of nanofluid past a nonlinear stretching sheet. They observed that positive correlations exist between the radiative heat flux and the rate of heat transfer. A negative sensitivity to the rate of heat transfer is shown by the exponential heat source. Akram et al. (Akram et al., 2022) discussed the concepts of non-linear stretching and EHS/SHS to describe the heat transfer through the stretchable cylinder. It is discovered that the temperature distribution in the fluid region is being disrupted by the non-linear stretching rate and the source of heat.

The scrutiny of second law analysis or, EG (entropy generation) in fluid flow and heat transport is a prominent area of study. Energy losses resulting from chemical processes, diffusion, solid-surface friction, and the viscosity of fluids all contribute to the production of entropy in thermodynamic systems. As a result, the formation of entropy generation (EG) has drawn a lot of attention to applications including heat exchangers, turbo machinery, and electronics cooling. Aiboud and Sauoli (Aiboud and Saouli, 2010) scrutinized the entropy optimized in viscoelastic flow through a flexible surface subjected to the magnetic field. The effect of slip across a heated vertical surface in entropy-optimized flow was calculated by But et al. (Butt et al., 2012). Slips have been found to allow for the control and adjustment of entropy formation in thermal systems. Abolbashari et al. (Abolbashari et al., 2014) employed HAM to examine, EG in magneto nanofluid flow near an unsteady stretched surface with H₂O-base liquid and several nanoparticle types. Tlau and Ontela (Tlau and Ontela, 2019) examined the role of magnetohydrodynamics on nanomaterial entropy-optimized flow from an inclined channel with a heat source/sink embedded in a porous media. Entropy optimization of nanomaterials flow across two stretchable rotating disks with effects of bioconvection was examined by Khan et al. (Khan et al., 2020). Mondal et al. (Mondal et al., 2021) used trapezoidal lid-driven enclosures with, EG to study the Lorentz forces on the constant buoyant flow of Al₂O₃ nanoparticles. According to the calculations, the average Nusselt and Sherwood numbers, and aspect ratio all decrease with increasing the percentage of nanoparticle volume.

The literature that is currently accessible indicates that no exploration has been performed on the entropy generation of the radiative cross-flow stimulated with nanofluid from a porous expanding/contracting sheet with an asymmetrical heat source/sink. The present investigation fills a research gap by demonstrating cross-flow and heat transfer towards a stagnation point of nanofluid via an expanding/contracting sheet with asymmetrical heat sink/source and thermal radiation. The proposed model is originally described via a highly nonlinear system of PDEs. The PDEs are rehabilitated into ODEs by using the proper similarity variables and then solved by employing a bvp4c solver. Multiple solutions are provided for certain values of the parameters such as mass suction and shrinking sheet. The study described in this paper was driven by the following research questions.

• What impact does have on the skin friction and Nusselt number by raising the suction parameter against the stretching/shrinking parameter?

The nanofluids’ stagnation points radiative cross flow and thermal system characteristics past a contracting/expanding surface with the mutual influence of mass suction and EHSE/EHSK are taken in this study. As portrayed in Figure 1, x a , y a are Cartesian coordinates measured along the horizontal and vertical surface of the sheet, respectively, with the flow occurring in the domain y a ≥ 0 . Nanofluid is a mixture of regular fluid (water) and single titanium dioxide (TiO₂) nanoparticles. The investigational features of the (water/TiO₂) nanomaterials are taken to be uniform. In addition, the horizontal surface of the sheet is assumed to have a variable velocity of u a = ε b 2 / 3 ν f 1 / 3 x a 1 / 3 γ b = u w a x a γ b , where γ b refers to the expanding/contracting factor with γ b < 0 , γ b > 0 and γ b = 0 signify the particular cases of the shrinking, stretching, and stationary/static sheet, respectively. In the meantime, ε b indicates the positive constant and υ f is the kinematic viscosity of the regular (water) fluid. The ambient or far-field (nanofluid) is also supposed to have a linear velocity of u e x y a = ε b y a , see Weidman et al. (Weidman, 2017). Moreover, the mass suction/injection or transpiration velocity at the surface of the sheet is v w a x a with v w a x a < 0 refers to the case of injection and v w a x a > 0 refers to the case of suction while v w a x a = 0 indicating the impermeable surface of the sheet. It is also supposed that the constant temperature of the sheet is T w a , while T ∞ represents the free stream temperature (inviscid fluid). With the help of these aforesaid assumptions, the governing equations in the Cartesian form are written as (Weidman, 2017; Roşca et al., 2021): ∂ u a ∂ x a + ∂ v a ∂ y a = 0 , (1) ρ n f u a ∂ u a ∂ x a + v a ∂ u a ∂ y a = μ n f ∂ 2 u a ∂ y a 2 , (2) ρ c p n f u a ∂ T a ∂ x a + v a ∂ T a ∂ y a = k n f ∂ 2 T a ∂ y a 2 − ∂ Q r a d ∂ y a + k n f u w a x a x a v n f A b * T w a − T ∞ e − η + B b * T a − T ∞ , (3) with boundary conditions (BCs) are: u a = γ b u w a x a , v a = v w a x a , T a = T w a at y a = 0 , ∂ u a ∂ y a → ∂ u e x ∂ y a = ε b , T a → T ∞ as y a → ∞ . (4)

In Eqs 1–4, u a and v a are the nanofluid velocities in the corresponding x a , and y a directions, T a the temperature of the nanofluid, A b * the exponentially decaying space coefficients, and B b * the temperature-dependent heat source/sink. Therefore, the heat source or absorption phenomenon is produced due to the positive value of A b * and B b * while the phenomenon of heat generation or sink is found by the negative value of both A b * and B b * .

The radiation heat flux Q r a d is expressed by the Rosseland approach as: Q r a d = − 4 σ a 3 k a ∂ T a 4 ∂ y a . (5)

Here, the Stefan Boltzmann constant and the mean absorption coefficient are denoted by k a and σ a , respectively. Moreover, the term T a 4 is simplified further by using the Taylor series at T ∞ and overlooking the power of higher-order yields: T a 4 ≅ 4 T a T ∞ 3 − 3 T ∞ 4 . (6)

In addition, executing Eqs 5, 6 into Eq 3 yields the final form: ρ c p n f u a ∂ T a ∂ x a + v a ∂ T a ∂ y a = k f k n f k f + 4 3 N r ∂ 2 T a ∂ y a 2 + k n f u w a x a x a v n f A b * T w a − T ∞ e − η + B b * T a − T ∞ , (7) where N r = 4 σ a T ∞ 3 k a k f signifies the thermal radiation parameter.

Furthermore, k n f indicates the thermal conductivity of the essential posited nanofluid (NFD), ρ c p n f indicates the heat capacitance of the NFD, ρ n f indicates the density of the NFD, and μ n f indicates the absolute viscosity of the NFD. The correlation of these NFDs is written as follows: k n f = k T i O 2 + 2 k f − 2 φ T i O 2 k f − k T i O 2 k T i O 2 + 2 k f + φ T i O 2 k f − k T i O 2 , ρ c p n f ρ c p f = φ T i O 2 ρ c p T i O 2 ρ c p f + 1 − φ T i O 2 , ρ n f ρ f = φ T i O 2 ρ T i O 2 ρ f + 1 − φ T i O 2 , μ n f ′ μ f = 1 − φ T i O 2 − 2.5 . (8)

Here, k f , ρ f , and μ f refer to the thermal conductivity, the density, and the absolute viscosity of the base (water) fluid while the heat capacity at constant pressure is represented via c p . Therefore, φ T i O 2 symbolizes the titania nanoparticles volume fraction and the special case φ T i O 2 = 0 reduces the Eqs 8 to a normal or a regular fluid (water). In addition, Table 1 displays the physical data of the titania (TiO₂) nanoparticles and the regular fluid (water).

For the considered model, the similarity transformations that can be expressed to further simplify the procedure for mathematical analysis are as follows: η = ε b / υ f 1 / 3 y a x a 1 / 3 , u a = ε b 2 / 3 υ f 1 / 3 x a 1 / 3 G ′ η , H η = T a − T ∞ T w a − T ∞ , v a = − ε b 1 / 3 υ f 2 / 3 3 x a 1 / 3 2 G η − η G ′ η , (9) which provides the opportunity to describe flows extremely broadly, regardless of the system size. Also, the prime corresponds to the derivative with respect to η , H is the non-dimensional temperature distribution, G describes the non-dimensional quantities, and G ′ is the dimensionless velocity profile. However, the mass suction/injection velocity at the surface of the sheet is written as: v w a x a = − 2 3 ε b υ f 2 x a 1 3 f w a . (10)

In Eq. 10, f w a is the constant mass suction/blowing constraint with f w a = 0 , f w a < 0 and f w a > 0 describe the phenomena of impermeable, blowing, and suction, respectively.

With the help of similarity transformations (9), the continuity Eq 1 of the governing model is satisfied while the rest of Eqs 2, 7 equations change to the resulting known ordinary (similarity) differential equations (ODEs) as: μ n f / μ f ρ n f / ρ f G ‴ + 2 3 G G ″ − 1 3 G ′ 2 = 0 , (11) k n f k f + 4 3 N r H ″ + 2 3 Pr ρ c p n f ρ c p f G H ′ + k n f k f ρ n f ρ f μ n f μ f A b * e − η + B b * H = 0 , (12) along with border conditions: G 0 = f w a , G ′ 0 = γ b , H   0 = 1  at  η = 0 , G ″ η → 1 , H η → 0  as  η → ∞ . (13)

In addition, Eq 11 for the case of φ T i O 2 = 0 is the same as equation (6.1, see Weidman et al. (Weidman, 2017)) when α = 1 , but Eq 12 with some special effects has been not taken in the same reference paper. Moreover, the dimensionless model comprised the following distinct parameters, the suction/injection f w a , the expanding/contracting γ b , and the Prandtl number Pr ⁡ = υ f / α f .

Entropy generation (E.G.,), also known as the second law analysis, is a necessary instrument for measuring energy loss and depreciation in the effectiveness of engineering and industrial systems, such as rate and transport operations. As a result, the systems expend less energy, making, E.G., analysis and comprehension crucial. Taking into account the scenario of, E.G., for the viscous Newtonian liquid with the inclusion of nanoparticles (Abolbashari et al., 2014; Tlau and Ontela, 2019). E G = 1 T ∞ 2 k n f + 16 3 σ a T ∞ 3 k a ∂ T a ∂ y a 2 + μ n f T ∞ ∂ u a ∂ y a 2 . (16)

Two fundamental elements are principally responsible for the, E.G., in the contemplated cross-flow of viscous Newtonian nanofluid. The first term in the statement, which is on the right-hand side, denotes the ensuing local heat transfer, and the final term, the consequent fluid friction or viscous dissipation. The second law analysis is defined as follows in the dimensionless form: N G * = x a 2 T ∞ 2 k f T w a − T ∞ 2 E G . (17)

The following formulas are obtained by incorporating the similarity transformations from Eq 9 into Eq 16. Hence, N G * = k n f k f + 4 3 N r Re x a H ′ 2 + μ n f μ f Re x a B r b Ω b G ′ ′ 2 , (18) where Ω b the dimensionless temperature difference and B r b the Brinkman number. Consequently, they are represented mathematically as: Ω b = T w a − T ∞ T ∞ , B r b = μ f u w a 2 k f T w a − T ∞ . (19)

This section demonstrates the analysis of the assumed cross-flow and the suspension of nanofluid for heat transfer. The requisite model equations are expressed as highly nonlinear ODEs (11) and (12) along with BCs (13) using similarity variables (9). A built-in function named bvp4c included in the MATLAB software is used to work out these equations numerically. It ought to be noted that the scheme of the finite difference is the foundation for the bvp4c package, which is further highlighted by the 3-stage Lobatto IIIA procedure. To instigate the bvp4c method, the transmuted ODEs are modified into a first-order system by launching new-fangled variables. By establishing this process, let G = A 1 , G ′ = A 2 , G ″ = A 3 , H = A 4 , H ′ = A 5 . (20)

Substituting Eq. 20 into Eqs 11, 12 along with BCs Eq. 13 we obtain the set of first-order ODEs as follows: d d η A 1 A 2 A 3 A 4 A 5 = A 2 A 3 ρ n f / ρ f μ n f / μ f 1 3 A 2 2 − 2 3 A 1 A 3 A 5 1 k n f k f + 4 3 N r − 2 3 Pr ρ c p n f ρ c p f A 4 A 2 − k n f / k f ρ n f / ρ f μ n f / μ f A b * e − η + B b * A 4 , (21) with BCs A 1 0 = f w a , A 2 0 = γ b , A 4 0 = 1 , A 3 ∞ = 1 , A 4 ∞ = 0 . (22)

The code desired initial estimations at the posited mesh point to solve Eq 21 and the corresponding conditions Eq. 22. The polynomial used in the collective type yields a continuous result. A fourth order accuracy set that is equally distributed over the spatial intervals where the function is integrated provides the result. The limitation at a distance η → ∞ is replaced by the value η = η ∞ = 8 in many successful boundary layer theory applications, and the relative tolerance error is pre-defined as 10⁻⁶. Additionally, the residual of the smooth output serves as the basis for both the mesh selection and error control. The initial mesh comprises four identical discredited points on the range 0 , η ∞ = 8 , and the mesh selection is then automatically changed using the bvp4c package. The problem now has two possible solutions, which means that the bvp4c package needed two alternative guesses for the unstable and stable solutions. The early initial presumes for the upper solution is moderately straightforward, whilst selecting a guess for the second solution is fairly challenging. Merkin et al. (Merkin, 1986) and Weidman et al. (Weidman et al., 2006) claim that the upper solution is physically stable and reliable but the lower solution is unstable and not physically dependable since the outcome only exists for a specific range of shrinking sheets.

The current portion of the work describes the binary (stable and unstable) solutions for a certain region or area of the dimensionless mass suction parameter as well as the contracting parameter γ b due to the variations in one physical parameter at the time of computations while the other factors are taken to be fixed. To ease analysis, the values of the basic physical influential parameters are φ T i O 2 = 0.025 , γ b = − 2.0 , f w a = 1.5 , A b * = 0.1 , N r = 0.05 , and B b * = 0.1 . Table 1 reveals the experimental data of titania (TiO₂) nanoparticles. Meanwhile, the comparison of the upshots for the unusual case is demonstrated in Table 2. Besides, the outcomes of the friction factor, heat transfer, and entropy generation are captured in various graphs (see Figures 2–8) of the nanofluid for the unstable and stable branches owed to the influence of the several comprised factors while their quantitative outputs are shown in Tables 3, 4. However, the branches of stable solutions (SBES), as well as solutions of unstable (USBES) are categorized by the black solid and black dash lines, respectively. The position in the graph or picture where both (SBES and USBES) curves meet at a single point is called the bifurcation or critical point. In this study, the SBES and USBES are invented only for the case of shrinking parameters.

The numerical data of the gradients (skin friction and heat transfer rate) with the impression of the several distinguished factors corresponding to water-based Titania nanofluid are illustrated in Tables 3, 4 for the SBSE and USBES, respectively. Upshots divulge that the friction factor upsurges for the SBSE owed to the superior values of φ T i O 2 and f w a while the branch of USBES behaves distinctly with variations in the mass suction parameter f w a but similarly with higher impacts of φ T i O 2 . Notably, the shear stress of the nanofluid is highest and lowest for the SBES and USBES with mass suction f w a . In contrast, the heat transfer escalates for both (SBES and USBES) results with superior consequences of φ T i O 2 . Therefore, owing to the rise in the heat source parameter A b * , B b * > 0 , the heat transfer shrinkages in the SBES as well as the USBES are endlessly enriching due to the higher role of the heat sink parameter A b * , B b * < 0 . In addition, the lowest and highest heat transfer approximations are perceived for the branch of SBES and USBES with the following selected values A b * , B b * = 0.9 , and φ T i O 2 = 0.035 .

Figures 2, 3 show the impact of f w a on the shear stress and heat transfer corresponding to (TiO₂/water) nanofluid for the SBES as well as the USBES, respectively. In this study, the dual (SBES and USBES) outcomes are possible to occur for a certain domain of the specific set of physical parameters. Therefore, the non-unique outcomes in either pictures or graphs exist for a posited shrinkable sheet. Moreover, it is clear from the above graphs that the position where both solution curves meet is at a point called the critical point. Mathematically, this point is expressed by the symbol γ b C , where the solutions are unique γ b = γ b C . The non-unique solutions (dual) and no solutions are possible to exist for the range γ b C < γ b < ∞ and − ∞ < γ b < γ b C , respectively. Besides, the outcomes refer to the shear stress and heat transfer escalating for the SBES due to the larger impact of f w a while they are declined for the USBES. Physically, the motion of the nanofluid stops due to the inspiration of f w a shifting the particles of the liquid moving toward the surface of the sheet and sticking with it. Hence, the friction and motion/velocity of the nanofluid hold the inverse relations, as a result, the shear stress is enhanced. Furthermore, the next eight distinct critical values −2.1903, −2.2941, −2.4016, −2.5125, −2.6274, −2.7463, −2.8686, and −2.9951 are obtained for the respective change value of f w a . With the rise of f w a causes an increment in the absolute value of γ b C . This behavior corresponds to the superior inclusion of f w a decelerating the boundary layer separation.

The impressions of γ b on C f Re x a 1 / 2 and N u x a Re x a − 1 / 2 versus f w a of (TiO₂/water) nanofluid for both (SBES and USBES) results are typified in Figures 4, 5, respectively. The dual (SBES and USBES) results are shown in both graphs for the case of f w a due to the variations in the shrinking parameter. In both graphs, it is seen that the SBES and USBES curves congregate at a point called the critical point which is denoted by f w a C . Meanwhile, the outcomes are unique for the case when f w a = f w a C , but the phenomena f w a C < f w a < ∞ , and − ∞ < f w a < f w a C indicate the non-unique and no solutions, respectively. Besides, the shear stress decays and rises for the SBES due to the higher values of γ b while it shrinks for the USBES. Alternatively, with the increase of γ b , the thermal transport phenomenon uplifts for the SBES and declines for the USBES. More significantly, it is understood from the diagrams that bifurcation values like 2.1615, 2.0760, 1.9888, 1.8986, 1.8057, 1.7096, 1.6098, and 1.5066 are found due to the several values of γ b . Also, it is noted here that the magnitude of the critical values f w a C is weakened once the shrinking parameter is boosted. This further specifies that the growth in the impacts of γ b hastens the boundary layer (BL) separations.

With the assistance or support of entropy generation or second law analysis (SLA), the thermodynamic system performance of water-based Titania nanoparticles can be improved. Figures 6–8 describe the consequence of parameters B r b , φ T i O 2 , and Ω b on SLA (second law analysis) corresponding to the posited nanofluid for the SBES and USBES, respectively. As premeditated in Figures 6, 7, an improvement in both parameters of B r b and φ T i O 2 results in an improvement of SLA for the branch of stable as well as unstable solutions. The second law analysis is also more susceptible to changes in this parameter B r b value at the place close to the wall surface of the sheet. As increases η , susceptibility rapidly decreases. Therefore, the physical data and prior research findings are compatible with this occurrence. Conversely, Figure 8 is designed to inspect the parameter Ω b on second law analysis of the (TiO₂/water) nanofluid for both (SBES and USBES). In general, the lower temperature difference between the wall and its surroundings is generally caused by larger values of Ω b , which reduces second law analysis or, EG. The aforementioned discussion leads to the conclusion that by altering the related parameters, the system’s second law analysis value can be decreased to increase the solar radiation utilization system’s effectiveness.

The theoretical inspection on radiative cross flow and heat transfer incorporated water-based TiO₂ nanofluid through a permeable stretching/shrinking sheet with irregular heat source/sink have been explored. The entropy generation was used to analyze the heat transfer process after the development of a computational model. Combined impacts of pertaining governing parameters like suction, expanding/contracting parameter, irregular heat source/sink parameter, radiation parameter, and volume fraction of the nanoparticles on shear stress and the heat transfer have been analyzed. The important conclusions of our research can be summed up as follows.

• The results suggest that the shear stress enhances due to f w a and φ T i O 2 , whilst the heat transfer accelerates due to φ T i O 2 .

Further, this problem can also be studied by incorporating unsteady flow or mixed convection flow in addition to incorporating various physical elements such as the slip effect, chemical reaction viscous dissipation, etc. Also, this model can be extended to any non-Newtonian models of the problem along with several other impacts like shape factors of the nanoparticles, waste discharge concentration, and thermophoresis particle deposition effects.