A Multilevel SPH Solver with Unified Solid Boundary Handling

Takahashi, Tetsuya

doi:10.1111/cgf.13048

Cited by 9 publications

(5 citation statements)

References 57 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Two‐way coupled simulation methods that combine rigid body contact handling with fluid dynamics have been heavily investigated in many contexts including Lattice–Boltzmann fluids [TIR06], height‐field fluids [TMFSG07, CM10, SBC*11], Lagrangian vortex methods [VHLL14], MPM fluids [YLCH18, HFG*18], Eulerian fluids [CMT04, GSLF05, BBB07, EWC*10] and Lagrangian fluids [MST*04, BTT09, AIA*12, MMCK14, TL16, KB17, BKWK19, GPB*19]. For a detailed discussion on the collision handling used by the mentioned simulation methods, we refer to Gissler et al.…”

Section: Related Workmentioning

confidence: 99%

Monolithic Friction and Contact Handling for Rigid Bodies and Fluids Using SPH

Probst

Teschner

2023

Computer Graphics Forum

View full text Add to dashboard Cite

We propose a novel monolithic pure SPH formulation to simulate fluids strongly coupled with rigid bodies. This includes fluid incompressibility, fluid-rigid interface handling and rigid-rigid contact handling with a viable implicit particle-based dry friction formulation. The resulting global system is solved using a new accelerated solver implementation that outperforms existing fluid and coupled rigid-fluid simulation approaches. We compare results of our simulation method to analytical solutions, show performance evaluations of our solver and present a variety of new and challenging simulation scenarios.

show abstract

Section: Related Workmentioning

confidence: 99%

Monolithic Friction and Contact Handling for Rigid Bodies and Fluids Using SPH

Probst

Teschner

2023

Computer Graphics Forum

View full text Add to dashboard Cite

show abstract

“…Dick et al conducted a thorough convergence analysis of multigrid methods as a direct solver and a preconditioner and proposed a topology‐aware multigrid solver that yields high convergence rates for complex scenes . Takahashi and Lin solved linear systems arising from particle‐based methods using the multigrid CG method . Aanjaneya et al solved the pressure Poisson equation using the multigrid CG method adapted to the power diagram discretization …”

Section: Previous Workmentioning

confidence: 99%

Sawtooth cycle revisited

Tsuruga

Iwasaki

2018

Computer Animation & Virtual

View full text Add to dashboard Cite

Solving the pressure Poisson equation dominates a large portion of computational time for incompressible fluid flow simulations. To solve the pressure Poisson equation efficiently, geometric multigrid methods are used directly or used as the preconditioner for the conjugate gradient (CG) method. Conventionally, the V-cycle multigrid method is widely employed, and little attentionhas been paid to other cycles. In this paper, we introduce the sawtooth cycle multigrid method and its simple extension called N-cycle, which provides better convergence in equal time comparison with the V-cycle as a direct solver of the pressure Poisson equation. We also apply the N-cycle to the preconditioner of the CG method and show that the N-cycle multigrid CG method can provide better convergence than the V-cycle multigrid CG method.

show abstract

“…Incorporating these ideas into our solver ensures accurate pressure solutions and avoids voxelization artifacts. While we focus on adaptive approaches for Eulerian methods, Lagrangian SPH and vortex methods have also been applied in the context of adaptivity [Adams et al 2007;Golas et al 2012;Solenthaler and Gross 2011;Takahashi and Lin 2016].…”

Section: Related Workmentioning

confidence: 99%

Power diagrams and sparse paged grids for high resolution adaptive liquids

et al. 2017

View full text Add to dashboard Cite

Fig. 1. (Left)Water filling a river bed surrounded by a canyon, with effective resolution 512 2 × 1024. Three refinement levels are used, based on proximity to the terrain. (Right) Sources inject water into a container and collide to form a thin sheet, with effective resolution 512 3 . Adaptivity pattern shown on background.We present an efficient and scalable octree-inspired fluid simulation framework with the flexibility to leverage adaptivity in any part of the computational domain, even when resolution transitions reach the free surface. Our methodology ensures symmetry, definiteness and second order accuracy of the discrete Poisson operator, and eliminates numerical and visual artifacts of prior octree schemes. This is achieved by adapting the operators acting on the octree's simulation variables to reflect the structure and connectivity of a power diagram, which recovers primal-dual mesh orthogonality and eliminates problematic T-junction configurations. We show how such operators can be efficiently implemented using a pyramid of sparsely populated uniform grids, enhancing the regularity of operations and facilitating parallelization. A novel scheme is proposed for encoding the topology of the power diagram in the neighborhood of each octree cell, allowing us to locally reconstruct it on the fly via a lookup table, rather than resorting to costly explicit meshing. The pressure Poisson equation is solved via a highly efficient, matrix-free multigrid preconditioner for Conjugate Gradient, adapted to the power diagram discretization. We use another sparsely † M. Aanjaneya and M. Gao are joint first authors. * M. Aanjaneya was with the University of Wisconsin -Madison during this work. This work was supported in part by National Science Foundation grants IIS-1253598, CCF-1423064, CCF-1533885 and by the Natural Sciences and Engineering Research Council of Canada under grant RGPIN-04360-2014. The authors are grateful to Nathan Mitchell for his indispensable help with modeling and rendering of examples. C. Batty would like to thank Ted Ying for carrying out preliminary explorations on quadtrees. Permission to make digital or hard copies of all or part of this work for personal or classroom use is granted without fee provided that copies are not made or distributed for profit or commercial advantage and that copies bear this notice and the full citation on the first page. Copyrights for components of this work owned by others than ACM must be honored. Abstracting with credit is permitted. To copy otherwise, or republish, to post on servers or to redistribute to lists, requires prior specific permission and/or a fee. Request permissions from permissions@acm.org. populated uniform grid for high resolution interface tracking with a narrow band level set representation. Using the recently introduced SPGrid data structure, sparse uniform grids in both the power diagram discretization and our narrow band level set can be compactly stored and efficiently updated via streaming operations. Additionally, we present enhancem...

show abstract

A Multilevel SPH Solver with Unified Solid Boundary Handling

Cited by 9 publications

References 57 publications

Monolithic Friction and Contact Handling for Rigid Bodies and Fluids Using SPH

Monolithic Friction and Contact Handling for Rigid Bodies and Fluids Using SPH

Sawtooth cycle revisited

Power diagrams and sparse paged grids for high resolution adaptive liquids

Contact Info

Product

Resources

About