Modified Newton-EHS method for solving nonlinear problems with complex symmetric Jacobian matrices

1.   Introduction

Fluids can be categorized into Newtonian and non-Newtonian categories based on shear stress, according to the rheological behavior presented earlier ^[1]. The linear relationship between stress and strain identifies the Newtonian type. Physical problems involving Newtonian fluids have piqued the interest of scientists and engineers in recent years ^[2,3,4,5,6]. Unlike Newtonian fluids, the non-Newtonian fluid's viscosity does not remain constant with the shear rate and can be affected by a variety of factors. Slurries, foams, polymer melts, emulsions, and solutions are examples of non-Newtonian fluids. Non-Newtonian fluids have been extensively researched due of their numerous industrial applications. Furthermore, given the wide range of fluids of interest, heat transfer with non-Newtonian fluids is a vast topic that cannot be treated in its entirety in this study. Paint and adhesives industries, nuclear reactors, cooling systems, and drilling rigs are all examples of non-Newtonian fluid applications. To discover applications in new sectors, non-Newtonian fluid flows need a full investigation in terms of analytical, experimental, and numerical features. A variety of models have been used to explain non-Newtonian fluids in a number of works ^{[7,8,9,10,11,12,13,14,15,16,17,18,19,20]}. Non-Newtonian fluids were examined further by Megahed ^[7] and Ahmed and Iqba ^[8] for the power-law model; Cortell ^[9], Midya ^[10] and Megahed ^[11] for the viscoelastic model; Ibrahim and Hindebu ^[12], Bilal and Ashbar ^[13] and Abbas and Megahed ^[14] for the Powell-Eyring model; Pramanik ^[15], Rana et al. ^[16] and Alali and Megahed ^[17] for the Casson model; and Hayat et al. ^[18], Prasad et al. ^[19] and Megahed ^[20] for the Maxwell model. It is important to note that this study considers the viscous dissipation phenomenon. Only if the friction increases to the point where the fluid temperature field is noticeably warmed should this phenomenon be taken into account. There is a vast amount of information accessible on the viscous dissipation phenomenon and many fluid flow models ^[21,22,23].

The Cattaneo-Christov heat transfer model is a modified form of the Fourier law that is used to calculate the characteristics of a heat flux model while considering the relaxation time for heat flux dispersion across the physical model ^[24]. As a result, several researchers ^{[25,26,27,28,29,30]} have conducted more research on this model's most essential characteristics. In many industrial and technical sectors, including natural processes, biomedical engineering, chemical engineering and petroleum engineering, the Cattaneo-Christov heat transfer model makes the heat transfer mechanism and fluid flow in porous media crucially important ^[31]. So, the cooling process, which is one of the benefits, is one of the most important real life applications of our research, especially in light of the influence of the double diffusive Cattaneo-Christov model on the non-Newtonian Williamson model. Also, this benefit can lower the cost of the finished product and serve as the primary component in preventing product fault. In addition, the novelty of the current work is to examine how Cattaneo-Christov heat mass fluxes affect non-Newtonian dissipative Williamson fluids due to a slippery stretching sheet under the presence of thermal radiation, chemical reactions, and changing thermal characteristics.

2.   Basic governing equations

Consider the flow of a non-Newtonian Williamson fluid with the Cattaneo-Christov phenomenon in two dimensions. It is assumed that the stretched sheet, which is embedded in a porous medium with permeability $k$, is what causes the fluid to flow. Likewise, it is considered that the fluid concentration along the sheet is $C_{w}$ and that the sheet at $y = 0$ is heated with temperature $T_{w}$. The fluid is flowing in a streamlined pattern, and its viscosity is low but not negligible. The modifying impact appears to be restricted within a narrow layer adjacent to the sheet surface; this is called the boundary layer region. Within such a layer, the fluid velocity rapidly changes from its starting value to its mainstream value. We chose this model because it can accurately describe a large number of non-Newtonian fluids, possibly the majority, over a vast scope of shear rates. The x-axis is chosen parallel to the surface of the sheet in the flow direction, with the origin at the sheet's leading edge, and the y-axis is perpendicular to it (Figure 1). Let $T$ be the Williamson fluid temperature and $u$ and $v$ be the $x$-axial and perpendicular velocity components, respectively.

