Description
- Abstract:
- In this Thesis we focus on simulations of blood flow in three-dimensional patient-specific arterial networks. We employ high-order spectral/hp-element spatial discretization and concentrate on computational efficiency in solving multi-million degrees of freedom (DOF) flow problems on petaflop computers. We develop new two-level domain decomposition method and multilevel communicating interface for ultra-parallel flow simulations. Specifically, at the coarse level the computational domain is subdivided into several big patches. Within each patch a spectral element discretization (fine level) is employed. New interface conditions for the Navier-Stokes equations are developed. The proposed numerical approach has been tested in arterial flow simulations with up to 147 arteries. Solution of 2.87B DOF problem was computed on 18576 processors in less than one second at each time step. A scalable and fast parallel low-energy bases preconditioner (LEBP) in conjunction with coarse-space linear vertex solver is developed. We provide details on optimization, parallel performance and implementation of the coarse-space solver and show scalability of LEBP on thousands processors of the IBM BlueGene/L and the Cray XT3. An embarrassingly parallel but extremely efficient accelerator for iterative solver has been proposed. The new approach reduces the number of conjugate gradient iterations and exhibits grid independent scaling, while the computational overhead is negligible. A new type of outflow boundary condition for networks with multiple outlets has been developed. The method is based on a time-dependent resistance-capacitance model, where the resistance values are related to the measured flowrates. The aforementioned methods have been employed to study unsteady flow in patient-specific intracranial arterial trees. Results of a comparative study of 3D rigid wall and 1D flexible wall modeling of flow in complex arterial networks are presented. Transient turbulent flow in a stenosed carotid artery has been studied in details. To analyze the intermittent in time and space laminar-turbulent flow a new methodology based on time-space-window Proper Orthogonal Decomposition (POD) is proposed. A simplified version of the POD analysis that utilizes 2D slices only - more appropriate in the clinical setting - is also investigated.
- Notes:
- Thesis (Ph.D.) -- Brown University (2009)
Access Conditions
- Rights
- In Copyright
- Restrictions on Use
- Collection is open for research.
Citation
Grinberg, Leopold,
"Topics in Ultrascale Scientific Computing with Application in Biomedical Modeling"
(2009).
Applied Mathematics Theses and Dissertations.
Brown Digital Repository. Brown University Library.
https://doi.org/10.7301/Z0X0659T
Relations
Collection:
-
Applied Mathematics Theses and Dissertations
Theses and Dissertations for the Applied Mathematics department....