2,108 research outputs found
Variational Multiscale Stabilization and the Exponential Decay of Fine-scale Correctors
This paper addresses the variational multiscale stabilization of standard
finite element methods for linear partial differential equations that exhibit
multiscale features. The stabilization is of Petrov-Galerkin type with a
standard finite element trial space and a problem-dependent test space based on
pre-computed fine-scale correctors. The exponential decay of these correctors
and their localisation to local cell problems is rigorously justified. The
stabilization eliminates scale-dependent pre-asymptotic effects as they appear
for standard finite element discretizations of highly oscillatory problems,
e.g., the poor approximation in homogenization problems or the pollution
effect in high-frequency acoustic scattering
On stability of discretizations of the Helmholtz equation (extended version)
We review the stability properties of several discretizations of the
Helmholtz equation at large wavenumbers. For a model problem in a polygon, a
complete -explicit stability (including -explicit stability of the
continuous problem) and convergence theory for high order finite element
methods is developed. In particular, quasi-optimality is shown for a fixed
number of degrees of freedom per wavelength if the mesh size and the
approximation order are selected such that is sufficiently small and
, and, additionally, appropriate mesh refinement is used near
the vertices. We also review the stability properties of two classes of
numerical schemes that use piecewise solutions of the homogeneous Helmholtz
equation, namely, Least Squares methods and Discontinuous Galerkin (DG)
methods. The latter includes the Ultra Weak Variational Formulation
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
A multiscale method for heterogeneous bulk-surface coupling
In this paper, we construct and analyze a multiscale (finite element) method
for parabolic problems with heterogeneous dynamic boundary conditions. As
origin, we consider a reformulation of the system in order to decouple the
discretization of bulk and surface dynamics. This allows us to combine
multiscale methods on the boundary with standard Lagrangian schemes in the
interior. We prove convergence and quantify explicit rates for low-regularity
solutions, independent of the oscillatory behavior of the heterogeneities. As a
result, coarse discretization parameters, which do not resolve the fine scales,
can be considered. The theoretical findings are justified by a number of
numerical experiments including dynamic boundary conditions with random
diffusion coefficients
Stabilised finite element methods for ill-posed problems with conditional stability
In this paper we discuss the adjoint stabilised finite element method
introduced in, E. Burman, Stabilized finite element methods for nonsymmetric,
noncoercive and ill-posed problems. Part I: elliptic equations, SIAM Journal on
Scientific Computing, and how it may be used for the computation of solutions
to problems for which the standard stability theory given by the Lax-Milgram
Lemma or the Babuska-Brezzi Theorem fails. We pay particular attention to
ill-posed problems that have some conditional stability property and prove
(conditional) error estimates in an abstract framework. As a model problem we
consider the elliptic Cauchy problem and provide a complete numerical analysis
for this case. Some numerical examples are given to illustrate the theory.Comment: Accepted in the proceedings from the EPSRC Durham Symposium Building
Bridges: Connections and Challenges in Modern Approaches to Numerical Partial
Differential Equation
An asymptotic preserving method for linear systems of balance laws based on Galerkin's method
We apply the concept of Asymptotic Preserving (AP) schemes to the linearized
p-system and discretize the resulting elliptic equation using standard
continuous Finite Elements instead of Finite Differences. The fully discrete
method is analyzed with respect to consistency, and we compare it numerically
with more traditional methods such as Implicit Euler's method. Numerical
results indicate that the AP method is indeed superior to more traditional
methods.Comment: Journal of Scientific Computing, 201
- …