Turbulent fluid-structure simulation of human vocal folds, heart valves, and marine renewable energy devices

SNIC 2017/1-110


SNAC Medium

Principal Investigator:

Johan Jansson


Kungliga Tekniska högskolan

Start Date:


End Date:


Primary Classification:

10105: Beräkningsmatematik

Secondary Classification:

20306: Strömningsmekanik och akustik



We have led the work package on "Simulation of Phonation" in the EUNISON ( EU FP7 project, with a 3 year duration and 3.9M EUR budget ending mid-2016. In this project we will carry out simulations to publish several papers on the methods and software developed in the project which represent a novel direction in human voice simulation. The goal of the project is to simulate human voice generation by solving the unified domain equations for the air flow and elastic deformation in the vocal folds, leading to a turbulent fluid-structure interaction (FSI) with vibrating and contacting vocal folds, generating the base sound of the human voice which is then modulated by the vocal tract and mouth. We will simulate the fluid-structure interaction using our adaptive FEM framework Unicorn/FEniCS-HPC [2] which has shown good results for turbulent flow and FSI problems [1, 2], and also good results for the vocal folds [1]. We have verified optimal strong scaling on both Lindren and Beskow up to ca. 5000 cores with a pure MPI backend using PETSc [2] and up to 12000 cores on Beskow using the experimental hybrid MPI+PGAS JANPACK backend [3]. We will specifically simulate the full realistic geometry of the human voice apparatus including aeroacoustics in close collaboration with the voice group at KTH headed by Sten Ternström and the acoustics group at La Salle in Barcelona headed by Oriol Guasch. We are now able to produce sound close to the spectrum of a human voice, a critical step forward in the voice community. We will also apply the developed FSI methodology [4] to heart simulation, to marine structures in recently started projects with Energimyndigheten, a collaboration with Tecnalia in the Basque Country funded by the ELKARTEK framework, and the H2020 MSO4SC project. The methodology will also be applied to the simulation of a Rolls Royce jet engine, funded by the GENTALVE Basque government ELKARTEK project. We will present our results at acoustics and computational mechanics conferences and in journals. [1] J. Jansson, A. Holmberg, R Vilela de Abreu, C. Degirmenci, J. Hoffman, M. K arlsson, M. Abom “Adaptive parallel finite element framework for simulation of vocal fold turbulent fluid-structure interaction”, Proceedings of Meetings on Acoustics, 2013 [2] J. Hoffman, J. Jansson, R. V. de Abreu, N. C. Degirmenci, N. Jansson, K. Muller, M. Nazarov, and J. H. Spuhler. Unicorn: Parallel adaptive finite element simulation of turbulent flow and fluid-structure interaction for deforming domains and complex geometry. Computers & Fluids, 2012 [3] N. Jansson. Optimizing Sparse Matrix Assembly in Finite Element Solvers with One-sided Communication. In High Performance Computing for Computational Science { VECPAR 2012, volume 7851 of Lecture Notes in Computer Science , pages 128{139. Springer Berlin Heidelberg, 2013. [4] Towards a parameter-free method for high Reynolds number turbulent flow simulation based on adaptive finite element approximation, J Hoffman, J Jansson, N Jansson, RV De Abreu, Computer Methods in Applied Mechanics and Engineering, 2015