2018/06/26 by Jakob M. Maljaars, Maljaars, Jakob M., Robert Jan Labeur +5
Engineering · #Advanced Numerical Methods in Computational Mathematics #FOS: Mathematics #Fluid Dynamics Simulations and Interactions #Lattice Boltzmann Simulation Studies #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.1806.09916
openalex publication_date 2018/06/26 · openalex created_date 2018/07/10 · openalex updated_date 2026/07/28
By combining concepts from particle-in-cell (PIC) and hybridized discontinuous Galerkin (HDG) methods, we present a particle-mesh scheme which allows for diffusion-free advection, satisfies mass and momentum conservation principles in a local sense, and allows the extension to high-order spatial accuracy. To achieve this, we propose a novel particle-mesh projection operator required for the exchange of information between the particles and the mesh. Key is to cast these projections as a PDE-constrained ℓ2-optimization problem to allow the advective field naturally located on Lagrangian particles to be expressed as a mesh quantity. By expressing the control variable in terms of single-valued functions at cell interfaces, this optimization problem seamlessly fits in a HDG framework. Owing to this framework, the resulting scheme can be implemented efficiently via static condensation. The performance of the scheme is demonstrated by means of various numerical examples for the linear advection-diffusion equation and the incompressible Navier-Stokes equations. The results show that optimal spatial accuracy can be achieved, and given the particular time-stepping strategy, second-order time accuracy is confirmed. The robustness of the scheme is illustrated by considering benchmarks for advection of discontinuous fields and the Taylor-Green vortex instability in the high Reynolds number regime.