Homotopy perturbation method for special nonlinear partial differential equations

Homotopy perturbation method (HPM) is a semi-analytical technique for solving linear as well as nonlinear ordinary/partial differential equations. The method may also be used to solve a system of coupled linear and nonlinear differential equations. This chapter first illustrates the basic idea of the HPM. It uses the present method to solve a nonlinear ordinary differential equation, a nonlinear partial differential equation, and a coupled nonlinear system of partial differential equations. It is worth mentioning that the HPM can also be used for handling linear ordinary and partial differential equations. The test problems demonstrated in the chapter confirm that the HPM is an efficient method for solving linear/nonlinear ordinary, partial, and coupled differential equations. The series usually converges with an increase in the number of terms, but one may not expect the compact form solution always.

: Complex physical phenomena involving chemically reacting systems or the transport of heat or fluids are often modeled by coupled systems of time-dependent, nonlinear partial differential equations. The difficulties in understanding the stability of the differential equation systems and in designing efficient, accurate numerical methods for their solution are widely recognized and were the focus of this research. We have worked on four general aspects of the analysis and numerical approximation of systems of partial differential equations. These areas of research are: (1) modeling aspects and stability analysis for nonlinear time-dependent partial differential equations; (2) use and analysis of finite element or finite difference methods to discretize coupled systems of nonlinear differential equations; (3) development of adaptive or local grid refinement capabilities to resolve local phenomena in large-scale applications; and (4) development of data structures, preconditioners, and efficient solution algorithms for large-scale problems on new computer architectures. Emphasis has been placed upon multiphase or multicomponent, transport-dominated flow processes with dynamic local phenomena. The research also involved a mix of analysis, algorithm development, and large-scale computation using newer computer architectures.

A new treatment for homotopy perturbation method is introduced. The new treatment is called He‐Laplace method which is the coupling of the Laplace transform and the homotopy perturbation method using He’s polynomials. The nonlinear terms can be easily handled by the use of He’s polynomials. The method is implemented on linear and nonlinear partial differential equations. It is found that the proposed scheme provides the solution without any discretization or restrictive assumptions and avoids the round‐off errors.

We apply the homotopy perturbation method to obtain the solution of partial differential equations of fractional order. This method is powerful tool to find exact and approximate solution of many linear and nonlinear partial differential equations of fractional order. Convergence of the method is proved and the convergence analysis is reliable enough to estimate the maximum absolute truncated error of the series solution. The fractional derivatives are described in the Caputo sense. Some examples are presented to verify convergence hypothesis and simplicity of the method.

In this paper, the Laplace transform homotopy analysis method (LHAM) is employed to obtain approximate analytical solutions of the linear and nonlinear differential equations.This method is a combined form of the Laplace transform method and the homotopy analysis method.The proposed scheme finds the solutions without any discretization or restrictive assumptions and is free from round-off errors and therefore, reduces the numerical computations to a great extent.Some illustrative examples are presented and the numerical results show that the solutions of the LHAM are in good agreement with those obtained by exact solution.

An elegant and powerful technique is Homotopy Perturbation Method (HPM) to solve linear and nonlinear ordinary and partial differential equations. The method, which is a coupling of the traditional perturbation method and homotopy in topology, deforms continuously to a simple problem which can be solved easily. The method does not depend upon a small parameter in the equation. Using the initial conditions this method provides an analytical or exact solution. From the calculation and its graphical representation it is clear that how the solution of the original equation and its behavior depends on the initial conditions. Therefore there have been attempts to develop new techniques for obtaining analytical solutions which reasonably approximate the exact solutions. Many problems in natural and engineering sciences are modeled by nonlinear partial differential equations (NPDEs). The theory of nonlinear problem has recently undergone much study. Nonlinear phenomena have important applications in applied mathematics, physics, and issues related to engineering. In this paper we have applied this method to Burger’s equation and an example of highly nonlinear partial differential equation to get the most accurate solutions. The final results tell us that the proposed method is more efficient and easier to handle when is compared with the exact solutions or Adomian Decomposition Method (ADM).

A new higher order finite element method for elliptic partial differential equations on a stationary smooth surface $\Gamma$ is introduced and analyzed. We assume that $\Gamma$ is characterized as the zero level of a level set function $\phi$ and only a finite element approximation $\phi_h$ (of degree $k \geq 1$) of $\phi$ is known. For the discretization of the partial differential equation, finite elements (of degree $m\geq 1$) on a piecewise linear approximation of $\Gamma$ are used. The discretization is lifted to $\Gamma_h$, which denotes the zero level of $\phi_h$, using a quasi-orthogonal coordinate system that is constructed by applying a gradient recovery technique to $\phi_h$. A complete discretization error analysis is presented in which the error is split into a geometric error, a quadrature error, and a finite element approximation error. The main result is a $H^1(\Gamma)$-error bound of the form $c(h^m + h^{k+1})$. Results of numerical experiments illustrate the higher order convergence of t...

In the paper, two new analytic methods using the decomposition and linearization technique on nonlinear differential/integral equations are developed, namely, the decomposition–linearization–sequential method (DLSM) and the linearized homotopy perturbation method (LHPM). The DLSM is realized by an integrating factor and the integral of certain function obtained at the previous step for obtaining a sequential analytic solution of nonlinear differential equation, which provides quite accurate analytic solution. Some first- and second-order nonlinear differential equations display the fast convergence and accuracy of the DLSM. An analytic approximation for the Volterra differential–integral equation model of the population growth of a species is obtained by using the LHPM. In addition, the LHPM is also applied to the first-, second-, and third-order nonlinear ordinary differential equations. To reduce the cost of computation of He’s homotopy perturbation method and enhance the accuracy for solving cubically nonlinear jerk equations, the LHPM is implemented by invoking a linearization technique in advance is developed. A generalization of the LHPM to the nth-order nonlinear differential equation is involved, which can greatly simplify the work to find an analytic solution by solving a set of second-order linear differential equations. A remarkable feature of those new analytic methods is that just a few steps and lower-order approximations are sufficient for producing reasonably accurate analytic solutions. For all examples, the second-order analytic solution x2(t) is found to be a good approximation of the real solution. The accuracy of the obtained approximate solutions are identified by the fourth-order Runge–Kutta method. The major objection is to unify the analytic solution methods of different nonlinear differential equations by simply solving a set of first-order or second-order linear differential equations. It is clear that the new technique considerably saves computational costs and converges faster than other analytical solution techniques existing in the literature, including the Picard iteration method. Moreover, the accuracy of the obtained analytic solution is raised.

We present a control design method for nonlinear partial differential equations (PDEs) based on a combination of gain scheduling and backstepping theory for linear PDEs. A benchmark first-order hyperbolic system with an in-domain nonlinearity is considered first. For this system a nonlinear feedback law, based on gain scheduling, is derived explicitly, and a proof of local exponential stability, with an estimate of the region of attraction, is presented for the closed-loop system. Control designs (without proofs) are then presented for a string PDE and a shear beam PDE, both with Kelvin–Voigt (KV) damping and free-end nonlinearities of a potentially destabilizing kind. String and beam simulation results illustrate the merits of the gain scheduling approach over the linearization based design.

In this paper, we extend the invariant subspace method to coupled fractional higher order nonlinear partial differential equations. In particular, we consider a two coupled fractional partial differential equations of the form [Formula: see text] where [Formula: see text], [Formula: see text] are third-order nonlinear differential operators with quadratic nonlinearities and study the four-dimensional subspaces admitted by them. Using the obtained invariant subspaces we find explicit solutions to many two coupled fractional partial differential equations including the time-fractional equivalent of potential Korteweg–de Vries equation, reduced Super Korteweg–de Vries equation, supersymmetric extension of Korteweg–de Vries equation and generalized Drinfeld–Sokolov equation.

ABSTRACTIn this article, we have established the homotopy analysis method (HAM) for solving a few partial differential equations arising in engineering. This technique provides the solutions in rapid convergence series with computable terms for the problems with high degree of nonlinear terms appearing in the governing differential equations. The convergence analysis of the proposed method is also discussed. Finally, we have given some illustrative examples to demonstrate the validity and applicability of the proposed method.

In this paper, the weighted average-based differential quadrature method is presented for solving one-dimensional homogeneous first-order non-linear parabolic partial differential equation. First, the given solution domain is discretized by using uniform discretization grid point. Next, by using Taylor series expansion we obtain central finite difference discretization of the partial differential equation involving with temporal variable associated with weighted average of partial derivative concerning spatial variable. From this, we obtain the system of nonlinear ordinary differential equations and it is linearized by using the quasilinearization method. Then by using the polynomial-based differential quadrature method for approximating derivative involving with spatial variable at specified grid point, we obtain the system of linear equation. Then they obtained linear system equation is solved by using the LU matrix decomposition method. To validate the applicability of the proposed method, two model examples are considered and solved at each specific grid point on its solution domain. The stability and convergent analysis of the present method is worked by supported the theoretical and mathematical statements and the accuracy of the solution is obtained. The accuracy of the present method has been shown in the sense of root mean square error norm and maximum absolute error norm and the local behavior of the solution is captured exactly. Numerical versus exact solutions and behavior of maximum absolute error between them have been presented in terms of graphs and the corresponding root mean square error norm and maximum absolute error norm presented in tables. The present method approximates the exact solution very well and it is quite efficient and practically well suited for solving the non-linear parabolic equation. The numerical result presented in tables and graphs indicates that the approximate solution is in good agreement with the exact solution.

In this paper, the homotopy analysis method is applied to solve linear fractional problems. Based on this method, a scheme is developed to obtain approximation solution of fractional wave, Burgers, Korteweg-de Vries (KdV), KdV-Burgers, and Klein-Gordon equations with initial conditions, which are introduced by replacing some integer-order time derivatives by fractional derivatives. The fractional derivatives are described in the Caputo sense. So the homotopy analysis method for partial differential equations of integer order is directly extended to derive explicit and numerical solutions of the fractional partial differential equations. The solutions are calculated in the form of convergent series with easily computable components. The results of applying this procedure to the studied cases show the high accuracy and efficiency of the new technique.

Variable separation method for a nonlinear time fractional partial differential equation with forcing term

The magnetohydrodynamic (MHD) flow of a third grade fluid through a rectangular channel, considering the effect of aspect ratio, has been investigated. The flow considered is steady, laminar, incompressible and hydro-dynamically fully developed. The equation, describing the flow, is a highly non-linear partial differential equation (PDE) with remote possibility of having an exact solution and even numerical solution also is very difficult to obtain. A combination of the homotopy perturbation method (HPM) and integral method (IM) has been employed to solve the non-linear PDE which is scarce in open literature. The results of the present study are compared with the results obtained by the least square method (LSM) of the MHD third grade fluid flow through a rectangular channel, without the effect of aspect ratio and are found to be in close agreement. The results indicate that the flow field is significantly affected by the aspect ratio which should be considered for practical applications. In all the available literatures of the third grade fluid flow, the aspect ratio effect is neglected and this simplifying assumption reduces the highly complicated non-linear PDE to a non-linear ordinary differential equation (ODE). The novelty of the subject work lies in the inclusion of the effects of aspect ratio in the governing equation describing the flow of a third grade fluid through a channel and solving this by a combined analytical method (HPM and IM). Further, the effects of the Hartmann number and non-Newtonian third grade fluid parameter on the flow filed are discussed.