In this study, a simple and high accurate series-based method called Differential Transformation Method (DTM) is used for solving the coupled nonlinear differential equations in fluids mechanic problems. The concept of the DTM is briefly introduced, and its application on two different cases, natural convection of a non-Newtonian nanofluid between two vertical plates and Newtonian nanofluid flow between two horizontal plates, has been studied. DTM results are compared with those obtained by a numerical solution (Fourth-order Runge–Kutta) to show the accuracy of the proposed method. Results reveal that DTM is very effective and convenient which can achieve more reliable results compared to other analytical methods in solving some engineering and sciences problems.
Coupled nonlinear differential equations; Differential Transformation Method; Natural convection; Newtonian and non-Newtonian nanofluids
The nonlinear problems and phenomena can be modeled by ordinary or partial nonlinear differential equations to find their physical behavior in certain cases. Most of these described physical and mechanical problems have a set of coupled nonlinear differential equations. For example, heat transfer by natural convection which commonly occurs in many physical and engineering applications such as geothermal systems, chemical catalytic reactors, and heat exchangers has a system of coupled nonlinear differential equations. One of the simple and reliable methods for solving the system of coupled nonlinear differential equations is analytical solution which recently is widely used in solving many problems [1], [2], [3], [4], [5], [6], [7] and [8]. Ziabakhsh and Domairry [9] have studied the natural convection of a non-Newtonian fluid between two infinite parallel vertical flat plates by Homotopy Analysis Method (HAM) which is a general form of the HPM.
An analytical method called Differential Transformation Method (DTM) has also been introduced for nonlinear problem solution. DTM is an analytical method based on Taylor series. It constructs an analytical solution in the form of a polynomial. Actually, DTM is different from the traditional high order Taylor series method, which requires symbolic calculation of the necessary derivatives of the data functions. This method was firstly applied in the engineering field by Zhou [10] and its abilities attracted many authors to use this method for solving the nonlinear equations. DTM or improved shape forms of it (Ms-DTM and Hybrid-DTM) can solve any coupled nonlinear equations set. They can also be used for three dimensional problems. Ghafoori et al. [11] solved a nonlinear oscillation equation by DTM and they revealed that DTM results are more accurate than HPM and VIM (Variation Iteration Method). Yahyazadeh et al. [12] evaluated the natural convection flow of a nanofluid over a linearly stretching sheet in the presence of magnetic field using DTM. Also, Biazar and Eslami [13] considered DTM to solve the quadratic Riccati differential equation and the derived results by DTM were compared with the results of HAM and ADM (adomian decomposition method). They showed that DTM results are more effective and favorable than HAM and ADM results. Gokdogan et al. [14] acquired an approximate analytical solution of the chaotic Genesio system by a modified DTM. Ayaz [15] studied two-dimensional differential transform method for solving an initial value problem in partial differential equations (PDEs). In another application of DTM, the natural frequencies and critical flow velocities of pipes conveying fluid with several typical boundary conditions were obtained by Ni et al. [16]. Although many analytical methods have been presented in the literature, DTM is expected to find more valuable applications due to the following advantages:
through h-curves against Homotopy Analysis method (HAM).
This paper aimed to apply DTM for system of coupled nonlinear equations and testify the reliability of DTM by comparing its results with HPM. For showing the accuracy of DTM in described problems, two cases of problems, natural convection of a non-Newtonian fluid between two vertical plates and Newtonian fluid flow analysis between two horizontal plates, were solved. As an important result, it was found that the DTM results are more accurate than those obtained by HPM in some areas. After this verification, the effects of some physical parameters were analyzed to show the efficiency of DTM for this type of the problems.
The study of the fluid flow between parallel plates has many applications in industries which motivated the researchers to investigate. Ziabakhsh and Domairry [9] investigated the natural convection of non-Newtonian fluid between vertical parallel plates using HAM, and based on their work, Hatami and Ganji [18] studied the effect of sodium alginate (SA) non-Newtonian nanofluid between parallel plates. Also, Rajagopal and Na [19] presented a numerical solution for the natural convection of non-Newtonian fluids between vertical plates. In the following section, two cases of problems in fluid flow between parallel plates with coupled nonlinear differential equations will be introduced. It is aimed to investigate the efficiency of DTM in comparison with numerical solution.
In first case the heat transfer analysis in the unsteady two-dimensional squeezing nanofluid flow between the infinite parallel plates (Fig. 1, left). The two plates are placed at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle z=\pm l{\left(1-\alpha t\right)}^{1/2}=\pm h(t)} . For Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha >0} , the two plates are squeezed until they touch where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle t=1/\alpha }
and for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \alpha <0} the two plates are separated. The viscous dissipation effect, the generation of heat due to friction caused by shear in the flow, is retained. This effect is quite important in the case when the fluid is largely viscous or flowing at a high speed. This behavior occurs at high Eckert number (≫1). The symmetric nature of the flow is adopted. The fluid is a water based nanofluid containing Cu (copper) nanoparticles. The nanofluid is a two component mixture with the following assumptions: incompressible; no-chemical reaction; negligible viscous dissipation; negligible radiative heat transfer; nano-solid-particles and the base fluid are in thermal equilibrium and no slip occurs between them. The governing equations for momentum and energy in unsteady two dimensional flow of a nanofluid are as follows [17]:
|
(1) |
|
(2) |
|
(3) |
|
(4) |
Here u and v are the velocities in the x and y directions respectively, T is the temperature, P is the pressure, and effective density Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \left({\rho }_{nf}\right)} , the effective dynamic viscosity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \left({\mu }_{nf}\right)} , the effective heat capacity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\left(\rho C_p\right)}_{nf}}
and the effective thermal conductivity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle k_{nf}} of the nanofluid are defined as follows [7]:
|
(5) |
The relevant boundary conditions are as follows:
|
(6) |
We introduce these parameters:
|
(7) |
Substituting the above variables into (2) and (3) and then eliminating the pressure gradient from the resulting equations give
|
(8) |
Using (7), Eqs. (3) and (4) reduces to the following differential equations:
|
(9) |
Here Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle A_2}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle A_3} are constants given by
|
(10) |
With these boundary conditions:
|
(11) |
where S is the squeeze number, Pr is the Prandtl number and Ec is the Eckert number, which are defined as follows:
|
(12) |
|
Fig. 1. Schematic of the two cases of problems with coupled nonlinear equations in fluids mechanic. |
The second case (Fig. 1, right) consists of two vertical flat plates separated by a distance of 2b . A non-Newtonian fluid flows between them due to natural convection. The walls at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=+b}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=-b} are held at constant temperatures Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle T_2} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle T_1}
, respectively, where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle T_1>T_2} . This difference in temperature causes the fluid near the wall at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=-b}
to rise and the fluid near the wall at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x=+b} to fall. The fluid is a non-Newtonian Sodium Alginate (SA) based nanofluid containing Cu and Ag nanoparticles. It is assumed that the base fluid and the nanoparticles are in thermal equilibrium and no slip occurs between them. The effective density Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\rho }_{nf}}
, the effective dynamic viscosity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\mu }_{nf}} , the heat capacitance Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle {\left(\rho C_p\right)}_{nf}}
and the thermal conductivity Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle k_{nf}} of the nanofluid are given as follows [7]:
|
(13) |
|
(14) |
|
(15) |
|
(16) |
Here, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi }
is the solid volume fraction. By definition following are the similarity variables:
|
(17) |
Under these assumptions and following the nanofluid model proposed by Maxwell–Garnett (MG) model [9], the Navier–Stokes and energy equations can be reduced to the following pair of ordinary differential equations:
|
(18) |
|
(19) |
where Prandtl number (Pr), Eckert number (Ec ), dimensionless non-Newtonian viscosity (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta } ) and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle A_1}
have the following forms:
|
(20) |
|
(21) |
The appropriate boundary conditions are as follows:
|
(22) |
For understanding the concept of DTM, we suppose that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)}
is an analytic function in domain D , and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle t=t_i} represents any point in the domain. The function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)} is then represented by one power series whose center is located at Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle t_i}
. The Taylor series expansion function of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)}
is in form of
|
(23) |
The Maclaurin series of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)}
can be obtained by taking Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle t_i=0} in Eq. (23) as follows:
|
(24) |
As explained in Ref. [10], the differential transformation of the function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)}
is defined as follows:
|
(25) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)}
represents the transformed function and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)} is the original function. The differential spectrum of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)} is confined within the interval Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle t\in [0\mbox{,}H]}
, where H is a constant value and it can be assumed as unity in common DTM, but for multi-step DTM it should be considered as the length of t steps. When Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle H\rightarrow \infty }
the methods no longer work and Pade approximation should be applied. More information about H can be found in Ref. [20]. The differential inverse transform of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)} is defined as follows:
|
(26) |
It is clear that the concept of differential transformation is based upon the Taylor series expansion. The values of function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)}
at values of argument k are referred to as discrete, i.e. Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(0)} is known as the zero discrete, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(1)} as the first discrete, etc. The more the discrete available, the more precise it is possible to restore the unknown function. The function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)} consists of the T -function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)}
, and its value is given by the sum of the T -function with Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle (t/H)k}
as its coefficient. In real applications, with the right choice of constant H, the larger values of argument k will lead to rapid reduction of the discrete of spectrum. The function Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)} is expressed by a finite series and Eq. (26) can be written as follows:
|
(27) |
where n is the number of statements of the DTM. Generally, by increasing n, DTM accuracy will be increased, but it depends on H and time step’s length [20]. Actually, H depends on the time steps which for classical DTM is considered unity, but for Multi step DTM or Ms-DTM it is equal to the time step value. Some important mathematical operations performed by differential transform method are listed in Table 1. As described before, one of the most advantages of DTM is its independence of small parameter, linearization or perturbation. Also, it does not need to determine auxiliary function, auxiliary parameter or suitable initial guess against other analytical methods.
Origin function | Transformed function |
---|---|
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)=\alpha f(x)\pm \beta g(t)} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)=\alpha F(k)\pm \beta G(k)} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)=\frac{d^mf(t)}{{dt}^m}} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)=\frac{\left(k+m)!F(k+m\right)}{k!}} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)=f(t)g(t)} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)={\sum }_{l=0}^kF(l)G(k-l)} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)=t^m} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)=\delta (k-m)=\begin{array}{ll} 1\mbox{,} & \mbox{if}\quad k=m\\ 0\mbox{,} & \mbox{if}\quad k\quad \not =\quad m\mbox{.} \end{array}} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)=exp(t)} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)=\frac{1}{k!}} |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)=sin(\omega t+\alpha )} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)=\frac{{\omega }^k}{k!}sin\left(\frac{k\pi }{2}+\right. }
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. \alpha \right) |
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle x(t)=cos(\omega t+\alpha )} | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle X(k)=\frac{{\omega }^k}{k!}cos\left(\frac{k\pi }{2}+\right. }
Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): \left. \alpha \right) |
Now Differential Transformation Method has been applied into the governing equations. Taking the differential transforms of Eqs. (8), (9) and (10) with respect to Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \chi }
and considering H = 1 gives:
|
(28) |
|
(29) |
|
(30) |
|
(31) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle F[k]}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Theta [k]} are the differential transforms of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f(\eta )\mbox{,}\theta (\eta )} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_1\mbox{,}a_2\mbox{,}a_3} are constants which can be obtained through boundary condition. This problem can be solved as follows:
|
(32) |
|
(33) |
The above process is continuous. By substituting Eqs. (32) and (33) into the main equation based on DTM, it can be obtained that the closed form of the solutions is as follows:
|
(34) |
|
(35) |
by substituting the boundary condition from Eq. (11) into Eqs. (34) and (35) in point Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \eta =1}
the values of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_1\mbox{,}a_2\mbox{,}a_3} can be obtained. By substituting obtained Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a_1\mbox{,}a_2\mbox{,}a_3} into Eqs. (34) and (35), the expression of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle F(\eta )} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Theta (\eta )} can be obtained. For example for Cu–water nanofluid when Pr = 6.2, Ec = 0.05, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta =0.1}
, S = 0.1 and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi =0.01}
following equations will be obtained:
|
(36) |
|
(37) |
Similarly, DTM is applied to Eqs. (18) and (19). Their transformed form will be as follows:
|
(38) |
The boundary conditions can be written as follows:
|
(39) |
In ith step, by solving the set of coupled Eq. (38), Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \overline{V}(i+2)}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Theta (i+2)} can be determined. For example using boundary condition and Eq. (22) second terms will be
|
(40) |
where Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \Theta }
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \overline{V}} represent the DTM transformed form of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta } and V respectively and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a\mbox{,}b\mbox{,}c} and d are unknown coefficients that after specifying Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta (X)} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle V(X)} and applying boundary condition (Eq. (22)) into it, will be determined. For example when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pr=\delta =Ec=1} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi =0.01} for Cu–SA following values were determined for Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle a\mbox{,}b\mbox{,}c} and d coefficients.
|
(41) |
Finally, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta (X)}
and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle V(X)} when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pr=\delta =Ec=1} and Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi =0.01} can be defined as follows:
|
(42) |
As seen in terms of above equation, although for increasing the accuracy, the number of statements has been increased, convergence of DTM is completely evident.
As has been discussed, Newtonian and non-Newtonian nanofluids can be found in many industrial applications [21], [22], [23], [24], [25], [26], [27] and [28]. Their analysis and treatment can be performed by well-known analytical and numerical methods. Some of these methods have been employed for solving the engineering problems [29], [30], [31] and [32]. In our study, efficiency of DTM for solving these kinds of problems is examined through two different problems as introduced in the previous sections. Our DTM solutions will be compared to numerical solution presented in the literature [19]. As mentioned before, in first case, heat transfer and nanofluid flow analysis in the unsteady squeezing nanofluids between parallel plates are studied using Differential Transformation Method (DTM) and results are compared with forth-order Runge–Kutta numerical method. Fig. 2 shows the results of DTM for solving the Eqs. (8) and (9) in different Eckert and squeeze numbers. As seen in these figures, DTM has a good agreement with numerical method in wide range of Ec and S numbers. Effect of Nanoparticles volume fraction on velocity and temperature profiles is shown in Fig. 3. Adding nanoparticles into the base fluid leads to increase in thermal boundary layer thickness while it has no significant effect on velocity boundary layer thickness. It is obvious that when nanoparticles were added into base fluid, heat transfer will increase due to their high thermal conductivity, and so temperature profiles will be decreased.
|
Fig. 2. Comparison of DTM and numerical results for (a) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta (\eta )} when Pr = 6.2, S = 0.1, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta =0.1} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi =0.01} and (b) Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle f(\eta )} when Pr = 6.2, Ec = 0.05, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta =0.1} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi =0.01} . |
|
Fig. 3. Effect of nanoparticles volume fraction (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi } ) on (a) temperature profile and (b) velocity profile when Pr = 6.2, Ec = 0.5, Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta =0.1} , S = 1.0. |
For the second case, the base fluid is considered as a non-Newtonian fluid containing sodium alginate (SA) and two types of nanoparticles namely silver (Ag) and copper (Cu) are added. For showing the efficiency of analytical applied method (DTM), Table 2 and Table 3 presented values and errors for Ag–SA, respectively. These tables also confirm the accuracy of DTM for solving such kinds of problems. The effect of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta }
number on non-dimensional velocity and temperature of the nanofluid is investigated through Fig. 4(a) and (b). These figures confirm that Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta } has no effect on temperature variations but it reduces velocity values. The same phenomenon was also observed by Ziabakhsh and Domairry [9]. Fig. 5 confirms that when nanoparticles volume fraction increases, velocity profiles increase but temperature values decrease due to increase in heat transfer. Fig. 5 shows the effect of nanoparticles volume fraction (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi }
) on velocity profile and temperature distribution for SA–TiO2 nanofluid when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Ec=\delta =1} , while Fig. 6 depictesthe effect of volume fraction of copper nanoparticles.
Table 2.
Velocity and temperature profile values by applied methods for Ag–SA nanofluid profiles when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pr=Ec=\delta =1} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi =0.01} .
X | V(X) | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta }
(X) | ||
---|---|---|---|---|
Numerical | DTM | Numerical | DTM | |
−1.0 | 0.00 | 2 × 10−11 | 0.50 | 0.50 |
0.45059686 | −0.9 | 0.013888702 | 0.01402428 | 0.45043148 |
−0.8 | 0.023620556 | 0.02368182 | 0.40071897 | 0.40113389 |
−0.7 | 0.029489465 | 0.02945437 | 0.35094362 | 0.35161056 |
−0.6 | 0.031888322 | 0.03182368 | 0.30114965 | 0.30202640 |
−0.5 | 0.031279886 | 0.03127150 | 0.25135346 | 0.25238090 |
−0.4 | 0.028164333 | 0.02827959 | 0.20155266 | 0.20267356 |
−0.3 | 0.023053394 | 0.02332968 | 0.15173412 | 0.15290389 |
−0.2 | 0.016454874 | 0.01690353 | 0.10188071 | 0.10307139 |
−0.1 | 0.008866496 | 0.00948289 | 0.05197647 | 0.05317557 |
0.0 | 0.000776279 | 0.00154951 | 0.002010267 | 0.00321593 |
0.1 | −0.007333249 | −0.00641485 | −0.04802206 | −0.04680801 |
0.2 | −0.014979179 | −0.01392846 | −0.09811697 | −0.09689681 |
0.3 | −0.021672314 | −0.02050956 | −0.14826378 | −0.14705093 |
0.4 | −0.026912911 | −0.02567640 | −0.19844679 | −0.19727080 |
0.5 | −0.030190188 | −0.02894723 | −0.24864875 | −0.24755700 |
0.6 | −0.030988255 | −0.02984030 | −0.29885599 | −0.29791001 |
0.7 | −0.028801308 | −0.02787389 | −0.34906512 | −0.34833031 |
0.8 | −0.023159274 | −0.02256617 | −0.39929111 | −0.39881842 |
0.9 | −0.013660270 | −0.01343546 | −0.44957637 | −0.44937485 |
1.0 | 0.00 | −2 × 10−11 | −0.50 | −0.50 |
X | V(X) | Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta }
(X) |
---|---|---|
DTM | DTM | |
−1.0 | 0.00 | 0.00 |
−0.9 | 0.00976 | 0.00037 |
−0.8 | 0.00259 | 0.00104 |
−0.7 | 0.00119 | 0.0019 |
−0.6 | 0.002027 | 0.00291 |
−0.5 | 0.000268 | 0.00409 |
−0.4 | 0.00409 | 0.00556 |
−0.3 | 0.01198 | 0.00771 |
−0.2 | 0.02727 | 0.01169 |
−0.1 | 0.06952 | 0.02307 |
0.0 | 0.99607 | 0.59975 |
0.1 | 0.125238 | 0.025281 |
0.2 | 0.070145 | 0.012436 |
0.3 | 0.053652 | 0.00818 |
0.4 | 0.045945 | 0.005926 |
0.5 | 0.041171 | 0.004391 |
0.6 | 0.037045 | 0.003165 |
0.7 | 0.032201 | 0.002105 |
0.8 | 0.02561 | 0.001184 |
0.9 | 0.016457 | 0.000448 |
1.0 | 0.00 | 0.00 |
|
Fig. 4. Effect of Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \delta } number on (a) velocity profile (V(X )) and (b) temperature profile (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta } (X )) for Cu–SA nanofluid when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pr=Ec=0.5} , Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \phi =0.05} . |
|
Fig. 5. Effect of nanoparticles volume fraction (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \varphi } ) on (a) velocity profile (V(X )) and (b) temperature distribution (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta } (X)) for SA–TiO2 nanofluid when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Ec=\delta =1} . |
|
Fig. 6. Effect of nanoparticles volume fraction on (a) velocity profile (V(X )) and (b) temperature profile (Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle \theta (X)} ) for Cu–SA nanofluid when Failed to parse (MathML with SVG or PNG fallback (recommended for modern browsers and accessibility tools): Invalid response ("Math extension cannot connect to Restbase.") from server "https://test.scipedia.com:8081/localhost/v1/":): {\textstyle Pr=Ec=\delta =0.5} . |
In present work, Differential Transformation Method (DTM) is applied for solving the coupled nonlinear differential equations in fluids mechanic to show the validity and simplicity of this method. For this aim, two cases of problems in nanofluid mechanics which have coupled nonlinear differential equations were selected and DTM was applied on them. For illustrating the accuracy of described method, a numerical method was also used to solve the problems. A very good agreement between these two methods was observed in both problems. Compared to previous works, DTM has a better agreement with numerical method compared to other analytical and perturbation methods. This DTM accuracy is due to independency to small parameter (p) compared to perturbation methods and direct solving of the equation without linearization or discretization.
D. Jing gratefully acknowledges the financial support from the National Natural Science Foundation of China (Nos. 51422604 and 21276206) and China Fundamental Research Funds for the Central Universities.
Published on 11/04/17
Licence: Other
Are you one of the authors of this document?