Physics-Based Animation

Christopher Batty, Ben Houston

We present the first spatially adaptive Eulerian fluid animation method to support challenging viscous liquid effects such as folding, coiling, and variable viscosity. We propose a tetrahedral node-based embedded finite volume method for fluid viscosity, adapted from popular techniques for Lagrangian deformable objects. Applied in an Eulerian fashion with implicit integration, this scheme stably and efficiently supports high viscosity fluids while yielding symmetric positive definite linear systems. To integrate this scheme into standard tetrahedral mesh-based fluid simulators, which store normal velocities on faces rather than velocity vectors at nodes, we offer two methods to reconcile these representations. The first incorporates a mapping between different degrees of freedom into the viscosity solve itself. The second uses a FLIP-like approach to transfer velocity data between nodes and faces before and after the linear solve. The former offers tighter coupling by enabling the linear solver to act directly on the face velocities of the staggered mesh, while the latter provides a sparser linear system and a simpler implementation. We demonstrate the effectiveness of our approach with animations of spatially varying viscosity, realistic rotational motion, and viscous liquid buckling and coiling.

A Simple Finite Volume Method for Adaptive Viscous Liquids

Raphael Diziol, Jan Bender, Daniel Bayer

We introduce an efficient technique for robustly simulating incompressible objects with thousands of elements in real-time. Instead of considering a tetrahedral model, commonly used to simulate volumetric bodies, we simply use their surfaces. Not requiring hundreds or even thousands of elements in the interior of the object enables us to simulate more elements on the surface, resulting in high quality deformations at low computation costs. The elasticity of the objects is robustly simulated with a geometrically motivated shape matching approach which is extended by a fast summation technique for arbitrary triangle meshes suitable for an efficient parallel computation on the GPU. Moreover, we present an oscillation-free and collision-aware volume constraint, purely based on the surface of the incompressible body. The novel heuristic we propose in our approach enables us to conserve the volume, both globally and locally. Our volume constraint is not limited to the shape matching method and can be used with any method simulating the elasticity of an object. We present several examples which demonstrate high quality volume conserving deformations and compare the run-times of our CPU implementation, as well as our GPU implementation with similar methods.

Robust Real-Time Deformation of Incompressible Surface Meshes

Yubo Zhang, Carlos Correa, Kwan-Liu Ma

We present a novel graph-based data-driven technique for cost-effective fire modeling. This technique allows composing long animation sequences using a small number of short simulations. While traditional techniques such as motion graphs and motion blending work well for character motion synthesis, they cannot be trivially applied to fluids to produce results with physically consistent properties which are crucial to the visual appearance of fluids. Motivated by the motion graph technique used in character animations, we introduce a new type of graph which can be applied to create various fire phenomena. Each graph node consists of a group of compact spatialtemporal flow pathlines instead of a set of volumetric state fields. Consequently, achieving smooth transitions between discontinuous graph nodes for modeling turbulent fires becomes feasible and computationally efficient.The synthesized particle flow results allow direct particle controls which is much more flexible than a full volumetric representation of the simulation output. The accompanying video shows the versatility and potential power of this new technique for synthesizing realtime complex fire at the quality comparable to production animations.

Graph-based Fire Synthesis

The draft program for SCA 2011 is online here. Ke-Sen Huang maintains links to the full set of papers here.  Many of the papers involve physical simulation, including:

• Physics-based Character Skinning using Multi-Domain Subspace Deformations
• A Multigrid Fluid Pressure Solver Handling Separating Solid Boundary Conditions
• Mass and Momentum Conservation for Fluid Simulation
• Mathematical Foundation of the Optimization-Based Fluid Animation Method
• A Simple Finite Volume Method for Adaptive Viscous Liquids
• Procedural Modelling of Explosion Phenomena Based on Physical Properties
• Preview-based Sampling for Controlling Gaseous Simulations
• Graph-based Fire Synthesis
• A Particle-based Method for Preserving Fluid Sheets
• A Level-set Method for Skinning Animated Particle Data
• SPH Granular Flow with Friction and Cohesion
• Hybrid Smoothed Particle Hydrodynamics
• Optimization for Sag-Free Simulations
• Robust Real-Time Deformation of Incompressible Surface Meshes
• Asynchronous Integration with Phantom Meshes
• Element-Wise Mixed Implicit-Explicit Integration for Stable Dynamic Simulation of Deformable Objects

