An almost L‐stable BDF‐type method for the numerical solution of stiff ODEs arising from the method of lines

A new BDF‐type scheme is proposed for the numerical integration of the system of ordinary differential equations that arises in the Method of Lines solution of time‐dependent partial differential equations. This system is usually stiff, so it is desirable for the numerical method to solve it to have good properties concerning stability. The method proposed in this article is almost L‐stable and of algebraic order three. Numerical experiments illustrate the performance of the new method on different stiff systems of ODEs after discretizing in the space variable some PDE problems. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

The fully Sinc‐Galerkin method for time‐dependent boundary conditions

The fully Sinc‐Galerkin method is developed for a family of complex‐valued partial differential equations with time‐dependent boundary conditions. The Sinc‐Galerkin discrete system is formulated and represented by a Kronecker product form of those equations. The numerical solution is efficiently calculated and the method exhibits an exponential convergence rate. Several examples, some with a real‐valued solution and some with a complex‐valued solution, are used to demonstrate the performance of this method. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

Variational iteration method for solving the space‐ and time‐fractional KdV equation

This paper presents numerical solutions for the space‐ and time‐fractional Korteweg–de Vries equation (KdV for short) using the variational iteration method. The space‐ and time‐fractional derivatives are described in the Caputo sense. In this method, general Lagrange multipliers are introduced to construct correction functionals for the problems. The multipliers in the functionals can be identified optimally via variational theory. The iteration method, which produces the solutions in terms of convergent series with easily computable components, requiring no linearization or small perturbation. The numerical results show that the approach is easy to implement and accurate when applied to space‐ and time‐fractional KdV equations. The method introduces a promising tool for solving many space–time fractional partial differential equations. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

On the solution of an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

Numerical solution of hyperbolic partial differential equation with an integral condition continues to be a major research area with widespread applications in modern physics and technology. Many physical phenomena are modeled by nonclassical hyperbolic boundary value problems with nonlocal boundary conditions. In place of the classical specification of boundary data, we impose a nonlocal boundary condition. Partial differential equations with nonlocal boundary specifications have received much attention in last 20 years. However, most of the articles were directed to the second‐order parabolic equation, particularly to heat conduction equation. We will deal here with new type of nonlocal boundary value problem that is the solution of hyperbolic partial differential equations with nonlocal boundary specifications. These nonlocal conditions arise mainly when the data on the boundary can not be measured directly. Several finite difference methods have been proposed for the numerical solution of this one‐dimensional nonclassic boundary value problem. These computational techniques are compared using the largest error terms in the resulting modified equivalent partial differential equation. Numerical results supporting theoretical expectations are given. Restrictions on using higher order computational techniques for the studied problem are discussed. Suitable references on various physical applications and the theoretical aspects of solutions are introduced at the end of this article. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A Cartesian‐grid collocation method based on radial‐basis‐function networks for solving PDEs in irregular domains

This paper reports a new Cartesian‐grid collocation method based on radial‐basis‐function networks (RBFNs) for numerically solving elliptic partial differential equations in irregular domains. The domain of interest is embedded in a Cartesian grid, and the governing equation is discretized by using a collocation approach. The new features here are (a) one‐dimensional integrated RBFNs are employed to represent the variable along each line of the grid, resulting in a significant improvement of computational efficiency, (b) the present method does not require complicated interpolation techniques for the treatment of Dirichlet boundary conditions in order to achieve a high level of accuracy, and (c) normal derivative boundary conditions are imposed by means of integration constants. The method is verified through the solution of second‐ and fourth‐order PDEs; accurate results and fast convergence rates are obtained. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Radial basis collocation method and quasi‐Newton iteration for nonlinear elliptic problems

This work presents a radial basis collocation method combined with the quasi‐Newton iteration method for solving semilinear elliptic partial differential equations. The main result in this study is that there exists an exponential convergence rate in the radial basis collocation discretization and a superlinear convergence rate in the quasi‐Newton iteration of the nonlinear partial differential equations. In this work, the numerical error associated with the employed quadrature rule is considered. It is shown that the errors in Sobolev norms for linear elliptic partial differential equations using radial basis collocation method are bounded by the truncation error of the RBF. The combined errors due to radial basis approximation, quadrature rules, and quasi‐Newton and Newton iterations are also presented. This result can be extended to finite element or finite difference method combined with any iteration methods discussed in this work. The numerical example demonstrates a good agreement between numerical results and analytical predictions. The numerical results also show that although the convergence rate of order 1.62 of the quasi‐Newton iteration scheme is slightly slower than rate of order 2 in the Newton iteration scheme, the former is more stable and less sensitive to the initial guess. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Solving hyperbolic conservation laws using multiquadric quasi‐interpolation

