Conservative, high-order particle–mesh scheme with applications to advection-dominated flows
Jakob M. Maljaars (TU Delft - Rivers, Ports, Waterways and Dredging Engineering)
Robert Jan Labeur (TU Delft - Environmental Fluid Mechanics)
Nathaniel A. Trask (Sandia National Laboratories, New Mexico)
Deborah L. Sulsky (University of New Mexico)
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
By combining concepts from particle-in-cell (PIC) and hybridized discontinuous Galerkin (HDG) methods, we present a particle–mesh scheme for flow and transport problems which allows for diffusion-free advection while satisfying mass and momentum conservation – locally and globally – and extending to high-order spatial accuracy. This is achieved via the introduction of a novel particle–mesh projection operator which casts the particle–mesh data transfer as a PDE-constrained optimization problem, permitting advective flux functionals at cell boundaries to be inferred from particle trajectories. This optimization problem seamlessly fits in a HDG framework, whereby the control variables in the optimization problem serve as advective fluxes in the HDG scheme. The resulting algebraic problem can be solved efficiently using static condensation. The performance of the scheme is demonstrated by means of numerical examples for the linear advection–diffusion equation and the incompressible Navier–Stokes equations. The results demonstrate optimal spatial accuracy, and when combined with a θ time integration scheme, second-order temporal accuracy is shown.
help_outline