Homotopy perturbation-reproducing kernel method for nonlinear systems of second order boundary value problems

In this paper, based on homotopy perturbation method (HPM) and reproducing kernel method (RKM), a new method is presented for solving nonlinear systems of second order boundary value problems (BVPs). HPM is based on the use of traditional perturbation method and homotopy technique. The HPM can reduce a nonlinear problem to a sequence of linear problems and generate a rapid convergent series solution in most cases. RKM is also an analytical technique, which can solve powerfully linear BVPs. Homotopy perturbation-reproducing kernel method (HP-RKM) combines advantages of these two methods and therefore can be used to solve efficiently systems of nonlinear BVPs. Three numerical examples are presented to illustrate the strength of the method.     

New method based on the HPM and RKHSM for solving forced Duffing equations with integral boundary conditions

This paper investigates the forced Duffing equation with integral boundary conditions. Its approximate solution is developed by combining the homotopy perturbation method (HPM) and the reproducing kernel Hilbert space method (RKHSM). HPM is based on the use of the traditional perturbation method and the homotopy technique. The HPM can reduce nonlinear problems to some linear problems and generate a rapid convergent series solution in most cases. RKHSM is also an analytical technique, which can solve powerfully linear boundary value problems. Therefore, the forced Duffing equation with integral boundary conditions can be solved using advantages of these two methods. Two numerical examples are presented to illustrate the strength of the method.     

Computing curve intersection by homotopy methods

Intersection problems are fundamental in computational geometry, geometric modeling and design and manufacturing applications, and can be reduced to solving polynomial systems. This paper introduces two homotopy methods, i.e. polyhedral homotopy method and linear homotopy method, to compute the intersections of two plane rational parametric curves. Extensive numerical examples show that computing curve intersection by homotopy methods has better accuracy, efficiency and robustness than by the Ehrlich-Aberth iteration method. Finally, some other applications of homotopy methods are also presented.     

Computing curve intersection by homotopy methods

Intersection problems are fundamental in computational geometry, geometric modeling and design and manufacturing applications, and can be reduced to solving polynomial systems. This paper introduces two homotopy methods, i.e. polyhedral homotopy method and linear homotopy method, to compute the intersections of two plane rational parametric curves. Extensive numerical examples show that computing curve intersection by homotopy methods has better accuracy, efficiency and robustness than by the Ehrlich–Aberth iteration method. Finally, some other applications of homotopy methods are also presented.     

A linearization‐based approach of homotopy analysis method for non‐linear time‐fractional parabolic PDEs

In this paper, a novel approach, namely, the linearization‐based approach of homotopy analysis method, to analytically treat non‐linear time‐fractional PDEs is proposed. The presented approach suggests a new optimized structure of the homotopy series solution based on a linear approximation of the non‐linear problem. A comparative study between the proposed approach and standard homotopy analysis approach is illustrated by solving two examples involving non‐linear time‐fractional parabolic PDEs. The performed numerical simulations demonstrate that the linearization‐based approach reduces the computational complexity and improves the performance of the homotopy analysis method.     

求解混合三角多项式方程组的直接GBQ方法

许多科学与工程领域,我们经常需要求混合三角多项式方程组的全部解.一般来说,混合三角多项式方程组可以通过变量替换及增加二次多项式转化为多项式方程组,进而利用数值方法进行求解,但这种转化会增大问题的规模从而增加计算量.在本文中,我们不将问题转化,考虑利用直接同伦方法求解,并给出基于GBQ方法构造的初始方程组及同伦定理的证明.数值实验结果表明我们构造的直接同伦方法较已有的直接同伦方法更加有效.     

A PL homotopy continuation method with the use of an odd map for the artificial level

