“Resolving fluid boundary layers with particle strength exchange and weak adaptivity”

  • ©Xinxin Zhang, Minchen Li, and Robert Bridson




    Resolving fluid boundary layers with particle strength exchange and weak adaptivity





    Most fluid scenarios in graphics have a high Reynolds number, where viscosity is dominated by inertial effects, thus most solvers drop viscosity altogether: numerical damping from coarse grids is generally stronger than physical viscosity while resembling it in character. However, viscosity remains crucial near solid boundaries, in the boundary layer, to a large extent determining the look of the flow as a function of Reynolds number. Typical graphics simulations do not resolve boundary layer dynamics, so their look is determined mostly by numerical errors with the given grid size and time step, rather than physical parameters. We introduce two complementary techniques to capture boundary layer dynamics, bringing more physical control and predictability. We extend the FLIP particle-grid method with viscous particle strength exchange[Rivoalen and Huberson 2001] to better transfer momentum at solid boundaries, dubbed VFLIP. We also introduce Weakly Higher Resolution Regional Projection (WHIRP), a cheap and simple way to increase grid resolution where important by overlaying high resolution grids on the global coarse grid.


    1. Ando, R., Thürey, N., and Wojtan, C. 2013. Highly adaptive liquid simulations on tetrahedral meshes. ACM Trans. Graph. (Proc. SIGGRAPH 2013) (July). Google ScholarDigital Library
    2. Ando, R., Thürey, N., and Wojtan, C. 2015. A dimension-reduced pressure solver for liquid simulations. EUROGRAPHICS 2015. Google ScholarDigital Library
    3. Batty, C., Bertails, F., and Bridson, R. 2007. A fast variational framework for accurate solid-fluid coupling. ACM Trans. Graph. (Proc. SIGGRAPH) 26, 3, 100. Google ScholarDigital Library
    4. Bridson, R. 2008. Fluid Simulation for Computer Graphics. A K Peters / CRC Press. Google ScholarDigital Library
    5. Chorin, A. J. 1973. Numerical study of slightly viscous flow. Journal of Fluid Mechanics Digital Archive 57, 04, 785–796.Google ScholarCross Ref
    6. Edwards, E., and Bridson, R. 2014. Detailed water with coarse grids: combining surface meshes and adaptive discontinuous Galerkin. ACM Trans. Graph. (Proc. SIGGRAPH) 33, 4, 136:1–9. Google ScholarDigital Library
    7. English, R. E., Qiu, L., Yu, Y., and Fedkiw, R. 2013. Chimera grids for water simulation. In Proceedings of the 12th ACM SIGGRAPH/Eurographics Symposium on Computer Animation, ACM, New York, NY, USA, SCA ’13, 85–94. Google ScholarDigital Library
    8. Fedkiw, R., Stam, J., and Jensen, H. W. 2001. Visual simulation of smoke. In Proc. ACM SIGGRAPH, 15–22. Google ScholarDigital Library
    9. Ferstl, F., Westermann, R., and Dick, C. 2014. Large-scale liquid simulation on adaptive hexahedral grids. Visualization and Computer Graphics, IEEE Transactions on PP, 99, 1–1.Google Scholar
    10. Gamito, M. N., Lopes, P. F., and Gomes, M. R. 1995. Two-dimensional simulation of gaseous phenomena using vortex particles. In Computer Animation and Simulation ’95, Eurographics. Springer Vienna, 3–15.Google Scholar
    11. Golas, A., Narain, R., Sewall, J., Krajcevski, P., Dubey, P., and Lin, M. 2012. Large-scale fluid simulation using velocity-vorticity domain decomposition. ACM Transactions on Graphics (TOG) 31, 6, 148. Google ScholarDigital Library
    12. Jiang, C., Schroeder, C., Selle, A., Teran, J., and Stomakhin, A. 2015. The affine particle-in-cell method. ACM Trans. Graph. 34, 4 (July), 51:1–51:10. Google ScholarDigital Library
    13. Kim, T., Thurey, N., James, D., and Gross, M. H. 2008. Wavelet turbulence for fluid simulation. ACM Trans. Graph. (Proc. SIGGRAPH) 27, 3, 50. Google ScholarDigital Library
    14. Lentine, M., Zheng, W., and Fedkiw, R. 2010. A novel algorithm for incompressible flow using only a coarse grid projection. ACM Trans. Graph. 29, 4 (July), 114:1–114:9. Google ScholarDigital Library
    15. Lienhard, J. 1966. Synopsis of Lift, Drag, and Vortex Frequency Data for Rigid Circular Cylinders. Bulletin (Washington State University. College of Engineering. Research Division). Technical Extension Service, Washington State University.Google Scholar
    16. Losasso, F., Gibou, F., and Fedkiw, R. 2004. Simulating water and smoke with an octree data structure. ACM Trans. Graph. (Proc. SIGGRAPH) 23, 3, 457–462. Google ScholarDigital Library
    17. Mas-Gallic, S. 1995. Particle approximation of a linear convection-diffusion problem with Neumann boundary conditions. SIAM J. Numer. Anal. 32, 4 (Aug.), 1098–1125. Google ScholarDigital Library
    18. McAdams, A., Sifakis, E., and Teran, J. 2010. A parallel multigrid poisson solver for fluids simulation on large grids. In Proc. ACM SIGGRAPH / Eurographics Symp. Comp. Anim., 65–74. Google ScholarDigital Library
    19. Molemaker, J., Cohen, J. M., Patel, S., and Noh, J. 2008. Low viscosity flow simulations for animation. In SCA ’08: Proceedings of the 2007 ACM SIGGRAPH/Eurographics symposium on Computer animation, Eurographics Association, 9–18. Google ScholarDigital Library
    20. Park, S. I., and Kim, M.-J. 2005. Vortex fluid for gaseous phenomena. In Proc. ACM SIGGRAPH / Eurographics Symp. Comp. Animation, 261–270. Google ScholarDigital Library
    21. Pfaff, T., Thuerey, N., Selle, A., and Gross, M. 2009. Synthetic turbulence using artificial boundary layers. ACM Trans. Graph. 28, 5, 121. Google ScholarDigital Library
    22. Ralston, A. 1962. Runge-Kutta methods with minimum error bounds. Mathematics of Computation 16, 80, 431–437.Google ScholarCross Ref
    23. Rivoalen, E., and Huberson, S. 2001. The particle strength exchange method applied to axisymmetric viscous flows. J. Comput. Phys. 168, 2 (Apr.), 519–526. Google ScholarDigital Library
    24. Schechter, H., and Bridson, R. 2008. Evolving sub-grid turbulence for smoke animation. In Proceedings of the 2008 ACM/Eurographics Symposium on Computer Animation. Google ScholarDigital Library
    25. Setaluri, R., Aanjaneya, M., Bauer, S., and Sifakis, E. 2014. SPGrid: A sparse paged grid structure applied to adaptive smoke simulation. ACM Trans. Graph. (Proc. SIGGRAPH Asia) 33, 6, 205:1–205:12. Google ScholarDigital Library
    26. Song, O.-Y., Kim, D., and Ko, H.-S. 2007. Derivative particles for simulating detailed movements of fluids. IEEE Trans. Vis. Comp. Graph. 13, 4, 711–719. Google ScholarDigital Library
    27. Steinhoff, J., and Underhill, D. 1994. Modification of the Euler equations for “vorticity confinement”: Application to the computation of interacting vortex rings. Physics of Fluids 6, 2738–2744.Google ScholarCross Ref
    28. Stock, M. J., and Gharakhani, A. 2010. A GPU-accelerated boundary element method and vortex particle method. In AIAA 40th Fluid Dynamics Conference and Exhibit (July 2010), 1.Google ScholarCross Ref
    29. Stock, M. J., Gharakhani, A., and Stone, C. P. 2010. Modeling rotor wakes with a hybrid OVERFLOW-vortex method on a GPU cluster. In 28th AIAA Applied Aerodynamics Conference, vol. 20.Google Scholar
    30. Van Dyke, M. 1982. An album of fluid motion. Parabolic Press.Google Scholar
    31. Zhang, X., and Bridson, R. 2014. A PPPM fast summation method for fluids and beyond. ACM Trans. Graph. (Proc. SIGGRAPH Asia) 33, 6, 206:1–11. Google ScholarDigital Library
    32. Zhang, X., Bridson, R., and Greif, C. 2015. Restoring the missing vorticity in advection-projection fluid solvers. ACM Trans. Graph. 34, 4 (July), 52:1–52:8. Google ScholarDigital Library
    33. Zhu, Y., and Bridson, R. 2005. Animating sand as a fluid. ACM Trans. Graph. (Proc. SIGGRAPH) 24, 3, 965–972. Google ScholarDigital Library

ACM Digital Library Publication: