Engineers and scientists have held a longstanding fascination with researching non-Newtonian fluids, primarily because of their behavior on stretched surfaces with heat and mass transfer. This is because of its numerous practicalities, such as paper manufacturing, wire drawing, paper manufacturing, petroleum production, and polymer sheet manufacturing. Because shear rate and shear stress are not linearly related in non-Newtonian fluids, one constitutive equation is insufficient for representing the relationship between shear rates and shear stresses. Owing to this, various frameworks have been presented in the literature, including the Maxwell model, Reiner–Rivlin model, generalized Burgers, Oldroyd-A model, Oldroyd-B model, Oldroyd-8 constants, Carreau model, and the Carreau–Yasuda model. In the current investigation, we took non-Newtonian tangent hyperbolic fluid under consideration. Effects of shear thinning are described by a four-constant pseudoplastic model of fluid, i.e., tangent hyperbolic fluid. When flow decreases, owing to increased shear stress, it is useful to evaluate the inter-particle friction in the fluid. Examples of such liquids include ketchup, paint, blood, nail polish, and whipped cream. [1] presented a numerical approach for the flow of rate-type fluid across the stretchable cylinder using MHD flow. [2] focused on thermal radiation, which is non-linear in nature, to probe the flow by a stretchable surface with Brownian diffusion and thermophoresis. The flow of the two-dimensional rate-type fluid toward a stretchable sheet with the applied magnetic effect was investigated [3]. [4] found solutions in the series form for a tangent hyperbolic fluid that is inclined uniformly in a tube. [5] analyzed the tangent hyperbolic fluid using a stretchable sheet while considering both the impacts of Soret and Dufour in MHD stagnation point flow. [6–13] provided more research articles on non-Newtonian fluids.

Recently, some energy-efficient mechanisms have been designed to minimize the quantity of wasted energy and boost energy utilization capabilities of liquid systems. The transport of heat and viscous dissipation results in energy waste, which aids the irreversibility of the liquid system. Electrochemical and chemical, evaporative cooling, natural convection, solar thermal, microchannels, air separators, chillers, gas turbines, fuel cells, and functionally graded materials are all examples of systems where entropy generation minimization (EGM) is beneficial. [14] recently analyzed the generation of entropy in nanofluids using copper and silver nanoparticles. [15] studied the increase in entropy in a burner model with various energy sources. The optimization of entropy in a nonlinear radiative flow while accounting for Joule effects and dissipation is investigated by [16].

The change in density over a region induces additive rapid movement in self-propelled microorganisms; as a result of this, there is macroscopic movement of the fluid, which is termed as bioconvection. Self-propelled microorganisms have a tendency to enhance the base fluid to create a bio-convective stream in a specific direction. The self-propelled microorganisms are categorized into different types, i.e., gyrotactic, negative gravitaxis, and chemotaxis or oxytactic microorganisms. Compared to self-propelled microorganisms, nanoparticles do not move on their own, and their movement is driven by the effect of Brownian motion and thermophoresis. If there is minimal nanoparticle concentration, then in the nanofluids, bioconvection is supposed to be feasible and it will not result in an essential increase in the thickness of the base fluid. Bioconvection with nanoparticles was originally discussed in [17, 18]. Afterward, combining nanoparticles with gyrotactic microorganisms was proposed by [19] using Buongiornos theory.[20] studied the flow of nanofluids with bioconvection in porous medium numerically. Bioconvective flow with nanoparticles through a symmetric channel was studied by [21], and a bio nanoengineering model was presented. In rapidly changing magnetic lines of force, the clot blood model was also analyzed by [22] using the Jeffrey fluid model with microorganisms and nanoparticles. A magnetized laminar flow of nanofluid along with phototactic microorganisms in a medium that is non-Darcy with pores was observed by [23]. The external magnetic effect and bioconvective flow were researched by [24] with convective boundary conditions and nanoparticles.

