Next Article in Journal
High-Speed Erosion Behavior of Hydrophobic Micro/Nanostructured Titanium Surfaces
Next Article in Special Issue
The Theoretical Study of Kink Deformation in Graphite Based on Differential Geometric Method
Previous Article in Journal
An Operational DNA Strand Displacement Encryption Approach
Previous Article in Special Issue
Virtual Vibrational Spectrometry of Stable Radicals—Necklaced Graphene Molecules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Quantitative Features Analysis of Water Carrying Nanoparticles of Alumina over a Uniform Surface

by
Abdul Hamid Ganie
1,*,
Fazlullah Fazal
2,
Carlos Andrés Tavera Romero
3 and
Muhammad Sulaiman
2,*
1
Basic Science Department, College of Science and Theoretical Studies, Saudi Electronic University, Abha Male 61421, Saudi Arabia
2
Department of Mathematics, Abdul Wali Khan University, Mardan 23200, Pakistan
3
COMBA R&D Laboratory, Faculty of Engineering, Universidad Santiago de Cali, Cali 76001, Colombia
*
Authors to whom correspondence should be addressed.
Nanomaterials 2022, 12(5), 878; https://doi.org/10.3390/nano12050878
Submission received: 23 December 2021 / Revised: 23 February 2022 / Accepted: 1 March 2022 / Published: 6 March 2022
(This article belongs to the Special Issue Computational Study of Nanomaterials)

Abstract

:
Little is known about the rising impacts of Coriolis force and volume fraction of nanoparticles in industrial, mechanical, and biological domains, with an emphasis on water conveying 47 nm nanoparticles of alumina nanoparticles. We explored the impact of the volume fraction and rotation parameter on water conveying 47 nm of alumina nanoparticles across a uniform surface in this study. The Levenberg–Marquardt backpropagated neural network (LMB-NN) architecture was used to examine the transport phenomena of 47 nm conveying nanoparticles. The partial differential equations (PDEs) are converted into a system of Ordinary Differential Equations (ODEs). To assess our soft-computing process, we used the RK4 method to acquire reference solutions. The problem is investigated using two situations, each with three sub-cases for the change of the rotation parameter K and the volume fraction ϕ . Our simulation results are compared to the reference solutions. It has been proven that our technique is superior to the current state-of-the-art. For further explanation, error histograms, regression graphs, and fitness values are graphically displayed.

1. Introduction

The Coriolis force causes gases and liquids in the northern hemisphere to bend to the right. The Coriolis force plays a significant role in fundamental flow equations as an inertial force, magnetohydrodynamic force, and viscous force. Each fluid that moves on the earth’s surface is subjected to gravitational force, pressure gradient force, centrifugal force, and frictional force. While the Coriolis force does not affect every flow, it does have the power to shift the direction of transportation phenomena on the earth’s surface in the atmosphere and ocean. As a result, accepting that the Coriolis force has no influence on any non-static transport phenomena on the earth’s surface is impractical.
In addition to the Coriolis effects, the earth’s rotation has a major role in global climate change. According to [1,2], the Coriolis force may be used to characterize the whirling directions of typhoons, hurricanes, strong cyclonic storms, and strong cyclonic storm surges. In [3], researchers present an analysis of layer flow with an integral boundary of an incompressible, constant momentum in three dimensions, with a focus on the importance of rotation as it applies to the blades of a wind turbine. One of the findings was that, in addition to centrifugal forces, Coriolis forces also play a role in transport phenomena. The radial flow observed in the boundary layer is caused by centrifugal acceleration, whereas the flow in the clockwise direction is caused by the Coriolis force.
The Coriolis force is employed in businesses such as photobioreactors, sewage treatment, and bioreactors for particular tissues; see [1,2]. In [3] researchers provide a three-dimensional analysis of an incompressible constant momentum applied to wind turbine blades, with an emphasis on the relevance of rotation. The Coriolis force and the centrifugal force are both responsible for radial flow in the boundary layer. On the other hand, only the Coriolis force is responsible for flow in the clockwise direction. The Coriolis effects on the dynamics in the centrifugal bioreactor were investigated by [4], and it was observed that as the rotation of the flow increases, vortices develop and that the dynamics are regulated by Coriolis force. Coriolis forces strongly impact smaller microbial cultures because they become locked in vortices. On the other hand, large cells have larger rates of sedimentation, resulting in streamline deflection and continuing Coriolis pressures. However, this results in a stably fluidized bed of cells. In the field of atmospheric science, researchers are experimenting with several forms of nanoparticles to remediate contaminated waterways [5]. The influence of various fluids with different thicknesses on the dynamics of the Coriolis force on an object, such as the sharp edge of an airplane or rotating automobiles, was recently explored by [6]. Internal shear stress may be reduced by increasing the flow’s maximum velocity, whereas buoyant forces can raise the flow’s total velocity. The shear stress will be calculated less as the rotation across the flow increases.
Term nanofluid is formed by [7] to designate fluid substances suitably combined along with particles whose diameter falls within the nanoscale range. Particles with a diameter of less than 100 nm have greater thermal, chemical, mechanical, optical, electrical, and magnetic capabilities than regular solids due to their large ratio of surface-area-to-volume. Industries confront challenges such as cooling concerns and product maintenance due to rising heat generation. However, scientists have come to appreciate the thermal and particular nature properties of various fluids generated by solid particles addition (inn micrometer or millimeter scales) to protect energy and save procedure time.
Nanofluids are employed in a wide range of applications in industries, comprising coolants for car engines, cancer therapies, delivery of nano-drugs, syphilis testing, and detergent having nanofluid [8]. The effect of size of the particles on the dynamic viscosity of water–47 nm Al 2 O 3 , water–29 nm CuO and water–36 nm Al 2 O 3 is given by [9]. When particle volume fractions are less than 4%, then the viscosities of both nanofluids containing Al 2 O 3 are equivalent. There are various small or no facts about heat sinks’ influence on the nanofluid movement in thermal management. Moreover, how to retain the temperature at the heat sink’s base as low as possible while the heat transfer rate increases are discussed in [10].
In comparison to nanofluid CuO–water and distilled water, results demonstrate that nanofluid Al 2 O 3 –water has a higher transfer rate of heat. Furthermore, the heat sink can reduce generated heat by 89.6% between the mini-channels. The dynamics of water transporting titania copper and nanoparticles of alumina through a cylinder that is shaped like a wave is investigated by [11]. Volume fraction at all levels, this was revealed that friction at the wall and the local skin friction are negligibly proportional during the motion of water transporting alumina nanoparticles. Nanoparticles addition to the wall can dramatically limit the rate of heat transfer; see [12]. The rising property of random movement of particles having nano-sized mixed thoroughly with base fluid is the distribution of temperature through the nanofluid dynamics; see [13]. In [14,15], the effect of increasing Coriolis force of non-Newtonian Casson fluid on dynamics in the absence and presence of Lorentz force was investigated. The impacts on dynamics of water transporting 47 nm nanoparticles of alumina, Coriolis force, volume fraction, and heat source/sink over uniform surfaces have not been studied. Following up the literature on the findings mentioned above, it is critical to find scientifically valid solutions to the following questions.
  • What effect does Coriolis force have on layer flow at the boundary of a nanofluid alumina-water subjected to heat?
  • How does the growing volume percentage of nanoparticles impact the wall friction, as well as the mass and heat transfer rates in a nanofluid of alumina–water?
  • How does the volume percentage of 47 nm alumina nanoparticles affect heat mass transfer rate and local skin friction?
  • How can numerical solutions to the model be calculated with accuracy and efficiency. Nowadays, in the modern world, we need a method that requires less effort while giving good efficiency.
  • How can we get a simple and easy method which gives accurate results?
  • How can we save our time?
The rest of the manuscript is organized such that in Section 2 a detailed problem formulation is elaborated. In Section 3, the solution methodology is presented. Results and discussions are given in Section 4. At the end of this manuscript, we present the conclusions of this study.

2. Formulation of the Principal Equations

The derivation is presented for transport phenomena principal equations and the inertial forces of a frame of reference that is rotating.

2.1. Forces for a Frame of Reference Which Is Rotating

The planned path and specious route are two visual pathways in a frame of reference with rotation. In the true/intended path, the trajectory is effortlessly warped by a fake force recognized as the Coriolis force, resulting in the specious path. Suppose the angular velocity Ω , position vector r = ( r 1 , r 2 , r 3 ) and in the rotating frame the absolute time derivative:
D r D t = d r d t + Ω × r , d r d t = ( d r 1 d t , d r 2 d t , d r 3 d t ) .
According to the second law of motion of Newton, force is:
F = m a r d Ω d t m × r 2 m Ω × d r d t m Ω × ( r × Ω ) .
Euler force are the three acting principal forces on a body travelling in such a frame which is rotating in Equations (1) and (2). F e = m d Ω d t × r , Coriolis force F c = 2 m Ω × d r d t , and Centrifugal force F g = m Ω × ( r × Ω ) . Azimuthal force (also known as Euler’s force) is required to keep the body in orbit because it acts parallel to but against the angular velocity. The centrifugal force acts radially outward from the axis of rotation of the rotating frame. The angular velocity and the Coriolis force are perpendiculars, and the Coriolis force acts outward. After combining the Navier–Stokes equation published by [4,6], the magnitude order was executed on d r d t for boundary layer flow. The buoyancy and Coriolis forces on the flow of fluid in a uniformly thick horizontal surface in x and z axes are referred to as the body force term.
f b x = g β ( T T ) + g β * ( C C ) 2 Ω w ,
f b z = g β ( T T ) + g β * ( C C ) + 2 Ω u .
The rotating frame’s angular velocity is Ω , volumetric thermal expansion is β , gravity is g, and  β is the nanofluid’s thermal expansion coefficient. Given that the length of characteristics is measured in meters, the unit of the body forces fbx and fbz is m s 2 .

2.2. Mathematical Formulation

