Numerical and intelligent neuro-computational modelling with Fourier’s energy and Fick’s mass flux theory of 3D fluid flow through a stretchable surface

Current work explores the intricacies of magnetohydrodynamic and mix convectional boundary-layer flow concerning couple stress Casson nanofluid (CSCNF) dynamics via a 3D stretchable surface. The addition of active and passive control mechanisms for nanoscales creates an innovative dimension to the exploration. Remarkably, the analysis incorporates the influence of non-Fourier and non-Fickian heat and mass flux, alongside the effects of thermophoresis and Brownian diffusion, to systematically investigate the heat and mass transportation phenomena. The governing equations (PDEs) describing the MHD-3DCSCNF model are converted into a set of ordinary differential equations to facilitate the ANN analysis. Employing the bvp4c technique, a dataset is systematically generated for the back propagation artificial neural network with Levenberg–Marquardt Algorithm (BANN-LMA) through four distinct scenarios. Via accurate testing, validation, and training, the BANN-LMA produces estimated results for the MHD-3DCSCNF problem. The performance validation of BANN-LMA is executed through several metrics, involving the mean squared error, error histogram and regression analysis. The training process, characterized by minimizing the MSE through a gradient descent methodology with optimized weights, exhibits a compelling correlation R = 1, between the target and network output. Furthermore, the consistent convergence observed highlights the method robustness and reliability.