The sheet on which the slip velocity phenomenon occurs is supposed to be rough and exposed to thermal radiation with a radiative heat flux of $q_{r}$. The concept of the slip phenomenon was put forth by Mahmoud ^[32], who proved that it linearly depends on the shear stress $\tau_{w}$ at the surface, i.e.,

where $\tau_{w} = \left[\mu\frac{\partial u}{\partial y}+\mu\frac{\Gamma}{\sqrt{2}}\left(\frac{\partial u}{\partial y}\right)^{2}\right]_{y = 0}$ and $U_{w}$ is the fluid velocity at the sheet. Also, chemical reactions and the viscous dissipation phenomenon are taken into consideration in this study. Similarly, the temperature $T_{w}$ and concentration $C_{w}$ are constant on the sheet. We continue to assume that the Williamson fluid thermal conductivity varies linearly with temperature $\kappa = \kappa_{\infty}(1+\varepsilon \theta)$ ^[33], although its viscosity varies nonlinearly with temperature as $\mu = \mu_{\infty}e^{-\alpha\theta}$ ^[33], where $\mu_{\infty}$ is the Williamson fluid's viscosity at ambient temperature, $\alpha$ is the viscosity parameter, $\kappa_{\infty}$ is the liquid thermal conductivity at ambient temperature, and $\varepsilon$ is the thermal conductivity parameter. According to the aforementioned hypotheses, the fundamental equations governing the flow are

where $\underline {\rm{U}} = (u, v, 0)$, $p$ is the pressure, $\underline {\rm{r}} = \frac{-\mu}{k} \underline {\rm{U}}$ is the Darcy impedance for a Williamson fluid, $\Phi$ is the function of dissipation, $q$ is the heat flux, and $j$ is the mass flux. Also, the heat flux $q$ and the mass flux $j$ satisfy the following relations, respectively ^[34]:

Then, using the conventional boundary layer approximation, we can show that the equations below control the flow and heat mass transfer in our physical description ^[35].

In terms of boundary conditions, the fluid flow model is subjected to the following ^[36]:

where

In addition, $\rho_{\infty}$ represents the ambient density, and $\mu$ represents the Williamson viscosity. $\Gamma$ stands for a Williamson time constant, $k$ is the permeability of the porous medium, $c_{p}$ signifies the specific heat at constant pressure, $k_{1}$ is the reaction rate, $D$ is the diffusion coefficient, $\gamma_{1}$ is the time it takes for the heat flux to relax, $\gamma_{2}$ is the time it takes for the mass flux to relax, and $C_{\infty}$ is the ambient fluid concentration. Now, we utilize suitable similarity transforms, such as ^[35]:

When the last postulate is taken into account, Equation (2.8) is consequently confirmed, and the other equations are reduced to the following form:

The relevant physical boundary conditions are also modified as a result of the invocation of the prior dimensionless transformations as

The dimensionless controlling factors that have emerged are the local Weissenberg number $W_{e} = \Gamma x\sqrt{\frac{2a^{3}}{\nu_{\infty}}}$, the porous parameter $\Delta = \frac{\nu_{\infty}}{a k}$, the radiation parameter $R = \frac{16 \sigma^{*}T_{\infty}^{3}}{3k^{*}\kappa_{\infty}}$, the slip velocity parameter $\lambda = \lambda_{1}\sqrt{\frac{a}{\nu_{\infty}}}$, the thermal Deborah number $D e_{1} = \gamma_{1} a$, the Eckert number $Ec = \frac{u_{\mathrm{w}}^{2}}{c_{p}(T_{\mathrm{w}}-T_{\infty})}$, the Deborah number which related to the mass transfer field $D e_{2} = \gamma_{2} a$, the chemical reaction parameter $K = \frac{k_{1}}{a}$ and the Prandtl number $\operatorname{Pr} = \frac{\mu_{\infty} c_{p}}{\kappa_{\infty}}$. It is interesting to see that when $W_{e} = 0$, the nature of the non-Newtonian Williamson fluid model transforms into a Newtonian fluid model.

Furthermore, the local Nusselt number $Nu_{x}$, the local Sherwood number $Sh_{x}$, and the drag force coefficient in terms of $Cf_{x}$ are determined by ^[37]:

where $Re = \frac{u_{\mathrm{w}} x}{\nu_{\infty}}$ is the local Reynolds number.

3.   Outcomes with discussion

