140 research outputs found
Singular solutions, momentum maps and computational anatomy
This paper describes the variational formulation of template matching
problems of computational anatomy (CA); introduces the EPDiff evolution
equation in the context of an analogy between CA and fluid dynamics; discusses
the singular solutions for the EPDiff equation and explains why these singular
solutions exist (singular momentum map). Then it draws the consequences of
EPDiff for outline matching problem in CA and gives numerical examples
Recommended from our members
Computational modes and grid imprinting on five quasi-uniform spherical C-grids
Currently, most operational forecasting models use latitude-longitude grids, whose convergence of meridians towards the poles limits parallel scaling. Quasi-uniform grids might avoid this limitation. Thuburn et al, JCP, 2009 and Ringler et al, JCP, 2010 have developed a method for arbitrarily-structured, orthogonal C-grids (TRiSK), which has many of the desirable properties of the C-grid on latitude-longitude grids but which works on a variety of quasi-uniform grids. Here, five quasi-uniform, orthogonal grids of the sphere are investigated using TRiSK to solve the shallow-water equations.
We demonstrate some of the advantages and disadvantages of the hexagonal and triangular icosahedra, a Voronoi-ised cubed sphere, a Voronoi-ised skipped latitude-longitude grid and a grid of kites in comparison to a full latitude-longitude grid. We will show that the hexagonal-icosahedron gives the most accurate results (for least computational cost). All of the grids suffer from spurious computational modes; this is especially true of the kite grid,
despite it having exactly twice as many velocity degrees of freedom as height degrees of freedom. However, the computational modes are easiest to control on the hexagonal icosahedron since they consist of vorticity oscillations on the dual grid which can be controlled using a diffusive advection scheme for potential vorticity
Compatible finite element methods for geophysical fluid dynamics
This article surveys research on the application of compatible finite element
methods to large scale atmosphere and ocean simulation. Compatible finite
element methods extend Arakawa's C-grid finite difference scheme to the finite
element world. They are constructed from a discrete de Rham complex, which is a
sequence of finite element spaces which are linked by the operators of
differential calculus. The use of discrete de Rham complexes to solve partial
differential equations is well established, but in this article we focus on the
specifics of dynamical cores for simulating weather, oceans and climate. The
most important consequence of the discrete de Rham complex is the
Hodge-Helmholtz decomposition, which has been used to exclude the possibility
of several types of spurious oscillations from linear equations of geophysical
flow. This means that compatible finite element spaces provide a useful
framework for building dynamical cores. In this article we introduce the main
concepts of compatible finite element spaces, and discuss their wave
propagation properties. We survey some methods for discretising the transport
terms that arise in dynamical core equation systems, and provide some example
discretisations, briefly discussing their iterative solution. Then we focus on
the recent use of compatible finite element spaces in designing structure
preserving methods, surveying variational discretisations, Poisson bracket
discretisations, and consistent vorticity transport.Comment: correction of some typo
- …