On Criticality for a Generalized Couette Flow of a Branch-Chain Thermal Reactive Third-Grade Fluid with Reynold's Viscosity Model

This research considers the third-grade liquid flow and criticality branched-chain of a thermal reaction in a Couette generalized medium with a nonlinear viscosity model. A dimensionless transformation of the system momentum and heat equations are carried out. Compared with the diffusion coefficient, the flow is stimulated by initiation reaction rate, reaction branch-chain order, non-Newtonian term, thermal Grashof number, and pressure gradient. The reactive fluid is fully exothermic with consumption of the material, and the heat exchange in the system is greater than the exchange of heat with the ambient. A semianalytical collocation weighted residual scheme is employed for the branch-chain slice bifurcation, dimensionless energy, and momentum solutions. The results show that exponential decreases in the thermal fluid viscosity can help in controlling the boundless heat produced by the Frank-Kamenetskii term and initiation reaction rate. Therefore, the results will help in stimulating positive combustion processes.


Introduction
Limited industrial usage of Newtonian fluid has substantially encouraged the interest in non-Newtonian fluid and its applications in science, technology, and manufacturing processes such as polymer film production and fiberglass. e flow of liquids under the action pressure gradient in a device with the moving part where lubrication takes place described Couette generalized flow [1,2]. Non-Newtonian fluid with viscoelastic characteristics can help in improving lubricants of industrial machines and technology devices' efficiency. As such, Nayak et al. [3] examined entropy optimization of the hydromagnetic non-Newtonian nanomaterial with a joint approach for heat transfer intensification and solar energy absorber. e problem was solved numerically, and the result shows that a rise in the Weissenberg number enhanced fluid velocity. Wang et al. [4] reported on the non-Newtonian reaction catalyst of the heterogeneous-homogeneous Oldroyd fluid with radiation and heat absorption. From the study, it was revealed that radiation and Biot number have a direct relation to the temperature distribution. With heat transfer, it is known that shear rate and friction between the fluid and moving devices can lead to significant heat production in the system that may affect the device and fluid properties [5][6][7][8]. Among the liquids that exhibit viscoelastic characteristics is the third-grade liquid. ird-grade fluid model is a non-Newtonian liquid model that predicts shear thickening or thinning properties over a definite boundary.
Due to its applications, Yilbas et al. [9] examined entropy production of the third-grade liquid with unvarying viscosity in an annular pipe using analytical approximate solutions for the entropy generation and other flow properties. It was reported that improved non-Newtonian material terms will assist in reducing the irreversibility process. Yurusoy et al. [10] improved on the work of Yilbas et al. [9] by considering Vogel model viscosity and solved the steady third-grade fluid flow in a concentric rigid cylinder by the perturbation method. It was noticed that heat-dependent viscosity supported increasing the velocity field. In a convective cooling channel, the irreversibility of the flow of the third-grade reactive liquid with dissipative viscous heating was carried out by Salawu and Adesanya et al. [11,12]. e solutions to the equations were obtained using semianalytical techniques, and it was reported that the non-Newtonian term decreases the flow rate. Okoya [13] studied computationally the nonlinear viscosity of the third-grade reactive fluid with heat effect on axial annular flow. It was revealed that the reactive third-grade and Newtonian fluids are the same qualitatively in an annulus pipe except for the level of inward translocation of the peak temperature and axial momentum. Khan et al. [14] considered modeling and simulation of the micropolar ferrofluid for slip velocity of second order with saturated permeable media. It was found that the fluid material terms decrease the fluid velocity due to simulation in the liquid viscosity. However, in a reactive combustion heat transfer process without heat diffusion, exothermic reaction will reach the state of criticality no matter the starting temperature [15,16]. ermal criticality performs a significant role in processing and handling of the non-Newtonian fluid. It exists when the heat generation rate in a reactive flow system transcends heat dissipation to the environments [17,18].
is occurrence is the basis for thermal ignition or runaway in a flow system, Frank-Kamenetskii [19]. e main objective of analyzing thermal criticality is to predict the unsafe or critical state of an exothermic reaction flow condition in a combustion process. Combustion is absolute and essential for the chemically reacting flow system with applications in pollution control, power generation, processing material industries, and so on, Balakrishnan et al. [20]. As a result of its importance, Salawu et al. [21] investigated the irreversibility and thermal criticality of Powell-Eyring reactive liquid flow with radiation and variable conductivity in permeable media. e authors revealed that thermal explosion can be avoided in a reactive system if heat source terms are minimized. Makinde and Maserumule [5] examined thermal ignition and the second law of Couette fluid flow with variable viscosity. Analytical solutions of the problem were provided using perturbation techniques along with the Hermite approximation method. e results obtained proved clearly that the Frank-Kamenetskii term is a strong internal heat generation term that leads to thermal criticality or ignition. In a porous-filled channel, analytical solution of viscous reactive thermal ignition in a slab was carried out by Makinde [22]. In the study, essential properties of heat distribution and criticality slice bifurcation were reported. Also, with series of solution technique, thermal explosion of exothermic reaction in a slab was examined by Makinde [23]. e study presented that energy reactive parameters can assist in controlling exothermic reaction combustion explosion.
is study is built based on Okoya [16] which is an extension of the work done by Varatharajan and Williams [24] in which branch-chain transition and criticality in a slab are considered. e indisputable results obtained in their studies and suggested further extension have motivated this investigation. is study presents coupling fluid flow and heat transfer to investigate the bifurcation slice for the thermal criticality of the third-grade reactive Couette fluid with nonlinear viscosity and heat distribution. erefore, it is necessary to study the criticality of the branched-chain thermal-reactive diffusion problem in order to obtain competing effects of emerging parameters and enhance its utilization in different manufacturing industries. is study is essential in determining combustion system capability under different operating conditions in order to improve safety; it will also help in decreasing combustion product pollution.
e considered problem is solved by weighted residual techniques coupled with the collocation scheme. Results from the method used are observed to quantitatively and qualitatively agree with the computed results.