In the recent past, as a result of advancements in nanotechnology, the researchers’ focus has shifted to improving heat transport through the nanoparticle interaction. Because of their fascinating thermo-physical aspects, the nanoparticles, which include metallic particles of micro-size, are useful in biological, mechanical, chemical, and other engineering fields. The conductivity of heat energy and the viscosity of the fluid improve not just with the addition of nanoparticles but also with the structure of particle size, and slip features in nanoparticles play an important role as well. [25] developed nanoparticles to improve the product efficiency and fluid thermal conductivity. Heat transfer near a stagnation point in the Maxwell nanofluid flow over a porous rotating disk, highlighting the impact of shear stress and flow field alterations on heat transfer, has been investigated by [26]. Some interesting studies on nanofluids have been discussed in [27–30]. [31] analyzed the magnetohydrodynamic stagnation point flow of Williamson hybrid nanofluid over a porous sheet, considering chemical reactions, energy generation, and boundary conditions. [32] examined the unsteady 2D laminar flow of Williamson hybrid nanoliquid over a convectively heated stretching sheet, focusing on maximizing energy and mass transfer rates. [33] investigated entropy generation in the mixed convection time-dependent flow of cross-hybrid nanoliquid at a stagnation point using engine oil with C u O   a n d   T i O 2 nanoparticles.

Activation energy and chemical reactions have been used in mass transportation phenomena, which have been carried out in the history of food sciences, chemistry, chemical, and mechanical industries. In 1889, Svante Arrhenius engaged in the key observations on activation energy and defined it in terms of the least quantity of energy that is required to begin a chemical process. Geothermal sciences, petroleum engineering applications, oil reserves, heat reactions, and a lot of other fields are part of the activation energy phenomenon. Some major contributions to this subject are discussed in [34–36]. [37] numerically investigated the natural convection boundary-layer flow over a truncated cone with gyrotactic microorganisms in water-based nanofluid, revealing enhanced transport properties and boundary conditions’ effects. [38] explored the bioconvection phenomenon of a nanofluid with gyrotactic microorganisms around a horizontal circular cylinder, utilizing Buongiorno’s model and finite-difference algorithm for the comprehensive analysis of transport phenomena. [39] investigated mixed convection of nanofluid with gyrotactic microorganisms over a permeable vertical surface in a porous medium, which reveals suction/injection effects, with significant impacts on heat and mass transfer. [40] also examined the mixed bio-convective stagnation point flow of a power-law nanofluid over a stretchable surface with passively controlled boundary conditions. [41] explored non-Newtonian nanofluid flow induced by a rotating stretchable disk, analyzing heat and mass transfer with Stefan blowing and Cattaneo–Christov fluxes. He also examined the flow of a non-Newtonian nanofluid in a conical gap between rotating/stationary surfaces, considering heat and mass transfer and a uniform magnetic field [42]. Some of the similar results are discussed in [43–45]. The study of the fluid flow in porous medium is of great importance in fluid mechanics. [46] analyzed the motion of Jeffrey nanofluid in porous medium with motile microorganisms between two revolving stretching disks. The study of different types of fluids in between porous medium is given in [47–50]. Furthermore, heat and mass transfer analysis of nanofluid is numerically scrutinized under the influence of magnetic field subject to marangoni boundary layer [51, 52].

The major goal of this work is to look at the optimization of entropy in the tangent hyperbolic nanofluid flow on a stretchable sheet. The uniqueness of the work under discussion is that it investigates the optimization of entropy generation in a non-Newtonian fluid model with non-linear thermal radiation, non-linear mixed convection, viscous dissipation, bioconvection, and activation energy as additional effects. To characterize activation energy’s effects, we used a unique two-step chemical reaction model. Various flow variables are utilized to compute the rate of total entropy generation. We have used a method called the shooting technique to deal with the nonlinear ordinary system. To examine the behavior of the included variables on fluid properties, graphs and tables are created. Analysis on various methodologies is shown in [53–56].

