Cattaneo–Christov Double Diffusion (CCDD) on Sutterby Nanofluid with Irreversibility Analysis and Motile Microbes Due to a RIGA Plate

In this article, a Riga plate is exhibited with an electric magnetization actuator consisting of permanent magnets and electrodes assembled alternatively. This Riga plate creates an electric and magnetic field, where a transverse Lorentz force is generated that contributes to the flow along the plate. A new study field has been created by Sutterby nanofluid flows down the Riga plate, which is crucial to the creation of several industrial advancements, including thermal nuclear reactors, flow metres, and nuclear reactor design. This article addresses the second law analysis of MHD Sutter by nanofluid over a stretching sheet with the Riga plate. The Cattaneo–Christov Double Diffusion heat and mass flux have been created to examine the behaviour of relaxation time. The bioconvection of motile microorganisms and chemical reactions are taken into consideration. Similarity transformations are used to make the governing equations non-dimensional ordinary differential equations (ODE’s) that are subsequently solved through an efficient and powerful analytic technique, the homotopy analysis method (HAM). The effect of pertained variables on velocity, temperature, concentration, and motile microorganism distributions are elaborated through the plot in detail. Further, the velocity distribution enhances and reduces for greater value Deborah number and Reynold number for the two cases of pseudoplastic and dilatant flow. Microorganism distribution decreases with the augmented magnitude of Peclet number (Pe), Bioconvection Lewis number (Lb), and microorganism concentration difference number (ϖ). The entropy production distribution is increased for the greater estimations of the Reynolds number (ReL) and Brinkman parameter (Br). Two sets of graphical outputs are presented for the Sutterby fluid parameter. Finally, for the justification of these outcomes, tables of comparison are made with various variables.


Introduction
The rate of heat transport characteristics is a subject of enduring interest for various researchers and scientists owing to its tremendous industrial applications, for example, in mechanical, optical, electrical, and cooling instruments. The rate of heat transport increment is very crucial in depositing energy. Thus, researchers are focused on the investigation of a new fluid substance that is a mixture of nanoparticles 100 nm in size and more thermophysical properties than ordinary fluid, known as nanoliquid. Nanoliquid is the colloidal suspension of the nanoparticles' thermal behaviour of the ordinary fluid. The first endeavour was performed by Choi et al. [1] in 1995. They exposed the thermal conductivity of nanoliquid by adding nanosized particles. Later, Buongiorno [2] used this notation of nanofluid to make a mathematical form by adding terms known as Brownian and Thermophoretic. The mixed convection nanofluid flows with different surfaces and geometry were conducted by Hussain et al. [3,4], and Haq et al. [5] deliberated the second law analysis on a cross nanofluid. The impact of microorganisms on the Darcy flow of Sutterby nanofluid, taking into account microorganisms over movable cylinders, was examined by Aldabesh et al. [6]. Mankiw et al. [7] analysed the MHD time-dependent flow of nanofluid with variable properties due to an inclined stretching sheet under thermal radiation. Shahid [8] demonstrated the effect of upper convective flow of Maxwell liquid over a permeable surface near the stagnation point. Rafique et al. [9] addressed the stratified micropolar nanofluid flow past an exponentially stretchable plate with the Riga surface. The unsteady viscous flow of nanofluid over the Riga plate using a rotating plate was investigated by Parvine et al. [10]. Abbas et al. [11] deliberated the entropy production over the Riga plate and the suction case. Sannad et al. [12] using a non-homogeneous dynamic model, which is more accurate in physically describing nanofluids than homogeneous ones, wherein free convective flow in a cubical cavity filled with copper-water nanofluid was explored numerically. Some of the latest developments in nanofluids were obtained from references [13][14][15][16][17][18][19].
Different fluid forms, including polymer melts, colloidal suspensions, organic chain mixes, etc., are used in a wide range of industries and production processes. The rheological behaviour of these fluids cannot be well described by the Naiver-Stokes equation alone. Therefore, several nonlinear fluid models are proposed to represent the rheological characteristics of complicated fluids. One non-Newtonian fluid model used to examine key characteristics of pseudoplastic and dilatant fluids is the Sutterby fluid model. Numerous experts have extensively studied the flow of Sutterby liquid. Waqas et al. [20] inspected Sutterby nanofluid using two rotating disks. Yahya et al. [21] investigated Williamson Sutterby nanoparticles under the Cattaneo-Christov heat flux. The effect of MHD on Sutterby nanoparticles due to porous movable sheets was discovered by Fayydh et al. [22]. Gowda et al. [23] examined the Cattaneo-Christof theory of heat diffusion in Sutterby nanofluid. The impact of thermal radiation and heat source/sink on Sutterby nanofluid has been studied by Ali et al. [24]. Hayat et al. [25] investigated the Sutterby fluid with thermal radiation due to a rotating disk. Recently, Fujii et al. [26] addressed Sutterby fluid with natural convection flow due to a vertical plate. The Darcy surface with MHD flow of Sutterby fluid was reported by Bilal et al. [27]. The bioconvection flow of Sutterby nanofluid due to a rotating disk was described by Khan et al. [28]. Sohail et al. [29] designed the free convection flow of Sutterby fluid with heat under the Cattaneo-Christov theory. The heat generation/absorption in the thermally stratified flow of stagnant Sutterby fluid through a linearly stretched plate is analysed in this article by Saif et al. [30]. Usman et al. [31] investigated the two-dimensional stagnant flow of Sutterby nanofluid across a stretching wedge placed in a porous media. Song et al. [32] discussed the two-dimensional flow of Sutterby nanofluids over a stretching cylinder, scrutinized in the presence of bioconvection and swimming microorganisms. The influence of homogeneous heterogeneous reaction on Sutterby fluid flow by a disk with Cattaneo-Christov heat flux is studied by Khan et al. [33].
In the modern period, the research on bio-convection exists due to the motion of microorganisms upwards because microorganisms are denser than water. The upward surface of the fluid develops thickness as a result of the collection of the microorganism. For this reason, the upper surface becomes disturbed, and microorganisms sink, which develops bioconvection. Bio-convection continues to be the subject of research broadly due to its many applications in the clinical area, manufacturing process, and biofuel production. Bio-convection could be organised as consisting of the motion of direction with a vast number of microorganism species. In this way, gyrotactic microorganisms are among those whose swimming direction is based on viscous torque and gravitational force. Kuznetsov et al. [34,35] reported the early investigation of bio-convection in the mixed suspension of nanoparticles with gyrotactic microorganisms. Kotha et al. [36] examined the MHD flow of nanofluid with motile gyrotactic microorganisms over a vertical plate. Siddiq et al. [37] numerically analysed the bioconvection of micropolar nanofluid flow restricted through the stretchable disk using the bvp4c method. Bagh et al. [38] studied the effect of bioconvection and Cattaneo-Christov with a vertical stretching sheet. Azam et al. [39] investigated the effect of bioconvection flow for the Sutterby nanoliquid with nonlinear radiation. Khashi'ie et al. [40] tested a hybrid nanofluid that had bioconvection with gyrotactic microorganisms. Azam [41] explored the time-dependent flow of the chemically reactive Sutterby nanofluid and the influence of gyrotactic microorganisms. Hayat et al. [42] operated the bio-convection flow of nanomaterial subject to melting effect. They addressed thermal nonlinear radiation and Joule heating for heat distribution characteristics. Reddy et al. [43] analysed the time-dependent flow of Cross nanofluid comprising the gyrostatic microorganisms due to slip velocity. Sarkar et al. [44] defined Sutterby nanofluid flow as having motile gyrotactic microorganisms over the Riga plate. Shah et al. [45] described the heat transfer properties of a magnetohydrodynamic Prandtl hybrid nanofluid over a stretched surface in the presence of bioconvection and chemical reaction effects. Ali et al. [46] explored the fluctuating temperature with Cattaneo-Christov features and selfmotivated bioconvective microbes immersed in the water-based nanofluid, observed with the excision/accretion of the leading edge.
Riga plate is an electromagnetic actuator, which is made by the combination of a span-wise aligned array of alternating electrodes and permanent magnets. This actuator is responsible for producing a Lorentz force that decreases exponentially. It is important to note that when an external magnetic field is applied to a fluid with high electrical conductivity, it has a provocative effect. However, this device has been used to reduce pressure drag by preventing boundary layer separation and reducing the generation of turbulence requirements. Riga plate was first introduced by Gailitis and Leilausis [47], and they generated a wall paralleled by Lorentz force to control the fluid flow. Further, in contrast to this, the Grinberg term was introduced with a most important feature that the boundary layer of the momentum equation is fully separate from the flow, and exponentially decreases in the direction of the normal to the plate. Subsequently, the Lorentz force was detected over a Riga plate, and then the researchers regained interest in the Gailitis-Lielausis actuator. Islam and Nasrin [48] described the one-dimensional unsteady micropolar fluid flow set in a porous medium along with an inclined infinite Riga plate. Nasrin et al. [49] explored the impulsively started horizontal Riga plate in two-dimensional unsteady Casson fluid flows with rotation. Anjum et al. [50] concentrated on thermal stratification in the flow of second-grade fluid past a Riga plate with linear stretching toward a stagnation region. Hayat et al. [51] investigated the convective heat transfer of electron magnetohydrodynamic squeezed flow past a Riga plate. Wahidunnisa et al. [52] stduied a mixed convective nanofluid flow along a heated Riga plate with viscous dissipation and heat source. Iqbal et al. [53,54] have found the stagnation point flow of the melting heat, thermal radiation, and viscous dissipation effects on the Riga plate using erratic thickness, and considered the Casson fluid stagnation point flow along with a Riga plate. Vishnu et al. [55] discussed the γ Al 2 O 3 -Water/Ethylene Glycol over a Gailitis and Lielausis device with an effective Prandtl number.
Newtonian fluids have been mostly discussed as base fluids in the prior work. We were motivated to examine Sutterby as a base fluid along with the Riga plate and chemical reaction effects with Cattaneo-Christove double diffusion for microorganism flow of nanofluids due to stretching surface because there have been so few studies on them. These facets of the issue, as far as the authors can tell, are not taken into consideration in the inves-tigations that have already been done. The improvement of heat and mass transportation is the main goal of this extensive investigation. The current investigation was encouraged to express the Sutterby nanofluid inserted on a porous surface with Cattaneo-Christof double diffusion due to the Riga plate containing gyrotactic microorganisms. The chemical reaction and heat source-sink are considered worthy of attention. The leading aim of this work is an inclusive analysis of this flow problem. The coupled equations of the current flow problem are altered into the nonlinear system by applying suitable correspondence transformations. Further, the solution of ordinary differential equations is utilized via the Homotopy Analysis Method (HAM). The novel outcomes of the current work are obtained through different parameters and explained with the assistance of graphs and tables.

