Numerical simulation of a class of fractional subdiffusion equations via the alternating direction implicit method

In this article, a new numerical technique is proposed for solving the two‐dimensional time fractional subdiffusion equation with nonhomogeneous terms. After a transformation of the original problem, standard central difference approximation is used for the spatial discretization. For the time step, a new fractional alternating direction implicit (FADI) scheme based on the L₁ approximation is considered. This FADI scheme is constructed by adding a small term, so it is different from standard FADI methods. The solvability, unconditional stability and H¹ norm convergence are proved. Finally, numerical examples show the effectiveness and accuracy of our proposed method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 531–547, 2016     

Supraconvergence and Supercloseness of a Scheme for Elliptic Equations on Nonuniform Grids

In this paper, we study the convergence of a finite difference scheme on nonuniform grids for the solution of second-order elliptic equations with mixed derivatives and variable coefficients in polygonal domains subjected to Dirichlet boundary conditions. We show that the scheme is equivalent to a fully discrete linear finite element approximation with quadrature. It exhibits the phenomenon of supraconvergence, more precisely, for s ∈ [1,2] order O(h ^s )-convergence of the finite difference solution, and its gradient is shown if the exact solution is in the Sobolev space H ^1+s (Ω). In the case of an equation with mixed derivatives in a domain containing oblique boundary sections, the convergence order is reduced to O(h ^3/2?ε) with ε > 0 if u ∈ H ³(Ω). The second-order accuracy of the finite difference gradient is in the finite element context nothing else than the supercloseness of the gradient. For s ∈ {1,2}, the given error estimates are strictly local.     

An Itô Formula of Generalized Functionals and Local Time for Fractional Brownian Sheet

Abstract

By using the white noise theory for a fractional Brownian sheet, we derive an Itô formula for the generalized functionals for the fractional Brownian sheet with arbitrary Hurst parameters H ₁, H ₂ ∈ (0,1). As an application, we give the integral representations for two versions of local times of a fractional Brownian sheet, respectively.     

Wick-Itô Formula for Gaussian Processes

Abstract

A Wick-Itô formula for Gaussian processes is obtained. This is a change of variables formula, which is to Wick-Itô integrals what the usual Itô formula is to Itô integrals. The conditions are weak enough to allow processes with infinite quadratic variation. They are satisfied by fractional Brownian motion with parameter 1/4 < H < 1.     

Stochastic Differential Equations Driven by Fractional Brownian Motion and Standard Brownian Motion

Abstract

We prove an existence and uniqueness theorem for solutions of multidimensional, time dependent, stochastic differential equations driven simultaneously by a multidimensional fractional Brownian motion with Hurst parameter H > 1/2 and a multidimensional standard Brownian motion. The proof relies on some a priori estimates, which are obtained using the methods of fractional integration and the classical Itô stochastic calculus. The existence result is based on the Yamada–Watanabe theorem.     

Exact Rate of Convergence of Some Approximation Schemes Associated to SDEs Driven by a Fractional Brownian Motion

In this article, we derive the exact rate of convergence of some approximation schemes associated to scalar stochastic differential equations driven by a fractional Brownian motion with Hurst index H. We consider two cases. If H>1/2, the exact rate of convergence of the Euler scheme is determined. We show that the error of the Euler scheme converges almost surely to a random variable, which in particular depends on the Malliavin derivative of the solution. This result extends those contained in J. Complex. 22(4), 459–474, 2006 and C.R. Acad. Sci. Paris, Ser. I 340(8), 611–614, 2005. When 1/6<H<1/2, the exact rate of convergence of the Crank-Nicholson scheme is determined for a particular equation. Here we show convergence in law of the error to a random variable, which depends on the solution of the equation and an independent Gaussian random variable.     

CONVERGENCE RATE OF AN APPROXIMATION TO MULTIPLE INTEGRAL OF FBM

In this article, we study the rate of convergence of the polygonal approximation to multiple stochastic integral Sp （f） of fractional Brownian motion of Hurst parameter H 〈 1/2 when the fractional Brownian motion is replaced by its polygonal approximation. Under different conditions on f and for different p, we obtain different rates.     

Variational formulation for the stationary fractional advection dispersion equation

In this article a theoretical framework for the Galerkin finite element approximation to the steady state fractional advection dispersion equation is presented. Appropriate fractional derivative spaces are defined and shown to be equivalent to the usual fractional dimension Sobolev spaces H^s. Existence and uniqueness results are proven, and error estimates for the Galerkin approximation derived. Numerical results are included that confirm the theoretical estimates. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

Extension Problem and Harnack's Inequality for Some Fractional Operators

The fractional Laplacian can be obtained as a Dirichlet-to-Neumann map via an extension problem to the upper half space. In this paper we prove the same type of characterization for the fractional powers of second order partial differential operators in some class. We also get a Poisson formula and a system of Cauchy–Riemann equations for the extension. The method is applied to the fractional harmonic oscillator H ^σ = (? Δ + |x|²)^σ to deduce a Harnack's inequality. A pointwise formula for H ^σ f(x) and some maximum and comparison principles are derived.     

Microlocal properties of scattering matrices

We consider the scattering theory for a pair of operators H₀ and H = H₀ + V on L²(M, m), where M is a Riemannian manifold, H₀ is a multiplication operator on M, and V is a pseudodifferential operator of order ? μ, μ > 1. We show that a time-dependent scattering theory can be constructed, and the scattering matrix is a pseudodifferential operator on each energy surface. Moreover, the principal symbol of the scattering matrix is given by a Born approximation type function. The main motivation of the study comes from applications to discrete Schrödigner operators, and it also applies to various differential operators with constant coefficients and short-range perturbations on Euclidean spaces.     

