Keenan Crane
Energy-Preserving Integrators for Fluid Animation
SIGGRAPH 2009 / ACM Transactions on Graphics
Patrick Mullen Keenan Crane Dmitry Pavlov
Caltech Caltech Caltech
Yiying Tong Mathieu Desbrun
Michigan State University Caltech
Numerical viscosity has long been a problem in fluid animation. Existing methods suffer from intrinsic artificial dissipation and often apply complicated computational mechanisms to combat such effects. Consequently, dissipative behavior cannot be controlled or modeled explicitly in a manner independent of time step size, complicating the use of coarse previews and adaptive time-stepping methods. This paper proposes simple, unconditionally stable, fully Eulerian integration schemes with no numerical viscosity that are capable of maintaining the liveliness of fluid motion without recourse to corrective devices. Pressure and fluxes are solved efficiently and simultaneously in a time-reversible manner on simplicial grids, and the energy is preserved exactly over long time scales in the case of inviscid fluids. These integrators can be viewed as an extension of the classical energy-preserving Harlow-Welch / Crank-Nicolson scheme to simplicial grids.
PDF, 443KB
The authors wish to thank our collaborators Eva Kanso and Jerrold E. Marsden, as well as Kevin Beason for his renderer Pane. This research was partially funded by the NSF (CCF-0811373/0811313, DMS-0453145, and CMMI-0757106/0757123), the DOE (DE-FG02-04ER25657), and Pixar.
@article{Mullen:2009:EIF, author = {Mullen, Patrick and Crane, Keenan and Pavlov, Dmitry and Tong, Yiying and Desbrun, Mathieu}, title = {Energy-preserving integrators for fluid animation}, journal = {ACM Trans. Graph.}, volume = {28}, issue = {3}, month = {July}, year = {2009}, }
Initial and final vorticity from teapot animation.
Targa, 164KB
Additional Images