Home Pub.s Proj.s   CV  
[Hide abstracts] [All in Bibtex format]
Sections: [Refereed Papers] [All others]

Journal / Refereed Papers:
[Back to top...]

- 2017 -


Data-Driven Synthesis of Smoke Flows with CNN-based Feature Descriptors   [CITE]
Mengyu Chu, Nils Thuerey
Transaction on Graphics (SIGGRAPH); ACM, 2017

Abstract: We present a novel data-driven algorithm to synthesize high resolution flow simulations with reusable repositories of space-time flow data. In our work, we employ a descriptor learning approach to encode the similarity between fluid regions with differences in resolution and numerical viscosity. We use convolutional neural networks to generate the descriptors from fluid data such as smoke density and flow velocity. At the same time, we present a deformation limiting patch advection method which allows us to robustly track deformable fluid regions. With the help of this patch advection, we generate stable space-time data sets from detailed fluids for our repositories. We can then use our learned descriptors to quickly localize a suitable data set when running a new simulation. This makes our approach very efficient, and resolution independent. We will demonstrate with several examples that our method yields volumes with very high effective resolutions, and non-dissipative small scale details that naturally integrate into the motions of the underlying flow.

Perceptual Evaluation of Liquid Simulation Methods   [CITE]
Kiwon Um, Xiangyu Hu, Nils Thuerey
Transaction on Graphics (SIGGRAPH); ACM, 2017

Abstract: This paper proposes a novel framework to evaluate fluid simulation methods based on crowd-sourced user studies in order to robustly gather large numbers of opinions. The key idea for a robust and reliable evaluation is to use a reference video from a carefully selected real-world setup in the user study. By conducting a series of controlled user studies and comparing their evaluation results, we observe various factors that affect the perceptual evaluation. Our data show that the availability of a reference video makes the evaluation consistent. We introduce this approach for computing scores of simulation methods as visual accuracy metric. As an application of the proposed framework, a variety of popular simulation methods are evaluated.

Robust eXtended Finite Elements for Complex Cutting of Deformables   [CITE]
Dan Koschier, Jan Bender, Nils Thuerey
Transaction on Graphics (SIGGRAPH); ACM, 2017

Abstract: In this paper we present a novel approach to simulate cutting of deformable solids in virtual environments. On the basis of the eXtended Finite Element Method (XFEM) we develop a framework that can accurately handle complex cut paths. A particular strength of our method is that there is no requirement to modify either topology or geometry of the underlying discretization mesh. Our approach makes three-dimensional XFEM simulations possible with the help of a modified Moment Fitting scheme. Our moment fitting approach allows us to construct highly accurate quadrature schemes for discontinuous integrands on polyhedral domains. Our approach is very efficient, and does not require an explicit discretization of the cut surface within each element. Moreover, we arrive at a well-conditioned equation system with our approach which allows for stable, implicit time integration with large time steps. We demonstrate that our method yields results that are significantly more accurate than approaches based on node duplication. In contrast to methods employing geometric refinement, we realize very detailed cutting effects without changing the computational mesh in any way. Finally, we verify the consistency and robustness of our approach on simple and reproducible academic examples, while we demonstrate scalability, stability and applicability of our method with complex cuts of highly detailed structures in large-scale scenarios.

Pre-computed Liquid Spaces with Generative Neural Networks   [CITE]
Boris Bonev, Lukas Prantl, Nils Thuerey
arXiv; , 2017
Preprint: [PDF]    App: [WWW]   

Abstract: Liquids exhibit complex non-linear behavior under changing simulation conditions such as user interactions. We propose a method to map this complex behavior over a parameter range onto reduced representation based on space-time deformations. In order to represent the complexity of the full space of inputs, we leverage the power of generative neural networks to learn a reduced representation. We introduce a novel deformation-aware loss function, which enables optimization in the highly non-linear space of multiple deformations. To demonstrate the effectiveness of our approach, we showcase the method with several complex examples in two and four dimensions. Our representation makes it possible to generate implicit surfaces of liquids very efficiently, which makes it possible to display the scene from any angle, and to add secondary effects such as particle systems. We have implemented a mobile application for our full output pipeline to demonstrate that real-time interaction is possible with our approach.

Liquid Splash Modeling with Neural Networks   [CITE]
Kiwon Um, Xiangyu Hu, Nils Thuerey
arXiv; , 2017
Preprint: [PDF]   

Abstract: This paper proposes a new data-driven approach for modeling detailed splashes for liquid simulations with neural networks. Our model learns to generate small-scale splash detail for fluid-implicit-particle methods using training data acquired from physically accurate, high-resolution simulations. We use neural networks to model the regression of splash formation using a classifier together with a velocity modification term. More specifically, we employ a heteroscedastic model for the velocity updates. Our simulation results demonstrate that our model significantly improves visual fidelity with a large amount of realistic droplet formation and yields splash detail much more efficiently than finer discretizations. We show this for two different spatial scales and simulation setups.

Primal-Dual Optimization for Fluids   [CITE]
Tiffany Inglis, Marie-Lena Eckert, James Gregson, Nils Thuerey
Computer Graphics Forum; Wiley, 2017
Preprint: [PDF]    Paper: [WWW]    Video: [WWW]   

Abstract: We apply a novel optimization scheme from the image processing and machine learning areas, a fast Primal-Dual method, to achieve controllable and realistic fluid simulations. While our method is generally applicable to many problems in fluid simulations, we focus on the two topics of fluid guiding and separating solid-wall boundary conditions. Each problem is posed as an optimization problem and solved using our method, which contains acceleration schemes tailored to each problem. In fluid guiding, we are interested in partially guiding fluid motion to exert control while preserving fluid characteristics. With our method, we achieve explicit control over both large-scale motions and small-scale details which is valuable for many applications, such as level-of-detail adjustment (after running the coarse simulation), spatially varying guiding strength, domain modification, and resimulation with different fluid parameters. For the separating solid-wall boundary conditions problem, our method effectively eliminates unrealistic artefacts of fluid crawling up solid walls and sticking to ceilings, requiring few changes to existing implementations. We demonstrate the fast convergence of our Primal-Dual method with a variety of test cases for both model problems.

- 2016 -


SMASH: Physics-guided Reconstruction of Collisions from Videos   [CITE]
Aron Monszpart, Nils Thuerey, Niloy Mitra
Transactions on Graphics (SIGGRAPH Asia); ACM, 2016
Paper: [PDF]    Video: [WWW]    Video: [AVI]    Project: [WWW]   

Abstract: Collision sequences are commonly used in games and entertainment to add drama and excitement. Authoring even two body collisions in real world can be difficult, as one has to get timing and the object trajectories to be correctly synchronized. After tedious trial-and-error iterations, when objects can actually be made to collide, then they are difficult to capture in 3D. In contrast, synthetically generating plausible collisions is difficult as it requires adjusting different collision parameters (e.g., object mass ratio, coefficient of restitution, etc.) and appropriate initial parameters. We present SMASH to directly "read off" appropriate collision parameters simply based on input video recordings. Technically we enable this by utilizing laws of rigid body collision to regularize the problem of lifting 2D trajectories to a 3D reconstruction of the collision. The reconstructed sequences can then be modified and combined to easily author novel and plausible collisions. We evaluate our system on a range of synthetic scenes and demonstrate the effectiveness of our method by accurately reconstructing several complex real world collision events.

