An Introduction to Numerical Methods for the Solutions of Partial Differential Equations

Partial differential equations arise in formulations of problems involving functions of several variables such as the propagation of sound or heat, electrostatics, electrodynamics, fluid flow, and elasticity, etc. The present paper deals with a general introduction and classification of partial differential equations and the numerical methods available in the literature for the solution of partial differential equations.

Comparative Studies between Picard’s and Taylor’s Methods of Numerical Solutions of First Ordinary Order Differential Equations Arising from Real-Life Problems

To solve the first-order differential equation derived from the problem of a free-falling object and the problem arising from Newton’s law of cooling, the study compares the numerical solutions obtained from Picard’s and Taylor’s series methods. We have carried out a descriptive analysis using the MATLAB software. Picard’s and Taylor’s techniques for deriving numerical solutions are both strong mathematical instruments that behave similarly. All first-order differential equations in standard form that have a constant function on the right-hand side share this similarity. As a result, we can conclude that Taylor’s approach is simpler to use, more effective, and more accurate. We will contrast Rung Kutta and Taylor’s methods in more detail in the following section.

THE NUMERICAL STABILITY OF THE BLOCK θ-METHODS FOR DELAY DIFFERENTIAL EQUATIONS

This paper focuses on the numerical stability of the block θ methods adapted to differential equations with a delay argument. For the block θ methods, an interpolation procedure is introduced which leads to the numerical processes that satisfy an important asymptotic stability condition related to the class of test problems y′(t)=ay(t)+by(t-τ) with a,b∈C, Re(a)<-|b| and τ>0. We prove that the block θ method is GP stable if and only if the method is A stable for ordinary differential equations. Furthermore, it is proved that the P and GP stability are equivalent for the block θ method.

Crank-Nicolson ADI Galerkin Finite Element Methods for Two Classes of Riesz Space Fractional Partial Differential Equations

In this paper,two classes of Riesz space fractional partial differential equations including space-fractional and space-time-fractional ones are considered.These two models can be regarded as the generalization of the classical wave equation in two space dimensions.Combining with the Crank-Nicolson method in temporal direction,efficient alternating direction implicit Galerkin finite element methods for solving these two fractional models are developed,respectively.The corresponding stability and convergence analysis of the numerical methods are discussed.Numerical results are provided to verify the theoretical analysis.

Modified Laguerre spectral and pseudospectral methods for nonlinear partial differential equations in multiple dimensions

The Laguerre spectral and pseudospectral methods are investigated for multidimensional nonlinear partial differential equations. Some results on the modified Laguerre orthogonal approximation and interpolation are established, which play important roles in the related numerical methods for unbounded domains. As an example, the modified Laguerre spectral and pseudospectral methods are proposed for two-dimensional Logistic equation. The stability and convergence of the suggested schemes are proved. Numerical results demonstrate the high accuracy of these approaches.

THE STABILITY OF LINEAR MULTISTEP METHODS FOR SYSTEMS OF DELAY DIFFERENTIAL EQUATIONS

This paper deals with the numerical solution of initial value problems for systems of differential equations with a delay argument. The numerical stability of a linear multistep method is investigated by analysing the solution of the lest equation y’(t)=Ay(t) + By(1-t),where A,B denote constant complex N×N-matrices,and t】0.We investigate carefully the characterization of the stability region.

Fourth-Order Splitting Methods for Time-Dependant Differential Equations

This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.

Method of Lines for Third Order Partial Differential Equations

The method of lines is applied to the boundary-value problem for third order partial differential equation. Explicit expression and order of convergence for the approximate solution are obtained.

A lattice Boltzmann model with an amending function for simulating nonlinear partial differential equations

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations （NPDEs） in the form ut ＋αuux ＋βu^nuz ＋γuxx ＋δuzxx ＋ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.

Symbolic computation and exact traveling solutions for nonlinear partial differential equations

In this paper, with the aid of the symbolic computation, a further extended tanh function method was presented. Based on the new general ansatz, many nonlinear partial differential equation（s）（NPDE（s）） can he solved. Especially, as applications, a compound KdV-mKdV equation and the Broer-Kaup equations are considered successfully, and many solutions including periodic solutions, triangle solutions, and rational solutions are obtained. The method can also be applied to other NPDEs.

New Implementation of Legendre Polynomials for Solving Partial Differential Equations

In this paper we present a proposal using Legendre polynomials approximation for the solution of the second order linear partial differential equations. Our approach consists of reducing the problem to a set of linear equations by expanding the approximate solution in terms of shifted Legendre polynomials with unknown coefficients. The performance of presented method has been compared with other methods, namely Sinc-Galerkin, quadratic spline collocation and LiuLin method. Numerical examples show better accuracy of the proposed method. Moreover, the computation cost decreases at least by a factor of 6 in this method.

Mixture of a New Integral Transform and Homotopy Perturbation Method for Solving Nonlinear Partial Differential Equations

