Intracranial hemodynamics simulations: An efficient and accurate immersed boundary scheme
Authors:
D. S. Lampropoulos,
G. C. Bourantas,
B. F. Zwick,
G. C. Kagadis,
A. Wittek,
K. Miller,
V. C. Loukopoulos
Abstract:
Computational fluid dynamics (CFD) studies have been increasingly used for blood flow simulations in intracranial aneurysms (ICAs). However, despite the continuous progress of body-fitted CFD solvers, generating a high quality mesh is still the bottleneck of the CFD simulation, and strongly affects the accuracy of the numerical solution. To overcome this challenge, which will allow preforming CFD…
▽ More
Computational fluid dynamics (CFD) studies have been increasingly used for blood flow simulations in intracranial aneurysms (ICAs). However, despite the continuous progress of body-fitted CFD solvers, generating a high quality mesh is still the bottleneck of the CFD simulation, and strongly affects the accuracy of the numerical solution. To overcome this challenge, which will allow preforming CFD simulations efficiently for a large number of aneurysm cases we use an Immersed Boundary (IB) method. The proposed scheme relies on Cartesian grids to solve the incompressible Navier-Stokes (N-S) equations, using a finite element solver, and Lagrangian points to discretize the immersed object. All grid generations are conducted through automated algorithms which require no user input. Consequently, we verify the proposed method by comparing our numerical findings (velocity values) with published experimental results. Finally, we test the ability of the scheme to efficiently handle hemodynamic simulations on complex geometries on a sample of four patient-specific intracranial aneurysms.
△ Less
Submitted 27 July, 2020;
originally announced July 2020.
Immersed boundary finite element method for blood flow simulation
Authors:
G. C. Bourantas,
D. L. Lampropoulos,
B. F. Zwick,
V. C. Loukopoulos,
A. Wittek,
K. Miller
Abstract:
We present an efficient and accurate immersed boundary (IB) finite element (FE) solver for numerically solving incompressible Navier--Stokes equations. Particular emphasis is given to internal flows with complex geometries (blood flow in the vasculature system). IB methods are computationally costly for internal flows, mainly due to the large percentage of grid points that lie outside the flow dom…
▽ More
We present an efficient and accurate immersed boundary (IB) finite element (FE) solver for numerically solving incompressible Navier--Stokes equations. Particular emphasis is given to internal flows with complex geometries (blood flow in the vasculature system). IB methods are computationally costly for internal flows, mainly due to the large percentage of grid points that lie outside the flow domain. In this study, we apply a local refinement strategy, along with a domain reduction approach in order to reduce the grid that covers the flow domain and increase the percentage of the grid nodes that fall inside the flow domain. The proposed method utilizes an efficient and accurate FE solver with the incremental pressure correction scheme (IPCS), along with the boundary condition enforced IB method to numerically solve the transient, incompressible Navier--Stokes flow equations. We verify the accuracy of the numerical method using the analytical solution for Poiseuille flow in a cylinder. We further examine the accuracy and applicability of the proposed method by considering flow within complex geometries, such as blood flow in aneurysmal vessels and the aorta, flow configurations which would otherwise be extremely difficult to solve by most IB methods. Our method offers high accuracy, as demonstrated by the verification examples, and high efficiency, as demonstrated through the solution of blood flow within complex geometry on an off-the-shelf laptop computer.
△ Less
Submitted 4 July, 2020;
originally announced July 2020.