The intelligent networks for double-diffusion and MHD analysis of thin film flow over a stretched surface.

This study presents a novel application of soft-computing through intelligent, neural networks backpropagated by Levenberg-Marquardt scheme (NNs-BLMS) to solve the mathematical model of unsteady thin film flow of magnetized Maxwell fluid with thermo-diffusion effects and chemical reaction (TFFMFTDECR) over a horizontal rotating disk. The expression for thermophoretic velocity is accounted. Energy expression is deliberated with the addition of non-uniform heat source. The PDEs of mathematical model of TFFMFTDECR are transformed to ODEs by the application of similarity transformations. A dataset is generated through Adams method for the proposed NNs-BLMS in case of various scenarios of TFFMFTDECR model by variation of rotation parameter, magnetic parameter, space dependent heat sink/source parameter, temperature dependent heat sink/source parameter and chemical reaction controlling parameter. The designed computational solver NNs-BLMS is implemented by performing training, testing and validation for the solution of TFFMFTDECR system for different variants. Variation of various physical parameters are designed via plots and explain in details. It is depicted that thin film thickness increases for higher values of disk rotation parameter, while it diminishes for higher magnetic parameter. Furthermore, higher values of Dufour number and the corresponding diminishing values of Soret number causes enhancement in fluid temperature profile. Further the effectiveness of NNs-BLMS is validated by comparing the results of the proposed solver and the standard solution of TFFMFTDECR model through error analyses, histogram representations and regression analyses.

Thin film flow can be expressed as thin liquid layer flowing over substrate or surface due to gravity or shear stresses (the external forces) and the holding free boundary is called free surface. Free surface is that one, at which, the parallel shear stresses vanish. Liquid’s free surface is that one where it is in contact with the air. The applications of thin film liquid flows have very wide range uses in the mass and heat transfer devises, printing technology, spinning disks reactors, drop spreading processes, microelectronic industry, lubricating gears, distillation columns. Moreover, coating or spin coating (thin film flow) is necessary for improvement of durability, cooling or heat transfer process and work efficiency. Emslie et al.[1] initiated the work on liquid film flow in research field, by considering the flow over a rotating disk. The stated work was extended by Jenekhe[2] and Flack et al. [3] by including heat and mass transfer effects and taking non-Newtonian fluids for analysis. Washo[4] performed experimental findings of spin coating problem. Nath and Kumari[5] carried out research work on thin film unsteady flow with disk rotation effect. Thin film flow along with various effects, like activation energy, radiation, Joule heating, nanomaterials, unsteadiness, thin film thickness Brownian motion and thermophoresis, over rotating disk were studied by various researchers[6-11]. The impact of double-diffusion and flow of thin film liquid over rotating surface have been extensively analyzed by many researchers[12-18].

The artificial intelligent based neural networks (NNs) is an effective and authentic solver, which has been used extensively in diversified fields of science. A brief survey is provided here, which provides the effectiveness of NNs technique. For instance, the effective applicability can be seen with applications include nonlinear HIV infection with CD4+ T-cells based mathematical model[19], the model of Maxwell nanofluid thin film flow over a rotating and stretched surface[20], the integrated system of power plant including stochastic wind[21], nonlinear dynamics system of COVID-19[22], the model of squeezing flow[23], singular nonlinear pantograph differential models[24], the piezostage actuator Bouc–Wen hysteresis system[25]. All the above-mentioned facts imply that the artificial intelligent base techniques have a variety of applications in sundry fields of science and technology. The application of computational intelligent neural networks processes can be helpful for the solution nonlinear systems of thin nanoliquid film flow model of Maxwell fluid with thermo-diffusion effect and chemical reaction (TFFMFTHECR) over rotating disk.

The notable features and motivations of this analysis are highlighted as:

A novel design of intelligent, neural networks backpropagated by Levenberg–Marquardt scheme (NNs-BLMS) is incorporated for numerical interpretations of TFFMFTHECR mode.

The system of PDEs of TFFMFTDECR model are transformed to system of ODEs by the use of suitable similarity transformations.

A dataset is generated through Adams method for the proposed NNs-BLMS in case of sundry scenarios of TFFMFTDECR model by varying rotation parameter, magnetic parameter, space dependent heat source/sink parameter, temperature dependent heat source/sink parameter and chemical reaction parameter.