In this paper, we present a new method, a mixture of homotopy perturbation method and a new integral transform to solve some nonlinear partial differential equations. The proposed method introduces also He’s polynomials [1]. The analytical results of examples are calculated in terms of convergent series with easily computed components [2].

Double Elzaki Transform Decomposition Method for Solving Non-Linear Partial Differential Equations

In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Transform and Adomian Decomposition Method. This technique is hereafter provided and supported with necessary illustrations, together with some attached examples. The results reveal that the new method is very efficient, simple and can be applied to other non-linear problems.

The θ－Methods in Numerical Solution of Systems of Differential Equations with Two Delay Terms

This paper deals with the numerical solution of initial value problems for systems of differential equations with two delay terms. We investigate the stability of adaptations of the θ-methods in the numerical solution of test equations u'(t) = a 11 u(t) + a12v(t) + b11 u(t - τ1) + b12v(t-τ2,v'(t) = a21 u(t) + a22 v(t) + b21 u(t -τ1,) + b22 v(t -τ2), t>0,with initial conditionsu(t)=u0(t),v(t) =v0(t), t≤0.where aij, bij∈C, τj >0, i,j = 1,2,, and u0(t), v0(t)are continuous and complex valued. Sufficient conditions for the asymptotic stability of test equation are derived. Furthermore, with respect to an appropriate definition of stability for the numerical method, it is proved that the linear θ-method is stable if and only if 1/2≤θ≤1 and the one-leg θ-method is stable if and only if θ= 1.

A SOLVING METHOD FOR A SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS WITH AN APPLICATION TO THE BENDING PROBLEM OF A THICK PLATE

A theorem of solving a system of linear non-homogeneous differential equations through integrating and adding its basic solutions is put forward and proved, the mathematical role and physical nature of the theorem is interpreted briefly. As an example, the theorem is applied to solve the problem of thermo-force bending of a thick plate.

The Stability of Runge－Kutta Methods for Systems of Delay Differential Equations

ＴｈｅＳｔａｂｉｌｉｔｙｏｆＲｕｎｇｅ－ＫｕｔｔａＭｅｔｈｏｄｓｆｏｒＳｙｓｔｅｍｓｏｆＤｅｌａｙＤｉｆｆｅｒｅｎｔｉａｌＥｑｕａｔｉｏｎｓ￥ＷＡＮＧＸｉａｏｂｉａｏ；ＬＩＵＭｉｎｇｚｈｕ(王晓彪)(刘明珠)（Ｄｅｐｔ．ｏｆＭａｔｈｅｍａｔｉｃｓ，Ｈａｒ...

Numerical Solution for Fractional Partial Differential Equation with Bernstein Polynomials

A framework to obtain numerical solution of the fractional partial differential equation using Bernstein polynomials is presented. The main characteristic behind this approach is that a fractional order operational matrix of Bernstein polynomials is derived. With the operational matrix, the equation is transformed into the products of several dependent matrixes which can also be regarded as the system of linear equations after dispersing the variable. By solving the linear equations, the numerical solutions are acquired. Only a small number of Bernstein polynomials are needed to obtain a satisfactory result. Numerical examples are provided to show that the method is computationally efficient.

SINE TRANSFORM PRECONDITIONERS FOR SECOND-ORDER PARTIAL DIFFERENTIAL EQUATIONS

In this paper, we are concerned with the numerical solution of second-order partial differential equations. We analyse the use of the Sine Transform precondilioners for the solution of linear systems arising from the discretization of p.d.e. via the preconditioned conjugate gradient method. For the second-order partial differential equations with Dirichlel boundary conditions, we prove that the condition number of the preconditioned system is O(1) while the condition number of the original system is O(m 2) Here m is the number of interior gridpoints in each direction. Such condition number produces a linear convergence rale.

Element-free Galerkin (EFG) method for analysis of the time-fractional partial differential equations

The present paper deals with the numerical solution of time-fractional partial differential equations using the element-free Galerkin （EFG） method, which is based on the moving least-square approximation. Compared with numerical methods based on meshes, the EFG method for time-fractional partial differential equations needs only scattered nodes instead of meshing the domain of the problem. It neither requires element connectivity nor suffers much degradation in accuracy when nodal arrangements are very irregular. In this method, the first-order time derivative is replaced by the Caputo fractional derivative of order α（0 〈 α≤ 1）. The Galerkin weak form is used to obtain the discrete equations, and the essential boundary conditions are enforced by the penalty method. Several numerical examples are presented and the results we obtained are in good agreement with the exact solutions.

The Multi-scale Method for Solving Nonlinear Time Space Fractional Partial Differential Equations

In this paper, we present a new algorithm to solve a kind of nonlinear time space-fractional partial differential equations on a finite domain. The method is based on B-spline wavelets approximations, some of these functions are reshaped to satisfy on boundary conditions exactly. The Adams fractional method is used to reduce the problem to a system of equations. By multiscale method this system is divided into some smaller systems which have less computations. We get an approximated solution which is more accurate on some subdomains by combining the solutions of these systems. Illustrative examples are included to demonstrate the validity and applicability of our proposed technique, also the stability of the method is discussed.