Langevin Particle: A Self-Adaptive Lagrangian Primitive For Flow Simulation Enhancement

We develop a new Lagrangian primitive, named Langevin particle, to incorporate turbulent flow details in fluid simulation. A group of the particles are distributed inside the simulation domain based on a turbulence energy model with turbulence viscosity. A particle in particular moves obeying the generalized Langevin equation, a well-known stochastic differential equation that describes the particle’s motion as a random Markov process. The resultant particle trajectory shows self-adapted fluctuation in accordance to the turbulence energy, while following the global flow dynamics. We then feed back Langevin forces to the simulation based on the stochastic trajectory, which drive the flow with necessary turbulence. The new hybrid flow simulation method features nonrestricted particle evolution requiring minimal extra manipulation after initiation. The flow turbulence is easily controlled and the total computational overhead of enhancement is minimal based on typical fluid solvers.

Langevin Particle: A Self-Adaptive Lagrangian Primitive For Flow Simulation Enhancement

Posted in Uncategorized | Comments Off on Langevin Particle: A Self-Adaptive Lagrangian Primitive For Flow Simulation Enhancement

PhD Thesis

Cem Yuksel, Texas A&M, 2010: Real-Time Water Waves with Wave Particles

Posted in Uncategorized | Comments Off on PhD Thesis

Hybrid Multiresolution Wire

We describe a method for the visual interactive simulation of wires contacting with rigid multibodies. The physical model used is a hybrid combining lumped elements and massless quasistatic representations. The latter is based on a kinematic constraint preserving the total length of the wire along a segmented path which can involve multiple bodies simultaneously and dry frictional contact nodes used for roping, lassoing and fastening. These nodes provide stick and slide friction along edges of the contacting geometries. The lumped element resolution is adapted dynamically based on local stability criteria, becoming coarser as the tension increases, and up to the purely kinematic representation. Kinematic segments and contact nodes are added and deleted and propagated based on contact geometries and dry friction configurations. The method gives dramatic increase on both performance and robustness because it quickly decimates superfluous nodes without loosing stability, yet adapts to complex configurations with many contacts and high curvature, keeping a fixed, large integration time step. Numerical results demonstrating the performance and stability of the adaptive multiresolution scheme are presented along with an array of representative simulation examples illustrating the versatility of the frictional contact model.

Hybrid Multiresolution Wire

Posted in Uncategorized | Comments Off on Hybrid Multiresolution Wire

Constraint Fluids

We present a fluid simulation method based on Smoothed Particle Hydrodynamics (SPH) in which incompressibility and boundary conditions are enforced using holonomic kinematic constraints on the density. This formulation enables systematic multiphysics integration in which interactions are modeled via similar constraints between the fluid pseudo-particles and impenetrable surfaces of other bodies. These conditions embody Archimede’s principle for solids and thus buoyancy results as a direct consequence. We use a variational time stepping scheme suitable for general constrained multibody systems we call SPOOK. Each step requires the solution of only one Mixed Linear Complementarity Problem (MLCP) with very few inequalities, corresponding to solid boundary conditions. We solve this MLCP with a fast iterative method. Overall stability is vastly improved in comparison to the unconstrained version of SPH, and this allows much larger time steps, and an increase in overall performance by two orders of magnitude. Proof of concept is given for computer graphics applications and interactive simulations.

Constraint Fluids

Posted in Uncategorized | Comments Off on Constraint Fluids

A Hexahedral Multigrid Approach for Simulating Cuts in Deformable Objects