In tangent hyperbolic nanomaterial flow, we analyzed Bejan numbers and entropy. In the current model, the impacts of non-linear heat radiation, non-linear mixed convection, activation energy, and Joule effect are studied. The stretchable sheet with the stretching rate being a generates the flow. T ∼ w , C ∼ w , and N ∼ w signify the temperature, concentration, and density of self-propelling microorganisms at the sheet, respectively, while T ∼ ∞ , C ∼ ∞ , and N ∼ ∞ represent the ambient temperature, concentration, and density of self-propelling microorganisms, respectively. B 0 magnetic strength is in the y-direction (see Figure 1). In light of these considerations, the governing equations for the given problem are Equations 1–6 [57]: ∂ u ∼ 1 ∂ x + ∂ u ∼ 2 ∂ y = 0 , (1) u ∼ 1 ∂ u ∼ 1 ∂ x + u ∼ 2 ∂ u ∼ 1 ∂ y = ∂ 2 u ∼ 1 ∂ y 2 ν 1 − s + ∂ u ∼ 1 ∂ y ∂ 2 u ∼ 1 ∂ y 2 2 ν s Γ − σ ρ f B 0 2 u ∼ 1 + g β 1 − C ∞ T ∼ − T ∼ ∞ − g ρ m − ρ f ρ f C ∼ − C ∼ ∞ − γ 1 g ρ m − ρ f ρ f N ∼ − N ∼ ∞ , (2) u ∼ 1 ∂ T ∼ ∂ x + u ∼ 2 ∂ T ∼ ∂ y = ∂ 2 T ∼ ∂ y 2 k ρ c p f + D B ∂ C ∼ ∂ y ∂ C ∼ ∂ y + D T T ∼ ∞ ∂ T ∼ ∂ y 2 ρ c p s ρ c p f + 3 T ∼ 2 ∂ T ∼ ∂ y 2 + T ∼ 3 ∂ 2 T ∼ ∂ y 2 × 1 ρ c p f 16 σ * 3 k * + σ ρ c p f B 0 2 u ∼ 1 2 + ∂ u ∼ 1 ∂ y 2 μ ρ c p f 1 − s + ∂ u ∼ 1 ∂ y ∂ u ∼ 1 ∂ y 2 × s μ Γ ρ c p f 2 , (3) u ∼ 1 ∂ C ∼ ∂ x + u ∼ 2 ∂ C ∼ ∂ y = ∂ 2 C ∼ ∂ y 2 D B + ∂ 2 T ∼ ∂ y 2 D T T ∼ ∞ − k r 2 T ∼ T ∼ ∞ n exp − E a κ T ∼ C ∼ w − C ∼ ∞ , (4) u ∼ 1 ∂ N ∼ ∂ x + u ∼ 2 ∂ N ∼ ∂ y + b W c C ∼ − C ∼ ∞ ∂ ∂ y N ∼ ∂ C ∼ ∂ y = D m ∂ 2 N ∼ ∂ y 2 . (5)

Following are the flow constraints used with the current flow model: u ∼ 1 = U w = a x , u ∼ 2 = 0 , T ∼ = T ∼ w , C ∼ = C ∼ w , N ∼ = N ∼ w   a t   y = 0 , u ∼ → 0 , T → T ∞ , C → C ∞ , N → N ∞ a s   y → ∞ , (6) where u ∼ 1 , u ∼ 2 denote the velocity in the x , y direction, respectively; s being power law index, ν is the kinematic viscosity, Γ is notified as the Williamson parameter, σ denotes the electrical conductivity, ρ f represents the density of the fluid, C p is the specific heat, T ∼ is the temperature of the fluid, σ * represents the Stefan–Boltzmann constant, μ is the dynamic viscosity, k * is the coefficient of mean absorption, D B is the coefficient of Brownian diffusion, k is the thermal conductivity, D T represents the coefficient of thermophoretic diffusion, k r 2 is the chemical reaction rate constant n denotes the fitted rate constants, C ∼ is the fluid concentration, E a is activation energy, κ = 8.61105   e V / K is the Boltzmann constant, T ∼ T ∼ ∞ n exp E a κ T ∼ represents the modified Arrhenius parameter, N ∼ is the microorganism concentration, D m shows the diffusion coefficient of microorganisms, b denotes the chemotaxis constant, and W c is the maximum swimming speed of cell.

For present flow, suitable transformations Equation 7 are derived from [53] ξ = a ν , u ∼ 1 = a x f ′ ξ , u ∼ 2 = − a ν f ξ , θ = T ∼ − T ∼ ∞ T ∼ w − T ∼ ∞ , ϕ = C ∼ − C ∼ ∞ C ∼ w − C ∼ ∞ , χ = N ∼ − N ∼ ∞ N ∼ w − N ∼ ∞ . (7)

