Blog

Mar 03, 2024

A novel technique for solving unsteady three

Scientific Reports volume 13, Article number: 13241 (2023) Cite this article 164 Accesses 3 Altmetric Metrics details The motion of the fluid due to the swirling of a disk/sheet has many applications

Scientific Reports volume 13, Article number: 13241 (2023) Cite this article

164 Accesses

3 Altmetric

Metrics details

The motion of the fluid due to the swirling of a disk/sheet has many applications in engineering and industry. Investigating these types of problems is very difficult due to the non-linearity of the governing equations, especially when the governing equations are to be solved analytically. Time is also considered a challenge in problems, and times dependent problems are rare. This study aims to investigate the problem related to a transient rotating angled plate through two analytical techniques for the three-dimensional thin film nanomaterials flow. The geometry of research is a swirling sheet with a three-dimensional unsteady nanomaterial thin-film moment. The problem's governing equations of the conservation of mass, momentum, energy, and concentration are partial differential equations (PDEs). Solving PDEs, especially their analytical solution, is considered a serious challenge, but by using similar variables, they can be converted into ordinary differential equations (ODEs). The derived ODEs are still nonlinear, but it is possible to approximate them analytically with semi-analytical methods. This study transformed the governing PDEs into a set of nonlinear ODEs using appropriate similarity variables. The dimensionless parameters such as Prandtl number, Schmidt number, Brownian motion parameter, thermophoretic parameter, Nusselt, and Sherwood numbers are presented in ODEs, and the impact of these dimensionless parameters was considered in four cases. Every case that is considered in this problem was demonstrated with graphs. This study used modified AGM (Akbari–Ganji Method) and HAN (Hybrid analytical and numerical) methods to solve the ODEs, which are the novelty of the current study. The modified AGM is novel and has made the former AGM more complete. The second semi-analytical technique is the HAN method, and because it has been solved numerically in previous articles, this method has also been used. The new results were obtained using the modified AGM and HAN solutions. The validity of these two analytical solutions was proved when compared with the Runge–Kutta fourth-order (RK4) numerical solutions.