Problem Formulation
Consider the flow of the reactive mixture in a Couette generalized device ( Figure 1) bounded with isothermal walls. e flow has a unidirectional velocity field with Reynold's viscosity under the influence of gravity, pre-exponential factor m, and n order of reaction branching. e isothermal flow configuration is positioned at y � [−h; h] with y-axis normal to the flow direction in the x-axis, and the non-Newtonian formulation is used to create the viscoelastic effects. e branched-chain parameter is described in the form of the reaction rate generalized law as e consumption chemical reactant is taken to be very small with constant coefficients. e terms P 0 , Z, k, E, m, n, ], T 0 , and T correspond to the branch-chain rate, Planck's number, Boltzmann constant, activation energy, preexponential constant, reaction branch order, frequency of vibration, initial heat, and fluid heat. e momentum singlefluid balance for the dynamic species is given as follows (Truesdell [25]): Given that the chemical reaction does not create flow momentum, the kinematic-related variables for the thirdgrade fluid that are thermodynamically appropriate in the stress tensor form according to Fosdick and Rajagopal [26] can be expressed as where μ is the dynamic viscosity, I is the unit tensor, p is the pressure, and tr is the trace matrix. e terms r 1 , c 3 , and r 2 are the variable temperature of the material coefficients defined as follows: e kinematic tensors B 1 and B 2 can be expressed as e term ∇ is an operator (gradient), index T denotes transpose, v is the vector velocity, and the derivative material term is d/dt. e flow ensues in an erect device with a spontaneous velocity U as the chemical reaction is being prompted to raise proliferation of chain carriers. e reactive fluid accelerates in the medium close to the wall, while the other wall remains static. Following the assumptions stated above and [16], the nondimensional steady velocity and energy balance equations are written as Here, the parameters w and θ depict nondimensionless flow velocity and energy. e terms Q, m, ϵ, n, λ, Gr, G, a, and Λ individually represent the initiation rate, pre-exponential factor, activation energy (that is, the energy needed for a species reaction mixture to take place; thus, it determines the transition state of a chemical reaction), branch-chain order, Frank-Kamenetskii (that is, homogeneous mixture of thermal ignition species at walls' constant temperature. It determines the species reaction time to heat conducting time), heat Grashof number, pressure gradient, variable viscosity, and non-Newtonian.
e suitable nondimensional boundary conditions are e Scientific World Journal e heat-dependent dynamic viscosity is assumed to exponentially vary according to Salawu et al. [17,27]: e succeeding quantities are utilized to get nondimensional formulated models (6) to (8):    (6) and (7) correspond to the case of the nonlinear viscosity flow model with heat transfer investigated by Okoya [28], while the corresponding model with Gr ≠ 0 was considered by Salawu and Fatunmbi [11,29]. In the absence of momentum and reaction branch order, the exceptional instance of ϵ � Q � 0 corresponds to the Frank-Kamenetskii traditional model [19], while the case ϵ ≠ 0 and m ∈ 0.5, 0, −2 { } was described was described by Boddington et al. [30].

Method of the Solution
e solution technique for the criticality bifurcation slice, temperature, and momentum equations is performed using a semianalytical technique. In the scheme, it is taken that as defined in [31,32], and ϕ i (y) is an assigned function with the boundary conditions satisfied. e basis function u(y, d) is individually illustrated for the boundary conditions and their respective equations. Hence, for the arbitrary chosen values of d ′ s, the residual equation is obtained as For the function ϕ i in successful approximation, the differential equations are satisfied by u(y, d). e main aim is to minimize W(y, d) errors to zero (say) in between the domain, i.e., e number of weighted functions V i must be translated to the number of unknown constants d i ′ s in u. e integration collocation scheme is adopted for the solution in which the weight function is presented in Dirac delta as e method is employed on boundary conditions (8) and on dimensionless equations (6) and (7) to have the following residual equations:   w r � G + exp −a y 10 b 10 + y 9 b 9 + y 8 b 8 + y 7 b 7 + y 6 b 6 + y 5 b 5 + y 4 b 4 + y 3 b 3 + y 2 b 2 + yb 1 + b 0 × 90y 8 a 10 + 72y 7 a 9 + 56y 6 a 8 + 42y 5 a 7 + 30y 4 a 6 + 20y 3 a 5 + 12y 2 a 4 + 6ya 3 + 2a 2 − a × exp −a y 10 b 10 + y 9 b 9 + y 8 b 8 + y 7 b 7 + y 6 b 6 + y 5 b 5 + y 4 b 4 + y 3 b 3 + y 2 b 2 + yb 1 + b 0 × 10y 9 a 10 + 9y 8 a 9 + 8y 7 a 8 + 7y 6 a 7 + 6y 5 a 6 + 5y 4 a 5 + 4y 3 a 4 + 3y 2 a 3 + 2ya 2 + . . .

Discussion of Results
e solutions to the thermal criticality, velocity, and heat equations are obtained via the weighted residual collocation scheme.
e computational default values are taken according to [33,34]. e results obtained are presented as plots in Figures 2-11. e accuracy and consistency of the method are established by comparing it with the

e Scientific World Journal
Fehlberg-Runge-Kutta scheme along with the shooting method (numerical method) as depicted in the table. In Table 1, the results for the weighted residual technique and numerical scheme are compared as presented. e semianalytical method used gives good solutions when relating with the numerical method results. e collocation WRM contracted well with the other solution approaches with the order of absolute error 10 − 8 . Figure 2 illustrates the impact of Reynold's viscosity term a on the velocity profile. As the values of the heat-dependent viscosity term increase, the fluid viscosity diminishes, and the fluid bonding force is discouraged. As such, the fluid molecular diffusion is boosted which correspondingly diminishes the flow resistance force. Hence, the non-Newtonian fluid particle collision is enhanced that leads to an increasing fluid velocity distribution. In Figure 3, the reaction of the non-Newtonian flow liquid to rising in the Grashof number Gr is demonstrated. e term Gr is the ratio of buoyancy to the viscous force exerting on a liquid which is equivalent to the Reynolds number. A rise in the thermal Grashof number enhances heat source terms which leads to a breakdown in the fluid bonding forces. e breaking down in the liquid bonding forces causes the fluid particles to move freely and thereby increases the flow rate. Figure 4 shows the response of the third-grade liquid to variation in the non-Newtonian term Λ. Large shear rate occurs in a generalized Couette flow, which acts as a source term for the heat and velocity equations. In a third-grade liquid, the term Λ needs to be kept at small values to discourage flow opposition forces. erefore, increasing the values of the term Λ decreases the flow velocity field due to a decrease in the source terms. Figures 5-7 depict the temperature distribution for the rising values of the terms n, Q, and λ for the function θ(y) plotted against y. e impacts of branch reaction order n on the temperature profile are displayed in Figure 5. e branch reaction order is a basis for thermal ignition; it entails steps of the reaction chain in which each step serves as a reagent for the next step. erefore, rising the reaction branch order reduces the heat profile. is is as a result of heat diffusion from the device that causes general decreases in the system heat field. e influence of the rate of initiation Q and Frank-Kamenetskii term λ on the energy distribution is illustrated in Figures 6 and 7. Both parameters are high heat generation terms and must be carefully monitored to avert the blow up

Criticality Branch-Chain and Ignition
Solutions. e slice bifurcation thermal runaway plot for 0 < ϵ ≪ 1 and different pre-exponential factor m in the plane (λ, θ max ) is demonstrated in Figures 8 and 9. e diagrams depict the qualitative difference in the reactive non-Newtonian liquid as the term λ increases. In fact, for 0 ≤ ϵ ≪ 1 and m � −2, 0, 0.5, critical value λ cr occurs so that, for 0 < λ < λ cr , two solution branches are found. e lower solution branch is steady, while the upper branch solution diverges to infinity as λ ⟶ 0. When λ > λ cr , there is no real solution but a traditional form demonstrating temperature criticality. A rise in the activation energy ϵ raises the system thermal explosion but declines the system heat steadiness as the chemical kinetics m is encouraged as described in Figures 8  and 9. e considerable variations in the maximum heat (θ max ,cr ) and ignition criticality (λ cr ) for a third-grade reactive flow as the parameter ϵ rises are confirmed in Figures 10 and 11. An increase in the branching reaction order n causes complete reduction in a reactive exothermic diffusion flow system as seen in the figures.    branch-chain criticality bifurcation, momentum, and temperature field solutions are analytically obtained. It is revealed from the investigation that (i) e activation energy has no influence on the considered non-Newtonian fluid and temperature distribution in a Couette flow device (ii) It is obtained that, for engineering processes, the non-Newtonian reactive third-grade term enhances fluid viscosity by decreasing the flow rate (iii) From the aspect of industrial design, boundless heat production in the system due to the Frank-Kamenetskii term and initiation reaction rate can be controlled by exponential decrease in the fluid viscosity e results will assist in combustion processes and in improving safety conditions of thermal explosion system. e solution algorithm is an interesting tool to study different parameter-dependent boundary-value nonlinear problems in engineering and science.

Data Availability
No data were used to support this study.

Conflicts of Interest
e authors declare that they have no conflicts of interest.