Saved in:
| Main Authors: | , , , |
|---|---|
| Format: | Preprint |
| Published: |
2025
|
| Subjects: | |
| Online Access: | https://arxiv.org/abs/2503.22695 |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| _version_ | 1866911151961407488 |
|---|---|
| author | Ma, Pengfei Cai, Li Wang, Xuan Gao, Hao |
| author_facet | Ma, Pengfei Cai, Li Wang, Xuan Gao, Hao |
| contents | The immersed boundary (IB) method has become a leading approach in cardiac fluid-structure interaction (FSI) modeling due to its ability to handle large deformations and complex geometries without requiring mesh regeneration. However, the use of nonlinear, fiber-reinforced hyperelastic materials for modeling soft cardiac tissues introduces challenges in computational efficiency, particularly due to the additional projection steps required for stability in the IB framework. These steps often involve sparse matrix storage and computation, which can degrade GPU performance. In this work, we present a novel, fully GPU-accelerated, matrix-free IB method for FSI in anatomically realistic cardiac models. By employing nodal coupling, our method eliminates the need for projection operations in the finite element space. Additionally, we solve the Navier-Stokes equations using Chorin's projection method combined with a matrix-free geometric multigrid solver, ensuring the entire FSI algorithm remains matrix-free and highly compatible with GPU acceleration. Our implementation features several GPU-specific optimizations, including the use of constant memory to store values of nodal basis functions and their derivatives at quadrature points, and texture memory to efficiently implement the semi-Lagrangian discretization of convection terms. These innovations maximize GPU utilization while preserving the complex mechanical behavior of soft cardiac tissue. Benchmark tests demonstrate that our GPU-accelerated solver achieves a $50\times$-$100\times$ speedup compared to a 20-core CPU implementation, with comparable accuracy. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2503_22695 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Fully GPU-Accelerated, Matrix-Free Immersed Boundary Method for Complex Fiber-reinforced Hyperelastic Cardiac Models Ma, Pengfei Cai, Li Wang, Xuan Gao, Hao Computational Physics The immersed boundary (IB) method has become a leading approach in cardiac fluid-structure interaction (FSI) modeling due to its ability to handle large deformations and complex geometries without requiring mesh regeneration. However, the use of nonlinear, fiber-reinforced hyperelastic materials for modeling soft cardiac tissues introduces challenges in computational efficiency, particularly due to the additional projection steps required for stability in the IB framework. These steps often involve sparse matrix storage and computation, which can degrade GPU performance. In this work, we present a novel, fully GPU-accelerated, matrix-free IB method for FSI in anatomically realistic cardiac models. By employing nodal coupling, our method eliminates the need for projection operations in the finite element space. Additionally, we solve the Navier-Stokes equations using Chorin's projection method combined with a matrix-free geometric multigrid solver, ensuring the entire FSI algorithm remains matrix-free and highly compatible with GPU acceleration. Our implementation features several GPU-specific optimizations, including the use of constant memory to store values of nodal basis functions and their derivatives at quadrature points, and texture memory to efficiently implement the semi-Lagrangian discretization of convection terms. These innovations maximize GPU utilization while preserving the complex mechanical behavior of soft cardiac tissue. Benchmark tests demonstrate that our GPU-accelerated solver achieves a $50\times$-$100\times$ speedup compared to a 20-core CPU implementation, with comparable accuracy. |
| title | Fully GPU-Accelerated, Matrix-Free Immersed Boundary Method for Complex Fiber-reinforced Hyperelastic Cardiac Models |
| topic | Computational Physics |
| url | https://arxiv.org/abs/2503.22695 |