The governing equations with the aforementioned variables yield these Equations 8–12 in dimensionless forms: 1 − s f ‴ + f ″ f ‴ s W e − f ′ 2 + f f ‴ − f ′ M + λ θ − N r ϕ − R b χ = 0 , (8) θ ″ 1 P r + f θ ′ + 1 + θ θ w − 1 2 θ θ w − 1 + 3 θ ′ 2 θ w − 1 + 1 θ ″ R P r + 1 − s E c f ″ 2 s 2 W e E c f ″ 3 + θ ′ 2 N t + θ ′ ϕ ′ N b + f ′ 2 M E c = 0 , (9) ϕ ″ + N t N b θ ″ + ϕ ′ S c f − ϕ ⁡ exp − E 1 + δ θ S c σ 1 2 1 + δ θ n = 0 , (10) χ ″ + χ ′ L b f − P e χ ′ ϕ ′ + ϕ ″ χ + δ 1 = 0 , (11) with f ′ ∞ → 0 , f 0 = 0 , f ′ 0 = 1 , θ ∞ → 0 , θ 0 = 1 , ϕ ∞ → 0 , ϕ 0 = 1 , χ 0 = 1 , χ ∞ → 0 , (12) where the Hartmann number is denoted by M = σ B 0 2 ρ a , Prandtl number Pr = ρ c p f ν k , R = 16 σ * T ∼ ∞ 3 3 k k * shows the radiation parameter, W e = a 2 3 2 Γ U ν represents the Weissenberg number, E c = a 2 U 2 c p T f − T ∞ is the Eckert number, λ = G r R e x 2 denotes the variable for mixed convection, temperature Grashof number, N r = g β c C w − C ∞ g β T T w − T ∞ is the buoyancy force, R b = γ 1 ρ m − ρ f N w − N ∞ ρ f β 1 − C ∞ T w − T ∞ is the bioconvection Rayleigh number, the Schmidt number is represented by S c = ν D B , the chemical reaction rate constant is σ 1 = k r 2 a , δ = T ∼ w − T ∼ ∞ T ∼ ∞ is the temperature relative parameter, υ = E a κ T ∼ ∞ represents the activation energy parameter, N t = ρ c p s D T T ∼ w − T ∼ ∞ T ∼ ∞ ρ c p f ν is the parameter for the thermophoresis effect, N b = ρ c p s D T C ∼ w − C ∼ ∞ T ∼ ∞ ν ρ c p f is the Brownian motion parameter, L b is the bioconvection Lewis parameter, P e = b W c D m shows the Peclet number, δ 1 = N ∼ ∞ N ∼ w − N ∼ ∞ is the motile microorganism difference parameter, and θ w = T ∼ w T ∼ ∞ is the temperature ratio parameter.

The following are the mathematical expressions (Equations 13–16) for the motile density number, local Sherwood number, skin friction coefficient, and the local Nusselt number: C f = 2 τ w ρ U w 2 , N u x = x q w k T ∼ w − T ∼ ∞ , S h x = x q m k C ∼ w − C ∼ ∞ , N n x = x q n k N ∼ w − N ∼ ∞ , (13)

with τ w μ 1 − s ∂ u ∼ 1 ∂ y y = 0 + μ s Γ 2 ∂ u ∼ 1 ∂ y y = 0 3 , q w = − ∂ T ∼ ∂ y y = 0 k + q r w , q m = − ∂ C ∼ ∂ y y = 0 D B , q n = − ∂ N ∼ ∂ y y = 0 D m , (14) where q r w = − 16 σ * T ∼ 3 3 k * ∂ T ∼ ∂ y y = 0 . (15)

The above quantities are expressed as follows in non-dimensional forms: 1 2 R e 1 2 C f = f ″ 0 1 − s + s W e 2 f ″ 0 2 , N u x R e − 0.5 = − θ ′ 0 1 + R θ w 3 , S h x R e − 0.5 = − ϕ ′ 0 , N n x R e − 0.5 = − χ ′ 0 , (16) where C f , N u , S h , N n , and R e x are the coefficients of skin friction, the local Sherwood, the local Nusselt, the Reynolds number R e x = a 2 x ν , and motile density, respectively.