In science, especially chemistry, condensate production from a cooling and saturated vapor is very substantial. Many researchers investigated this phenomenon under various circumstances. Sparrow and Gregg1 analyzed film condensation on a rotating plate on pure saturated steam. The centrifugal field associated with the rotation moves the condensate outward along the disc's surface without requiring gravitational forces. In this problem, the governing equations were solved numerically, and finally, results were given for heat transfer and condensate layer thickness, torque, temperature, and velocity profiles. Beckett et al.2 investigated the problem of laminar condensation on a swirling disk in a large volume of static vapor for low and high cooling rates on the disk surface. The governing equations were converted into a set of ODEs using similarity transformation and solved numerically, and solutions were compared via previously published results. Chary and Sarma3 considered the problem of vapor-to-liquid transition in the presence of constant axial suction at a permeable condensing surface. The governing equations were reduced into a set of ODEs. The Runge–Kutta numerical method was used to calculate the heat transfer coefficient, and limiting solutions for very thin condensate films were obtained. They determined that the heat transfer coefficient can be increased to any desired level by correctly selecting the suction parameter value. Attia and Aboul-Hassan4 investigated the transient motion of a viscous conducting fluid due to the swirling of an infinite, non-conducting, porous disk with a uniform magnetic field and the Hall effect. The governing equations were solved numerically, and the solution showed that including injection or suction from the disk surface in addition to the Hall flow gives interesting results. Bachok et al.5 investigated the transient boundary layer of a nanofluid flow on a permeable stretching/shrinking sheet. The governing equations are reduced into nonlinear ODEs and solved numerically. Freidoonimehr et al.6 studied a nanofluid unsteady MHD laminar free convection flow on a perpendicular sheet. The governing equations are reduced into the system of ODEs by a suitable similarity transformation and solved numerically with the RK4 method. Makinde et al.7 investigated the combined effects of thermal radiation, thermophoresis, Brownian motion, magnetic field and variable viscosity on boundary layer flow, heat and mass transfer of an electrically conducting nanofluid on a radially stretching convectively heated sheet. The governing equations transformed into a system of ODEs by using suitable similarity variables and solved numerically with RK4 method. Akbar et al.8 studied the two-dimensional non-transient incompressible viscous nanofluid flow on a stretching/shrinking plate. The governing PDEs were transformed into a set of ODEs by similarity variables and solved numerically via shooting method. Ramzan et al.9 studied the non-transient incompressible MHD nanofluid flow due to an infinite swirling disk with constant angular velocity, and the various velocity slip conditions are considered either. The governing equations were transformed into a set of nonlinear ODEs and solved numerically via the RK4 method. Alshomrani and Gul10 studied the nanofluid flow of a liquid film in a porous medium on a stretching sheet via the presence of velocity slip and thermal slip. The governing equations were transformed into a set of ODEs via suitable similarity variables and solved via the Homotopy Analysis Method (HAM). Gul and Sohail11 investigated the various Marangoni convection over a thin film flow on a stretching cylinder. The suitable similarity variables transformed the governing equations of this study into a set of ODEs and solved numerically via RK4 method. Ellahi12 investigated the MHD non-Newtonian nanofluid flow inside a pipe with an assumption that temperature of the pipe was higher than fluid temperature and also considered two particular temperature dependent viscosity models. The governing equations were transformed into a set of ODEs via suitable similarity variables and solved by the HAM. The analytical solutions of velocity field, the temperature distribution, and nano concentration have been derived. Khan and Pop13 investigated the steady two-dimensional laminar nano fluid flow and heat transfer arising from the stretching of a sheet and Brownian motion and thermophoresis was also considered in the problem. The governing equations were solved numerically after transforming the governing PDEs into a set of ODEs. Mustafa et al.14 studied the incompressible nanofluid flow, heat and mass transfer in a channel with presence of Brownian motion and thermophoresis effects. The governing equations were converted from PDEs into ODEs using suitable similarity transformation and then solved with both numerical method of RK4 and analytically with HAM. Akbar and Nadeem15 studied the two dimensional incompressible steady peristaltic flow of a nanofluid flow, heat, and mass transfer in an endoscope. The governing equations were transformed into dimensionless form and solved analytically via the Homotopy Perturbation method (HPM). Lakshmisha et al.16 investigated the three-dimensional transient laminar motion of a viscous, incompressible MHD fluid flow and heat transfer caused by the stretching of an infinite flat surface. The fluid was stationary at infinity, and the no-slip condition was imposed at the stretching surface in two lateral directions, where suction or injection can be applied. The governing equations were reduced into ODEs and solved by three different numerical methods. Wang17 investigated three-dimensional fluid flow due to the stretching of a sheet in two directions. The governing equations were reduced into a set of ODEs via suitable similarity transformation and then solved by the numerical method of RK4. Ahmad et al.18 investigated the problem of forced convection boundary layer nanofluid flow and heat transfer from a stationary semi-infinite flattish sheet and another problem similar to the previous one, but this time the flat sheet was not stationary. The governing equations were converted into a set of ODEs by a transformation and then solved with the numerical method of RK4. Chamkha et al.19 investigated the problem of boundary-layer nanofluid flow, heat and mass transfer on a dynamic porous media in the presence of magnetic field, heat generation or absorption, thermophoresis, Brownian motion, and suction or injection effects. The governing equations were reduced into a system of ODEs and solved numerically via the finite difference method (FDM). Kandasamy et al.20 studied the problem of three-dimensional unsteady laminar nanofluid flow, heat, and mass transfer due to the stretchy perpendicular sheet with changing stream conditions in the presence of Brownian motion and thermophoresis effects. The governing equations were reduced into a system of coupled nonlinear ODEs and solved numerically with the Oberbeck–Boussinesq approximation. Berkan et al.21 studied the problem of intransient three-dimensional condensation film over an angled swirling disk. The governing equations were reduced into a set of ODEs via transformation and solved analytically with AGM. The results were compared with the previously published studies. Mirgolbabaee et al.22 studied a two-dimensional intransient MHD laminar flow of fluid along parallel porous walls in which fluid is uniformly injected or removed. The governing equations were reduced into a set of ODEs via a similarity transformation and solved analytically. Jalili et al.23 studied the impacts of angled Lorentz body force and changing viscosity for the flow of non-Newtonian Williamson nanofluid over a stretchy sheet. The governing equations were transformed into ODEs via similarity variables and solved analytically. Jalili et al.24 studied the flow of an intransient two-dimensional MHD nanofluid over a semi-infinite stretchy flat plate. The governing equations were reduced into a set of ODEs and solved analytically. Jalili et al.25 investigated the problem of two-dimensional steady boundary layer micropolar ferrofluid flow and heat transfer due to the constricting plate with presence of thermal radiation and transverse magnetic field. The governing equations reduced into system of ODEs and solved analytically and numerically. Jalili et al.26 proposed the Hybrid Analytical and Numerical method (the HAN method) for solving a the problem of viscous, incompressible, laminar axisymmetric flow of a micropolar fluid with presence of magnetic field between two stretchable disks. The governing equations were reduced into ODEs by similarity variables and solved analytically. Jalili et al.27,28 used the same method of HAN in two other studied either. Many problems29,30,31,32,33,34,35,36 related to fluid mechanics were studied and used the similarity transformation to convert the PDEs into ODEs but they solved numerically, Meanwhile, the modified HAN or AGM method had the potential to solve these problems analytically. The novelty of this article is that these two methods were used and the analytical answer was obtained.