Description of the Physical Model
Consider the incompressible and steady flow of second law analysis in Sutterby nanofluid due to the Riga Plate containing gyrotactic microorganisms displayed in Figure 1. Cattaneo-Christov with heat and mass flux has also been pondering the temperature and concentration equation. The x-axis is considered along with the sheet and the y-axis is taken perpendicular to the sheet. Moreover, the sheet of velocity is U w = ax. At the surface of the sheet, the temperature of the surface, the concentration of the surface, and the microorganisms of the surface are symbolized by T w , C w , χ w . The swimming route and velocity of microorganisms do not alter when nanoparticles are present. However, if the volume percentage of nanoparticles exceeds 1%, the motility of microorganisms is impacted. As a result, the base liquid is blended with microorganisms and solid nanoparticles to provide the necessary bioconvection stability. Furthermore, it has been assumed that the fluid contains gyrotactic bacteria. The fluid microorganisms gravitate towards the light. Gyrotactic phenomena, or movement against gravity, are made possible by the "bottom heavy" bulk microorganism, which orients their bodies. The existence of microorganisms is advantageous for the suspension of the nanoparticles. The motion of microorganisms has been taken, irrespective of that of the nanoparticles, to ensure the stability of convection. The flow of a double-diffusive fluid across a stretching sheet containing gyrotactic microorganisms has not yet been investigated, and this study aims to fill that gap with the simplifying of boundary layer approximations of the leading expressions given by [23]:

Fluid Flow Problem
Cauchy Stress tensor  for the Sutterby fluid [56,57]   1 pI, Sutterby viscosity model is represented as

