# Application of the Homotopy Perturbation Method for Solving the Foam Drainage Equation

Authors: Abdolhosein  Fereidoon; Hessameddin Yaghoobi; Mohammadreza Davoudabadi
Journal: International Journal of Differential Equations
 (2011)

---

## Abstract

We use homotopy perturbation method (HPM) to handle the foam drainage equation. Foaming occurs in many distillation and absorption processes. The drainage of liquid foams involves the interplay of gravity, surface tension, and viscous forces. The concept of He's homotopy perturbation method is introduced briefly for applying this method for problem solving. The results of HPM as an analytical solution are then compared with those derived from Adomian's decomposition method (ADM) and the variational iteration method (VIM). The results reveal that the HPM is very effective and convenient in predicting the solution of such problems, and it is predicted that HPM can find a wide application in new engineering problems.


---

## Body

## 1. Introduction

Most scientific problems and physical phenomena occur nonlinearly. Except in a limited number of these problems, finding the exact analytical solutions of such problems are rather difficult. Therefore, there have been attempts to develop new techniques for obtaining analytical solutions which reasonably approximate the exact solutions [1]. In recent years, several such techniques have drawn special attention, such as Hirota’s bilinear method [2], the homogeneous balance method [3, 4], inverse scattering method [5], Adomian’s decomposition method (ADM) [6, 7], the variational iteration method (VIM) [8, 9], and homotopy analysis method (HAM) [10] as well as homotopy perturbation method (HPM).The homotopy perturbation method was established by He [11]. The method has been used by many authors to handle a wide variety of scientific and engineering applications to solve various functional equations. In this method, the solution is considered as the sum of an infinite series, which converges rapidly to accurate solutions. Recently, considerable research work has been conducted in applying this method to the linear and nonlinear equations [12–25].Foams are of great importance in many technological processes and applications, and their properties are subject of intensive studies from both practical and scientific points of view [26]. Liquid foam is an example of soft matter (or complex fluid) with a very well-defined structure that first clearly described by Joseph plateau in the 19th century. Weaire et al. [27] showed in their work simple answers to many such questions exist, but no going experiments continue to challenge our understanding. Foams and emulsions are wellknown to scientists and the general public alike because of their everyday occurrence [28, 29]. Foams are common in foods and personal care products such as creams and lotions, and foams often occur, even when not desired, during cleaning (clothes, dishes, scrubbing) and dispensing processes [30]. They have important applications in the food and chemical industries, firefighting, mineral processing, and structural material science [31]. Less obviously, they appear in acoustic cladding, lightweight mechanical components, and impact absorbing parts on cars, heat exchangers, and textured wallpapers (incorporated as foaming inks) and even have an analogy in cosmology. The packing of bubbles or cells can form both random and symmetrical arrays, such as sea foam and bees’ honeycomb. History connects foams with a number of eminent scientists, and foams continue to excite imaginations [32]. There are now many applications of polymeric foams [33] and more recently metallic foams, which are foams made of metals such as aluminum [34]. Some commonly mentioned applications include the use of foams for reducing the impact of explosions and for cleaning up oil spills. In addition, industrial applications of polymeric foams and porous metals include their use for structural purposes and as heat exchange media analogous to common “finned” structures [35]. Polymeric foams are used in cushions and packing and structural materials [36]. Glass, ceramic, and metal foams [37] can also be made and find an increasing number of new applications. In addition, mineral processing utilizes foam to separate valuable products by flotation. Finally, foams enter geophysical studies of the mechanics of volcanic eruptions [30]. Recent research in foams and emulsions has centered on three topics which are often treated separately but are, in fact, interdependent: drainage, coarsening, and rheology; see Figure 1. We focus here on a quantitative description of the coupling of drainage and coarsening. Foam drainage is the flow of liquid through channels (plateau borders) and nodes (intersections of four channels) between the bubbles, driven by gravity and capillarity [38–40]. During foam production, the material is in the liquid state, and fluid can rearrange while the bubble structure stays relatively unchanged. The flow of liquid relative to the bubbles is called drainage. Generally, drainage is driven by gravity and/or capillary (surface tension) forces and is resisted by viscous forces [30]. Because of their limited time stability and despite the numerous studies reported in the literature, many of their properties are still not well understood, in particular the drainage of the liquid in between the bubbles under the influence of gravity [41, 42]. Drainage plays an important role in foam stability. Indeed, when foam dries, its structure becomes more fragile; the liquid films between adjacent bubbles being thinner, then can break, leading to foam collapse. In the case of aqueous foams, surfactant is added into water, and it adsorbs at the surface of the films, protecting them against rupture [43]. Most of the basic rules that explain the stability of liquid gas foams were introduced over 100 years ago by the Belgian Joseph Plateau who was blind before he completed his important book on the subject. This modern-day book by Weaire and Hutzler provides valuable summaries of plateaus work on the laws of equilibrium of soap films, and it is especially useful since the original 1873 French text does not appear to be in a fully translated English version. Weaire and Hutzler note that Sir W. Thompson (Lord Kelvin) was simulated by Plateau’s book to examine the optimum packing of free space by regular geometrical cells. His solution to the problem remained the best until quite recently. Why does this area of theoretical research, still active today, have connections with the apparently frivolous theme of bubbles? It is because the packing of free space involves the minimization of the surface energy of the structure (i.e., least amount of boundary material). Thus, one might ask why such an often-observed medium as a foam has not provided the optimum solution to this problem much earlier; perhaps, this shows that observation is often biased towards what one expects to see, rather than to the unexpected. Also, in nature, there are packing problems, such as the bees’ honeycomb. Its shaped ends provide a nice example of Plateau’s rules in a natural environment [32]. Recent theoretical studies by Verbist and Weaire describe the main features of both free drainage [44, 45], where liquid drains out of a foam due to gravity, and forced drainage [46], where liquid is introduced to the top of a column of foam. In the latter case, a solitary wave of constant velocity is generated when liquid is added at a constant rate [47]. Forced foam drainage may well be the best prototype for certain general phenomena described by nonlinear differential equations, particularly the type of solitary wave which is most familiar in tidal bores. The model developed by Verbist and Weaire idealizes the network of Plateau borders, through which the majority of liquid is assumed to drain, as a set of N identical pipes of cross section A, which is a function of position and time [48](1.1)∂A∂t+∂∂x(A2-A2∂A∂x)=0,
where x and t are scaled position and time coordinates, respectively. In the case of forced drainage, the solution can be expressed as [48](1.2)A(x,t)={ctanh2⁡(c(x-ct))x≤ct,0x≤ct,
where c is the velocity of the wave front [46].Figure 1
Schematic illustration of the interdependence of drainage, coarsening, and rheology of foams [31].The pursuit of analytical solutions for foam drainage equation is of intrinsic scientific interest. To the best of the authors’ knowledge, there is no paper that has solved the nonlinear foam drainage equation by HPM. In this paper, the basic idea of HPM is described, and then, it is applied to study the following nonlinear foam drainage equation [49]. Finally, the results of HPM as an analytical solution are then compared with those derived from Adomian’s decomposition method [49] and the variational iteration method [50].


## 2. Basic Idea of Homotopy Perturbation Method

To explain this method, let us consider the following function:(2.1)A(u)-f(r)=0,r∈Ω,
with the boundary conditions of(2.2)B(u,∂u∂n)=0,r∈Γ,
where A is a general differential operator f(r) is a known analytic function, B is a boundary operator, and Γ is the boundary of the domain Ω. The operator A can be generally divided into two operators, L and N, where L is a linear and N a nonlinear operator. Equation (2.1) can be, therefore, written as follows:(2.3)L(u)+N(u)-f(r)=0.
Using the homotopy technique, we constructed a homotopy v(r,p):Ω×[0,1]→R, which satisfies(2.4)H(v,p)=(1-p)[L(v)-L(u0)]+p[A(v)-f(r)]=0,
or(2.5)H(v,p)=L(v)-L(u0)+pL(u0)+p[N(v)-f(r)]=0,
where p∈[0,1] is called homotopy parameter and u0 is an initial approximation for the solution of (2.1), which satisfies the boundary conditions. Obviously, from (2.4) or (2.5), we will have(2.6)H(v,0)=L(v)-L(u0)=0,H(v,1)=A(v)-f(r)=0.
We can assume that the solution of (2.4) or (2.5) can be expressed as a series in p as follows:(2.7)v=v0+pv1+p2v2+⋯.
Setting p=1 results in the approximate solution of  (2.1)(2.8)u=lim⁡p→1⁡v=v0+v1+v2+⋯.


## 3. Implementation of HPM

In this section, we obtain an analytical solution of (1.1). We first use the transformation(3.1)A(x,t)=u2(x,t),
to convert (1.1) to(3.2)ut+2u2ux-(ux)2-12uxxu=0,
with initial condition(3.3)u(x,0)=-ctanh⁡(cx).
To solve (3.2) with the initial condition (3.3), according to the homotopy perturbation method, we construct the following He’s polynomials corresponding to (2.5):(3.4)L(v)=vt,N(v)=2v2vx-(vx)2-12vxxv,H(v,p)=vt-u0t+pu0t+p[2v2vx-(vx)2-12vxxv].
Substituting v=v0+pv1+p2v2+⋯ into (3.2) and rearranging the resultant equation based on powers of p-terms, one has(3.5)p0:v0t-u0t=0,p1:v1t+u0t+[2v02v0x-(v0x)2-12v0xxv0]=0,p2:v2t+4v0v1v0x+2v02v1x-2v0xv1x-12v0xxv1-12v1xxv0=0,
with the following conditions:(3.6)v0(x,0)=-ctanh(cx),(3.7)vi(x,0)=0fori=1,2,….
With the effective initial approximation for v0 from the condition (3.6), the solutions of (3.5) are obtained as follows:(3.8)v0(x,t)=-ctanh⁡(cx),v1(x,t)=2c2tcosh⁡⁡(2cx)2,v2(x,t)=c7/2t2sinh⁡⁡(cx)cosh⁡⁡(cx)3.
In the same manner, the rest of components were obtained using the Maple package. According to the HPM, we can conclude that(3.9)u(x,t)=v0(x,t)+v1(x,t)+v2(x,t).
Therefore, we will have(3.10)u(x,t)=-ctanh(cx)+2c2tcosh⁡⁡(2cx)2+c7/2t2sinh⁡⁡(cx)cosh⁡⁡(cx)3.
By substituting (3.10) into (3.1), we can find the solution of (1.1).


## 4. Results and Discussion

In this section, we present the results with HPM to show the efficiency of the method, described in the previous section for solving (1.1). By Figures 2, 3, and 4, we may simply compare the HPM solution and exact solution of (1.1) for c=1, respectively. It is easy to verify the accuracy of the obtained results if we graphically compare HPM solutions with the exact ones. The absolute error of HPM is drawn in Figure 5 at t=0.01. It can be seen from this figure that the absolute error is very small and the mentioned method is very accurate. The effect of velocity of the wave, c, is demonstrated in Figure 6.The solution of (1.1) at c=1 (a) HPM (b) Exact.
(a)(b)The solution of (1.1) at c=2 (a) HPM (b) Exact.
(a)(b)The solution of (1.1) at c=3 (a) HPM (b) Exact.
(a)(b)The absolute error of HPM att=0.01 (a) c=1 (b) c=2.
(a)(b)Figure 6
The effect ofc at t=0.01.It is important to see the difference between the results obtained by ADM [49] and VIM [50] and the results which are obtained by the analytic solution using by the HPM. These can be seen by comparing Tables 1, 2, and 3. The results which are obtained by the HPM and the exact solutions are quite similar. It is shown the accuracy of the HPM. The solution of the HPM is more accurate than the ADM. While using HPM, the difficulty of calculating Adomian's polynomials does not occur. Another benefit of using the HPM is that the time consumption of the numerical calculation is less than the time consumption of ADM.Table 1
Comparison between absolute errors of ADM, VIM, and HPM fort=0.001, c=3.
xAexact-AADM [49]Aexact-AVIM [50]Aexact-AHPM−104.44089E-1600−82.24265E-1302.4E-18−62.29754E-1002.1059E-15−42.34498E-702.14918E-12−20.00023665602.17229E-9−10.0052383405.10302E-805.2479E-801.45797E-9Table 2
Comparison between absolute errors of ADM, VIM, and HPM fort=0.01, c=3.
xAexact-AADM [49]Aexact-AVIM [50]Aexact-AHPM−101.77636E-1502.00000E-18−81.86962E-1202.01600E-15−61.9087E-902.05744E-12−41.94811E-66E-082.09993E-9−20.001972966.0881E-050.000002123−10.0485679—0.00005043300.000515928.08544E-30.000014557Table 3
Comparison between absolute errors of ADM, VIM, and HPM fort=0.1, c=3.
xAexact-AADM [49]Aexact-AVIM [50]Aexact-AHPM−101.42109E-1401.59300E-15−81.40941E-1101.62593E-12−61.43842E-85E-091.65952E-9−40.00001467274.527E-061.69379E-6−20.00642184.59239E-030.001716375−10.094494—0.04376768903.713676.8374E-010.126259125


## 5. Conclusion

In this work, He’s homotopy perturbation method has been successfully utilized to derive approximate explicit analytical solution for the nonlinear foam drainage equation. The results show that this perturbation scheme provides excellent approximations to the solution of this nonlinear equation with high accuracy and avoids linearization and physically unrealistic assumptions. This new method is extremely simple, easy to apply, needs less computation, and accelerates the convergence to the solutions. The results obtained here are compared with results of exact solution. The current work illustrates that the HPM is indeed a powerful analytical technique for most types of nonlinear problems and several such problems in scientific studies and engineering may be solved by this method.



---

*Source: 864023-2011-08-03.xml*