In this article, we apply the univariate multiquadric (MQ) quasi‐interpolation to solve the hyperbolic conservation laws. At first we construct the MQ quasi‐interpolation corresponding to periodic and inflow‐outflow boundary conditions respectively. Next we obtain the numerical schemes to solve the partial differential equations, by using the derivative of the quasi‐interpolation to approximate the spatial derivative of the differential equation and a low‐order explicit difference to approximate the temporal derivative of the differential equation. Then we verify our scheme for the one‐dimensional Burgers' equation (without viscosity). We can see that the numerical results are very close to the exact solution and the computational accuracy of the scheme is ??(τ), where τ is the temporal step. We can improve the accuracy by using the high‐order quasi‐interpolation. Moreover the methods can be generalized to the other equations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

An efficient numerical approach to solve a class of variable‐order fractional integro‐partial differential equations

The main purpose of this work is to investigate an initial boundary value problem related to a suitable class of variable order fractional integro‐partial differential equations with a weakly singular kernel. To discretize the problem in the time direction, a finite difference method will be used. Then, the Sinc‐collocation approach combined with the double exponential transformation is employed to solve the problem in each time level. The proposed numerical algorithm is completely described and the convergence analysis of the numerical solution is presented. Finally, some illustrative examples are given to demonstrate the pertinent features of the proposed algorithm.     

Using empirical Eigenfunctions and Galerkin method to two‐phase transport models

In this article, we consider a nonlinear partial differential system describing two‐phase transports and try to recover the source term and the nonlinear diffusion term when the state variable is known at different profile times. To this end, we use a POD‐Galerkin procedure in which the proper orthogonal decomposition technique is applied to the ensemble of solutions to derive empirical eigenfunctions. These empirical eigenfunctions are then used as basis functions within a Galerkin method to transform the partial differential equation into a set of ordinary differential equations. Finally, the validation of the used method has been evaluated by some numerical examples. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 456–474, 2007     

Adaptive pseudo‐transient‐continuation‐Galerkin methods for semilinear elliptic partial differential equations

In this article, we investigate the application of pseudo‐transient‐continuation (PTC) schemes for the numerical solution of semilinear elliptic partial differential equations, with possible singular perturbations. We will outline a residual reduction analysis within the framework of general Hilbert spaces, and, subsequently, use the PTC‐methodology in the context of finite element discretizations of semilinear boundary value problems. Our approach combines both a prediction‐type PTC‐method (for infinite dimensional problems) and an adaptive finite element discretization (based on a robust a posteriori residual analysis), thereby leading to a fully adaptive PTC ‐Galerkin scheme. Numerical experiments underline the robustness and reliability of the proposed approach for different examples.© 2017 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 2005–2022, 2017     

Numerical scheme for solving system of fractional partial differential equations with Volterra‐type integral term through two‐dimensional block‐pulse functions

In the current study, an approximate scheme is established for solving the fractional partial differential equations (FPDEs) with Volterra integral terms via two‐dimensional block‐pulse functions (2D‐BPFs). According to the definitions and properties of 2D‐BPFs, the original problem is transformed into a system of linear algebra equations. By dispersing the unknown variables for these algebraic equations, the numerical solutions can be obtained. Besides, the proof of the convergence of this system is given. Finally, several numerical experiments are presented to test the feasibility and effectiveness of the proposed method.     

Use of radial basis functions for solving the second‐order parabolic equation with nonlocal boundary conditions

Nonlocal mathematical models appear in various problems of physics and engineering. In these models the integral term may appear in the boundary conditions. In this paper the problem of solving the one‐dimensional parabolic partial differential equation subject to given initial and nonlocal boundary conditions is considered. These kinds of problems have certainly been one of the fastest growing areas in various application fields. The presence of an integral term in a boundary condition can greatly complicate the application of standard numerical techniques. As a well‐known class of meshless methods, the radial basis functions are used for finding an approximation of the solution of the present problem. Numerical examples are given at the end of the paper to compare the efficiency of the radial basis functions with famous finite‐difference methods. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

Riesz basis and stabilization for the flexible structure of a symmetric tree‐shaped beam network

The stabilization of a symmetric tree‐shaped network of Euler–Bernoulli beams described by a system of partial differential equations is considered. The boundary controllers are designed based on passivity principle. The eigenfrequencies are analysed in detail and the asymptotic expansion of eigenvalues are presented. It is shown that there is a set of generalized eigenfunctions for the closed‐loop system, which forms a Riesz basis with parentheses for the energy state space. This concludes the spectrum‐determined growth condition and the exponential stability of the closed‐loop system. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical solution of the one‐dimensional wave equation with an integral condition

