Misztal, Marek Krzysztof2; Erleben, Kenny3; Bargteil, Adam4; Fursund, Jens5; Christensen, Brian Bunch6; Bærentzen, Jakob Andreas2; Bridson, Robert7
1 Department of Informatics and Mathematical Modeling, Technical University of Denmark2 Department of Applied Mathematics and Computer Science, Technical University of Denmark3 University of Copenhagen4 University of Utah5 Imagination Technologies6 Alexandra Institute7 University of British Columbia
In this paper, we present a method for animating multiphase flow of immiscible fluids using unstructured moving meshes. Our underlying discretization is an unstructured tetrahedral mesh, the deformable simplicial complex (DSC), that moves with the flow in a Lagrangian manner. Mesh optimization operations improve element quality and avoid element inversion. In the context of multiphase flow, we guarantee that every element is occupied by a single fluid and, consequently, the interface between fluids is represented by a set of faces in the simplicial complex. This approach ensures that the underlying discretization matches the physics and avoids the additional book-keeping required in grid-based methods where multiple fluids may occupy the same cell. Our Lagrangian approach naturally leads us to adopt a finite element approach to simulation, in contrast to the finite volume approaches adopted by a majority of fluid simulation techniques that use tetrahedral meshes. We characterize fluid simulation as an optimization problem allowing for full coupling of the pressure and velocity fields and the incorporation of a second-order surface energy. We introduce a preconditioner based on the diagonal Schur complement and solve our optimization on the GPU. We provide the results of parameter studies as well as a performance analysis of our method, together with suggestions for performance optimization.
I E E E Transactions on Visualization and Computer Graphics, 2013, Vol 20, Issue 1, p. 4-16