Error analysis of a FVEM for fractional order evolution equations with nonsmooth initial data

Samir Karaa, Amiya K. Pani

Research output: Contribution to journalArticle

6 Citations (Scopus)

Abstract

In this paper, a finite volume element (FVE) method is considered for spatial approximations of time fractional diffusion equations involving a Riemann-Liouville fractional derivative of order α ∈ (0, 1) in time. Improving upon earlier results [Karaa et al., IMA J. Numer. Anal. 37 (2017) 945-964], error estimates in L2(Ω)- and H1(Ω)-norms for the semidiscrete problem with smooth and mildly smooth initial data, i.e., v ∈ H2(Ω) ∩ H1 0 (Ω) and v ∈ H1 0 (Ω) are established. For nonsmooth data, that is, v ∈ L2(Ω), the optimal L2(Ω)-error estimate is shown to hold only under an additional assumption on the triangulation, which is known to be satisfied for symmetric triangulations. Super-convergence result is also proved and as a consequence, a quasi-optimal error estimate is established in the L(Ω)-norm. Further, two fully discrete schemes using convolution quadrature in time generated by the backward Euler and the second-order backward difference methods are analyzed, and error estimates are derived for both smooth and nonsmooth initial data. Based on a comparison of the standard Galerkin finite element solution with the FVE solution and exploiting tools for Laplace transforms with semigroup type properties of the FVE solution operator, our analysis is then extended in a unified manner to several time fractional order evolution problems. Finally, several numerical experiments are conducted to confirm our theoretical findings.

Original languageEnglish
Pages (from-to)773-801
Number of pages29
JournalESAIM: Mathematical Modelling and Numerical Analysis
Volume52
Issue number2
DOIs
Publication statusPublished - Mar 1 2018

Fingerprint

Fractional Order
Error Analysis
Error analysis
Evolution Equation
Finite Volume Element
Optimal Error Estimates
Triangulation
Error Estimates
Finite Volume Element Method
Riemann-Liouville Fractional Derivative
Norm
Fractional Diffusion Equation
Evolution Problems
Superconvergence
Laplace transforms
Finite Element Solution
Convolution
Quadrature
Galerkin
Laplace transform

Keywords

  • Backward Euler and second-order backward difference methods
  • Convolution quadrature
  • Finite volume element method
  • Fractional order evolution equation
  • Laplace transform
  • Optimal error estimate
  • Smooth and nonsmooth data
  • Subdiffusion

ASJC Scopus subject areas

  • Analysis
  • Numerical Analysis
  • Modelling and Simulation
  • Applied Mathematics

Cite this

Error analysis of a FVEM for fractional order evolution equations with nonsmooth initial data. / Karaa, Samir; Pani, Amiya K.

In: ESAIM: Mathematical Modelling and Numerical Analysis, Vol. 52, No. 2, 01.03.2018, p. 773-801.

Research output: Contribution to journalArticle

@article{5646b8d9fed6468ab40982dc139ac6a4,
title = "Error analysis of a FVEM for fractional order evolution equations with nonsmooth initial data",
abstract = "In this paper, a finite volume element (FVE) method is considered for spatial approximations of time fractional diffusion equations involving a Riemann-Liouville fractional derivative of order α ∈ (0, 1) in time. Improving upon earlier results [Karaa et al., IMA J. Numer. Anal. 37 (2017) 945-964], error estimates in L2(Ω)- and H1(Ω)-norms for the semidiscrete problem with smooth and mildly smooth initial data, i.e., v ∈ H2(Ω) ∩ H1 0 (Ω) and v ∈ H1 0 (Ω) are established. For nonsmooth data, that is, v ∈ L2(Ω), the optimal L2(Ω)-error estimate is shown to hold only under an additional assumption on the triangulation, which is known to be satisfied for symmetric triangulations. Super-convergence result is also proved and as a consequence, a quasi-optimal error estimate is established in the L∞(Ω)-norm. Further, two fully discrete schemes using convolution quadrature in time generated by the backward Euler and the second-order backward difference methods are analyzed, and error estimates are derived for both smooth and nonsmooth initial data. Based on a comparison of the standard Galerkin finite element solution with the FVE solution and exploiting tools for Laplace transforms with semigroup type properties of the FVE solution operator, our analysis is then extended in a unified manner to several time fractional order evolution problems. Finally, several numerical experiments are conducted to confirm our theoretical findings.",
keywords = "Backward Euler and second-order backward difference methods, Convolution quadrature, Finite volume element method, Fractional order evolution equation, Laplace transform, Optimal error estimate, Smooth and nonsmooth data, Subdiffusion",
author = "Samir Karaa and Pani, {Amiya K.}",
year = "2018",
month = "3",
day = "1",
doi = "10.1051/m2an/2018029",
language = "English",
volume = "52",
pages = "773--801",
journal = "ESAIM: Mathematical Modelling and Numerical Analysis",
issn = "0764-583X",
publisher = "EDP Sciences",
number = "2",

}