Alumina–water nanofluid dynamics on a rotating horizontal surface at a fixed time with varying thicknesses are investigated. Figure 1 shows how the surface rotates counter-clockwise. The dynamics of an alumina–water nanofluid on a rotating horizontal surface at a fixed time and with varying thicknesses are investigated [16]. Figure 1 shows how the surface rotates counter-clockwise. As a result, regulating of flow by boundary layer equation in the presence of Coriolis force are shown in Figure 1 is:
u x + v y = 0 ,
u u x + v u y = μ n f ρ n f 2 u y 2 + ( T T ) g β + ( C C ) g β * 2 Ω w ,
u w x + v w y = μ n f ρ n f 2 w y 2 + g β ( T T ) + g β * ( C C ) + 2 Ω u ,
u T x + v T y = α n f 2 T y 2 + τ [ D B Δ C C y T y + D T T T y T y ] + Q o [ T w ( x ) T ] ( ρ C p ) n f Exp [ n y a ϑ b f ] ,
u C x + v C y = D B 2 C y 2 + D T Δ C T 2 T y 2 ,
the following are the boundary conditions:
v = 0 , u = a x , w = 0 , C w = C , T = T w , at y = 0 ,
u 0 , w 0 , T = T , C = C when y .
The components of velocity are u in the x-direction and v in the y-direction, the concentration of expansion coefficient is β * , free stream temperature T , nano fluid temprature T, nanofluid kinematic viscosity is ϑ n f , the temperature of wall is T w , nanofluid concentration is C, the concentration at free stream C , wall concentration C w , the volume fraction of nanoparticle ϕ , coefficient of thermophoretic diffusion D T , coefficient of Brownian diffusion D B , angular velocity Ω , nanofluid thermal diffusivity is α n f , the thermal conductivity of nanofluid k n f , nanofluid specific heat capacity ( ρ c ρ ) n f , and the ratio of specific heat capacity of the base fluid to the specific heat capacity of the nanoparticle is τ = ( ρ c p ) n p ( ρ c p ) b f . Following [17,18,19,20], the density, thermal conductivity, similarity variables, and specific heat capacity of the nanofluid are specified as:
α n f = κ n f ( ρ c p ) n f , κ n f = κ b f [ κ n p + 2 κ b f 2 ϕ ( κ b f κ n p ) κ n p + 2 κ b f + ϕ ( κ b f κ n p ) ] , ρ n f = ( 1 ϕ ) ρ b f + ϕ ρ n p , ( ρ c p ) n f = ( 1 ϕ ) ( ρ c p ) b f + ϕ ( ρ c p ) n p , μ n f ρ n f ν b f = 0.904 e 0.148 ϕ 1 ϕ + ρ n p ρ b f ϕ , v = ( a ϑ b f ) 1 / 2 F ( t ) , u = a x d F d t , w = axH ( t ) , t = y a ϑ b f , Φ = C C C w C Θ = T T T w T .
Obtaining the dimensionless equations as:
μ n f ρ n f ν b f F + F F F F K H + G r t Θ + G r n p Φ = 0 ,
μ n f ρ n f ν b f H + H F H F + K F + G r t Θ + G r n p Φ = 0 ,
[ k n p + 2 k b f 2 ϕ ( k b f k n p ) k n p + 2 k b f + ϕ ( k b f k n p ) ] ( 1 ϕ ) + ϕ ( ρ c p ) n p ( ρ c p ) b f Θ + N b Θ Φ + N t Θ Θ + PrF Θ + Pr γ Exp ( n η ) ( 1 ϕ + ϕ ( ρ C p ) s p ( ρ C p ) b f ) = 0 ,
Φ + N t N b Θ + S c F Φ = 0 ,
where
γ = Q o ( ρ C p ) b f a , K = 2 Ω a , G r t = g β ( T w T ) a 2 x , P r = ϑ b f α b f , N b = τ D B α b f , N t = τ D T ( T b f T ) T α b f , G r n p = g β * ( C w C ) a 2 x , S c = ϑ b f D B ,
are the thermal Grashof number, rotation parameter, Brownian parameter, Schmidt number, thermophoretic parameter, nanoparticle Grashof number, and Prandtl number. The following are the non-dimensionalized boundary conditions:
F = 0 , F = 1 , H = 0 , Θ = 1 , Φ = 1 at t = 0 ,
F 0 , H 0 , Θ 0 , Φ 0 , as t .
The thermophysical properties of the water and 47 nm alumina nanoparticles are essential to remember, which are given as:
k n p = 40 W m 1 K 1 , k b f = 0.613 W m 1 K 1 , ( ρ c p ) n p = 3,037,050 J K 1 m 3 , ( ρ c p ) b f = 4,166,880.90 J K 1 m 3 , ρ n p = 3970 K g m 3 , and ρ b f = 997.10 K g m 3
see Nehad et al [21]. First, we get the solutions of the problem by Mathematica using Rk4 command, and then we apply neural on that data and get the outputs of the neural using
G r t = 1 , G r n p = 1 , N b = N t = 0.1 , n = 1 , S c = 0.62 .
Through transformation first, we obtained the Ordinary Differential Equations (ODEs) system of the considered problem and solved the system. In this paper, we solve the system of obtained Ordinary Differential Equations (ODEs) and then get the solution of the system at 1001 different points. Then, we discuss the effect of two other parameters i.e., rotation parameter K and the volume fraction ϕ , on the given problem. We also solve the system through neural networks and compare the values with the reference solution using different values of different parameters. We discuss the effects of rotation and volume fraction parameters on the system using a neural network. The heat transfer rate is proportional to the Nusselt number, and the mass transfer rate of the nanofluid of the problem is proportional to the Sherwood number.

3. Solution Methodology

We use artificial neural networks in this article due to the important individual qualities of this method. Artificial neural networks avoid singularity; if there is any singularity present in the equations, the neural networks avoid it and give good results [22,23]. A neural network provides a general solution to the problem, due to which we can find values of the problem outside the interval as well [24]. The significance of the neural network approach to mainly non-linear modeling capability, likely interpolations, forecasting, extrapolations, ill-conditioning, robustness to noise and insufficient data, and having ease of use, are discussed in [25]. It also has good fault tolerance ability [26]. NN can train the machine itself; that is, ANNs learn from the events and apply those properties on similar events [27]. NN has the property to perform different functions at the same time [28]. Neural networks (NN) are deep learning algorithms. Neural networks are machine learning networks, and the human brain inspires these networks. Neural networks contain one input layer, hidden neurons, and one or more than one output layer. The hidden neurons are connected and pass the signals between one another. The artificial neural network (ANN) is an algorithm used for non-linear training. It is essential due to its fast converging speed and accurate results. Performance achieved by this method is outstanding.
The Levenberg–Marquardt neural network consists of structure layer, the number of hidden neurons, network topology, and selection of input data and output data arbitrary for training, validation, and testing samples. All these things need to be determined for the solution of a problem. In this system of ODEs, we used 1001 points obtained by using the “NDSolve” command in Mathematica by using the Runge-Kutta order-4 method. We set the step size 0.004 for getting the 1001 points between 0 and 4. Then, we use those data set points in the neural toolbox by using the Matlab command window and suggest the LMB-NN method for the best result, and then we set input and output in LMB-NN. We choose the following set of data in this problem.
  • Data for training is 90 percent.
  • Data for testing is 5 percent.
  • Data for validation is 5 percent.
Choosing as in Figure 2 the number of hidden neurons is arbitrary, and we choose in our problem the number of hidden neurons is 40, which gives us good results.
Figure 2 shows us the basic working phenomenon of a neural network that how ANN works. We see that the neural network consists of inputs and the weights assigned to these inputs, then the hidden neuron, activation function, and in the last, it gives us the outputs.
The model we obtained for our problem of the neural networks using the nftool command in the Matlab window is shown in Figure 3. This problem contains 40 neurons. You may change the number of neurons, but the 40 neurons give us the best results in this problem. That is why we use 40 neurons in our problem. The inputs and outputs can also be clearly seen from Figure 3.
First, the partial differential equations are transformed into the system of ordinary differential equations. Then, one must solve those systems of ordinary differential equations in Mathematica using RK4 technique and obtain the system’s solution on 1001 different points by using 0.004 step size. Then, this data at 1001 points are copied to an Excel sheet, and from the MS Excel sheet, the data are transferred to Matlab for using the neural networks. In the neural networks, LMB-NN chooses for the solution of the system. In this system, there are four input values, four hidden layers, and four output values, as shown in Figure 3. The number of hidden neurons chosen is 40. Ninety percent values are selected for training, five percent for testing, and five percent is for validation. Run the command for several iterations. If the output values obtained are excellent and accurate, the process is ended; otherwise, it is necessary to retrain the network and run again until it reaches the minimum mean square error value (MSE). Depending which iteration the network gives better outputs, that will be the stopping iteration.

4. Results and Discussion

This study examined the effects of volume fraction, heat sink, and Coriolis force on the behavior of water conveying 47 nm alumina nanoparticles over a uniform surface using a soft computing technique. To generate nanofluids, researchers used Al 2 O 3 nanoparticles ranging in size from 13 nm to 302 nm, with a 2 percent to 36 percent increase in thermal conductivity (see [29]). As a result of ever-increasing heat production, industries face cooling issues and product maintenance obstacles. However, scientists appreciate the precise nature and thermal properties of various fluids generated by adding solid particles (on the micrometer and millimeter scales) to reduce energy consumption and processing time. Nanofluids are employed in a variety of applications, including automotive engine coolants, cancer therapies, nano-drug delivery, syphilis diagnosis, and detergent with nanofluids (see [8]). Particles with less than 100 nanometers diameter have better mechanical, thermal, optical, magnetic, chemical, and electrical properties than ordinary solids (see [16]). Alumina nanoparticles with a diameter of 47 nanometers are suitable in the above applications because they have a substantial surface-area-to-volume ratio. Furthermore, Ali et al. (see [10]) discussed how to keep the temperature at the heat sink’s base as low as possible as the heat transfer rate increases. In comparison to CuO–water nanofluid and distilled water situations, the results demonstrate that Al 2 O 3 –water nanofluid has a higher heat transfer rate. Furthermore, the heat sink can reduce generated heat by 89.6% between the mini-channels. At all levels of volume fraction, it was revealed that local skin friction proportional to friction at the wall is negligible during the motion of water transporting alumina nanoparticles. The effect of particle size on the dynamic viscosity of Nguyen et al. [9] described the effect of particle size on the dynamic viscosity of:
(a)
Water–36 nm Al 2 O 3 ;
(b)
Water–47 nm Al 2 O 3 ; and
(c)
Water–29 nm CuO nanoparticles experimentally.
When the particle volume fraction is less than 4%, the viscosities of both nanofluids with Al 2 O 3 are similar, according to the findings. Water has a viscosity of 47 nm above this region. The viscosity of Al 2 O 3 nanofluid is significantly higher than that of water–36 nm Al 2 O 3 nanofluid. As nanoparticles, Al 2 O 3 and TiO 2 were employed, using thermal oil as the base fluid. The needed amount of nanoparticles and base fluid were mixed. The diameter of alumina Al 2 O 3 nanoparticles in spherical shape ranged from 5 nm to 250 nm, with a mean diameter of 47 nm according to the manufacturer (Sigma-Aldrich Co., Burlington, MA, USA) (see [30]).
This means that nanoparticles with a higher viscosity have better thermal, electrical, chemical, mechanical, magnetic, and optical capabilities. The viscosity of Al 2 O 3 nanofluid is significantly higher than water–36 nm. Here, the particle volume fraction is more significant than 4%, why is why we are solely interested in studying 47 nm alumina particles. Additionally, Al 2 O 3 and TiO 2 were employed as nanoparticles, with thermal oil serving as the base fluid. The needed amount of nanoparticles and base fluid were mixed. The diameter of alumina Al 2 O 3 nanoparticles in spherical shape ranged from 5 nm to 250 nm, with a mean diameter of 47 nm, according to the manufacturer (Sigma-Aldrich Co.) (see [30]). This is also the reason that we are talking only about the 47 nm alumina nanoparticles in our study.
The method we used in this paper and all the results obtained are given in this section. All 1001 points were obtained using the “NDSolve” command in Mathematica using the Runge-Kutta order-4 method. We set the step size 0.004 for receiving the 1001 points between 0 and 4. Then, we used those data set points in the neural toolbox using Matlab and suggest the LMB-NN method for the best result.
Figure 4 shows us the basic working phenomenon of a neural network that how ANN works. We see that the neural network consists of inputs and the weights assigned to these inputs, then the hidden neuron, activation function, and, lastly, the outputs.

Tabular and Statistical Analysis

Two scenarios are discussed in this paper; see Table 1. In scenario number 1, we change the value of rotation parameter K (0.1, 0.2, 0.3) and discuss the behavior of the solution of the system of ODEs and train a neural network for each case individually. Similarly, in scenario number 2, we change the value of volume fraction ϕ (0.1, 0.2, 0.3) and observe the behavior of the solution of the given system of ODEs and then train the neural network for each case and compare the results with RK4 data obtained by using the NDSolve command in Mathematica. We discussed two effects: the rotation parameter K and the volume fraction ϕ on the water containing 47 nm alumina nanoparticles. In Algorithm 1, all the procedures of the LMBNN are given in detail.
In Table 2 we present the numerical results obtained from the solution of LMB-NN. We discuss all the cases of both scenarios and observe the time taken to get the mean square error of the training, testing, and validation samples. Also shown are the performance, gradient, Mu values, and the epoch at which we got the good values. Table 3, Table 4, Table 5 and Table 6 show the comparison of RK4 and LMB-NN in all the three cases of scenario 1 in the solution of all the ODEs i.e., F, H, Θ , and Φ , respectively. As mentioned above, the first column in all the tables present the input value. The second column shows the results obtained by RK4, while the third column shows the results obtained through LMB-NN when K = 0.1. Similarly, the fourth and sixth columns represent the values of the results of RK4, while the fifth and seventh columns show the importance of LMB-NN when the value of K = 0.2 and K = 0.3, respectively. We can observe the data set of the methods, and it is clear that there is a tiny error between the data set of both approaches.
Algorithm 1All the process of LMB-NN is given in the pseudocode
  • Starting of LMB-NN
  • Step 1: Construction
  • Construct input and data set
  • Step 2: Selection of Data
  • Target data and input data is chosen in non-linear form i.e., matrix form.
  • Step 3: Startup
  • Startup the ratio of Neuron numbers, testing, validation and training.
  • ▸ 90 Percent is for training
  • ▸ 5 percent is for validation
  • ▸ 5 percent is for testing
  • ▸ Number of hidden neurons is 40
  • ▸ Number of hidden layers is 4
  • Step 4: Weights for training
  • The selected data is trained from the activation function in LMB-NN
  • Step 5: Stopping criteria
  • Step 4 will stop automatically if the following conditions are satisfied.
  • ⋇ Reaching Mu to the maximum value
  • ⋇ Performance value reaches to minimum
  • ⋇ Maximum number of epoch achieved
  • ⋇ Performance of validation is less than maximum fail
  • ⋇ Gradient of performance less than minimum gradient
  • Testing data help us determine that the network is generalized. If the outputs are good and useful forward to step 7, and if the outputs are not desirable, retrain the network.
  • Step 6: Retraining
  • For retraining the hidden neurons, the ratio of testing, training, and validation is changed. Then move again to step 4 and do the same procedure.
  • Phase 7: Output saving
  • The process is ended by saving the output simulation of data statistically as well as numerically.
  • Ending of LMB-NN
Figure 5, Figure 6, Figure 7 and Figure 8 show the solution and absolute error graphs of the variation of rotation parameter K of F, H, Θ , and Φ , respectively. In these graphs, we see for the given system of equations of ODEs, when we change the value of rotation parameter K, how the solution of the system of ODEs changes. The error graph shows the error between the RK4 solution and ANN solutions. The error between both methods is minimal, and both methods’ graphs overlap. The error graphs show us the residual error in the solution of both methods.
Figure 5a shows the numerical solutions of F of all the cases of scenario 1, while Figure 5b–d show the residual errors of F in case 1, case 2, and case 3 of scenario 1, respectively in the solutions of RK4 and LMB-NN. Figure 6a shows the numerical solutions of H of all the cases of scenario 1, while Figure 6b–d show the residual errors of H in case 1, case 2, and case 3 of scenario 1, respectively in the solutions of RK4 and LMB-NN. Figure 7a shows the numerical solutions of Θ of all the cases of scenario 1, while Figure 7b–d show the residual errors of Θ in case 1, case 2, and case 3 of scenario 1, respectively in the solutions of RK4 and LMB-NN. Figure 8a shows the numerical solutions of Φ of all the cases of scenario 1, while Figure 8b–d show the residual errors of Φ case 1, case 2, and case 3 of scenario 1, respectively in the solutions of RK4 and LMB-NN.
Scenario 1 is the variation of rotation parameter (K). We change the value of K and apply ANN for every individual value of the rotation parameter (K). In each case, we obtain the mean square error graphs, training graphs, error histogram graphs, regression graphs, and fitness graphs of the system of ODEs. We note how the values of the above change when we change the rotation parameter (K). The effect of the variation of K can be observed clearly from the values shown in the graphs. We train a neural network for each case by using the “nftool” command in the Matlab environment, setting the input and output data for the network, and using the LMB-NN method to find the best solution for the present system of ODEs. In this problem, we have one input and four hidden output layers. Moreover, we used 40 hidden neurons for a good result.
The results obtained from LMB-NN of case 1, 2, and 3 of scenario 1 are shown in Figure 9, Figure 10 and Figure 11 respectively. These figures contain performance, training, error histogram, regression, and fitness graphs. Figure 9, Figure 10 and Figure 11 show the solutions obtained by artificial neural networks (ANNs). In these figures, we show the performance, training, error histogram, regression, and fitness graphs of the ANNs of all the cases of scenario 1. We use the input data of RK4 from Mathematica by using the “NDSolve” command. Figure 9a shows the performance of case 1 of scenario 1, in which we show the mean square error. Figure 9b shows the training samples of the ANNs of case 1 of scenario 1. Figure 9c shows the error histogram graph of case 1 of scenario 1. The regression of case 1 of scenario 1 is shown in Figure 9d and the fitness graphs of case 1 scenario 1 are shown in Figure 9e. Similarly, we show the performance, training, error histogram, regression, and the fitness figures in subfigures (a–e) of the remaining two cases of scenario 1 in Figure 10 and Figure 11, respectively.
When the value of MSE of a system is smaller, the system is stable. Among all the above cases discussed in the paper, the MSE of scenario 1 case 1 is smaller than all the other cases. The mean square error (MSE) of all the cases discussed in the paper ranges from 10 11 to 10 12 . The Mu and gradient show a better rate of convergence. Auto-correlation shows the relationship between two variables. The 10 11 and 10 12 value of Mu shows us the better convergence. The histogram indicates the reliability of a technique, and the regression of all the above-discussed cases gives good results. The linear relationship between target data and output data are shown by regression. The data gives an accurate solution and is well trained investigated by fitness plots.
We obtained the RK4 outputs from the Mathematica window and then applied the LMBNN method in the neural network to get the outcomes of the technique. The errors between the results of both the methods are given in the Table 7, Table 8, Table 9 and Table 10. The residual error in the variation of the rotation parameter K in all the four ODEs can be easily observed from the tables. Column 1 of the above tables represents the input values, while columns 2, 3, and 4 illustrate the residual error present between the methods RK4 and LMB-NN obtained from Matlab when the rotation parameter K varies. The error in the tables is minimal, which shows that this technique has greater accuracy; moreover, the method is so simple and easy to implement on many problems of various fields and obtains fruitful results.
Table 11, Table 12, Table 13 and Table 14 show us the solution of all the ODEs (F, H, Θ , Φ ) of the system of RK4 and LMB-NN methods of the variation of the volume fraction ϕ respectively. The comparison of RK4 and LMB-NN outputs are shown. The first column of Table 11, Table 12, Table 13 and Table 14 shows the input values. The second column shows the results obtained by the RK4 method, while the third column shows the results generated by LMB-NN at different input values. Similarly, the fourth and sixth columns represent the results generated by RK4, and the fifth and seventh columns of the tables show the results generated by LMB-NN at different input points, respectively. For each ODE, a separate table is given for comparison. All three cases of both scenarios are present in the tables.
Figure 12, Figure 13, Figure 14 and Figure 15 show the solution of the variation of volume fraction ϕ of all the ODEs along with the figures of absolute error. The effect of the volume fraction ϕ on F is given in Figure 12a. The errors in the variation of volume fraction ϕ (i.e., 0.1, 0.2, 0.3) are given in Figure 12b–d, respectively. Similarly, the effect of volume fraction ϕ on H is shown in Figure 13a and the absolute errors in the variation of volume fraction ϕ are shown in Figure 13b–d. In the same manner, the effect of ϕ and errors on the rest of the ODEs are shown in Figure 14 and Figure 15.
Figure 16, Figure 17 and Figure 18 show the solutions obtained by artificial neural networks (ANNs). In these graphs, we show the performance, training, error histogram, regression, and fitness graphs of the ANNs of all the cases of scenario 2. We use the input data of RK4 from Mathematica by using the “NDSolve” command. Figure 16a shows the performance of case 1 of scenario 2, in which we show the mean square error. Figure 16b shows the training samples of the ANNs of case 1 of scenario 2. Figure 16c shows the error histogram graph of case 1 of scenario 2. The regression of case of scenario 2 is shown in Figure 16d and the fitness graphs of case 1 scenario 2 are shown in Figure 16e. Similarly, we show the performance, training, error histogram, regression, and the fitness graphs in subfigures (a–e) of the remaining two cases of scenario 2 in Figure 17 and Figure 18, respectively.
In this paper, first, the outputs are obtained from Mathematica using the RK4 command, and then these outputs are used, and new outcomes are obtained through Matlab by using the “nftool” command for both the parameters. The error in both the outputs of all the four ODEs of the system when ϕ varies is given in Table 15, Table 16, Table 17 and Table 18. These tables show the residual errors in the outputs of both methods when different values of volume fraction ϕ are used. The first column of Table 15, Table 16, Table 17 and Table 18 represents the input values. The second column of the tables represents the residual error in both the methods when ϕ = 0.1. Similarly, the last two columns of the table represent the residual error in both methods when ϕ = 0.2 and ϕ = 0.3, respectively. The residual error is obtained from Matlab.
In our problem, we compare the solutions of artificial neural networks with the RK4 method, and it is clear from all the analyses that the solution of ANNs are entirely overlapping with the solutions of the RK4 method. Our method is best because we used the ANNs method having four hidden output layers, and the result is good to fast and accurate. Another thing is that our method is a computer-based method and an advanced method in which the chances of mistakes in the calculation are about negligible, and we can completely trust the results and require less time, and provide good results. The results of both scenarios’ cases are accurate, and significantly fewer errors mean that this method is very consistent. Moreover, the method is less time-consuming. This method may be used for various problems in the various fields to obtain accurate, fruitful, and valuable results, such as in the water present over a rotating surface in medicine, biodiesel, and maintaining cooling.