The efficiency and performance of NNs-BLMS is verified by comparing the results with the standard solution of TFFMFTDECR model through analyses of error, graphical outcomes of histograms and regression.

The outcomes depicted that thin film thickness increases for higher values of disk rotation parameter, while it diminishes for higher magnetic parameter. Furthermore, higher values of Dufour number and the corresponding diminishing values of Soret number causes enhancement in fluid temperature profile.

The rest of the article is presented as: the problem TFFMFTDECR model formulation is given in “Mathematical model formulation and interpretation of the outcomes” section, the results of NNs-BLMS in case of involved parameters are elaborated in “Solution methodology and results” section, “conclusion” of the work is placed in the very last.

Mathematical model formulation and interpretation of the outcomes

The Maxwell fluid thin film unsteady MHD[26-28] flow over a rotating stretchable disk in cylindrical coordinate system with Soret and Dufour effects is considered. A first order chemical reaction is included. The liquid thin film flow is due to rotation and stretching behavior of the disk having velocities in directions, where is rotation rate, c is stretching rate and is a constant as shown in Fig. 1. The thickness of the Maxwell nanofluid film over rotating and stretching disk is . The disk surface concentration and temperature are mathematically defined as:  = , where concentration and temperature at origin are , while the reference constant concentration and temperature are . The flow is influenced by magnetic field in the z-direction. The mathematical representation of the above-mentioned physical problem is given by the governing equations[29,30] as:where non-uniform, space and temperature dependent, heat source/sink [31], thermophoretic velocity [32] and thermophoretic coefficient are given by:

Figure 1

The geometry of the transient thin film flow problem.

The boundary conditions for the thin film flow of Maxwell nanofluid problem are:where the triplet represents the velocity components corresponding to . The expressions appearing in the above equations are fluid density, thermal conductivity, temperature, concentration, time relaxation parameter, specific heat at constant pressure, temperature at surface, concentration at surface, thermal diffusion, diffusion coefficient for species, chemical reaction coefficient, kinematic viscosity, concentration susceptibility, temperature at origin, concentration at origin, space dependent heat source/sink, temperature dependent heat source/sink, thermal conductivity of fluid and fluid mean temperature, respectively. The positive values of and i.e. and corresponds to internal heat generation, while negative values of these parameters i.e. and represents absorption parameters, respectively.

The similarity variables are defined as:

Substituting Eq. (9) into Eqs. (1)–(5) and Eq. (8) the following system of ODEs with corresponding boundary conditions are obtained:

The physical parameters appearing in above equations are given by:where is unsteadiness parameter, is Magnetic parameter, is Deborah number, is Prandtl number, is Dufour number, is thermophoretic number, is chemical reaction parameter, is Soret number, is Schmidt number and is rotation parameter.

The value of parameter represents at free surface. Therefore, Eq. (9) gives:

The local Nusselt number, which represents heat transfer rate and Sherwood number, which is mass transfer rate are given by:

The expressions in dimensionless form of these quantities are:where, stands for local Reynolds number.

Solution methodology and results

The application of computational numerical solver, intelligent neural networks backpropagated by Levenberg–Marquardt scheme (NNs-BLMS), is applied to numerically investigate thin film flow of Maxwell fluid with thermo-diffusion effects TFFMFTDECR model. The solver NNs-BLMS is applied to process the dataset generated by Adams numerical technique. The process of neural networks is presented in the form of single biological neuron as well as mathematical equivalent structure is given in Fig. 2. The overall process of this study is presented as single flow chart in Fig. 3. Eight variants of TFFMFTDECR model are processed through ‘nftool’, using NNs-BLMS for the solution.

Figure 2

A single neural structure.

Figure 3

Work flow diagram of NNs- BLMS for TFFMFTDECR.

The TFFMFTDECR model, which is given in Eqs. (10)–(14) is treated numerically for eight variations , , , , , , and , where each variation represents a scenario. The eight scenarios are presented in detail with the help of Table 1 along with three cases for each scenario, while Table 2 represents the detail of non-varying parameters during the analysis.

Table 1

Illustration of eight scenarios (1–8) along with cases (1–3) for TFFMFTDECR model.

Scenario	Case	Description of variation of physical parameters
		β	ω	M	Variation of A and B for positive and negative values				Sr	Du	K_R
					A(-)	B(-)	A(+)	B(+)
