21,924 research outputs found
On round-off error for adaptive finite element methods
Round-off error analysis has been historically studied by analyzing the condition number of the associated matrix. By controlling the size of the condition number, it is possible to guarantee a prescribed round-off error tolerance. However, the opposite is not true, since it is possible to have a system of linear equations with an arbitrarily large condition number that still delivers a small round-off error. In this paper, we perform a round-off error analysis in context of 1D and 2D hp-adaptive Finite Element simulations for the case of Poisson equation. We conclude that boundary conditions play a fundamental role on the round-off error analysis, specially for the so-called 'radical meshes'. Moreover, we illustrate the importance of the right-hand side when analyzing the round-off error, which is independent of the condition number of the matrix. © 2012 Published by Elsevier Ltd
Geostrophic balance preserving interpolation in mesh adaptive shallow-water ocean modelling
The accurate representation of geostrophic balance is an essential
requirement for numerical modelling of geophysical flows. Significant effort is
often put into the selection of accurate or optimal balance representation by
the discretisation of the fundamental equations. The issue of accurate balance
representation is particularly challenging when applying dynamic mesh
adaptivity, where there is potential for additional imbalance injection when
interpolating to new, optimised meshes.
In the context of shallow-water modelling, we present a new method for
preservation of geostrophic balance when applying dynamic mesh adaptivity. This
approach is based upon interpolation of the Helmholtz decomposition of the
Coriolis acceleration. We apply this in combination with a discretisation for
which states in geostrophic balance are exactly steady solutions of the
linearised equations on an f-plane; this method guarantees that a balanced and
steady flow on a donor mesh remains balanced and steady after interpolation
onto an arbitrary target mesh, to within machine precision. We further
demonstrate the utility of this interpolant for states close to geostrophic
balance, and show that it prevents pollution of the resulting solutions by
imbalanced perturbations introduced by the interpolation
A cloudy Vlasov solution
We propose to integrate the Vlasov-Poisson equations giving the evolution of
a dynamical system in phase-space using a continuous set of local basis
functions. In practice, the method decomposes the density in phase-space into
small smooth units having compact support. We call these small units ``clouds''
and choose them to be Gaussians of elliptical support. Fortunately, the
evolution of these clouds in the local potential has an analytical solution,
that can be used to evolve the whole system during a significant fraction of
dynamical time. In the process, the clouds, initially round, change shape and
get elongated. At some point, the system needs to be remapped on round clouds
once again. This remapping can be performed optimally using a small number of
Lucy iterations. The remapped solution can be evolved again with the cloud
method, and the process can be iterated a large number of times without showing
significant diffusion. Our numerical experiments show that it is possible to
follow the 2 dimensional phase space distribution during a large number of
dynamical times with excellent accuracy. The main limitation to this accuracy
is the finite size of the clouds, which results in coarse graining the
structures smaller than the clouds and induces small aliasing effects at these
scales. However, it is shown in this paper that this method is consistent with
an adaptive refinement algorithm which allows one to track the evolution of the
finer structure in phase space. It is also shown that the generalization of the
cloud method to the 4 dimensional and the 6 dimensional phase space is quite
natural.Comment: 46 pages, 25 figures, submitted to MNRA
- …