title: Constraint fluids.
abstract: 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 pseudoparticles 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.
authors.affiliation: HPC2N/UMIT, Umeå University, Umeå, Sweden. bodin@hpc2n.umu.se
authors.forename: Kenneth
authors.initials: K
authors.lastname: Bodin
authors.forename: Claude
authors.initials: C
authors.lastname: Lacoursière
authors.forename: Martin
authors.initials: M
authors.lastname: Servin
mesh_terms: D000465:Algorithms; D003196:Computer Graphics; D003198:Computer Simulation; D057446:Hydrodynamics; D008962:Models, Theoretical