Affine particle in cell method for MAC grids and fluid simulation
| Autor: | Ounan Ding, Craig Schroeder, Tamar Shinar |
|---|---|
| Rok vydání: | 2020 |
| Předmět: |
Angular momentum
Physics and Astronomy (miscellaneous) FLIP MPM Basis function 010103 numerical & computational mathematics Computational fluid dynamics 01 natural sciences Hybrid Lagrangian/Eulerian Mathematical Sciences Quantitative Biology::Cell Behavior Momentum Engineering 0101 mathematics Physics Numerical Analysis APIC business.industry PIC Applied Mathematics Mathematical analysis Dissipation Computer Science Applications 010101 applied mathematics Computational Mathematics Modeling and Simulation Physical Sciences Dissipative system Affine transformation business Interpolation |
| Popis: | We present a new technique for transferring momentum and velocity between particles and MAC grids based on the Affine-Particle-In-Cell (APIC) framework [1] , [2] previously developed for co-located grids. APIC represents particle velocities as locally affine, rather than locally constant as in traditional PIC. These previous APIC schemes were designed primarily as an improvement on Particle-in-Cell (PIC) transfers, which tend to be heavily dissipative, and as an alternative to Fluid Implicit Particle (FLIP) transfers, which tend to be noisy. The original APIC paper [1] proposed APIC-style transfers for MAC grids, based on a limit for multilinear interpolation. We extend these to the case of smooth basis functions and show that the proposed transfers satisfy all of the original APIC properties. In particular, we achieve conservation of angular momentum across our transfers, something which eluded [1] . Early indications in [1] suggested that APIC might be suitable for simulating high Reynolds fluids due to favorable retention of vortices, but these properties were not studied further and the question was left unresolved. One significant drawback of APIC relative to FLIP is that energy is dissipated even when Δ t = 0 . We use two dimensional Fourier analysis to investigate dissipation in this important limit. We investigate dissipation and vortex retention numerically in the general case to quantify the effectiveness of APIC compared with PIC, FLIP, and XPIC. |
| Databáze: | OpenAIRE |
| Externí odkaz: |
|