5. Conclusions

This study considers the industrial uses of nanofluids, knowing the effects of different parameters on the fluid, and how other parameters affect the heat transfer rate and the mass transfer rate in the flow of a nanofluid. The uses of 47 nm nanoparticles in the water present over a rotating surface can include medicine, biodiesel, maintaining cooling in different industries, etc. Variations of these nanoparticles still need to be investigated. The LMB-NN is a machine learning-based technique that is fast converging and accurate. All the results given by this method for all cases of both scenarios have minimal residual errors, which means that this method is consistent and reliable.
We have calculated solutions of problems under consideration with less mean squared errors (MSE); thus, it is established that our methodology is stable. In all test cases, which we have discussed in the paper, the MSE of scenario 1 case 1 is smaller than all the other cases. The mean square error (MSE) of all the cases discussed in the paper ranges from 10 11 to 10 12 . The Mu and gradient show a better rate of convergence. Auto-correlation shows the relationship between two variables and the 10 11 to 10 12 value of Mu shows us the better convergence performance of our proposed strategy. The reliability of our technique is indicated by the histogram and regression plots, which show promising results in all scenarios. The linear relationship between target data and output data showed by regression plots is a further strength of our soft computing technique. Our solution technique gives an accurate solution, and fitness plots are presented for better clarity of our claims.
We have discussed variations in two parameters i.e., rotation parameter K and volume fraction ϕ , and compared our results with RK4. Our solutions are in strong agreement with the reference points. Our technique helps solve problems with unknown landscapes and is a robust procedure. From the outputs of our simulations, we see that it is very straightforward, easy, and simple to implement, and the method is computer-based and the outputs are matching with the reference inputs. The mean square error in the outputs of this method is very small; see Table 2. The mean square error range in this paper is 10 11 to 10 12 . In the future, the new variants of the artificial intelligent networks based on the integrated intelligence will be used for many problems in different fields to obtain highly accurate and good results.

Author Contributions

Data curation, M.S. and F.F.; Formal analysis, M.S. and F.F.; Funding acquisition, C.A.T.R. and A.H.G.; Investigation, M.S. and F.F.; Methodology, M.S. and F.F.; Project administration, M.S.; Resources, M.S., C.A.T.R. and A.H.G.; Software, M.S.; Supervision, M.S.; Visualization, F.F.; Writing—original draft preparation, F.F. and M.S.; Writing—review and editing, M.S., C.A.T.R. and A.H.G. All authors have read and agreed to the published version of the manuscript.

Funding

The APC was funded by Dirección General de Investigaciones of Universidad Santiago de Cali under call No. 01-2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

This research has been funded by Dirección General de Investigaciones of Universidad Santiago de Cali under call No. 01-2021.

Conflicts of Interest

The authors declare that they have no competing interest.

Abbreviations

The following abbreviations are used in this manuscript:
ANNArtificial neural network
RK4Runge-kutta order four technique
LMBNNLevenberg–Marquardt backpropagation neural network
mseMean squared error

References

  1. Debnath, L. On unsteady magnetohydrodynamic boundary layers in a rotating flow. ZAMM-J. Appl. Math. Mech. Für Angew. Math. Und Mech. 1972, 52, 623–626. [Google Scholar] [CrossRef]
  2. Deng, L.; Li, T.; Bi, M.; Liu, J.; Peng, M. Dependence of tropical cyclone development on Coriolis parameter: A theoretical model. Dyn. Atmos. Ocean. 2018, 81, 51–62. [Google Scholar] [CrossRef]
  3. Du, Z.; Selig, M. The effect of rotation on the boundary layer of a wind turbine blade. Renew. Energy 2000, 20, 167–181. [Google Scholar] [CrossRef]
  4. Detzel, C.J.; Thorson, M.R.; Van Wie, B.J.; Ivory, C.F. A study of the Coriolis effect on the fluid flow profile in a centrifugal bioreactor. Biotechnol. Prog. 2009, 25, 1025–1034. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Punia, P.; Bharti, M.K.; Chalia, S.; Dhar, R.; Ravelo, B.; Thakur, P.; Thakur, A. Recent advances in synthesis, characterization, and applications of nanoparticles for contaminated water treatment—A review. Ceram. Int. 2021, 47, 1526–1550. [Google Scholar] [CrossRef]
  6. Koriko, O.; Adegbie, K.; Oke, A.; Animasaun, I. Corrigendum: Exploration of Coriolis force on motion of air over the upper horizontal surface of a paraboloid of revolution. (2020 Phys. Scr. 95 035210). Phys. Scr. 2020, 95, 119501. [Google Scholar] [CrossRef]
  7. Choi, S.U.; Eastman, J.A. Enhancing Thermal Conductivity of Fluids with Nanoparticles; Technical Report; Argonne National Lab.: Lemont, IL, USA, 1995. [Google Scholar]
  8. Kaufui, V.; Omar, D. Applications of nanofluids: Current and future. Adv. Mech. Eng. 2010, 11, 105–132. [Google Scholar]
  9. Nguyen, C.; Desgranges, F.; Roy, G.; Galanis, N.; Maré, T.; Boucher, E.; Mintsa, H.A. Temperature and particle-size dependent viscosity data for water-based nanofluids–hysteresis phenomenon. Int. J. Heat Fluid Flow 2007, 28, 1492–1506. [Google Scholar] [CrossRef]
  10. Ali, M.; Shoukat, A.A.; Tariq, H.A.; Anwar, M.; Ali, H. Header design optimization of mini-channel heat sinks using CuO–H2O and Al2O3–H2 O nanofluids for thermal management. Arab. J. Sci. Eng. 2019, 44, 10327–10338. [Google Scholar] [CrossRef]
  11. Ullah, N.; Nadeem, S.; Khan, A.; Haq, R.U.; Tlili, I. Influence of metallic nanoparticles in water driven along a wavy circular cylinder. Chin. J. Phys. 2020, 63, 168–185. [Google Scholar] [CrossRef]
  12. Patil, P.; Kulkarni, M.; Hiremath, P. Effects of surface roughness on mixed convective nanofluid flow past an exponentially stretching permeable surface. Chin. J. Phys. 2020, 64, 203–218. [Google Scholar] [CrossRef]
  13. Animasaun, I.; Ibraheem, R.; Mahanthesh, B.; Babatunde, H. A meta-analysis on the effects of haphazard motion of tiny/nano-sized particles on the dynamics and other physical properties of some fluids. Chin. J. Phys. 2019, 60, 676–687. [Google Scholar] [CrossRef]
  14. Oke, A.S.; Mutuku, W.N.; Kimathi, M.; Animasaun, I.L. Insight into the dynamics of non-Newtonian Casson fluid over a rotating non-uniform surface subject to Coriolis force. Nonlinear Eng. 2020, 9, 398–411. [Google Scholar] [CrossRef]
  15. Oke, A.S.; Mutuku, W.N.; Kimathi, M.; Animasaun, I.L. Coriolis effects on MHD newtonian flow over a rotating non-uniform surface. Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci. 2021, 235, 3875–3887. [Google Scholar] [CrossRef]
  16. Oke, A.; Animasaun, I.; Mutuku, W.; Kimathi, M.; Shah, N.A.; Saleem, S. Significance of Coriolis force, volume fraction, and heat source/sink on the dynamics of water conveying 47 nm alumina nanoparticles over a uniform surface. Chin. J. Phys. 2021, 71, 716–727. [Google Scholar] [CrossRef]
  17. Wasp, E.J.; Kenny, J.P.; Gandhi, R.L. Solid–liquid flow: Slurry pipeline transportation [Pumps, valves, mechanical equipment, economics]. Ser. Bulk Mater. Handl. U. S. 1977, 1, 224. [Google Scholar]
  18. Pak, B.C.; Cho, Y.I. Hydrodynamic and heat transfer study of dispersed fluids with submicron metallic oxide particles. Exp. Heat Trans. 1998, 11, 151170. [Google Scholar] [CrossRef]
  19. Xuan, Y.; Roetzel, W. Conceptions for heat transfer correlation of nanofluids. Int. J. Heat Mass Transf. 2000, 43, 3701–3707. [Google Scholar] [CrossRef]
  20. Yang, R.Q. Infrared laser based on intersubband transitions in quantum wells. Superlattices Microstruct. 1995, 17, 77–83. [Google Scholar] [CrossRef]
  21. Shah, N.A.; Animasaun, I.; Wakif, A.; Koriko, O.; Sivaraj, R.; Adegbie, K.; Abdelmalek, Z.; Vaidyaa, H.; Ijirimoye, A.; Prasad, K. Significance of suction and dual stretching on the dynamics of various hybrid nanofluids: Comparative analysis between type I and type II models. Phys. Scr. 2020, 95, 095205. [Google Scholar] [CrossRef]
  22. Yahya, S.; Moghavvemi, M.; Mohamed, H.A. Artificial neural networks aided solution to the problem of geometrically bounded singularities and joint limits prevention of a three dimensional planar redundant manipulator. Neurocomputing 2014, 137, 34–46. [Google Scholar] [CrossRef]
  23. Mall, S.; Chakraverty, S. Numerical solution of nonlinear singular initial value problems of Emden–Fowler type using Chebyshev Neural Network method. Neurocomputing 2015, 149, 975–982. [Google Scholar] [CrossRef]
  24. Saritha, M.; Joseph, K.P.; Mathew, A.T. Classification of MRI brain images using combined wavelet entropy based spider web plots and probabilistic neural network. Pattern Recognit. Lett. 2013, 34, 2151–2156. [Google Scholar] [CrossRef]
  25. Hoptroff, R.G. The principles and practice of time series forecasting and business modelling using neural nets. Neural Comput. Appl. 1993, 1, 59–66. [Google Scholar] [CrossRef]
  26. Patan, K. Neural network-based model predictive control: Fault tolerance and stability. IEEE Trans. Control Syst. Technol. 2014, 23, 1147–1155. [Google Scholar] [CrossRef]
  27. Li, B.; Wang, Y.; Wang, Y.; Chen, Y.; Yang, H. Training itself: Mixed-signal training acceleration for memristor-based neural network. In Proceedings of the 2014 19th Asia and South Pacific Design Automation Conference (ASP-DAC), Singapore, 20–23 January 2014; pp. 361–366. [Google Scholar]
  28. Nusbaum, M.P.; Blitz, D.M.; Swensen, A.M.; Wood, D.; Marder, E. The roles of co-transmission in neural network modulation. Trends Neurosci. 2001, 24, 146–154. [Google Scholar] [CrossRef]
  29. Sridhara, V.; Satapathy, L.N. Al2O3-based nanofluids: A review. Nanoscale Res. Lett. 2011, 6, 456. [Google Scholar] [CrossRef] [Green Version]
  30. Cieśliński, J.T.; Ronewicz, K.; Smoleń, S. Measurement of temperature-dependent viscosity and thermal conductivity of alumina and titania thermal oil nanofluids. Arch. Thermodyn. 2015, 36, 35–47. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Graphical representation of the water having 47 nm nanoparticles of alumina of the transport phenomenon.
Figure 1. Graphical representation of the water having 47 nm nanoparticles of alumina of the transport phenomenon.
Nanomaterials 12 00878 g001
Figure 2. Basic design of a neural network.
Figure 2. Basic design of a neural network.
Nanomaterials 12 00878 g002
Figure 3. Showing the model which is generated through nftool command in matlab window.
Figure 3. Showing the model which is generated through nftool command in matlab window.
Nanomaterials 12 00878 g003
Figure 4. Flow chart of the proposed methodology.
Figure 4. Flow chart of the proposed methodology.
Nanomaterials 12 00878 g004
Figure 5. Graphical representation of F and the residual error graphs when rotation parameter K varies. (a) Graphical illustration of F when rotation parameter K varies. (b) Error in F when K = 0.1. (c) Error in F when K = 0.2. (d) Error in F when K = 0.3.
Figure 5. Graphical representation of F and the residual error graphs when rotation parameter K varies. (a) Graphical illustration of F when rotation parameter K varies. (b) Error in F when K = 0.1. (c) Error in F when K = 0.2. (d) Error in F when K = 0.3.
Nanomaterials 12 00878 g005aNanomaterials 12 00878 g005b
Figure 6. Graphical illustration of the effect of rotation parameter K on H and graphical illustration of error in the output data obtained through RK4 and LMBNN of H. (a) Graphs of H when K varies. (b) Error in H when K = 0.1. (c) Error in H when K = 0.2. (d) Error in H when K = 0.3.
Figure 6. Graphical illustration of the effect of rotation parameter K on H and graphical illustration of error in the output data obtained through RK4 and LMBNN of H. (a) Graphs of H when K varies. (b) Error in H when K = 0.1. (c) Error in H when K = 0.2. (d) Error in H when K = 0.3.
Nanomaterials 12 00878 g006
Figure 7. Figures of Θ and error graphs in the solution obtained from RK4 and LMBNN when K varies. (a) Graph of Θ when K varies. (b) Error in Θ K = 0.1. (c) Error in Θ K = 0.2. (d) Error in Θ K = 0.3.
Figure 7. Figures of Θ and error graphs in the solution obtained from RK4 and LMBNN when K varies. (a) Graph of Θ when K varies. (b) Error in Θ K = 0.1. (c) Error in Θ K = 0.2. (d) Error in Θ K = 0.3.
Nanomaterials 12 00878 g007
Figure 8. Effect of rotation parameter K on Φ and the graphs of error in numerical outputs of RK4 and LMBNN. (a) Numerical solutions of Φ when K varies. (b) Error in Φ when K = 0.1. (c) Error in Φ when K = 0.2. (d) Error in Φ when K = 0.3.
Figure 8. Effect of rotation parameter K on Φ and the graphs of error in numerical outputs of RK4 and LMBNN. (a) Numerical solutions of Φ when K varies. (b) Error in Φ when K = 0.1. (c) Error in Φ when K = 0.2. (d) Error in Φ when K = 0.3.
Nanomaterials 12 00878 g008
Figure 9. Performance, training, error histogram, regression, and fitness graphs of our system were obtained through LMBNN using the value of rotation parameter K = 0.1. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Figure 9. Performance, training, error histogram, regression, and fitness graphs of our system were obtained through LMBNN using the value of rotation parameter K = 0.1. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Nanomaterials 12 00878 g009
Figure 10. Figures of performance, training, error histogram, regression, and fitness of our system were obtained from LMBNN while using the value of rotation parameter K = 0.2. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Figure 10. Figures of performance, training, error histogram, regression, and fitness of our system were obtained from LMBNN while using the value of rotation parameter K = 0.2. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Nanomaterials 12 00878 g010
Figure 11. Performance, training, error histogram, regression, and fitness graphs of our system were obtained through LMBNN while using the value of rotation parameter K = 0.3. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Figure 11. Performance, training, error histogram, regression, and fitness graphs of our system were obtained through LMBNN while using the value of rotation parameter K = 0.3. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Nanomaterials 12 00878 g011
Figure 12. Graphical illustration of F and the residual error graphs in the outputs of RK4 and LMBNN when ϕ varies. (a) Numerical outputs of F when K varies. (b) Error in F when ϕ = 0.1. (c) Error in F when ϕ = 0.2. (d) Error in F when ϕ = 0.3.
Figure 12. Graphical illustration of F and the residual error graphs in the outputs of RK4 and LMBNN when ϕ varies. (a) Numerical outputs of F when K varies. (b) Error in F when ϕ = 0.1. (c) Error in F when ϕ = 0.2. (d) Error in F when ϕ = 0.3.
Nanomaterials 12 00878 g012
Figure 13. Effects of volume fraction ϕ on H and residual error graphs in the variance of ϕ . (a) Numerical outputs of H through RK4 and LMBNN for different values of ϕ . (b) Error in H when ϕ = 0.1. (c) Error in H when ϕ = 0.2. (d) Error in H when ϕ = 0.3.
Figure 13. Effects of volume fraction ϕ on H and residual error graphs in the variance of ϕ . (a) Numerical outputs of H through RK4 and LMBNN for different values of ϕ . (b) Error in H when ϕ = 0.1. (c) Error in H when ϕ = 0.2. (d) Error in H when ϕ = 0.3.
Nanomaterials 12 00878 g013
Figure 14. Graph of Θ when ϕ varies with the comparison of RK4 and LMBNN techniques and the graphs of the residual error between both the techniques. (a) Graphs of outputs of Θ obtained from RK4 and LMBNN technique when ϕ varies. (b) Error in Θ when ϕ = 0.1. (c) Error in Θ when ϕ = 0.2. (d) Error in Θ when ϕ = 0.3.
Figure 14. Graph of Θ when ϕ varies with the comparison of RK4 and LMBNN techniques and the graphs of the residual error between both the techniques. (a) Graphs of outputs of Θ obtained from RK4 and LMBNN technique when ϕ varies. (b) Error in Θ when ϕ = 0.1. (c) Error in Θ when ϕ = 0.2. (d) Error in Θ when ϕ = 0.3.
Nanomaterials 12 00878 g014
Figure 15. Graphs of Φ representing the outputs obtained from the RK4 and LMB-NN methods along with the graphs of the residual error between both the outputs. (a) Graphical representation of the outputs of Φ through the RK4 and LMBNN when ϕ varies. (b) Error graph of Φ when ϕ = 0.1. (c) Error graph of Φ when ϕ = 0.2. (d) Error graph of Φ when ϕ = 0.3.
Figure 15. Graphs of Φ representing the outputs obtained from the RK4 and LMB-NN methods along with the graphs of the residual error between both the outputs. (a) Graphical representation of the outputs of Φ through the RK4 and LMBNN when ϕ varies. (b) Error graph of Φ when ϕ = 0.1. (c) Error graph of Φ when ϕ = 0.2. (d) Error graph of Φ when ϕ = 0.3.
Nanomaterials 12 00878 g015
Figure 16. Figures showing the performance, training, error histogram, regression, and fitness graphs of the LMB-NN when value of the volume fraction ϕ = 0.1. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Figure 16. Figures showing the performance, training, error histogram, regression, and fitness graphs of the LMB-NN when value of the volume fraction ϕ = 0.1. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Nanomaterials 12 00878 g016
Figure 17. Figures showing the performance, training, error histogram, regression, and fitness graphs of the LMB-NN when value of the volume fraction ϕ = 0.2. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Figure 17. Figures showing the performance, training, error histogram, regression, and fitness graphs of the LMB-NN when value of the volume fraction ϕ = 0.2. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Nanomaterials 12 00878 g017
Figure 18. The figures clearly show the mean square error (mse), error histogram, gradient, Mu, validation checks, training, regression, and the epoch takes to reach the best results when ϕ = 0.3. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Figure 18. The figures clearly show the mean square error (mse), error histogram, gradient, Mu, validation checks, training, regression, and the epoch takes to reach the best results when ϕ = 0.3. (a) Performance. (b) Training. (c) Error Histogram. (d) Regression. (e) Fitness.
Nanomaterials 12 00878 g018
Table 1. All the cases and scenarios discussed in the paper.
Table 1. All the cases and scenarios discussed in the paper.
Senarios CasesParameters
S1C 1K = 0.1
C 2K = 0.2
C 3K = 0.3
S2C 1 ϕ = 0.1
C 2 ϕ = 0.2
C 3 ϕ = 0.3
Table 2. Numerical results of ANN of different cases of the rotation parameter K and the volume fraction ϕ showing the time taken, training, validation, testing, mean square error, gradient, and epoch to reach the best solution.
Table 2. Numerical results of ANN of different cases of the rotation parameter K and the volume fraction ϕ showing the time taken, training, validation, testing, mean square error, gradient, and epoch to reach the best solution.
ScenarioCasesTimeMSEPerformanceGradiantMuEpoch
TrainingValidationTesting
KK = 0.1179.57 × 10 12 1.27 × 10 11 1.43 × 10 11 9.57 × 10 12 9.76 × 10 8 1.00 × 10 12 139
K = 0.2111.00 × 10 11 1.29 × 10 11 1.22 × 10 11 1.00 × 10 11 9.95 × 10 8 1.00 × 10 12 132
K = 0.3149.92 × 10 12 1.02 × 10 11 1.51 × 10 10 9.92 × 10 12 9.97 × 10 8 1.00 × 10 12 155
ϕ ϕ = 0.1191.69 × 10 11 2.12 × 10 11 2.34 × 10 11 1.69 × 10 11 9.99 × 10 8 1.00 × 10 11 721
ϕ = 0.2102.50 × 10 11 3.61 × 10 11 3.57 × 10 11 2.50 × 10 11 9.99 × 10 8 1.00 × 10 11 365
ϕ = 0.3153.56 × 10 11 4.45 × 10 11 4.73 × 10 11 3.56 × 10 11 9.98 × 10 8 1.00 × 10 11 294
Table 3. Different values of F were obtained by RK4 and LMBNN methods when rotation parameter K varies (from 0.1–0.3).
Table 3. Different values of F were obtained by RK4 and LMBNN methods when rotation parameter K varies (from 0.1–0.3).
F[t]
tRK4LMBNNRK4LMBNNRK4LMBNN
K = 0.1K = 0.1K = 0.2K = 0.2K = 0.3K = 0.3
t = 0.0000.00064300.0010400.001605
t = 0.20.1837940.1839330.1834570.1835480.1830980.18311
t = 0.40.3403170.3404330.3389880.3390630.3375680.337585
t = 0.60.475990.4759810.4730620.4730390.4699360.46993
t = 0.80.5956330.5956290.5905910.5905930.58520.5852
t = 1.000.7026530.7026560.6950960.6950820.6870080.687014
t = 1.20.7992720.7992680.7889410.7889110.7778640.777864
t = 1.40.8867990.8868070.8735740.8735740.8593720.859354
t = 1.60.9658990.9659020.9497980.949790.9324780.932481
t = 1.81.036851.0368481.018011.0180050.9977070.997713
t = 2.01.099741.099741.078391.0783951.055341.055339
t = 2.21.154621.154621.131061.1310621.105581.105586
t = 2.41.201621.2016111.176181.176191.148621.148618
t = 2.61.240951.2409451.213971.2139691.18471.184696
t = 2.81.272981.2729831.24481.2447911.214181.214175
t = 3.01.29821.2981991.269121.2691011.237481.237475
t = 3.21.317211.3172081.287491.2874651.255131.255136
t = 3.41.330671.3306761.300541.3005271.26771.267698
t = 3.61.339331.3393341.308961.3089521.275841.275842
t = 3.81.343961.3439641.313481.3134671.280221.280214
t = 4.01.345331.345331.314821.3148161.281521.281517
Table 4. The values of H by RK4 and LMB-NN methods at 21 different points when rotation parameter K varies.
Table 4. The values of H by RK4 and LMB-NN methods at 21 different points when rotation parameter K varies.
H[t]
tRK4LMBNNRK4LMBNNRK4LMBNN
K = 0.1K = 0.1K = 0.2K = 0.2K = 0.3K = 0.3
t = 0.0000.00022600.00036500.000492
t = 0.20.0590890.0591380.0599550.0599550.0609120.060922
t = 0.40.1158410.1158540.1175840.1176030.1195110.119524
t = 0.60.1675630.1675620.1701890.1701770.1730970.173067
t = 0.80.2115340.2115560.2150340.2150330.2189110.218908
t = 1.000.2455020.2454650.2498350.249840.254640.254616
t = 1.20.2680670.2680370.2731550.2731390.2788040.278815
t = 1.40.2788410.2788680.2845650.2846010.2909270.291016
t = 1.60.2784010.2783430.28460.2846390.2915010.291523
t = 1.80.2680830.2680850.2745650.2745720.2817930.281789
t = 2.00.2497180.2497130.2562710.2562710.2635930.263593
t = 2.20.2253630.225370.2317690.2317720.2389430.238933
t = 2.40.1970770.1970740.2031310.2031230.2099270.209928
t = 2.60.1667650.1667660.1722860.172290.1784990.178504
t = 2.80.1360780.1360780.1409180.1409350.146380.146388
t = 3.00.1063660.1063670.110420.1104530.1150090.11502
t = 3.20.0786790.0786830.0818860.0819270.0855260.085514
t = 3.40.0537820.0537760.0561220.0561530.0587870.058794
t = 3.60.0321930.0321940.0336870.0337280.0353940.035389
t = 3.80.0142190.0142120.0149230.0149670.0157310.01574
t = 4.0−6.71 × 10 16 0.000214−7.27 × 10 16 0.000337−9.45 × 10 16 0.000264
Table 5. Comparison of RK4 and LM-BNN in Θ using different values of rotation parameter K.
Table 5. Comparison of RK4 and LM-BNN in Θ using different values of rotation parameter K.
Θ [ t ]
tRK4LMBNNRK4LMBNNRK4LMBNN
K = 0.1K = 0.1K = 0.2K = 0.2K = 0.3K = 0.3
t = 0.0000.00191500.00348100.005383
t = 0.20.5510840.5510720.5543310.554340.557940.558028
t = 0.40.8778230.8779170.8841180.8840220.8911150.891085
t = 0.61.01631.0157751.025261.0248871.035231.035112
t = 0.81.0191.0190421.030091.0301031.042431.042399
t = 1.000.9381670.9381050.9506970.950670.9646480.964616
t = 1.20.8160220.8159740.8292120.8292390.8439310.843901
t = 1.40.6818570.681810.6949480.6949270.70960.70959
t = 1.60.5533370.5533120.5656820.5656870.5795540.579547
t = 1.80.4393930.4393980.4505250.4505080.4630920.463091
t = 2.00.3431810.3431760.3528310.3528160.3637810.363787
t = 2.20.2644870.264480.2725650.2725730.2817780.281775
t = 2.40.2014180.2014160.2079690.2079530.2154790.21548
t = 2.60.1514730.1514720.156630.1566360.1625710.162578
t = 2.80.1121460.1121450.1160860.1161040.1206450.12065
t = 3.00.0812120.0812110.0841230.0841540.0875040.087512
t = 3.20.0568290.0568280.0588890.0589230.0612890.061276
t = 3.40.0375290.0375260.0388960.0389190.0404930.040497
t = 3.60.0221720.0221720.0229810.0230080.0239270.023923
t = 3.80.0098850.009880.0102440.0102740.0106660.010673
t = 4.0−1.58 × 10 16 0.0001542.22 × 10 16 0.00024−1.39 × 10 17 0.000185
Table 6. Different outputs of Φ using different values of rotation parameter K by solving on RK4 and LMBNN.
Table 6. Different outputs of Φ using different values of rotation parameter K by solving on RK4 and LMBNN.
Φ [t]
tRK4LMBNNRK4LMBNNRK4LMBNN
K = 0.1K = 0.1K = 0.2K = 0.2K = 0.3K = 0.3
t = 0.0000.000011200.00015300.000766
t = 0.2−0.50429−0.50421−0.50798−0.50803−0.51209−0.5121
t = 0.4−0.77634−0.77635−0.78346−0.78354−0.7914−0.79146
t = 0.6−0.85358−0.85254−0.86364−0.86241−0.87485−0.87385
t = 0.8−0.79515−0.79523−0.80745−0.80747−0.82119−0.82111
t = 1.00−0.6608−0.66077−0.67447−0.67443−0.68974−0.68981
t = 1.2−0.49841−0.49835−0.51244−0.51247−0.52817−0.52822
t = 1.4−0.34007−0.34004−0.35351−0.35349−0.36865−0.36868
t = 1.6−0.20365−0.20363−0.21571−0.21572−0.2294−0.22939
t = 1.8−0.09637−0.09637−0.10655−0.10655−0.1182−0.11819
t = 2.0−0.01865−0.01865−0.02672−0.02673−0.03606−0.03607
t = 2.20.0328360.0328350.0268520.0268350.0198260.019831
t = 2.40.062940.0629430.0588350.058840.0539180.053924
t = 2.60.0767260.0767310.0741840.0741570.0710430.071049
t = 2.80.0787410.0787620.0774020.0773380.0756490.075584
t = 3.00.0727280.0727240.0722360.0722620.0714850.071493
t = 3.20.061580.0615830.0616190.0616360.0615130.061517
t = 3.40.0474390.0474360.0477410.0477520.0479780.047979
t = 3.60.0318270.0318280.0321780.0322070.032510.032504
t = 3.80.0157860.0157840.0160210.0160630.0162580.016264
t = 4.01.54 × 10 16 −2.74 × 10 5 3.57 × 10 16 −1.2 × 10 4 6.26 × 10 16 −4.94 × 10 5
Table 7. Residual error in F when rotation parameter K varies.
Table 7. Residual error in F when rotation parameter K varies.
Residual Error in F
tResidual Error When K = 0.1Residual Error When K = 0.2Residual Error When K = 0.3
0−6.43 × 10 4 −1.04 × 10 3 −1.60 × 10 3
0.2−1.39 × 10 4 −9.10 × 10 5 −1.20 × 10 5
0.4−1.16 × 10 4 −7.53 × 10 5 −1.70 × 10 5
0.69.43 × 10 6 2.34 × 10 5 6.19 × 10 6
0.84.25 × 10 6 −1.68 × 10 6 4.68 × 10 7
1−2.63 × 10 6 1.36 × 10 5 −5.88 × 10 6
1.24.02 × 10 6 3.05 × 10 5 −2.81 × 10 7
1.4−7.89 × 10 6 −2.06 × 10 7 1.84 × 10 5
1.6−3.00 × 10 6 7.54 × 10 6 −3.28 × 10 6
1.81.94 × 10 6 4.84 × 10 6 −5.69 × 10 6
2−2.32 × 10 7 −5.38 × 10 6 1.34 × 10 6
2.2−2.44 × 10 8 −1.55 × 10 6 −6.47 × 10 6
2.48.79 × 10 6 −1.00 × 10 5 2.30 × 10 6
2.65.18 × 10 6 1.21 × 10 6 3.56 × 10 6
2.8−2.60 × 10 6 9.06 × 10 6 5.21 × 10 6
31.43 × 10 6 1.87 × 10 5 4.93 × 10 6
3.22.30 × 10 6 2.50 × 10 5 −5.97 × 10 6
3.4−5.54 × 10 6 1.29 × 10 5 2.36 × 10 6
3.6−4.40 × 10 6 7.97 × 10 6 −2.48 × 10 6
3.8−4.43 × 10 6 1.32 × 10 5 5.80 × 10 6
42.60 × 10 7 3.90 × 10 6 2.73 × 10 6
Table 8. Residual error in H between the Mathematica and neural outputs.
Table 8. Residual error in H between the Mathematica and neural outputs.
Residual Error in H
tResidual Error When K = 0.1Residual Error When K = 0.2Residual Error When K = 0.3
0−2.26 × 10 4 −3.65 × 10 4 −4.92 × 10 4
0.2−4.90 × 10 5 1.90 × 10 7 −9.72 × 10 6
0.4−1.28 × 10 5 −1.91 × 10 5 −1.25 × 10 5
0.61.46 × 10 6 1.15 × 10 5 2.98 × 10 5
0.8−2.17 × 10 5 1.22 × 10 6 3.04 × 10 6
13.68 × 10 5 −5.03 × 10 6 2.39 × 10 5
1.22.96 × 10 5 1.63 × 10 5 −1.13 × 10 5
1.4−2.69 × 10 5 −3.57 × 10 5 −8.88 × 10 5
1.65.77 × 10 5 −3.91 × 10 5 −2.24 × 10 5
1.8−2.29 × 10 6 −6.56 × 10 6 3.60 × 10 6
25.42 × 10 6 −1.87 × 10 7 2.34 × 10 7
2.2−6.81 × 10 6 −2.91 × 10 6 9.74 × 10 6
2.43.14 × 10 6 7.62 × 10 6 −9.27 × 10 7
2.6−7.15 × 10 7 −3.81 × 10 6 −4.54 × 10 6
2.8−2.39 × 10 7 −1.68 × 10 5 −7.62 × 10 6
3−7.62 × 10 7 −3.33 × 10 5 −1.07 × 10 5
3.2−3.75 × 10 6 −4.10 × 10 5 1.21 × 10 5
3.46.78 × 10 6 −3.06 × 10 5 −6.99 × 10 6
3.6−9.27 × 10 7 −4.06 × 10 5 5.52 × 10 6
3.86.95 × 10 6 −4.39 × 10 5 −9.77 × 10 6
4−2.14 × 10 4 −3.37 × 10 4 −2.64 × 10 4
Table 9. Residual error with the reference solution of RK4 in Θ . These residual errors are obtained from Matlab.
Table 9. Residual error with the reference solution of RK4 in Θ . These residual errors are obtained from Matlab.
Residual Error in Θ
tResidual Error When K = 0.1Residual Error When K = 0.2Residual Error When K = 0.3
0−1.92 × 10 3 −3.48 × 10 3 −5.38 × 10 3
0.21.25 × 10 5 −9.31 × 10 6 −8.81 × 10 5
0.4−9.39 × 10 5 9.61 × 10 5 3.02 × 10 5
0.65.25 × 10 4 3.73 × 10 4 1.18 × 10 4
0.8−4.17 × 10 5 −1.34 × 10 5 3.07 × 10 5
16.16 × 10 5 2.73 × 10 5 3.24 × 10 5
1.24.81 × 10 5 −2.70 × 10 5 3.02 × 10 5
1.44.71 × 10 5 2.07 × 10 5 1.01 × 10 5
1.62.50 × 10 5 −5.03 × 10 6 7.11 × 10 6
1.8−5.45 × 10 6 1.70 × 10 5 9.49 × 10 7
25.23 × 10 6 1.48 × 10 5 −6.01 × 10 6
2.27.39 × 10 6 −7.77 × 10 6 2.62 × 10 6
2.42.17 × 10 6 1.60 × 10 5 −6.64 × 10 7
2.61.14 × 10 6 −6.22 × 10 6 −7.10 × 10 6
2.85.46 × 10 7 −1.77 × 10 5 −4.74 × 10 6
35.71 × 10 7 −3.08 × 10 5 −8.59 × 10 6
3.25.45 × 10 7 −3.44 × 10 5 1.25 × 10 5
3.42.62 × 10 6 −2.29 × 10 5 −3.71 × 10 6
3.6−5.12 × 10 7 −2.79 × 10 5 4.18 × 10 6
3.84.56 × 10 6 −2.98 × 10 5 −6.62 × 10 6
4−1.54 × 10 4 −2.40 × 10 4 −1.85 × 10 4
Table 10. Residual error in RK4 and LMB-NN in Φ .
Table 10. Residual error in RK4 and LMB-NN in Φ .
Residual Error in Φ
tResidual Error When K = 0.1Residual Error When K = 0.2Residual Error When K = 0.3
0−1.12 × 10 5 −1.53 × 10 4 −7.66 × 10 4
0.2−7.59 × 10 5 4.73 × 10 5 6.38 × 10 6
0.48.83 × 10 6 7.49 × 10 5 6.41 × 10 5
0.6−1.04 × 10 3 −1.23 × 10 3 −1.00 × 10 3
0.88.72 × 10 5 1.97 × 10 5 −7.58 × 10 5
1−3.62 × 10 5 −3.18 × 10 5 6.76 × 10 5
1.2−5.75 × 10 5 3.16 × 10 5 5.20 × 10 5
1.4−3.38 × 10 5 −1.60 × 10 5 2.92 × 10 5
1.6−1.34 × 10 5 6.64 × 10 6 −7.06 × 10 6
1.83.54 × 10 6 −5.05 × 10 6 −1.06 × 10 5
2−1.39 × 10 7 8.68 × 10 6 8.21 × 10 6
2.28.88 × 10 7 1.73 × 10 5 −5.32 × 10 6
2.4−3.96 × 10 6 −5.11 × 10 6 −6.21 × 10 6
2.6−4.35 × 10 6 2.76 × 10 5 −5.77 × 10 6
2.8−2.01 × 10 5 6.31 × 10 5 6.54 × 10 5
33.46 × 10 6 −2.58 × 10 5 −8.42 × 10 6
3.2−2.51 × 10 6 −1.68 × 10 5 −4.01 × 10 6
3.43.31 × 10 6 −1.12 × 10 5 −1.36 × 10 6
3.6−6.46 × 10 7 −2.97 × 10 5 6.47 × 10 6
3.81.78 × 10 6 −4.25 × 10 5 −5.95 × 10 6
42.74 × 10 5 1.19 × 10 4 4.94 × 10 5
Table 11. Numerical values of F obtained through RK4 and LMBNN using different values of volume fraction ϕ .
Table 11. Numerical values of F obtained through RK4 and LMBNN using different values of volume fraction ϕ .
F[t]
tRK4LMBNNRK4LMBNNRK4LMBNN
ϕ = 0.1 ϕ = 0.1 ϕ = 0.2 ϕ = 0.2 ϕ = 0.3 ϕ = 0.3
t = 0.0006.89 × 10 8 05.75 × 10 6 0−1.04 × 10 6
t = 0.20.1870680.1870680.1868210.1868270.1864960.18649
t = 0.40.3532670.3532680.3523190.3523090.351060.35106
t = 0.60.5044780.5044770.5025140.5025110.4998430.499852
t = 0.80.6445060.6444980.6414360.6414330.6371040.637107
t = 1.000.775330.7753320.7713280.7713360.7653710.765375
t = 1.20.897540.8975410.8930.8930040.8857340.885733
t = 1.41.010851.0108491.006291.0062880.9982390.998239
t = 1.61.114561.1145571.110531.1105281.102311.102315
t = 1.81.207941.2079341.20491.2048991.197121.197119
t = 2.01.290441.2904351.288731.2887321.281911.281915
t = 2.21.361821.3618241.361631.3616271.356121.356118
t = 2.41.42221.4222041.423561.4235691.419561.419565
t = 2.61.471981.4719861.474821.4748281.472361.472357
t = 2.81.511851.5118511.5161.5160061.5151.514998
t = 3.01.542681.542681.547931.5479251.548181.548188
t = 3.21.565471.5654751.571581.5715761.572861.572856
t = 3.41.58131.5813041.588031.5880271.590061.590061
t = 3.61.591281.591281.59841.5984031.600941.600943
t = 3.81.596491.5964891.603821.6038251.606631.606631
t = 4.01.5981.5980011.605391.6053961.608281.608284
Table 12. Comparison of output data obtained through RK4 and LMBNN methods when volume fraction ϕ varies.
Table 12. Comparison of output data obtained through RK4 and LMBNN methods when volume fraction ϕ varies.
H[t]
tRK4LMBNNRK4LMBNNRK4LMBNN
ϕ = 0.1 ϕ = 0.1 ϕ = 0.2 ϕ = 0.2 ϕ = 0.3 ϕ = 0.3
t = 0.0001.38 × 10 6 03.57 × 10 6 0−1.06 × 10 6
t = 0.20.1120120.1120130.1091360.109140.1051570.105154
t = 0.40.2180490.2180520.2128540.2128480.2054690.205469
t = 0.60.3116040.3116040.3052110.305210.2955920.295595
t = 0.80.3867180.3867130.3806040.3806030.3703670.370366
t = 1.000.4393330.4393340.4349190.4349220.4257810.42578
t = 1.20.4680080.4680080.4662860.4662860.459670.459669
t = 1.40.4738380.4738380.4751870.4751890.4719540.471955
t = 1.60.4598370.4598360.4640520.4640550.4644110.464408
t = 1.80.4301220.4301240.4365870.4365880.440180.44018
t = 2.00.3891720.3891740.3970670.3970660.403150.403147
t = 2.20.3412690.3412690.3497570.349760.3574050.357411
t = 2.40.290160.2901560.2985090.2985030.3067990.306794
t = 2.60.2388960.2388960.2465450.2465380.2546740.254683
t = 2.80.1898070.1898120.1963840.1963820.203740.203741
t = 3.00.1445560.1445620.1498650.1498730.1560450.156043
t = 3.20.1042360.1042330.1082260.108230.1130220.113023
t = 3.40.0694750.0694780.0722080.0722110.0755810.075591
t = 3.60.0405430.0405380.0421620.0421510.0442030.044191
t = 3.80.0174460.0174530.0181430.0181380.019040.019047
t = 4.0−2.14 × 10 15 −5.78 × 10 6 9.35 × 10 16 1.09 × 10 5 2.39 × 10 15 1.10 × 10 5
Table 13. Comparison of different values of Θ between RK4 and LMB-NN when ϕ varies.
Table 13. Comparison of different values of Θ between RK4 and LMB-NN when ϕ varies.
Θ [t]
tRK4LMBNNRK4LMBNNRK4LMBNN
ϕ = 0.1 ϕ = 0.1 ϕ = 0.2 ϕ = 0.2 ϕ = 0.3 ϕ = 0.3
t = 0.0009.21 × 10 6 01.17 × 10 5 01.52 × 10 5
t = 0.21.555341.555341.294231.2942321.074131.074122
t = 0.42.244182.2441791.955941.9559311.683271.683266
t = 0.62.309332.3093272.125252.1252431.912411.912411
t = 0.82.046562.0465631.989271.9892741.878471.878473
t = 1.001.674461.6744581.709411.7094131.692981.692973
t = 1.21.31121.3111961.393021.3930241.441561.441564
t = 1.41.003581.0035781.097751.0977511.179941.17994
t = 1.60.7604520.7604470.8473380.8473390.9390350.939026
t = 1.80.5746020.5746050.6461970.6462020.7324190.732429
t = 2.00.4343260.4343330.489540.4895430.5630260.563016
t = 2.20.3284660.328460.3693860.3693920.4281080.428115
t = 2.40.2479910.2479860.2776260.2776190.322490.322486
t = 2.60.1861960.1861960.2073570.2073530.2405130.240525
t = 2.80.1382290.1382310.1531630.153160.1770070.177006
t = 3.00.1005960.1006010.1109770.1109830.1276680.127665
t = 3.20.0708040.0708040.0778280.0778320.0891140.089114
t = 3.40.0470350.0470360.0515490.0515510.0587680.058775
t = 3.60.0279380.0279340.0305510.0305430.0346970.034688
t = 3.80.0125140.0125180.013660.0136560.0154650.015471
t = 4.06.07 × 10 17 −5.54 × 10 6 4.81 × 10 16 1.18 × 10 5 5.64 × 10 16 1.21 × 10 5
Table 14. Numerical outputs of Φ at different values of ϕ using RK4 and LMBNN techniques.
Table 14. Numerical outputs of Φ at different values of ϕ using RK4 and LMBNN techniques.
Φ [t]
tRK4LMBNNRK4LMBNNRK4LMBNN
ϕ = 0.1 ϕ = 0.1 ϕ = 0.2 ϕ = 0.2 ϕ = 0.3 ϕ = 0.3
t = 0.000−8.62 × 10 6 0−2.56 × 10 6 0−1.50 × 10 5
t = 0.2−1.43028−1.43028−1.18136−1.18136−0.97385−0.97384
t = 0.4−1.97489−1.97489−1.71308−1.71309−1.46774−1.46773
t = 0.6−1.88905−1.88906−1.74254−1.74253−1.57038−1.57038
t = 0.8−1.49282−1.49283−1.47602−1.47602−1.41311−1.41312
t = 1.00−1.02319−1.02319−1.09275−1.09275−1.12306−1.12306
t = 1.2−0.60554−0.60553−0.71023−0.71023−0.79701−0.79701
t = 1.4−0.28453−0.28453−0.38801−0.388−0.49553−0.49553
t = 1.6−0.06174−0.06174−0.14572−0.14572−0.24874−0.24873
t = 1.80.0790540.0790550.0195180.019516−0.06549−0.0655
t = 2.00.1581810.1581770.1206930.1206880.0579840.057989
t = 2.20.1939920.1939990.1733630.1733630.1315670.131568
t = 2.40.2010220.2010210.1918670.1918670.1669140.166912
t = 2.60.1898870.1898850.1877850.1877790.1749840.174985
t = 2.80.1680360.1680420.1697720.1697730.1649060.164909
t = 3.00.1405440.1405470.1439550.143960.1437030.143704
t = 3.20.1107340.110730.1144620.1144640.1164360.116437
t = 3.40.0807310.0807360.0839750.0839780.0865760.086585
t = 3.60.0518440.0518380.0541650.0541530.0563880.056375
t = 3.80.024810.0248190.0260020.0259940.0272650.027273
t = 4.06.51 × 10 16 −1.35 × 10 5 −1.59 × 10 15 2.75 × 10 5 −1.23 × 10 15 2.64 × 10 5
Table 15. The residual error in F when ϕ varies are shown in the table.
Table 15. The residual error in F when ϕ varies are shown in the table.
Residual Error in F
tResidual Error When ϕ = 0.1Residual Error When ϕ = 0.2Residual Error When ϕ = 0.3
0−6.89 × 10 8 −5.75 × 10 6 1.04 × 10 6
0.22.26 × 10 7 −5.66 × 10 6 5.87 × 10 6
0.4−1.44 × 10 6 9.77 × 10 6 −3.97 × 10 7
0.69.52 × 10 7 2.77 × 10 6 −9.14 × 10 6
0.87.88 × 10 6 2.80 × 10 6 −2.94 × 10 6
1−2.36 × 10 6 −7.82 × 10 6 −4.30 × 10 6
1.2−8.33 × 10 7 −3.78 × 10 6 5.70 × 10 7
1.48.21 × 10 7 2.15 × 10 6 −3.74 × 10 7
1.63.36 × 10 6 1.96 × 10 6 −4.57 × 10 6
1.85.80 × 10 6 1.32 × 10 6 5.16 × 10 7
24.61 × 10 6 −1.86 × 10 6 −4.65 × 10 6
2.2−4.33 × 10 6 2.70 × 10 6 2.36 × 10 6
2.4−3.78 × 10 6 −8.82 × 10 6 −5.48 × 10 6
2.6−6.32 × 10 6 −7.79 × 10 6 3.19 × 10 6
2.8−5.86 × 10 7 −6.04 × 10 6 2.47 × 10 6
3−7.84 × 10 9 5.48 × 10 6 −7.87 × 10 6
3.2−4.60 × 10 6 3.94 × 10 6 3.69 × 10 6
3.4−3.76 × 10 6 3.12 × 10 6 −1.11 × 10 6
3.63.55 × 10 7 −2.60 × 10 6 −3.37 × 10 6
3.81.17 × 10 6 −4.80 × 10 6 −1.41 × 10 6
4−8.76 × 10 7 −6.28 × 10 6 −4.47 × 10 6
Table 16. The residual error between RK4 and LMBNN in the outputs of H.
Table 16. The residual error between RK4 and LMBNN in the outputs of H.
Residual Error in H
tResidual Error When ϕ = 0.1Residual Error When ϕ = 0.2Residual Error When ϕ = 0.3
0−1.38 × 10 6 −3.57 × 10 6 1.06 × 10 6
0.2−5.09 × 10 7 −4.04 × 10 6 3.49 × 10 6
0.4−2.59 × 10 6 6.17 × 10 6 4.49 × 10 7
0.62.84 × 10 7 7.50 × 10 7 −2.63 × 10 6
0.85.39 × 10 6 1.44 × 10 6 5.63 × 10 7
1−8.18 × 10 7 −2.79 × 10 6 6.30 × 10 7
1.21.46 × 10 8 −3.01 × 10 7 7.21 × 10 7
1.4−2.48 × 10 8 −1.67 × 10 6 −8.33 × 10 7
1.65.24 × 10 7 −3.24 × 10 6 2.79 × 10 6
1.8−2.16 × 10 6 −1.28 × 10 6 1.95 × 10 7
2−1.68 × 10 6 1.22 × 10 6 3.40 × 10 6
2.2−1.12 × 10 7 −3.34 × 10 6 −6.20 × 10 6
2.43.98 × 10 6 6.37 × 10 6 5.24 × 10 6
2.62.05 × 10 7 7.08 × 10 6 −9.05 × 10 6
2.8−5.21 × 10 6 1.68 × 10 6 −1.07 × 10 6
3−5.63 × 10 6 −8.12 × 10 6 2.08 × 10 6
3.22.52 × 10 6 −4.18 × 10 6 −9.71 × 10 7
3.4−3.59 × 10 6 −2.92 × 10 6 −1.03 × 10 5
3.65.45 × 10 6 1.09 × 10 5 1.18 × 10 5
3.8−6.66 × 10 6 5.54 × 10 6 −7.04 × 10 6
45.78 × 10 6 −1.09 × 10 5 −1.10 × 10 5
Table 17. The residual error in Θ when volume fraction ϕ varies.
Table 17. The residual error in Θ when volume fraction ϕ varies.
Residual Error in Θ
tResidual Error When ϕ = 0.1Residual Error When ϕ = 0.2Residual Error When ϕ = 0.3
0−9.21 × 10 6 −1.17 × 10 5 −1.52 × 10 5
0.2−2.88 × 10 7 −2.39 × 10 6 8.44 × 10 6
0.48.56 × 10 7 8.62 × 10 6 3.85 × 10 6
0.62.76 × 10 6 6.57 × 10 6 −1.01 × 10 6
0.8−2.58 × 10 6 −3.67 × 10 6 −2.56 × 10 6
12.08 × 10 6 −2.62 × 10 6 6.56 × 10 6
1.24.22 × 10 6 −3.68 × 10 6 −3.66 × 10 6
1.41.99 × 10 6 −1.32 × 10 6 −3.82 × 10 7
1.64.59 × 10 6 −1.18 × 10 6 8.73 × 10 6
1.8−2.58 × 10 6 −5.09 × 10 6 −9.82 × 10 6
2−6.68 × 10 6 −2.64 × 10 6 9.86 × 10 6
2.26.25 × 10 6 −5.51 × 10 6 −7.37 × 10 6
2.45.30 × 10 6 7.41 × 10 6 4.37 × 10 6
2.6−4.81 × 10 7 4.50 × 10 6 −1.19 × 10 5
2.8−2.07 × 10 6 3.35 × 10 6 5.52 × 10 7
3−5.27 × 10 6 −6.20 × 10 6 3.29 × 10 6
3.25.32 × 10 8 −3.95 × 10 6 −5.06 × 10 7
3.4−1.81 × 10 6 −1.84 × 10 6 −7.95 × 10 6
3.63.75 × 10 6 7.71 × 10 6 9.03 × 10 6
3.8−4.51 × 10 6 3.90 × 10 6 −5.38 × 10 6
45.54 × 10 6 −1.18 × 10 5 −1.21 × 10 5
Table 18. The residual error in Φ are shown in the following the table.
Table 18. The residual error in Φ are shown in the following the table.
Residual Error in Φ
tResidual Error When ϕ = 0.1Residual Error When ϕ = 0.2Residual Error When ϕ = 0.3
08.62 × 10 6 2.56 × 10 6 1.50 × 10 5
0.2−3.20 × 10 6 1.63 × 10 6 −3.56 × 10 6
0.41.09 × 10 7 5.54 × 10 6 −1.17 × 10 5
0.65.06 × 10 6 −6.33 × 10 6 2.49 × 10 6
0.85.29 × 10 6 2.26 × 10 6 1.00 × 10 5
13.24 × 10 6 −3.88 × 10 6 −3.43 × 10 6
1.2−3.15 × 10 6 −1.04 × 10 6 2.45 × 10 6
1.41.54 × 10 6 −5.99 × 10 6 −5.66 × 10 6
1.6−4.82 × 10 6 −3.24 × 10 6 −4.33 × 10 6
1.8−7.96 × 10 7 2.23 × 10 6 1.03 × 10 5
24.41 × 10 6 4.83 × 10 6 −5.61 × 10 6
2.2−6.79 × 10 6 2.33 × 10 7 −8.62 × 10 7
2.41.39 × 10 6 4.53 × 10 7 2.45 × 10 6
2.61.56 × 10 6 5.67 × 10 6 −5.85 × 10 7
2.8−5.62 × 10 6 −7.14 × 10 7 −2.98 × 10 6
3−2.52 × 10 6 −4.74 × 10 6 −1.38 × 10 6
3.24.40 × 10 6 −1.72 × 10 6 −1.43 × 10 6
3.4−4.35 × 10 6 −3.41 × 10 6 −9.23 × 10 6
3.66.31 × 10 6 1.20 × 10 5 1.28 × 10 5
3.8−8.53 × 10 6 7.89 × 10 6 −8.32 × 10 6
41.35 × 10 5 −2.75 × 10 5 −2.64 × 10 5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ganie, A.H.; Fazal, F.; Tavera Romero, C.A.; Sulaiman, M. Quantitative Features Analysis of Water Carrying Nanoparticles of Alumina over a Uniform Surface. Nanomaterials 2022, 12, 878. https://doi.org/10.3390/nano12050878

AMA Style

Ganie AH, Fazal F, Tavera Romero CA, Sulaiman M. Quantitative Features Analysis of Water Carrying Nanoparticles of Alumina over a Uniform Surface. Nanomaterials. 2022; 12(5):878. https://doi.org/10.3390/nano12050878

Chicago/Turabian Style

Ganie, Abdul Hamid, Fazlullah Fazal, Carlos Andrés Tavera Romero, and Muhammad Sulaiman. 2022. "Quantitative Features Analysis of Water Carrying Nanoparticles of Alumina over a Uniform Surface" Nanomaterials 12, no. 5: 878. https://doi.org/10.3390/nano12050878

APA Style

Ganie, A. H., Fazal, F., Tavera Romero, C. A., & Sulaiman, M. (2022). Quantitative Features Analysis of Water Carrying Nanoparticles of Alumina over a Uniform Surface. Nanomaterials, 12(5), 878. https://doi.org/10.3390/nano12050878

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop