Discrete Lie Advection of Differential Forms

In this paper, we present a numerical technique for performing Lie advection of arbitrary differential forms. Leveraging advances in high-resolution finite volume methods for scalar hyperbolic conservation laws, we first discretize the interior product (also called contraction) through integrals over Eulerian approximations of extrusions. This, along with Cartan's homotopy formula and a discrete exterior derivative, can then be used to derive a discrete Lie derivative. The usefulness of this operator is demonstrated through the numerical advection of scalar fields and 1-forms on regular grids.Comment: Accepted version; to be published in J. FoC

Poisson integrators for Volterra lattice equations

The Volterra lattice equations are completely integrable and possess bi-Hamiltonian structure. They are integrated using partitioned Lobatto IIIA-B methods which preserve the Poisson structure. Modified equations are derived for the symplectic Euler and second order Lobatto IIIA-B method. Numerical results confirm preservation of the corresponding Hamiltonians, Casimirs, quadratic and cubic integrals in the long-term with different orders of accuracy.Comment: 9 pages, 2 figure

On Weak Tractability of the Clenshaw-Curtis Smolyak Algorithm

We consider the problem of integration of d-variate analytic functions defined on the unit cube with directional derivatives of all orders bounded by 1. We prove that the Clenshaw Curtis Smolyak algorithm leads to weak tractability of the problem. This seems to be the first positive tractability result for the Smolyak algorithm for a normalized and unweighted problem. The space of integrands is not a tensor product space and therefore we have to develop a different proof technique. We use the polynomial exactness of the algorithm as well as an explicit bound on the operator norm of the algorithm.Comment: 18 page

Finite element methods for surface PDEs

In this article we consider finite element methods for approximating the solution of partial differential equations on surfaces. We focus on surface finite elements on triangulated surfaces, implicit surface methods using level set descriptions of the surface, unfitted finite element methods and diffuse interface methods. In order to formulate the methods we present the necessary geometric analysis and, in the context of evolving surfaces, the necessary transport formulae. A wide variety of equations and applications are covered. Some ideas of the numerical analysis are presented along with illustrative numerical examples