This paper investigates heat and mass transfer in a transient swirling angled sheet analytically with two techniques for a 3D thin nanomaterial film flow. The semi-analytical methods used in this study are modified AGM and HAN methods. The modified AGM is novel and has made the former AGM more complete. The second semi-analytical technique is the HAN method, which is applied because Zeeshan et al. already solved the numerical solution to this problem37. The results from these two semi-analytical solutions were compared with the previously published RK4 solution.

The geometry of the study is a swirling sheet with a three-dimensional transient nanomaterial thin-film moment, as illustrated in Fig. 1. The plate swirls with the angular velocity of \(\Omega ,\) and the inclined plate has an angle of \(\beta\) with the horizon. The nanomaterial thickness of the sheet is indicated by \(h\), and the speed of the sprayed fluid is denoted by \(W\). The terminal effect is neglected because the thickness of the fluid film in comparison to the radius of the disc is not thick enough. The gravitation force exists, and it is denoted by \(\overline{g}\) , and its direction is illustrated in the following figure. The temperature of the film surface is denoted by \({T}_{0}\). The temperature of the inclined swirling surface is denoted by \({T}_{w}\). The concentration of the film surface is denoted by \({C}_{0}\). The concentration of the inclined swirling surface is denoted by \({C}_{w}\).

The geometry of the problem.

The thickness of the fluid film is very thin, and the pressure at the surface of the surface is denoted by \({p}_{0}\) and it is just a function of \(z\). The viscous dissipation function in the energy equation is negligible. The governing equations of the problem are as follows2,3,5,6,8,37:

The equation of conservation of mass:

The equation of conservation of momentum in \(x\) direction:

The equation of conservation of momentum in \(y\) direction:

The equation of conservation of momentum in \(z\) direction:

The equation of conservation of energy:

The equation of conservation of concentration:

where \(D/Dt\) denotes the total derivative to the variable of time, \(\nabla\) is the gradient operator, \(u\), \(v,\) and \(z\) are the velocities in the \(x\), \(y,\) and \(z\) directions, respectively, \({\nabla }^{2}\) is the Laplacian operator, \(T\) is the temperature, \(C\) is the concentration, \({\rho }_{f}\), is the density of the base fluid, \(\mu\) is the dynamic viscosity, \(\alpha\) is the thermal diffusivity, \({c}_{p}\), is the specific heat capacity at a constant pressure of nanofluid, \({\left(\rho {c}_{p}\right)}_{p}/{\left(\rho {c}_{p}\right)}_{f}\) , is the ratio of nanoparticles' heat capacity to the base fluid heat capacity, \({D}_{B}\) is the Brownian diffusion coefficient, and \({D}_{T}\), is the thermophoretic diffusion coefficient.

The corresponding boundary conditions of Eqs. (1)–(6) are as follows:

The similarity transformations are considered as follows8,11,37:

Here, \(\nu\) is the kinematic viscosity, \(\theta\) is the dimensionless temperature, and \(\phi\) is the dimensionless concentration. The similarity variables of Eq. (8) can be substituted in Eqs. (2)–(6) for converting a system of nonlinear PDEs into a system of nonlinear dimensionless coupled ODEs:

Substituting the similarity variables of Eq. (8) in Eq. (7) will be as follows:

where \(Pr\) is the Prandtl number, \(Sc\) is the Schmidt number, \(Nb\) is the Brownian motion parameter, \(S\) is the parameter that depends on the angular velocity of the rotating surface, and \(Nt\) is the thermophoretic parameter, which is defined as37,38:

The constant normalized thickness of \(\delta\) is as follows37:

The dimensionless Nusselt and Sherwood numbers are as follows37:

Jalili et al.26,27,28 developed the HAN method for approximating an analytical solution for a differential equation. In this part, the explanation of the HAN method is as follows:

Consider an ODE of the m \(th\) order as follow:

Equation (20) is a nonlinear differential equation, and \(\Gamma\) is the function of \(\zeta\) and its derivatives to \(\upxi\). The parameter \(\zeta\) is the function of the independent variable \(\upxi\). The derivatives of the function \(\zeta \left(\upxi \right)\) with respect to \(\upxi\) at \(\upxi =0\) and \(\upxi =L\) are denoted as follows:

The solution of Eq. (20) is considered as follows:

Here, \({a}_{0}\), \({a}_{1}\), …, \({a}_{n}\) are \(n+1\) constant coefficients which \(n>m\). By solving a system of \(n+1\) unknowns and \(n+1\) equations, constant coefficients will be determined. The boundary conditions of the problem can construct some of these equations as follows:

The constructed equations from boundary conditions of the problem as they can be seen in Eqs. (23), (24) are limited because we assume the value of \(n\) is higher than \(m\) earlier in this methodology. But more boundary equations are needed, and a numerical method (no matter which numerical method and no matter what kind of software package) can approximate these additional boundary conditions for making the remaining needed equations. So, the new approximated boundary conditions are as follows:

For instance, the following equations are constructed from approximated boundary conditions of Eq. (25):

According to Eqs. (26)–(29), it can be derived as many equations as possible are needed to create a system with \(n+1\) equations and \(n+1\) unknowns. The limitation of the HAN method is just in the numerical method that is used, and this means that if no numerical method could solve a problem, the HAN method could not be used because this method seriously needs a numerical solution. To summarize the mentioned method in a more compact form, the following Fig. 2, the flow chart is presented for the HAN method:

The flow chart of the HAN method.

For applying the HAN method, let us assume the following functions are the semi-analytical solutions of Eqs. (9)–(14):

Based on Eq. (30), there are 43 unknown coefficients, and 43 equations are needed to obtain them. Equation (15) makes only 13 equations, and the remaining 30 must be made numerically. This study used the numerical solution of Zeeshan et al.37. Finally, according to Table 1, the system of ODEs of Eqs. (9)–(14) for 4 cases can be solved by calculating the system of 46 equations and 46 unknowns and the solutions of Eqs. (9)–(14) for all available cases in Table 1 are as follows:

Solutions of case 1 where \(Pr=6.6\), \(Nt=0.2\), \(Nb=0.2\), \(Sc=2.0\), \(S=0.0\), \(\delta =1.0\) are are demostrated in Eqs. (31)–(36) as follows:

Solutions of case 2 where \(Pr=6.7\), \(Nt=0.4\), \(Nb=0.4\), \(Sc=4.0\), \(S=0.3\), \(\delta =1.0\) are are demostrated in Eqs. (37)–(42) as follows:

Solutions of case 3 where \(Pr=7.1\), \(Nt=0.6\), \(Nb=0.6\), \(Sc=6.0\), \(S=0.5\), \(\delta =1.0\) are are demostrated in Eqs. (43)–(48) as follows:

Solutions of case 4 where \(Pr=7.3\), \(Nt=0.8\), \(Nb=0.8\), \(Sc=8.0\), \(S=0.6\), \(\delta =1.0\) are are demostrated in Eqs. (49)–(54) as follows:

The Akbari–Ganji Method was developed for solving nonlinear differential equations analytically. This method has solved many problems21,22,23,24,25,39,40 for which no exact analytical method exists. This paper introduces the modification of this method due to needing more accurate solutions.

To explain the main idea of modified AGM, the general form of the m \(th\) order differential equation is assumed as:

With boundary conditions:

To solve Eq. (55), we can consider the answer as the following polynomial of degree \(n\) with unknown constant coefficients:

Here, \({a}_{0}\), \({a}_{1}\), …, \({a}_{n}\) are \(n+1\) constant coefficients which \(n>m\). By solving a system of \(n+1\) unknowns and \(n+1\) equations, constant coefficients will be determined. The boundary conditions of the problem can construct some of these equations as follows:

The constructed equations from boundary conditions of the problem as they can be seen in Eqs. (58), (59) are limited because we assume the value of \(n\) is higher than \(m\) earlier in this methodology. But more equations are needed to construct a system of \(n+1\) unknowns and \(n+1\) equations. So, the remaining equations can be made by substituting Eq. (57) in Eq. (55) as follows:

So, it can be derived as many equations as possible from Eqs. (60)–(62) to construct a system of \(n+1\) unknowns and \(n+1\) equations. Finally, series constant coefficients and, thus, the solution to the problem will be determined by solving the equations. Unlike the HAN method, AGM does not depend on the numerical solution and is more independent, but the limitation of this method is that the more nonlinear the problem, the more difficult it is to solve with the AGM method. To summarize the mentioned method in a more compact form, the following Fig. 3, the flow chart is presented for the modified AGM:

The flow chart of the modified AGM.

In this part, Eqs. (9)–(14) are solved with the modified Akbari–Ganji method for cases (1, 2) according to Table 1. For applying the Modified Akbari–Ganji Method, the following functions are to be assumed the semi-analytical solutions of Eqs. (9)–(14) for cases (1, 2):

Based on Eq. (63), there are 31 unknown coefficients, and 43 equations are needed to obtain them. Equation (15) makes only 13 equations, and the remaining 18 equations must be made through Eq. (60), and the solutions for cases (1, 2) are as follows:

Solutions of case 1 where \(Pr=6.6\), \(Nt=0.2\), \(Nb=0.2\), \(Sc=2.0\), \(S=0.0\), \(\delta =1.0\) are are demostrated in Eqs. (64)–(69) as follows:

Solutions of case 2 where \(Pr=6.7\), \(Nt=0.4\), \(Nb=0.4\), \(Sc=4.0\), \(S=0.3\), \(\delta =1.0\) are are demostrated in Eqs. (70)–(75) as follows:

But it can reach more accurate solutions by increasing \(n\) in assumed functions. The following functions are the semi-analytical solutions of Eqs. (9)–(14) for cases (3, 4):

Based on Eq. (76), 49 unknown coefficients and 49 equations are needed to obtain them. Since the number of constant coefficients has increased, the number of equations that we must create to make a system of \(n\) equations and \(n\) unknowns increases, and in addition to Eqs. (60), (61) should also be used. Equation (15) makes only 13 equations, and the remaining 36 equations must be made through Eqs. (60), (61) and the solutions for cases (3, 4) are as follows:

Solutions of case 3 where \(Pr=7.1\), \(Nt=0.6\), \(Nb=0.6\), \(Sc=6.0\), \(S=0.5\), \(\delta =1.0\) are are demostrated in Eqs. (77)–(82) as follows:

Solutions of case 4 where \(Pr=7.3\), \(Nt=0.8\), \(Nb=0.8\), \(Sc=8.0\), \(S=0.6\), \(\delta =1.0\) are demostrated in Eqs. (83)–(88) as follows:

Heat and mass transfer in an unsteady rotating inclined plane has been investigated using 3D thin film nanomaterials flow. The solutions were obtained using the modified AGM and HAN methods. These two analytical solutions' validity was proved when compared with the Zeeshan et al.37 Runge–Kutta fourth-order (RK4) numerical solutions. Figures 4, 5, 6, 7, 8, 9, 10, 11 and 12 show the accuracy of the modified AGM and HAN results. The following table relates to the four cases considered in this study.

The impact of different cases on \(f\left(\upxi \right)\).

The impact of different cases on \(\mathrm{g}\left(\upxi \right)\).

The impact of different cases on \(h\left(\upxi \right)\).

The impact of different cases on \(k\left(\upxi \right)\).

The impact of different cases on \(\uptheta \left(\upxi \right)\).

The impact of different cases on \(\upphi \left(\upxi \right)\).

The impact of different cases on \(f^{\prime}\left(\upxi \right)\).

The impact of different cases on \(\uptheta ^{\prime}\left(\upxi \right)\).

The impact of different cases on \(\upphi ^{\prime}\left(\upxi \right)\).

As the impact of four cases is shown in Figs. 4, 5, 6, 7, 8, 9 and variations of Sherwood number, heat transmission, and radial velocity profiles are illustrated in Figs. 10, 11 and 12. In this study, the average value of some arbitrary scalar function of \(\upchi \left(\upxi \right)\), is defined as follows:

where the \(a\) and \(b\) are integer numbers. According to Eq. (89), the average values of \(f\left(\upxi \right)\), \(f^{\prime}\left(\upxi \right)\),\(\mathrm{g}\left(\upxi \right)\), \(k\left(\upxi \right)\), \(h\left(\upxi \right)\), \(\uptheta \left(\upxi \right)\), \(\upphi \left(\upxi \right)\), \(\uptheta {^{\prime}}\left(\upxi \right)\), and \(\upphi{ ^{\prime}}\left(\upxi \right)\) are denoted by \({\overline{f} }_{avg}\), \({\overline{f} }_{avg}^{\prime}\),\({\overline{\mathrm{g}} }_{avg}\), \({\overline{k} }_{avg}\), \({\overline{h} }_{avg}\), \({\overline{\uptheta } }_{avg}\), \({\overline{\upphi } }_{avg}\), \({\overline{\uptheta } }_{avg}^{\prime}\), and \({\overline{\upphi } }_{avg}^{\prime}\) , respectively. According to Figs. 4, 5, 6, 7, 8, 9, 10, 11 and 12, when cases (1–4) related to Table 1 occur, the average values of \({\overline{f} }_{avg}\), \({\overline{f} }_{avg}^{\prime}\), \({\overline{\mathrm{g}} }_{avg}\), \({\overline{k} }_{avg}\), \({\overline{h} }_{avg}\), \({\overline{\uptheta } }_{avg}\), \({\overline{\upphi } }_{avg}\), \({\overline{\uptheta } }_{avg}^{\prime}\), and \({\overline{\upphi } }_{avg}^{\prime}\) , are given in the following tables:

In Tables 2, 3 and 4, some average results increased when the conditions changed from case 1 to case 4, and some decreased. The decrease or increase of these values is calculated using the following relationship:

where \(\Omega\) is the amount of percentage increase or decrease values, \({\mathrm{ Z}}_{2}\) is the second value and \({\mathrm{ Z}}_{1}\) is the first value. According to Table 2, when the constant coefficients of cases in Table 1 change from case 1 to case 4, the average values of the results from Ref.37 change respectively. When the conditions change from case 1 to case 4, \({\overline{f} }_{avg}\) , will decrease by 34.71362167%, \({\overline{f} }_{avg}^{\prime}\) will decrease by 35.78465385%, \({\overline{\mathrm{g}} }_{avg}\), will decrease by 8.197818397%, \({\overline{k} }_{avg}\) will decrease by 9.504907967%, \({\overline{h} }_{avg}\) will increase by 23.38976383%, \({\overline{\uptheta } }_{avg}\) will increase by 9.441723369%, \({\overline{\uptheta } }_{avg}^{\prime}\) will increase by 7.742384880%, and \({\overline{\upphi } }_{avg}^{\prime}\) will decrease by 9.974419462% but according to Table 2, \({\overline{\upphi } }_{avg}\) will increase by 1.541136126% when it changes from case 1 to case 2 but, \({\overline{\upphi } }_{avg}\) , will decrease by 1.839033024% when it changes from case 2 to case 4. According to Table 3, when the constant coefficients of cases in Table 1 change from case 1 to case 4, the average values of the results from Modified AGM change respectively. When the conditions change from case 1 to case 4, \({\overline{f} }_{avg}\) , will decrease by 34.60665497%, \({\overline{f} }_{avg}^{\prime}\) will decrease by 35.68188954%, \({\overline{\mathrm{g}} }_{avg}\) , will decrease by 8.140379566%, \({\overline{k} }_{avg}\) will decrease by 9.388178325%, \({\overline{h} }_{avg}\) will increase by 23.82760391%, \({\overline{\uptheta } }_{avg}\) will increase by 9.259738964%, \({\overline{\uptheta } }_{avg}^{\prime}\) will increase by 7.367258115%, and \({\overline{\upphi } }_{avg}^{\prime}\) will decrease by 9.528160070% but according to Table 2, \({\overline{\upphi } }_{avg}\) will increase by 1.313715539% when it changes from case 1 to case 2 but, \({\overline{\upphi } }_{avg}\) , will decrease by 1.109559162% when it changes from case 2 to case 4. According to Table 4, when the constant coefficients of cases in Table 1 change from case 1 to case 4, the average values of the results from the HAN Method change respectively. When the conditions change from case 1 to case 4, \({\overline{f} }_{avg}\) , will decrease by 34.71362166%, \({\overline{f} }_{avg}^{\prime}\) will decrease by 35.78461613%, \({\overline{\mathrm{g}} }_{avg}\) , will decrease by 8.197818386%, \({\overline{k} }_{avg}\) will decrease by 9.504907931%, \({\overline{h} }_{avg}\) will increase by 23.38976385%, \({\overline{\uptheta } }_{avg}\) will increase by 9.441723405%, \({\overline{\uptheta } }_{avg}^{\prime}\) will increase by 7.592028101%, and \({\overline{\upphi } }_{avg}^{\prime}\) will decrease by 9.661365531% but according to Table 2, \({\overline{\upphi } }_{avg}\) will increase by 1.541136126% when it changes from case 1 to case 2 but, \({\overline{\upphi } }_{avg}\) , will decrease by 1.839033024% when it changes from case 2 to case 4.

This study investigates the problem of heat and mass transfer in an unsteady rotating inclined plane using 3D thin film nanomaterial flow. The governing equations were set PDEs, and by using suitable similarity transformation, the PDEs were reduced into a set of nonlinear ODEs. The ODEs in four cases were solved with two semi-analytical techniques of Modified AGM and HAN. The Modified AGM that is used in this study is a novel technique, and the novelty of current work is related to solving this problem analytically. Unlike the former AGM, the Modified Agm has solved the previous issues and can replace the previous method of AGM. The HAN Method is another semi-analytical method that transforms a numerical solution into an analytical one. Technically, if the numerical solution exists for some problem, then HAN Method can be applied to obtain an analytic solution. The results of the HAN solution are very close to the Numerical solutions of Zeeshan et al.37 when compared with the modified AGM, but at the same time, the modified AGM is not dependent on any numerical methods for approximating analytical solutions. So, this paper is concluded that:

A new semi-analytical is introduced by modifying the former AGM technique.

The exact analytic solutions were obtained through HAN Method.

The solutions of both analytical solutions were compared with previously published papers.

The results of both analytical solutions were presented quantitatively.

The Sherwood number of the film surface and inclined swirling surface will decrease as the Schmidt number increases and the angular velocity of the rotating surface decreases.

The Nusselt number of inclined swirling surfaces will increase as the Prandtl number increases and the angular velocity of the rotating surface decreases.

The Nusselt number of film surfaces will decrease as the Prandtl number increases and the angular velocity of the rotating surface decreases.

The datasets used and/or analysed during the current study available from the corresponding author on reasonable request.

Sparrow, E. M. & Gregg, J. L. A theory of rotating condensation. J. Heat Transf. 81(2), 113–119 (1959).

Article CAS Google Scholar

Beckett, P. M., Hudson, P. C. & Poots, G. Laminar film condensation due to a rotating disk. J. Eng. Math. 7(1), 63–73 (1973).

Article MATH Google Scholar

