2016/10/17 by Daniel Duque, Duque, Daniel, Pep Español +1 · 1 citation
Engineering · Mathematics · Physics and Astronomy · #Advection #Applied mathematics #Classical mechanics #Computational Physics (physics.comp-ph) #Computer science #Context (archaeology) #Convergence (economics) #Eulerian path #FOS: Physical sciences #Finite element method #Fluid Dynamics (physics.flu-dyn) #Fluid Dynamics Simulations and Interactions #Fluid Dynamics and Vibration Analysis #Interpolation (computer graphics) #Lagrangian #Lattice Boltzmann Simulation Studies #Mathematics #Mechanics #Physics #physics.comp-ph #physics.flu-dyn
paper · pdf · doi:10.48550/arxiv.1610.05266
published in arXiv (Cornell University) (Cornell University) · Draft article
openalex publication_date 2016/10/17 · arxiv created 2017/01/26 · arxiv updated 2017/01/27 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
The advantage of particle Lagrangian methods in computational fluid dynamics is that advection is accurately modeled. However, this complicates the calculation of space derivatives. If a mesh is employed, it must be updated at each time step. On the other hand, fixed mesh, Eulerian, formulations benefit from the mesh being defined at the beginning of the simulation, but feature non-linear advection terms. It therefore seems natural to combine the two approaches, using a fixed mesh to perform calculations related to space derivatives, and using the particles to advect the information with time. The idea of combining Lagrangian particles and a fixed mesh goes back to Particle-in-Cell methods, and is here considered within the context of the finite element method (FEM) for the fixed mesh, and the particle FEM (pFEM) for the particles. Our results, in agreement with recent works, show that interpolation ("projection") errors, especially from particles to mesh, are the culprits of slow convergence of the method if standard, linear, FEs, are employed. By including quadratically consistent shape functions for the particles, the rate of convergence is restored to values competitive with purely Lagrangian simulations. The procedure is validated on the Zalesak's disk problem, the Taylor-Green vortex sheet flow, and the Rayleigh-Taylor instability.