A nonuniform L2 formula of Caputo derivative and its application to a fractional Benjamin–Bona–Mahony-type equation with nonsmooth solutions

To recover the full accuracy of discretized fractional derivatives, nonuniform mesh technique is a natural and simple approach to efficiently resolve the initial singularities that always appear in the solutions of time-fractional linear and nonlinear differential equations. We first construct a nonuniform L2 approximation for the fractional Caputo's derivative of order 1 < α < 2 and present a global consistency analysis under some reasonable regularity assumptions. The temporal nonuniform L2 formula is then utilized to develop a linearized difference scheme for a time-fractional Benjamin–Bona–Mahony-type equation. The unconditional convergence of our scheme on both uniform and nonuniform (graded) time meshes are proven with respect to the discrete H¹-norm. Numerical examples are provided to justify the accuracy.     

Fourier-Integral-Operator Approximation of Solutions to First-Order Hyperbolic Pseudodifferential Equations I: Convergence in Sobolev Spaces

An approximation Ansatz for the operator solution, U(z′,z), of a hyperbolic first-order pseudodifferential equation, ? _z  + a(z,x,D _x ) with Re (a) ≥ 0, is constructed as the composition of global Fourier integral operators with complex phases. An estimate of the operator norm in L(H ^(s),H ^(s)) of these operators is provided, which yields a convergence result for the Ansatz to U(z′,z) in some Sobolev space as the number of operators in the composition goes to ∞.     

A new approximation for the generalized fractional derivative and its application to generalized fractional diffusion equation

In this paper, we derive a fourth order approximation for the generalized fractional derivative that is characterized by a scale function z(t) and a weight function w(t) . Combining the new approximation with compact finite difference method, we develop a numerical scheme for a generalized fractional diffusion problem. The stability and convergence of the numerical scheme are proved by the energy method, and it is shown that the temporal and spatial convergence orders are both 4. Several numerical experiments are provided to illustrate the efficiency of our scheme.     

Two temporal second-order <Emphasis Type="Italic">H</Emphasis><Superscript>1</Superscript>-Galerkin mixed finite element schemes for distributed-order fractional sub-diffusion equations

In this paper, a special point is found for the interpolation approximation of the distributed order fractional derivatives to achieve at least second-order accuracy. Then, two H¹-Galerkin mixed finite element schemes combined with the higher accurate interpolation approximation are introduced and analyzed to solve the distributed order fractional sub-diffusion equations. The stable results, which just depend on initial value and source item, are derived. Some a priori estimates with optimal order of convergence both for the unknown function and its flux are established rigorously. It is shown that the H¹-Galerkin mixed finite element approximations have the same rates of the convergence as in the classical mixed finite element method, but without LBB consistency condition and quasiuniformity requirement on the finite element mesh. Finally, some numerical experiments are presented to show the efficiency and accuracy of H¹-Galerkin mixed finite element schemes.     

A General Approximation Scheme for Attractors of Abstract Parabolic Problems

We consider semilinear problems of the form u′ = Au + f(u), where A generates an exponentially decaying compact analytic C ₀-semigroup in a Banach space E, f:E ^α → E is differentiable globally Lipschitz and bounded (E ^α = D((?A)^α) with the graph norm). Under a very general approximation scheme, we prove that attractors for such problems behave upper semicontinuously. If all equilibrium points are hyperbolic, then there is an odd number of them. If, in addition, all global solutions converge as t → ±∞, then the attractors behave lower semicontinuously. This general approximation scheme includes finite element method, projection and finite difference methods. The main assumption on the approximation is the compact convergence of resolvents, which may be applied to many other problems not related to discretization.     

A Discontinuous Galerkin Method for Fractional Differential Equations

We investigate the numerical treatment of fractional derivatives using a discontinuous Galerkin approach. We consider fractional differential equations in one dimension. The most common definitions, the Riemann-Liouville and Caputo derivatives, as well as a third definition are considered. Different types of initial condition statements are investigated for the extended discontinuous Galerkin method. We obtain numerical results on the order of convergence in the L₂ and H¹ norms. We compare the classical discontinuous Galerkin method with our adapted approach and show results for the different types of fractional derivatives. Further work was carried out on parabolic differential equations with fractional time derivatives, which is to be published cf. [2]. A keystone in the further research is the development of analytical tools for theoretical convergence results. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

The p-version of the boundary element method for hypersingular operators on piecewise plane open surfaces

Summary. We prove an optimal a priori error estimate for the p-version of the boundary element method with hypersingular operators on piecewise plane open surfaces. The solutions of problems on open surfaces typically exhibit a singular behavior at the edges and corners of the surface which prevent an approximation analysis in H¹. We analyze the approximation by polynomials of typical singular functions in fractional order Sobolev spaces thus giving, as an application, the optimal rate of convergence of the p-version of the boundary element method. This paper extends the results of [C. Schwab, M. Suri, The optimal p-version approximation of singularities on polyhedra in the boundary element method, SIAM J. Numer. Anal., 33 (1996), pp. 729–759] who only considered closed surfaces where the solution is in H¹.Mathematics Subject Classification (2000): 41A10, 65N15, 65N38Financed by the FONDAP Program in Applied Mathematics, Chile.Supported by the FONDAP Program in Applied Mathematics and Fondecyt project no. 1010220, both Chile.     

Reaction-Diffusion Equations with Polynomial Drifts Driven by Fractional Brownian Motions

A reaction-diffusion equation on [0, 1] ^d with the heat conductivity κ > 0, a polynomial drift term and an additive noise, fractional in time with H > 1/2, and colored in space, is considered. We have shown the existence, uniqueness and uniform boundedness of solution with respect to κ. Also we show that if κ tends to infinity, then the corresponding solutions of the equation converge to a process satisfying a stochastic ordinary differential equation.