Analytical and numerical solutions of a two‐dimensional multi‐term time‐fractional Oldroyd‐B model

In this paper, we consider a two‐dimensional multi‐term time‐fractional Oldroyd‐B equation on a rectangular domain. Its analytical solution is obtained by the method of separation of variables. We employ the finite difference method with a discretization of the Caputo time‐fractional derivative to obtain an implicit difference approximation for the equation. Stability and convergence of the approximation scheme are established in the L_∞ ‐norm. Two examples are given to illustrate the theoretical analysis and analytical solution. The results indicate that the present numerical method is effective for this general two‐dimensional multi‐term time‐fractional Oldroyd‐B model.     

Numerical solution for (2 + 1)‐dimensional time‐fractional coupled Burger equations using fractional natural decomposition method

The aim of the present work is to find the numerical solutions for time‐fractional coupled Burgers equations using a new novel technique, called fractional natural decomposition method (FNDM). Two examples are considered in order to illustrate and validate the efficiency of the proposed algorithm. The numerical simulation has been conducted to ensure the exactness of the present method, and the obtained solutions are offered graphically to reveal the applicability and reliability of the FNDM. The outcomes of the study reveal that the FNDM is computationally very effective and accurate to study the (2 + 1)‐dimensional coupled Burger equations of arbitrary order.     

Numerical scheme with convergence for a generalized time‐fractional Telegraph‐type equation

In this work, we design and analyze a numerical scheme for solving the generalized time‐fractional Telegraph‐type equation (GTFTTE) which is defined using the generalized time fractional derivative (GTFD) proposed recently by Agrawal. The GTFD involves the scale and the weight functions, and reduces to the traditional Caputo derivative for a particular choice of the weight and the scale functions. The scale and the weight functions play an important role in describing the behavior of real‐life physical systems and thus we study the solution behavior of the GTFTTE by varying the weight and the scale functions in the GTFD. We investigate the solution profile of the GTFTTE under some of these choices. We also provide the stability and the convergence analysis of the proposed numerical scheme for the GTFTTE. We consider two test examples to perform numerical simulations.     

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.     

分数阶Burgers方程的有限元计算(英文)

建立了一维和二维分数阶Burgers方程的有限元格式.时间分数阶导数使用L1方法离散,空间方向使用有限元方法离散.通过选择合适的基函数,将离散后的方程转化成一个非线性代数方程组,并应用牛顿迭代方法求解.数值实验显示出了方法的有效性.     

Gegenbauer spectral method for time‐fractional convection–diffusion equations with variable coefficients

In this paper, we study the numerical solution to time‐fractional partial differential equations with variable coefficients that involve temporal Caputo derivative. A spectral method based on Gegenbauer polynomials is taken for approximating the solution of the given time‐fractional partial differential equation in time and a collocation method in space. The suggested method reduces this type of equation to the solution of a linear algebraic system. Finally, some numerical examples are presented to illustrate the efficiency and accuracy of the proposed method. Copyright © 2014 John Wiley & Sons, Ltd.     

Local discontinuous Galerkin approximations to fractional Bagley-Torvik equation

In this research, numerical approximation to fractional Bagley-Torvik equation as an important model arising in fluid mechanics is investigated. Our discretization algorithm is based on the local discontinuous Galerkin (LDG) schemes along with using the natural upwind fluxes, which enables us to solve the model problem element by element. This means that we require to solve a low-order system of equations in each subinterval, hence avoiding the need for a full global solution. The proposed schemes are tested on a range of initial- and boundary-value problems including a variable coefficient example, a nonsmooth problem, and some oscillatory test cases with exact solutions. Also, the validation of the proposed methods was compared with those obtained available existing computational procedures. Overall, it was found that LDG methods indicated highly satisfactory performance with comparatively lower degree of polynomials and number of elements compared with other numerical models.     

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     

Weak solutions of fractional differential equations in non cylindrical domains

We study a time fractional heat equation in a non cylindrical domain. The problem is one-dimensional. We prove existence of properly defined weak solutions by means of the Galerkin approximation.     

Well‐posedness of time‐space fractional stochastic evolution equations driven by α‐stable noise

This paper is devoted to the well‐posedness for time‐space fractional Ginzburg‐Landau equation and time‐space fractional Navier‐Stokes equations by α‐stable noise. The spatial regularity and the temporal regularity of the nonlocal stochastic convolution are firstly established, and then the existence and uniqueness of the global mild solution are obtained by the Banach fixed point theorem and Mittag‐Leffler functions, respectively. Numerical simulations for time‐space fractional Ginzburg‐Landau equation are provided to verify the analysis results.     

The dual reciprocity boundary elements method for the linear and nonlinear two‐dimensional time‐fractional partial differential equations

