Hemodynamic simulation of aneurysmal aorta using a parallel algorithm and a flow-rate based resistance splitting method.

Hemodynamic analysis of the aortic aneurysm can provide clinically useful information for the treatment and management of the disease. However, it is still a big challenge to accurately resolve the flow field within a clinically acceptable computing time. In this work, we adopt a two-level domain decomposition-based preconditioned parallel nonlinear solver to solve the unsteady Navier-Stokes equations with three-element Windkessel outlet boundary conditions, where model parameters are estimated by a novel method based on the blood flow rate ratio obtained from the ultrasound. The simulated results are verified by using clinically measured data. Several numerical experiments are performed to show the robustness of the proposed algorithm against different parameters. The proposed algorithm scales up to 2400 processor cores with 50% parallel efficiency and greatly reduces the computing time of the full-size aortic blood flow simulation.