Lifeng Wang¹, Yunpeng Ma^1,2, Yongqiang Yang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.101, No.2, pp. 97-111, 2014, DOI:10.3970/cmes.2014.101.097

Abstract In this paper, a numerical method based on the Legendre polynomials is presented for a class of variable order fractional differential equation. We adopt the Coimbra variable order fractional operator, which can be viewed as a Caputo-type definition. Three different kinds of operational matrixes with Legendre polynomials are derived. A truncated the Legendre polynomials series together with the products of several dependent matrixes are utilized to reduce the variable order fractional differential equation to a system of algebraic equations. The solution of this system gives the approximation solution for the truncated limited n. An error analysis technique is also given.… More >

Abdul-Majid Wazwaz¹, R,olph Rach², Lazhar Bougoffa³, Jun-Sheng Duan^{4, 5}

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.6, pp. 507-529, 2014, DOI:10.3970/cmes.2014.100.507

Abstract In this paper, we construct the Lane–Emden–Fowler type equations of higher orders. We study the linear and the nonlinear Lane–Emden–Fowler type equations of the third and fourth orders, where other forms can be treated in a similar manner. We use the systematic Adomian decomposition method to handle these types of equations with specified initial conditions. We confirm that the Adomian decomposition method provides an efficient algorithm for exact and approximate analytic solutions of these equations. We corroborate this study by investigating several numerical examples that emphasize initial value problems. More >

Lei Lu^1,2,3, Jun-Sheng Duan², Long-Zhen Fan¹

CMES-Computer Modeling in Engineering & Sciences, Vol.100, No.6, pp. 463-475, 2014, DOI:10.3970/cmes.2014.100.463

Abstract In this paper, the modified decomposition method (MDM) for solving the nonlinear two-dimensional viscous flow equations is presented. This study investigates the problem of laminar, isothermal, incompressible and viscous flow in a rectangular domain bounded by two moving porous walls, which enable the fluid to enter or exit during successive expansions or contractions. We first transform the original two-dimensional viscous flow problem into an equivalent fourth-order boundary value problem (BVP), then solve the problem by the MDM. The figures and tables clearly show high accuracy of the method to solve two-dimensional viscous flow. More >

Mingxu Yi¹, Jun Huang¹, Lifeng Wang¹

CMES-Computer Modeling in Engineering & Sciences, Vol.96, No.5, pp. 361-377, 2013, DOI:10.3970/cmes.2013.096.361

Abstract In this paper, operational matrix method based upon the Bernstein polynomials is proposed to solve the variable order fractional integro-differential equations in the Caputo derivative sense. We derive the Bernstein polynomials operational matrix of fractional order integration and introduce the product operational matrix of Bernstein polynomials. A truncated the Bernstein polynomials series together with the polynomials operational matrix are utilized to reduce the variable order fractional integro-differential equations to a system of algebraic equations. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Some examples are included to demonstrate the validity and applicability of the… More >

Jun-Sheng Duan^1,2, Randolph Rach³, Abdul-Majid Wazwaz⁴

CMES-Computer Modeling in Engineering & Sciences, Vol.94, No.1, pp. 77-118, 2013, DOI:10.3970/cmes.2013.094.077

Abstract In this paper, we propose a new modification of the Adomian decomposition method for solution of higher-order nonlinear initial value problems with variable system coefficients and solutions of systems of coupled nonlinear initial value problems. We consider various algorithms for the Adomian decomposition series and the series of Adomian polynomials to calculate the solutions of canonical first- and second-order nonlinear initial value problems in order to derive a systematic algorithm for the general case of higher-order nonlinear initial value problems and systems of coupled higher-order nonlinear initial value problems. Our new modified recursion scheme is designed to decelerate the Adomian… More >

Shou-Zhong Fu¹, Zhong Wang¹, Jun-Sheng Duan^1,2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.92, No.4, pp. 369-385, 2013, DOI:10.3970/cmes.2013.092.369

Abstract Quadratic integral equations are a class of nonlinear integral equations having many important uses in engineering and sciences. In this work we display an efficient application of the Adomian decomposition method to the quadratic integral equations of Volterra type. The analytical approximate solution obtained can be directly inserted into the original equation to verify the accuracy and estimate the error with a computing software. Four numerical examples demonstrate the efficiency of the method. More >

Yiming Chen, Mingxu Yi, Chen Chen, Chunxiao Yu

CMES-Computer Modeling in Engineering & Sciences, Vol.83, No.6, pp. 639-654, 2012, DOI:10.3970/cmes.2012.083.639

Abstract In this paper, Bernstein polynomials method is proposed for the numerical solution of a class of space-time fractional convection-diffusion equation with variable coefficients. This method combines the definition of fractional derivatives with some properties of Bernstein polynomials and are dispersed the coefficients efficaciously. The main characteristic behind this method is that the original problem is translated into a Sylvester equation. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples show that the method is effective. More >

Jun-Sheng Duan¹, Ai-Ping Guo²

CMES-Computer Modeling in Engineering & Sciences, Vol.60, No.2, pp. 139-150, 2010, DOI:10.3970/cmes.2010.060.139

Abstract Adomian polynomials are constituted of reduced polynomials and derivatives of nonlinear operator. The reduced polynomials are independent of the form of the nonlinear operator. A recursive algorithm of the reduced polynomials is discovered and its symbolic implementation by the software Mathematica is given. As a result, a new and convenient algorithm for the Adomian polynomials is obtained. More >

Chia-Cheng Tsai¹, Chein-Shan Liu², Wei-Chung Yeih³

CMES-Computer Modeling in Engineering & Sciences, Vol.56, No.2, pp. 131-152, 2010, DOI:10.3970/cmes.2010.056.131

Abstract The fictitious time integration method (FTIM) previously developed by Liu and Atluri (2008a) is combined with the method of fundamental solutions and the Chebyshev polynomials to solve Poisson-type nonlinear PDEs. The method of fundamental solutions with Chebyshev polynomials (MFS-CP) is an exponentially-convergent meshless numerical method which is able to solving nonhomogeneous partial differential equations if the fundamental solution and the analytical particular solutions of the considered operator are known. In this study, the MFS-CP is extended to solve Poisson-type nonlinear PDEs by using the FTIM. In the solution procedure, the FTIM is introduced to convert a Poisson-type nonlinear PDE into… More >

M. A. Kelmanson¹ and M. C. Tenwick¹

CMES-Computer Modeling in Engineering & Sciences, Vol.55, No.2, pp. 191-210, 2010, DOI:10.3970/cmes.2010.055.191

Abstract A method is presented for improving the accuracy of the widely used Gauss-Legendre Nyström method for determining approximate solutions of Fredholm integral equations of the second kind on finite intervals. The authors' recent continuous-kernel approach is generalised in order to accommodate kernels that are either singular or of limited continuous differentiability at a finite number of points within the interval of integration. This is achieved by developing a Gauss-Jacobi Nyström method that moreover includes a mean-value estimate of the truncation error of the Hermite interpolation on which the quadrature rule is based, making it particularly accurate at low orders. A… More >