A discontinuous Galerkin finite element method for time dependent partial differential equations with higher order derivatives

In this paper, we develop a new discontinuous Galerkin (DG) finite element method for solving time dependent partial differential equations (PDEs) with higher order spatial derivatives. Unlike the traditional local discontinuous Galerkin (LDG) method, the method in this paper can be applied without introducing any auxiliary variables or rewriting the original equation into a larger system. Stability is ensured by a careful choice of interface numerical fluxes. The method can be designed for quite general nonlinear PDEs and we prove stability and give error estimates for a few representative classes of PDEs up to fifth order. Numerical examples show that our scheme attains the optimal -th order of accuracy when using piecewise -th degree polynomials, under the condition that is greater than or equal to the order of the equation.

     

Bernstein modal basis: Application to the spectral Petrov‐Galerkin method for fractional partial differential equations

In the spectral Petrov‐Galerkin methods, the trial and test functions are required to satisfy particular boundary conditions. By a suitable linear combination of orthogonal polynomials, a basis, that is called the modal basis, is obtained. In this paper, we extend this idea to the nonorthogonal dual Bernstein polynomials. A compact general formula is derived for the modal basis functions based on dual Bernstein polynomials. Then, we present a Bernstein‐spectral Petrov‐Galerkin method for a class of time fractional partial differential equations with Caputo derivative. It is shown that the method leads to banded sparse linear systems for problems with constant coefficients. Some numerical examples are provided to show the efficiency and the spectral accuracy of the method.     

Stability and error analysis for a spectral Galerkin method for the Navier‐Stokes equations with H2 or H1 initial data

In this article we consider a spectral Galerkin method with a semi‐implicit Euler scheme for the two‐dimensional Navier‐Stokes equations with H² or H¹ initial data. The H²‐stability analysis of this spectral Galerkin method shows that for the smooth initial data the semi‐implicit Euler scheme admits a large time step. The L²‐error analysis of the spectral Galerkin method shows that for the smoother initial data the numerical solution u exhibits faster convergence on the time interval [0, 1] and retains the same convergence rate on the time interval [1, ∞). © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005.     

Least-squares spectral element method for non-linear hyperbolic differential equations

The least-squares spectral element method has been applied to the one-dimensional inviscid Burgers equation which allows for discontinuous solutions. In order to achieve high order accuracy both in space and in time a space–time formulation has been applied. The Burgers equation has been discretized in three different ways: a non-conservative formulation, a conservative system with two variables and two equations: one first order linear PDE and one linearized algebraic equation, and finally a variant on this conservative formulation applied to a direct minimization with a QR-decomposition at elemental level. For all three formulations an h/p-convergence study has been performed and the results are discussed in this paper.     

Legendre–Gauss–Lobatto spectral collocation method for nonlinear delay differential equations

A Legendre–Gauss–Lobatto spectral collocation method is introduced for the numerical solutions of a class of nonlinear delay differential equations. An efficient algorithm is designed for the single‐step scheme and applied to the multiple‐domain case. As a theoretical result, we obtain a general convergence theorem for the single‐step case. Numerical results show that the suggested algorithm enjoys high‐order accuracy both in time and in the delayed argument and can be implemented in a robust and efficient manner. Copyright © 2013 John Wiley & Sons, Ltd.     

Legendre spectral Galerkin method for second-kind Volterra integral equations

The Legendre spectral Galerkin method for the Volterra integral equations of the second kind is proposed in this paper. We provide a rigorous error analysis for the proposed method, which indicates that the numerical errors (in the L ₂ norm) will decay exponentially provided that the kernel function and the source function are sufficiently smooth. Numerical examples are given to illustrate the theoretical results.      

Galerkin finite element approximation of symmetric space-fractional partial differential equations

In this paper, symmetric space-fractional partial differential equations (SSFPDE) with the Riesz fractional operator are considered. The SSFPDE is obtained from the standard advection-dispersion equation by replacing the first-order and second-order space derivatives with the Riesz fractional derivatives of order 2β ∈ (0, 1) and 2α ∈ (1, 2], respectively. We prove that the variational solution of the SSFPDE exists and is unique. Using the Galerkin finite element method and a backward difference technique, a fully discrete approximating system is obtained, which has a unique solution according to the Lax-Milgram theorem. The stability and convergence of the fully discrete schemes are derived. Finally, some numerical experiments are given to confirm our theoretical analysis.     

High-order finite element methods for time-fractional partial differential equations

The aim of this paper is to develop high-order methods for solving time-fractional partial differential equations. The proposed high-order method is based on high-order finite element method for space and finite difference method for time. Optimal convergence rate O((Δt)^2−α+N^−r) is proved for the (r−1)th-order finite element method (r≥2).     

On discontinuous Galerkin finite element method for singularly perturbed delay differential equations

A nonsymmetric discontinuous Galerkin FEM with interior penalties has been applied to one-dimensional singularly perturbed problem with a constant negative shift. Using higher order polynomials on Shishkin-type layer-adapted meshes, a robust convergence has been proved in the corresponding energy norm. Numerical experiments support theoretical findings.     

Stability and error analysis for spectral Galerkin method for the Navier–Stokes equations with L2 initial data

In this article we consider the spectral Galerkin method with the implicit/explicit Euler scheme for the two‐dimensional Navier–Stokes equations with the L² initial data. Due to the poor smoothness of the solution on [0,1), we use the the spectral Galerkin method based on high‐dimensional spectral space H_M and small time step Δt² on this interval. While on [1,∞), we use the spectral Galerkin method based on low‐dimensional spectral space H_m(m = O(M^1/2)) and large time step Δt. For the spectral Galerkin method, we provide the standard H²‐stability and the L²‐error analysis. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2007     

Exact solutions of nonlinear time fractional partial differential equations by sub‐equation method