The hyperbolic partial differential equation with an integral condition arises in many physical phenomena. In this research a numerical technique is developed for the one‐dimensional hyperbolic equation that combine classical and integral boundary conditions. The proposed method is based on shifted Legendre tau technique. Illustrative examples are included to demonstrate the validity and applicability of the presented technique. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 282–292, 2007     

A new approach by two‐dimensional wavelets operational matrix method for solving variable‐order fractional partial integro‐differential equations

In this paper, a new computational scheme based on operational matrices (OMs) of two‐dimensional wavelets is proposed for the solution of variable‐order (VO) fractional partial integro‐differential equations (PIDEs). To accomplish this method, first OMs of integration and VO fractional derivative (FD) have been derived using two‐dimensional Legendre wavelets. By implementing two‐dimensional wavelets approximations and the OMs of integration and variable‐order fractional derivative (VO‐FD) along with collocation points, the VO fractional partial PIDEs are reduced into the system of algebraic equations. In addition to this, some useful theorems are discussed to establish the convergence analysis and error estimate of the proposed numerical technique. Furthermore, computational efficiency and applicability are examined through some illustrative examples.     

Two‐dimensional shifted Legendre polynomial collocation method for electromagnetic waves in dielectric media via almost operational matrices

In this paper, a numerical solution of fractional partial differential equations (FPDEs) for electromagnetic waves in dielectric media will be discussed. For the solution of FPDEs, we developed a numerical collocation method using an algorithm based on two‐dimensional shifted Legendre polynomials approximation, which is proposed for electromagnetic waves in dielectric media. By implementing the partial Riemann–Liouville fractional derivative operators, two‐dimensional shifted Legendre polynomials approximation and its operational matrix along with collocation method are used to convert FPDEs first into weakly singular fractional partial integro‐differential equations and then converted weakly singular fractional partial integro‐differential equations into system of algebraic equation. Some results concerning the convergence analysis and error analysis are obtained. Illustrative examples are included to demonstrate the validity and applicability of the technique. Copyright © 2017 John Wiley & Sons, Ltd.     

Minimum time‐step criteria for the Galerkin finite element methods applied to one‐dimensional parabolic partial differential equations

The finite element method has been well established for numerically solving parabolic partial differential equations (PDEs). Also it is well known that a too large time step should not be chosen in order to obtain a stable and accurate numerical solution. In this article, accuracy analysis shows that a too small time step should not be chosen either for some time‐stepping schemes. Otherwise, the accuracy of the numerical solution cannot be improved or can even be worsened in some cases. Furthermore, the so‐called minimum time step criteria are established for the Crank‐Nicolson scheme, the Galerkin‐time scheme, and the backward‐difference scheme used in the temporal discretization. For the forward‐difference scheme, no minimum time step exists as far as the accuracy is concerned. In the accuracy analysis, no specific initial and boundary conditions are invoked so that such established criteria can be applied to the parabolic PDEs subject to any initial and boundary conditions. These minimum time step criteria are verified in a series of numerical experiments for a one‐dimensional transient field problem with a known analytical solution. The minimum time step criteria developed in this study are useful for choosing appropriate time steps in numerical simulations of practical engineering problems. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

An Eulerian‐Lagrangian single‐node collocation method for transient advection‐diffusion equations in multiple space dimensions

We developed a nonconventional Eulerian‐Lagrangian single‐node collocation method for transient advection‐diffusion transport partial differential equations in multiple space dimensions. This method greatly reduces the number of unknowns in conventional collocation method, generates accurate numerical solutions, and allows large time steps to be used in numerical simulations. We perform numerical experiments to show the strong potential of the method. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 20: 284–301, 2004     

Bernstein Ritz‐Galerkin method for solving an initial‐boundary value problem that combines Neumann and integral condition for the wave equation

In this article, the Ritz‐Galerkin method in Bernstein polynomial basis is implemented to give an approximate solution of a hyperbolic partial differential equation with an integral condition. We will deal here with a type of nonlocal boundary value problem, that is, the solution of a hyperbolic partial differential equation with a nonlocal boundary specification. The nonlocal conditions arise mainly when the data on the boundary cannot be measured directly. The properties of Bernstein polynomial and Ritz‐Galerkin method are first presented, then Ritz‐Galerkin method is used to reduce the given hyperbolic partial differential equation to the solution of algebraic equations. Illustrative examples are included to demonstrate the validity and applicability of the technique presented in this article. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

A second‐order linearized difference scheme for a strongly coupled reaction‐diffusion system

This article deals with the numerical solution to some models described by the system of strongly coupled reaction–diffusion equations with the Neumann boundary value conditions. A linearized three‐level scheme is derived by the method of reduction of order. The uniquely solvability and second‐order convergence in L₂‐norm are proved by the energy method. A numerical example is presented to demonstrate the accuracy and efficiency of the proposed method. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007