Introduction
An important and cutting-edge technique to artificial intelligence is exemplified by artificial neural networks (ANNs).ANNs are remarkably adaptable, dynamically responding to the data that flows through the network while learning, whether it comes from within or outside the system.Artificial neural networks make use of the Back-Propagation approach for concurrent training to improve the performance of Multilayer Perceptron networks.This approach stands out as the most common, effective, and understandable framework for complex multi-layered networks.Backpropagation is a supervised learning technique that follows the gradient of the error curve to reduce network errors.Notably, Paul Werbos first introduced this algorithm in 1974.Rumelhart and Parker just revived it.In feedforward multilayer neural networks, backpropagation is now a frequently used learning technique.The Levenberg-Marquardt methodology offers computational solutions to a wide range of fluid flow problems in various technical disciplines by introducing a novel and convergent stability method for ANNs.Newtonian and non-Newtonian fluid systems have been the subject of recent research, which used the Levenberg-Marquardt artificial neural network (LMA-NN) for experimentation.To examine the impact of radiative heat transport on the Williamson fluid flow model, (Shoaib, Raja, et al., 2022) reported a groundbreaking use of the LMA-NN.By considering mixed convectional and slip impact between rotating disks, (Zhao et al., 2023) used LMA-NN to systematically examine the dynamics of Ree-Eyring fluid flow.Furthermore, led to a deeper knowledge of the fluid, (Asif Zahoor Raja et al., 2022) used the computational intelligence algorithms based on LMA-NN to learn more about ferrofluid flow, incorporating the effects of both magnetic forces and nonlinear thermal radiation effects.Recent studies (Dai et al., 2023;Karimmaslak et al., 2021;Z. Shah et al., 2022;B. K. Sharma et al., 2023;Shoaib, Kainat, et al., 2022;X. Zhang et al., 2022;Zuhra et al., 2022) that use the ANN approach to analyse heat transport mechanisms in various geometries while taking a variety of characteristics are presented in the literature.
In thermodynamic processes, irreversibility is caused by several phenomena, including fluid opposition, radiative heating, biological repetition, energy transmission, chemical changes, and electric and magnetic propagation.Various mechanisms such as energy transportation and movement of liquid inherently demonstrate irreversibility, supporting enhanced entropy and the loss of favourable properties.In these processes, the mechanism of entropy production remains constant during closed isothermal processes, but in case of irreversible, the entropy production increases.Such distinction is helpful in identifying reversible and irreversible mechanisms.When liquid flow experiences a drop in pressure, entropy is produced in any thermodynamic system.Bejan (1980Bejan ( , 2013) ) first investigated the idea of irreversibility in convection liquid flow because of the viscous effects of thermal transport.Tayebi et al. (2020) highlighted the numerical analysis of entropy production in the flow field of a nanofluid in presence of convectional heating.Khan et al. (2020) explained on the statistical modelling of slip flow in a micropolar liquid.Hayat et al. (2020) evaluated entropy-optimized flow with melting heat exchange.The rate of entropy production and nonlinear radiated thermal fluxes over a deformable surface with changeable properties were investigated by Nasir, Berrouk, Aamir, et al. (2023).Entropy improvement in high-temperature surfaces, including nanoparticles, was studied by Boujelbene et al. (2023).In the domain of literature, intriguing and contemporary flow-related challenges linked to the generation of entropy are explored in references (Mandal & Pal, 2023;Rafique et al., 2023;Sakkaravarthi & Reddy, 2023).
The quality of energy transmission becomes fundamental as distinct areas inside the same body undergo thermal variations.Fourier first proposed the idea of coolant temperature conduction, and Catteneo later improved it by adding a relaxation time to address the paradox of heat conduction.The Cattaneo-Christov heat flux model (CCHFM), which incorporates the notion of relaxation time, was created because of Christov's further development of the Fourier model.The CCHFM was used by Nayak et al. (2022) to investigate heat transmission in hybrid nanofluids and the subsequent entropy formation.The CCHFM was used by Gowda et al. (2021) to examine heat transmission in a nanofluid flow.To study the movement of chemically reactive micropolar carbon nanotube-blood nanofluid through a constriction channel, (Shaw, 2021) applied the idea of CCHFM mechanism.Using the CCHFM, (Gowda et al., 2023) investigated the flow of liquid nanomaterial over a surface.The effect of the CCHFM on the flow of ferromagnetic Maxwell nanofluids across an extended geometry was explained by F. Shah et al, (2023).
The phenomenon of thermal radiation refers to the transmission of thermal energy between two bodies in the absence of physical contact.Thermal radiation is used by pyrometers and infrared cameras to measure and regulate temperature in a variety of sensing processes.Thermal radiation plays a dynamic role in a wide range of engineering applications, especially at high temperatures.Convective heat transfer was examined by Jayaprakash et al. (2022) together with radiative heat transmission.The impact of a homogenous heat source/sink on the flow of hybrid nanofluids was explored by Soumya et al. (2022).In order to account for heat radiation (Nasir, Berrouk, Gul, et al., 2023), described the behaviour of nanofluids over a surface.The effects of heat sources and radiation on the convective flow of a nanoliquid in a cavity were discussed by Gul et al. (2023).Thermal radiation has great uses in the field of medical sciences in addition to industrial and engineering domains.As a result, (Hsu et al., 2017) evaluated potential application in the field of medicine, providing illumination on how thermal radiation from the human body might result in both heating and cooling effects.The vital importance of maintaining a healthy body temperature as a basic requirement for living is emphasized by their research.In a similar fashion, (Wehinger & Flaischlen, 2019) established frameworks for thermal radiation in computational fluid dynamics, (Petela, 2010) concentrated on industrial thermodynamics, (Navarro et al., 2021;Wang et al., 2023) investigated thermal radiation within fluid dynamics.All the authors researched the novel and advanced uses of thermal radiation in various contexts.
The properties of mixed convective flow across an extended surface have numerous technological and industrial applications.When gravitational forces are at work, the effect of free convection is more prominent.The mechanisms of fluid flow and energy transmission are impacted by both the strained and buoyant forces.Thermal buoyancy, which affects the rate of heat movement in manufacturing environments, results from variations in the heat on the expanding surface.Applications of this phenomenon include precautionary stoppage of nuclear cooler systems, environmental fluxes, heaters, photovoltaic panels, conditioning of electrical gadgets, oceanic mechanics, lowvelocity heat transfer mechanisms, precooler mechanisms, vehicle demisters, and more.According to (Patil et al., 2020), mixed convective nanoparticle transport along an exponentially expanded porous surface was influenced by suction/injection.Mixed convective radiative flow over an inclined permeability surface was explored by Moradi et al. (2013).Time-dependent magneto Williamson nanoparticle flow through a stretchy surface was evaluated by Das et al. (2022) considering velocity slips, nonlinear mixed convection, and radiative heat effects.In contrast to the focus on their Newtonian behaviour, the impact of nanofluids' non-Newtonian properties on their heat transfer rate has not been adequately investigated.Many nanofluids display non-Newtonian characteristics, even though a lot of study has been devoted to employing Newtonian nanofluids to improve energy transmission (Elboughdiri et al., 2023;Rasool et al., 2023;J. Sharma et al., 2023;Wakif, 2023;Wakif et al., 2022;K. Zhang et al., 2023).
To the author's knowledge, most of the literature in the area of Casson couple-stressed nanofluids concentrates on convective boundary conditions rather than active and passive control nanoparticle boundary conditions for examining thermal and concentration distributions.The analysis of the double-diffusive Cattaneo-Christov energy flux in couple stress nanofluids using the theory from reference (Hayat et al., 2017) while simultaneously considering passive and active control conditions has also not been the subject of any published study.Additionally, this theoretical and numerical studies have the potential to advance the field of bio-fluid dynamics and have applications in several engineering fields, including photovoltaic conversion, nuclear reactor implementation, cooling systems, and electronic generators.So, according to the authors, it helps to investigate thermal radiation, magnetohydrodynamic effects, and Darcy-analysis Forchheimer's of couple stress nanofluids' three-dimensional mixed convection flow.Because of this, authors feel it is appropriate to investigate the effects of passive and active control on magnetohydrodynamics, thermal radiation, Darcy-analysis Forchheimer's of three-dimensional mixed convection flow of couple stress nanofluid with non-heat Fourier's flux, and nonmass Fick's flux theory.The couple stress nanofluid flow of magnetic, mix convectional, and double-diffusive Cattaneo-Christov energy flux within porous media is investigated in this research using the BANN-LMA approach.The importance of the constructed ANN procedure is emphasized through a comparison of convergence and accuracy measures, employing performance analysis, regression, correlation, histograms, and the Mean Squared Error.The reference dataset for the ANNs is established using the bvp4c technique for various model parameter variations.Additionally, using graphs and tables, the effects of various parameters on velocity, temperature, and concentration fields are thoroughly explored and explained.

