Fluid-Structure Interaction Based on HPC Multicode Coupling

Abstract

\bfA \bfb \bfs \bft \bfr \bfa \bfc \bft. The fluid-structure interaction (FSI) problem has received great attention in the last few years, mainly because it is present in many physical systems, industrial applications, and almost every biological system. In the parallel computational field, outstanding advances have been achieved for the individual components of the problem, allowing, for instance, simulations around complex geometries at very high Reynolds numbers or simulations of the contraction of a beating heart. However, it is not an easy task to combine the advances of both fields, given that they have followed development paths in a rather independent way, and also because physical and numerical instabilities arise when dealing with two highly nonlinear partial differential equations. Nonetheless, in the last few years great advances in the coupled FSI field have been achieved, recognizing the most challenging problems to tackle and enabling a new generation of numerical simulations in aerodynamics, biological systems, and complex industrial devices. Keeping in mind that efficient parallel codes for the individual components already exist, this paper presents a framework to build a massively parallel FSI solver in a multicode coupling partitioned approach, with strong focus in the parallel implementation aspects and the parallel performance of the resulting application. The problem is casted in an algebraic form, and the main points of interest are the parallel environment needed to be able to transfer data among the codes, the location of the exchange surface, and the exchange of information among the parallel applications. The proposed framework has been implemented in the HPC multiphysics code Alya, and the multicode coupling is carried out running separated instances of this code. Two coupling algorithms with different acceleration schemes are revised, and three representative cases of different areas of interest showing the reach of the proposed framework are solved. Good agreement with literature and experiments is obtained. In addition to the numerical validation of the FSI solver, an assessment of the parallel performance of the proposed multicode strategy is done. In particular, a special distribution of the fluid code and solid code MPI processes on the supercomputer nodes based on computing cores overloading is investigated. Finally, a strong scalability test, running up to a 30 million elements case using 1280 MPI processes, is done. \bfK \bfe \bfy \bfw \bfo \bfr \bfd \bfs. fluid-structure interaction, HPC, multicode coupling \bfA \bfM \bfS \bfs \bfu \bfb \bfj \bfe \bfc \bft \bfc \bfl \bfa \bfs \bfs \bfi fi\bfc \bfa \bft \bfi \bfo \bfn \bfs. 74F10, 68U20, 6804 \bfD \bfO \bfI. 10.1137/17M1138868 1. Introduction. Fluid-structure interaction (FSI) governs so many physical problems, industrial devices, and biological systems that its relevance is out of the question. FSI combines the challenges of computational fluid dynamics (CFD) with those of computational solid dynamics (CSD).

Fluid–structure interaction simulations outperform computational fluid dynamics in the description of thoracic aorta haemodynamics and in the differentiation of progressive dilation in Marfan syndrome patients

“…The efficacy of such an approach was remarked by the substantial error reduction in the estimation of aortic haemodynamics. The algorithm used for FSI and the solid and fluid mechanics parts individually, and its implementation in our code ALYA, have already been deeply analysed in previous works [40,48]. There, important matters such as validation, mesh dependency, accuracy, efficiency and implementation aspects, especially those related to the numeric and parallel issues, are widely described.…”

The Cooperative Parallel: A Discussion About Run-Time Schedulers for Nested Parallelism

“…Coupled runs, consisting in simulations that solve different physics for a single run, are very common in HPC environments [10]. They can be found in a variety of examples from earth science, where some processes simulate the earth while other the ocean, to biological ones.…”

A Parallel Partitioned Approach On Fluid-Structure Interaction Simulations Using the Multiscale Universal Interface Coupling Library

“…This phenomenon occurs in a wide range of scientific research and engineering applications at different spatial and temporal scales. Examples include the hydroelasticity of a hydrofoil, the aeroelasticity of turbine blades and the flow of blood in an artery [1,2]. To accurately and efficiently tackle large-scale FSI problems involving billions of degrees of freedom, such a tool needs to show good scalability.…”