Parallel adaptive finite element methods for simulation of turbulent flow and fluid-structure interaction

Dnr:

SNIC 2019/3-354

Type:

SNIC Medium Compute

Principal Investigator:

Johan Hoffman

Affiliation:

Kungliga Tekniska högskolan

Start Date:

2019-07-01

End Date:

2020-07-01

Primary Classification:

10105: Computational Mathematics

Secondary Classification:

20301: Applied Mechanics

Tertiary Classification:

20603: Medical Image Processing

Webpage:

- Centre Storage at NSC: 500 GiB
- Beskow at PDC: 200 x 1000 core-h/month
- Tetralith at NSC: 100 x 1000 core-h/month

The research objectives of this project is to develop parallel adaptive finite element methods for simulation of turbulent flow and fluid-structure interaction (FSI).
The underlying mathematical equations are the fundamental conservation laws of mass, momentum and energy. Simulation of turbulent flow is challenging since the range of turbulent scales requires a level of computational resolution that is beyond the capacity of the methods, software and hardware of today; simulation of FSI is challenging since it represents a multi-physics problem where two (or more) models are coupled.
We have pioneered (i) a method for turbulence simulation that avoids explicit parameterization of unresolved turbulent scales [1]; (ii) a monolithic approach to FSI that circumvents the coupling problem, by expressing the FSI problem as a unified continuum model [2]. Both (i) and (ii) are imbedded into the underlying methodology of adaptive finite element methods, implemented in the open source FEniCS-HPC framework [3] that we develop in our research group. In benchmark workshops organized by NASA, Boeing and others, we have proven our technology [4], and since 2008 SNIC has supported our research and the development of open source software on Sweden’s most powerful supercomputing resources. In an ongoing project funded by the Swedish Research Council we have developed a clinical pathway for patient-specific simulation of the blood flow in the left ventricle of the human heart [5,6].
The objectives for this project are: (i) to investigate fundamental mathematical questions regarding existence of solutions to the mathematical equations that describe turbulent flow (with University College Dublin and Imperial College London); (ii) patient-specific simulation of the blood flow in the left ventricle of the human heart (with Karolinska hospital, Karolinska Institute and CNRS Paris); (iii) fluid-structure interaction simulation of blood flow in the brain (with University of Illinois), and (iv) fluid-structure interaction of human cells (with Linköping University).
[1] J. Hoffman et al., Towards a parameter-free method for high Reynolds number turbulent flow simulation based on adaptive finite element approximation, Comput. Meth. Appl. Mech. Engrg., Vol.288, pp.60-74, 2015.
[2] J. Hoffman et al., Unified continuum modeling of fluid-structure interaction, Mathematical Models and Methods in Applied Sciences, Vol.21(3), pp.491-513, 2011.
[3] J.Hoffman et al., FEniCS-HPC: Coupled Multiphysics in Computational Fluid Dynamics, Jülich Aachen Research Alliance (JARA) High-Performance Computing Symposium, Springer, pp.58—69, 2016.
[4] Johan Jansson et al., Time-resolved Adaptive Direct FEM Simulation of High-lift Aircraft Configurations, in Numerical Simulation of the Aerodynamics of High-Lift Configurations, pp.67-92, Springer, 2018.
[5] D. Larsson et al., "Patient-Specific Left Ventricular Flow Simulations From Transthoracic Echocardiography : Robustness Evaluation and Validation Against Ultrasound Doppler and Magnetic Resonance Imaging," IEEE Transactions on Medical Imaging, vol. 36, no. 11, s. 2261-2275, 2017.
[6] J. H. Spühler et al., "3D Fluid-Structure Interaction Simulation of Aortic Valves Using a Unified Continuum ALE FEM Model," Frontiers in Physiology, vol. 9, 2018.