Fluid Flow Problem
Cauchy Stress tensor τ for the Sutterby fluid [56,57] Sutterby viscosity model is represented as where n, µ • , . β, are described as power law index, zero share rate viscosity, and time material constant.
Using Equation (2) in Equation (1), we have The governing equation of the following form [23] ∂u ∂x where u and v are velocity components in the x and ydirections. ν represents kinematic viscosity of the fluid, υ represents kinematic viscosity of the fluid, ρ represents the density of the fluid, α represents thermal diffusivity, C is the concentration, D B , D t represent Brownian diffusion and thermophoretic diffusion, C p denotes volumetric expansion, D m represents the microorganism coefficient, Ω E and Ω C fluid relaxation time, In the above equations, the terms Ω E Ω C are stated as The relevant boundary conditions are assumed to be the form where T 4 can be expanded as follows: Replacing Equation (12) into Equation (13), Introducing the similarity variables [24] Using Equation (15), Equations (4)- (6) become The corresponding boundary conditions are where is the Modified Hartmann number; Re L = ax 2 v is the Reynolds number; a is the Biot number; Pe = bW c D m is the Peclet number; Lb = v Dm is the Bio convection Lewis number; • ρa is the Magnetic number. The thermo-fluidic quantities of engineering interest in this study are skin friction C f x , heat transfer rate Nu x , mass transfer rate Sh x , and motile density Wh x .
where τ w surface shear stress, q w surface heat flux, q m surface mass flux, and q n motile density are given by the following expressions: The dimensionless form of the above parameters is expressed as Here, Re x = xU w v is the local Reynolds number.

Entropy Generation Analysis
Entropy generation with Sutterby fluid is communicated as [39]: From Equation (23) the entropy production consists of some terms. The first term represents heat transfer, the second-term fluid friction, and the third to sixth terms represent diffusive irreversibility. The significance of entropy production can be written as Using Equation (15) the rate of entropy production with Sutter by Equation (23) can be converted The Bejan number Be is defined as the ratio over the entropy generation with heat transport S T and the total entropy production N G , and it can be written as:

Homotopy Expression
HAM was used to investigate the explicit and analytic solutions of Equations (16)-(19) [51]. This analytical technique is quite effective and flexible in solving the parabolic-type boundary value problems of any order, is unconditionally stable, and attains remarkable accuracy. Taking the initial guesses and the linear operators as with the propertyL in which B i , (i = 1 − 9) are the constants. For s = 0 and s = 1, the results are achieved

The mth Order Formulation
The mth order deformation can be presented in the following formŝ Boundary conditions are where The general solutions are in which f * m , θ * m , φ * m , W * m are the special solution.

Convergence of Homotopy Solutions
These parameters f , θ , φ and W are the converging control of the desired series solution. For the function f (0), θ (0), φ (0), W (0) to seek the permissible values to obtain the 25th and 30th order. Figures 1-4

Results and Discussion
The system of Equations (16)- (19) subject to the boundary conditions (20) has been tackled through the analytical technique known as the Homotopy analysis method (HAM) for distinct parameters that emerged in the analytical simulation. To discuss the

Results and Discussion
The system of Equations (16)- (19) subject to the boundary conditions (20) has been tackled through the analytical technique known as the Homotopy analysis method (HAM) for distinct parameters that emerged in the analytical simulation. To discuss the perfor-mance of physical significance against the velocity field f (η), temperature distribution θ(η), concentration field φ(η), motile microorganism profile W(η), entropy production, Bejan number, as well as skin friction, Nusselt number, Sherwood number and motile density microorganism, are delineated in Figures 3-9. Table 2 verifies −θ(0) with Ali and Zaib [55], and found a good agreement.  [55], and found a good agreement.   [55], and found a good agreement.    , , Pe Lb .

Results and Discussion
The system of Equations (16)- (19) subject to the boundary conditions (20) has been tackled through the analytical technique known as the Homotopy analysis method (HAM) for distinct parameters that emerged in the analytical simulation. To discuss the performance of physical significance against the velocity field   [58], and found a good agreement.

Velocity Field
The effects of different numerous parameters on velocity distribution Figure 3a-d. Figure 3a shows that for pseudoplastic, the velocity of fluid lessens with the larger value of Deborah number  , and the velocity field is enhanced for the rising value of Deborah number  for the two cases of when 0   and 0   . This result is shown to be the same as in F. Ali et al. [56].