This note presents a new piecewise linear homotopy continuation method for solving a system of nonlinear equations. The important feature of the method is the use of an odd map for the artificial level of the homotopy. Some sufficient conditions for the global convergence of the method are given. They are different from the known conditions for the global convergence of the existing homotopy continuation methods. Specifically, they cover all the systems of nondegenerate linear equations. Partially supported by the Sakkokai Foundation and by the Grand-in-Aid for Scientific Research of the Ministry of Education and Culture No. 58750265. Partially supported by the Grant-in-Aid for Scientific Research of the Ministry of Education and Culture No. 56460103.     

A spline smoothing homotopy method for nonconvex nonlinear programming

Homotopy methods are globally convergent under weak conditions and robust; however, the efficiency of a homotopy method is closely related with the construction of the homotopy map and the path tracing algorithm. Different homotopies may behave very different in performance even though they are all theoretically convergent. In this paper, a spline smoothing homotopy method for nonconvex nonlinear programming is developed using cubic spline to smooth the max function of the constraints of nonlinear programming. Some properties of spline smoothing function are discussed and the global convergence of spline smoothing homotopy under the weak normal cone condition is proven. The spline smoothing technique uses a smooth constraint instead of m constraints and acts also as an active set technique. So the spline smoothing homotopy method is more efficient than previous homotopy methods like combined homotopy interior point method, aggregate constraint homotopy method and other probability one homotopy methods. Numerical tests with the comparisons to some other methods show that the new method is very efficient for nonlinear programming with large number of complicated constraints.     

A seminumeric approach for solution of the Eikonal partial differential equation and its applications

In this work, a partial differential equation, which has several important applications, is investigated, and some techniques based on semianalytic (or quasi‐numerical) approaches are developed to find its solution. In this article, the homotopy perturbation method (HPM), Adomian decomposition method, and the modified homotopy perturbation method are proposed to solve the Eikonal equation. HPM yields solution in convergent series form with easily computable terms, and in some case, yields exact solutions in one iteration. In other hand, in Adomian decomposition method, the approximate solution is considered as an infinite series usually converges to the accurate solution. Moreover, these methods do not require any discretization, linearization, or small perturbation, and therefore reduce the numerical computation a lot. Several test problems are given and results are compared with the variational iteration method. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2010     

Comparison of approximate symmetry and approximate homotopy symmetry to the Cahn–Hilliard equation

Complete infinite order approximate symmetry and approximate homotopy symmetry classifications of the Cahn–Hilliard equation are performed and the reductions are constructed by an optimal system of one-dimensional subalgebras. Zero order similarity reduced equations are nonlinear ordinary differential equations while higher order similarity solutions can be obtained by solving linear variable coefficient ordinary differential equations. The relationship between two methods for different order are studied and the results show that the approximate homotopy symmetry method is more effective to control the convergence of series solutions than the approximate symmetry one.     

The homotopy analysis method for solving the time-fractional Fornberg–Whitham equation and comparison with Adomian’s decomposition method

In this paper, we applied relatively new analytical techniques, the homotopy analysis method (HAM) and the Adomian’s decomposition method (ADM) for solving time-fractional Fornberg–Whitham equation. The homotopy analysis method contains the auxiliary parameter, which provides us with a simple way to adjust and control the convergence region of solution series. The fractional derivatives are described in the Caputo sense. A comparison is made the between HAM and ADM results. The present methods performs extremely well in terms of efficiency and simplicity. Numerical results for different particular cases of the problem are presented.     

On homotopy analysis method applied to linear optimal control problems

In this paper, the homotopy analysis method (HAM) is employed to solve the linear optimal control problems (OCPs), which have a quadratic performance index. The study examines the application of the homotopy analysis method in obtaining the solution of equations that have previously been obtained using the Pontryagin’s maximum principle (PMP). The HAM approach is also applied in obtaining the solution of the matrix Riccati equation. Numerical results are presented for several test examples involving scalar and 2nd-order systems to demonstrate the applicability and efficiency of the method.     

Method for solving nonlinear initial value problems by combining homotopy perturbation and reproducing kernel Hilbert space methods

