504 research outputs found
Numerical methods for time-fractional evolution equations with nonsmooth data: a concise overview
Over the past few decades, there has been substantial interest in evolution
equations that involving a fractional-order derivative of order
in time, due to their many successful applications in
engineering, physics, biology and finance. Thus, it is of paramount importance
to develop and to analyze efficient and accurate numerical methods for reliably
simulating such models, and the literature on the topic is vast and fast
growing. The present paper gives a concise overview on numerical schemes for
the subdiffusion model with nonsmooth problem data, which are important for the
numerical analysis of many problems arising in optimal control, inverse
problems and stochastic analysis. We focus on the following aspects of the
subdiffusion model: regularity theory, Galerkin finite element discretization
in space, time-stepping schemes (including convolution quadrature and L1 type
schemes), and space-time variational formulations, and compare the results with
that for standard parabolic problems. Further, these aspects are showcased with
illustrative numerical experiments and complemented with perspectives and
pointers to relevant literature.Comment: 24 pages, 3 figure
Recommended from our members
Space-time Methods for Time-dependent Partial Differential Equations
Modern discretizations of time-dependent PDEs consider the full problem in the space-time cylinder and aim to overcome limitations of classical approaches such as the method of lines (first discretize in space and then solve the resulting ODE) and the Rothe method (first discretize in time and then solve the PDE). A main advantage of a holistic space-time method is the direct access to space-time adaptivity and to the backward problem (required for the dual problem in optimization or error control). Moreover, this allows for parallel solution strategies simultaneously in time and space.
Several space-time concepts where proposed (different conforming and nonconforming space-time finite elements, the parareal method, wavefront relaxation etc.) but this topic has become a rapidly growing field in numerical analysis and scientific computing. In this workshop the focus is the development of adaptive and flexible space-time discretization methods for solving parabolic and hyperbolic space-time partial differential equations
An Analysis of Galerkin Proper Orthogonal Decomposition for Subdiffusion
In this work, we develop a novel Galerkin-L1-POD scheme for the subdiffusion
model with a Caputo fractional derivative of order in time,
which is often used to describe anomalous diffusion processes in heterogeneous
media. The nonlocality of the fractional derivative requires storing all the
solutions from time zero. The proposed scheme is based on continuous piecewise
linear finite elements, L1 time stepping, and proper orthogonal decomposition
(POD). By constructing an effective reduced-order scheme using problem-adapted
basis functions, it can significantly reduce the computational complexity and
storage requirement. We shall provide a complete error analysis of the scheme
under realistic regularity assumptions by means of a novel energy argument.
Extensive numerical experiments are presented to verify the convergence
analysis and the efficiency of the proposed scheme.Comment: 25 pp, 5 figure
Discontinuous Galerkin for the wave equation: a simplified a priori error analysis
Standard discontinuous Galerkin methods, based on piecewise polynomials of
degree \qq\geq 0, are considered for temporal semi-discretization for second
order hyperbolic equations. The main goal of this paper is to present a simple
and straight forward a priori error analysis of optimal order with minimal
regularity requirement on the solution. Uniform norm in time error estimates
are also proved for the constant and linear cases. To this end, energy
identities and stability estimates of the discrete problem are proved for a
slightly more general problem. These are used to prove optimal order a priori
error estimates with minimal regularity requirement on the solution. The
combination with the classic continuous Galerkin finite element discretization
in space variable is used, to formulate a full-discrete scheme. The a priori
error analysis is presented. Numerical experiments are performed to verify the
theoretical rate of convergence.Comment: 23 pages, 1 figur
Recommended from our members
Mini-Workshop: Numerical Analysis for Non-Smooth PDE-Constrained Optimal Control Problems
This mini-workshop brought together leading experts working on various aspects of numerical analysis for optimal control problems with nonsmoothness. Fifteen extended abstracts summarize the presentations at this mini-workshop
Error Estimates for Approximations of Distributed Order Time Fractional Diffusion with Nonsmooth Data
In this work, we consider the numerical solution of an initial boundary value
problem for the distributed order time fractional diffusion equation. The model
arises in the mathematical modeling of ultra-slow diffusion processes observed
in some physical problems, whose solution decays only logarithmically as the
time tends to infinity. We develop a space semidiscrete scheme based on the
standard Galerkin finite element method, and establish error estimates optimal
with respect to data regularity in and norms for both smooth
and nonsmooth initial data. Further, we propose two fully discrete schemes,
based on the Laplace transform and convolution quadrature generated by the
backward Euler method, respectively, and provide optimal convergence rates in
the norm, which exhibits exponential convergence and first-order
convergence in time, respectively. Extensive numerical experiments are provided
to verify the error estimates for both smooth and nonsmooth initial data, and
to examine the asymptotic behavior of the solution.Comment: 25 pages, 2 figure
- …