We propose a two-scale method for particle-based fluids that allocates computing resources to regions of the fluid where complex flow behavior emerges. Our method uses a low- and a high-resolution simulation that run at the same time. While in the coarse simulation the whole fluid is represented by large particles, the fine level simulates only a subset of the fluid with small particles. The subset can be arbitrarily defined and also dynamically change over time to capture complex flows and small-scale surface details. The low- and high-resolution simulations are coupled by including feedback forces and defining appropriate boundary conditions. Our method offers the benefit that particles are of the same size within each simulation level. This avoids particle splitting and merging processes, and allows the simulation of very large resolution differences without any stability problems. The model is easy to implement, and we show how it can be integrated into a standard SPH simulation as well as into the incompressible PCISPH solver. Compared to the single-resolution simulation, our method produces similar surface details while improving the efficiency linearly to the achieved reduction rate of the particle number.

Two-Scale Particle Simulation

We present a new algorithm for near-interactive simulation of skeleton driven, high resolution elasticity models. Our methodology is used for soft tissue deformation in character animation. The algorithm is based on a novel discretization of corotational elasticity over a hexahedral lattice. Within this framework we enforce positive definiteness of the stiffness matrix to allow efficient quasistatics and dynamics. In addition, we present a multigrid method that converges with very high efficiency. Our design targets performance through parallelism using a fully vectorized and branch-free SVD algorithm as well as a stable one-point quadrature scheme. Since body collisions, self collisions and soft-constraints are necessary for real-world examples, we present a simple framework for enforcing them. The whole approach is demonstrated in an end-to-end production-level character skinning system.

Efficient Elasticity for Character Skinning with Contact and Collisions

This paper shows a method to extend 3D nonlinear elasticity model reduction to open-loop multi-level reduced deformable structures. Given a volumetric mesh, we decompose the mesh into several subdomains, build a reduced deformable model for each domain, and connect the domains using inertia coupling. This makes model reduction deformable simulations much more versatile: localized deformations can be supported without prohibitive computational costs, parts can be re-used and precomputation times shortened. Our method does not use constraints, and can handle large domain rigid body motion in addition to large deformations, due to our derivation of the gradient and Hessian of the rotation matrix in polar decomposition. We show real-time examples with multi-level domain hierarchies and hundreds of reduced degrees of freedom.

Real-time Large-deformation Substructuring

We present a new Eulerian fluid simulation method, which allows real-time simulations of large scale three dimensional liquids. Such scenarios have hitherto been restricted to the domain of off-line computation. To reduce computation time we use a hybrid grid representation composed of regular cubic cells on top of a layer of tall cells. With this layout water above an arbitrary terrain can be represented without consuming an excessive amount of memory and compute power, while focusing effort on the area near the surface where it most matters. Additionally, we optimized the grid representation for a GPU implementation of the fluid solver. To further accelerate the simulation, we introduce a specialized multigrid algorithm for solving the Poisson equation and propose solver modifications to keep the simulation stable for large time steps. We demonstrate the efficiency of our approach in several real-world scenarios, all running above 30 frames per second on a modern GPU. Some scenes include additional features such as two-way rigid body coupling as well as particle representations of sub-grid detail.

Real-Time Eulerian Water Simulation Using a Restricted Tall Cell Grid

We propose a new fast and robust method to simulate various types of solid including rigid, plastic and soft bodies as well as one, two and three dimensional structures such as ropes, cloth and volumetric objects. The underlying idea is to use oriented particles that store rotation and spin, along with the usual linear attributes, i.e. position and velocity. This additional information adds substantially to traditional particle methods. First, particles can be represented by anisotropic shapes such as ellipsoids, which approximate surfaces more accurately than spheres. Second, shape matching becomes robust for sparse structures such as chains of particles or even single particles because the undefined degrees of freedom are captured in the rotational states of the particles. Third, the full transformation stored in the particles, including translation and rotation, can be used for robust skinning of graphical meshes and for transforming plastic deformations back into the rest state.

Solid Simulation with Oriented Particles

We introduce Hodge-optimized triangulations (HOT), a family of well-shaped primal-dual pairs of complexes designed for fast and accurate computations in computer graphics. Previous work most commonly employs barycentric or circumcentric duals; while barycentric duals guarantee that the dual of each simplex lies within the simplex, circumcentric duals are often preferred due to the induced orthogonality between primal and dual complexes. We instead promote the use of weighted duals (“power diagrams”). They allow greater flexibility in the location of dual vertices while keeping primal-dual orthogonality, thus providing a valuable extension to the usual choices of dual by only adding one additional scalar per primal vertex. Furthermore, we introduce a family of functionals on pairs of complexes that we derive from bounds on the errors induced by diagonal Hodge stars, commonly used in discrete computations. The minimizers of these functionals, called HOT meshes, are shown to be generalizations of Centroidal Voronoi Tesselations and Optimal Delaunay Triangulations, and to provide increased accuracy and flexibility for a variety of computational purposes.

HOT: Hodge-Optimized Triangulations