The nonlinear singular initial value problems including generalized Lane–Emden-type equations are investigated by combining homotopy perturbation method (HPM) and reproducing kernel Hilbert space method (RKHSM). He’s HPM is based on the use of traditional perturbation method and homotopy technique and can reduce a nonlinear problem to some linear problems and generate a rapid convergent series solution in most cases. RKHSM is also an analytical technique, which can overcome the difficulty at the singular point of non-homogeneous, linear singular initial value problems; especially when the singularity appears on the right-hand side of this type of equations, so it can solve powerfully linear singular initial value problems. Therefore, using advantages of these two methods, more general nonlinear singular initial value problems can be solved powerfully. Some numerical examples are presented to illustrate the strength of the method.     

Application of He’s Variational Iteration Method for Solving Nonlinear BBMB Equations and Free Vibration of Systems

A new analytical method called He’s variational iteration method (VIM) is introduced to be applied to solve nonlinear Benjamin-Bona-Mahony-Burgers (BBMB) equations and free vibration of a nonlinear system having combined linear and nonlinear springs in series in this article. In this method, general Lagrange multipliers are introduced to construct correction functionals for the problems. The multipliers can be identified optimally via the variational theory. The results are compared with the results of the homotopy analysis method and also with the exact solution. He’s Variational iteration method in this problem functions so better than the homotopy analysis method and exact solutions one of them in per section.     

An effective modification of the homotopy perturbation method for stiff systems of ordinary differential equations

The present work deals with employing a new form of the homotopy perturbation method (NHPM) for solving stiff systems of linear and nonlinear ordinary differential equations (ODEs). In this scheme, the solution is considered as an infinite series that converges rapidly to the exact solution. Two problems are chosen as illustrative examples to show the effectiveness of the present method. In obtaining the exact solution for each case, the capability and the simplicity of the proposed technique is clarified.     

函数因子法——非线性方程组求解中处理奇异问题的一种新方法

§1.引言 非线性方程组 F(x)=0,F:D?R~n→R~n (1.1)嵌入参数t,构成同伦H:[0,T]×D?R~(n 1)→R~n,使得 H(0,x~0)=0,H(T,x)=F(x),(1.2)这里T可以是有限的或 ∞,当T为 ∞时以极限过程代替求值.若 H(t.x)=0(1.3)存在连续解x(t):[0,T]→D,则非线性方程组(1.1)的解x~*=x(T).若(1.3)的解     

Comparison of the Adomian decomposition method with homotopy perturbation method for the solutions of seventh order boundary value problems

The aim of this paper is to compare the Adomian decomposition method and the homotopy perturbation method for solving the linear and nonlinear seventh order boundary value problems. The approximate solutions of the problems obtained with a small amount of computation in both methods. Two numerical examples have been considered to illustrate the accuracy and implementation of the methods.     

On the numerical solution of singular two‐point boundary value problems: A domain decomposition homotopy perturbation approach

This paper reports a modified homotopy perturbation algorithm, called the domain decomposition homotopy perturbation method (DDHPM), for solving two‐point singular boundary value problems arising in science and engineering. The essence of the approach is to split the domain of the problem into a number of nonoverlapping subdomains. In each subdomain, a method based on a combination of HPM and integral equation formalism is implemented. The boundary condition at the right endpoint of each inner subdomain is established before deriving an iterative scheme for the components of the solution series. The accuracy and efficiency of the DDHPM are demonstrated by 4 examples (2 nonlinear and 2 linear). In comparison with the traditional HPM, the proposed domain decomposition HPM is highly accurate.     

Laplace transform for the solution of higher order deformation equations arising in the homotopy analysis method

A new analytic approach for solving nonlinear ordinary differential equations with initial conditions is proposed. First, the homotopy analysis method is used to transform a nonlinear differential equation into a system of linear differential equations; then, the Laplace transform method is applied to solve the resulting linear initial value problems; finally, the solutions to the linear initial value problems are employed to form a convergent series solution to the given problem. The main advantage of the new approach is that it provides an effective way to solve the higher order deformation equations arising in the homotopy analysis method.