Integration of the EPDiff equation by particle methods

The purpose of this paper is to apply particle methods to the numerical solution of the EPDiff equation. The weak solutions of EPDiff are contact discontinuities that carry momentum so that wavefront interactions represent collisions in which momentum is exchanged. This behavior allows for the description of many rich physical applications, but also introduces difficult numerical challenges. We present a particle method for the EPDiff equation that is well-suited for this class of solutions and for simulating collisions between wavefronts. Discretization by means of the particle method is shown to preserve the basic Hamiltonian, the weak and variational structure of the original problem, and to respect the conservation laws associated with symmetry under the Euclidean group. Numerical results illustrate that the particle method has superior features in both one and two dimensions, and can also be effectively implemented when the initial data of interest lies on a submanifold

Convergence of a linearly transformed particle method for aggregation equations

We study a linearly transformed particle method for the aggregation equation with smooth or singular interaction forces. For the smooth interaction forces, we provide convergence estimates in $L^1$ and $L^\infty$ norms depending on the regularity of the initial data. Moreover, we give convergence estimates in bounded Lipschitz distance for measure valued solutions. For singular interaction forces, we establish the convergence of the error between the approximated and exact flows up to the existence time of the solutions in $L^1 \cap L^p$ norm

The instanton method and its numerical implementation in fluid mechanics

A precise characterization of structures occurring in turbulent fluid flows at high Reynolds numbers is one of the last open problems of classical physics. In this review we discuss recent developments related to the application of instanton methods to turbulence. Instantons are saddle point configurations of the underlying path integrals. They are equivalent to minimizers of the related Freidlin-Wentzell action and known to be able to characterize rare events in such systems. While there is an impressive body of work concerning their analytical description, this review focuses on the question on how to compute these minimizers numerically. In a short introduction we present the relevant mathematical and physical background before we discuss the stochastic Burgers equation in detail. We present algorithms to compute instantons numerically by an efficient solution of the corresponding Euler-Lagrange equations. A second focus is the discussion of a recently developed numerical filtering technique that allows to extract instantons from direct numerical simulations. In the following we present modifications of the algorithms to make them efficient when applied to two- or three-dimensional fluid dynamical problems. We illustrate these ideas using the two-dimensional Burgers equation and the three-dimensional Navier-Stokes equations

Uniform convergence of a linearly transformed particle method for the Vlasov-Poisson system

International audienceA particle method with linear transformation of the particle shape functions is studied for the 1d-1v Vlasov-Poisson equation, and a priori error estimates are proven which show that the approximated densities converge in the uniform norm. When compared to standard fixed-shape particle methods, the present approach can be seen as a way to gain one order in the convergence rate of the particle trajectories at the cost of linearly transforming each particle shape. It also allows to compute strongly convergent densities with particles that overlap in a bounded way