Physics-Based Animation

In this paper, we present a point-based method for animating incompressible flow. The advection term is handled by moving the sample points through the flow in a Lagrangian fashion. However, unlike most previous approaches, the pressure term is handled by performing a projection onto a divergence-free field. To perform the pressure projection, we compute a Voronoi diagram with the sample points as input. Borrowing from Finite Volume Methods, we then invoke the divergence theorem and ensure that each Voronoi cell is divergence free. To handle complex boundary conditions, Voronoi cells are clipped against obstacle boundaries and free surfaces. The method is stable, flexible and combines many of the desirable features of point-based and grid-based methods. We demonstrate our approach on several examples of splashing and streaming liquid and swirling smoke.

A Point-based method for Animating Incompressible Flow

Recent Particle System models have evolved toward accurate representation of elastic stiffness based on continuum mechanics, converging to formulations that make them quite analogous to fast Finite Element methods. These formulations usually involve the linearization of tensors that help their formulation in the context of linear elasticity. Toward our objective of simulating the nonlinear properties of cloth accurately, we show through this work that this linearization can indeed be suppressed and replaced by adapted strain-stress laws relating precisely the nonlinear behavior of the material. This leads to very streamlined computations that are particularly efficient for simulating the nonlinear anisotropic tensile elasticity of deformable surfaces. Through a simple and efficient implementation using the Particle System formalism, we demonstrate the efficiency of this method with examples related to garment simulation.

Simple, yet accurate tensile stiffness

We present a practical approach to isotropic tetrahedral meshing of 3D domains bounded by piecewise smooth surfaces. Building upon recent theoretical and practical advances, our algorithm interleaves Delaunay refinement and mesh optimization to generate quality meshes that satisfy a set of user-defined criteria. This interleaving is shown to be more conservative in number of Steiner point insertions than refinement alone, and to produce higher quality meshes than optimization alone. A careful treatment of boundaries and their features is presented, offering a versatile framework for designing smoothly graded tetrahedral meshes.

Interleaving Delaunay Refinement and Optimization for Practical Isotropic Tetrahedron Mesh Generation

We develop a method for reliable simulation of elastica in complex contact scenarios. Our focus is on firmly establishing three parameter-independent guarantees: that simulations of well-posed problems (a) have no interpenetrations, (b) obey causality, momentum- and energy-conservation laws, and (c) complete in finite time. We achieve these guarantees through a novel synthesis of asynchronous variational integrators, kinetic data structures, and a discretization of the contact barrier potential by an infinite sum of nested quadratic potentials. In a series of two- and three-dimensional examples, we illustrate that this method more easily handles challenging problems involving complex contact geometries, sharp features, and sliding during extremely tight contact.

Asynchronous Contact Mechanics

This paper introduces a data-driven representation and modeling technique for simulating non-linear heterogeneous soft tissue. It simplifies the construction of convincing deformable models by avoiding complex selection and tuning of physical material parameters, yet retaining the richness of non-linear heterogeneous behavior. We acquire a set of example deformations of a real object, and represent each of them as a spatially varying stress-strain relationship in a finite-element model. We then model the material by non-linear interpolation of these stress-strain relationships in strain-space. Our method relies on a simple-to-build capture system and an efficient run-time simulation algorithm based on incremental loading, making it suitable for interactive computer graphics applications. We present the results of our approach for several nonlinear materials and biological soft tissue, with accurate agreement of our model to the measured data.

Capture and Modeling of Non-Linear Heterogeneous Soft Tissue

It’s that time of year again… Here’s the link to Ke-Sen’s more complete list.

Physics-related papers:

• Energy-Preserving Integrators for Fluid Animation
• Interactive Simulation of Surgical Needle Insertion and Steering
• Modular Bases for Fluid Dynamics
• Numerical Coarsening of Inhomogeneous Elastic Materials
• Capture and Modeling of Non-Linear Heterogeneous Soft Tissue
• Interleaving Delaunay Refinement and Optimization for Practical Isotropic Tetrahedron Mesh Generation
• Physically Guided Liquid Surface Modeling from Videos
• Deforming Meshes that Split and Merge
• Predictive-Corrective Incompressible SPH
• Preserving Topology and Elasticity for Embedded Deformable Models
• Deformable Object Animation Using Reduced Optimal Control
• Enrichment Textures for Detailed Cutting of Shells
• Directable, High-Resolution Simulation of Fire on the GPU
• Asynchronous Contact Mechanics
• Detail-Preserving Continuum Simulation of Straight Hair

We present a method for simulating highly detailed cutting and fracturing of thin shells using low-resolution simulation meshes. Instead of refining or remeshing the underlying simulation domain to resolve complex cut paths, we adapt the extended finite element method (XFEM) and enrich our approximation by custom designed basis functions, while keeping the simulation mesh unchanged. The enrichment functions are stored in enrichment textures, which allows for fracture and cutting discontinuities at a resolution much finer than the underlying mesh, similar to image textures for increased visual resolution. Furthermore, we propose harmonic enrichment functions to handle multiple, intersecting, arbitrarily shaped, progressive cuts per element in a simple and unified framework. Our underlying shell simulation is based on discontinuous Galerkin (DG) FEM, which relaxes the restrictive requirement of C1 continuous basis functions and thus allows for simpler, C0 continuous XFEM enrichment functions.

Enrichment Textures for Detailed Cutting of Shells

We present a novel, incompressible fluid simulation method based on the Lagrangian Smoothed Particle Hydrodynamics (SPH) model. In our method, incompressibility is enforced by using a prediction-correction scheme to determine the particle pressures. For this, the information about density fluctuations is actively propagated through the fluid and pressure values are updated until the targeted density is satisfied. With this approach, we avoid the computational expenses of solving a pressure Poisson equation, while still being able to use large time steps in the simulation. The achieved results show that our predictive-corrective incompressible SPH (PCISPH) method clearly outperforms the commonly used weakly compressible SPH (WCSPH) model by more than an order of magnitude while the computations are in good agreement with the WCSPH results.

Predictive-Corrective Incompressible SPH

We present a new approach to fluid simulation that balances the speed of model reduction with the flexibility of grid-based methods. We construct a set of composable reduced models, or tiles, which capture spatially localized fluid behavior. We then precompute coupling terms so that these models can be rearranged at runtime. To enforce consistency between tiles, we introduce constraint reduction. This technique modifies a reduced model so that a given set of linear constraints can be fulfilled. Because dynamics and constraints can be solved entirely in the reduced space, our method is extremely fast and scales to large domains.

Modular Bases for Fluid Dynamics

We present algorithms for simulating and visualizing the insertion and steering of needles through deformable tissues for surgical training and planning. Needle insertion is an essential component of manyclinical procedures such as biopsies, injections, neurosurgery, and brachytherapy cancer treatment. The success of these procedures dependson accurate guidanceofthe needletiptoaclinical target while avoiding vital tissues. Needle insertion deforms body tissues, making accurate placement difficult. Our interactive needle insertion simulator models the coupling between a steerable needle and deformable tissue. We introduce (1) a novel algorithm for local remeshing that quickly enforces the conformity of a tetrahedral mesh to a curvilinear needle path, enabling accurate computation of contact forces, (2) an efficient method for coupling a 3D finite element simulation with a 1D inextensible rod with stick-slip friction, and (3) optimizations that reduce the computation time for physically based simulations.We can realistically and interactively simulate needle insertion into a prostate mesh of 13,375 tetrahedra and 2,763 vertices at a 25 Hz frame rate on an 8-core 3.0 GHz Intel Xeon PC. The simulation models prostate brachytherapy with needles of varying stiffness, steering needles around obstacles, and supports motion planning for robotic needle insertion. We evaluate the accuracyof the simulation by comparing against real-world experiments in which flexible, steerable needles were inserted into gel tissue phantoms.

Interactive Simulation of Surgical Needle Insertion and Steering