DeepAI AI Chat
Log In Sign Up

LEoPart: a particle library for FEniCS

by   Jakob M. Maljaars, et al.

This paper introduces LEoPart, an add-on for the open source finite element software library FEniCS to seamlessly integrate Lagrangian particle functionality with (Eulerian) mesh-based finite element (FE) approaches. LEoPart - which is so much as to say: `Lagrangian-Eulerian on Particles' - contains tools for efficient, accurate and scalable advection of Lagrangian particles on arbitrary polyhedral meshes. In addition, LEoPart comes with several projection operators for exchanging information between the scattered particles and the mesh and vice versa. These projection operators are based on a variational framework, which allows extension to high-order accuracy. In particular, by implementing a dedicated PDE-constrained particle-mesh projection operator, LEoPart provides all the tools for diffusion-free advection, while simultaneously achieving optimal convergence and ensuring conservation of the projected particle quantities on the underlying mesh. A range of numerical examples that are prototypical to passive and active tracer methods highlight the properties and the parallel performance of the different tools in LEoPart. Finally, future developments are identified. The source code for LEoPart is actively maintained and available under an open source license at


page 24

page 25

page 30


MFEM: a modular finite element methods library

MFEM is an open-source, lightweight, flexible and scalable C++ library f...

Fully Parallel Mesh I/O using PETSc DMPlex with an Application to Waveform Modeling

Large-scale PDE simulations using high-order finite-element methods on u...

Veamy: an extensible object-oriented C++ library for the virtual element method

This paper summarizes the development of Veamy, an object-oriented C++ l...

Quality tetrahedral mesh generation with HXT

We proposed, in a recent paper (10.1002/nme.5987), a fast 3D parallel De...

Checkpoint/Restart for Lagrangian particle mesh with AMR in community code FLASH-X

In this work we present the design decisions and advantages for accompli...

Parcels v0.9: prototyping a Lagrangian Ocean Analysis framework for the petascale age

As Ocean General Circulation Models (OGCMs) move into the petascale age,...