Hani Cheik Sleiman,Kevin Moerman,Diana de Oliveira,Joseph Jacob,Nesrin Mogulkoc,Brian Davidson,Simon Walker-Samuel,Rebecca Shipley
Abstract:This paper presents a comprehensive and highly-automated open-source pipeline for simulating flow and flow-related processes in (embedded) tubular structures. Addressing a critical gap in computational fluid dynamics (CFD) and simulation sciences, it facilitates the transition from raw three-dimensional imaging, graph networks, or CAD models of tubular objects to refined, simulation-ready meshes. This transition, traditionally labor-intensive and challenging, is streamlined and highly-automated through a series of innovative steps that include surface mesh processing, centre-line construction, anisotropic mesh generation, and volumetric meshing, leading to Finite Element Method (FEM) simulations. The pipeline leverages a range of open-source software and libraries, notably GIBBON, FEniCS, and Paraview, to provide flexibility and broad applicability across different simulation scenarios, ranging from biomedical to industrial applications. We demonstrate the versatility of our approach through five distinct applications, including the mesh generation for soil-root systems, lung airways, microcirculation networks, and portal vein networks, each originating from a different data source. Moreover, for several of these cases, we incorporate Computational Fluid Dynamics (CFD) simulations and strategies for 3D-1D coupling between the embedding domain and the embedded structures. Finally, we outline some future perspectives aimed at enhancing accuracy, reducing computational time, and incorporating advanced modeling and boundary condition strategies to further refine the framework's capabilities.
What problem does this paper attempt to address?
### Problems Addressed by the Paper
The paper aims to address a critical issue in computational fluid dynamics (CFD) and simulation science: the process of converting raw 3D imaging, graph networks, or computer-aided design (CAD) models into fine meshes suitable for simulation. This process has traditionally been very time-consuming and challenging, especially when it comes to converting raw data into a refined, simulation-ready format. The paper proposes a comprehensive and highly automated open-source pipeline (Tube2FEM) that simplifies this process through a series of innovative steps, including surface mesh processing, centerline construction, anisotropic mesh generation, and volumetric mesh partitioning, ultimately enabling finite element method (FEM) simulation. The pipeline leverages a range of open-source software and libraries, such as GIBBON, FEniCS, and Paraview, providing flexibility and broad applicability across different simulation scenarios, covering multiple fields from biomedical to industrial applications.
### Key Points of the Solution
1. **High Automation**: Reduces manual intervention and increases efficiency through a series of automated steps.
2. **Flexible Input**: Supports various input forms, including 3D imaging data, synthetic/skeletonized networks, and CAD models.
3. **Fine Mesh Generation**: Ensures precise mesh partitioning in tubular structures of varying diameters through anisotropic mesh generation techniques.
4. **Multidimensional Coupling**: Supports 3D-1D coupling, enhancing simulation accuracy.
5. **Wide Applicability**: Suitable for various application scenarios, including soil-root systems, lung airways, microcirculation networks, and portal vein networks.
### Application Examples
The paper demonstrates the versatility and effectiveness of the framework through five different application cases, including:
- Mesh generation for soil-root systems
- Mesh generation for lung airways
- Mesh generation for microcirculation networks
- Mesh generation for portal vein networks
These cases not only validate the feasibility of the framework but also showcase its broad application potential across different data sources.
### Future Prospects
The paper also outlines future development directions, including improving simulation accuracy, reducing computation time, and introducing advanced modeling and boundary condition strategies to further enhance the framework's capabilities.