In this article, the sub‐equation method is presented for finding the exact solutions of a nonlinear fractional partial differential equations. For this, the fractional complex transformation method has been used to convert fractional‐order partial differential equation to ordinary differential equation. The fractional derivatives are described in Jumarie's the modified Riemann–Liouville sense. We apply to this method for the nonlinear time fractional differential equations. With the aid of symbolic computation, a variety of exact solutions for them are obtained. Copyright © 2014 John Wiley & Sons, Ltd.     

Reproducing kernel method to solve parabolic partial differential equations with nonlocal conditions

In this study, the parabolic partial differential equations with nonlocal conditions are solved. To this end, we use the reproducing kernel method (RKM) that is obtained from the combining fundamental concepts of the Galerkin method, and the complete system of reproducing kernel Hilbert space that was first introduced by Wang et al. who implemented RKM without Gram–Schmidt orthogonalization process. In this method, first the reproducing kernel spaces and their kernels such that satisfy the nonlocal conditions are constructed, and then the RKM without Gram–Schmidt orthogonalization process on the considered problem is implemented. Moreover, convergence theorem, error analysis theorems, and stability theorem are provided in detail. To show the high accuracy of the present method several numerical examples are solved.     

Time accurate fast wavelet‐Taylor Galerkin method for partial differential equations

We introduce the concept of fast wavelet‐Taylor Galerkin methods for the numerical solution of partial differential equations. In wavelet‐Taylor Galerkin method discretization in time is performed before the wavelet based spatial approximation by introducing accurate generalizations of the standard Euler, θ and leap‐frog time‐stepping scheme with the help of Taylor series expansions in the time step. We will present two different time‐accurate wavelet schemes to solve the PDEs. First, numerical schemes taking advantage of the wavelet bases capabilities to compress the operators and sparse representation of functions which are smooth, except for in localized regions, up to any given accuracy are presented. Here numerical experiments deal with advection equation with the spiky solution in one dimension, two dimensions, and nonlinear equation with a shock in solution in two dimensions. Second, our schemes deal with more regular class of problems where wavelets are not efficient procedure for data compression but we can use the good approximation properties of wavelet. Here time‐accurate schemes lead to consistent mass matrix in an explicit time stepping, which can be solved by approximate factorization techniques. Numerical experiment deals with more regular class of problems like heat equation as well as coupled linear system in two dimensions. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Waveform relaxation method for stochastic differential equations with constant delay

This paper extends the waveform relaxation method to stochastic differential equations with constant delay terms, gives sufficient conditions for the mean square convergence of the method. A lot of attention is paid to the rate of convergence of the method. The conditions of the superlinear convergence for a special case, which bases on the special splitting functions, are given. The theory is applied to a one-dimensional model problem and checked against results obtained by numerical experiments.     

Galerkin approximations with embedded boundary conditions for retarded delay differential equations

Finite-dimensional approximations are developed for retarded delay differential equations (DDEs). The DDE system is equivalently posed as an initial-boundary value problem consisting of hyperbolic partial differential equations (PDEs). By exploiting the equivalence of partial derivatives in space and time, we develop a new PDE representation for the DDEs that is devoid of boundary conditions. The resulting boundary condition-free PDEs are discretized using the Galerkin method with Legendre polynomials as the basis functions, whereupon we obtain a system of ordinary differential equations (ODEs) that is a finite-dimensional approximation of the original DDE system. We present several numerical examples comparing the solution obtained using the approximate ODEs to the direct numerical simulation of the original non-linear DDEs. Stability charts developed using our method are compared to existing results for linear DDEs. The presented results clearly demonstrate that the equivalent boundary condition-free PDE formulation accurately captures the dynamic behaviour of the original DDE system and facilitates the application of control theory developed for systems governed by ODEs.     

Quasilinear parabolic stochastic partial differential equations: Existence,uniqueness

We present a direct approach to existence and uniqueness of strong (in the probabilistic sense) and weak (in the PDE sense) solutions to quasilinear stochastic partial differential equations, which are neither monotone nor locally monotone. The proof of uniqueness is very elementary, based on a new method of applying Itô’s formula for the $L^1$-norm. The proof of existence relies on a recent regularity result and is direct in the sense that it does not rely on the stochastic compactness method.     

Flow invariance for solutions to nonlinear nonautonomous partial differential delay equations

We investigate the problem of existence and flow invariance of mild solutions to nonautonomous partial differential delay equations , t?s, us=φ, where B(t) is a family of nonlinear multivalued, α-accretive operators with D(B(t)) possibly depending on t, and the operators F(t,.) being defined—and Lipschitz continuous—possibly only on “thin” subsets of the initial history space E. The results are applied to population dynamics models. We also study the asymptotic behavior of solutions to this equation. Our analysis will be based on the evolution operator associated to the equation in the initial history space E.     

Efficient hybrid method for solving special type of nonlinear partial differential equations

In this article, an efficient hybrid method has been developed for solving some special type of nonlinear partial differential equations. Hybrid method is based on tanh–coth method, quasilinearization technique and Haar wavelet method. Nonlinear partial differential equations have been converted into a nonlinear ordinary differential equation by choosing some suitable variable transformations. Quasilinearization technique is used to linearize the nonlinear ordinary differential equation and then the Haar wavelet method is applied to linearized ordinary differential equation. A tanh–coth method has been used to obtain the exact solutions of nonlinear ordinary differential equations. It is easier to handle nonlinear ordinary differential equations in comparison to nonlinear partial differential equations. A distinct feature of the proposed method is their simple applicability in a variety of two‐ and three‐dimensional nonlinear partial differential equations. Numerical examples show better accuracy of the proposed method as compared with the methods described in past. Error analysis and stability of the proposed method have been discussed.     

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