Abstract
Hybrid Lagrangian/Eulerian simulation is commonplace in computer graphics for fluids and other materials undergoing large deformation. In these methods, particles are used to resolve transport and topological change, while a background Eulerian grid is used for computing mechanical forces and collision responses. Particle-in-Cell (PIC) techniques, particularly the Fluid Implicit Particle (FLIP) variants have become the norm in computer graphics calculations. While these approaches have proven very powerful, they do suffer from some well known limitations. The original PIC is stable, but highly dissipative, while FLIP, designed to remove this dissipation, is more noisy and at times, unstable. We present a novel technique designed to retain the stability of the original PIC, without suffering from the noise and instability of FLIP. Our primary observation is that the dissipation in the original PIC results from a loss of information when transferring between grid and particle representations. We prevent this loss of information by augmenting each particle with a locally affine, rather than locally constant, description of the velocity. We show that this not only stably removes the dissipation of PIC, but that it also allows for exact conservation of angular momentum across the transfers between particles and grid.
Supplemental Material
Available for Download
Supplemental files
- Ando, R., and Tsuruno, R. 2011. A particle-based method for preserving fluid sheets. In Proc ACM SIGGRAPH/Eurographics Symp Comp Anim, SCA '11, 7--16. Google ScholarDigital Library
- Ando, R., Thurey, N., and Tsuruno, R. 2012. Preserving fluid sheets with adaptively sampled anisotropic particles. IEEE Trans Vis Comp Graph 18, 8, 1202--1214. Google ScholarDigital Library
- Ando, R., Thurey, N., and Wojtan, C. 2013. Highly adaptive liquid simulations on tetrahedral meshes. ACM Trans Graph 32, 4, 103:1--103:10. Google ScholarDigital Library
- Bargteil, A., Wojtan, C., Hodgins, J., and Turk, G. 2007. A finite element method for animating large viscoplastic flow. ACM Trans Graph 26, 3. Google ScholarDigital Library
- Batty, C., and Bridson, R. 2008. Accurate viscous free surfaces for buckling, coiling, and rotating liquids. Proc ACM SIGGRAPH/ Eurograph Symp Comp Anim, 219--228. Google ScholarDigital Library
- Batty, C., Bertails, F., and Bridson, R. 2007. A fast variational framework for accurate solid-fluid coupling. ACM Trans Graph 26, 3. Google ScholarDigital Library
- Boyd, L., and Bridson, R. 2012. Multiflip for energetic two-phase fluid simulation. ACM Trans Graph 31, 2, 16:1--16:12. Google ScholarDigital Library
- Brackbill, J., and Ruppel, H. 1986. Flip: A method for adaptively zoned, particle-in-cell calculations of fluid flows in two dimensions. J Comp Phys 65, 314--343. Google ScholarDigital Library
- Brackbill, J., Kothe, D., and Ruppel, H. 1988. Flip: A low-dissipation, pic method for fluid flow. Comp Phys Comm 48, 25--38.Google ScholarCross Ref
- Brackbill, J. 1988. The ringing instability in particle-in-cell calculations of low-speed flow. J Comp Phys 75, 2, 469--492. Google ScholarDigital Library
- Bridson, R. 2008. Fluid simulation for computer graphics. Taylor & Francis. Google ScholarDigital Library
- Chentanez, N., and Muller, M. 2010. Real-time simulation of large bodies of water with small scale details. In Proc ACM SIGGRAPH/Eurograph Symp Comp Anim, SCA '10, 197--206. Google ScholarDigital Library
- Chentanez, N., and Muller, M. 2011. Real-time eulerian water simulation using a restricted tall cell grid. ACM Trans Graph 30, 4, 82:1--82:10. Google ScholarDigital Library
- Chentanez, N., and Muller, M. 2014. Coupling 3d eulerian, height field and particle methods for the simulation of large scale liquid phenomena. In Proc ACM SIGGRAPH/Eurograph Symp Comp Anim, SCA '14.Google Scholar
- Cornelis, J., Ihmsen, M., Peer, A., and Teschner, M. 2014. Iisph-flip for incompressible fluids. Comp Graph Forum 33, 2, 255--262. Google ScholarDigital Library
- Edwards, E., and Bridson, R. 2012. A high-order accurate particle-in-cell method. Int J Numer Meth Eng 90, 1073--1088.Google ScholarCross Ref
- Edwards, E., and Bridson, R. 2014. Detailed water with coarse grids: combining surface meshes and adaptive discontinuous galerkin. ACM Trans Graph 33, 4, 136:1--136:9. Google ScholarDigital Library
- Enright, D., Marschner, S., and Fedkiw, R. 2002. Animation and rendering of complex water surfaces. ACM Trans Graph 21, 3, 736--744. Google ScholarDigital Library
- Feldman, B., O'Brien, J., and Arikan, O. 2003. Animating suspended particle explosions. SIGGRAPH '03 22, 3, 708--715. Google ScholarDigital Library
- Foster, N., and Metaxas, D. 1996. Realistic animation of liquids. Graph Mod Imag Proc 58, 471--483. Google ScholarDigital Library
- Gao, Y., li, C., Hu, S., and Barsky, B. 2009. Simulating gaseous fluids with low and high speeds. Comp Graph Forum 28, 28, 1845--1852.Google ScholarCross Ref
- Gerszewski, D., and Bargteil, A. 2013. Physics-based animation of large-scale splashing liquids. ACM Trans Graph 32, 6, 185:1--185:6. Google ScholarDigital Library
- Harlow, F., and Welch, E. 1965. Numerical calculation of time dependent viscous flow of fluid with a free surface. Phys Fluid 8, 12, 2182--2189.Google ScholarCross Ref
- Harlow, F. 1964. The particle-in-cell method for numerical solution of problems in fluid dynamics. Meth Comp Phys 3, 319--343.Google Scholar
- Hong, J., Lee, H., Yoon, J., and Kim, C. 2008. Bubbles alive. ACM Trans Graph 27, 3, 48:1--48:4. Google ScholarDigital Library
- Hong, W., House, D., and Keyser, J. 2008. Adaptive particles for incompressible fluid simulation. Vis Comp 24, 7, 535--543. Google ScholarDigital Library
- Hong, W., House, D., and Keyser, J. 2009. An adaptive sampling approach to incompressible particle-based fluid. Theory Pract Comp Graph, 69--76.Google Scholar
- Ihmsen, M., Cornelis, J., Solenthaler, B., Horvath, C., and Teschner, M. 2013. Implicit incompressible sph. IEEE Trans Vis Comp Graph 20, 3, 426--435. Google ScholarDigital Library
- Kim, J., Cha, D., Chang, B., Koo, B., and Ihm, I. 2006. Practical animation of turbulent splashing water. In Proc ACM SIGGRAPH/Eurograph Symp Comp Anim, SCA '06, 335--344. Google ScholarDigital Library
- Lee, H., Hong, J., and Kim, C. 2009. Interchangeable sph and level set method in multiphase fluids. Vis Comp 25, 5, 713--718. Google ScholarDigital Library
- Losasso, F., Talton, J., Kwatra, N., and Fedkiw, R. 2008. Two-way coupled sph and particle level set fluid simulation. IEEE Trans Vis Comp Graph 14, 797--804. Google ScholarDigital Library
- Love, E., and Sulsky, D. 2006. An unconditionally stable, energy-momentum consistent implementation of the the material point method. Comp Meth App Mech Eng 195, 3903--3925.Google ScholarCross Ref
- Mihalef, V., Metaxas, D., and Sussman, M. 2007. Textured liquids based on the marker level set. Comp Graph Forum, 457--466.Google Scholar
- Muller, K., Fedosov, D., and Gompper, G. 2015. Smoothed dissipative particle dynamics with angular momentum conservation. J Comp Phys 281, 301--315.Google ScholarDigital Library
- Narain, R., Golas, A., and Lin, M. 2013. Free-flowing granular materials with two-way solid coupling. ACM Trans Graph 29, 6, 173:1--173:10. Google ScholarDigital Library
- Patkar, S., Aanjaneya, M., Karpman, D., and Fedkiw, R. 2013. A hybrid lagrangian-eulerian formulation for bubble generation and dynamics. In Proc ACM SIGGRAPH/Eurograp Symp Comp Anim, SCA '13, 105--114. Google ScholarDigital Library
- Raveendran, K., Wojtan, C., and Turk, G. 2011. Hybrid sph. In Proc 2011 ACM SIGGRAPH/Eurograp Symp Comp Anim, SCA '11, 33--42. Google ScholarDigital Library
- Sifakis, E., Shinar, T., Irving, G., and Fedkiw, R. 2007. Hybrid simulation of deformable solids. In Proc ACM SIGGRAPH/Eurograph Symp Comp Anim, 81--90. Google ScholarDigital Library
- Sin, F., Bargteil, A., and Hodgins, J. 2009. A point-based method for animating incompressible flow. In Proc ACM SIGGRAPH/Eurograph Symp Comp Anim, 247--255. Google ScholarDigital Library
- Song, O., Kim, D., and Ko, H. 2009. Derivative particles for simulating detailed movements of fluids. IEEE Trans Vis Comp Graph, 247--255. Google ScholarDigital Library
- Stomakhin, A., Schroeder, C., Chai, L., Teran, J., and Selle, A. 2013. A material point method for snow simulation. ACM Trans Graph 32, 4, 102:1--102:10. Google ScholarDigital Library
- Stomakhin, A., Schroeder, C., Jiang, C., Chai, L., Teran, J., and Selle, A. 2014. Augmented mpm for phasechange and varied materials. ACM Trans Graph 33, 4, 138:1--138:11. Google ScholarDigital Library
- Sulsky, D., Zhou, S., and Schreyer, H. 1995. Application of a pic method to solid mechanics. Comp Phys Comm 87, 1, 236--252.Google ScholarCross Ref
- Um, K., Baek, S., and Han, J. 2014. Advanced hybrid particle-grid method with sub-grid particle correction. Comp Graph Forum 33, 209--218. Google ScholarDigital Library
- Yabe, T., Xiao, F., and Utsumi, T. 2001. The constrained interpolation profile method for multiphase analysis. J Comp Phys 169, 556--593. Google ScholarDigital Library
- Zhu, Y., and Bridson, R. 2005. Animating sand as a fluid. ACM Trans Graph 24, 3, 965--972. Google ScholarDigital Library
- Zhu, B., Yang, X., and Fan, Y. 2010. Creating and preserving vortical details in sph fluid. Comp Graph Forum 29, 7, 2207--2214.Google ScholarCross Ref
Index Terms
- The affine particle-in-cell method
Recommendations
An angular momentum conserving affine-particle-in-cell method
We present a new technique for transferring momentum and velocity between particles and grid with Particle-In-Cell (PIC) [1] calculations which we call Affine-Particle-In-Cell (APIC). APIC represents particle velocities as locally affine, rather than ...
A polynomial particle-in-cell method
Recently the Affine Particle-In-Cell (APIC) Method was proposed by Jiang et al.[2015; 2017b] to improve the accuracy of the transfers in Particle-In-Cell (PIC) [Harlow 1964] techniques by augmenting each particle with a locally affine, rather than ...
IQ-MPM: an interface quadrature material point method for non-sticky strongly two-way coupled nonlinear solids and fluids
We propose a novel scheme for simulating two-way coupled interactions between nonlinear elastic solids and incompressible fluids. The key ingredient of this approach is a ghost matrix operator-splitting scheme for strongly coupled nonlinear elastica and ...
Comments