An angular momentum conserving Affine-Particle-In-Cell method

Abstract

We present a new technique for transferring momentum and velocity between particles and grid with Particle-In-Cell (PIC) calculations which we call Affine-Particle-In-Cell (APIC). APIC represents particle velocities as locally affine, rather than locally constant as in traditional PIC. We show that this representation allows APIC to conserve linear and angular momentum across transfers while also dramatically reducing numerical diffusion usually associated with PIC. Notably, conservation is achieved with lumped mass, as opposed to the more commonly used Fluid Implicit Particle (FLIP) transfers which require a 'full' mass matrix for exact conservation. Furthermore, unlike FLIP, APIC retains a filtering property of the original PIC and thus does not accumulate velocity modes on particles as FLIP does. In particular, we demonstrate that APIC does not experience velocity instabilities that are characteristic of FLIP in a number of Material Point Method (MPM) hyperelasticity calculations. Lastly, we demonstrate that when combined with the midpoint rule for implicit update of grid momentum that linear and angular momentum are exactly conserved.

26 pages, 9 figures

Cited by in corpus (13)

• Stabilized material point methods for coupled large deformation and fluid flow in porous materials
• Circumventing volumetric locking in explicit material point methods: A simple, efficient, and general approach
• Hierarchical Optimization Time Integration for CFL-rate MPM Stepping
• BFEMP: Interpenetration-Free MPM-FEM Coupling with Barrier Contact
• Hybrid continuum-discrete simulation of granular impact dynamics
• Data-Driven Multi-Grid Solver for Accelerated Pressure Projection
• Turbulent Details Simulation for SPH Fluids via Vorticity Refinement
• Unstructured Moving Least Squares Material Point Methods: A Stable Kernel Approach With Continuous Gradient Reconstruction on General Unstructured Tessellations
• CK-MPM: A Compact-Kernel Material Point Method
• A hybrid material-point spheropolygon-element method for solid and granular material interaction
• A-ULMPM: An Arbitrary Updated Lagrangian Material Point Method for Efficient Simulation of Solids and Fluids
• A momentum preserving frictional contact algorithm based on affine particle-in-cell grid transfers
• Constrained particle-mesh projections in a hybridized discontinuous Galerkin framework with applications to advection-dominated flows