AFSI: Automated Fluid-Structure Interaction Solver Development for Nonlinear Solid Mechanics

Abstract: AFSI is a novel, open-source fluid-structure interaction (FSI) solver that extends the capabilities of the FEniCS finite element library through an immersed boundary (IB) framework. Designed to simulate large deformations in hyperelastic materials (such as cardiac tissue), AFSI avoids the need for expensive remeshing by coupling a Lagrangian representation of the solid with an Eulerian description of the surrounding fluid. This approach retains the full expressiveness of FEniCS's variational formulations, function spaces, and time integration schemes. Implemented in a hybrid Python/C++ architecture, AFSI allows users to define geometries, constitutive models (e.g., the Holzapfel-Ogden law for myocardium), and strain energy functions directly in Python, while delegating performance-critical tasks such as assembly and linear solvers to optimized C++ backends. Its concise and modular Python API facilitates the setup of FSI simulations, enabling users to easily modify discretization strategies or analyze results using standard FEniCS post-processing tools. By combining the flexibility of FEniCS with a robust immersed boundary formulation, AFSI empowers rapid prototyping of complex nonlinear solid-fluid interaction problems, making it a powerful tool for simulating biomechanical systems and other applications involving highly deformable structures in flow.