Mathematical formulation
Here, a mathematical model for steady, incompressible Couple stress nanofluid flowing through a nonlinear extending sheet submerged in a porous media created by Darcy and Forchheimer is represented.In this analysis the surface is traced at z = 0 and the fluid flow is produced as the extending of sheet surface in x− and y−directions.Also, we assume that the starching velocity along x−direction is U w = a (x + y) n and along y−direction is V w = b(x + y) n (Gowda et al., 2023;Ramzan et al., 2013).Wherea, b and n are positive constants.Its also suggested that T w −constant surface and T ∞ −ambient temperature such that T w − T ∞ > 0. Furthermore, we settled that C w −constant surface and C ∞ −ambient temperature such that C w − C ∞ > 0. In this effort, the Buongiorno nanofluid model and Cattaneo-Christov heat and mass fluxes are also used (Figure 1).

System assumptions and limitations
In the current flow system, we have implemented the following assumptions and constraints:

and mass fluxes models
The velocity, temperature, and concentration fields for three-dimensional steady incompressible fluid flow can be represented as V = [u(x, y, z), v(x, y, z), w(x, y, z)], T = T(x, y, z) and C = C(x, y, z).For an incompressible couple-stress nanofluid, the fundamental equations of conservation of mass, motion, heat, and concentration are (Gowda et al., 2023;Ramzan et al., 2013;Zuhra et al., 2022): The heat flux ( q) and mass flux ( J) model is presented as follows (F.Shah et al., 2023): where E −relaxation time heat flux and C − relaxation time mass flux.Equation reduces classical Fouriers and Fick's law for E = 0 = C (Shaw, 2021;Zuhra et al., 2022).Since the fluid flow in this suggested model is mass-conserving, or steady, at ∇. Ṽ = 0, Equation ( 5) become.
For the isotropic and incompressible flow of a Casson fluid, the rheological equation of state is as follows (J.Sharma et al., 2023;Wang et al., 2023;Zuhra et al., 2022).
Here π − is the product of the component of deformation rate with itself, where , the dynamic viscosity could be described as We obtain by inserting the value of py The high order system of partial differential equations formed from the controlling problem are as follows, based on the above-mentioned explanations (Gowda et al., 2023;Ibrahim et al., 2022;Ramzan et al., 2013;Shaw, 2021;Zuhra et al., 2022): where u and v (velocities components), ), E and C (thermal and concentration relaxation time coefficient), α f = k/(ρc p ) f , k(thermal diffusivity and conductivity), τ = (ρc p ) p /(ρc p ) f (heat capacities ratio), (ρc p ) f (effective heat capacity of fluid), D B (Brownian diffusivity), D T (thermophoretic diffusion coefficient).

