header

Profile


photo

Tassilo Kugelstadt, M.Sc.
Email: kugelstadt@cs.rwth-aachen.de



Publications


Implicit Density Projection for Volume Conserving Liquids


Tassilo Kugelstadt, Andreas Longva, Nils Thuerey, Jan Bender
IEEE Transactions on Visualization and Computer Graphics
pubimg

We propose a novel implicit density projection approach for hybrid Eulerian/Lagrangian methods like FLIP and APIC to enforce volume conservation of incompressible liquids. Our approach is able to robustly recover from highly degenerate configurations and incorporates volume-conserving boundary handling. A problem of the standard divergence-free pressure solver is that it only has a differential view on density changes. Numerical volume errors, which occur due to large time steps and the limited accuracy of pressure projections, are invisible to the solver and cannot be corrected. Moreover, these errors accumulate over time and can lead to drastic volume changes, especially in long-running simulations or interactive scenarios. Therefore, we introduce a novel method that enforces constant density throughout the fluid. The density itself is tracked via the particles of the hybrid Eulerian/Lagrangian simulation algorithm. To achieve constant density, we use the continuous mass conservation law to derive a pressure Poisson equation which also takes density deviations into account. It can be discretized with standard approaches and easily implemented into existing code by extending the regular pressure solver. Our method enables us to relax the strict time step and solver accuracy requirements of a regular solver, leading to significantly higher performance. Moreover, our approach is able to push fluid particles out of solid obstacles without losing volume and generates more uniform particle distributions, which makes frequent particle resampling unnecessary. We compare the proposed method to standard FLIP and APIC and to previous volume correction approaches in several simulations and demonstrate significant improvements in terms of incompressibility, visual realism and computational performance.

» Show Videos
» Show BibTeX

@Article{BKKW21,
author = {Tassilo Kugelstadt and Andreas Longva and Nils Thuerey and Jan Bender},
title = {Implicit Density Projection for Volume Conserving Liquids},
journal = {IEEE Transactions on Visualization and Computer Graphics},
year = {2021},
publisher = {IEEE},
volume = {27},
number = {4},
doi={ 10.1109/TVCG.2019.2947437},
}





Accurately Solving Rod Dynamics with Graph Learning


Han Shao, Tassilo Kugelstadt, Torsten Hädrich, Wojciech Palubicki, Jan Bender, Sören Pirk, Dominik L. Michels
Advances in Neural Information Processing Systems (NeurIPS 2021)
pubimg

Iterative solvers are widely used to accurately simulate physical systems. These solvers require initial guesses to generate a sequence of improving approximate solutions. In this contribution, we introduce a novel method to accelerate iterative solvers for rod dynamics with graph networks (GNs) by predicting the initial guesses to reduce the number of iterations. Unlike existing methods that aim to learn physical systems in an end-to-end manner, our approach guarantees long-term stability and therefore leads to more accurate solutions. Furthermore, our method improves the run time performance of traditional iterative solvers for rod dynamics. To explore our method we make use of position-based dynamics (PBD) as a common solver for physical systems and evaluate it by simulating the dynamics of elastic rods. Our approach is able to generalize across different initial conditions, discretizations, and realistic material properties. We demonstrate that it also performs well when taking discontinuous effects into account such as collisions between individual rods. Finally, to illustrate the scalability of our approach, we simulate complex 3D tree models composed of over a thousand individual branch segments swaying in wind fields.

» Show Videos
» Show BibTeX

@inproceedings{Shao:2021:GraphLearning,
title={Accurately Solving Rod Dynamics with Graph Learning},
author={Han Shao and Tassilo Kugelstadt and Torsten H\"{a}drich and Wojciech Pa\l{}ubicki and Jan Bender and S\"{o}ren Pirk and Dominik L. Michels},
year={2021},
booktitle={Advances in Neural Information Processing Systems (NeurIPS)},
URL={http://computationalsciences.org/publications/shao-2021-physical-systems-graph-learning.html}
}





Fast Corotated Elastic SPH Solids with Implicit Zero-Energy Mode Control


Tassilo Kugelstadt, Jan Bender, José Antonio Fernández-Fernández, Stefan Rhys Jeske, Fabian Löschner, Andreas Longva
Proceedings of the ACM on Computer Graphics and Interactive Techniques (Best Paper Award at SCA)
pubimg

We develop a new operator splitting formulation for the simulation of corotated linearly elastic solids with Smoothed Particle Hydrodynamics (SPH). Based on the technique of Kugelstadt et al. [KKB2018] originally developed for the Finite Element Method (FEM), we split the elastic energy into two separate terms corresponding to stretching and volume conservation, and based on this principle, we design a splitting scheme compatible with SPH. The operator splitting scheme enables us to treat the two terms separately, and because the stretching forces lead to a stiffness matrix that is constant in time, we are able to prefactor the system matrix for the implicit integration step. Solid-solid contact and fluid-solid interaction is achieved through a unified pressure solve. We demonstrate more than an order of magnitude improvement in computation time compared to a state-of-the-art SPH simulator for elastic solids.

We further improve the stability and reliability of the simulation through several additional contributions. We introduce a new implicit penalty mechanism that suppresses zero-energy modes inherent in the SPH formulation for elastic solids, and present a new, physics-inspired sampling algorithm for generating high-quality particle distributions for the rest shape of an elastic solid. We finally also devise an efficient method for interpolating vertex positions of a high-resolution surface mesh based on the SPH particle positions for use in high-fidelity visualization.

» Show Videos
» Show BibTeX

@article{KBF+21,
author = {Kugelstadt, Tassilo and Bender, Jan and Fern{\'{a}}ndez-Fern{\'{a}}ndez, Jos{\'{e}} Antonio and Jeske, Stefan Rhys and L{\"{o}}schner, Fabian and Longva, Andreas},
title = {Fast Corotated Elastic SPH Solids with Implicit Zero-Energy Mode Control},
year = {2021},
issue_date = {September 2021},
publisher = {Association for Computing Machinery},
address = {New York, NY, USA},
volume = {4},
number = {3},
url = {https://doi.org/10.1145/3480142},
doi = {10.1145/3480142},
journal = {Proc. ACM Comput. Graph. Interact. Tech.},
month = sep,
articleno = {33},
numpages = {21},
keywords = {Smoothed Particle Hydrodynamics, fluid simulation, deformable solids, solid-fluid coupling}
}





Higher-Order Finite Elements for Embedded Simulation


Andreas Longva, Fabian Löschner, Tassilo Kugelstadt, José Antonio Fernández-Fernández, Jan Bender
ACM Transactions on Graphics (SIGGRAPH Asia 2020)
pubimg

As demands for high-fidelity physics-based animations increase, the need for accurate methods for simulating deformable solids grows. While higher-order finite elements are commonplace in engineering due to their superior approximation properties for many problems, they have gained little traction in the computer graphics community. This may partially be explained by the need for finite element meshes to approximate the highly complex geometry of models used in graphics applications. Due to the additional per-element computational expense of higher-order elements, larger elements are needed, and the error incurred due to the geometry mismatch eradicates the benefits of higher-order discretizations. One solution to this problem is the embedding of the geometry into a coarser finite element mesh. However, to date there is no adequate, practical computational framework that permits the accurate embedding into higher-order elements.

We develop a novel, robust quadrature generation method that generates theoretically guaranteed high-quality sub-cell integration rules of arbitrary polynomial accuracy. The number of quadrature points generated is bounded only by the desired degree of the polynomial, independent of the embedded geometry. Additionally, we build on recent work in the Finite Cell Method (FCM) community so as to tackle the severe ill-conditioning caused by partially filled elements by adapting an Additive-Schwarz-based preconditioner so that it is suitable for use with state-of-the-art non-linear material models from the graphics literature. Together these two contributions constitute a general-purpose framework for embedded simulation with higher-order finite elements.

We finally demonstrate the benefits of our framework in several scenarios, in which second-order hexahedra and tetrahedra clearly outperform their first-order counterparts.

» Show Videos
» Show BibTeX

@ARTICLE{ LLKFB20,
author= {Andreas Longva and Fabian L{\"{o}}schner and Tassilo Kugelstadt and Jos{\'{e}} Antonio Fern{\'{a}}ndez-Fern{\'{a}}ndez and Jan Bender },
title= {{Higher-Order Finite Elements for Embedded Simulation}},
year= {2020},
journal= {ACM Transactions on Graphics (SIGGRAPH Asia)},
publisher= {ACM},
volume = {39},
number = {6},
pages= {14}
}





Implicit Frictional Boundary Handling for SPH


Jan Bender, Tassilo Kugelstadt, Marcel Weiler, Dan Koschier
IEEE Transactions on Visualization and Computer Graphics
pubimg

In this paper, we present a novel method for the robust handling of static and dynamic rigid boundaries in Smoothed Particle Hydrodynamics (SPH) simulations. We build upon the ideas of the density maps approach which has been introduced recently by Koschier and Bender. They precompute the density contributions of solid boundaries and store them on a spatial grid which can be efficiently queried during runtime. This alleviates the problems of commonly used boundary particles, like bumpy surfaces and inaccurate pressure forces near boundaries. Our method is based on a similar concept but we precompute the volume contribution of the boundary geometry. This maintains all benefits of density maps but offers a variety of advantages which are demonstrated in several experiments. Firstly, in contrast to the density maps method we can compute derivatives in the standard SPH manner by differentiating the kernel function. This results in smooth pressure forces, even for lower map resolutions, such that precomputation times and memory requirements are reduced by more than two orders of magnitude compared to density maps. Furthermore, this directly fits into the SPH concept so that volume maps can be seamlessly combined with existing SPH methods. Finally, the kernel function is not baked into the map such that the same volume map can be used with different kernels. This is especially useful when we want to incorporate common surface tension or viscosity methods that use different kernels than the fluid simulation.

» Show Videos
» Show BibTeX

@Article{BKWK2020,
author = {Jan Bender and Tassilo Kugelstadt and Marcel Weiler and Dan Koschier },
title = {Implicit Frictional Boundary Handling for SPH},
journal = {IEEE Transactions on Visualization and Computer Graphics},
year = {2020},
publisher = {IEEE},
volume={26},
number={10},
pages={2982-2993},
doi={10.1109/TVCG.2020.3004245},
}





Higher-Order Time Integration for Deformable Solids


Fabian Löschner, Andreas Longva, Stefan Rhys Jeske, Tassilo Kugelstadt, Jan Bender
Computer Graphics Forum
pubimg

Visually appealing and vivid simulations of deformable solids represent an important aspect of physically based computer animation. For the temporal discretization, it is customary in computer animation to use first-order accurate integration methods, such as Backward Euler, due to their simplicity and robustness. Although there is notable research on second-order methods, their use is not widespread. Many of these well-known methods have significant drawbacks such as severe numerical damping or scene-dependent time step restrictions to ensure stability. In this paper, we discuss the most relevant requirements on such methods in computer animation and motivate the interest beyond first-order accuracy. Keeping these requirements in mind, we investigate several promising methods from the families of diagonally implicit Runge-Kutta (DIRK) and Rosenbrock methods which currently do not appear to have considerable popularity in this field. We show that the usage of such methods improves the visual quality of physical animations. In addition, we demonstrate that they allow distinctly more control over damping at lower computational cost than classical methods. As part of our theoretical contribution, we review aspects of simulations that are often considered more intricate with higher-order methods, such as contact handling. To this end, we derive an implicit linearized contact model based on a predictor-corrector approach that leads to consistent behavior with higher-order integrators as predictors. Our contact model is well suited for the simulation of stiff, nonlinear materials with the integration methods presented in this paper and more common methods such as Backward Euler alike.

» Show Videos
» Show BibTeX

@article{LLJKB20,
author = {Fabian L{\"{o}}schner and Andreas Longva and Stefan Jeske and Tassilo Kugelstadt and Jan Bender},
title = {Higher-Order Time Integration for Deformable Solids},
year = {2020},
journal = {Computer Graphics Forum},
volume = {39},
number = {8}
}





Turbulent Micropolar SPH Fluids with Foam


Jan Bender, Dan Koschier, Tassilo Kugelstadt, Marcel Weiler
IEEE Transactions on Visualization and Computer Graphics
pubimg

In this paper we introduce a novel micropolar material model for the simulation of turbulent inviscid fluids. The governing equations are solved by using the concept of Smoothed Particle Hydrodynamics (SPH). As already investigated in previous works, SPH fluid simulations suffer from numerical diffusion which leads to a lower vorticity, a loss in turbulent details and finally in less realistic results. To solve this problem we propose a micropolar fluid model. The micropolar fluid model is a generalization of the classical Navier-Stokes equations, which are typically used in computer graphics to simulate fluids. In contrast to the classical Navier-Stokes model, micropolar fluids have a microstructure and therefore consider the rotational motion of fluid particles. In addition to the linear velocity field these fluids also have a field of microrotation which represents existing vortices and provides a source for new ones. However, classical micropolar materials are viscous and the translational and the rotational motion are coupled in a dissipative way. Since our goal is to simulate turbulent fluids, we introduce a novel modified micropolar material for inviscid fluids with a non-dissipative coupling. Our model can generate realistic turbulences, is linear and angular momentum conserving, can be easily integrated in existing SPH simulation methods and its computational overhead is negligible. Another important visual feature of turbulent liquids is foam. Therefore, we present a post-processing method which considers microrotation in the foam particle generation. It works completely automatic and requires only one user-defined parameter to control the amount of foam.

» Show Videos
» Show BibTeX

@Article{BKKW19,
author = {Bender, Jan and Koschier, Dan and Kugelstadt, Tassilo and Weiler, Marcel},
title = {Turbulent Micropolar SPH Fluids with Foam},
journal = {IEEE Transactions on Visualization and Computer Graphics},
year = {2019},
publisher = {IEEE},
volume={25},
number={6},
pages={2284-2295},
doi={10.1109/TVCG.2018.2832080},
ISSN={1077-2626},
month={June},
}





Volume Maps: An Implicit Boundary Representation for SPH


Jan Bender, Tassilo Kugelstadt, Marcel Weiler, Dan Koschier
ACM SIGGRAPH Motion, Interaction and Games (Best Paper Award)
pubimg

In this paper, we present a novel method for the robust handling of static and dynamic rigid boundaries in Smoothed Particle Hydrodynamics (SPH) simulations. We build upon the ideas of the density maps approach which has been introduced recently by Koschier and Bender. They precompute the density contributions of solid boundaries and store them on a spatial grid which can be efficiently queried during runtime. This alleviates the problems of commonly used boundary particles, like bumpy surfaces and inaccurate pressure forces near boundaries. Our method is based on a similar concept but we precompute the volume contribution of the boundary geometry and store it on a grid. This maintains all benefits of density maps but offers a variety of advantages which are demonstrated in several experiments. Firstly, in contrast to the density maps method we can compute derivatives in the standard SPH manner by differentiating the kernel function. This results in smooth pressure forces, even for lower map resolutions, such that precomputation times and memory requirements are reduced by more than two orders of magnitude compared to density maps. Furthermore, this directly fits into the SPH concept so that volume maps can be seamlessly combined with existing SPH methods. Finally, the kernel function is not baked into the map such that the same volume map can be used with different kernels. This is especially useful when we want to incorporate common surface tension or viscosity methods that use different kernels than the fluid simulation.

» Show Videos
» Show BibTeX

@inproceedings{Bender2019,
author = {Jan Bender and Tassilo Kugelstadt and Marcel Weiler and Dan Koschier},
title = {Volume Maps: An Implicit Boundary Representation for SPH},
booktitle = {Proceedings of ACM SIGGRAPH Conference on Motion, Interaction and Games},
series = {MIG '19},
year = {2019},
publisher = {ACM}
}





Fast Corotated FEM using Operator Splitting


Tassilo Kugelstadt, Dan Koschier, Jan Bender
Computer Graphics Forum (SCA)
pubimg

In this paper we present a novel operator splitting approach for corotated FEM simulations. The deformation energy of the corotated linear material model consists of two additive terms. The first term models stretching in the individual spatial directions and the second term describes resistance to volume changes. By formulating the backward Euler time integration scheme as an optimization problem, we show that the first term is invariant to rotations. This allows us to use an operator splitting approach and to solve both terms individually with different numerical methods. The stretching part is solved accurately with an optimization integrator, which can be done very efficiently because the system matrix is constant over time such that its Cholesky factorization can be precomputed. The volume term is solved approximately by using the compliant constraints method and Gauss-Seidel iterations. Further, we introduce the analytic polar decomposition which allows us to speed up the extraction of the rotational part of the deformation gradient and to recover inverted elements. Finally, this results in an extremely fast and robust simulation method with high visual quality that outperforms standard corotated FEMs by more than two orders of magnitude and even the fast but inaccurate PBD and shape matching methods by more than one order of magnitude without having their typical drawbacks. This enables a very efficient simulation of complex scenes containing more than a million elements.

» Show Videos
» Show BibTeX

@article{KKB2018,
author = {Tassilo Kugelstadt and Dan Koschier and Jan Bender},
title = {Fast Corotated FEM using Operator Splitting},
year = {2018},
journal = {Computer Graphics Forum (SCA)},
volume = {37},
number = {8}
}





Direct Position-Based Solver for Stiff Rods


Crispin Deul, Tassilo Kugelstadt, Marcel Weiler, Jan Bender
Computer Graphics Forum
pubimg

In this paper, we present a novel direct solver for the efficient simulation of stiff, inextensible elastic rods within the Position-Based Dynamics (PBD) framework. It is based on the XPBD algorithm, which extends PBD to simulate elastic objects with physically meaningful material parameters. XPBD approximates an implicit Euler integration and solves the system of non-linear equations using a non-linear Gauss-Seidel solver. However, this solver requires many iterations to converge for complex models and if convergence is not reached, the material becomes too soft. In contrast we use Newton iterations in combination with our direct solver to solve the non-linear equations which significantly improves convergence by solving all constraints of an acyclic structure (tree), simultaneously. Our solver only requires a few Newton iterations to achieve high stiffness and inextensibility. We model inextensible rods and trees using rigid segments connected by constraints. Bending and twisting constraints are derived from the well-established Cosserat model. The high performance of our solver is demonstrated in highly realistic simulations of rods consisting of multiple ten-thousand segments. In summary, our method allows the efficient simulation of stiff rods in the Position-Based Dynamics framework with a speedup of two orders of magnitude compared to the original XPBD approach.

» Show Videos
» Show BibTeX

@article{DKWB2018,
author = {Crispin Deul and Tassilo Kugelstadt and Marcel Weiler and Jan Bender},
title = {Direct Position-Based Solver for Stiff Rods},
year = {2018},
journal = {Computer Graphics Forum},
volume = {37},
number = {6},
pages = {313-324},
keywords = {physically based animation, animation, Computing methodologies → Physical simulation},
doi = {10.1111/cgf.13326},
url = {https://onlinelibrary.wiley.com/doi/abs/10.1111/cgf.13326},
eprint = {https://onlinelibrary.wiley.com/doi/pdf/10.1111/cgf.13326},
}





A Micropolar Material Model for Turbulent SPH Fluids


Jan Bender, Dan Koschier, Tassilo Kugelstadt, Marcel Weiler
ACM SIGGRAPH / EUROGRAPHICS Symposium on Computer Animation (Best Paper Award)
pubimg

In this paper we introduce a novel micropolar material model for the simulation of turbulent inviscid fluids. The governing equations are solved by using the concept of Smoothed Particle Hydrodynamics (SPH). As already investigated in previous works, SPH fluid simulations suffer from numerical diffusion which leads to a lower vorticity, a loss in turbulent details and finally in less realistic results. To solve this problem we propose a micropolar fluid model. The micropolar fluid model is a generalization of the classical Navier-Stokes equations, which are typically used in computer graphics to simulate fluids. In contrast to the classical Navier-Stokes model, micropolar fluids have a microstructure and therefore consider the rotational motion of fluid particles. In addition to the linear velocity field these fluids also have a field of microrotation which represents existing vortices and provides a source for new ones. However, classical micropolar materials are viscous and the translational and the rotational motion are coupled in a dissipative way. Since our goal is to simulate turbulent fluids, we introduce a novel modified micropolar material for inviscid fluids with a non-dissipative coupling. Our model can generate realistic turbulences, is linear and angular momentum conserving, can be easily integrated in existing SPH simulation methods and its computational overhead is negligible.

» Show Videos
» Show BibTeX

@INPROCEEDINGS{Bender2017,
author = {Jan Bender and Dan Koschier and Tassilo Kugelstadt and Marcel Weiler},
title = {A Micropolar Material Model for Turbulent SPH Fluids},
booktitle = {Proceedings of the 2017 ACM SIGGRAPH/Eurographics Symposium on Computer
Animation},
year = {2017},
publisher = {ACM}
}





Position and Orientation Based Cosserat Rods


Tassilo Kugelstadt, Elmar Schoemer
ACM SIGGRAPH / EUROGRAPHICS Symposium on Computer Animation
pubimg

We present a novel method to simulate bending and torsion of elastic rods within the position-based dynamics (PBD) framework. The main challenge is that torsion effects of Cosserat rods are described in terms of material frames which are attached to the centerline of the rod. But frames or orientations do not fit into the classical position-based dynamics formulation. To solve this problem we introduce new types of constraints to couple orientations which are represented by unit quaternions. For constraint projection quaternions are treated in the exact same way as positions. Unit length is enforced with an additional constraint. This allows us to use the strain measures form Cosserat theory directly as constraints in PBD. It leads to very simple algebraic expressions for the correction displacements which only contain quaternion products and additions. Our results show that our method is very robust and accurately produces the complex bending and torsion effects of rods. Due to its simplicity our method is very efficient and more than one order of magnitude faster than existing position-based rod simulation methods. It even achieves the same performance as position-based simulations without torsion effects.

» Show Videos
» Show BibTeX

@INPROCEEDINGS{Kugelstadt2016,
author = {Tassilo Kugelstadt and Elmar Schoemer},
title = {Position and Orientation Based Cosserat Rods},
booktitle = {Proceedings of the 2016 ACM SIGGRAPH/Eurographics Symposium on Computer
Animation},
year = {2016},
publisher = {Eurographics Association},
location = {Zurich, Switzerland}
}





Disclaimer Home Visual Computing institute RWTH Aachen University