Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces a heterogeneous CPU-GPU framework with a pipelined Parareal architecture to accelerate long-time simulations of filamentous microswimmers in viscous fluid. The framework leverages GPU kernels for spatial interactions using the Method of Regularized Stokeslets (MRS) and MPI-GPU pipelining for temporal concurrency, overlapping coarse and fine propagators. Numerical experiments demonstrate order-of-magnitude speedups compared to CPU-only methods, enabling scalable simulations of complex biological and physical systems.
Simulating microswimmers just got an order of magnitude faster thanks to a new heterogeneous pipelined parallel-in-time framework.
Simulating large-scale microswimmer dynamics in viscous fluid poses significant challenges due to the coupled high spatial and temporal complexity. Conventional high-performance computing (HPC) methods often address these two dimensions in isolation, leaving a critical gap for synergistic acceleration. This paper introduces a heterogeneous CPU--GPU computing framework specifically optimized for the long-time simulation of filamentous microswimmers in viscous fluid. We propose a two-level parallelization strategy: (1) high-intensity GPU kernels to resolve the quadratic spatial interactions given by the Method of Regularized Stokeslets (MRS), and (2) a distributed MPI-GPU pipelined Parareal architecture to exploit temporal concurrency. By mapping the asynchronous pipeline onto multiple GPU devices, our framework effectively overlaps coarse and fine propagators, overcoming the serial bottlenecks of traditional Parareal method. Furthermore, we employ a GPU-optimized numerical routine for computing the matrix square root arising in the numerical scheme of the filamentous microswimmer simulations. Theoretical analysis of the efficiency improvement of the pipelined Parareal is presented. Numerical experiments demonstrate that the proposed framework achieves order-of-magnitude speedups over CPU-only methods, providing a scalable pathway for simulating complex emergent behaviors in large-scale biology and physics systems.