Neural networking analysis for MHD mixed convection Casson flow past a multiple surfaces: A numerical solution

1.   Introduction

Non-Newtonian fluids are used in a variety of products, including culinary items, personal protective equipment, printing technologies, damping and braking systems, and drag-reducing agents ^[1]. Fluids that are not Newtonian exhibit complex rheological behavior ^[2]. Unlike Newtonian fluids, these materials do not behave linearly between the shear rate or shear and viscosity. Non-Newtonian fluids ^[3,4] are categorized into many groups based on how they react to changes in shear stress or rate. Based on their behavior, non-Newtonian fluids can be classified as shear thinning ^[5,6], shear thickening ^[7,8], thixotropic, visco-plastic, and rheopectic fluids. Following such characteristics, various studies were given by researchers to explore the type of non-Newtonian fluid models and the corresponding daily life applications. In this direction, the Casson fluid model is one of the non-Newtonian fluid models that best capture yield stress and has noteworthy applications in the biomechanics and polymer sectors. Several different goods, including synthetic lubricants, medicinal compounds, paints, coal, china clay, and many others, can be prepared using Casson fluid. Human blood can also be referred to as Casson fluid since it contains a variety of components including human red blood cells, fibrinogen, protein, and globulin in aqueous base plasma. The operative use of Casson fluid in biological treatments, drilling processes, bio-engineering and food processing has captured the attention of researchers like the analysis of diffusion phenomena in stenosed capillary-tissue exchange systems provided by Siddiqui and Mishra ^[9]. In order to depict the blood flow they utilize a modified Casson's fluid description. Through the alteration of a measure known as the retention parameter, the severity of the condition could be evaluated. Another significant aspect of this analysis was the increase of these variables together with the progression of the stenosis. For the investigation of the normal and sick states, the concentration profiles were identified. In the recent past, according to Abolbashari et al. ^[10], engineering and industrial applications have found usage for nano-fluids because of their fantastic thermal conductivity increase. Due to their significance, analyze the fluid non-Newtonian flow with nanoparticles towards the sheet. In special circumstances, the current analytical solution shows a very strong correlation with those of the earlier research that was published. The impact of various flow physical factors were thoroughly explored. Bhattacharyya et al. ^[11] conducted an examination of Casson fluid flow past a shrinking surface with different effects. Self-similar nonlinear ordinary differential equations (ODEs) were obtained. The translated energy and velocity equations both had two exact answers. According to the depicted data, the temperature drops for greater values of radiation parameter, Prandtl number, and power-law exponent. Additionally, the thickness of the thermal boundary layer increases with direct variation in wall temperature. In some instances, temperature overshoot was seen in the temperature field's graphical depiction. Therefore, under some circumstances, heat transfer from the surface happens instead of heat absorption at the surface. In the presence of fluctuating suction, Abd El-Aziz and Yahya ^[12] investigated the impact of Hall current on the Casson fluid flow over a sheet. To obtain non-dimensional flow fields, the resulting equations were analytically solved. The temperature and velocity profiles are studied for various parameter values that entered the problem and were displayed graphically. Physically relevant Nusselt number was sorted out and displayed as graphics. For an axisymmetric Casson stagnation point, Nawaz et al. ^[13] looked at the effects of Joule heating and viscous dissipation. The ensuing boundary value issue was resolved using the homotopy analysis method. The homotopy solutions and the numerical outcomes were compared. The behavior of various parameters was examined. The Casson fluid parameter's resultant effect seems to accelerate the fluid's velocity. The actions of dissipation and Joule heating cause the system to warm up. The thermal and velocity slip effects on a Magnetohydrodynamics (MHD) Casson fluid were studied by Usman et al. ^[14] over a cylinder with Buongiorno's model. The Casson nanofluid model's equations were constructed and reduced by using similarity transformation. The collocation approach was used to arrive at the numerical solution. The graphs that compare the nanofluid fluid flow to the emerging physical parameters provide an analysis of the physical quantities of interest. The effects of Newtonian heating on the Casson fluid flow were examined by Tassaddiq et al. ^[15]. In this study, porous effects for such fluids and MHD effects were also taken into account. Equations involving partial differentials were used to model the fundamental issue. The analytical tool is known as the Laplace transform was used to achieve the "Velocity" and "Temperature" functions. In order to analyze the modeling parameters that were used, graphical representations were used. To verify the facts, numerical calculations were made. The graphical results show that velocity clearly declines for magnetic parameter is intensified and increases as the porosity parameter is increased (conjugate parameter). For all conceivable Casson parameter values, the fluid flow might be controlled. A comparative examination of multiphase flow over a steep channel was carried out by Hussain et al. ^[16]. The combination of microscopic gold particles with Casson fluid creates an intuitive non-Newtonian particulate suspension. Consideration of a uniformly slanted conduit allowed for the consideration of gravitational and magnetic influences. The shear-thinning phenomena is further influenced by heating effects at the channel's edge. For two-phase fluid flow with heat transfer, an exaction was found. In addition, Newtonian-gold particulate flow and Casson-gold particulate flow, which was a non-Newtonian suspension of particle flow, are contrasted. It was discovered that the shear-thinning effects of magnetized Newtonian particle suspension across the inclined channel are more pronounced and experience reduced skin friction. However, it was discovered that Casson fluid was a very helpful suspension for the fabrication and coatings processes, among other things. Additionally, the magnetic field affects the mobility of both phases as a resistive force. Due to the system's significant viscous dissipation, more energy was added. Akolade and Tijani ^[17] discussed the effect of various transport parameters on Casson- nanofluids through a Riga surface. After imposing sufficient similarity factors, the flow equations were derived. To account for all relevant flow characteristics, the spectral quasi-linearization technique (SQLM) was used. Aside from great agreement with previously published work, the approach performance (numerical estimation tools) was examined in terms of CPU time and mistakes. For the current challenge, engineering dimensionless parameters were also supplied. The outcomes demonstrated that the Casson fluid had a lower flow resistance than the Williamson fluid. Furthermore, Casson fluid diffuses and conducts less than Williamson fluid. Siddiqui et al. ^[18] investigated the Casson nanofluid between two spinning and stretchy disks in this study. The effects of thermal radiation, magnetic field, Darcy Forchheimer flow, and heat source were investigated. The shooting strategy was owned to solve the existing equations with the help of MATLAB packages. The temperature profile when the Prandtl number was varied was also discussed. Furthermore, the thermal profile towards the thermophoresis, radiation, and Brownian parameter was investigated. Graphs were used to evaluate the impact of nanoparticle concentration on microorganisms with rising Peclet and Lewis numbers. References ^{[19,20,21,22]} provide an evaluation of the most recent attempts at Casson fluid flow in a variety of configurations.