In this paper, we apply the dual reciprocity boundary elements method for the numerical solution of two‐dimensional linear and nonlinear time‐fractional modified anomalous subdiffusion equations and time‐fractional convection–diffusion equation. The fractional derivative of problems is described in the Riemann–Liouville and Caputo senses. We employ the linear radial basis function for interpolation of the nonlinear, inhomogeneous and time derivative terms. This method is improved by using a predictor–corrector scheme to overcome the nonlinearity which appears in the nonlinear problems under consideration. The accuracy and efficiency of the proposed schemes are checked by five test problems. The proposed method is employed for solving some examples in two dimensions on unit square and also in complex regions to demonstrate the efficiency of the new technique. Copyright © 2016 John Wiley & Sons, Ltd.     

Time discontinuous Galerkin space-time finite element method for nonlinear Sobolev equations

This article presents a complete discretization of a nonlinear Sobolev equation using space-time discontinuous Galerkin method that is discontinuous in time and continuous in space. The scheme is formulated by introducing the equivalent integral equation of the primal equation. The proposed scheme does not explicitly include the jump terms in time, which represent the discontinuity characteristics of approximate solution. And then the complexity of the theoretical analysis is reduced. The existence and uniqueness of the approximate solution and the stability of the scheme are proved. The optimalorder error estimates in L²(H¹) and L²(L²) norms are derived. These estimates are valid under weak restrictions on the space-time mesh, namely, without the condition k_n≥ch², which is necessary in traditional space-time discontinuous Galerkin methods. Numerical experiments are presented to verify the theoretical results.     

Convergence of the interpolated coefficient finite element method for the two‐dimensional elliptic sine‐Gordon equations

An interpolated coefficient finite element method is presented and analyzed for the two‐dimensional elliptic sine‐Gordon equations with Dirichlet boundary conditions. It is proved that the discretization scheme admits at least one solution, and that a subsequence of the approximation solutions converges to an exact solution in L²‐norm as the mesh size tends to zero. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

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     

A new method for exact solutions of variant types of time‐fractional Korteweg‐de Vries equations in shallow water waves

The current article devoted on the new method for finding the exact solutions of some time‐fractional Korteweg–de Vries (KdV) type equations appearing in shallow water waves. We employ the new method here for time‐fractional equations viz. time‐fractional KdV‐Burgers and KdV‐mKdV equations for finding the exact solutions. We use here the fractional complex transform accompanied by properties of local fractional calculus for reduction of fractional partial differential equations to ordinary differential equations. The obtained results are demonstrated by graphs for the new solutions. Copyright © 2016 John Wiley & Sons, Ltd.     

The finite element method for fractional diffusion with spectral fractional Laplacian

This paper focuses on the finite element method for Caputo-type parabolic equation with spectral fractional Laplacian, where the time derivative is in the sense of Caputo with order in (0,1) and the spatial derivative is the spectral fractional Laplacian. The time discretization is based on the Hadamard finite-part integral (or the finite-part integral in the sense of Hadamard), where the piecewise linear interpolation polynomials are used. The spatial fractional Laplacian is lifted to the local spacial derivative by using the Caffarelli–Silvestre extension, where the finite element method is used. Full-discretization scheme is constructed. The convergence and error estimates are obtained. Finally, numerical experiments are presented which support the theoretical results.     

A backward euler orthogonal spline collocation method for the time‐fractional Fokker–Planck equation

We formulate and analyze a novel numerical method for solving a time‐fractional Fokker–Planck equation which models an anomalous subdiffusion process. In this method, orthogonal spline collocation is used for the spatial discretization and the time‐stepping is done using a backward Euler method based on the L1 approximation to the Caputo derivative. The stability and convergence of the method are considered, and the theoretical results are supported by numerical examples, which also exhibit superconvergence. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1534–1550, 2015     

A fractional‐order Jacobi Tau method for a class of time‐fractional PDEs with variable coefficients

This paper presents a shifted fractional‐order Jacobi orthogonal function (SFJF) based on the definition of the classical Jacobi polynomial. A new fractional integral operational matrix of the SFJF is presented and derived. We propose the spectral Tau method, in conjunction with the operational matrices of the Riemann–Liouville fractional integral for SFJF and derivative for Jacobi polynomial, to solve a class of time‐fractional partial differential equations with variable coefficients. In this algorithm, the approximate solution is expanded by means of both SFJFs for temporal discretization and Jacobi polynomials for spatial discretization. The proposed tau scheme, both in temporal and spatial discretizations, successfully reduced such problem into a system of algebraic equations, which is far easier to be solved. Numerical results are provided to demonstrate the high accuracy and superiority of the proposed algorithm over existing ones. Copyright © 2015 John Wiley & Sons, Ltd.