Parallel Simulation of High Reynolds Number Vascular Flows
|Title||Parallel Simulation of High Reynolds Number Vascular Flows|
|Year of Publication||2006|
|Authors||Fischer, PF, Loth, F, Lee, S-W, Smith, D, Tufo, HM, Bassiouny, HS|
The simulation of turbulent vascular flows presents significant numerical challenges. Because such flows are only weakly turbulent (i.e., transitional), they lack an inertial subrange that is amenable to subgrid-scale (SGS) modeling required for large-eddy or Reynolds-averaged Navier-Stokes simulations. The only reliable approach at present is to directly resolve all scales of motion. While the Reynolds number is not high (Re=1000-2000, typ.), the physical dissipation is small, and high-order methods are essential for efficiency. Weakly turbulent blood flow, such as occurs in post-stenotic regions or subsequent to graft implantation, exhibits a much broader range of scales than does its laminar (healthy) counterpart and thus requires an order of magnitude increase in spatial and temporal resolution, making fast iterative solvers and parallel computing necessities. The paper is organized as follows. Section 2 provides a brief overview of the governing equations, time advancement scheme, and spectral element method. Section 3 describes boundary condition treatment for simulating transition in bifurcation geometries. Section 4 presents parallel considerations and performance results, and Section 5 gives results for transitional flow in an arteriovenous graft model.