For tangent hyperbolic fluid, entropy generation is [58] S G = k T ∼ ∞ 2 1 + 16 σ * T ∼ 3 3 k k * ⏟ Thermal irreversibility + μ T ∼ ∞ Φ ⏟ Fluid Friction irreversibility + σ T ∼ ^ ∞ B 0 2 u 1 2 ⏟ Joule dissipation irreversibility + R D T ∼ ∞ ∂ T ∼ ∂ y ∂ C ∼ ∂ y + R D C ∼ ∞ ∂ C ∼ ∂ y 2 ⏟ Concentration irreversibility + R D N ∼ ∞ ∂ N ∼ ∂ y 2 + R D T ∼ ∞ ∂ N ∼ ∂ y ∂ T ∼ ∂ y , ⏟ Microorganisms concentration irreversibility (17)

where Φ = ∂ u 1 ∂ y 2 1 − s + s Γ 2 ∂ u 1 ∂ y ∂ u 1 ∂ y 2 . (18)

Now Equation 17 becomes using Equation 18 S G = k T ∼ ∞ 2 1 + 16 σ * T ∼ 3 3 k k * ⏟ Thermal irreversibility + μ T ∼ ∞ 1 − s ∂ u 1 ∂ y 2 + s Γ 2 ∂ u 1 ∂ y ∂ u 1 ∂ y ⏟ Fluid Friction irreversibility + σ T ∼ ∞ ^ B 0 2 u 1 2 ⏟ Joule dissipation irreversibility + R D C ∼ ∞ ∂ C ∼ ∂ y 2 + ∂ C ∼ ∂ y ∂ T ∼ ∂ y R D T ∼ ∞ ⏟ Nanoparticle concentration irreversibility + ∂ N ∼ ∂ y 2 R D N ∼ ∞ + R D T ∼ ∞ ∂ T ∼ ∂ y ∂ N ∼ ∂ y ⏟ . Microorganisms concentration irreversibility (19)

Equation 19 contains five factors: i) viscosity dissipation in the entropy generation for tangent hyperbolic fluid, ii) thermal irreversibility, iii) nanoparticle concentration irreversibility, iv) Joule effect irreversibility, and v) microorganism concentration irreversibility. The dimensionless form are in Equation 20 N g = R θ θ w − 1 + 1 3 + 1 α 1 θ ′ 2 + 1 − s f ″ 2 B r + f ″ 3 B r W e 2 + B M r f ′ 2 + L α 2 α 1 ϕ ′ 2 + L θ ′ ϕ ′ + L * α 3 α 1 χ ′ 2 + L * θ ′ χ ′ , (20) α 1 = Δ T ∼ T ∼ ∞ = T ∼ w − T ∼ ∞ T ∼ ∞ , α 2 = Δ C ∼ C ∼ ∞ = C ∼ w − C ∼ ∞ C ∼ ∞ , α 3 = N ∼ w − N ∼ ∞ N ∼ ∞ = Δ N ∼ N ∼ ∞ , B r = μ a 2 k Δ T ∼ , N g = T ∼ ∞ S G ν a k Δ T ∼ , L = R D C ∼ w − C ∼ ∞ k , L * = R D N ∼ w − N ∼ ∞ k , (21) where Equation 21 represents α 1 , α 2 , α 3 , B r , L , L * , and  N g differences in the temperature variable, Brinkman number, concentration difference parameter, diffusion parameter due to nanoparticles and microorganisms’ concentration, and local entropy generation. The dimensionless Bejan number is expressed in Equation 22 and mathematical form is in Equation 23 B e = Entropy   generation   due   to   mass   and   heat   transfer Total   entropy   generation , (22) B e = 1 + R θ θ w − 1 + 1 3 θ ′ 2 α 1 L θ ′ ϕ ′ + L α 2 α 1 ϕ ′ 2 + L θ ′ ϕ ′ + L * α 3 α 1 χ ′ 2 + L * θ ′ χ ′ 1 − s B r f ″ 2 + 1 + R θ θ w − 1 + 1 3 θ ′ 2 α 1 + W e 2 f ″ 3 B r + B M r f ′ 2 + α 2 α 1 L ϕ ′ 2 + θ ′ ϕ ′ L + L * α 3 α 1 χ ′ 2 + L * θ ′ χ ′ . (23)