We present a hexahedral finite element method for simulating cuts in deformable bodies using the corotational formulation of strain at high computational efficiency. Key to our approach is a novel embedding of adaptive element refinements and topological changes of the simulation grid into a geometric multigrid solver. Starting with a coarse hexahedral simulation grid, this grid is adaptively refined at the surface of a cutting tool until a finest resolution level, and the cut is modeled by separating elements along the cell faces at this level. To represent the induced discontinuities on successive multigrid levels, the affected coarse grid cells are duplicated and the resulting connectivity components are distributed to either side of the cut. Drawing upon recent work on octree and multigrid schemes for the numerical solution of partial differential equations, we develop efficient algorithms for updating the systems of equations of the adaptive finite element discretization and the multigrid hierarchy. To construct a surface that accurately aligns with the cuts, we adapt the splitting cubes algorithm to the specific linked voxel representation of the simulation domain we use. The paper is completed by a convergence analysis of the finite element solver and a performance comparison to alternative numerical solution methods. These investigations show that our approach offers high computational efficiency and physical accuracy, and that it enables cutting of deformable bodies at very high resolutions.

A Hexahedral Multigrid Approach for Simulating Cuts in Deformable Objects

Posted in Uncategorized | Comments Off on A Hexahedral Multigrid Approach for Simulating Cuts in Deformable Objects

PhD Thesis

A thesis from 2009 that I missed:

Barbara Solenthaler, University of Zurich: Incompressible Fluid Simulation and Advanced Surface Handling with SPH

Posted in Uncategorized | 1 Comment

Physics-Animation Forum?

I was recently asked whether there exists any internet forum dedicated to physics-based animation. The two that came to mind are Gamedev’s Math & Physics forum and Bullet’s physics simulation forums (both linked in the right column of the page), though my impression is that these naturally tend to be focused towards games, and with somewhat greater emphasis on rigid bodies and collision detection. So I wanted to pose some questions to you…

  1. Are there other forum sites that focus on physical simulation for computer graphics? If so please post a link in the comments.
  2. Would a dedicated physics-based animation forum be of genuine use to you? (Alternatively, do existing forums already serve this purpose, or do you not see a need for such a forum?)

Feel free to share any other thoughts you have on the subject.

Posted in Uncategorized | 3 Comments

Iso-geometric Analysis Based on Catmull-Clark Subdivision Solids

We present a volumetric iso-geometric finite element analysis based on Catmull-Clark solids. This concept allows one to use the same representation for the modeling, the physical simulation, and the visualization, which optimizes the design process and narrows the gap between CAD and CAE. In our method the boundary of the solid model is a Catmull-Clark surface with optional corners and creases to support the modeling phase. The crucial point in the simulation phase is the need to perform efficient  integration for the elements. We propose a method similar to the standard subdivision surface evaluation technique, such that numerical quadrature can be used.
Experiments show that our approach converges faster than methods based on tri-linear and tri-quadratic elements. However, the topological structure of Catmull-Clark elements is as simple as the structure of linear elements. Furthermore, the Catmull-Clark elements we use are C2-continuous on the boundary and in the interior except for irregular vertices and edges.

Iso-geometric Analysis Based on Catmull-Clark Subdivision Solids

Posted in Uncategorized | Comments Off on Iso-geometric Analysis Based on Catmull-Clark Subdivision Solids

Fast and Scalable CPU/GPU Collision Detection for Rigid and Deformable Surfaces

We present a new hybrid CPU/GPU collision detection technique for rigid and deformable objects based on spatial subdivision. Our approach efficiently exploits the massive computational capabilities of modern CPUs and GPUs commonly found in off-the-shelf computer systems. The algorithm is specifically tailored to be highly scalable on both the CPU and the GPU sides. We can compute discrete and continuous external and self-collisions of non-penetrating rigid and deformable objects consisting of many tens of thousands of triangles in few milliseconds on a modern PC. Our approach is orders of magnitude faster than earlier CPU-based approaches and up to twice as fast as the most recent GPU-based techniques.

Fast and Scalable CPU/GPU Collision Detection for Rigid and Deformable Surfaces

Posted in Uncategorized | Comments Off on Fast and Scalable CPU/GPU Collision Detection for Rigid and Deformable Surfaces

PhD Theses

Chris Wojtan,  Georgia Tech: Animating Physical Phenomena with Embedded Surface Meshes

Andreas Söderström, Linköping University: Memory Efficient Methods for Eulerian Free Surface Fluid Animation

Posted in Uncategorized | Comments Off on PhD Theses