Boundary conditions
The relevant boundary conditions are as follows (Elboughdiri et al., 2023;Ibrahim et al., 2022;Rasool et al., 2023): The arising terms χ E and χ C from Equations ( 13) and ( 14) are estimated as: Considering the appropriate similarity transformations as (Ibrahim et al., 2022): For active control of φ, where the velocity components for present investigations are represented as The aforementioned equations can be substituted in governing partial differential equations which transform them into couplet ODEs of the following form: With transform boundary conditions: where the term δ * 1 and δ * 2 are where the several model parameters obtained after transformation are presented in Table 1.
The engineering aspects for the current analysis are Nu and Sh are given by Gowda et al. (2023), Ibrahim et al. (2022), Ramzan et al. (2013) and Zuhra et al. (2022): Here terms q w → surface heat flux and p w → mass flux mathematically stated as: In non-dimensionless form is When the mass flux is equal to zero, then the Sherwood number will also equal zero.
Equation ( 29) signifies the Local Reynolds number.

Formation of data set
Assuming and organizing the relevant information for training the BANN-LMA or for resolving fluid mechanics related problems is referred to as the process of constructing a dataset for projected solutions.Such data set operates as the foundation for training and evaluating the performance of several algorithms, containing neural networks.The choice of representative models, data preparation, and ensuring data integrity and quality are all carefully considered.The dataset's extensiveness and structure play a crucial role in assisting the model to generalize and produce dependable results.Consequently, the creation of a well-crafted dataset stands as a critical initial phase when utilizing neural networks for successful solutions.With regards to addressing the nonlinear couple stress Casson fluid model (CSCFM), the strategies and techniques presented here are characterized into the following three parts.The classification of present nonlinear CSCFM fluid model contains three categories, specifically velocity field, temperature field, and concentration field.Within each classification, there are two distinct cases.The first phase offers the principal data for the numerical solver bvp4c.This solver is working to establish the dataset based on final couplet system of equations.This consists of transforming higherorder nonlinear ODEs into first-order ODEs while containing the boundary conditions.The succeeding phase presents a suggested computational framework, BANN-LMA, which involves elements like the backpropagation procedure, layer layout, and network plan.The concluding phase offers guidelines for applying the suggested BANN-LMA methodology (Asif Zahoor Raja et al., 2022; Z. Shah et al., 2022;Shoaib, Raja, et al., 2022) to address the tasks of the nonlinear CSCFM.The resulting higherorder ODEs after transformation are (Ibrahim et al., 2022): With transform boundary conditions: π 13 (0) = 1, → For active control, Nb.π 14 (0) + Nt.π 12 (0) = 0, → For passive control

Creating artificial neural network
In programming courses, artificial neural networks (ANNs) are a crucial tool.Artificial neural networks, which were originally utilized in the middle of the 20th century, gained popularity because of their improved ability to handle complex nonlinear functions (Dai et al., 2023;Shoaib, Raja, et al., 2022).The most well-known multilayer perceptron (MLP) systems are used in artificial neural network (ANN) models (Zhao et al., 2023).The neural network is made up of (i) weights, (ii) input/output, (iii) activation mechanism, and hidden layer.Backpropagation feed-forward (BPFF) is the most popular approach used in multilayer perceptron networks (Shoaib, Kainat, et al., 2022).The network then creates a relationship with at least one hidden layer, where the weighted connection's structure is doing critical computations.Towards the end of training, error eventually and θ(η) using 410 data sets from 8 distinct data components.349 samples were used for training in this analysis, whereas 61 data sets were used for testing and validation.

Time relaxation parameter of concentration
Hartman number Porosity parameter → For passive control of φ

Results and discussion
Using MATLAB 'bvp4c' numerical solver, the nonlinear differential equations and their related boundary conditions are solved numerically.The original differential equations are changed into first-order form to facilitate this procedure.By changing the values of the physical model parameters, as shown in Table 2, several scenarios are investigated.Figure 3 depicts the full computing process of current analysis.

Effect of model parameters on heat and concentration profile
Following the implementation of the MATLAB 'bvp4c' numerical solver to acquire the solution of the nonlinear ODEs, the solution is drawn to show the presence of numerous significant factors, including the M, γ , K * , β, Nt, Nb, Pr and ω, on the (ξ ) and φ(ξ ) fields.Results are obtained for both passive and active control of nanoparticles. Figure 4(a,b) displays the variation in (ξ ) and φ(ξ ) for varying values of M.An enhancement in M signifies an augmentation in the (ξ ) and φ(ξ ) profiles for both the active and passive control cases.The magnetic field, which lowers the liquid's movement and increases the heat and concentration profiles, causes the body force to diminish as would be expected.Additionally, from Figure 4(a) the (ξ ) distribution for passive control is smaller than it is for active control.Figure 4(b) illustrates how the M affect the φ(ξ ) profile.
By increasing the M, the φ(ξ ) profile is enhanced for both active and passive controller scenarios.Figure 4(c) shows the effect of the γ stimulation on (ξ ) field.To produce stronger convection, the γ can be elevated.As a result, (ξ ) field and thermal layer thickness rise under control situations.Figure 4(d) is developed to demonstrate how the γ impacts the φ(ξ ) profile.In this situation, the φ(ξ ) profile shows a propensity to rise in both active and passive nanoparticle control conditions.The effect of Nb on the (ξ ) and φ(ξ ) field is shown in Figure 4(e,f) for both cases.The collision of fluid atoms or molecules will result in an arbitrary motion known as the Brownian motion of suspended (pendulous) particles, which will enlarge the boundary layer in both active and passive control scenarios, as shown in Figure 4(e).Figure 4(f) demonstrates that choosing a higher Nb parameter value results in a weaker φ(ξ ) field for both types of nanoparticles (Table 3).

Analysis of arterial neural network model results
To implement the suggested solution methodology, reference solutions in the form of the numerical data set produced via bvp4c against the F (ξ ), (ξ ) and φ(ξ ) are delivered to the neural network input.The neural network is trained with 75% of the points, with the remaining 15% being used for testing and validating the     in total time 17s with MSE is 2.6915×10 −10 , 8(d) after 78 epochs in total time 21s with MSE is 3.7068×10 −11 ) respectively.