Table 1 compares the values of the local skin-friction coefficient ($Cf_{x}Re^{\frac{1}{2}}$) with those from Andersson's earlier work ^[38] to validate the current findings acquired by the shooting technique. This comparison is carried out for various slip velocity parameter values. We can confidently state that our results are in good accord with those references based on this comparison. Furthermore, the collected results show that the proposed strategy is both reliable and efficient.

The simulation studies of a non-Newtonian Williamson fluid over a stretching surface are presented in this section. The momentum field considers the slip velocity phenomenon, the energy equation includes the viscous dissipation and thermal radiation effects, and the mass transport equation includes the chemical process. Via the shooting approach, the governing equations are numerically solved using supported dimensionless transformation. All the governing physical parameters are utilized in the following ranges: $0.0\leq \Delta\leq0.5, 0.0\leq \alpha\leq0.5$, $0.0\leq \lambda\leq0.4, 0.0\leq W_{e}\leq0.5, 0.0\leq \varepsilon\leq0.5$,$0.0\leq De_{1}\leq0.4, 0.0\leq De_{2}\leq1.0$,$0.0\leq K\leq0.5$ and $0.0\leq Ec\leq0.5$. Thus, the values of the parameters with fixed values which are used for graphical display can be chosen as $\Delta = 0.2, W_{e} = 0.4, Sc = 4.0$, $R = 0.5, \alpha = 0.2, K = 0.2, \varepsilon = 0.2, \lambda = 0.2, De_{1} = 0.2, Pr = 7.0, De_{2} = 0.2$ and $Ec = 0.2$. This section discusses the graphical effects of physical dimensionless amounts on complicated profiles. The effect of porous parameter $\Delta$ on the Williamson fluid velocity is tested and introduced by means of Figure 2(a). When we assume $\Delta = 0.0, 0.2$ and $0.5$, both the velocity of the Williamson fluid and the thickness of the boundary layer are dramatically reduced. The porous parameter exhibits this behavior due to its direct presence in the velocity field. Figure 2(b) shows the impact of the porous parameter on both the fluid temperature and the fluid concentration. We discovered that the Williamson fluid temperature $\theta(\eta)$ rises as the porous parameter rises. Likewise, the Williamson fluid concentration $\phi(\eta)$ is subject to the same phenomenon. Physically, the porous parameter causes a restriction in the flow of the fluid, which reduces the fluid's velocity and raises the temperature and concentration distributions.

The velocity, temperature and concentration profiles for the viscosity parameter $\alpha$ are shown in Figure 3. It's worth noting that $\alpha$ has a greater impact on velocity profiles than on temperature or concentration profiles. Increases in the viscosity parameter $\alpha$ result in decreases in the velocity profile $f'(\eta)$ and associated boundary layer thickness, whereas the temperature $\theta(\eta)$ and concentration $\phi(\eta)$ fields show the opposite tendency. Physically, a barrier type of force will be produced in the Williamson flow due to the fluid viscosity's reliance on temperature. The fluid velocity can be slowed down by this force. The fluid layers consequently acquire little thermal energy via the same force.

Figure 4 shows the velocity, temperature, and concentration for the slip velocity parameter $\lambda$. When the slip velocity parameter is increased, both the temperature $\theta(\eta)$ and concentration $\phi(\eta)$ profiles increase, and as a result, the thermal boundary layer thickness increases. With the same parameter, however, the reverse tendency is observed for both the velocity distribution $f'(\eta)$ and the sheet velocity $f'(0)$. In terms of physics, the existence of the slip phenomenon produces a resistance force between the fluid layers, which in turn reduces the fluid's velocity and enhances the fluid temperature.

The influence of the local Weissenberg number $W_{e}$ on Williamson fluid velocity, Williamson fluid temperature, and Williamson fluid concentration is depicted in Figure 5. The graph illustrates that raising the local Weissenberg number $W_{e}$ causes temperature and concentration distributions to rise, whereas the same value of $W_{e}$ causes fluid velocity to decrease. Physically, a high local Weissenberg number increases the viscous forces that hold the Williamson fluid layers together, which lowers the fluid's velocity and improves the fluid's heat distribution through the boundary layer.

The impact of the thermal conductivity parameter $\varepsilon$ has been depicted in Figure 6. The larger values of $\varepsilon$ correspond to a broader temperature distribution $\theta(\eta)$ and a marginal increase in the concentration field $\phi(\eta)$, but an increase in the same parameter results in a slight decrease in the velocity graphs $f'(\eta)$. Physically, higher values of the thermal conductivity characteristic indicate a fluid's ability to gain higher temperatures, which may cause the fluid's temperature to rise through the thermal layer.