It is a well-known observation that a number of physical models contain differential equations which are nonlinear depending on the included constraints that can be solved using numerical approaches. In the suggested numerical approaches, the shooting method is perhaps a more reliable strategy for achieving adequate accuracy. After generating dimensionless forms, the bvp4c built-in package is used to tackle numerically. This approach is reliable and efficient, having an error tolerance of 10⁻⁶. As a result, we used this strategy to target Equations 8–12. The collection of these mathematical equations is turned into first-order differential equations using dummy variables as in Equations 24–28: f = z 1 , f ′ = z 2 , f ″ = z 3 , f ″ = z 3 ′ , θ = z 4 , θ ′ = z 5 , θ ″ = z 5 ′ , ϕ = z 6 , ϕ ′ = z 7 , ϕ ″ = z 7 ′ , χ = z 8 , χ ′ = z 9 , χ ″ = z 9 ′ .

Then, z 3 ′ = z 2 2 − z 1 z 3 + M z 2 − λ z 4 1 + β t z 4 − λ z 6 N * 1 + β c z 6 − λ z 8 N * * 1 + β n z 8 1 − s + s W e z 3 , (24) z 5 ′ = − 3 R z 4 θ w − 1 + 1 2 z 5 2 θ w − 1 − Pr ⁡ z 1 z 5 − 1 − s E c 2 z 3 5 s 2 W e P r − N t P r z 5 2 − N b P r z 5 z 7 − M z 2 2 E c P r 1 + R z 4 θ w − 1 + 1 2 z 4 θ w − 1 + 1 , (25) z 7 ′ = − N t N b z 5 ′ − S c z 1 z 7 + S c σ 1 2 1 + δ z 4 n exp − υ 1 + δ z 4 z 6 , (26) z 9 ′ = − L b z 1 z 9 + P e z 7 ′ z 8 + δ 1 + z 7 z 9 . (27)

The transformed boundary conditions are z 1 0 = 0 , z 2 ∞ → 0 , − 1 + z 2 0 = 0 , − 1 + z 4 0 = 0 , z 4 ∞ → 0 , z 6 0 − 1 = 0 , z 6 ∞ → 0 , z 8 0 − 1 = 0 , z 8 ∞ → 0 . (28)

This part describes the visual representation of physical characteristics such as velocity distribution f ′ , temperature θ , nanoparticle concentration ϕ , and motile gyrotactic χ . In the process of modifying the related physical parameter, the quantities that are left have been allotted constant values such as M = 0.5 , s = λ = β c = W e = σ 1 = δ = E c = 0.1 , S c = N b = N t = N * = 0.01 , ⁡ Pr = 2 , and θ w = 1.1 . Figures 2–4 examine the behavior of a variety of physical factors, including the Hartmann number M , Weissenberg number W e , first-order velocity slip α , buoyancy ratio parameter N r , second-order slip parameter β , and mixed convection variable λ . The effects of the Hartmann number M and W e Weissenberg number on the velocity graph f ′ are depicted in Figure 2. A declining behavior is prominent for velocity distribution by increasing both M and W e . For higher values of (M = 0.1, 0.4, 8, 1.2), Figure 2 shows a decreasing trend of f ′ ξ . Because M depends on the Lorentz force, the inter-particle resistance increases as M increases, which reduces f ′ ξ . We observe that as W e increases, fluid motion gradually decreases. Given that the fraction of relaxation time by a particular processing time is denoted by W e , larger W e will give rise to a longer relaxation period, which will lessen the velocity profile. Figure 3 discusses the first-order slip parameter α . In addition, the mixed convection parameter λ for the profile of velocity f ′ ξ . First, it should be noted that the first order is the feature of the medium that is linked to the pathways between fluid flows. The values for the velocity declines as we change first-order slip parameters. Moreover, it is noticed that the thickness of the boundary, which is entitled to first-order slip factor, is decreasing. The velocity profile increases with the increase in the mixed convection parameter λ . Physically, the buoyant force overcomes the force, owing to the increase in inertia as the parameter for mixed convection, and as a result, a leading increase is noticed for f ′ ξ .

