Filament based smoke with vortex shedding and variational reconnection

Simulating fluids based on vortex filaments is highly attractive for the creation of special effects because it gives artists full control over the simulation using familiar tools like curve editors or the scripted generation of new vortex filaments over time. Because filaments offer a very compact description of fluid flow, real time applications like games or virtual reality are also possible. We present a complete model that includes moving obstacles with vortex shedding, all represented as filaments. Due to variational reconnection the long-time behavior of our method is excellent: Energy and momentum stay constant within reasonable bounds and computational complexity does not increase over time.

A Practical Simulation of Dispersed Bubble Flow

In this paper, we propose a simple and efficient framework for simulating dispersed bubble flow. Instead of modeling the complex hydrodynamics of numerous small bubbles explicitly, our method approximates the average motion of these bubbles using a continuum multiphase solver. Then, the subgrid interactions among bubbles are computed using our new stochastic solver. Using the proposed scheme, we can efficiently simulate complex scenes with millions of bubbles.

A Practical Simulation of Dispersed Bubble Flow

A Multiscale Approach to Mesh-based Surface Tension

We present an approach to simulate flows driven by surface tension based on triangle meshes. Our method consists of two simulation layers: the first layer is an Eulerian method for simulating surface tension forces that is free from typical strict time step constraints. The second simulation layer is a Lagrangian finite element method that simulates sub-grid scale wave details on the fluid surface. The surface wave simulation employs an unconditionally stable, symplectic time integration method that allows for a high propagation speed due to strong surface tension. Our approach can naturally separate the grid- and sub-grid scales based on a volume-preserving mean curvature flow. As our model for the sub-grid dynamics enforces a local conservation of mass, it leads to realistic pinch off and merging effects. In addition to this method for simulating dynamic surface tension effects, we also present an efficient non-oscillatory approximation for capturing damped surface tension behavior. These approaches allow us to efficiently simulate complex phenomena associated with strong surface tension, such as Rayleigh-Plateau instabilities and crown splashes, in a short amount of time.

A Multiscale Approach to Mesh-based Surface Tension

Matching Fluid Simulation Elements to Surface Geometry and Topology

We introduce an Eulerian liquid simulation framework based on the Voronoi diagram of a potentially unorganized collection of pressure samples. Constructing the simulation mesh in this way allows us to place samples anywhere in the computational domain; we exploit this by choosing samples that accurately capture the geometry and topology of the liquid surface. When combined with high-resolution explicit surface tracking this allows us to simulate nearly arbitrarily thin features, while eliminating noise and other artifacts that arise when there is a resolution mismatch between the simulation and the surface—and allowing a precise inclusion of surface tension based directly on and at the same resolution as the surface mesh. In addition, we present a simplified Voronoi/Delaunay mesh velocity interpolation scheme, and a direct extension of embedded free surfaces and solid boundaries to Voronoi meshes.

Matching Fluid Simulation Elements to Surface Geometry and Topology

Fluid Simulation with Articulated Bodies

We present an algorithm for creating realistic animations of characters that are swimming through fluids. Our approach combines dynamic simulation with data-driven kinematic motions (motion capture data) to produce realistic animation in a fluid. The interaction of the articulated body with the fluid is performed by incorporating joint constraints with rigid animation and by extending a solid/fluid coupling method to handle articulated chains. Our solver takes as input the current state of the simulation and calculates the angular and linear accelerations of the connected bodies needed to match a particular motion sequence for the articulated body. These accelerations are used to estimate the forces and torques that are then applied to each joint. Based on this approach, we demonstrate simulated swimming results for a variety of different strokes, including crawl, backstroke, breaststroke and butterfly. The ability to have articulated bodies interact with fluids also allows us to generate simulations of simple water creatures that are driven by simple controllers.

Fluid Simulation with Articulated Bodies

Accurate Tangential Velocities for Solid-Fluid Coupling

We propose a novel method for obtaining more accurate tangential velocities for solid fluid coupling. Our method works for both rigid and deformable objects as well as both volumetric objects and thin shells. The fluid can be either one phase such as smoke or two phase such as water with a free surface. The coupling between the solid and the fluid can either be one-way with kinematic solids or fully two-way coupled. The only previous scheme that was general enough to handle both two-way coupling and thin shells required a mass lumping strategy that did not allow for freely flowing tangential velocities. Similar to that previous work, our method prevents leaking of fluid across a thin shell, however unlike that work our method does not couple the tangential velocities in any fashion, allowing for the proper slip independently on each side of the body. Moreover, since it accurately and directly treats the tangential velocity, it does not rely on grid refinement to obtain a reasonable solution. Therefore, it gives a highly improved result on coarse meshes.

Accurate Tangential Velocities for Solid-Fluid Coupling

A Point-based Method for Animating Elastoplastic Solids

In this paper we describe a point-based approach for animating elastoplastic materials. Our primary contribution is a simple method for computing the deformation gradient for each particle in the simulation. The deformation gradient is computed for each particle by finding the affine transformation that best approximates the motion of neighboring particles over a single timestep. These transformations are then composed to compute the total deformation gradient that describes the deformation around a particle over the course of the simulation. Given the deformation gradient we can apply arbitrary constitutive models and compute the resulting elastic forces. Our method has two primary advantages: we do not store or compare to an initial rest configuration and we work directly with the deformation gradient. The first advantage avoids poor numerical conditioning and the second naturally leads to a multiplicative model of deformation appropriate for finite deformations. We demonstrate our approach on a number of examples that exhibit a wide range of material behaviors.

A Point-based Method for Animating Elastoplastic Solids

A Point-based Method for Animating Incompressible Flow

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

Simulation of Two-Phase Flow with Sub-Scale Droplets and Bubble Effects

We present a new Eulerian-Lagrangian method for physics-based simulation of fluid flow, which includes automatic generation of sub-scale spray and bubbles. The Marker Level Set method is used to provide a simple geometric criterion for free marker generation. A filtering method, inspired from Weber number thresholding, further controls the free marker generation (in a physics-based manner). Two separate models are used, one for sub-scale droplets, the other for sub-scale bubbles. Droplets are evolved in a Newtonian manner, using a density extension drag force field, while bubbles are evolved using a model based on Stokes’ Law. We show that our model for sub-scale droplet and bubble dynamics is simple to couple with a full (macro-scale) Navier-Stokes two-phase flow model and is quite powerful in its applications. Our animations include coarse grained multiphase features interacting with fine scale multiphase features.

Simulation of Two-Phase Flow with Sub-Scale Droplets and Bubble Effects

Real-Time Fluid Simulation Using Discrete Sine/Cosine Transforms

Recent advances in fluid simulations have yielded exceptionally realistic imagery. However, most algorithms have computational requirements that are prohibitive for real-time simulations. Using Fourier based solutions mitigates this issue, although due to wraparound, boundary conditions are not naturally available, leading to inconsistencies near the boundary. We show that slip boundary conditions can be imposed by solving the mass conservation step using cosine and sine transforms instead of the Fourier transform.
Further, we show that measures against density dissipation can be computed using cosine transforms and we describe a new method to compute surface tension in the same domain. This combination of related algorithms leads to real-time simulations with boundary conditions.

Real-Time Fluid Simulation Using Discrete Sine/Cosine Transforms