Abstract
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.
1. 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. The 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. The 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–8]. Among the liquids that exhibit viscoelastic characteristics is the third-grade liquid. Third-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]. The 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].
Thermal 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]. This occurrence is the basis for thermal ignition or runaway in a flow system, Frank-Kamenetskii [19]. The 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. The 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. The 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]. The study presented that energy reactive parameters can assist in controlling exothermic reaction combustion explosion.
This 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. The indisputable results obtained in their studies and suggested further extension have motivated this investigation. This 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. Therefore, 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. This 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. The 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.
2. Problem Formulation
Consider the flow of the reactive mixture in a Couette generalized device (Figure 1) bounded with isothermal walls. The flow has a unidirectional velocity field with Reynold’s viscosity under the influence of gravity, pre-exponential factor , and order of reaction branching. The isothermal flow configuration is positioned at with y-axis normal to the flow direction in the x-axis, and the non-Newtonian formulation is used to create the viscoelastic effects. The branched-chain parameter is described in the form of the reaction rate generalized law as
The consumption chemical reactant is taken to be very small with constant coefficients. The terms , , , , , , , , and correspond to the branch-chain rate, Planck’s number, Boltzmann constant, activation energy, pre-exponential constant, reaction branch order, frequency of vibration, initial heat, and fluid heat. The momentum single-fluid 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 third-grade fluid that are thermodynamically appropriate in the stress tensor form according to Fosdick and Rajagopal [26] can be expressed aswhere is the dynamic viscosity, is the unit tensor, is the pressure, and is the trace matrix. The terms , , and are the variable temperature of the material coefficients defined as follows:
The kinematic tensors and can be expressed as
The term is an operator (gradient), index denotes transpose, is the vector velocity, and the derivative material term is . The flow ensues in an erect device with a spontaneous velocity as the chemical reaction is being prompted to raise proliferation of chain carriers. The 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 and depict nondimensionless flow velocity and energy. The terms , , , , , , , , 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. The suitable nondimensional boundary conditions are
The heat-dependent dynamic viscosity is assumed to exponentially vary according to Salawu et al. [17, 27]:
The succeeding quantities are utilized to get nondimensional formulated models (6) to (8):
2.1. Limiting Cases
When , equations (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 was considered by Salawu and Fatunmbi [11, 29]. In the absence of momentum and reaction branch order, the exceptional instance of corresponds to the Frank-Kamenetskii traditional model [19], while the case and was described was described by Boddington et al. [30].
3. Method of the Solution
The solution technique for the criticality bifurcation slice, temperature, and momentum equations is performed using a semianalytical technique. In the scheme, it is taken thatas defined in [31, 32], and is an assigned function with the boundary conditions satisfied. The basis function is individually illustrated for the boundary conditions and their respective equations. Hence, for the arbitrary chosen values of , the residual equation is obtained as
For the function in successful approximation, the differential equations are satisfied by . The main aim is to minimize errors to zero (say) in between the domain, i.e.,
The number of weighted functions must be translated to the number of unknown constants in . The integration collocation scheme is adopted for the solution in which the weight function is presented in Dirac delta as such that .
The method is employed on boundary conditions (8) and on dimensionless equations (6) and (7) to have the following residual equations:
Within the domain, the collocation techniques are applied on equations (14) and (15) which are then solved together with the boundary conditions to get the coefficients and . The solution algorithms are repeated for various values of parameters. Maple software is used to determine the constant coefficients and solve the equations completely.
4. Discussion of Results
The solutions to the thermal criticality, velocity, and heat equations are obtained via the weighted residual collocation scheme. The computational default values are taken according to [33, 34]. The results obtained are presented as plots in Figures 2–11. The accuracy and consistency of the method are established by comparing it with the 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. The semianalytical method used gives good solutions when relating with the numerical method results. The collocation WRM contracted well with the other solution approaches with the order of absolute error .
4.1. Velocity Profile for Parameter-Dependent Solutions
Figure 2 illustrates the impact of Reynold’s viscosity term 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 is demonstrated. The term 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. The 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. Therefore, increasing the values of the term decreases the flow velocity field due to a decrease in the source terms.
4.2. Temperature Field for Parameter-Dependent Solutions
Figures 5–7 depict the temperature distribution for the rising values of the terms , , and for the function plotted against . The impacts of branch reaction order on the temperature profile are displayed in Figure 5. The 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. Therefore, rising the reaction branch order reduces the heat profile. This is as a result of heat diffusion from the device that causes general decreases in the system heat field. The influence of the rate of initiation 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 of the reactive solution. The temperature profile is enhanced at the fixed wall and consistently rises until it gets to the highest at the channel center. It then diminishes continuously close to the moving surface. Continuous decrease in the heat profiles is observed as the viscosity of the liquid reduces in a reactive exothermic third-grade liquid. From the viewpoint of technology design, total heat generation can be minimized if the viscosity of the liquid as a result of rises in heat can be exponentially reduced.
4.3. Criticality Branch-Chain and Ignition Solutions
The slice bifurcation thermal runaway plot for and different pre-exponential factor in the plane is demonstrated in Figures 8 and 9. The diagrams depict the qualitative difference in the reactive non-Newtonian liquid as the term increases. In fact, for and , critical value occurs so that, for , two solution branches are found. The lower solution branch is steady, while the upper branch solution diverges to infinity as . When , 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 is encouraged as described in Figures 8 and 9. The considerable variations in the maximum heat and ignition criticality for a third-grade reactive flow as the parameter rises are confirmed in Figures 10 and 11. An increase in the branching reaction order causes complete reduction in a reactive exothermic diffusion flow system as seen in the figures.
5. Conclusion
The criticality of a reactive generalized third-grade Couette fluid flow with Reynold’s viscosity model is analyzed using the collocation weighted residual scheme. The thermal branch-chain criticality bifurcation, momentum, and temperature field solutions are analytically obtained. It is revealed from the investigation that(i)The 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
The results will assist in combustion processes and in improving safety conditions of thermal explosion system. The 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
The authors declare that they have no conflicts of interest.