Results and Discussion
The system of Equations (16)- (19) subject to the boundary conditions (20) has been tackled through the analytical technique known as the Homotopy analysis method (HAM) for distinct parameters that emerged in the analytical simulation. To discuss the performance of physical significance against the velocity field f (η), temperature distribution θ(η), concentration field φ(η), motile microorganism profile W(η), entropy production, Bejan number, as well as skin friction, Nusselt number, Sherwood number and motile density microorganism, are displayed graphically and in tabular form. Table 2 verifies −θ(0) with Ali and Zaib [58], and found a good agreement.

Velocity Field
The effects of different numerous parameters on velocity distribution f (η) are discussed in Figure 3a-d. Figure 3a shows that for pseudoplastic, the velocity of fluid lessens with the larger value of Deborah number δ, and the velocity field is enhanced for the rising value of Deborah number δ for the two cases of when β < 0 and β > 0. This result is shown to be the same as in F. Ali et al. [56]. Figure 3b shows that the velocity curve f (η) shortens with increasing values of the Reynold number Re L in the case of β < 0; a higher Re L , the more diminished the viscous force and hence fluid depreciates for the pseudoplastic fluid. Thus, for shear thickening fluid, the velocity field enhances as Re L increasing the impact of the values of β in the two cases. As shown in the figure, the velocity for both values β shows mixed behaviour with increasing fluid parameters. The Reynold number described the proportion of inertial to viscous forces due to this inertial force being prominent for the two cases of velocity magnitude when (β = 2.5) and increases when (β = −2.5). Figure 3c manifests the effect of the enhanced power-law index β parameter, resulting in a better velocity profile for shear thickening fluid. Figure 3d demonstrates the significance of the Hartmann number Z on velocity fluid field for the two cases of β < 0 and β > 0. It has been revealed that with changes in the strength of Z the velocity of the fluid escalates in both cases. Physically, an increment in Z corresponds to enhancing the external electric field that constructs wall-parallel Lorentz force. Therefore f (η) rises. The velocity is more pronounced for pseudo-plastic fluid than a fluid dilatant fluid with increasing Z for the shear-thinning as well as shear-thickening fluid.

Temperature Field
Figure 4a-f plotted the consequences of temperature θ(η) against different values of various involved parameters. The thermal relaxation time parameter causes are demonstrated in Figure 4a. The thermal relaxation parameter tends to decrease the temperature curve in both cases of dilatants and pseudoplastic. It impacts the strength of heat transfer. So, the dissimilar decayed impact of λ 1 on the temperature field takes place away from the surface. Figure 4b reveals the inclination of θ(η) for specific values of thermal radiation parameters fo β > 0 and β < 0 The temperature of the fluid increases due to the enlargement in the radiation parameter Rd in both cases. Physically, by increasing the radiation parameter more energy is added to the system, which uplifts the kinetic energy of molecules and also causes a rise in the temperature profile. Figure 4c shows that temperature increases with the increments in the values of Nt. Thermophoresis is a phenomenon where small particles diffuse under the effect of a temperature gradient. An increase in the value of the thermophoresis parameter sets the nanoparticle in motion, and thus on the hotter side, the momentum of the nanoparticles rises. Due to the rise in the momentum, nanoparticles transfer their kinetic energy in the direction of the cooler side, and thus the region gets warmed up quickly; i.e., the temperature of the fluid is increased. Figure 4d also suggests that the nanofluid temperature rises with the increasing value of the Brownian motion parameter Nb, and the opposite effect is examined in the Figure 4e. The reason behind this is that Brownian motion is the irregular movement of the particles suspended in the fluid. The temperature of the fluid increases as a result of the random collision of particles suspended in the liquid, which further leads to an expected improvement in the temperature profile θ(η) Figure 4e depicts the impact of the Biot number θ(η). These figures point out that the temperature field is boosted by enhancing the value of the Biot number. It seems that the value of the Biot number is directly proportional to the temperature, as the higher value of Bi corresponds to greater convective heating along the sheet, which increases the temperature gradient on the sheet. Therefore, the thickness of the boundary layer and temperature are increasing functions of the Biot number. Figure 4f display variations of the Prandtl number against θ(η). Temperature is restrained against higher Pr. Basically, the Prandtl number refers to thermal diffusivity. Larger Pr indicates lower thermal diffusivity, which causes the temperature to decompose.

