Novel Paradigms in Physics-Based Animation: Pointwise Divergence-Free Fluid Advection and Mixed-Dimensional Elastic Object Simulation

This thesis explores important but so far less studied aspects of physics-based animation: a simulation method for mixed-dimensional and/or non-manifold elastic objects, and a pointwise divergence-free velocity interpolation method applied to fluid simulation. Considering the popularity of single-type models e.g., hair, cloths, soft bodies, etc., in deformable body simulations, more complicated coupled models have gained less attention in graphics research, despite their relative ubiquity in daily life. This thesis presents a unified method to simulate such models: elastic bodies consisting of mixed-dimensional components represented with potentially non-manifold simplicial meshes. Building on well-known simplicial rod, shell, and solid models, this thesis categorizes and defines a comprehensive palette expressing all possible constraints and elastic energies for stiff and flexible connections between the 1D, 2D, and 3D components of a single conforming simplicial mesh. For fluid animation, this thesis proposes a novel methodology to enhance grid-based fluid animation with pointwise divergence-free velocity interpolation. Unlike previous methods which interpolate discrete velocity values directly for advection, this thesis proposes using intermediate steps involving vector potentials: first build a discrete vector potential field, interpolate these values to form a pointwise potential, and apply the continuous curl to recover a pointwise divergence-free flow field. Particles under these pointwise divergence-free flows exhibit significantly better particle distributions than divergent flows over time. To accelerate the use of vector potentials, this thesis proposes an efficient method that provides boundary-satisfying and smooth discrete potential fields on uniform and cut-cell grids. This thesis also introduces an improved ramping strategy for the “Curl-Noise” method of Bridson et al. (2007), which enforces exact no-normal-flow on the exterior domain boundaries and solid surfaces. The ramping method in the thesis effectively reduces the incidence of particles colliding with obstacles or creating erroneous gaps around the obstacles, while significantly alleviating the artifacts the original ramping strategy produces

Curl-Flow: Boundary-Respecting Pointwise Incompressible Velocity Interpolation for Grid-Based Fluids

We propose to augment standard grid-based fluid solvers with pointwise divergence-free velocity interpolation, thereby ensuring exact incompressibility down to the sub-cell level. Our method takes as input a discretely divergence-free velocity field generated by a staggered grid pressure projection, and first recovers a corresponding discrete vector potential. Instead of solving a costly vector Poisson problem for the potential, we develop a fast parallel sweeping strategy to find a candidate potential and apply a gauge transformation to enforce the Coulomb gauge condition and thereby make it numerically smooth. Interpolating this discrete potential generates a pointwise vector potential whose analytical curl is a pointwise incompressible velocity field. Our method further supports irregular solid geometry through the use of level set-based cut-cells and a novel Curl-Noise-inspired potential ramping procedure that simultaneously offers strictly non-penetrating velocities and incompressibility. Experimental comparisons demonstrate that the vector potential reconstruction procedure at the heart of our approach is consistently faster than prior such reconstruction schemes, especially those that solve vector Poisson problems. Moreover, in exchange for its modest extra cost, our overall Curl-Flow framework produces significantly improved particle trajectories that closely respect irregular obstacles, do not suffer from spurious sources or sinks, and yield superior particle distributions over time

A Unified Simplicial Model for Mixed-Dimensional and Non-Manifold Deformable Elastic Objects

Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from [email protected] present a unified method to simulate deformable elastic bodies consisting of mixed-dimensional components represented with potentially non-manifold simplicial meshes. Building on well-known simplicial rod, shell, and solid models for elastic continua, we categorize and define a comprehensive palette expressing all possible constraints and elastic energies for stiff and flexible connections between the 1D, 2D, and 3D components of a single conforming simplicial mesh. This palette consists of three categories: point connections, in which simplices meet at a single vertex around which they may twist and bend; curve connections in which simplices share an edge around which they may rotate (bend) relative to one another; and surface connections, in which a shell is embedded on or into a solid. To define elastic behaviors across non-manifold point connections, we adapt and apply parallel transport concepts from elastic rods. To address discontinuous forces that would otherwise arise when large accumulated relative rotations wrap around in the space of angles, we develop an incremental angle-update strategy. Our method provides a conceptually simple, flexible, and highly expressive framework for designing complex elastic objects, by modeling the geometry with a single simplicial mesh and decorating its elements with appropriate physical models (rod, shell, solid) and connection types (point, curve, surface). We demonstrate a diverse set of possible interactions achievable with our method, through technical and application examples, including scenes featuring complex aquatic creatures, children's toys, and umbrellas.This work was supported in part by the Natural Sciences and Engineering Research Council of Canada (RGPIN-04360-2014