Chary, S. P. & Sarma, P. K. Condensation on a rotating disk with constant axial suction. ASME J. Heat Transf. 98, 682–684 (1976).

Article ADS Google Scholar

Attia, H. A. & Aboul-Hassan, A. L. Effect of Hall current on the unsteady MHD flow due to a rotating disk with uniform suction or injection. Appl. Math. Model. 25(12), 1089–1098 (2001).

Article MATH Google Scholar

Bachok, N., Ishak, A. & Pop, I. Unsteady boundary-layer flow and heat transfer of a nanofluid over a permeable stretching/shrinking sheet. Int. J. Heat Mass Transf. 55(7), 2102–2109 (2012).

Article CAS Google Scholar

Freidoonimehr, N., Rashidi, M. M. & Mahmud, S. Unsteady MHD free convective flow past a permeable stretching vertical surface in a nano-fluid. Int. J. Therm. Sci. 87, 136–145 (2015).

Article CAS Google Scholar

Makinde, O. D. et al. MHD flow of a variable viscosity nanofluid over a radially stretching convective surface with radiative heat. J. Mol. Liq. 219, 624–630 (2016).

Article CAS Google Scholar

Akbar, T. et al. Magnetohydrodynamics flow of nanofluid due to stretching/shrinking surface with slip effect. Adv. Mech. Eng. 9(12), 1687814017740266 (2017).

Article Google Scholar

Ramzan, M., Chung, J. D. & Ullah, N. Partial slip effect in the flow of MHD micropolar nanofluid flow due to a rotating disk—A numerical approach. Results Phys. 7, 3557–3566 (2017).

Article ADS Google Scholar

Alshomrani, A. S. & Gul, T. A convective study of Al2O3–H2O and Cu–H2O nano-liquid films sprayed over a stretching cylinder with viscous dissipation. Eur. Phys. J. Plus 132(11), 495 (2017).

Article Google Scholar

Gul, T. & Sohail, M. Marangoni liquid film scattering over an extending cylinder. Theor. Appl. Mech. Lett. 9(2), 106–112 (2019).

Article Google Scholar

Ellahi, R. The effects of MHD and temperature dependent viscosity on the flow of non-Newtonian nanofluid in a pipe: Analytical solutions. Appl. Math. Model. 37(3), 1451–1467 (2013).

Article MathSciNet MATH Google Scholar

Khan, W. A. & Pop, I. Boundary-layer flow of a nanofluid past a stretching sheet. Int. J. Heat Mass Transf. 53(11), 2477–2483 (2010).

Article CAS MATH Google Scholar

Mustafa, M. et al. Influence of wall properties on the peristaltic flow of a nanofluid: Analytic and numerical solutions. Int. J. Heat Mass Transf. 55(17), 4871–4877 (2012).

Article CAS Google Scholar

Akbar, N. S. & Nadeem, S. Endoscopic effects on peristaltic flow of a nanofluid. Commun. Theor. Phys. 56(4), 761 (2011).

Article ADS MATH Google Scholar

Lakshmisha, K.N., Venkateswaran, S., & Nath, G. Three-dimensional unsteady flow with heat and mass transfer over a continuous stretching surface. (1988).

Wang, C. The three-dimensional flow due to a stretching flat surface. Phys. Fluids 27(8), 1915–1917 (1984).

Article ADS MathSciNet MATH Google Scholar

Ahmad, S., Rohni, A. M. & Pop, I. Blasius and Sakiadis problems in nanofluids. Acta Mech. 218(3), 195–204 (2011).

Article MATH Google Scholar

Chamkha, A. J., Aly, A. M. & Al-Mudhaf, H. Laminar MHD mixed convection flow of a nanofluid along a stretching permeable surface in the presence of heat generation or absorption effects. Int. J. Microscale Nanoscale Therm. Fluid Trans. Phenomena 2(1), 51–70 (2011).

Google Scholar

Kandasamy, R., Loganathan, P. & Arasu, P. P. Scaling group transformation for MHD boundary-layer flow of a nanofluid past a vertical stretching surface in the presence of suction/injection. Nucl. Eng. Des. 241(6), 2053–2059 (2011).

Article CAS Google Scholar

Berkan, S., Hoseini, S. R. & Ganji, D. D. Analytical investigation of steady three-dimensional problem of condensation film on inclined rotating disk by Akbari–Ganji’s method. Propuls. Power Res. 6(4), 277–284 (2017).

Article Google Scholar

Mirgolbabaee, H., Ledari, S. T. & Ganji, D. D. Semi-analytical investigation on micropolar fluid flow and heat transfer in a permeable channel using AGM. J. Assoc. Arab Univ. Basic Appl. Sci. 24, 213–222 (2017).

Google Scholar

Jalili, B. et al. Thermal analysis of Williamson fluid flow with Lorentz force on the stretching plate. Case Stud. Therm. Eng 39, 102374 (2022).