Analysis of comparison of suggested and ref solutions with absolute errors
In comparison to the reference dataset of CSCFM, specifically the parameter K * , Figure 9 shows the resultant graphs and evaluation of the variation in F (ξ ) using BANN-LMA.We investigate the effect of the K * on F (ξ )in Figure 9(a).The behaviour of K * depends on the couple stress viscosity, which acts as a barrier to flow and increases fluid density, which in turn increases fluid motion.Figure 9(a) demonstrates clearly how fluid velocity inside the flow zone accelerates as value of K * increases.An assessment of F (ξ ) s absolute errors is shown in Figure 9(b), with the absolute errors lying within the range 10 −4 and 10 −8 .While observing Figure 10(a,b), it becomes evident that F (ξ ) profile experience a reduction as the value of Mincreases, and this is accompanied by absolute error ranges between 10 −3 and 10 −7 .This decline in F (ξ ) profile is mostly because of the Lorentz force produced by the applied magnetic field.It's significant that the Lorentz force resists the motion of liquid, preceding to a reduction in its overall motion.Furthermore, such force creates additional heat, identified as Joule heating, within the system.Therefore, the employment of higher value of M in the system results in an expansion in (ξ ) profile within the flow region, which is spotted in both active and passive control cases, as indicated in Figure 11(a,b) with absolute error ranges ranging from 10 −5 to 10 −8 .In Figure 11(a,b), we display the visual interpretations of how Rd impact the (ξ ) in both scenarios, while simultaneously operating an AE analysis within the range of 10 −5 to 10 −8 .Rd plays an extensive role in governing the transfer of heat energy within and between the liquid and its surroundings.The Rd serves as a quantification of the impact of radiation in the situation of heat transport.A higher value Rd specifies an elevated status of radiation effects, indicating to higher rates of heat transport and alterations in thermal distribution.As the Rd escalates, the surface temperature also experiences an increase, then ensuing a rise in both the thickness of the thermal boundary layer and thermal resistance.Also, observing Figure 12(a,b), it becomes evident that φ(ξ ) profile experience a reduction as the value of γ increases in case of active control while for passive control present opposite behaviour, and this is accompanied by AE ranges between 10 −5 and 10 −7 .Similarly, the employment of higher value of Nb in the system results in an expansion in φ(ξ ) profile within the flow region, which is spotted in both active and passive control cases, as indicated in Figure 13(a,b) with absolute error ranges ranging from 10 −3 to 10 −8 .
Figure 14 (a) shows comparisons between the current study and existing literature (Ramzan et al., 2013), whereas Figure 14(b) emphasizes the numerical comparison between numerical and ANN outcomes.
Tables 4-7 show the significance of each variation within the suggested methodology and were individually prepared for each of the four scenarios (1-4).Figure 15(a) illustrates the contrast between our present work and the existing literature (Ramzan et al., 2013), whereas Figure 15(b) demonstrates the validation of our findings by comparing the ANN results with the numerical dataset.Both figures exhibit a high level of concordance, thereby guaranteeing the precision of our current work and preventing any potential confusion or misinterpretation.

Conclusions
The present investigation aims to explore the impacts of coupling stress Casson fluid models with convectional properties under active and passive control using advanced AI-based computational methodologies.These approaches power the neural networks enhanced with the Levenberg-Marquardt Algorithm (NN-LMA).The investigation incorporates the Cattaneo-Christov  heat and mass flux models and employs the Rosseland approximation to account for thermal radiation effects.The basic partial differential equations (PDEs) explaining flow patterns are transformed into ODEs through appropriate conversion manner.This system of equations are then solved using the robust bvp4c       demonstrating a strong agreement between expected and actual values.Additionally, the regression value constantly equals 1 in each case, highlighting the precision of the results.The main objective of the current analysis is described below, focusing on how the physical characteristics of couple stress Casson fluid models change.
• The M parameter reduces F (ξ ) profile and results in a thinner momentum boundary layer, while enhances the fluid's thermal behaviour.• Increasing the K * and λ parameter leads to an acceleration of fluid motion.
• (ξ ) distribution can be controlled by adjusting the Rd, Nt, Pr and γ to finer values.
Both temperature and concentration distributions are more pronounced in cases of active control compared to passive control.
• As the Nt increases, the mass transfer rate declines, while the Nb exhibits the opposite trend.• A higher γ signifies increased convective effects, which enhance the temperature.• In the future, one could employ the remarkable capabilities of artificial intelligence-based methodologies,   utilizing both supervised and unsupervised neural networks, in various applications within the field of fluid mechanics such as (Nasir, Alghamdi, et al., 2023;Ramzan et al., 2013;Wang et al., 2023;K. Zhang et al., 2023).Future study may also expand on the current research to include non-Newtonian liquids and unsteady behaviour under diverse physical conditions.
To ensure that the study continues to make significant contributions and remains current in its specific field, these suggested areas for extension might serve as a road map for future research directions.

Nomenclature
List of symbols

Figure 1 .
Figure 1.Physical model of the problem.
of this product based on the non-Newtonian model, μB − is the plastic dynamic viscosity of the non-Newtonian fluid and py − is the yield stress of the fluid (Sakkaravarthi & Reddy, 2023), but for Casson fluid we supposed π c < π and py = √ 2πμ B β

Figure 2 .
Figure 2. (a) The structure of signal neural network model, (b) Construction of neural network for CSCFM.
i , w i is the weighted sum of input connection has been included in the present research as the transformation function /activation function of the MLP network.Among the most popular transformation functions is the sigmoid function.Its precise performance on nonlinear regression issues in various applications leads to its adoption in MLP networks.During the training phase, errors were evaluated using a performance measure.In the present research, the determination MSE and R 2 (correlation coefficient) have been used to assess the network's accuracy.Which is define as (Asif Zahoor Raja

Figure 3 .
Figure 3. Working flowchart for current analysis.

Figure 5 .
Figure 5. Performance results variations of Case 1, scenario I.
Figure 5(a) displays the transient stability for case (1) in scenario (I).An error histogram, which breaks down the full error range into 20 bins, is used to analyse the dispersion of errors.The total number of errors for the network's training, validation, and testing stages is shown as a vertical line for each error range.The reliability and correctness of the solution outcomes are indicated by the presence of maximum error values near to the zero-error line is −7.5 × 10 −6 .The error distribution for case (1) in scenario (I) is displayed in Figure 5(b).Another helpful technique for evaluating the consistency between desired and actual values for the training, validation, and testing phases separately is regression analysis.This evaluation can be done mathematically or visually.A straight line is used to represent the reference solution in graphical representations, while BANN-LMA results are represented by little rings or spots.An R = 0 denotes a lack of correlation between the targeted and obtained results, whereas an R = 1 indicates outstanding correlation.A regression plot for case (1) in scenario (I) is shown in Figure 5(c).When analysing the quality and reliability of the machine learning model, the MSE (mean square error) is an essential aspect.MSE is calculated numerically by averaging the squares of the discrepancies between the values that were produced and a benchmark solution.A lower MSE value indicates greater accuracy and reliability in the findings, demonstrating the dependability of the solution technique.The MSE plot for case (1) in scenario (I) over the training, validation, and testing data sets is shown in Figure 5(d) after 445 epochs in total time 28s with MSE is 7.7741 × 10 −10 .For scenarios II, III, and IV, the subfigures of Figures 6-8 contain time-series response plots (Figures 6(a), 7(a), 8(a)), error histograms (Figure 6(b) with maximum error values near to the zero-error line is −7.11×10 −6 , 7(b) with −1.5×10 −6 , 8(b) with −9.6×10 −7 ), regression graphs (Figures 6(c), 7(c) and 8(c)) and MSE plots (Figure 6(d) after 336 epochs in total time 22s with MSE is 1.804×10 −9 , 7(d) after 180 epochs

Figure 6 .
Figure 6.Performance results variations of Case 1, scenario II.

Figure 9 .
Figure 9.Comparison of suggested and ref solutions and absolute errors for various values of K * .

Figure 10 .
Figure 10.Comparison of suggested and ref solutions and absolute errors for various values of M.

Figure 11 .
Figure 11.Comparison of suggested and ref solutions and absolute errors for various values of M.

Figure 12 .
Figure 12.Comparison of suggested and ref solutions and absolute errors for various values of Rd.

Figure 13 .
Figure 13.Comparison of suggested and ref solutions and absolute errors for various values of γ .

Figure 14 .
Figure 14.Comparison of suggested and ref solutions and absolute errors for various values of Nt.

Table 1 .
Details regarding the dimensions and arrangement of the proposed network.

Table 3 .
Variation in different parameters of ANN for various scenarios.
solver to produce datasets for ANN analysis.Subsequently, the BANN-LMA method is applied to the reference dataset to approximate solutions for the CSCFM models.As shown by a minimal MSE of about 10 −10 in the training charts, the findings of this study show a high degree of agreement between the approximations and the real data.In error histograms, around 90% of the output dynamics are centred around the zero line,
Radiation parameter C f , C w , C ∞ Fluid, wall, free space concentration (C p ) f EThermal relaxation time coefficientEConcentration relaxation time coefficient βCasson parameter