4,231 research outputs found
Energy-corrected FEM and explicit time-stepping for parabolic problems
The presence of corners in the computational domain, in general, reduces the
regularity of solutions of parabolic problems and diminishes the convergence
properties of the finite element approximation introducing a so-called
"pollution effect". Standard remedies based on mesh refinement around the
singular corner result in very restrictive stability requirements on the
time-step size when explicit time integration is applied. In this article, we
introduce and analyse the energy-corrected finite element method for parabolic
problems, which works on quasi-uniform meshes, and, based on it, create fast
explicit time discretisation. We illustrate these results with extensive
numerical investigations not only confirming the theoretical results but also
showing the flexibility of the method, which can be applied in the presence of
multiple singular corners and a three-dimensional setting. We also propose a
fast explicit time-stepping scheme based on a piecewise cubic energy-corrected
discretisation in space completed with mass-lumping techniques and numerically
verify its efficiency
A-stable Runge-Kutta methods for semilinear evolution equations
We consider semilinear evolution equations for which the linear part
generates a strongly continuous semigroup and the nonlinear part is
sufficiently smooth on a scale of Hilbert spaces. In this setting, we prove the
existence of solutions which are temporally smooth in the norm of the lowest
rung of the scale for an open set of initial data on the highest rung of the
scale. Under the same assumptions, we prove that a class of implicit,
-stable Runge--Kutta semidiscretizations in time of such equations are
smooth as maps from open subsets of the highest rung into the lowest rung of
the scale. Under the additional assumption that the linear part of the
evolution equation is normal or sectorial, we prove full order convergence of
the semidiscretization in time for initial data on open sets. Our results
apply, in particular, to the semilinear wave equation and to the nonlinear
Schr\"odinger equation
Numerical Methods for the Fractional Laplacian: a Finite Difference-quadrature Approach
The fractional Laplacian is a non-local operator which
depends on the parameter and recovers the usual Laplacian as . A numerical method for the fractional Laplacian is proposed, based on
the singular integral representation for the operator. The method combines
finite difference with numerical quadrature, to obtain a discrete convolution
operator with positive weights. The accuracy of the method is shown to be
. Convergence of the method is proven. The treatment of far
field boundary conditions using an asymptotic approximation to the integral is
used to obtain an accurate method. Numerical experiments on known exact
solutions validate the predicted convergence rates. Computational examples
include exponentially and algebraically decaying solution with varying
regularity. The generalization to nonlinear equations involving the operator is
discussed: the obstacle problem for the fractional Laplacian is computed.Comment: 29 pages, 9 figure
Numerical analysis of nonlinear subdiffusion equations
We present a general framework for the rigorous numerical analysis of
time-fractional nonlinear parabolic partial differential equations, with a
fractional derivative of order in time. The framework relies
on three technical tools: a fractional version of the discrete Gr\"onwall-type
inequality, discrete maximal regularity, and regularity theory of nonlinear
equations. We establish a general criterion for showing the fractional discrete
Gr\"onwall inequality, and verify it for the L1 scheme and convolution
quadrature generated by BDFs. Further, we provide a complete solution theory,
e.g., existence, uniqueness and regularity, for a time-fractional diffusion
equation with a Lipschitz nonlinear source term. Together with the known
results of discrete maximal regularity, we derive pointwise norm
error estimates for semidiscrete Galerkin finite element solutions and fully
discrete solutions, which are of order (up to a logarithmic factor)
and , respectively, without any extra regularity assumption on
the solution or compatibility condition on the problem data. The sharpness of
the convergence rates is supported by the numerical experiments
Regularity of mean curvature flow of graphs on Lie groups free up to step 2
We consider (smooth) solutions of the mean curvature flow of graphs over
bounded domains in a Lie group free up to step two (and not necessarily
nilpotent), endowed with a one parameter family of Riemannian metrics
\sigma_\e collapsing to a subRiemannian metric as \e\to 0. We
establish estimates for this flow, that are uniform as \e\to 0
and as a consequence prove long time existence for the subRiemannian mean
curvature flow of the graph. Our proof extend to the setting of every step two
Carnot group (not necessarily free) and can be adapted following our previous
work in \cite{CCM3} to the total variation flow.Comment: arXiv admin note: text overlap with arXiv:1212.666
- …