The parameter for the buoyancy ratio N r and the parameter for the second-order slip β for the dimensionless distribution of velocity f ′ ξ are graphed in Figure 4. The profile for velocity decreases as both of these factors increase.

The impacts of the Prandtl number Pr , Biot number B i , heat sink/source Q , parameter for thermophoresis N t , temperature ratio parameter θ w , and Eckert number E c on θ temperature field are depicted in Figures 5–7. Figure 5 shows the physical impact caused by Pr and B i on temperature field θ . It is easy to observe that θ is the decreasing function of ( Pr = 2.0, 3.0, 4.0, 5.0). Furthermore, the thermal Biot number B i has a connection to the coefficient of heat transfer that causes the increase in the temperature distribution of nanoparticles. The Eckert number E c and the characteristics of parameters for the temperature ratio θ w and on the profile for temperature θ are demonstrated in Figure 6. In fluid temperature θ , an increasing trend is prominent for the increase in the values of parameters for the temperature ratio ( θ w = 1.5, 1.6, 1.7, 1.8). For increasing values of E c , it displays that the fluid temperature shows an increase. The reason for it being higher E c fluid friction results in an interconversion of mechanical energy and thermal energy, which peaks the temperature θ of the fluid.

The effects of the Q < 0 heat sink or Q > 0 heat source and thermophoresis parameter N t are shown in Figure 7. In the boundary layer, when a heat source is present, it increases energy, which results in the increase in fluid’s temperature. The existence of a heat sink in the boundary layer absorbs energy, which, in turn, lowers the temperature of the fluid. These behaviors are depicted in Figure 7. For larger N t , particles of fluid increase from hot to cold areas of the system. The temperature profile increases as a result of an increase in thermophoresis force.

To examine the changes in concentration distribution ϕ versus different values of the E activation energy, Pr Prandtl number, L e Lewis number, and N b Brownian movement, Figures 8, 9 are formulated. Figure 8 presents the graph of the Prandtl number Pr and Lewis number L e on ϕ . Increases in Pr and L e result in decreases in the profile of concentration of nanoparticles. The ramifications of yet two other crucial parameters, the parameters for Brownian motion N b and activation energy E , are shown in Figure 9. Although the concentration of nanoparticles decreases with N b , but for E , the profile for the concentration improves.

The Peclet number P e and Lewis number L b for bioconvection effects on the self-propelled microorganism profile χ are shown in Figure 10. By altering Pe, a retarding behavior is discovered; this tendency establishes because P e develops a reverse relationship with microorganism diffusivity. χ reduced as a result of this inverse relationship. A distinguishing motile microorganism profile is found when the value of L b increases. The physical reasoning is that a bigger change in L b is connected with the reduced diffusion of self-propelled microorganisms, and χ decays as a result of this. Figures 11, 12 examine the behavior of a variety of physical factors, including the Hartmann number M , power law index s , and N b parameters for Brownian motion on entropy generation N g . Salient features of the Hartmann number M and Weissenberg number W e for entropy generation N g are displayed in Figure 11. The increase in resistance for increasing M entropy generation shows an increasing trend. Relaxation time increases for higher W e , which indicates increased heat loss due to increased resistance between liquid particles. Consequently, entropy generation increases. Figure 12 is plotted to examine the entropy generation for greater values of s , which is the power law index and Brownian motion parameter N b . For the increase in s and N b , entropy generation N g shows decreasing behavior.

The streamlines of the flow problem visualize the flow pattern and movement of the molecules, as displayed in Figures 13A–D. The laminar pattern is visible among all profiles, and no intersection is found. These streamlines describe the movement, velocity distribution, heat transport efficiency, and interaction with the surface body.

In the tabulated values in Table 1, we see the comparison of our computed results with the already in the tabulated values we see the comparison of our computed results with the already results shows the great accord with available literature.

The aim of this contribution is to study the flow of tangent hyperbolic nanomaterials for entropy optimization, owing to a stretchable surface. Some key findings are as follows:

• For the first-order velocity slip factor, the velocity distribution increases, but the Hartmann number, Weissenberg number, buoyancy ratio constant, and second-order slip parameter show a declining trend.

In the future, this approach may be extended to include the various heat and mass flux models, as well as the observation of changes in transport rates due to geometrical changes.