Concentration Field
The outcomes of different leading parameters φ(η) are presented in Figure 5a-d. Figure 5a is designed to show the characteristics of Nb on φ(η). The concentration distribution is depleted by rising Nb. The reason behind the decrease in the concentration field is the relation of Brownian motion with the Brownian diffusion coefficient which is responsible for decreasing the concentration field. The influence of the thermophoresis variable φ(η) is rendered in Figure 5b. It can be seen that increasing values of the thermophoresis parameter enhances mass concentration. This is because the increase in Nt corresponds to the increase in thermophoretic diffusion coefficient enhancing mass concentration. The influence of chemical reaction on the profile of concentration is shown in Figure 5c. The enhanced values Cr result in a fluid particle breakthrough near the surface, which reduces the concentration and the corresponding boundary layer thickness. The further rate of inclination is slower in the presence of the chemical reaction rate parameter. The Schmidt number outcomes for concentration are displayed in Figure 5d. A depreciation in concentration arises with greater Sc due to the reduction of mass diffusion. In the physical phenomenon of sight, an increase in the Sc lessens the molecular diffusivity and results in the declination of the concentration gradient.

Microorganism Field
The effect of different influential variables on the microorganism's field is shown in Figure 6a-c. Variations in motile microorganisms against the Lewis number Lb for various values are seen in Figure 6a. The Lb reduces the density of microorganisms. Actually, Lb have an opposite trend with thermal diffusivity, as an escalation in Lb decreases the thermal diffusivity with regard to a decline in motile density. Figure 6b points out the behaviour of the Peclet number in the microorganism field. The motile density profile decreases with an increase in Pe. Actually, the Peclet number is the quotient of heat transfer by fluid motion to the heat transfer through thermal conduction, and the microorganism's diffusivity is in an inverse relationship with the Peclet number. That is why the motile density field is on the decline for large estimates of Pe The relationship between the microorganisms' difference parameter and the motile density is present in Figure 6c. It has been noticed that motile density shrinks for larger . In fact, improving the values escalates the concentration of microorganisms in the ambient concentration, and so W(η) declines.

Entropy Generation
Figure 7a-c examined the performance of numerous variable parameters δ Re L and Br entropy production N G . Figure 7a sketched the effect of Deborah's number δ versus the entropy generation number N G . This shows an enhancement in entropy production close to the wall for dilatant β < 0, and deduction close to the wall for pseudo-plastic fluid β < 0. The increasing influence of the Reynolds number Re γ in entropy generation studied for dilatant and pseudoplastic fluids is plotted in Figure 7b. The disparity Br has been plotted in Figure 7c. Entropy generation is boosted with increasing values of Br in both cases. Subsequently, Br attributes a proportion of the free heat through viscous heating to the molecular condition. Therefore, heat is created in the system for increased values of Br while the disorderliness of Re also rises in the system, which increases the entropy of the system.

Bejan Number
The performance Be with the variations in the variables δ, Re γ and Br, Be plotted in Figure 8a-c. Figure 8a,b signify the behaviour of the physical parameters Deborah's number and the Reynold number on the Bejan number. It is shown that the Bejan number declines as the larger values of Deborah's and the Reynold number for shear-thickening and upsurge in both numbers for shear-thinning fluid. Furthermore, Figure 8c shows that the influence of the Bejan number Be is reduced β > 0 and β < 0 with the growing values of Br. Actually, the thermal conductivity of the fluid declines with the rise in Br, so that a greater amount of heat is transmitted through the fluid. The parameter Br characterizes the heat generated by viscous dissipation. So as a result, the Bejan number is reduced via the uplift value of the Brickman number Br. Figure 9a shows that the skin friction coefficient Re 1/2