The effect of the thermal Deborah number $D e_{1}$ on velocity distribution, temperature distribution and concentration distribution is depicted in Figure 7. The Deborah number is a rheological term that describes the fluidity of materials under specified flow circumstances. Low-relaxation-time materials flow freely and exhibit quick stress decay as a result. The fluid motion $f'(\eta)$ scarcely improves with the bigger $D e_{1}$ parameter, but the concentration profiles $\phi(\eta)$ are slightly reduced. Furthermore, extended values of the thermal Deborah number $D e_{1}$ reduce both the thickness of the thermal region and the profiles of temperature $\theta(\eta)$.

Figure 8 depicts the dimensionless velocity, dimensionless concentration, and thermal profiles for different Eckert number $Ec$ estimates. It is evident that as the Eckert number is increased, the dimensionless velocities $f'(\eta)$ diminish modestly, while the dimensionless concentration $\phi(\eta)$ increases slightly. Furthermore, the increase in Eckert number $Ec$ obviously increases both the temperature profile and the thickness of the thermal region. Physically, the Williamson fluid moves quickly, causing some kinetic energy to be transformed into thermal energy as a result of the viscous dissipation phenomenon. This promotes the fluid distributing heat more through the thermal layer.

Figure 9 shows the effects of the Deborah number $D e_{2}$, which is related to mass transfer, and the chemical reaction parameter $K$ on Williamson fluid concentration. Because both the concentration characteristics $\phi(\eta)$ and the concentration boundary thickness of the Williamson fluid decrease when both the chemical reaction parameter and the Deborah number increase, the mass transfer rate increases as well.

From an engineering standpoint, we now concentrate on the fluctuations of physical quantities of interest. For all regulating factors of our model, the local skin-friction $Cf_{x}(Re_{x})^{\frac{1}{2}}$, the local Nusselt number $\frac{Nu_{x}Re_{x}^{\frac{-1}{2}}}{1+R}$, and the local Sherwood number $Sh_{x}Re_{x}^{\frac{-1}{2}}$ are introduced in Table 2. Skin friction coefficient values decrease when the viscosity parameter, slip velocity parameter, and local Weissenberg number increase, lowering both the local Nusselt number and the local Sherwood number. The skin friction coefficient and the local Nusselt number, respectively, grow and shrink uniformly with the thermal Deborah number and the Eckert number. Furthermore, increasing the chemical reaction parameter or the Deborah number increases the rate of mass transfer, while increasing the thermal conductivity parameter reduces it significantly.

4.   Conclusions

A description of the heat and mass transport features of a viscous non-Newtonian Williamson fluid across a stretching sheet that is embedded in a porous medium was attempted here. Heat and mass transfer in the Williamson fluid in the presence of thermal radiation, slip velocity, variable thermal characteristics, Cattaneo-Christov heat mass fluxes, chemical reaction and viscous dissipation are all investigated in this study. The governing equations have been converted via dimensionless transformations into a set of coupled nonlinear ordinary differential equations, which are then numerically solved using the shooting technique. The following are the major characteristics of the current work.

$1)$ A concentration field grows as the viscosity and slip parameters increase, but as the chemical reactions and Deborah number increase, the concentration field decreases.

$2)$ The Sherwood number decreases as the Eckert number rises, whereas the Deborah number and chemical parameter increase it.

$3)$ The velocity profile is reduced by the viscosity and slip velocity parameters, whereas a minor increase in velocity profile occurs due to the Deborah number being increased.

$4)$ The Deborah number has a lower temperature profile than the thermal conductivity parameter and the Eckert number.

$5)$ The heat transfer rate is reduced by both the thermal conductivity parameter and the viscous dissipation phenomenon.

$6)$ Skin-friction coefficient values increase with the porous parameter, but they decrease with an increase in either the slip velocity or the viscosity parameter.

$7)$ When the local Weissenberg number is higher, both the skin-friction coefficient and the local Nusselt number's magnitude fall.

Acknowledgments

The authors are thankful to the honorable editor and anonymous reviewers for their useful suggestions and comments which improved the quality of this paper.

Conflict of interest

The authors declare that they have no competing interests.