TY - JOUR

T1 - Error analysis of a FVEM for fractional order evolution equations with nonsmooth initial data

AU - Karaa, Samir

AU - Pani, Amiya K.

PY - 2018/3/1

Y1 - 2018/3/1

N2 - In this paper, a finite volume element (FVE) method is considered for spatial approximations of time fractional diffusion equations involving a Riemann-Liouville fractional derivative of order α ∈ (0, 1) in time. Improving upon earlier results [Karaa et al., IMA J. Numer. Anal. 37 (2017) 945-964], error estimates in L2(Ω)- and H1(Ω)-norms for the semidiscrete problem with smooth and mildly smooth initial data, i.e., v ∈ H2(Ω) ∩ H1 0 (Ω) and v ∈ H1 0 (Ω) are established. For nonsmooth data, that is, v ∈ L2(Ω), the optimal L2(Ω)-error estimate is shown to hold only under an additional assumption on the triangulation, which is known to be satisfied for symmetric triangulations. Super-convergence result is also proved and as a consequence, a quasi-optimal error estimate is established in the L∞(Ω)-norm. Further, two fully discrete schemes using convolution quadrature in time generated by the backward Euler and the second-order backward difference methods are analyzed, and error estimates are derived for both smooth and nonsmooth initial data. Based on a comparison of the standard Galerkin finite element solution with the FVE solution and exploiting tools for Laplace transforms with semigroup type properties of the FVE solution operator, our analysis is then extended in a unified manner to several time fractional order evolution problems. Finally, several numerical experiments are conducted to confirm our theoretical findings.

AB - In this paper, a finite volume element (FVE) method is considered for spatial approximations of time fractional diffusion equations involving a Riemann-Liouville fractional derivative of order α ∈ (0, 1) in time. Improving upon earlier results [Karaa et al., IMA J. Numer. Anal. 37 (2017) 945-964], error estimates in L2(Ω)- and H1(Ω)-norms for the semidiscrete problem with smooth and mildly smooth initial data, i.e., v ∈ H2(Ω) ∩ H1 0 (Ω) and v ∈ H1 0 (Ω) are established. For nonsmooth data, that is, v ∈ L2(Ω), the optimal L2(Ω)-error estimate is shown to hold only under an additional assumption on the triangulation, which is known to be satisfied for symmetric triangulations. Super-convergence result is also proved and as a consequence, a quasi-optimal error estimate is established in the L∞(Ω)-norm. Further, two fully discrete schemes using convolution quadrature in time generated by the backward Euler and the second-order backward difference methods are analyzed, and error estimates are derived for both smooth and nonsmooth initial data. Based on a comparison of the standard Galerkin finite element solution with the FVE solution and exploiting tools for Laplace transforms with semigroup type properties of the FVE solution operator, our analysis is then extended in a unified manner to several time fractional order evolution problems. Finally, several numerical experiments are conducted to confirm our theoretical findings.

KW - Backward Euler and second-order backward difference methods

KW - Convolution quadrature

KW - Finite volume element method

KW - Fractional order evolution equation

KW - Laplace transform

KW - Optimal error estimate

KW - Smooth and nonsmooth data

KW - Subdiffusion

UR - http://www.scopus.com/inward/record.url?scp=85049005291&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85049005291&partnerID=8YFLogxK

U2 - 10.1051/m2an/2018029

DO - 10.1051/m2an/2018029

M3 - Article

AN - SCOPUS:85049005291

VL - 52

SP - 773

EP - 801

JO - ESAIM: Mathematical Modelling and Numerical Analysis

JF - ESAIM: Mathematical Modelling and Numerical Analysis

SN - 0764-583X

IS - 2

ER -