Dispersion Kernels for Water Wave Simulation   [CITE]
Jose-Angel Canabal Delgado, David Miraut, Nils Thuerey, Theodore Kim, Javier Portilla, Miguel Otaduy
Transactions on Graphics (SIGGRAPH Asia); ACM, 2016
Paper: [PDF]    Video: [MP4]    Video: [WWW]    Project: [WWW]   

Abstract: We propose a method to simulate the rich, scale-dependent dynamics of water waves. Our method preserves the dispersion properties of real waves, yet it supports interactions with obstacles and is computationally efficient. Fundamentally, it computes wave accelerations by way of applying a dispersion kernel as a spatially variant filter, which we are able to compute efficiently using two core technical contributions. First, we design novel, accurate, and compact pyramid kernels which compensate for low-frequency truncation errors. Second, we design a shadowed convolution operation that efficiently accounts for obstacle interactions by modulating the application of the dispersion kernel. We demonstrate a wide range of behaviors, which include capillary waves, gravity waves, and interactions with static and dynamic obstacles, all from within a single simulation.

Interpolations of Smoke and Liquid Simulations   [CITE]
Nils Thuerey
Transactions on Graphics; ACM, 2016
Preprint: [PDF]    ACM: [PDF]    Source code: [WWW]    Video: [WWW]   

Abstract: We present a novel method to interpolate smoke and liquid simulations in order to perform data-driven fluid simulations. Our approach calculates a dense space-time deformation using grid-based signed-distance functions of the inputs. A key advantage of this implicit Eulerian representation is that it allows us to use powerful techniques from the optical flow area. We employ a five-dimensional optical flow solve. In combination with a projection algorithm, and residual iterations, we achieve a robust matching of the inputs. Once the match is computed, arbitrary in between variants can be created very efficiently. To concatenate multiple long-range deformations, we propose a novel alignment technique. Our approach has numerous advantages, including automatic matches without user input, volumetric deformations that can be applied to details around the surface, and the inherent handling of topology changes. As a result, we can interpolate swirling smoke clouds, and splashing liquid simulations. We can even match and interpolate phenomena with fundamentally different physics: a drop of liquid, and a blob of heavy smoke.

Example-based Plastic Deformation of Rigid Bodies   [CITE]
Ben Jones, Nils Thuerey, Tamar Shinar, Adam Bargteil
Transactions on Graphics (SIGGRAPH); ACM, 2016
Paper: [PDF]    Project: [WWW]    Video: [WWW]   

Abstract: Physics-based animation is often used to animate scenes containing destruction of near-rigid, man-made materials. For these applications, the most important visual features are plastic deformation and fracture. Methods based on continuum mechanics model these materials as elastoplastic, and must perform expensive elasticity computations even though elastic deformations are imperceptibly small for rigid materials. We introduce an example-based plasticity model based on linear blend skinning that allows artists to author simulation objects using familiar tools. Dynamics are computed using an unmodified rigid body simulator, making our method computationally efficient and easy to integrate into existing pipelines. We introduce a flexible technique for mapping impulses computed by the rigid body solver to local, example-based deformations. For completeness, our method also supports prescoring based fracture. We demonstrate the practicality of our method by animating a variety of destructive scenes.

Narrow Band FLIP for Liquid Simulations   [CITE]
Florian Ferstl, Ryoichi Ando, Chris Wojtan, Ruediger Westermann, Nils Thuerey
Computer Graphics Forum (Eurographics'16); Eurographics Association, 2016
Paper: [PDF]    Project: [WWW]    Video: [WWW]    Code: [TGZ]   

Abstract: The Fluid Implicit Particle method (FLIP) for liquid simulations uses particles to reduce numerical dissipation and provide important visual cues for events like complex splashes and small-scale features near the liquid surface. Unfortunately, FLIP simulations can be computationally expensive, because they require a dense sampling of particles to fill the entire liquid volume. Furthermore, the vast majority of these FLIP particles contribute nothing to the fluid's visual appearance, especially for larger volumes of liquid. We present a method that only uses FLIP particles within a narrow band of the liquid surface, while efficiently representing the remaining inner volume on a regular grid. We show that a naive realization of this idea introduces unstable and uncontrollable energy fluctuations, and we propose a novel coupling scheme between FLIP particles and regular grid which overcomes this problem. Our method drastically reduces the particle count and simulation times while yielding results that are nearly indistinguishable from regular FLIP simulations. Our approach is easy to integrate into any existing FLIP implementation.

- 2015 -


Surface Turbulence for Particle-Based Liquid Simulations   [CITE]
Olivier Mercier, Cynthia Beauchemin, Nils Thuerey, Theodore Kim, Derek Nowrouzezahrai
Transactions on Graphics (SIGGRAPH Asia); ACM, 2015
Paper: [PDF]    Video: [AVI]    Video: [WWW]   

Abstract: We present a method to increase the apparent resolution of particlebased liquid simulations. Our method first outputs a dense, temporally coherent, regularized point set from a coarse particle-based liquid simulation. We then apply a surface-only Lagrangian wave simulation to this high-resolution point set. We develop novel methods for seeding and simulating waves over surface points, and use them to generate high-resolution details. We avoid error-prone surface mesh processing, and robustly propagate waves without the need for explicit connectivity information. Our seeding strategy combines a robust curvature evaluation with multiple bands of seeding oscillators, injects waves with arbitrarily fine-scale structures, and properly handles obstacle boundaries. We generate detailed fluid surfaces from coarse simulations as an independent post-process that can be applied to most particle-based fluid solvers.

A Stream Function Solver for Liquid Simulations   [CITE]
Ryoichi Ando, Nils Thuerey, and Chris Wojtan
Transactions on Graphics (SIGGRAPH); ACM, 2015
Paper: [PDF]    Video: [WWW]    Supp.-Vid.: [WWW]    Project: [WWW]   

Abstract: This paper presents a liquid simulation technique that enforces the incompressibility condition using a stream function solve instead of a pressure projection. Previous methods have used stream function techniques for the simulation of detailed single-phase flows, but a formulation for liquid simulation has proved elusive in part due to the free surface boundary conditions. In this paper, we introduce a stream function approach to liquid simulations with novel boundary conditions for free surfaces, solid obstacles, and solid-fluid coupling. Although our approach increases the dimension of the linear system necessary to enforce incompressibility, it provides interesting and surprising benefits. First, the resulting flow is guaranteed to be divergence-free regardless of the accuracy of the solve. Second, our free-surface boundary conditions guarantee divergence-free motion even in the un-simulated air phase, which enables two-phase flow simulation by only computing a single phase. We implemented this method using a variant of FLIP simulation which only samples particles within a narrow band of the liquid surface, and we illustrate the effectiveness of our method for detailed two-phase flow simulations with complex boundaries, detailed bubble interactions, and two-way solid-fluid coupling.

A Dimension-reduced Pressure Solver for Liquid Simulations   [CITE]
Ryoichi Ando, Nils Thuerey, and Chris Wojtan
Computer Graphics Forum (Eurographics); Eurographics Association, 2015
Paper: [PDF]    Video: [WWW]    Project: [WWW]   

Abstract: This work presents a method for efficiently simplifying the pressure projection step in a liquid simulation. We first devise a straightforward dimension reduction technique that dramatically reduces the cost of solving the pressure projection. Next, we introduce a novel change of basis that satisfies free-surface boundary conditions exactly, regardless of the accuracy of the pressure solve. When combined, these ideas greatly reduce the computational complexity of the pressure solve without compromising free surface boundary conditions at the highest level of detail. Our techniques are easy to parallelize, and they effectively eliminate the computational bottleneck for large liquid simulations.

- 2014 -


Blending Liquids   [CITE]
Karthik Raveendran, Nils Thuerey, Chris Wojtan, Greg Turk
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2014
Paper: [PDF]    Video: [WWW]   

Abstract: We present a method for smoothly blending between existing liquid animations. We introduce a semi-automatic method for matching two existing liquid animations, which we use to create new fluid motion that plausibly interpolates the input. Our contributions include a new space-time non-rigid iterative closest point algorithm that incorporates user guidance, a subsampling technique for efficient registration of meshes with millions of vertices, and a fast surface extraction algorithm that produces 3D triangle meshes from a 4D space-time surface. Our technique can be used to instantly create hundreds of new simulations, or to interactively explore complex parameter spaces. Our method is guaranteed to produce output that does not deviate from the input animations, and it generalizes to multiple dimensions. Because our method runs at interactive rates after the initial precomputation step, it has potential applications in games and training simulations.

From Capture to Simulation - Connecting Forward and Inverse Problems in Fluids   [CITE]
James Gregson, Nils Thuerey, Ivo Ihrke, Wolfgang Heidrich
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2014
Paper: [PDF]    Video: [AVI]    Video: [WWW]    Project: [WWW]   

Abstract: We explore the connection between fluid capture, simulation and proximal methods, a class of algorithms commonly used for inverse problems in image processing and computer vision. Our key finding is that the proximal operator constraining fluid velocities to be divergence-free is directly equivalent to the pressure-projection methods commonly used in incompressible flow solvers. This observation lets us treat the inverse problem of fluid tracking as a constrained flow problem all while working in an efficient, modular framework. In addition it lets us tightly couple fluid simulation into flow tracking, providing a global prior that significantly increases tracking accuracy and temporal coherence as compared to previous techniques. We demonstrate how we can use these improved results for a variety of applications, such as re-simulation, detail enhancement, and domain modification. We furthermore give an outlook of the applications beyond fluid tracking that our proximal operator framework could enable by exploring the connection of deblurring and fluid guiding.

- 2013 -


Highly Adaptive Liquid Simulations on Tetrahedral Meshes   [CITE]
Ryoichi Ando, Nils Thuerey, Chris Wojtan
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2013
Paper: [PDF]    Video: [MOV]    Video: [WWW]    Project: [WWW]   

Abstract: We introduce a new method for efficiently simulating liquid with extreme amounts of spatial adaptivity. Our method combines several key components to drastically speed up the simulation of large-scale fluid phenomena: We leverage an alternative Eulerian tetrahedral mesh discretization to significantly reduce the complexity of the pressure solve while increasing the robustness with respect to element quality and removing the possibility of locking. Next, we enable subtle free-surface phenomena by deriving novel second-order boundary conditions consistent with our discretization. We couple this discretization with a spatially adaptive Fluid-Implicit Particle (FLIP) method, enabling efficient, robust, minimally-dissipative simulations that can undergo sharp changes in spatial resolution while minimizing artifacts. Along the way, we provide a new method for generating a smooth and detailed surface from a set of particles with variable sizes. Finally, we explore several new sizing functions for determining spatially adaptive simulation resolutions, and we show how to couple them to our simulator. We combine each of these elements to produce a simulation algorithm that is capable of creating animations at high maximum resolutions while avoiding common pitfalls like inaccurate boundary conditions and inefficient computation.

Closest-Point Turbulence for Liquid Surfaces   [CITE]
Theodore Kim, Jerry Tessendorf, Nils Thuerey
ACM Transactions on Graphics; ACM Press, 2013
Paper: [PDF]    Video: [MOV]    Video: [WWW]    Project: [WWW]   

Abstract: We propose a method of increasing the apparent spatial resolution of an existing liquid simulation. Previous approaches to this "up-resing" problem have focused on increasing the turbulence of the underlying velocity field. Motivated by measurements in the free surface turbulence literature, we observe that past certain frequencies, it is sufficient to perform a wave simulation directly on the liquid surface, and construct a reduced-dimensional surface-only simulation. We sidestep the considerable problem of generating a surface parameterization by employing an embedding technique known as the Closest Point Method (CPM) that operates directly on a 3D extension field. The CPM requires 3D operators, and we show that for surface operators with no natural 3D generalization, it is possible to construct a viable operator using the inverse Abel transform. We additionally propose a fast, frozen core closest point transform, and an advection method for the extension field that reduces smearing considerably. Finally, we propose two turbulence coupling methods that seed the high resolution wave simulation in visually expected regions.

- 2012 -


Lagrangian Vortex Sheets for Animating Fluids   [CITE]
Tobias Pfaff, Nils Thuerey, Markus Gross
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2012
Paper: [PDF]    Video: [MP4]    Video: [WWW]   

Abstract: Buoyant turbulent smoke plumes with a sharp smoke-air interface, such as volcanic plumes, are notoriously hard to simulate. The surface clearly shows small-scale turbulent structures which are costly to resolve. In addition, the turbulence onset is directly visible at the interface, and is not captured by commonly used turbulence models. We present a novel approach that employs a triangle mesh as a high-resolution surface representation combined with a coarse Eulerian solver. On the mesh, we solve the interfacial vortex sheet equations, which allows us to accurately simulate buoyancy induced turbulence. For complex boundary conditions we propose an orthogonal turbulence model that handles vortices caused by obstacle interaction. In addition, we demonstrate a re-sampling scheme to remove surfaces that are hidden inside the bulk volume. In this way we are able to achieve highly detailed simulations of turbulent plumes efficiently.

Controlling Fluids using Meshes   [CITE]
Karthik Raveendran, Nils Thuerey, Chris Wojtan, Greg Turk
SCA '12: Proceedings of the 2012 ACM SIGGRAPH/Eurographics Symposium on Computer Animation; Eurographics Association, 2012
Paper: [PDF]    Video: [WWW]   

Abstract: We present an approach for artist-directed animation of liquids using multiple levels of control over the simulation, ranging from the overall tracking of desired shapes to highly detailed secondary effects such as dripping streams, separating sheets of fluid, surface waves and ripples. The first portion of our technique is a volume preserving morph that allows the animator to produce a plausible fluid-like motion from a sparse set of control meshes. By rasterizing the resulting control meshes onto the simulation grid, the mesh velocities act as boundary conditions during the projection step of the fluid simulation. We can then blend this motion together with uncontrolled fluid velocities to achieve a more relaxed control over the fluid that captures natural inertial effects. Our method can produce highly detailed liquid surfaces with control over sub-grid details by using a mesh-based surface tracker on top of a coarse grid-based fluid simulation. We can create ripples and waves on the fluid surface attracting the surface mesh to the control mesh with spring-like forces and also by running a wave simulation over the surface mesh. Our video results demonstrate how our control scheme can be used to create animated characters and shapes that are made of water.

- 2011 -


Preserving Fluid Sheets with Adaptively Sampled Anisotropic Particles   [CITE]
Ryoichi Ando, Nils Thuerey, Reiji Tsuruno
IEEE Transactions on Visualization and Computer Graphics; IEEE, 2011
Preprint: [PDF]    Video: [WWW]    Project: [WWW]   

Abstract: This paper presents a particle-based model for preserving fluid sheets of animated liquids with an adaptively sampled Fluid-Implicit-Particle (FLIP) method. In our method, we preserve fluid sheets by filling the breaking sheets with particle splitting in the thin regions, and by collapsing them in the deep water. To identify the critically thin parts, we compute the anisotropy of the particle neighborhoods, and use this information as a resampling criterion to reconstruct thin liquid surfaces. Unlike previous approaches, our method does not suffer from diffusive surfaces or complex re-meshing operations, and robustly handles topology changes with the use of a meshless representation. We extend the underlying FLIP model with an anisotropic position correction to improve the particle spacing, and adaptive sampling to efficiently perform simulations of larger volumes. Due to the Lagrangian nature of our method, it can be easily implemented and efficiently parallelized. The results show that our method can produce visually complex liquid animations with thin structures and vivid motions.

- 2010 -


Scalable Fluid Simulation using Anisotropic Turbulence Particles   [CITE]
Tobias Pfaff, Nils Thuerey, Jonathan Cohen, Sarah Tariq, Markus Gross
ACM Transactions on Graphics (SIGGRAPH Asia); ACM Press, 2010
Paper: [PDF]    Video: [MP4]    Video: [WWW]    Project: [WWW]   

Abstract: We present a novel, scalable turbulence method for fluid simulations that simulates small-scale detail based on a particle representation, which does not require neighborhood information. We compute transport of turbulent energy using a full two-equation k-e model with extensions for stability and for realistically capturing anisotropic turbulence effects. Turbulence is synthesized using a curl noise function, and we only require a very low grid resolution to resolve an underlying base flow. As we offload complexity from the fluid solver to the particle system, we can control the level of detail of the simulation easily by adjusting the particle number, without changing the large-scale behavior. In addition, no computations are wasted on areas that are not visible. We demonstrate that due to the design of our algorithm it is highly suitable for massively parallel architectures, and is able to generate detailed turbulent simulations with millions of particles at high frame rates.

A Multiscale Approach to Mesh-based Surface Tension Flows   [CITE]
Nils Thuerey, Chris Wojtan, Markus Gross, Greg Turk
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2010
Paper: [PDF]    Video: [MOV]    Video: [WWW]   

Abstract: 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.

Physics-Inspired Topology Changes for Thin Fluid Features   [CITE]
Chris Wojtan, Nils Thuerey, Markus Gross, Greg Turk
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2010
Paper: [PDF]    Video: [MOV]    Video: [WWW]   

Abstract: We propose a mesh-based surface tracking method for fluid animation that both preserves fine surface details and robustly adjusts the topology of the surface in the presence of arbitrarily thin features like sheets and strands. We replace traditional re-sampling methods with a convex hull method for connecting surface features during topological changes. This technique permits arbitrarily thin fluid features with minimal re-sampling errors by reusing points from the original surface. We further reduce re-sampling artifacts with a subdivision-based mesh-stitching algorithm, and we use a higher order interpolating subdivision scheme to determine the location of any newly-created vertices. The resulting algorithm efficiently produces detailed fluid surfaces with arbitrarily thin features while maintaining a consistent topology with the underlying fluid simulation.

- 2009 -


Synthetic Turbulence using Artificial Boundary Layers   [CITE]
Tobias Pfaff, Nils Thuerey, Andrew Selle, Markus Gross
ACM Transactions on Graphics (SIGGRAPH Asia); ACM Press, 2009
Paper: [PDF]    Video: [MOV]    Video: [WWW]    Sources: [TGZ]   

Abstract: Turbulent vortices in fluid flows are crucial for a visually interesting appearance. Although there has been a significant amount of work on turbulence in graphics recently, these algorithms rely on the underlying simulation to resolve the flow around objects. We build upon work from classical fluid mechanics to design an algorithm that allows us to accurately precompute the turbulence being generated around an object immersed in a flow. This is made possible by separating the flow into a mean and a fluctuating component, and relying on universal laws describing the flow near a wall. We precompute the turbulence characteristics around an object, and inject corresponding vorticity during a fluid simulation run, giving regions of confined vorticity. Then, a turbulence model is used to identify areas where this vorticity will transition into actual turbulence. We sample these regions with vortex particles, and simulate the further dynamics of the vortices based on these particles. We will show how our method complements previous work on synthetic turbulence, and yields physically plausible results. In addition, we demonstrate that our method can efficiently compute turbulent flows around a variety of objects including cars, whisks, as well as boulders in a river flow. We can even apply our model to precomputed static flow fields, yielding turbulent dynamics without a costly simulation.

Deforming Meshes that Split and Merge   [CITE]
Chris Wojtan, Nils Thuerey, Markus Gross and Greg Turk
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2009
Paper: [PDF]    Video: [MOV]    Video: [WWW]   

Abstract: We present a method for accurately tracking the moving surface of deformable materials in a manner that gracefully handles topological changes. We employ a Lagrangian surface tracking method, and we utilize a triangle mesh for our surface representation so that fine features can be retained. We make topological changes to the mesh by first identifying merging or splitting events at a particular grid resolution, and then locally creating new pieces of the mesh in the affected cells using a standard isosurface creation method. We stitch the new, topologically simplified portion of the mesh to the rest of the mesh at the cell boundaries. Our method detects and treats topological events with an emphasis on the preservation of thin features, while simultaneously simplifying those portions of the material that are not visible. Our surface tracker is not tied to a particular method for simulating deformable materials. In particular, we show results from two significantly different simulators: a Lagrangian FEM simulator with tetrahedral elements, and an Eulerian grid-based fluid simulator. Although our surface tracking method is generic, it is particularly well-suited for simulations that exhibit thin features and numerous topological events. Highlights of our results include a taffy-pulling animation with many fold and merge events, the creation and separation of thin strands in the simulation of viscoelastic materials, and the retention of thin sheets and surface details in a splashing fluid animation.

Visibility Transition Planning for Real-Time Camera Control   [CITE]
Thomas Oskam, Robert W. Sumner, Nils Thuerey, Markus Gross
SCA '09: Proceedings of the 2009 ACM SIGGRAPH/Eurographics Symposium on Computer Animation; Eurographics Association, 2009
Paper: [PDF]    Video: [MOV]   

Abstract: We present a real-time camera control system that uses a global planning algorithm to compute large, occlusion free camera paths through complex environments. The algorithm incorporates the visibility of a focus point into the search strategy, so that a path is chosen along which the focus target will be in view. The efficiency of our algorithm comes from a visibility-aware roadmap data structure that permits the precomputation of a coarse representation of all collision-free paths through an environment, together with an estimate of the pair-wise visibility between all portions of the scene. Our runtime system executes a path planning algorithm using the precomputed roadmap values to find a coarse path, and then refines the path using a sequence of occlusion maps computed on-the-fly. An iterative smoothing algorithm, together with a physically-based camera model, ensures that the path followed by the camera is smooth in both space and time. Our global planning strategy on the visibility-aware roadmap enables large-scale camera transitions as well as a local third-person camera module that follows a player and avoids obstructed viewpoints. The data structure itself adapts at run-time to dynamic occluders that move in an environment. We demonstrate these capabilities in several realistic game environments.

Detail-Preserving Fluid Control   [CITE]
Nils Thuerey, Richard Keiser, Ulrich Ruede, Mark Pauly
Graphical Models; Elsevier, 2009
Paper-Online: [WWW]    Video: [AVI]    Video: [WWW]   

Abstract: We propose a new fluid control technique that uses scale-dependent force control to preserve small-scale fluid detail. Control particles define local force fields and can be generated automatically from either a physical simulation or a sequence of target shapes. We use a multi-scale decomposition of the velocity field and apply control forces only to the coarse-scale components of the flow. Small-scale detail is thus preserved in a natural way avoiding the artificial viscosity often introduced by force-based control methods. We demonstrate the effectiveness of our method for both Lagrangian and Eulerian fluid simulation environments.

The results of this paper can also be seen in the "Magic Fluid Control" animation, which can be found below.

- 2008 -


Wavelet Turbulence for Fluid Simulation   [CITE]
Theodore Kim, Nils Thuerey, Doug James, Markus Gross
ACM Transactions on Graphics (SIGGRAPH); ACM Press, 2008
Webpage: [WWW]    Source-code: [WWW]    Paper: [PDF]    Video: [WWW]   

Abstract: We present a novel wavelet method for the simulation of fluids at high spatial resolution. The algorithm enables large- and small-scale detail to be edited separately, allowing high-resolution detail to be added as a post-processing step. Instead of solving the Navier-Stokes equations over a highly refined mesh, we use the wavelet decomposition of a low-resolution simulation to determine the location and energy characteristics of missing high-frequency components. We then synthesize these missing components using a novel incompressible turbulence function, and provide a method to maintain the temporal coherence of the resulting structures. There is no linear system to solve, so the method parallelizes trivially and requires only a few auxiliary arrays. The method guarantees that the new frequencies will not interfere with existing frequencies, allowing animators to set up a low resolution simulation quickly and later add details without changing the overall fluid motion.

Robust and Efficient Wave Simulations on Deforming Meshes   [CITE]
Roland Angst, Nils Thuerey, Mario Botsch, Markus Gross
Computer Graphics Forum; Blackwell Publishing, 2008
Paper: [PDF]    Video: [MOV]   

Abstract: The goal of this paper is to enable the interactive simulation of phenomena such as animated fluid characters. While full 3D fluid solvers achieve this with control algorithms, these 3D simulations are usually too costly for real-time environments. In order to achieve our goal, we reduce the problem from a three- to a two-dimensional one, and make use of the shallow water equations to simulate surface waves that can be solved very efficiently. In addition to a low runtime cost, stability is likewise crucial for interactive applications. Hence, we make use of an implicit time integration scheme to obtain a robust solver. To ensure a low energy dissipation, we apply an Implicit Newmark time integration scheme. We propose a general formulation of the underlying equations that is tailored towards the use with an Implicit Newmark integrator. Furthermore, we gain efficiency by making use of a direct solver. Due to the generality of our formulation, the fluid simulation can be coupled interactively with arbitrary external forces, such as forces caused by inertia or collisions. We will discuss the properties of our algorithm, and demonstrate its robustness with simulations on strongly deforming meshes.

Stable free surface flows with the lattice Boltzmann method on adaptively coarsened grids   [CITE]
Nils Thuerey, Ulrich Ruede
Computing and Visualization in Science; Springer, 2009
Online: [WWW]    Preprint: [PDF]    Video: [AVI] [MPG]   

Abstract: In this paper we will present an algorithm to perform free surface flow simulations with the lattice Boltzmann method (LBM) on adaptive grids. This reduces the required computational time by more than a factor of three for simulations with large volumes of fluid. To achieve this, the simulation of large fluid regions is performed with coarser grid resolutions.We have developed a set of rules to dynamically adapt the coarse regions to the movement of the free surface, while ensuring the consistency of all grids. Furthermore, the free surface treatment is combined with a Smagorinsky turbulence model and a technique for adaptive time steps to ensure stable simulations. The method is validated by comparing the position of the free surface with an uncoarsened simulation. It yields speedup factors of up to 3.85 for a simulation with a resolution of 480^3 cells and three coarser grid levels, and thus enables efficient and stable simulations of free surface flows, e.g. for highly detailed physically based animations of fluids.

The ETH Game Programming Laboratory: A Capstone for Computer Science and Visual Computing   [CITE]
Robert Sumner, Nils Thuerey and Markus Gross
Game Development in Computer Science Education (GDCSE); ACM, 2008
Paper: [PDF]    Course: [WWW]   

Abstract: The Visual Computing bachelors/masters program at ETH Zurich provides an internationally renowned degree in computer science with a specialization track in computer graphics. A new project-based game development course serves as a capstone to the program by reinforcing core computer science concepts and specialized topics in Visual Computing. Additionally, students learn design principles and obtain a better understanding of the interplay between the desires of game design and the realities of technical implementation. Finally, students practice crucial ┐soft skills┐ such as team work, effective communication, time management, and leadership. This article details the course goals and structure, presents three case studies of student-made games and the effect of the class on the students, and evaluates the overall class design. We hope that this document presents a compelling argument in favor of game development as a capstone to computer science and also provides useful insights for other academics wishing to incorporate game development into the computer science curriculum.

- 2007 -


Hybrid Parallelization Techniques for Lattice Boltzmann Free Surface Flows   [CITE]
Nils Thuerey, Thomas Pohl and Ulrich Ruede
Proceedings of Parallel CFD 2007; -, 2007
Paper: [PDF]   

Abstract: In the following, we will present an algorithm to perform adaptive free surface simulations with the lattice Boltzmann method (LBM) on machines with shared and distributed memory architectures. Performance results for different test cases and architectures will be given. The algorithm for parallelization yields a high performance, and can be combined with the adaptive LBM simulations. Moreover, the effects of the adaptive simulation on the parallel performance will be evaluated.

Simulation of moving particles in 3D with the Lattice Boltzmann method   [CITE]
Klaus Iglberger, Nils Thuerey and Ulrich Ruede
Computers and Mathematics with Applications, Mesoscopic Methods in Engineering and Science; Elsevier, 2008
Webpage: [WWW]   

Abstract: The paper presents a way to simulate the behavior of particle agglomerates in a fluid flow by coupling the Lattice Boltzmann Method to a rigid body physics engine. By extending the basic algorithm by a fluid/particle force interaction method, the hydrodynamic forces acting on the particles can be calculated. By the use of this force interaction between the fluid and the particles and by the use of the rigid body physics engine, the movement and collision behavior of particles in a flow can be simulated. Additionally, this coupled simulation system is able to simulate the internal particle forces in the connections between sintered particles, which could break due to the forces and torques of a shear flow. This permits a prediction of possible break-ups or structural displacements.

Real-time Breaking Waves for Shallow Water Simulations   [CITE]
Nils Thuerey, Matthias Mueller-Fischer, Simon Schirm, Markus Gross
Proceedings of the Pacific Conference on Computer Graphics and Applications 2007; IEEE Computer Society, 2007
Paper: [PDF]    Video: [MOV]   

Abstract: We present a new method for enhancing shallow water simulations by the effect of overturning waves. While full 3D fluid simulations can capture the process of wave breaking, this is beyond the capabilities of a pure height field model. 3D simulations, however, are still too expensive for real-time applications, especially when large bodies of water need to be simulated. The extension we propose overcomes this problem and makes it possible to simulate scenes such as waves near a beach, and surf riding characters in real-time. In a first step, steep wave fronts in the height field are detected and marked by line segments. These segments then spawn sheets of fluid represented by connected particles. When the sheets impinge on the water surface, they are absorbed and result in the creation of particles representing drops and foam. To enable interesting applications, we furthermore present a two-way coupling of rigid bodies with the fluid simulation. The capabilities and efficiency of the method will be demonstrated with several scenes, which run in real-time on today┐s commodity hardware.

Real-time simulations of bubbles and foam within a shallow water framework   [CITE]
Nils Thuerey, Filip Sadlo, Simon Schirm, Matthias Mueller-Fischer and Markus Gross
SCA '07: Proceedings of the 2007 ACM SIGGRAPH/Eurographics Symposium on Computer Animation; Eurographics Association, 2007
Paper: [PDF]    Video: [MOV]   

Abstract: Bubbles and foam are important fluid phenomena on scales that we encounter in our lives every day. While different techniques to handle these effects were developed in the past years, they require a full 3D fluid solver with free surfaces and surface tension. We present a shallow water based particle model that is coupled with a smoothed particle hydrodynamics simulation to demonstrate that real-time simulations of bubble and foam effects are possible with high frame rates. A shallow water simulation is used to represent the overall water volume. It is coupled to a particle-based bubble simulation with a flow field of spherical vortices. This bubble simulation is interacting with a smoothed particle hydrodynamics simulation including surface tension to handle foam on the fluid surface. The realism and performance of our approach is demonstrated with several test cases that run with high frame rates on a standard PC.

- 2006 -


Detail-Preserving Fluid Control   [CITE]
Nils Thuerey, Richard Keiser, Ulrich Ruede and Mark Pauly
SCA '06: Proceedings of the 2006 ACM SIGGRAPH/Eurographics Symposium on Computer Animation; Eurographics Association, 2006
Paper: [PDF]    Video: [AVI]   

Abstract: We propose a new fluid control technique that uses scale-dependent force control to preserve small-scale fluid detail. Control particles define local force fields and can be generated automatically from either a physical simulation or a sequence of target shapes. We use a multi-scale decomposition of the velocity field and apply control forces only to the coarse-scale components of the flow. Small-scale detail is thus preserved in a natural way avoiding the artificial viscosity often introduced by force-based control methods. We demonstrate the effectiveness of our method for both Lagrangian and Eulerian fluid simulation environments.

The results of this paper can also be seen in the "Magic Fluid Control" animation, which can be found here.

Animation of Open water Phenomena with coupled Shallow Water and Free Surface Simulation   [CITE]
Nils Thuerey, U. Ruede and M. Stamminger
SCA '06: Proceedings of the 2006 ACM SIGGRAPH/Eurographics Symposium on Computer Animation; Eurographics Association, 2006
Paper: [PDF]    Video: [AVI]   

Abstract: The goal of this paper is to perform simulations that capture fluid effects from small drops up to the propagation of large waves. To achieve this, we present a hybrid simulation method, that couples a two-dimensional shallow water simulation with a full three-dimensional free surface fluid simulation. We explain the approximations imposed by the shallow water model, and how to parametrize it according to the parameters of a 3D simulation. Each simulation is used to initialize double layered boundary conditions for the other one. The area covered by the 2D region can be an order of magnitude larger than the 3D region without significantly effecting the overall computation time. The 3D region can furthermore be easily moved within the 2D region during the course of the simulation. To achieve realistic results we combine our simulation method with a physically based model to generate and animate drops. For their generation we make use of the fluid turbulence model, and animate them with a simplified drag calculation. This allows simulations with relatively low resolutions.

Free Surface Flows with Moving and Deforming Objects for LBM   [CITE]
Nils Thuerey, K. Iglberger and U. Ruede
Proceedings of Vision, Modeling and Visualization 2006; IOS Press, 2006
Paper: [PDF]    Video: [MPG]   

Abstract: We present an algorithm to handle moving and deforming objects for free surface fluid simulations with the lattice Boltzmann method (LBM). To achieve this we extend methods available for flows without a free surface to enables simulations of moving objects with varying surface roughness, two-way coupled interaction and improved mass conservation. We furthermore show how to efficiently initialize boundary conditions for the moving objects from an arbitrary triangle mesh. The triangulation of the free surface can be improved by removing the interface between fluid and an immersed object. We will demonstrate the capabilities of our approach with two simpler test cases, and a more complicated simulation using an animated character as obstacle.

Enhanced Motion Blur Calculation with Optical Flow   [CITE]
Yuan Zheng, Harald Koestler, Nils Thuerey and Ulrich Ruede
Proceedings of Vision, Modeling and Visualization 2006; IOS Press, 2006
Paper: [PDF]   

Abstract: Computing motion blur is a complex task for modern raytracers, usually temporal supersampling is used to compute motion blur. However, this leads to increased rendering time and high computational complexity. In the following we present an image based technique which is based on the optical flow method to achieve an approximated motion field in an image sequence. In contrast to vector fields generated by the raytracer itself, this enables us to also attain motion vectors for indirect movements, such as shadows or caustics. For the optical flow computation we use an extended Horn and Schunk model based on the assumption that changes in the illumination are ignored. Experiments to compare the results of motion blurred pictures, computed by a raytracer and the optical flow method, will be presented and evaluated with respect to visual effects and runtime.

Simulation of the Hydrodynamic Drag of Aggregated Particles   [CITE]
Christian Binder, Christian Feichtinger, Hans-Joachim Schmid, Nils Thuerey, Wolfgang Peukert and Ulrich Ruede
Journal of Colloid and Interface Science; Elsevier, 2006
Paper online: [WWW]   

Abstract: The drag force on aggregates and partially sintered agglomerates is assessed using the lattice Boltzmann method (LBM) and accelerated Stokesian dynamics (ASD). Both methods have been compared in terms of accuracy and computational effort. It is shown that they give comparable results if all numerical parameters are controlled carefully. LBM requires a much higher computational effort, however, in contrast to ASD it is able to simulate partially sintered agglomerates as well. The results show that even a very small amount of sintering leads to a significant reduction in the drag force. The analysis of the drag force on agglomerates as well as on aggregates shows that there is no simple geometric quantity which is uniquely related to the drag force. Moreover, there is a significant variation in drag force for single aggregates at different orientations or for the orientation averaged drag force of different aggregates of the same size. This is explained by the structural effects which may lead to a variation in the drag force up to 20%.

Parallel Lattice Boltzmann Methods for CFD Applications   [CITE]
Carolin Koerner, Thomas Pohl, Ulrich Ruede, Nils Thuerey and Thomas Zeiser
Numerical Solution of Partial Differential Equations on Parallel Computers; Springer, 2006
Book Chaper: [PDF]    Book: [Table of Contents]   

Abstract: The lattice Boltzmann method (LBM) has evolved to a promising alternative to the well-established methods based on finite elements/volumes for computational fluid dynamics simulations. Ease of implementation, extensibility, and computational efficiency are the major reasons for LBM's growing field of application and increasing popularity. In this paper we give a brief introduction to the involved theory and equations for LBM, present various techniques to increase the single-CPU performance, outline the parallelization of a standard LBM implementation, and show performance results. In order to demonstrate the straightforward extensibility of LBM, we then focus on an application in material science involving fluid flows with free surfaces. We discuss the required extensions to handle this complex scenario, and the impact on the parallelization technique.

Optimization and Stabilization of LBM Free Surface Flow Simulations using Adaptive Parameterization   [CITE]
Nils Thuerey, Thomas Pohl, Ulrich Ruede, Markus Oechsner and Carolin Koerner
Computers and Fluids; Elsevier, 2006
Webpage: [WWW]    Preprint: [PDF]   

Abstract: We present a method to speed up and stabilize free surface simulations with the lattice Boltzmann method (LBM). This is done by adaptively changing the parameterization of the simulation in a way that corresponds to a different size of the simulation time step. This means that the Mach number changes as well, and requires a rescaling of all distribution functions. Hence we only perform the rescaling when the velocities in the simulation become too large or small. We will demonstrate the effect of this procedure for two and three-dimensional test cases. In addition to a reduction of the necessary LBM steps, this method can also be used to stabilize gravity driven simulations, where the maximum velocities are not known a priori.

- 2005 -


Lattice Boltzmann Model for Free Surface Flow for Modeling Foaming   [CITE]
Carolin Koerner, Michael Thies, Thomas Hofmann, Nils Thuerey and Ulrich Ruede
Journal of Statistical Physics; Springer, 2005
Paper: [PDF]   

Abstract: We present a 2D- and 3D-lattice Boltzmann model for the treatment of free surface flows including gas diffusion. Interface advection and related boundary conditions are based on the idea of the lattice Boltzmann equation. The fluid dynamic boundary conditions are approximated by using the mass and momentum fluxes across the interface, which do not require explicit calculation of gradients. A similar procedure is applied to fulfill the diffusion boundary condition. Simple verification tests demonstrate the correctness of the algorithms. 2D- and 3D-foam evolution examples demonstrate the potential of the method.

- 2004 -


Free Surface Lattice-Boltzmann fluid simulations with and without level sets   [CITE]
Nils Thuerey and U. Ruede
Proc. of Vision, Modelling, and Visualization VMV; IOS Press, 2004
Paper: [PDF]    Poster: [PDF]   

Abstract: We present two variants of free surface Lattice-Boltzmann fluid simulations for the animation of liquids in computer graphics. The Lattice-Boltzmann method is an attractive alternative to conventional fluid solvers, due to its simplicity and flexibility, especially for changing geometries and topologies. While our first method directly calculates the mass fluxes between the cells of the computational grid, another variant of the method is explained, that uses level sets to track the fluid surface. This has advantages for the smoothness of the fluid surface and improves the representation of details in the free surface, but makes the conservation of mass more difficult. Several examples will be shown to highlight the differences between the two methods.

Performance Evaluation of Parallel Large-Scale Lattice Boltzmann Applications on Three Supercomputing Architectures   [CITE]
Thomas Pohl, Frank Deserno, Nils Thuerey, Ulrich Ruede, Peter Lammers, Gerhard Wellein and Thomas Zeiser
SC '04: Proceedings of the 2004 ACM/IEEE conference on Supercomputing; IEEE Computer Society, 2004
Paper: [PDF]    TechReport: [PDF] [PS.GZ]   

Abstract: Computationally intensive programs with moderate communication requirements such as CFD codes suffer from the standard slow interconnects of commodity off the shelf (COTS) hardware. We will introduce different large-scale applications of the Lattice Boltzmann Method (LBM) in fluid dynamics, material science, and chemical engineering and present results of the parallel performance on different architectures. It will be shown that a high speed communication network in combination with an efficient CPU is mandatory in order to achieve the required performance. An estimation of the necessary CPU count to meet the performance of 1 TFlop/s will be given as well as a prediction as to which architecture is the most suitable for LBM. Finally, ratios of costs to application performance for tailored HPC systems and COTS architectures will be presented.

- 2002 -


Performance Optimization of 3D Multigrid on Hierarchical Memory Architectures   [CITE]
Markus Kowarschik, Ulrich Ruede, Nils Thuerey and Christian Weiss
Proceedings of PARA'02; Springer, Lecture Notes in Computer Science, 2002
Paper: [PDF]   
Theses / Articles:
[Back to top...]

- 2011 -


The Design and Development of Computer Games   [CITE]
Markus Gross, Robert Sumner, Nils Thuerey
The Design of Material, Organism, and Minds (Editors: S. Lang, M. Hampe); Springer, 2011
Online: [WWW]    Preprint: [PDF]   

Abstract: The design of modern computer games is as much an art as painting, sculpting, music, or writing. Albeit relatively young and still evolving, game design is highly complex and requires a broad spectrum of artistic and technical skills. The following contribution reviews the process of designing and developing modern computer games. We will acquaint the reader with the most important fundamentals of game design, walk through the various stages of the design process, and highlight the specifics of each stage.

- 2008 -


Real-time Physics, SIGGRAPH 2008 Course
Matthias Mueller-Fischer, Doug James, Jos Stam, Nils Thuerey
Webpage: [WWW]    SWS-source-code: [TGZ]    LBM-source-code: [TGZ]    Course-notes: [PDF]   

Abstract: Physical simulations have become an important component of computer games. In next-generation games, players expect to see fully dynamic and destructible worlds, and this requires fast and stable simulation methods. In this class, lecturers who have made significant contributions in simulation methods present a wide spectrum of state-of-the-art methods for real-time simulation of rigid and deformable solids, and smoke and liquid simulation. In addition to the underlying physical equations, they present practical simulation methods and algorithms that will help physical-simulation developers and game developers apply these techniques properly.

Magic Fluid Control, Animation
Nils Thuerey, Mark Pauly, Richard Keiser, Ulrich Ruede
Video: [MPG]    Webpage: [WWW]    Video (uploaded by someone else): [WWW]   

Abstract: The movie consists of three clips to demonstrate the possibilities of controlled water simulations: 1. A magician pulls out water from a basin which forms a teapot by moving his hand upwards from the basin. 2. Water flows magically upwards a stair where it forms a human figure. This is achieved by using control particles from a reversed water simulation and from an invisible model. 3. The magician moves a duck out of water from one basin to another and vice versa. Finally, he pulls out two ducks from both basins which end up in a short "fight".

So far, the animation was shown at the following events:
Eurographics 2006, Animations & Interactive Media
Computer Graphics International 2007, Video Festival
SIGGRAPH 2007, Computer Animation Festival
Melbourne International Animation Festival 2008

- 2007 -


Physically based Animation of Free Surface Flows with the Lattice Boltzmann Method   [CITE]
N. Thuerey
PhD thesis; Dept. of Computer Science 10, University of Erlangen-Nuremberg, 2007
Thesis: [PDF]   

Abstract: The numerical simulation of fluids has become an established tool in many engineering applications. Free surface fluids represent a special case that is important for a variety of applications. For a free surface simulation, a two phase system, such as air and water, is described by a single fluid phase with a sharp interface and corresponding boundary conditions. This allows the efficient representation and simulation of complex problems. In this thesis, the main application for free surface flows will be the generation of animations of liquids. Additionally, engineering applications from material science and particle technology are considered.

- 2006 -


Fluid Simulation with Blender   [CITE]
Nils Thuerey
Dr. Dobbs Journal; CMP Media, 2006
Webpage: [WWW]   

Abstract: The goal of the Fluid Simulation with Blender project was the integration of the free surface fluid solver developed at the system simulation group (LSS) in Erlangen, Germany into the Blender graphics toolset (3D modeling, animation, rendering, and the like) and to publish it as open-source software.

- 2005 -


Simulation of moving Nano-Particles with the Lattice Boltzmann Method in 3D   [CITE]
Iglberger, Nils Thuerey, Ulrich Ruede, H.-J. Schmid and W. Peukert
Proceedings of ASIM '05; SCS Publishing House, 2005
Paper: [PDF]   

Abstract: The following work presents a way to simulate the nano-particle behavior in a flow by coupling the Lattice Boltzmann Method to a rigid body physics engine. By extending the basic fluid simulation for the treatment of curved particle surfaces and by a force interaction method, the fluid forces acting on the nano-particles can be calculated.With this force interaction between the fluid and the particles and the use of the rigid body physics engine, the movement and collision behavior of nano-particles in a flow can be simulated. Additionally, this coupled simulation system is able to calculate the internal particle forces in the connections between sintered particles, which could break due to the forces of a shear flow. Therefore a prediction of possible break-ups becomes possible.

Drag Force Simulations of Particle Agglomerates with the Lattice-Boltzmann Method   [CITE]
Christian Feichtinger, Nils Thuerey, Ulrich Ruede, Christian Binder, Hans Schmid, Wolfgang Peukert
Proceedings of ASIM '05; SCS Publishing House, 2005
Paper: [PDF]   

Abstract: In this work we show, that the Lattice Boltzmann method is capable of calculating the drag force on spheres with a sufficient accuracy. The main interest lies on laminar flows around nano particle agglomerates. To give an overview over the simulated force accuracy the results are compared with analytical solutions. Thereby the benefit from curved boundary treatments such as [YD03] and [PL03] is taken into account. Through the use of curved boundary treatments the lattice radius of the spheres can be reduced and thus smaller domain sizes can be used with the same force accuracy. Further more, results for complex agglomerates are presented and compared to solutions of the Stokesian Dynamics approach [AS01].

Optimized Free Surface Fluids on Adaptive Grids with the Lattice Boltzmann Method   [CITE]
Nils Thuerey and Ulrich Ruede
Poster; SIGGRAPH '05, 2005
Video: [AVI]    Poster: [PDF]    Abstract: [PDF]   

Abstract: The algorithm used to simulate the free surface fluid is based on the Lattice Boltzmann Method (LBM). The free surface tracking is similar to Volume-of-Fluid methods for conventional Navier-Stokes solvers. For the free surface three types of cells are distinguished - fluid cells that are treated with the normal LBM, empty cells, that do not require any computations, and interface cells. For these the amount of fluid mass is tracked for each time step, together with appropriate boundary conditions for the fluid gas interface. The algorithm is described in more detail in e.g. [Thuerey 2003]. To speed up the simulations we use a modied version of the grid-refinement described in [Filippova and Haenel 1998] together with a method to adapt the grid boundaries according to the free surface movement.

Interactive Free Surface Fluids with the Lattice Boltzmann Method   [CITE]
Nils Thuerey, Carolin Koerner and Ulrich Ruede
Technical Report 05-4; Department of Computer Science 10 System Simulation, 2005
Report: [PDF]    Video: [DivX] [MPG]   

Abstract: In this paper we present our algorithm for animating fluids with a free surface. It is based on the Lattice-Boltzmann Method, instead of a direct discretization of the Navier-Stokes equations. This allows a relatively simple treatment of the free surface boundary conditions at high computational efficiency, without sacrificing the underlying physics. We give a detailed description of our algorithm, focussing on details that are required to achieve a good visual appearance. Furthermore we describe how to implement our adaptive time stepping technique to achieve flexible and stable simulations. We will demonstrate the speed and capabilities of the method with animations from different interactive test cases. These run with, on average, more than 20 frames per second on a standard desktop PC.

- 2003 -


Simulation von Metallschaum mittels der Lattice-Boltzmann Methode   [CITE]
Nils Thuerey, Ulrich Ruede and Carolin Koerner
Konwihr Quartl; KONWIHR, 2003
Paper: [PDF]   

A Lattice Boltzmann method for single-phase free surface flows in 3D   [CITE]
Nils Thuerey
Master thesis; Dept. of Computer Science 10, University of Erlangen-Nuremberg, 2003
Thesis: [PDF] [PS.GZ]   

Abstract: The fluid solver presented here is capable of simulating a single fluid phase with a free surface, including surface tension, bubbles and coalescence. In contrast to the standard multi-phase models of [Gunstensen et al., 1991] and [Swift et al., 1996] it is not necessary to simulate the motion of the gas phase, which overcomes restrictions in difference of the viscosity of the two phases, and improves the computational time needed for the numerical simulation, as the gas volumes in the foam do not require any additional computations. The forces at the fluid interface, gas pressure and surface tension, are applied by reconstructing the missing information from the gas phase. For surface tension, points on the fluid surface are calculated with the marching cubes algorithm and used to retrieve an average curvature at the fluid-gas interface. As part of this thesis, the fluid solver is tested and validated with standard problems like the breaking dam problem and rising bubble problems.

- 2002 -


Cache Optimizations for Multigrid in 3D   [CITE]
Nils Thuerey
Semester thesis; Dept. of Computer Science 10, University of Erlangen-Nuremberg, 2002
Thesis: [PDF]