A Scalable Sharp-Interface Immersed Boundary Method for Moving Boundary Flow Simulation on Multi-GPU Systems
Abstract
Sharp-interface immersed boundary methods provide a versatile framework for simulating flows with complex moving boundaries, but their irregular data structures, dynamic sparsity, and discontinuous pressure operators present significant challenges for efficient multi-GPU implementation. This work presents a fully GPU-native, multi-GPU acceleration of the sharp-interface immersed boundary solver, ViCar3D. The implementation combines a hybrid OpenACC/CUDA programming model with a matrix decomposition that isolates dynamically changing immersed-boundary operators, GPU-optimized geometric kernels, asynchronous dual-stream sparse matrix--vector execution, and a symmetric-memory-based sparse communication strategy. A GPU-oriented iterative solver employing a preconditioned stabilized bi-conjugate gradient (BiCGStab) method is developed to efficiently treat the pressure Poisson equation on highly stretched Cartesian grids. The resulting framework preserves the accuracy and versatility of the original solver while delivering substantial acceleration and scalability for simulations involving complex moving and deforming immersed boundaries, thereby extending the applicability of sharp-interface immersed boundary methods to extreme-scale computational fluid dynamics.
Related
- Mechanism-Resolved Interface Momentum Transfer in Immersed-Boundary Lattice Boltzmann Simulations
- On the hydrodynamic behaviour of the immersed boundary-lattice Boltzmann method for wetting problems (a)
- High-order immersed boundary method for subsonic aeroacoustics and compressible Navier-Stokes equationsHigh-order immersed boundary method for subsonic aeroacoustics and compressible Navier-Stokes equations
- Handling geometric singularities efficiently with the immersed-boundary method: application to riblets
- A ghost-cell immersed boundary method for the shallow water equations on Cartesian grids: mass conservation and wetting and drying over variable topography
BibTeX
@misc{kumar2026scalable,
title = {{A Scalable Sharp-Interface Immersed Boundary Method for Moving Boundary Flow Simulation on Multi-GPU Systems}},
author = {Sushrut Kumar and Joshua Romero and Massimiliano Fatica and Jung-Hee Seo and Rajat Mittal},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.21032},
eprint = {2610.21032},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.21032},
note = {CC BY-NC-SA 4.0},
}
RIS
TY - UNPB
TI - A Scalable Sharp-Interface Immersed Boundary Method for Moving Boundary Flow Simulation on Multi-GPU Systems
AU - Sushrut Kumar
AU - Joshua Romero
AU - Massimiliano Fatica
AU - Jung-Hee Seo
AU - Rajat Mittal
PY - 2026
DA - 2026/10/15
PB - arXiv
DO - 10.48550/arxiv.2610.21032
UR - https://arxiv.org/abs/2610.21032
AB - Sharp-interface immersed boundary methods provide a versatile framework for simulating flows with complex moving boundaries, but their irregular data structures, dynamic sparsity, and discontinuous pressure operators present significant challenges for efficient multi-GPU implementation. This work presents a fully GPU-native, multi-GPU acceleration of the sharp-interface immersed boundary solver, ViCar3D. The implementation combines a hybrid OpenACC/CUDA programming model with a matrix decomposition that isolates dynamically changing immersed-boundary operators, GPU-optimized geometric kernels, asynchronous dual-stream sparse matrix--vector execution, and a symmetric-memory-based sparse communication strategy. A GPU-oriented iterative solver employing a preconditioned stabilized bi-conjugate gradient (BiCGStab) method is developed to efficiently treat the pressure Poisson equation on highly stretched Cartesian grids. The resulting framework preserves the accuracy and versatility of the original solver while delivering substantial acceleration and scalability for simulations involving complex moving and deforming immersed boundaries, thereby extending the applicability of sharp-interface immersed boundary methods to extreme-scale computational fluid dynamics.
ER -