Guide Shapes for High Resolution Naturalistic Liquid Simulation

Art direction of high resolution naturalistic liquid simulations is notoriously hard, due to both the chaotic nature of the physics and the computational resources required. Resimulating a scene at higher resolution often produces very different results, and is too expensive to allow many design cycles. We present a method of constraining or guiding a high resolution liquid simulation to stay close to a finalized low resolution version (either simulated or directly animated), restricting the solve to a thin outer shell of liquid around a guide shape. Our method is generally faster than an unconstrained simulation and can be integrated with a standard fluid simulator. We demonstrate several applications, with both simulated and hand-animated inputs.

Guide Shapes for High Resolution Naturalistic Liquid Simulation

SIGGRAPH 2011

Ke-Sen Huang’s  list of 2011 SIGGRAPH papers is here.

The physics-related SIGGRAPH papers so far include:

The  TOG papers being presented at SIGGRAPH related to physics:

A Nonsmooth Newton Solver for Capturing Exact Coulomb Friction in Fiber Assemblies

We focus on the challenging problem of simulating thin elastic rods in contact, in the presence of friction. Most previous approaches in computer graphics rely on a linear complementarity formulation for handling contact in a stable way, and approximate Coulombs’s friction law for making the problem tractable. In contrast, following the seminal work by Alart and Curnier in contact mechanics, we simultaneously model contact and exact Coulomb friction as a zero finding problem of a nonsmooth function. A semi-implicit time-stepping scheme is then employed to discretizethe dynamics of rods constrained by frictional contact: this leads to a set of linear equations subject to an equality constraint involving a non-differentiable function. To solve this one-step problem we introduce a simple and practical nonsmooth Newton algorithm, which proves to be reasonably efficient and robust for systems that are not over-constrained. We show that our method is able to finely capture the subtle effects that occur when thin elastic rods with various geometries enter into contact, such as stick-slip instabilities in free configurations, entangling curls, resting contacts in braid-like structures, or the formation of tight knots under large constraints. Our method can be viewed as a first step towards the accurate modeling of dynamic fibrous materials.

A Nonsmooth Newton Solver for Capturing Exact Coulomb Friction in Fiber Assemblies

Optimization-based Fluid Simulation on Unstructured Meshes

We present a novel approach to fluid simulation, allowing us to take into account the surface energy in a precise manner. This new approach combines a novel, topology-adaptive approach to deformable interface tracking, called the deformable simplicial complexes method (DSC) with an optimization-based, linear finite element method for solving the incompressible Euler equations. The deformable simplicial complexes track the surface of the fluid: the fluid-air interface is represented explicitly as a piecewise linear surface which is a subset of tetrahedralization of the space, such that the interface can be also represented implicitly as a set of faces separating tetrahedra marked as inside from the ones marked as outside. This representation introduces insignificant and controllable numerical diffusion, allows robust topological adaptivity and provides both a volumetric finite element mesh for solving the fluid dynamics equations as well as direct access to the interface geometry data, making inclusion of a new surface energy term feasible. Furthermore, using an unstructured mesh makes it straightforward to handle curved solid boundaries and gives us a possibility to explore several fluid-solid interaction scenarios.

Optimization-based Fluid Simulation on Unstructured Meshes

Langevin Particle: A Self-Adaptive Lagrangian Primitive For Flow Simulation Enhancement

We develop a new Lagrangian primitive, named Langevin particle, to incorporate turbulent flow details in fluid simulation. A group of the particles are distributed inside the simulation domain based on a turbulence energy model with turbulence viscosity. A particle in particular moves obeying the generalized Langevin equation, a well-known stochastic differential equation that describes the particle’s motion as a random Markov process. The resultant particle trajectory shows self-adapted fluctuation in accordance to the turbulence energy, while following the global flow dynamics. We then feed back Langevin forces to the simulation based on the stochastic trajectory, which drive the flow with necessary turbulence. The new hybrid flow simulation method features nonrestricted particle evolution requiring minimal extra manipulation after initiation. The flow turbulence is easily controlled and the total computational overhead of enhancement is minimal based on typical fluid solvers.

Langevin Particle: A Self-Adaptive Lagrangian Primitive For Flow Simulation Enhancement

Hybrid Multiresolution Wire

We describe a method for the visual interactive simulation of wires contacting with rigid multibodies. The physical model used is a hybrid combining lumped elements and massless quasistatic representations. The latter is based on a kinematic constraint preserving the total length of the wire along a segmented path which can involve multiple bodies simultaneously and dry frictional contact nodes used for roping, lassoing and fastening. These nodes provide stick and slide friction along edges of the contacting geometries. The lumped element resolution is adapted dynamically based on local stability criteria, becoming coarser as the tension increases, and up to the purely kinematic representation. Kinematic segments and contact nodes are added and deleted and propagated based on contact geometries and dry friction configurations. The method gives dramatic increase on both performance and robustness because it quickly decimates superfluous nodes without loosing stability, yet adapts to complex configurations with many contacts and high curvature, keeping a fixed, large integration time step. Numerical results demonstrating the performance and stability of the adaptive multiresolution scheme are presented along with an array of representative simulation examples illustrating the versatility of the frictional contact model.

Hybrid Multiresolution Wire

Constraint Fluids

We present a fluid simulation method based on Smoothed Particle Hydrodynamics (SPH) in which incompressibility and boundary conditions are enforced using holonomic kinematic constraints on the density. This formulation enables systematic multiphysics integration in which interactions are modeled via similar constraints between the fluid pseudo-particles and impenetrable surfaces of other bodies. These conditions embody Archimede’s principle for solids and thus buoyancy results as a direct consequence. We use a variational time stepping scheme suitable for general constrained multibody systems we call SPOOK. Each step requires the solution of only one Mixed Linear Complementarity Problem (MLCP) with very few inequalities, corresponding to solid boundary conditions. We solve this MLCP with a fast iterative method. Overall stability is vastly improved in comparison to the unconstrained version of SPH, and this allows much larger time steps, and an increase in overall performance by two orders of magnitude. Proof of concept is given for computer graphics applications and interactive simulations.

Constraint Fluids

A Hexahedral Multigrid Approach for Simulating Cuts in Deformable Objects

We present a hexahedral finite element method for simulating cuts in deformable bodies using the corotational formulation of strain at high computational efficiency. Key to our approach is a novel embedding of adaptive element refinements and topological changes of the simulation grid into a geometric multigrid solver. Starting with a coarse hexahedral simulation grid, this grid is adaptively refined at the surface of a cutting tool until a finest resolution level, and the cut is modeled by separating elements along the cell faces at this level. To represent the induced discontinuities on successive multigrid levels, the affected coarse grid cells are duplicated and the resulting connectivity components are distributed to either side of the cut. Drawing upon recent work on octree and multigrid schemes for the numerical solution of partial differential equations, we develop efficient algorithms for updating the systems of equations of the adaptive finite element discretization and the multigrid hierarchy. To construct a surface that accurately aligns with the cuts, we adapt the splitting cubes algorithm to the specific linked voxel representation of the simulation domain we use. The paper is completed by a convergence analysis of the finite element solver and a performance comparison to alternative numerical solution methods. These investigations show that our approach offers high computational efficiency and physical accuracy, and that it enables cutting of deformable bodies at very high resolutions.

A Hexahedral Multigrid Approach for Simulating Cuts in Deformable Objects