Article Google Scholar

Jalili, B. et al. Characteristics of ferrofluid flow over a stretching sheet with suction and injection. Case Stud. Therm. Eng. 14, 100470 (2019).

Article Google Scholar

Jalili, B. et al. effect of magnetic and boundary parameters on flow characteristics analysis of micropolar ferrofluid through the shrinking sheet with effective thermal conductivity. Chin. J. Phys. 71, 136–150 (2021).

Article MathSciNet CAS Google Scholar

Jalili, P. et al. Heat transfer analysis in cylindrical polar system with magnetic field: A novel Hybrid Analytical and Numerical Technique. Case Stud. Therm. Eng. 40, 102524 (2022).

Article Google Scholar

Jalili, P. et al. Study of nonlinear radiative heat transfer with magnetic field for non-Newtonian Casson fluid flow in a porous medium. Results Phys. 48, 106371 (2023).

Article Google Scholar

Jalili, P. et al. The HAN method for a thermal analysis of forced non-Newtonian MHD Reiner-Rivlin viscoelastic fluid motion between two disks. Heliyon 9, e17535 (2023).

Article PubMed PubMed Central Google Scholar

Zeeshan, A. et al. Simultaneous effects of slip and wall stretching/shrinking on radiative flow of magneto nanofluid through porous medium. J. Magn. 23(4), 491–498 (2018).

Article Google Scholar

Yousif, M. A. et al. Numerical study of momentum and heat transfer of MHD Carreau nanofluid over an exponentially stretched plate with internal heat source/sink and radiation. Heat Transf. Res. 50(7), 649–658 (2019).

Article Google Scholar

Bhatti, M. et al. Simultaneous influence of thermo-diffusion and diffusion-thermo on non-Newtonian hyperbolic tangent magnetised nanofluid with Hall current through a nonlinear stretching surface. Pramana 93, 1–10 (2019).

Article CAS Google Scholar

Nazir, U. et al. Finite element analysis for thermal enhancement in power law hybrid nanofluid. Front. Phys. 10, 996174 (2022).

Article Google Scholar

Sohail, M. et al. A study of triple-mass diffusion species and energy transfer in Carreau–Yasuda material influenced by activation energy and heat source. Sci. Rep. 12(1), 10219 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Sohail, M. et al. Finite element analysis for ternary hybrid nanoparticles on thermal enhancement in pseudo-plastic liquid through porous stretching sheet. Sci. Rep. 12(1), 9219 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Nazir, U. et al. Significant production of thermal energy in partially ionized hyperbolic tangent material based on ternary hybrid nanomaterials. Energies 14(21), 6911 (2021).

Article CAS Google Scholar

Nazir, U. et al. Thermal and mass species transportation in tri-hybridized Sisko martial with heat source over vertical heated cylinder. Int. Commun. Heat Mass Transf. 134, 106003 (2022).

Article CAS Google Scholar

Zeeshan, A. et al. Numerical computation of 3D Brownian motion of thin film nanofluid flow of convective heat transfer over a stretchable rotating surface. Sci. Rep. 12, 2708 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Rehman, K. U., Ali Khan, A. & Malik, M. Magneto-nanofluid numerical modelling of chemically reactive Eyring-Powell fluid flow towards both flat and cylindrical an inclined surfaces: A comparative study. AIP Adv. 7(6), 065103 (2017).

Article ADS Google Scholar

Attar, M. A. et al. Analytical solution of fractional differential equations by Akbari–Ganji’s method. Partial Differ. Equ. Appl. Math. 6, 100450 (2022).

Article Google Scholar

Jalili, B. et al. Thermal analysis of boundary layer nanofluid flow over the movable plate with internal heat generation, radiation, and viscous dissipation. Case Stud. Therm. Eng. 49, 103203 (2023).

Article Google Scholar

Download references

Department of Mechanical Engineering, North Tehran Branch, Islamic Azad University, Tehran, Iran

Payam Jalili, Ali Ahmadi Azar & Bahram Jalili

Department of Mechanical Engineering, Babol Noshirvani University of Technology, P.O. Box 484, Babol, Iran

Davood Domiri Ganji

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

A.A.A. and P.J. wrote the main manuscript text and B.J. and D.D.G. prepared figures. All authors reviewed the manuscript.

Correspondence to Bahram Jalili or Davood Domiri Ganji.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Jalili, P., Ahmadi Azar, A., Jalili, B. et al. A novel technique for solving unsteady three-dimensional brownian motion of a thin film nanofluid flow over a rotating surface. Sci Rep 13, 13241 (2023). https://doi.org/10.1038/s41598-023-40410-3

Download citation

Received: 19 March 2023

Accepted: 09 August 2023

Published: 15 August 2023

DOI: https://doi.org/10.1038/s41598-023-40410-3

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.