In this paper, we used the prediction application of artificial intelligence to examine the thermal stagnation point magnetized flow field of non-Newtonian fluid flow due to inclined heated stretched surfaces with thermos-physical effects namely viscous dissipation, mixed convection, heat generation, and variable thermal conductivity. The flow is mathematically formulated for a cylindrical surface and later mathematical reduction is done for thermal flow over a flat surface. The ultimate results are obtained by using the shooting method and are offered as line graphs. The quantity of interest includes the temperature and NN at both flat and cylindrical surfaces. Two different ANN models are constructed to forecast the heat transfer coefficient at surfaces. We are confident that the thermal findings by the use of ANN models will helpful for investigators affiliated with the area of thermal engineering.

2.   Flow formulation

Stretching inclined surfaces, such as the cylinder and plate, are used to introduce non-Newtonian fluid flow. By taking into account the magnetic field ^[23,24,25], mixed convection, stagnation point flow ^[26], temperature-dependent conductivity ^[27], viscous dissipation ^[28], and heat generation ^[29] the novelty of the flow field is strengthened. By utilizing the concentration equation, the mass transfer is taken into account. It is assumed that the strength of temperature and concentration is considered higher at surfaces in comparison with strength far away from the surface. The concluding Casson fluid flow equations ^[26] for the given problem are expressed as

The radioactive flux and thermal conductivity relations are given as:

The flow conditions are:

Here, $\nu$ is kinematic viscosity, ${R}_{c}$ stands for radius of cylinder, ${B}_{T}$ is the thermal expansion coefficient, ${B}_{C}$ denote the coefficient of solutal expansion, $\rho, \widehat{T}, \widehat{C}, {Q}_{0}, \varepsilon , {\kappa }_{\infty }~~ \text{and}~~ {D}_{m}$ are fluid density, temperature, concentration, heat generation coefficient, small parameter, heat conductivity far from the surface and mass diffusivity respectively. By using NN-ANN model, we are interested in comparing how flow parameters affect the Nusselt number and temperature of Casson fluid flow toward both surfaces. Consequently, we must convert these partial differential equations (PDEs) into corresponding ODEs in order to solve Eq (1) through (4). For this reason, we have own the following set of variables ^[26,27]:

Owning Eq (7), Eqs (1)-(3) turns into

and flow conditions reduced as:

The parameters are given as:

The Casson fluid parameter, Prandtl number, Eckert number, thermal radiation parameter, Solutal Grashof number, thermal Grashof number, velocities ratio parameter, and Schmidt number are symbolized by ${\beta }_{p}, {Pr}, {E}_{n}, {R}_{p}, {G}_{C}, {G}_{T}, {A}_{p}\text{and}Sc$. The heat transfer at inclined surfaces is considered a key physical effect and to evaluate heat transfer we considered the Nusselt number. The mathematical relationship ^[26] of NN is:

3.   Numerical scheme

The lower order equations, Eqs (8)-(10), can be used to examine the non-Newtonian fluid flow regime. Due to their coupling and nonlinearity, they are challenging to solve analytically. The flow differential equations are reported using a number of various solution approaches ^{[30,31,32,33,34,35]}, however, the shooting method will be employed to get numerical solution for the present research problem. The name of the shooting technique is taken from an analogy with target shooting: we shoot at the target and watch where it lands; based on the misses, we can modify our aim and shoot again in the hopes that it would hit close to the target. To convert boundary value problems into similar initial value problems, the shooting method was developed. The supportive material in this regard can be accessed in Refs. ^[25,26,27]. Having an initial value system is a vital first step in this regard. We can achieve it by presuming:

under Eq (14), Eqs (8)-(10) reduces to

The flow equations turns down to:

The self-coding scheme of shooting is carried out in Matlab for the present problem and outcomes are shared in terms of graphs and numerical data.

4.   Artificial neural networking model

Mathematical modeling is done for the Casson flow over heated stretched surfaces. The equations for non-Newtonian flow are numerically solved. Both a flat plate and a cylindrical surface are used to evaluate the Nusselt number. The Nusselt number at surfaces is predicted using ANN models in both situations. Researchers firmly think that the multilayer perceptron (MLP) and ANN model ^{[36,37,38,39,40]} have excellent learning capabilities and can be utilized to anticipate a variety of physical occurrences. In MLP networks, three separate layers are utilized. The first layer uses the inputs, while the centrally important layer is referred to as the hidden layer. The output layer, which is the final layer, is where the prediction data is stored. The NN-ANN-I and II are two ANN models that we have created. Whereas NN-ANN-II is designed to provide predictions of Nusselt number values at cylindrical surfaces, NN-ANN-I is built to forecast the Nusselt number at flat surfaces. The five parameters are employed as inputs and the NN is taken into account as an output for both NN-ANN models. Variable thermal conductivity, heat generation, Prandtl number, the Casson fluid, and the thermal radiation parameters are the flow parameters.

Table 1 provides the symbolic information for models. Network architecture is given in Figure 1. We have gathered 80 samples for the five inputs, leading to 80 sample values for the Nusselt number. 15% of the data is used for validation while the remaining 70% is used to train the network. Testing also makes advantage of the 15%. Table 2 provides all available details in this regard. Ten (10) neurons are taken into account in the hidden layer, and the network is trained using the Levenberg-Marquardt method. The transfer functions used in hidden and output have the following expressions:

We have taken into consideration the mean squared error (MSE) and coefficient of determination (R) for performance analysis of ANN models. The mathematical formulas are expressed as:

5.   Results and discussion

Using ANN models, the Casson flow toward two distinct heated surfaces is studied. The heat transfer coefficient and Casson fluid temperature are included in the quantity of interest. A mathematical formulation and shooting method is used to solve flow equations.

By taking the curvature parameter's zero value into consideration, the Casson fluid flow over a flat surface is assessed. The most important flow parameters for studying the flow field are the thermal radiation parameter, Prandt number, heat generation parameter, and variable thermal conductivity. We have assessed these parameters' effects on the NN at flat and cylindrical surfaces see Tables 3-10. Table 3 shows the variation in the NN in the presence of positive values for the variable thermal conductivity parameter. Both plates and cylinders are subject to this kind of inspection. The heat transfer normal to the surface is measured as a Nusselt number. The NN is seen to decrease for variations in the variable thermal conductivity parameter. For cylindrical surfaces, the Nusselt number is larger in magnitude. Table 4 shows the influence of a varied thermal conductivity parameter on the NN in the absence of heat generation.

The NN permits the inverse relation with varied heat conductivity for both surfaces. It is significant to notice that the Nusselt number's strength is higher for flat surfaces than for cylindrical surfaces. At the plate and cylinder, the effect of the Prandtl number on the NN is explored (see Table 5). It is obvious that the NN displays levels of stimulation for higher Prandtl numbers. When the influence of heat generation is taken into account, the effect of the Prandtl number on NN is explored (see Table 6). We have seen, the NN rises as the Prandtl number does for both plates and cylinders. Both plate and cylindrical surfaces are considered in the evaluation of the Casson fluid parameter's impact on the Nusselt number. About the scenario of heat generation, Table 7 is provided. We found that the Nusselt number declines normally when Casson fluid parameter values increase. Moreover, cylindrical surfaces have a somewhat greater Nusselt number magnitude.

In Table 8, the NN at both surfaces are analyzed in relation to the Casson fluid parameter for the case when heat generation has no impact. We found that heat transfer normal to the surface tends to incline for larger Casson fluid parameter values.

Table 9 shows how the thermal radiation parameter affects the NN. NN increases as thermal radiation strength does. Also, in the case of cylindrical surfaces, the NN is larger. When the heat generation effect is absent, Table 10 shows how the thermal radiation parameter affects the NN at both surfaces. For both frames, we saw that greater radiation parameter values lead to larger NN values. It is significant to notice that in the case of flat plates, the Nusselt number strength is higher. Our interest is to predict the values of the NN by using ANN models. For the prediction of the Nusselt number at both surfaces, we have built two separate ANN models (NN-ANN-I and NN-ANN-II). At a flat surface, NN-ANN-I owns the Nusselt number values, whereas NN-ANN-II owns the Nusselt number values at a cylindrical surface. One crucial step in the development of an ANN model is network training. The network was trained using the Levenberg-Marquardt method. Figures 2 and 3 show the training results for NN-ANN models. Specifically, Figure 2 presents the training of the ANN model for the forecast of NN at a flat surface, whereas Figure 3 for the prediction of NN at a cylindrical surface. For NN-ANN-I, the finest validation (FV) is 0.029329 at epoch 5 and for NN-ANN-II, the FV is 0.0.0099961 at epoch 29. We noticed that the MSE values for both graphs start out high and then start to decline as the stages go-on. The error histogram (EH) of the NN-ANN models being built to forecast Nusselt number values is shown in Figures 4 and 5. The EH of model-I is shown in Figure 4, whereas the EH of the model-II for a cylindrical surface is shown in Figure 5. The EHs in both cases have 20 Bins.

Since both NN-ANN models' error values are quite low, the training stages of NN-ANN-I and NN-ANN-II to estimate Nusselt numbers at both surfaces are successfully finished. The mean square error graphs for the Nusselt number are shown in Figures 6 and 7. Here, Figure 6 shows the MSE for the NN at the plat surface, and Figure 7 shows the MSE for the NN at cylinder. Both figures show that NN-ANN-I and NN-ANN-II have successfully completed their learning stages. For use of 80 data samples, the MSE values are displayed. The MSE values' proximity to zero indicates that the NN-ANN models were trained with fewer mistakes to predict the Nusselt number. While MSE is recorded as 0.03261 for flat surface, and the average MSE for cylindrical surfaces is 0.008562777. The regression outcomes of Nusselt number data for both surfaces are given in Figures 8 and 9. It is key to remember that the regression value provides a correlation between the objectives and forecasts values; if this value is near to one (1), we can state there is a close link.

The developed NN-ANN-I and NN-ANN-II are the best models to forecast the NN at both surfaces, according to Figures 8 and 9.

Figures 10 and 11 compare the Nusselt number's actual values with those predicted by an artificial neural networking model. Figure 10 compares NN-ANN and the target data set for plate whereas Figure 11 compares NN-ANN and the target data set for cylinder. We can note from figures that the bulk of NN-ANN model outputs is in good agreement with the Nusselt number target values. Due to this overlap, we may conclude that the generated models, NN-ANN-I and NN-ANN-II, are highly accurate in predicting Nusselt numbers.

The effect of a positive modification in the ${H}_{p}$ on temperature is shown in Figures 12 and 13. Figure 12 shows the impact of heat generation on temperature for a plate, while Figure 13 shows the impact of heat generation on temperature for a cylindrical surface. It can be demonstrated that the temperature of the Casson fluid and the heat generation parameter are directly related. Positive values of the ${H}_{p}$ cause energy generation, which raises the temperature. The effect of the Pr on the temperature is seen in Figures 14 and 15. In more detail, Figures 14 and 15 show how the Prandt number affects the temperature of a fluid over a flat surface and a cylindrical surface, respectively. We have seen that the Prandt number enables an indirect relation to Casson fluid temperature.

As the Prandtl number admits an inverse relationship to thermal conductivity, larger Prandtl number values result in a decrease in thermal conductivity. The fluid temperature drops as a result. For both flat and cylindrical surfaces, the effect of varying thermal conductivity on Casson fluid temperature is examined. Figures 16 and 17 provide evidence in this regard. Figure 16 describes the effect of variable thermal conductivity on the temperature in more detail, and we found that temperature rises as the ${\varepsilon }_{T}$ rises. Similar to Figure 16, Figure 17 demonstrates how ${\varepsilon }_{T}$ affects temperature over a cylinder. We found that fluid temperature exhibits tendencies that are indicative of variations in ${\varepsilon }_{T}$. Table 11 provides performance details for both ANN models, and it is clear that the constructed models are capable of accurately forecasting the values of the heat transfer coefficient at both flat and cylindrical surfaces.

6.   Conclusions

Thermally magnetized Casson stagnation flow over an inclined surfaces is studied numerically. Mixed convection, thermal radiations, heat generation, viscous dissipation, and temperature-dependent thermal conductivity all contribute to the flow. Using artificial neural networking models, the NN is evaluated at both surfaces toward five various flow factors. The following are the key outcomes:

• For an inclined cylindrical surface, an average Mean Square Error (MSE) is recorded as 0.008562777 while for a flat surface it is recorded as 0.03261.

• The values of the coefficient of determination suggest artificial neural networks as the best method for predicting the NN at both surfaces.

• For both surfaces, the temperature admits direct relation for heat generation and thermal conductivity parameters.

• For higher Prandtl number, the Casson fluid temperature shows decline values. The such trend holds for both surfaces.

• According to the predictions made by NN-ANN models, we found that for cylindrical surface, the strength of the NN is greater than the flat surface.

• The Nusselt number admits indirect relation for variable thermal conductivity, and Casson parameters while the opposite is the case for thermal radiation parameter and Prandtl number.

• A neural networking approach can be used to predict the surface quantities being involved for the fluid flow over non-linear stretched surfaces having engineering standpoints.

Acknowledgments

The authors would like to thank Prince Sultan University, Saudi Arabia, for the technical support through the TAS research lab. Author (Haitham M. S. Bahaidara) is thankful for the financial support provided by the Interdisciplinary Research Center for Renewable Energy and Power Systems, King Fahd University of Petroleum and Minerals, for funding this research through Project no: INRE2212.

Conflict of interest

The authors declare no conflict of interest.