Central moments multiple relaxation time LBM for hemodynamic simulations in intracranial aneurysms: An in-vitro validation study using PIV and PC-MRI. (April 2021)

• Record Type: Journal Article
• Authors: Hosseini, Seyed Ali; Berg, Philipp; Huang, Feng; Roloff, Christoph; Janiga, Gábor; Thévenin, Dominique
• Is Part Of: Computers in biology and medicine. Volume 131(2021)
• Publication Date: 2021-04
• Subjects: Lattice Boltzmann method -- Particle image velocimetry -- Intracranial aneurysm -- Magnetic resonance imaging -- Validation -- Single relaxation time -- Central hermite multiple relaxation time -- Computational fluid dynamics; 00–01 -- 99-00; Medicine -- Data processing -- Periodicals; Biology -- Data processing -- Periodicals; 610.285
• DOI: 10.1016/j.compbiomed.2021.104251
• Languages: English
• ISSNs: 0010-4825

Abstract

Abstract: The lattice Boltzmann method (LBM) has recently emerged as an efficient alternative to classical Navier-Stokes solvers. This is particularly true for hemodynamics in complex geometries. However, in its most basic formulation, i.e. with the so-called single relaxation time (SRT) collision operator, it has been observed to have a limited stability domain in the Courant/Fourier space, strongly constraining the minimum time-step and grid size. The development of improved collision models such as the multiple relaxation time (MRT) operator in central moments space has tremendously widened the stability domain, while allowing to overcome a number of other well-documented artifacts, therefore opening the door for simulations over a wider range of grid and time-step sizes. The present work focuses on implementing and validating a specific collision operator, the central Hermite moments multiple relaxation time model with the full expansion of the equilibrium distribution function, to simulate blood flows in intracranial aneurysms. The study further proceeds with a validation of the numerical model through different test-cases and against experimental measurements obtained via stereoscopic particle image velocimetry (PIV) and phase-contrast magnetic resonance imaging (PC-MRI). For a patient-specific aneurysm both PIV and PC-MRI agree fairly well with the simulation. Finally, low-resolution simulations were shown to be able to capture blood flow information with sufficient accuracy, as demonstrated through both qualitative and quantitative analysis of the flow field while leading to strongly reduced computation times. For instance in the case of the patient-specific configuration, increasing the grid-size by a factor of two led to a reduction of computation time by a factor of 14 with very good similarity indices still ranging from 0.83 to 0.88. Highlights: Lattice Boltzmann simulation of intracranial aneurysm flows. Use of Galilean invariant collision operator based on central Hermite moments. Cross-validation with stereoscopic particle image velocimetry and phase-contrast magnetic resonance imaging. Under-resolved simulations via extended numerical properties of collision operator.