Speaker
Description
Modern Fortran offers powerful features for building extensible and efficient solvers in computational science. We present a fully Fortran based framework for simulating coupled fluid–structure interactions involving arrays of flexible cilia and deformable particles immersed in flow. The solver employs the immersed boundary method to couple a Navier–Stokes fluid solver with Kirchhoff rod models for elastic cilia and finite element models for deformable particles.
Key design elements include the use of derived types, allocatable arrays, and modular programming to represent cilia and particles as objects with geometry, forces, and update procedures. This object-oriented approach enables handling large arrays of interacting structures with clarity and efficiency. The entire solver, including fluid, solid, and coupling modules, is implemented in Fortran, demonstrating its
continuing relevance for multiphysics simulations. GPU acceleration of the pressure Poisson equation is optionally achieved through Fortran–C interoperability with NVIDIA’s AmgX library.
This work illustrates how modern Fortran facilitates high-performance scientific computing while supporting modular, maintainable code design. The presented cilia–particle interaction study show-cases Fortran’s capability to address complex biological and physical problems in computational fluid dynamics.