1	1	2.21	1.0	1.0	–	–	0.2	0.2	0.2	0.2	0.3
	2	2.57	1.3
	3	3.14	1.9
2	1	2.73	1.0	0.0	–	–	0.2	0.2	0.2	0.2	0.3
	2	2.21		1.0
	3	1.70		3.0
3	1	2.21	1.0	1.0	− 0.5	0.1	0.2	0.2	0.2	0.2	0.3
	2	2.21			− 1.5
	3	2.21			− 2.0
4	1	2.21	1.0	1.0	–	− 0.5	0.2	0.2	0.2	0.2	0.3
	2	2.21				− 1.5
	3	2.21				− 2.0
5	1	2.21	1.0	1.0	–	0.1	0.5	0.2	0.2	0.2	0.3
	2	2.21					1.5
	3	2.21					2.0
6	1	2.21	1.0	1.0	–	–	1.0	0.5	0.2	0.2	0.3
	2	2.21						1.5
	3	2.21						2.0
7	1	2.21	1.0	1.0	–	–	0.2	0.2	0.07	1.0	0.3
	2	2.21							0.4	0.175
	3	2.21							1.0	0.07
8	1	2.21	1.0	1.0	–	–	0.2	0.2	0.2	0.2	0.0
	2	2.21									0.3
	3	2.21									0.9

Table 2

Description of physical quantities, which are non-varying during the study.

S	β_1	Pr	Sc	τ
1.0	0.6	1.7	1.0	0.9

The ‘NDSolve[33-36]’ tool is used in Mathematica environment to solve the system of ODEs of TFFMFTDECR for eight variations as described in Table 1. After generating the numerical reference dataset through Adams technique, it was exported to Matlab (Matlab 2021a, with Licence 40727596, and URL is as: https://www.mathworks.com/products/matlab.html) environment NNs-BLMS. The inputs used for cases (1–3) of scenarios (1–8) for step size 0.01, the detail is given in Table 3. The input values are distributed as 70%, 15%, and 15% for train, validation and testing, respectively. The designed NNs-BLMS with of two layered structure is given in Fig. 4.

Table 3

Illustration of eight scenarios (1–8) along with cases (1–3) for TFFMFTDECR model.

Scenario	Case	β	ω		Starting value		Step size		End value		Input Values in dataset for NNs-BLMS
1 (ω)	1	2.21	1.0		0		0.01		2.21		222
	2	2.57	1.3		0		0.01		2.57		258
	3	3.14	1.9		0		0.01		3.14		315
			M
2 (M)	1	2.73	0.0		0		0.01		2.73		274
	2	2.21	1.0		0		0.01		2.21		222
	3	1.70	3.0		0		0.01		1.70		171
			A(-)
3 (A(-))	1	2.21	− 0.5		0		0.01		2.21		222
	2	2.21	− 1.5		0		0.01		2.21		222
	3	2.21	− 2.0		0		0.01		2.21		222
			B(-)
4 (B(-))	1	2.21	− 0.5		0		0.01		2.21		222
	2	2.21	− 1.5		0		0.01		2.21		222
	3	2.21	− 2.0		0		0.01		2.21		222
			A(+)
5 (A(+))	1	2.21	0.5		0		0.01		2.21		222
	2	2.21	1.5		0		0.01		2.21		222
	3	2.21	2.0		0		0.01		2.21		222
			B(+)
6 (B(+))	1	2.21	0.5		0		0.01		2.21		222
	2	2.21	1.5		0		0.01		2.21		222
	3	2.21	2.0		0		0.01		2.21		222
			Sr	Du
7 (Sr, Du)	1	2.21	0.07	1.0	0		0.01		2.21		222
	2	2.21	0.4	0.175	0		0.01		2.21		222
	3	2.21	1.0	0.07	0		0.01		2.21		222
			K_R
8 (K_R)	1	2.21	0.0		0		0.01		2.21		222
	2	2.21	0.3		0		0.01		2.21		222
	3	2.21	0.9		0		0.01		2.21		222

Figure 4

Architecture of neural networks.

The Figs. 5a–h and 6a–h illustrate the results of case 3 of all the eight scenarios (1–8) showing performance as mean square error (MSE) analysis and the training states detail, respectively. Figure 7a–h explains Fitting graphs for case 3 of scenarios 1–8, respectively. In Fig. 8a–h, the histograms based on error analysis are given, while regression analysis, for each case 3 of all eight scenarios of TFFMFTDECR model, is presented in Fig. 9. The solution graphs are given in Figs. 10p1–p4, 12p5–p8, 14p9–p12 and 16p13–p15, while graphs of absolute error (AE) of the outcomes of the proposed NNs-BLMS and that of the reference solutions are provided in Figs. 11e1–e4, 13e5–e8, 15e9–e12 and 17e13–e15, respectively. Moreover, the Table 4 is provided to illustrate the convergence based on mean square error, epochs, time of execution, performance and the measures of backpropagation scheme for cases 1–3 of scenarios 1–8 for TFFMFTDECR model.

Figure 5

Performance illustration through graphs of C3 of all eight scenarios.

Figure 6

Graphical illustration of Training states of C3 of all eight scenarios.

Figure 7

Fitness graphs of C3 of all eight scenarios.

Figure 8

Illistration via Error histograms(EH) for C3 of all eight scenarios.

Figure 9

Illustrations via regression plots for C3 of all eight scenarios.

Figure 10

Solution plots of , G , F and by using the proposed scheme.

Figure 12

Solution plots of , G , F and by using the proposed scheme.

Figure 14

Solution plots of by using the proposed scheme.

Figure 16

Solution plots of and by using the proposed scheme.

Figure 11

AE plots for outputs by proposed scheme and reference data.

Figure 13

AE plots for outputs by proposed scheme and reference data.

Figure 15

AE plots for outputs by proposed scheme and reference data.

Figure 17

AE plots for outputs by proposed scheme and reference data.

Table 4

Various statistics of scenarios (1–8) for cases (1–3) of TFFMFTDECR model.

Scenario	Case	Mean square error			Performance	Gradient	Mu	Epoch	Time
		Training	Validation	Testing
1	1	3.51E−10	1.58E−9	5.56E−10	3.52E−10	7.06E−8	1.00E−10	96	1
	2	5.16E−10	1.03E−9	6.65E−10	5.16E−10	9.96E−8	1.00E−8	728	17
	3	2.18E−9	4.73E−9	2.14E−8	2.19E−9	9.99E−8	1.00E−8	809	14
2	1	3.49E−11	6.75E−10	1.01E−10	3.5E−8	2.56E−8	1.00E−10	198	3
	2	6.91E−10	7.21E−9	1.06E−8	6.91E−10	9.92E−8	1.00E−11	128	2
	3	9.27E−10	4.93E−9	4.08E−8	9.27E−10	8.74E−8	1.00E−10	81	1
3	1	1.07E−9	5.67E−9	4.59E−9	1.089E−9	8.30E−8	1.00E−10	67	1
	2	2.65E−10	5.35E−10	4.30E−10	7.99E−11	9.47E−8	1.00E−11	38	1
	3	1.25E−9	2.96E−9	6.20E−9	1.25E−9	2.30E−8	1.00E−10	23	1
4	1	7.53E−10	5.92E−9	2.95E−9	7.53E−10	4.41E−8	1.00E−10	86	1
	2	2.59E−10	9.20E−10	1.21E−9	2.07E−10	5.89E−8	1.00E−10	88	1
	3	1.46E−9	2.07E−9	2.44E−9	9.59E−10	2.12E−6	1.00E−11	60	1
5	1	1.29E−9	5.27E−9	4.74E−9	1.30E−9	9.83E−8	1.00E−10	48	1
	2	1.81E−9	6.98E−9	6.67E−9	1.67E−9	9.74E−8	1.00E−10	58	1
	3	4.45E−10	4.08E−9	3.01E−9	3.00E−10	3.20E−7	1.00E−11	131	2
6	1	8.64E−10	2.16E−8	3.30E−9	8.34E−10	1.45E−6	1.00E−9	851	12
	2	2.16E−10	2.80E−10	3.45E−10	9.65E–11	2.78E−6	1.00E−12	44	1
	3	2.54E−9	2.14E−8	1.59E−8	2.55E−9	9.69E−8	1.00E−9	87	1
7	1	1.16E−9	4.49E−9	4.20E−9	9.10E−10	5.39E−7	1.00E−11	63	1
	2	1.21E−9	9.59E−9	5.90E−9	1.22E−9	7.23E−8	1.00E−9	131	2
	3	1.19E−10	3.38E−9	5.45E−10	1.20E−10	7.24E−8	1.00E−10	108	1
8	1	1.74E−9	8.01E−9	1.09E−8	1.66E−9	9.28E−8	1.00E−9	116	1
	2	4.05E−10	2.20E−9	2.40E−9	3.61E−10	2.00E−7	1.00E−11	116	1
	3	1.74E−9	8.01E−9	1.09E−8	1.66E−9	9.28E−8	1.00E−9	116	1

The eight sub-figures (a)–(h) of Fig. 5, gives the convergence of the solution in terms of MSE for train, validation and testing of case 3 of scenarios 1–8 of TFFMFTDECR model. In Table 4, it can be seen that best performance achieved at 809, 81, 23, 60, 131, 87, 108 and 116 epochs with MSE around 10–9–10–8, 10–10–10–8, 10–9, 10–9, 10–10–10–9, 10–9–10–8, 10–10–10–9 and 10–9–10–8 , respectively. Corresponding values of gradients are [9.9 × 10–08, 8.74 × 10–08, 2.30 × 10–08, 2.12 × 10–06, 3.20 × 10–07, 9.69 × 10–08, 7.24 × 10–08 , 9.28 × 10–08], while step size for case 3 of all the eight scenarios is [10–8, 10–10, 10–10, 10–11, 10–11 , 10–9, 10–10, 10–09] and time of execution in seconds for the corresponding case 3 of eight scenarios is [14, 1, 1, 1, 2, 1, 1, 1], which shows the convergence and accuracy of the proposed solver NNs-BLMS.

The states of convergence for case 3 of all eight scenarios are given in subgraphs (a)–(h) of Fig. 6, while the comparison of the results of designed solver NNs-BLMS and the reference solution of the corresponding cases is presented in subgraphs (a)–(h) of Fig. 7. The above mentioned-cases are further studied in the form of error histograms (EH) in subfigures (a)–(h) of Fig. 8. It can be seen that the bins near to desired zero line having values about 5.7 × 10–06, 3.7 × 10–06, 1.1 × 10–05, − 8.5 × 10–06, 5.5 × 10–06, 3.0 × 10–05, − 1.6 × 10–05 and 1.2 × 10–05. The results of regression analysis are provided graphically for the stated cases. The numerical value of correlation coefficient is consistently nearly equal to 1, for testing, validation and train datasets, which justify the worthy performance of designed solver NNs-BLMS to solve the eight scenarios of TFFMFTDECR model.

The variation of physical parameters , , , , , , and represents eight scenarios, where and represents the negative values of parameters and respectively. Similarly and the positive values of parameters and respectively. The solution plots of velocity component in radial direction , in axial direction , azimuthal direction and temperature profile are displayed for the variation of i.e. scenario 1 in subplots (p1)–(p4) of Fig. 10, while there corresponding Absolut Error (AE) plots are presented by subplots (e1)–(e4) of Fig. 11. It is obvious from Fig. 10 that and increases, while and decreases for higher . The rotation parameter is the ratio of stretching to swirling rates. Thus, the progression in the rotation of disc results a higher stretching behavior in comparison with rotation effect, which causes higher radial flow velocity. In case of moving away from the disc surface, the radial velocity becomes lower. As the stretching effect is dominant at the surface, therefore, moving away from the surface results a decrease in the radial velocity. The decreasing trend of axial velocity component is because of fluid particles, which are moving towards radial direction are compensated, due to rotation, by moving downward at the expense of diminishing axial velocity profile. The results are in line with the published work of Mushtaq and Mustafa[37]. It is depicted in Fig. 11 that AE for , , and are 10–7 to 10–3, 10–8 to 10–3, 10–8 to 10–3 and 10–8 to 10–3, respectively. The solution plots of , , and are shown variation of i.e. scenario 2 in subplots (p5)–(p8) of Fig. 12, while there corresponding AE plots are given by subplots (e5)–(e8) of Fig. 13. It can be seen from Fig. 12 that , and decrease, while enhances for higher . A drag-like force is induced by the implementation of magnetic field in the axial direction. The consequences of this force are in the form of dinishing the velocities in all the directions at the cost of enhancement of the temperature profile. The same results are observed by Maleque[38]. Its is notable that higher magnetic parameter in axial direction enforces the fluid particles a higher resistance, reduces the reduction in fluid velocity and enhances temperature profile of the fluid. It is seen in Fig. 13 that AE for , , and are 10–8 to 10–3, 10–8 to 10–3, 10–8 to 10–3 and 10–7 to 10–3, respectively.

The solution plots of are shown variation of , , and i.e. scenarios 3–6 in subplots (p9)–(p12) of Fig. 14, while there corresponding AE plots are given by subplots (e9)–(e12) of Fig. 15. It can be seen from Fig. 14 that decreases for decreasing values of and , while it increases for higher values of and . The temperature/space dependent heat source/sink parameters have considerable effect on the dynamics of temperature profile. The decreasing values of , which represent the heat sink parameters, cause reduction in temperature profile of the fluid. This result is physically acceptable because it implies absorption of heat in the fluid flow, which consequently, results the reduction in the thickness of boundary layer of Maxwell fluid. On the other hand higher values of mean heat generation of space/temperature dependent parameters. The enhancement in these parameters depict the enhancement of the temperature boundary layer thickness, which results mounting in the temperature profile. The results of this study have an agreement with the outcomes obtained by Irfan et al. [31]. It is seen in Fig. 15 that AE for are 10–7 to 10–3, 10–8 to 10–3, 10–7 to 10–3 and 10–8 to 10–3, corresponding to , , and , respectively.

The solution plots of , for increasing Soret number and consequently decreasing Dufour number, while the solution graph of for increasing values of reaction rate parameter are presented in subgraphs (p13)–(p15) of Fig. 16, whereas the corresponding error graphs are provided in subgraphs (e13)–(e15) of Fig. 17. The stated variation of and causes an increase in concentration profile, while a decrease in temperature profile (see Fig. 16p13, p15), whereas the concentration profile diminishes for the stated variation of (see Fig. 16p14). The boosting chemical reaction rate has influence on concentration profile to reduce because of the reason that during chemical reaction, the consumption of chemical species takes place, which causes the reduction in solutal concentration profile. The mass flux caused by a temperature difference termed as Soret effect, while the energy flux due to concentration gradient is called Dufour effect. The higher Soret and diminishing Dufour effects causes the higher concentration boundary layer thickness of the fluid. This is because of the higher mass flux causes an increase in the concentration profile. A declining trend in temperature profile is seen for higher Soret and simultaneous declining Dufour number. This is due to decreasing temperature flux due to less concentration difference. The AE plots Fig. 17e13,e15 have AE 10–7 to 10–3 and 10–6 to 10–3, while Fig. 17e14 have AE of 10–6 to 10–3.

Conclusions

The designed numerical computing paradigm through intelligent neural networks with backpropagation scheme of Levenberg–Marquardt, is applied to solve the Maxwell nanofluid model of thin film flow with thermo-diffusion effect and chemical reaction. The reference dataset for NNs-BLMS was taken by using Adams numerical solver to solve TFFMFTDECR system in case of eight variants based on rotation parameter, magnetic parameter, space dependent heat source/sink, heat dependent heat source/sink, combine effects of Soret and Dufour numbers, and reaction rate parameter. The 70% data was used for train, 15% for validation and 15% for testing processes through NNs-BLMS and matching level of accuracy around 10–11 to 10–3 between the proposed and reference solutions show the precision and consistency. The performance validated through inferences from presented numerical tables, graphical descriptions of convergence plots, MSE, regression illustrations and error-histogram analyses. The outcomes of physical parameters of interest for TFFMFTDECR model revealed that:

The decreasing trend of axial velocity component is because of fluid particles, which are moving towards radial direction are compensated, due to rotation, by moving downward at the expense of diminishing axial velocity profile.

The higher rotation parameter enhances fluid film thickness and diminishes the temperature profile.

Higher magnetic field diminishes liquid film thickness, while enhances temperature profile.

Higher Soret number and corresponding diminishing Dufour number enhances the concentration profile, while declining trend in temperature profile is seen for higher Soret and simultaneous declining Dufour number. This is due to decreasing temperature flux due to less concentration difference.

Higher reaction rate parameter causes decreasing concentration profile.