Physical Entitles
x C f x is deformed in both cases for the larger values of the parameter when β > 0 and β < 0. The influence of the non-Newtonian nanofluid parameter on the Nusselt Re −1/2 x Nu against thermal radiation Rd are highlighted in Figure 9b. The heat transport gradient upsurges with rises in the thermal radiation Rd and heat source/sink parameters Υ. The significance of Nt and Nb on the Sherwood number is intimated in Figure 9c. It is determined that there is an amplification in the Sherwood number Re −1/2 x Sh for raised values of fluid parameters. Figure 9d elucidates the substantial rescaled density number of motile microorganisms. The rescaled density number of motile microorganisms is voluminous for higher variations of Pe and Lb. Figure 10a-d shows the 3D representation of the Skin, Nusselt, Sherwood and Motile density respectively.
The numerical estimation of drag friction , , Re , L Z A  has been given in Table 3. The larger the magnitude of , , Re , L Z A  , the lower the drag friction. Tables 4 and 5 produce the data of heat and mass related to the diverse value of , , Pr, , Rd Bi Cr Sc . Here it is illustrated that the magnitude of heat and mass transport is increased via a large value of  for the Nusselt number Nt and for the Sherwood number. Table 6 shows the variation of , , Pe Lb  on the motile density profile. In Table 6, it is shown the motile density profile is the leading function Lb .  The numerical estimation of drag friction Z, δ, Re L , A has been given in Table 3. The larger the magnitude of Z, δ, Re L , A, the lower the drag friction. Tables 4 and 5 produce the data of heat and mass related to the diverse value of Rd, Bi, Pr, Cr, Sc. Here it is illustrated that the magnitude of heat and mass transport is increased via a large value of Υ for the Nusselt number Nt and for the Sherwood number. Table 6 shows the variation of Pe, Lb, on the motile density profile. In Table 6, it is shown the motile density profile is the leading function Lb.    Figure 11a,b exhibits the behaviour of the stream function for the current flow. The patterns show that the streamlines are more obscured and split into two sections pseudoplastic β < 0 and dilatant β > 0; the shape is modest and fills the floe field. Figure 12a This investigation examines MHD Sutterby nanofluid with Cattaneo-Christov Double Diffusion (CCDD) dual diffusion theory due to the Riga Plate, as well as the bioconvection of motile microorganisms; a chemical reaction is also analysed. Appropriately modified transformations are invoked to get a non-linear system of differential problems. The model under consideration consists of both gyrotactic microorganisms and nanoparticles. By using bio-convective flow, which is produced by the combined impacts of nanoparticles and buoyancy forces, microorganisms maintain the suspension of nanoparticles. Refs. [59][60][61][62][63] are extended in future with different flow assumptions. Along with Newtonian heating, the mechanics of Brownian motion and thermophoresis are also considered. The main findings of the current study are the following: This investigation examines MHD Sutterby nanofluid with Cattaneo-Christov Double Diffusion (CCDD) dual diffusion theory due to the Riga Plate, as well as the bioconvection of motile microorganisms; a chemical reaction is also analysed. Appropriately modified transformations are invoked to get a non-linear system of differential problems. The model under consideration consists of both gyrotactic microorganisms and nanoparticles. By using bio-convective flow, which is produced by the combined impacts of nanoparticles and buoyancy forces, microorganisms maintain the suspension of nanoparticles. Refs. [59][60][61][62][63] are extended in future with different flow assumptions. Along with Newtonian heating, the mechanics of Brownian motion and thermophoresis are also considered. The main findings of the current study are the following: This investigation examines MHD Sutterby nanofluid with Cattaneo-Christov Double Diffusion (CCDD) dual diffusion theory due to the Riga Plate, as well as the bioconvection of motile microorganisms; a chemical reaction is also analysed. Appropriately modified transformations are invoked to get a non-linear system of differential problems. The model under consideration consists of both gyrotactic microorganisms and nanoparticles. By using bio-convective flow, which is produced by the combined impacts of nanoparticles and buoyancy forces, microorganisms maintain the suspension of nanoparticles. Refs. [59][60][61][62][63] are extended in future with different flow assumptions. Along with Newtonian heating, the mechanics of Brownian motion and thermophoresis are also considered. The main findings of the current study are the following:

•
The convergence of HAM solutions is ensured up to the 25th iteration. • Deterioration and elevation behaviour in the momentum of fluid is manifested in relation to Deborah's number and Reynold's number when β = −2.5 and β = 2.5.

•
The velocity shows continuous improvement with increases in the Hartman number in cases of dilatant and pseudoplastic fluid.

•
The growing estimate of the variable leads to increases in the dimensionless temperature field. • Declining aptitude expressed in the concentration field against the Schmidt number and Prandtl number. • A larger chemical reaction Cr portrays a decline in the concentration, while the Biot number Bi lead to the expansion in concentration.

•
The microorganism field has deteriorated for the higher value of Pe and microorganism difference parameter.