1
|
Nahas H, Ishii T, Yiu BYS, Yu ACH. A GPU-Based, Real-Time Dealiasing Framework for High-Frame-Rate Vector Doppler Imaging. IEEE TRANSACTIONS ON ULTRASONICS, FERROELECTRICS, AND FREQUENCY CONTROL 2023; 70:1384-1400. [PMID: 37549086 DOI: 10.1109/tuffc.2023.3303349] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/09/2023]
Abstract
Vector Doppler is well regarded as a potential way of deriving flow vectors to intuitively visualize complex flow profiles, especially when it is implemented at high frame rates. However, this technique's performance is known to suffer from aliasing artifacts. There is a dire need to devise real-time dealiasing solutions for vector Doppler. In this article, we present a new methodological framework for achieving aliasing-resistant flow vector estimation at real-time throughput from precalculated Doppler frequencies. Our framework comprises a series of compute kernels that have synergized: 1) an extended least squares vector Doppler (ELS-VD) algorithm; 2) single-instruction, multiple-thread (SIMT) processing principles; and 3) implementation on a graphical processing unit (GPU). Results show that this new framework, when executed on an RTX-2080 GPU, can effectively generate aliasing-free flow vector maps using high-frame-rate imaging datasets acquired from multiple transmit-receive angle pairs in a carotid phantom imaging scenario. Over the entire cardiac cycle, the frame processing time for aliasing-resistant vector estimation was measured to be less than 16 ms, which corresponds to a minimum processing throughput of 62.5 frames/s. In a human femoral bifurcation imaging trial with fast flow (150 cm/s), our framework was found to be effective in resolving two-cycle aliasing artifacts at a minimum throughput of 53 frames/s. The framework's processing throughput was generally in the real-time range for practical combinations of ELS-VD algorithmic parameters. Overall, this work represents the first demonstration of real-time, GPU-based aliasing-resistant vector flow imaging using vector Doppler estimation principles.
Collapse
|
2
|
Waldron C, Zafar MA, Ziganshin BA, Weininger G, Grewal N, Elefteriades JA. Evidence Accumulates: Patients with Ascending Aneurysms Are Strongly Protected from Atherosclerotic Disease. Int J Mol Sci 2023; 24:15640. [PMID: 37958625 PMCID: PMC10650782 DOI: 10.3390/ijms242115640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Revised: 10/23/2023] [Accepted: 10/24/2023] [Indexed: 11/15/2023] Open
Abstract
Ascending thoracic aortic aneurysms may be fatal upon rupture or dissection and remain a leading cause of death in the developed world. Understanding the pathophysiology of the development of ascending thoracic aortic aneurysms may help reduce the morbidity and mortality of this disease. In this review, we will discuss our current understanding of the protective relationship between ascending thoracic aortic aneurysms and the development of atherosclerosis, including decreased carotid intima-media thickness, low-density lipoprotein levels, coronary and aortic calcification, and incidence of myocardial infarction. We also propose several possible mechanisms driving this relationship, including matrix metalloproteinase proteins and transforming growth factor-β.
Collapse
Affiliation(s)
- Christina Waldron
- Aortic Institute at Yale-New Haven Hospital, Yale University School of Medicine, New Haven, CT 06519, USA; (C.W.); (M.A.Z.); (B.A.Z.)
| | - Mohammad A. Zafar
- Aortic Institute at Yale-New Haven Hospital, Yale University School of Medicine, New Haven, CT 06519, USA; (C.W.); (M.A.Z.); (B.A.Z.)
| | - Bulat A. Ziganshin
- Aortic Institute at Yale-New Haven Hospital, Yale University School of Medicine, New Haven, CT 06519, USA; (C.W.); (M.A.Z.); (B.A.Z.)
- Department of Cardiovascular and Endovascular Surgery, Kazan State Medical University, 420012 Kazan, Russia
| | - Gabe Weininger
- Aortic Institute at Yale-New Haven Hospital, Yale University School of Medicine, New Haven, CT 06519, USA; (C.W.); (M.A.Z.); (B.A.Z.)
| | - Nimrat Grewal
- Department of Cardiothoracic Surgery, Amsterdam University Medical Center, 1105 AZ Amsterdam, The Netherlands;
| | - John A. Elefteriades
- Aortic Institute at Yale-New Haven Hospital, Yale University School of Medicine, New Haven, CT 06519, USA; (C.W.); (M.A.Z.); (B.A.Z.)
| |
Collapse
|
3
|
Wei H, Amlani F, Pahlevan NM. Direct 0D-3D coupling of a lattice Boltzmann methodology for fluid-structure aortic flow simulations. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2023; 39:e3683. [PMID: 36629353 DOI: 10.1002/cnm.3683] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/13/2022] [Revised: 11/29/2022] [Accepted: 01/06/2023] [Indexed: 05/05/2023]
Abstract
This work introduces a numerical approach and implementation for the direct coupling of arbitrary complex ordinary differential equation- (ODE-)governed zero-dimensional (0D) boundary conditions to three-dimensional (3D) lattice Boltzmann-based fluid-structure systems for hemodynamics studies. In particular, a most complex configuration is treated by considering a dynamic left ventricle- (LV-)elastance heart model which is governed by (and applied as) a nonlinear, non-stationary hybrid ODE-Dirichlet system. Other ODE-based boundary conditions, such as lumped parameter Windkessel models for truncated vasculature, are also considered. Performance studies of the complete 0D-3D solver, including its treatment of the lattice Boltzmann fluid equations and elastodynamics equations as well as their interactions, is conducted through a variety of benchmark and convergence studies that demonstrate the ability of the coupled 0D-3D methodology in generating physiological pressure and flow waveforms-ultimately enabling the exploration of various physical and physiological parameters for hemodynamics studies of the coupled LV-arterial system. The methods proposed in this paper can be easily applied to other ODE-based boundary conditions as well as to other fluid problems that are modeled by 3D lattice Boltzmann equations and that require direct coupling of dynamic 0D boundary conditions.
Collapse
Affiliation(s)
- Heng Wei
- Department of Aerospace and Mechanical Engineering, University of Southern California, Los Angeles, California, USA
| | - Faisal Amlani
- Department of Aerospace and Mechanical Engineering, University of Southern California, Los Angeles, California, USA
- Université Paris-Saclay, CentraleSupélec, ENS Paris-Saclay, CNRS, LMPS - Laboratoire de Mécanique Paris-Saclay, Gif-sur-Yvette, France
| | - Niema M Pahlevan
- Department of Aerospace and Mechanical Engineering, University of Southern California, Los Angeles, California, USA
- School of Medicine, University of Southern California, Los Angeles, California, USA
| |
Collapse
|
4
|
Du P, Wang JX. Reducing Geometric Uncertainty in Computational Hemodynamics by Deep Learning-Assisted Parallel-Chain MCMC. J Biomech Eng 2022; 144:121009. [PMID: 36166284 PMCID: PMC9632478 DOI: 10.1115/1.4055809] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2022] [Revised: 09/21/2022] [Indexed: 11/08/2022]
Abstract
Computational hemodynamic modeling has been widely used in cardiovascular research and healthcare. However, the reliability of model predictions is largely dependent on the uncertainties of modeling parameters and boundary conditions, which should be carefully quantified and further reduced with available measurements. In this work, we focus on propagating and reducing the uncertainty of vascular geometries within a Bayesian framework. A novel deep learning (DL)-assisted parallel Markov chain Monte Carlo (MCMC) method is presented to enable efficient Bayesian posterior sampling and geometric uncertainty reduction. A DL model is built to approximate the geometry-to-hemodynamic map, which is trained actively using online data collected from parallel MCMC chains and utilized for early rejection of unlikely proposals to facilitate convergence with less expensive full-order model evaluations. Numerical studies on two-dimensional aortic flows are conducted to demonstrate the effectiveness and merit of the proposed method.
Collapse
Affiliation(s)
- Pan Du
- Department of Aerospace and Mechanical Engineering, University of Notre Dame, Notre Dame, IN 46556
| | - Jian-Xun Wang
- Department of Aerospace and Mechanical Engineering, University of Notre Dame, Notre Dame, IN 46556
| |
Collapse
|
5
|
Non-invasive detection of coronary artery disease from photoplethysmograph using lumped parameter modelling. Biomed Signal Process Control 2022. [DOI: 10.1016/j.bspc.2022.103781] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022]
|
6
|
Boumpouli M, Sauvage EL, Capelli C, Schievano S, Kazakidi A. Characterization of Flow Dynamics in the Pulmonary Bifurcation of Patients With Repaired Tetralogy of Fallot: A Computational Approach. Front Cardiovasc Med 2021; 8:703717. [PMID: 34660711 PMCID: PMC8514754 DOI: 10.3389/fcvm.2021.703717] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2021] [Accepted: 09/07/2021] [Indexed: 11/18/2022] Open
Abstract
The hemodynamic environment of the pulmonary bifurcation is of great importance for adult patients with repaired tetralogy of Fallot (rTOF) due to possible complications in the pulmonary valve and narrowing of the left pulmonary artery (LPA). The aim of this study was to computationally investigate the effect of geometrical variability and flow split on blood flow characteristics in the pulmonary trunk of patient-specific models. Data from a cohort of seven patients was used retrospectively and the pulmonary hemodynamics was investigated using averaged and MRI-derived patient-specific boundary conditions on the individualized models, as well as a statistical mean geometry. Geometrical analysis showed that curvature and tortuosity are higher in the LPA branch, compared to the right pulmonary artery (RPA), resulting in complex flow patterns in the LPA. The computational analysis also demonstrated high time-averaged wall shear stress (TAWSS) at the outer wall of the LPA and the wall of the RPA proximal to the junction. Similar TAWSS patterns were observed for averaged boundary conditions, except for a significantly modified flow split assigned at the outlets. Overall, this study enhances our understanding about the flow development in the pulmonary bifurcation of rTOF patients and associates some morphological characteristics with hemodynamic parameters, highlighting the importance of patient-specificity in the models. To confirm these findings, further studies are required with a bigger cohort of patients.
Collapse
Affiliation(s)
- Maria Boumpouli
- Department of Biomedical Engineering, University of Strathclyde, Glasgow, United Kingdom
| | - Emilie L. Sauvage
- Institute of Cardiovascular Science and Great Ormond Street Hospital for Children, NHS Foundation Trust, University College London, London, United Kingdom
| | - Claudio Capelli
- Institute of Cardiovascular Science and Great Ormond Street Hospital for Children, NHS Foundation Trust, University College London, London, United Kingdom
| | - Silvia Schievano
- Institute of Cardiovascular Science and Great Ormond Street Hospital for Children, NHS Foundation Trust, University College London, London, United Kingdom
| | - Asimina Kazakidi
- Department of Biomedical Engineering, University of Strathclyde, Glasgow, United Kingdom
| |
Collapse
|
7
|
Seven Mathematical Models of Hemorrhagic Shock. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2021; 2021:6640638. [PMID: 34188690 PMCID: PMC8195646 DOI: 10.1155/2021/6640638] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Accepted: 04/02/2021] [Indexed: 11/17/2022]
Abstract
Although mathematical modelling of pressure-flow dynamics in the cardiocirculatory system has a lengthy history, readily finding the appropriate model for the experimental situation at hand is often a challenge in and of itself. An ideal model would be relatively easy to use and reliable, besides being ethically acceptable. Furthermore, it would address the pathogenic features of the cardiovascular disease that one seeks to investigate. No universally valid model has been identified, even though a host of models have been developed. The object of this review is to describe several of the most relevant mathematical models of the cardiovascular system: the physiological features of circulatory dynamics are explained, and their mathematical formulations are compared. The focus is on the whole-body scale mathematical models that portray the subject's responses to hypovolemic shock. The models contained in this review differ from one another, both in the mathematical methodology adopted and in the physiological or pathological aspects described. Each model, in fact, mimics different aspects of cardiocirculatory physiology and pathophysiology to varying degrees: some of these models are geared to better understand the mechanisms of vascular hemodynamics, whereas others focus more on disease states so as to develop therapeutic standards of care or to test novel approaches. We will elucidate key issues involved in the modeling of cardiovascular system and its control by reviewing seven of these models developed to address these specific purposes.
Collapse
|
8
|
Fathi MF, Perez-Raya I, Baghaie A, Berg P, Janiga G, Arzani A, D'Souza RM. Super-resolution and denoising of 4D-Flow MRI using physics-Informed deep neural nets. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2020; 197:105729. [PMID: 33007592 DOI: 10.1016/j.cmpb.2020.105729] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/01/2020] [Accepted: 08/23/2020] [Indexed: 06/11/2023]
Abstract
BACKGROUND AND OBJECTIVE Time resolved three-dimensional phase contrast magnetic resonance imaging (4D-Flow MRI) has been used to non-invasively measure blood velocities in the human vascular system. However, issues such as low spatio-temporal resolution, acquisition noise, velocity aliasing, and phase-offset artifacts have hampered its clinical application. In this research, we developed a purely data-driven method for super-resolution and denoising of 4D-Flow MRI. METHODS The flow velocities, pressure, and the MRI image magnitude are modeled as a patient-specific deep neural net (DNN). For training, 4D-Flow MRI images in the complex Cartesian space are used to impose data-fidelity. Physics of fluid flow is imposed through regularization. Creative loss function terms have been introduced to handle noise and super-resolution. The trained patient-specific DNN can be sampled to generate noise-free high-resolution flow images. The proposed method has been implemented using the TensorFlow DNN library and tested on numerical phantoms and validated in-vitro using high-resolution particle image velocitmetry (PIV) and 4D-Flow MRI experiments on transparent models subjected to pulsatile flow conditions. RESULTS In case of numerical phantoms, we were able to increase spatial resolution by a factor of 100 and temporal resolution by a factor of 5 compared to the simulated 4D-Flow MRI. There is an order of magnitude reduction of velocity normalized root mean square error (vNRMSE). In case of the in-vitro validation tests with PIV as reference, there is similar improvement in spatio-temporal resolution. Although the vNRMSE is reduced by 50%, the method is unable to negate a systematic bias with respect to the reference PIV that is introduced by the 4D-Flow MRI measurement. CONCLUSIONS This work has demonstrated the feasibility of using the readily available machinery of deep learning to enhance 4D-Flow MRI using a purely data-driven method. Unlike current state-of-the-art methods, the proposed method is agnostic to geometry and boundary conditions and therefore eliminates the need for tedious tasks such as accurate image segmentation for geometry, image registration, and estimation of boundary flow conditions. Arbitrary regions of interest can be selected for processing. This work will lead to user-friendly analysis tools that will enable quantitative hemodynamic analysis of vascular diseases in a clinical setting.
Collapse
Affiliation(s)
- Mojtaba F Fathi
- Dept. of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, WI, USA
| | - Isaac Perez-Raya
- Dept. of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, WI, USA
| | - Ahmadreza Baghaie
- Dept. of Electrical and Computer Engineering, New York Institute of Technology, Long Island, NY, USA
| | - Philipp Berg
- Lab. of Fluid Dynamics and Technical Flows, University of Magdeburg, Germany; Research Campus STIMULATE, University of Magdeburg, Magdeburg, Germany
| | - Gabor Janiga
- Lab. of Fluid Dynamics and Technical Flows, University of Magdeburg, Germany; Research Campus STIMULATE, University of Magdeburg, Magdeburg, Germany
| | - Amirhossein Arzani
- Dept. of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, USA
| | - Roshan M D'Souza
- Dept. of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, WI, USA.
| |
Collapse
|
9
|
Perez-Raya I, Fathi MF, Baghaie A, Sacho RH, Koch KM, D'Souza RM. Towards multi-modal data fusion for super-resolution and denoising of 4D-Flow MRI. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2020; 36:e3381. [PMID: 32627366 DOI: 10.1002/cnm.3381] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/06/2020] [Revised: 05/26/2020] [Accepted: 05/27/2020] [Indexed: 06/11/2023]
Abstract
4D-Flow magnetic resonance imaging (MRI) has enabled in vivo time-resolved measurement of three-dimensional blood flow velocities in the human vascular system. However, its clinical use has been hampered by two main issues, namely, low spatio-temporal resolution and acquisition noise. While patient-specific computational fluid dynamics (CFD) simulations can address the resolution and noise issues, its fidelity is impacted by accuracy of estimation of boundary conditions, model parameters, vascular geometry, and flow model assumptions. In this paper a scheme to address limitations of both modalities through data-fusion is presented. The solutions of the patient-specific CFD simulation are characterized using proper orthogonal decomposition (POD). Next, a process of projecting the 4D-Flow MRI data onto the POD basis and projection coefficient mapping using generalized dynamic mode decomposition (DMD) enables simultaneous super-resolution and denoising of 4D-Flow MRI. The method has been tested using numerical phantoms derived from patient-specific aneurysmal geometries and applied to in vivo 4D-Flow MRI data.
Collapse
Affiliation(s)
- Isaac Perez-Raya
- Department of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, USA
| | - Mojtaba F Fathi
- Department of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, USA
| | - Ahmadreza Baghaie
- Department of Electrical and Computer Engineering, New York Institute of Technology, Long Island, New York, USA
| | - Raphael H Sacho
- Department of Neurosurgery, Medical College of Wisconsin, Wauwatosa, Wisconsin, USA
| | - Kevin M Koch
- Department of Radiology, Medical College of Wisconsin, Wauwatosa, Wisconsin, USA
| | - Roshan M D'Souza
- Department of Mechanical Engineering, University of Wisconsin-Milwaukee, Milwaukee, Wisconsin, USA
| |
Collapse
|
10
|
Wang J, Zhang S. Fluid shear stress modulates endothelial inflammation by targeting LIMS2. Exp Biol Med (Maywood) 2020; 245:1656-1663. [PMID: 32752897 DOI: 10.1177/1535370220943837] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Abstract
Mechanosensitive genes regulate multiple cardiovascular pathophysiological processes and disorders; however, the role of flow-sensitive genes in atherosclerosis is still unknown. In this study, we identify LIM Zinc Finger Domain Containing 2 (LIMS2) that acts as a mechanosensitive gene downregulated by disturbed flow (d-flow) both in human endothelial cells (ECs) in vitro and in mice in vivo. Mechanistically, d-flow suppresses LIMS2 expression, which leads to endothelial inflammation by upregulating typical inflammatory factors, VCAM-1, and ICAM-1 in human ECs. The findings indicate that LIMS2, the new flow-sensitive gene, may help us to find a new insight to explain how d-flow caused endothelial inflammation and provide a new therapeutic approach for atherosclerosis in the future.
Collapse
Affiliation(s)
- Junyao Wang
- Department of Cardiovascular Medicine, Xiangya Hospital, Central South University, Changsha, Hunan 410008, China
| | - Shiyanjin Zhang
- Xiangya School of Medicine, Central South University, Changsha, Hunan 410013, China
| |
Collapse
|
11
|
Correlation-based Doppler-angle estimation with plane-wave excitation. INFORMATICS IN MEDICINE UNLOCKED 2020. [DOI: 10.1016/j.imu.2020.100315] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022] Open
|
12
|
Kumar S, Williams D, Sur S, Wang JY, Jo H. Role of flow-sensitive microRNAs and long noncoding RNAs in vascular dysfunction and atherosclerosis. Vascul Pharmacol 2019; 114:76-92. [PMID: 30300747 PMCID: PMC6905428 DOI: 10.1016/j.vph.2018.10.001] [Citation(s) in RCA: 62] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2018] [Revised: 09/19/2018] [Accepted: 10/05/2018] [Indexed: 02/07/2023]
Abstract
Atherosclerosis is the primary underlying cause of myocardial infarction, ischemic stroke, and peripheral artery disease. The disease preferentially occurs in arterial regions exposed to disturbed blood flow, in part, by altering expression of flow-sensitive coding- and non-coding genes. In this review, we summarize the role of noncoding RNAs, [microRNAs (miRNAs) and long noncoding RNAs(lncRNAs)], as regulators of gene expression and outline their relationship to the pathogenesis of atherosclerosis. While miRNAs are small noncoding genes that post-transcriptionally regulate gene expression by targeting mRNA transcripts, the lncRNAs regulate gene expression by diverse mechanisms, which are still emerging and incompletely understood. We focused on multiple flow-sensitive miRNAs such as, miR-10a, -19a, -23b, -17~92, -21, -663, -92a, -143/145, -101, -126, -712, -205, and -155 that play a critical role in endothelial function and atherosclerosis by targeting inflammation, cell cycle, proliferation, migration, apoptosis, and nitric oxide signaling. Flow-dependent regulation of lncRNAs is just emerging, and their role in vascular dysfunction and atherosclerosis is unknown. Here, we discuss the flow-sensitive lncRNA STEEL along with other lncRNAs studied in the context of vascular pathophysiology and atherosclerosis such as MALAT1, MIAT1, ANRIL, MYOSLID, MEG3, SENCR, SMILR, LISPR1, and H19. Also discussed is the use of these noncoding RNAs as potential biomarkers and therapeutics to reduce and regress atherosclerosis.
Collapse
Affiliation(s)
- Sandeep Kumar
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, USA
| | - Darian Williams
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, USA
| | - Sanjoli Sur
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, USA
| | - Jun-Yao Wang
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, USA
| | - Hanjoong Jo
- Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology and Emory University, Atlanta, USA; Division of Cardiology, Emory University, Atlanta, USA.
| |
Collapse
|
13
|
Wain RAJ, Gaskell NJ, Fsadni AM, Francis J, Whitty JPM. Finite Element Predictions of Sutured and Coupled Microarterial Anastomoses. ADVANCED BIOMEDICAL ENGINEERING 2019. [DOI: 10.14326/abe.8.63] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022] Open
Affiliation(s)
- Richard AJ Wain
- John Tyndall Institute, School of Engineering, University of Central Lancashire
- School of Medicine and Dentistry, University of Central Lancashire
- Department of Plastic & Reconstructive Surgery, Royal Preston Hospital, Lancashire Teaching Hospitals NHS Foundation Trust
| | - Nicolas J Gaskell
- John Tyndall Institute, School of Engineering, University of Central Lancashire
| | - Andrew M Fsadni
- John Tyndall Institute, School of Engineering, University of Central Lancashire
| | - Jonathan Francis
- John Tyndall Institute, School of Engineering, University of Central Lancashire
| | - Justin PM Whitty
- John Tyndall Institute, School of Engineering, University of Central Lancashire
| |
Collapse
|
14
|
Melis A, Clayton RH, Marzo A. Bayesian sensitivity analysis of a 1D vascular model with Gaussian process emulators. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2017; 33. [PMID: 28337862 DOI: 10.1002/cnm.2882] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/15/2016] [Revised: 01/18/2017] [Accepted: 03/10/2017] [Indexed: 06/06/2023]
Abstract
One-dimensional models of the cardiovascular system can capture the physics of pulse waves but involve many parameters. Since these may vary among individuals, patient-specific models are difficult to construct. Sensitivity analysis can be used to rank model parameters by their effect on outputs and to quantify how uncertainty in parameters influences output uncertainty. This type of analysis is often conducted with a Monte Carlo method, where large numbers of model runs are used to assess input-output relations. The aim of this study was to demonstrate the computational efficiency of variance-based sensitivity analysis of 1D vascular models using Gaussian process emulators, compared to a standard Monte Carlo approach. The methodology was tested on four vascular networks of increasing complexity to analyse its scalability. The computational time needed to perform the sensitivity analysis with an emulator was reduced by the 99.96% compared to a Monte Carlo approach. Despite the reduced computational time, sensitivity indices obtained using the two approaches were comparable. The scalability study showed that the number of mechanistic simulations needed to train a Gaussian process for sensitivity analysis was of the order O(d), rather than O(d×103) needed for Monte Carlo analysis (where d is the number of parameters in the model). The efficiency of this approach, combined with capacity to estimate the impact of uncertain parameters on model outputs, will enable development of patient-specific models of the vascular system, and has the potential to produce results with clinical relevance.
Collapse
Affiliation(s)
- Alessandro Melis
- INSIGNEO Institute for in Silico Medicine, The University of Sheffield, Sheffield, U.K
- Department of Mechanical Engineering, The University of Sheffield, Sheffield, U.K
| | - Richard H Clayton
- INSIGNEO Institute for in Silico Medicine, The University of Sheffield, Sheffield, U.K
- Department of Computer Science, The University of Sheffield, Sheffield, U.K
| | - Alberto Marzo
- INSIGNEO Institute for in Silico Medicine, The University of Sheffield, Sheffield, U.K
- Department of Mechanical Engineering, The University of Sheffield, Sheffield, U.K
| |
Collapse
|
15
|
Chen Z, Zhou Y, Wang J, Liu X, Ge S, He Y. Modeling of coarctation of aorta in human fetuses using 3D/4D fetal echocardiography and computational fluid dynamics. Echocardiography 2017; 34:1858-1866. [PMID: 28833523 DOI: 10.1111/echo.13644] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
OBJECTIVES We sought to develop a hemodynamic model of aortic and ductal arches using computational fluid dynamics (CFD) and 3D/4D spatio-temporal image correlation (STIC) fetal echocardiography and to investigate the hemodynamics of coarctation of aorta (CoA) in human fetuses using this approach. METHODS We obtained 3D/4D STIC fetal echocardiographic images of the aortic and ductal arches (DA) in five normal fetuses. Based on these images, we simulated the hemodynamics in the two arches using CFD. Subsequently, we reduced the dimensions of aortic isthmus from 100% to 85%, 70%, 55%, 40%, and 25% of the original dimension digitally. Numerical simulation was repeated in each condition, and flow profile, velocity, pressure, and wall shear stress (WSS) were compared with those of the baseline normal aortic and ductal arches. RESULTS With the progressive narrowing in the aortic isthmus, there were alterations in the flow profile, velocity, pressure, and WSS. The downstream vortexes disappeared, and the double helix profile became single helix. When the aortic isthmus reduced by 55% in dimension, there was an exponential increase in velocity and WSS and decrease in pressure. CONCLUSIONS The aortic and ductal arch geometry and flow lead to the alterations in flow profile, velocity, pressure, and WSS in the aortic isthmus in normal and CoA models, which are conductive of ductal issue migration into these areas. A 55% reduction in the dimension of aortic isthmus is associated with exponential change in velocity, pressure, and WSS, a probable threshold for hemodynamically significant CoA.
Collapse
Affiliation(s)
- Zhuo Chen
- Department of Ultrasound, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
| | - Yue Zhou
- School of Aeronautic Sciences and Engineering, BEIHANG Unversity, Beijing, China
| | - Jingying Wang
- School of Aeronautic Sciences and Engineering, BEIHANG Unversity, Beijing, China
| | - Xiaowei Liu
- Department of Ultrasound, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
| | - Shuping Ge
- Section of Cardiology, St. Christopher's Hospital for Children and Drexel University College of Medicine, Philadelphia, PA, USA
| | - Yihua He
- Department of Ultrasound, Beijing Anzhen Hospital, Capital Medical University, Beijing, China
| |
Collapse
|
16
|
Abstract
This review examines four imaging modalities; ultrasound (US), digital subtraction angiography (DSA), magnetic resonance imaging (MRI) and computed tomography (CT), that have common or potential applications in vascular access (VA). The four modalities are reviewed under their primary uses, techniques, advantages and disadvantages, and future directions that are specific to VA. Currently, US is the most commonly used modality in VA because it is cheaper (relative to other modalities), accessible, non-ionising, and does not require the use of contrast agents. DSA is predominantly only performed when an intervention is indicated. MRI is limited by its cost and the time required for image acquisition that mainly confines it to the realm of research where high resolution is required. CT’s short acquisition times and high resolution make it useful as a problem-solving tool in complex cases, although accessibility can be an issue. All four imaging modalities have advantages and disadvantages that limit their use in this particular patient cohort. Current imaging in VA comprises an integrated approach with each modality providing particular uses dependent on their capabilities. MRI and CT, which currently have limited use, may have increasingly important future roles in complex cases where detailed analysis is required.
Collapse
|
17
|
Bakhshinejad A, Baghaie A, Vali A, Saloner D, Rayz VL, D'Souza RM. Merging computational fluid dynamics and 4D Flow MRI using proper orthogonal decomposition and ridge regression. J Biomech 2017; 58:162-173. [PMID: 28577904 PMCID: PMC5527690 DOI: 10.1016/j.jbiomech.2017.05.004] [Citation(s) in RCA: 36] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2016] [Revised: 04/24/2017] [Accepted: 05/05/2017] [Indexed: 10/19/2022]
Abstract
Time resolved phase-contrast magnetic resonance imaging 4D-PCMR (also called 4D Flow MRI) data while capable of non-invasively measuring blood velocities, can be affected by acquisition noise, flow artifacts, and resolution limits. In this paper, we present a novel method for merging 4D Flow MRI with computational fluid dynamics (CFD) to address these limitations and to reconstruct de-noised, divergence-free high-resolution flow-fields. Proper orthogonal decomposition (POD) is used to construct the orthonormal basis of the local sampling of the space of all possible solutions to the flow equations both at the low-resolution level of the 4D Flow MRI grid and the high-level resolution of the CFD mesh. Low-resolution, de-noised flow is obtained by projecting in vivo 4D Flow MRI data onto the low-resolution basis vectors. Ridge regression is then used to reconstruct high-resolution de-noised divergence-free solution. The effects of 4D Flow MRI grid resolution, and noise levels on the resulting velocity fields are further investigated. A numerical phantom of the flow through a cerebral aneurysm was used to compare the results obtained using the POD method with those obtained with the state-of-the-art de-noising methods. At the 4D Flow MRI grid resolution, the POD method was shown to preserve the small flow structures better than the other methods, while eliminating noise. Furthermore, the method was shown to successfully reconstruct details at the CFD mesh resolution not discernible at the 4D Flow MRI grid resolution. This method will improve the accuracy of the clinically relevant flow-derived parameters, such as pressure gradients and wall shear stresses, computed from in vivo 4D Flow MRI data.
Collapse
Affiliation(s)
- Ali Bakhshinejad
- Department of Mechanical Engineering, University of Wisconsin-Milwaukee, United States.
| | - Ahmadreza Baghaie
- Department of Biomedical Engineering, Purdue University, United States
| | - Alireza Vali
- Department of Radiology, Northwestern University, United States
| | - David Saloner
- Department of Radiology, College of Medicine, University of California, San Francisco, United States
| | - Vitaliy L Rayz
- Department of Biomedical Engineering, Purdue University, United States
| | - Roshan M D'Souza
- Department of Mechanical Engineering, University of Wisconsin-Milwaukee, United States
| |
Collapse
|
18
|
Zhou X, Xia C, Stephen G, Khan F, Corner GA, Hoskins PR, Huang Z. Investigation of Ultrasound-Measured Flow Velocity, Flow Rate and Wall Shear Rate in Radial and Ulnar Arteries Using Simulation. ULTRASOUND IN MEDICINE & BIOLOGY 2017; 43:981-992. [PMID: 28236532 DOI: 10.1016/j.ultrasmedbio.2016.12.015] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/28/2016] [Revised: 12/02/2016] [Accepted: 12/31/2016] [Indexed: 06/06/2023]
Abstract
Parameters of blood flow measured by ultrasound in radial and ulnar arteries, such as flow velocity, flow rate and wall shear rate, are widely used in clinical practice and clinical research. Investigation of these measurements is useful for evaluating accuracy and providing knowledge of error sources. A method for simulating the spectral Doppler ultrasound measurement process was developed with computational fluid dynamics providing flow-field data. Specific scanning factors were adjusted to investigate their influence on estimation of the maximum velocity waveform, and flow rate and wall shear rate were derived using the Womersley equation. The overestimation in maximum velocity increases greatly (peak systolic from about 10% to 30%, time-averaged from about 30% to 50%) when the beam-vessel angle is changed from 30° to 70°. The Womersley equation was able to estimate flow rate in both arteries with less than 3% error, but performed better in the radial artery (2.3% overestimation) than the ulnar artery (15.4% underestimation) in estimating wall shear rate. It is concluded that measurements of flow parameters in the radial and ulnar arteries with clinical ultrasound scanners are prone to clinically significant errors.
Collapse
Affiliation(s)
- Xiaowei Zhou
- School of Science and Engineering, University of Dundee, Dundee, United Kingdom
| | - Chunming Xia
- School of Mechanical and Power Engineering, East China University of Science and Technology, Shanghai, China; NHS Tayside Medical Physics, Ninewells Hospital, Dundee, United Kingdom
| | - Gandy Stephen
- School of Mechanical and Power Engineering, East China University of Science and Technology, Shanghai, China; NHS Tayside Medical Physics, Ninewells Hospital, Dundee, United Kingdom
| | - Faisel Khan
- Cardiovascular and Diabetes Medicine, Ninewells Hospital and Medical School, University of Dundee, Dundee, United Kingdom
| | - George A Corner
- School of Science and Engineering, University of Dundee, Dundee, United Kingdom
| | - Peter R Hoskins
- Centre for Cardiovascular Science, University of Edinburgh, Edinburgh, United Kingdom
| | - Zhihong Huang
- School of Science and Engineering, University of Dundee, Dundee, United Kingdom.
| |
Collapse
|
19
|
Botar CC, Tóth ÁÁ, Klisurić OR, Nićiforović DD, Vučaj Ćirilović VA, Till VE. Dynamic simulation and Doppler Ultrasonography validation of blood flow behavior in Abdominal Aortic Aneurysm. Phys Med 2017; 37:1-8. [DOI: 10.1016/j.ejmp.2017.03.021] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/13/2016] [Revised: 03/20/2017] [Accepted: 03/27/2017] [Indexed: 10/19/2022] Open
|
20
|
Dumas L, El Bouti T, Lucor D. A Robust and Subject-Specific Hemodynamic Model of the Lower Limb Based on Noninvasive Arterial Measurements. J Biomech Eng 2017; 139:2565258. [PMID: 27685182 DOI: 10.1115/1.4034833] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2016] [Indexed: 11/08/2022]
Abstract
Cardiovascular diseases are currently the leading cause of mortality in the population of developed countries, due to the constant increase in cardiovascular risk factors, such as high blood pressure, cholesterol, overweight, tobacco use, lack of physical activity, etc. Numerous prospective and retrospective studies have shown that arterial stiffening is a relevant predictor of these diseases. Unfortunately, the arterial stiffness distribution across the human body is difficult to measure experimentally. We propose a numerical approach to determine the arterial stiffness distribution of an arterial network using a subject-specific one-dimensional model. The proposed approach calibrates the optimal parameters of the reduced-order model, including the arterial stiffness, by solving an inverse problem associated with the noninvasive in vivo measurements. An uncertainty quantification analysis has also been carried out to measure the contribution of the model input parameters variability, alone or by interaction with other inputs, to the variation of clinically relevant hemodynamic indices, here the arterial pulse pressure. The results obtained for a lower limb model, demonstrate that the numerical approach presented here can provide a robust and subject-specific tool to the practitioner, allowing an early and reliable diagnosis of cardiovascular diseases based on a noninvasive clinical examination.
Collapse
Affiliation(s)
- Laurent Dumas
- Professor Lab. de Mathématiques de Versailles, CNRS, UVSQ, Université Paris-Saclay, Versailles 78035, France e-mail:
| | - Tamara El Bouti
- MINES ParisTech PSL, Centre D'efficacité Énergétique des Systèmes, 5 rue Léon Blum, Palaiseau 91120, France e-mail:
| | - Didier Lucor
- LIMSI, CNRS, Université Paris-Saclay, Campus Universitaire, bât 508, Rue John von Neumann, Orsay Cedex F-91405, France e-mail:
| |
Collapse
|
21
|
Daughter Sac Formation Related to Blood Inflow Jet in an Intracranial Aneurysm. World Neurosurg 2016; 96:396-402. [DOI: 10.1016/j.wneu.2016.09.040] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2016] [Revised: 09/08/2016] [Accepted: 09/08/2016] [Indexed: 10/21/2022]
|
22
|
Computational Hemodynamic Analysis for the Diagnosis of Atherosclerotic Changes in Intracranial Aneurysms: A Proof-of-Concept Study Using 3 Cases Harboring Atherosclerotic and Nonatherosclerotic Aneurysms Simultaneously. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2016; 2016:2386031. [PMID: 27703491 PMCID: PMC5040129 DOI: 10.1155/2016/2386031] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/06/2016] [Accepted: 08/04/2016] [Indexed: 11/17/2022]
Abstract
This was a proof-of-concept computational fluid dynamics (CFD) study designed to identify atherosclerotic changes in intracranial aneurysms. We selected 3 patients with multiple unruptured aneurysms including at least one with atherosclerotic changes and investigated whether an image-based CFD study could provide useful information for discriminating the atherosclerotic aneurysms. Patient-specific geometries were constructed from three-dimensional data obtained using rotational angiography. Transient simulations were conducted under patient-specific inlet flow rates measured by phase-contrast magnetic resonance velocimetry. In the postanalyses, we calculated time-averaged wall shear stress (WSS), oscillatory shear index, and relative residence time (RRT). The volume of blood flow entering aneurysms through the neck and the mean velocity of blood flow inside aneurysms were examined. We applied the age-of-fluid method to quantitatively assess the residence of blood inside aneurysms. Atherosclerotic changes coincided with regions exposed to disturbed blood flow, as indicated by low WSS and long RRT. Blood entered aneurysms in phase with inlet flow rates. The mean velocities of blood inside atherosclerotic aneurysms were lower than those inside nonatherosclerotic aneurysms. Blood in atherosclerotic aneurysms was older than that in nonatherosclerotic aneurysms, especially near the wall. This proof-of-concept study demonstrated that CFD analysis provided detailed information on the exchange and residence of blood that is useful for the diagnosis of atherosclerotic changes in intracranial aneurysms.
Collapse
|
23
|
Simmons RD, Kumar S, Thabet SR, Sur S, Jo H. Omics-based approaches to understand mechanosensitive endothelial biology and atherosclerosis. WILEY INTERDISCIPLINARY REVIEWS-SYSTEMS BIOLOGY AND MEDICINE 2016; 8:378-401. [PMID: 27341633 DOI: 10.1002/wsbm.1344] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/10/2016] [Revised: 04/11/2016] [Accepted: 04/12/2016] [Indexed: 12/16/2022]
Abstract
Atherosclerosis is a multifactorial disease that preferentially occurs in arterial regions exposed to d-flow can be used to indicate disturbed flow or disturbed blood flow. The mechanisms by which d-flow induces atherosclerosis involve changes in the transcriptome, methylome, proteome, and metabolome of multiple vascular cells, especially endothelial cells. Initially, we begin with the pathogenesis of atherosclerosis and the changes that occur at multiple levels owing to d-flow, especially in the endothelium. Also, there are a variety of strategies used for the global profiling of the genome, transcriptome, miRNA-ome, DNA methylome, and metabolome that are important to define the biological and pathophysiological mechanisms of endothelial dysfunction and atherosclerosis. Finally, systems biology can be used to integrate these 'omics' datasets, especially those that derive data based on a single animal model, in order to better understand the pathophysiology of atherosclerosis development in a holistic manner and how this integrative approach could be used to identify novel molecular diagnostics and therapeutic targets to prevent or treat atherosclerosis. WIREs Syst Biol Med 2016, 8:378-401. doi: 10.1002/wsbm.1344 For further resources related to this article, please visit the WIREs website.
Collapse
Affiliation(s)
- Rachel D Simmons
- The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Sandeep Kumar
- The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Salim Raid Thabet
- The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Sanjoli Sur
- The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| | - Hanjoong Jo
- The Wallace H. Coulter Department of Biomedical Engineering, Georgia Institute of Technology, Atlanta, GA, USA
| |
Collapse
|
24
|
The role of endothelial mechanosensitive genes in atherosclerosis and omics approaches. Arch Biochem Biophys 2015; 591:111-31. [PMID: 26686737 DOI: 10.1016/j.abb.2015.11.005] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2015] [Revised: 10/29/2015] [Accepted: 11/04/2015] [Indexed: 12/24/2022]
Abstract
Atherosclerosis is the leading cause of morbidity and mortality in the U.S., and is a multifactorial disease that preferentially occurs in regions of the arterial tree exposed to disturbed blood flow. The detailed mechanisms by which d-flow induces atherosclerosis involve changes in the expression of genes, epigenetic patterns, and metabolites of multiple vascular cells, especially endothelial cells. This review presents an overview of endothelial mechanobiology and its relation to the pathogenesis of atherosclerosis with special reference to the anatomy of the artery and the underlying fluid mechanics, followed by a discussion of a variety of experimental models to study the role of fluid mechanics and atherosclerosis. Various in vitro and in vivo models to study the role of flow in endothelial biology and pathobiology are discussed in this review. Furthermore, strategies used for the global profiling of the genome, transcriptome, miR-nome, DNA methylome, and metabolome, as they are important to define the biological and pathophysiological mechanisms of atherosclerosis. These "omics" approaches, especially those which derive data based on a single animal model, provide unprecedented opportunities to not only better understand the pathophysiology of atherosclerosis development in a holistic and integrative manner, but also to identify novel molecular and diagnostic targets.
Collapse
|
25
|
Rispoli VC, Nielsen JF, Nayak KS, Carvalho JLA. Computational fluid dynamics simulations of blood flow regularized by 3D phase contrast MRI. Biomed Eng Online 2015; 14:110. [PMID: 26611470 PMCID: PMC4661988 DOI: 10.1186/s12938-015-0104-7] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2015] [Accepted: 11/16/2015] [Indexed: 11/23/2022] Open
Abstract
Background Phase contrast magnetic resonance imaging (PC-MRI) is used clinically for quantitative assessment of cardiovascular flow and function, as it is capable of providing directly-measured 3D velocity maps. Alternatively, vascular flow can be estimated from model-based computation fluid dynamics (CFD) calculations. CFD provides arbitrarily high resolution, but its accuracy hinges on model assumptions, while velocity fields measured with PC-MRI generally do not satisfy the equations of fluid dynamics, provide limited resolution, and suffer from partial volume effects. The purpose of this study is to develop a proof-of-concept numerical procedure for constructing a simulated flow field that is influenced by both direct PC-MRI measurements and a fluid physics model, thereby taking advantage of both the accuracy of PC-MRI and the high spatial resolution of CFD. The use of the proposed approach in regularizing 3D flow fields is evaluated. Methods The proposed algorithm incorporates both a Newtonian fluid physics model and a linear PC-MRI signal model. The model equations are solved numerically using a modified CFD algorithm. The numerical solution corresponds to the optimal solution of a generalized Tikhonov regularization, which provides a flow field that satisfies the flow physics equations, while being close enough to the measured PC-MRI velocity profile. The feasibility of the proposed approach is demonstrated on data from the carotid bifurcation of one healthy volunteer, and also from a pulsatile carotid flow phantom. Results The proposed solver produces flow fields that are in better agreement with direct PC-MRI measurements than CFD alone, and converges faster, while closely satisfying the fluid dynamics equations. For the implementation that provided the best results, the signal-to-error ratio (with respect to the PC-MRI measurements) in the phantom experiment was 6.56 dB higher than that of conventional CFD; in the in vivo experiment, it was 2.15 dB higher. Conclusions The proposed approach allows partial or complete measurements to be incorporated into a modified CFD solver, for improving the accuracy of the resulting flow fields estimates. This can be used for reducing scan time, increasing the spatial resolution, and/or denoising the PC-MRI measurements.
Collapse
Affiliation(s)
- Vinicius C Rispoli
- Department of Electrical Engineering, University of Brasilia, Brasília, Brazil. .,UnB Gama College, University of Brasilia, Brasília, Brazil.
| | - Jon F Nielsen
- fMRI Laboratory, Biomedical Engineering Department, University of Michigan, Ann Arbor, USA.
| | - Krishna S Nayak
- Magnetic Resonance Engineering Laboratory, Ming Hsieh Department of Electrical Engineering, University of Southern California, Los Angeles, USA.
| | - Joao L A Carvalho
- Department of Electrical Engineering, University of Brasilia, Brasília, Brazil.
| |
Collapse
|
26
|
Aristokleous N, Seimenis I, Georgiou GC, Nicolaides A, Anayiotos AS. The Effect of Head Rotation on the Geometry and Hemodynamics of Healthy Vertebral Arteries. Ann Biomed Eng 2015; 43:1287-97. [DOI: 10.1007/s10439-015-1340-5] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2015] [Accepted: 05/13/2015] [Indexed: 12/30/2022]
|
27
|
Du T, Hu D, Cai D. Outflow boundary conditions for blood flow in arterial trees. PLoS One 2015; 10:e0128597. [PMID: 26000782 PMCID: PMC4441455 DOI: 10.1371/journal.pone.0128597] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2014] [Accepted: 04/28/2015] [Indexed: 11/18/2022] Open
Abstract
In the modeling of the pulse wave in the systemic arterial tree, it is necessary to truncate small arterial crowns representing the networks of small arteries and arterioles. Appropriate boundary conditions at the truncation points are required to represent wave reflection effects of the truncated arterial crowns. In this work, we provide a systematic method to extract parameters of the three-element Windkessel model from the impedance of a truncated arterial tree or from experimental measurements of the blood pressure and flow rate at the inlet of the truncated arterial crown. In addition, we propose an improved three-element Windkessel model with a complex capacitance to accurately capture the fundamental-frequency time lag of the reflection wave with respect to the incident wave. Through our numerical simulations of blood flow in a single artery and in a large arterial tree, together with the analysis of the modeling error of the pulse wave in large arteries, we show that both a small truncation radius and the complex capacitance in the improved Windkessel model play an important role in reducing the modeling error, defined as the difference in dynamics induced by the structured tree model and the Windkessel models.
Collapse
Affiliation(s)
- Tao Du
- Department of Mathematics, Institute of Natural Sciences, and MOE-LSC, Shanghai Jiao Tong University, Shanghai, China
| | - Dan Hu
- Department of Mathematics, Institute of Natural Sciences, and MOE-LSC, Shanghai Jiao Tong University, Shanghai, China
- * E-mail:
| | - David Cai
- Department of Mathematics, Institute of Natural Sciences, and MOE-LSC, Shanghai Jiao Tong University, Shanghai, China
- Courant Institute of Mathematical Sciences and Center for Neural Science, New York University, New York, U.S.A.
- NYUAD Institute, New York University Abu Dhabi, Abu Dhabi, UAE
| |
Collapse
|
28
|
Yiu BYS, Lai SSM, Yu ACH. Vector projectile imaging: time-resolved dynamic visualization of complex flow patterns. ULTRASOUND IN MEDICINE & BIOLOGY 2014; 40:2295-309. [PMID: 24972498 DOI: 10.1016/j.ultrasmedbio.2014.03.014] [Citation(s) in RCA: 112] [Impact Index Per Article: 11.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/17/2013] [Revised: 03/02/2014] [Accepted: 03/10/2014] [Indexed: 05/22/2023]
Abstract
Achieving non-invasive, accurate and time-resolved imaging of vascular flow with spatiotemporal fluctuations is well acknowledged to be an ongoing challenge. In this article, we present a new ultrasound-based framework called vector projectile imaging (VPI) that can dynamically render complex flow patterns over an imaging view at millisecond time resolution. VPI is founded on three principles: (i) high-frame-rate broad-view data acquisition (based on steered plane wave firings); (ii) flow vector estimation derived from multi-angle Doppler analysis (coupled with data regularization and least-squares fitting); (iii) dynamic visualization of color-encoded vector projectiles (with flow speckles displayed as adjunct). Calibration results indicated that by using three transmit angles and three receive angles (-10°, 0°, +10° for both), VPI can consistently compute flow vectors in a multi-vessel phantom with three tubes positioned at different depths (1.5, 4, 6 cm), oriented at different angles (-10°, 0°, +10°) and of different sizes (dilated diameter: 2.2, 4.4 and 6.3 mm; steady flow rate: 2.5 mL/s). The practical merit of VPI was further illustrated through an anthropomorphic flow phantom investigation that considered both healthy and stenosed carotid bifurcation geometries. For the healthy bifurcation with 1.2-Hz carotid flow pulses, VPI was able to render multi-directional and spatiotemporally varying flow patterns (using a nominal frame rate of 416 fps or 2.4-ms time resolution). In the case of stenosed bifurcations (50% eccentric narrowing), VPI enabled dynamic visualization of flow jet and recirculation zones. These findings suggest that VPI holds promise as a new tool for complex flow analysis.
Collapse
Affiliation(s)
- Billy Y S Yiu
- Medical Engineering Program, University of Hong Kong, Pokfulam, Hong Kong
| | - Simon S M Lai
- Medical Engineering Program, University of Hong Kong, Pokfulam, Hong Kong
| | - Alfred C H Yu
- Medical Engineering Program, University of Hong Kong, Pokfulam, Hong Kong.
| |
Collapse
|
29
|
Accuracy and Reproducibility of Patient-Specific Hemodynamic Models of Stented Intracranial Aneurysms: Report on the Virtual Intracranial Stenting Challenge 2011. Ann Biomed Eng 2014; 43:154-67. [DOI: 10.1007/s10439-014-1082-9] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2014] [Accepted: 07/31/2014] [Indexed: 10/24/2022]
|
30
|
Aristokleous N, Seimenis I, Georgiou GC, Papaharilaou Y, Brott BC, Nicolaides A, Anayiotos AS. Impact of Head Rotation on the Individualized Common Carotid Flow and Carotid Bifurcation Hemodynamics. IEEE J Biomed Health Inform 2014; 18:783-9. [DOI: 10.1109/jbhi.2014.2305575] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
31
|
Sugiyama SI, Niizuma K, Nakayama T, Shimizu H, Endo H, Inoue T, Fujimura M, Ohta M, Takahashi A, Tominaga T. Relative Residence Time Prolongation in Intracranial Aneurysms. Neurosurgery 2013; 73:767-76. [DOI: 10.1227/neu.0000000000000096] [Citation(s) in RCA: 55] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
|
32
|
Marsden AL. Simulation based planning of surgical interventions in pediatric cardiology. PHYSICS OF FLUIDS (WOODBURY, N.Y. : 1994) 2013; 25:101303. [PMID: 24255590 PMCID: PMC3820639 DOI: 10.1063/1.4825031] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2013] [Accepted: 09/22/2013] [Indexed: 05/17/2023]
Abstract
Hemodynamics plays an essential role in the progression and treatment of cardiovascular disease. However, while medical imaging provides increasingly detailed anatomical information, clinicians often have limited access to hemodynamic data that may be crucial to patient risk assessment and treatment planning. Computational simulations can now provide detailed hemodynamic data to augment clinical knowledge in both adult and pediatric applications. There is a particular need for simulation tools in pediatric cardiology, due to the wide variation in anatomy and physiology in congenital heart disease patients, necessitating individualized treatment plans. Despite great strides in medical imaging, enabling extraction of flow information from magnetic resonance and ultrasound imaging, simulations offer predictive capabilities that imaging alone cannot provide. Patient specific simulations can be used for in silico testing of new surgical designs, treatment planning, device testing, and patient risk stratification. Furthermore, simulations can be performed at no direct risk to the patient. In this paper, we outline the current state of the art in methods for cardiovascular blood flow simulation and virtual surgery. We then step through pressing challenges in the field, including multiscale modeling, boundary condition selection, optimization, and uncertainty quantification. Finally, we summarize simulation results of two representative examples from pediatric cardiology: single ventricle physiology, and coronary aneurysms caused by Kawasaki disease. These examples illustrate the potential impact of computational modeling tools in the clinical setting.
Collapse
Affiliation(s)
- Alison L Marsden
- Mechanical and Aerospace Engineering Department, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
33
|
Melchionna S, Amati G, Bernaschi M, Bisson M, Succi S, Mitsouras D, Rybicki FJ. Risk assessment of atherosclerotic plaques based on global biomechanics. Med Eng Phys 2013; 35:1290-7; discussion 1290. [DOI: 10.1016/j.medengphy.2013.02.002] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2012] [Revised: 01/31/2013] [Accepted: 02/05/2013] [Indexed: 01/14/2023]
|
34
|
Russell JH, Kelson N, Barry M, Pearcy M, Fletcher DF, Winter CD. Computational Fluid Dynamic Analysis of Intracranial Aneurysmal Bleb Formation. Neurosurgery 2013; 73:1061-8; discussion 1068-9. [DOI: 10.1227/neu.0000000000000137] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022] Open
Abstract
Abstract
BACKGROUND:
The management of unruptured aneurysms is controversial, with the decision to treat influenced by aneurysm characteristics including size and morphology. Aneurysmal bleb formation is thought to be associated with an increased risk of rupture.
OBJECTIVE:
To correlate computational fluid dynamic (CFD) indices with bleb formation.
METHODS:
Anatomic models were constructed from 3-dimensional rotational angiography data in 27 patients with cerebral aneurysms harboring a single bleb. Additional models representing the aneurysm before bleb formation were constructed by digitally removing the bleb. We characterized hemodynamic features of models both with and without the blebs using CFDs. Flow structure, wall shear stress (WSS), pressure, and oscillatory shear index (OSI) were analyzed.
RESULTS:
There was a statistically significant association between bleb location at or adjacent to the point of maximal WSS (74%, P = .019), irrespective of rupture status. Aneurysmal blebs were related to the inflow or outflow jet in 89% of cases (P < .001), whereas 11% were unrelated. Maximal wall pressure and OSI were not significantly related to bleb location. The bleb region attained a lower WSS after its formation in 96% of cases (P < .001) and was also lower than the average aneurysm WSS in 86% of cases (P < .001).
CONCLUSION:
Cerebral aneurysm blebs generally form at or adjacent to the point of maximal WSS and are aligned with major flow structures. Wall pressure and OSI do not contribute to determining bleb location. The measurement of WSS using CFD models may potentially predict bleb formation and thus improve the assessment of rupture risk in unruptured aneurysms.
Collapse
Affiliation(s)
- Jeremy H. Russell
- Kenneth Jamieson Department of Neurosurgery, Royal Brisbane and Women's Hospital, Brisbane, Australia
| | - Neil Kelson
- Department of Information Technology Services (High Performance Computing)
| | - Mark Barry
- Department of Information Technology Services (High Performance Computing)
| | - Mark Pearcy
- School of Chemistry, Physics and Mechanical Engineering, Queensland University of Technology, Brisbane, Australia
| | - David F. Fletcher
- School of Chemical and Biomolecular Engineering, The University of Sydney, Sydney, Australia
| | - Craig D. Winter
- Kenneth Jamieson Department of Neurosurgery, Royal Brisbane and Women's Hospital, Brisbane, Australia
| |
Collapse
|
35
|
Castro MA. Understanding the role of hemodynamics in the initiation, progression, rupture, and treatment outcome of cerebral aneurysm from medical image-based computational studies. ISRN RADIOLOGY 2013; 2013:602707. [PMID: 24967285 PMCID: PMC4045510 DOI: 10.5402/2013/602707] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 05/15/2013] [Accepted: 06/19/2013] [Indexed: 12/31/2022]
Abstract
About a decade ago, the first image-based computational hemodynamic studies of cerebral aneurysms were presented. Their potential for clinical applications was the result of a right combination of medical image processing, vascular reconstruction, and grid generation techniques used to reconstruct personalized domains for computational fluid and solid dynamics solvers and data analysis and visualization techniques. A considerable number of studies have captivated the attention of clinicians, neurosurgeons, and neuroradiologists, who realized the ability of those tools to help in understanding the role played by hemodynamics in the natural history and management of intracranial aneurysms. This paper intends to summarize the most relevant results in the field reported during the last years.
Collapse
Affiliation(s)
- Marcelo A. Castro
- Grupo de Investigación y Desarrollo en Bioingeniería, Facultad Regional Buenos Aires, Universidad Tecnológica Nacional, CONICET, Medrano 951, CP 1179, Buenos Aires, Argentina
| |
Collapse
|
36
|
Chiastra C, Morlacchi S, Gallo D, Morbiducci U, Cárdenes R, Larrabide I, Migliavacca F. Computational fluid dynamic simulations of image-based stented coronary bifurcation models. J R Soc Interface 2013; 10:20130193. [PMID: 23676893 DOI: 10.1098/rsif.2013.0193] [Citation(s) in RCA: 86] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
One of the relevant phenomenon associated with in-stent restenosis in coronary arteries is an altered haemodynamics in the stented region. Computational fluid dynamics (CFD) offers the possibility to investigate the haemodynamics at a level of detail not always accessible within experimental techniques. CFD can quantify and correlate the local haemodynamics structures which might lead to in-stent restenosis. The aim of this work is to study the fluid dynamics of realistic stented coronary artery models which replicate the complete clinical procedure of stent implantation. Two cases of pathologic left anterior descending coronary arteries with their bifurcations are reconstructed from computed tomography angiography and conventional coronary angiography images. Results of wall shear stress and relative residence time show that the wall regions more prone to the risk of restenosis are located next to stent struts, to the bifurcations and to the stent overlapping zone for both investigated cases. Considering a bulk flow analysis, helical flow structures are generated by the curvature of the zone upstream from the stent and by the bifurcation regions. Helical recirculating microstructures are also visible downstream from the stent struts. This study demonstrates the feasibility to virtually investigate the haemodynamics of patient-specific coronary bifurcation geometries.
Collapse
Affiliation(s)
- Claudio Chiastra
- Chemistry, Materials and Chemical Engineering Department, Politecnico di Milano, Milan, Italy.
| | | | | | | | | | | | | |
Collapse
|
37
|
Papaharilaou Y, Aristokleous N, Seimenis I, Khozeymeh MI, Georgiou GC, Brott BC, Eracleous E, Anayiotos AS. Effect of head posture on the healthy human carotid bifurcation hemodynamics. Med Biol Eng Comput 2012; 51:207-18. [PMID: 23143389 DOI: 10.1007/s11517-012-0985-6] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2012] [Accepted: 10/29/2012] [Indexed: 11/25/2022]
Abstract
Head and neck postures may cause morphology changes to the geometry of the carotid bifurcation (CB) that alter the low and oscillating wall shear stress (WSS) regions previously reported as important in the development of atherosclerosis. Here the right and left CB were imaged by MRI in two healthy subjects in the neutral head posture with the subject in the supine position and in two other head postures with the subject in the prone position: (1) rightward rotation up to 80°, and (2) leftward rotation up to 80°. Image-based computational models were constructed to investigate the effect of posture on arterial geometry and local hemodynamics. The area exposure to unfavorable hemodynamics, based on thresholds set for oscillatory shear index (OSI), WSS and relative residence time, was used to quantify the hemodynamic impact on the wall. Torsion of the head was found to: (1) cause notable changes in the bifurcation and internal carotid artery angles and, in most cases, on cross-sectional area ratios for common, internal and external carotid artery, (2) change the spatial distribution of wall regions exposed to unfavorable hemodynamics, and (3) cause a marked change in the hemodynamic burden on the wall when the OSI was considered. These findings suggest that head posture may be associated with the genesis and development of atherosclerotic disease as well as complications in stenotic and stented vessels.
Collapse
Affiliation(s)
- Yannis Papaharilaou
- Institute of Applied and Computational Mathematics, Foundation for Research and Technology, Hellas, Heraklion, Crete, Greece
| | | | | | | | | | | | | | | |
Collapse
|
38
|
Sugiyama SI, Meng H, Funamoto K, Inoue T, Fujimura M, Nakayama T, Omodaka S, Shimizu H, Takahashi A, Tominaga T. Hemodynamic Analysis of Growing Intracranial Aneurysms Arising from a Posterior Inferior Cerebellar Artery. World Neurosurg 2012; 78:462-8. [DOI: 10.1016/j.wneu.2011.09.023] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/04/2011] [Accepted: 09/06/2011] [Indexed: 11/30/2022]
|
39
|
Capelli C, Biglino G, Petrini L, Migliavacca F, Cosentino D, Bonhoeffer P, Taylor AM, Schievano S. Finite Element Strategies to Satisfy Clinical and Engineering Requirements in the Field of Percutaneous Valves. Ann Biomed Eng 2012; 40:2663-73. [DOI: 10.1007/s10439-012-0617-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2012] [Accepted: 07/04/2012] [Indexed: 10/28/2022]
|
40
|
Renner J, Nadali Najafabadi H, Modin D, Länne T, Karlsson M. Subject-specific aortic wall shear stress estimations using semi-automatic segmentation. Clin Physiol Funct Imaging 2012; 32:481-91. [DOI: 10.1111/j.1475-097x.2012.01146.x] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2012] [Accepted: 05/13/2012] [Indexed: 11/28/2022]
Affiliation(s)
- J. Renner
- Department of Mechanical Engineering; Linköping University; Linköping; Sweden
| | | | - D. Modin
- Department of Medicine and Health Sciences; Linköping University; Linköping; Sweden
| | - T. Länne
- Department of Medicine and Health Sciences; Linköping University; Linköping; Sweden
| | - M. Karlsson
- Department of Mechanical Engineering; Linköping University; Linköping; Sweden
| |
Collapse
|
41
|
Sensitivity of the gradient oscillatory number to flow input waveform shapes. J Biomech 2012; 45:985-9. [DOI: 10.1016/j.jbiomech.2012.01.008] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2011] [Revised: 12/12/2011] [Accepted: 01/10/2012] [Indexed: 11/24/2022]
|
42
|
LaDisa JF, Alberto Figueroa C, Vignon-Clementel IE, Kim HJ, Xiao N, Ellwein LM, Chan FP, Feinstein JA, Taylor CA. Computational simulations for aortic coarctation: representative results from a sampling of patients. J Biomech Eng 2012; 133:091008. [PMID: 22010743 DOI: 10.1115/1.4004996] [Citation(s) in RCA: 85] [Impact Index Per Article: 7.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
Treatments for coarctation of the aorta (CoA) can alleviate blood pressure (BP) gradients (Δ), but long-term morbidity still exists that can be explained by altered indices of hemodynamics and biomechanics. We introduce a technique to increase our understanding of these indices for CoA under resting and nonresting conditions, quantify their contribution to morbidity, and evaluate treatment options. Patient-specific computational fluid dynamics (CFD) models were created from imaging and BP data for one normal and four CoA patients (moderate native CoA: Δ12 mmHg, severe native CoA: Δ25 mmHg and postoperative end-to-end and end-to-side patients: Δ0 mmHg). Simulations incorporated vessel deformation, downstream vascular resistance and compliance. Indices including cyclic strain, time-averaged wall shear stress (TAWSS), and oscillatory shear index (OSI) were quantified. Simulations replicated resting BP and blood flow data. BP during simulated exercise for the normal patient matched reported values. Greatest exercise-induced increases in systolic BP and mean and peak ΔBP occurred for the moderate native CoA patient (SBP: 115 to 154 mmHg; mean and peak ΔBP: 31 and 73 mmHg). Cyclic strain was elevated proximal to the coarctation for native CoA patients, but reduced throughout the aorta after treatment. A greater percentage of vessels was exposed to subnormal TAWSS or elevated OSI for CoA patients. Local patterns of these indices reported to correlate with atherosclerosis in normal patients were accentuated by CoA. These results apply CFD to a range of CoA patients for the first time and provide the foundation for future progress in this area.
Collapse
Affiliation(s)
- John F LaDisa
- Department of Biomedical Engineering, Marquette University, Milwaukee, WI 53233, USA.
| | | | | | | | | | | | | | | | | |
Collapse
|
43
|
Petersson S, Dyverfeldt P, Ebbers T. Assessment of the accuracy of MRI wall shear stress estimation using numerical simulations. J Magn Reson Imaging 2012; 36:128-38. [PMID: 22336966 DOI: 10.1002/jmri.23610] [Citation(s) in RCA: 97] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2011] [Accepted: 01/11/2012] [Indexed: 11/07/2022] Open
Abstract
PURPOSE To investigate the accuracy of wall shear stress (WSS) estimation using MRI. Specifically, to investigate the impact of different parameters and if MRI WSS estimates are monotonically related to actual WSS. MATERIALS AND METHODS The accuracy of WSS estimation using methods based on phase-contrast (PC) MRI velocity mapping, Fourier velocity encoding (FVE) and intravoxel velocity standard deviation mapping were studied using numerical simulations. The influence of spatial resolution, velocity encoding, wall segmentation, and voxel location were investigated over a range of WSS values. RESULTS WSS estimates were found to be sensitive to parameter settings in general and spatial resolution in particular. All methods underestimated WSS, except for the FVE-based method, which instead was extremely sensitive to voxel position relative to the wall. Methods using PC-based WSS estimation with wall segmentation showed to be accurate for low WSS, but were sensitive to segmentation errors. CONCLUSION Even in the absence of noise and for relatively simple velocity profiles, MRI WSS estimates cannot always be assumed to be linearly or even monotonically related to actual WSS. High WSS values cannot be resolved and the estimates depend on parameter setting. Nevertheless, distinguishing areas of low and moderate WSS may be feasible.
Collapse
Affiliation(s)
- Sven Petersson
- Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, Linköping, Sweden.
| | | | | |
Collapse
|
44
|
Xiong G, Choi G, Taylor CA. Virtual Interventions for Image-based Blood Flow Computation. COMPUTER AIDED DESIGN 2012; 44:3-14. [PMID: 22121255 PMCID: PMC3222596 DOI: 10.1016/j.cad.2011.01.004] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Image-based blood flow computation provides great promise for evaluation of vascular devices and assessment of surgical procedures. However, many previous studies employ idealized arterial and device models or only patient-specific models from the image data after device deployment, since the tools for model construction are unavailable or limited and tedious to use. Moreover, in contrast to retrospective studies from existing data, there is a pressing need for prospective analysis with the goal of surgical planning. Therefore, it is necessary to construct models with deployed devices in a fast, virtual and interactive fashion. The goal of this paper is to develop new geometric methods to deploy stents or stent grafts virtually to patient-specific geometric models constructed from a 3D segmentation of medical images. A triangular surface representing the vessel lumen boundary is extracted from the segmentation. The diseased portion is either clipped and replaced by the surface of a deployed device or rerouted in the case of a bypass graft. For diseased arteries close to bifurcations, bifurcated device models are generated. A method to map a 2D strut pattern on the surface of a device is also presented. We demonstrate three applications of our methods in personalized surgical planning for aortic aneurysms, aortic coarctation, and coronary artery stenosis using blood flow computation. Our approach enables prospective model construction and may help to expand the throughput required by routine clinical uses in the future.
Collapse
Affiliation(s)
- Guanglei Xiong
- Biomedical Informatics Program, Stanford University, Stanford, CA 94305, USA
| | - Gilwoo Choi
- Department of Mechanical Engineering, Stanford University, Stanford, CA 94305, USA
| | - Charles A. Taylor
- Departments of Bioengineering and Surgery, Stanford University, Stanford, CA 94305
| |
Collapse
|
45
|
Frydrychowicz A, Berger A, Munoz Del Rio A, Russe MF, Bock J, Harloff A, Markl M. Interdependencies of aortic arch secondary flow patterns, geometry, and age analysed by 4-dimensional phase contrast magnetic resonance imaging at 3 Tesla. Eur Radiol 2011; 22:1122-30. [PMID: 22207269 DOI: 10.1007/s00330-011-2353-6] [Citation(s) in RCA: 67] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2011] [Revised: 10/25/2011] [Accepted: 10/29/2011] [Indexed: 10/14/2022]
Abstract
OBJECTIVE It was the aim to analyse the impact of age, aortic arch geometry, and size on secondary flow patterns such as helix and vortex flow derived from flow-sensitive magnetic resonance imaging (4D PC-MRI). METHODS 62 subjects (age range = 20-80 years) without circumscribed pathologies of the thoracic aorta (ascending aortic (AAo) diameter: 3.2 ± 0.6 cm [range 2.2-5.1]) were examined by 4D PC-MRI after IRB-approval and written informed consent. Blood flow visualisation based on streamlines and time-resolved 3D particle traces was performed. Aortic diameter, shape (gothic, crook-shaped, cubic), angle, and age were correlated with existence and extent of secondary flow patterns (helicity, vortices); statistical modelling was performed. RESULTS Helical flow was the typical pattern in standard crook-shaped aortic arches. With altered shapes and increasing age, helicity was less common. AAo diameter and age had the highest correlation (r = 0.69 and 0.68, respectively) with number of detected vortices. None of the other arch geometric or demographic variables (for all, P ≥ 0.177) improved statistical modelling. CONCLUSION Substantially different secondary flow patterns can be observed in the normal thoracic aorta. Age and the AAo diameter were the parameters correlating best with presence and amount of vortices. Findings underline the importance of age- and geometry-matched control groups for haemodynamic studies. KEY POINTS • Secondary blood flow patterns (helices, vortices) are commonly observed in the aorta • Secondary flow patterns predominantly depend on patient age and aortic diameter • Geometric factors show a lesser impact on blood flow patterns than age and diameter • Future analyses of flow patterns should incorporate age- and diameter dependencies.
Collapse
Affiliation(s)
- Alex Frydrychowicz
- Clinic for Radiology and Nuclear Medicine, University Hospital Schleswig-Holstein, Campus Lübeck, Ratzeburger Allee 160, Haus 40, 23538 Lübeck, Germany.
| | | | | | | | | | | | | |
Collapse
|
46
|
Hunter KS, Lammers SR, Shandas R. Pulmonary vascular stiffness: measurement, modeling, and implications in normal and hypertensive pulmonary circulations. Compr Physiol 2011; 1:1413-35. [PMID: 23733649 PMCID: PMC4113421 DOI: 10.1002/cphy.c100005] [Citation(s) in RCA: 35] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
This article introduces the concept of pulmonary vascular stiffness, discusses its increasingly recognized importance as a diagnostic marker in the evaluation of pulmonary vascular disease, and describes methods to measure and model it clinically, experimentally, and computationally. It begins with a description of systems-level methods to evaluate pulmonary vascular compliance and recent clinical efforts in applying such techniques to better predict patient outcomes in pulmonary arterial hypertension. It then progresses from the systems-level to the local level, discusses proposed methods by which upstream pulmonary vessels increase in stiffness, introduces concepts around vascular mechanics, and concludes by describing recent work incorporating advanced numerical methods to more thoroughly evaluate changes in local mechanical properties of pulmonary arteries.
Collapse
Affiliation(s)
- Kendall S. Hunter
- Department of Bioengineering, University of Colorado at Denver Anschutz Medical Campus (UCD-AMC), Aurora, Colorado
- Division of Cardiology, Department of Pediatrics, The Children’s Hospital of Denver, UCD-AMC, Aurora, Colorado
| | - Steven R. Lammers
- Department of Bioengineering, University of Colorado at Denver Anschutz Medical Campus (UCD-AMC), Aurora, Colorado
- Cardiovascular Pulmonary (CVP) Research Laboratory, UCD-AMC, Aurora, Colorado
| | - Robin Shandas
- Department of Bioengineering, University of Colorado at Denver Anschutz Medical Campus (UCD-AMC), Aurora, Colorado
- Division of Cardiology, Department of Pediatrics, The Children’s Hospital of Denver, UCD-AMC, Aurora, Colorado
- Department of Surgery, UCD-AMC, Aurora, Colorado
| |
Collapse
|
47
|
Xiong G, Figueroa CA, Xiao N, Taylor CA. Simulation of blood flow in deformable vessels using subject-specific geometry and spatially varying wall properties. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2011; 27:1000-1016. [PMID: 21765984 PMCID: PMC3137382 DOI: 10.1002/cnm.1404] [Citation(s) in RCA: 41] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Simulation of blood flow using image-based models and computational fluid dynamics has found widespread application to quantifying hemodynamic factors relevant to the initiation and progression of cardiovascular diseases and for planning interventions. Methods for creating subject-specific geometric models from medical imaging data have improved substantially in the last decade but for many problems, still require significant user interaction. In addition, while fluid-structure interaction methods are being employed to model blood flow and vessel wall dynamics, tissue properties are often assumed to be uniform. In this paper, we propose a novel workflow for simulating blood flow using subject-specific geometry and spatially varying wall properties. The geometric model construction is based on 3D segmentation and geometric processing. Variable wall properties are assigned to the model based on combining centerline-based and surface-based methods. We finally demonstrate these new methods using an idealized cylindrical model and two subject-specific vascular models with thoracic and cerebral aneurysms.
Collapse
|
48
|
Shi Y, Lawford P, Hose R. Review of zero-D and 1-D models of blood flow in the cardiovascular system. Biomed Eng Online 2011; 10:33. [PMID: 21521508 PMCID: PMC3103466 DOI: 10.1186/1475-925x-10-33] [Citation(s) in RCA: 173] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2010] [Accepted: 04/26/2011] [Indexed: 11/16/2022] Open
Abstract
Background Zero-dimensional (lumped parameter) and one dimensional models, based on simplified representations of the components of the cardiovascular system, can contribute strongly to our understanding of circulatory physiology. Zero-D models provide a concise way to evaluate the haemodynamic interactions among the cardiovascular organs, whilst one-D (distributed parameter) models add the facility to represent efficiently the effects of pulse wave transmission in the arterial network at greatly reduced computational expense compared to higher dimensional computational fluid dynamics studies. There is extensive literature on both types of models. Method and Results The purpose of this review article is to summarise published 0D and 1D models of the cardiovascular system, to explore their limitations and range of application, and to provide an indication of the physiological phenomena that can be included in these representations. The review on 0D models collects together in one place a description of the range of models that have been used to describe the various characteristics of cardiovascular response, together with the factors that influence it. Such models generally feature the major components of the system, such as the heart, the heart valves and the vasculature. The models are categorised in terms of the features of the system that they are able to represent, their complexity and range of application: representations of effects including pressure-dependent vessel properties, interaction between the heart chambers, neuro-regulation and auto-regulation are explored. The examination on 1D models covers various methods for the assembly, discretisation and solution of the governing equations, in conjunction with a report of the definition and treatment of boundary conditions. Increasingly, 0D and 1D models are used in multi-scale models, in which their primary role is to provide boundary conditions for sophisticate, and often patient-specific, 2D and 3D models, and this application is also addressed. As an example of 0D cardiovascular modelling, a small selection of simple models have been represented in the CellML mark-up language and uploaded to the CellML model repository http://models.cellml.org/. They are freely available to the research and education communities. Conclusion Each published cardiovascular model has merit for particular applications. This review categorises 0D and 1D models, highlights their advantages and disadvantages, and thus provides guidance on the selection of models to assist various cardiovascular modelling studies. It also identifies directions for further development, as well as current challenges in the wider use of these models including service to represent boundary conditions for local 3D models and translation to clinical application.
Collapse
Affiliation(s)
- Yubing Shi
- Medical Physics Group, Department of Cardiovascular Science, Faculty of Medicine, Dentistry and Health, University of Sheffield, Sheffield S10 2RX, UK
| | | | | |
Collapse
|
49
|
Orlowski P, Al-Senani F, Summers P, Byrne J, Noble JA, Ventikos Y. Towards treatment planning for the embolization of arteriovenous malformations of the brain: intranidal hemodynamics modeling. IEEE Trans Biomed Eng 2011; 58:1994-2001. [PMID: 21349788 DOI: 10.1109/tbme.2011.2119317] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
Abstract
This paper presents a patient-derived model for the simulation of the hemodynamics of arteriovenous malformations of the brain (BAVM). This new approach is a step toward the simulation of the outcome of the embolization of the BAVM during treatment planning. More specifically, two aspects of the planning are pursued: simulation of the change of blood flow in the brain vasculature after the blocking of the malformation and simulation of the transport of the embolic liquid. The method we propose is tested on 3 BAVM cases of varying complexity. Twenty two out of 24 main BAVM flow paths have been identified well by simulation.
Collapse
Affiliation(s)
- Piotr Orlowski
- Institute of Biomedical Engineering, Department of Engineering Science, University of Oxford, Oxford, OX3 7DQ, UK
| | | | | | | | | | | |
Collapse
|
50
|
A Rapid and Computationally Inexpensive Method to Virtually Implant Current and Next-Generation Stents into Subject-Specific Computational Fluid Dynamics Models. Ann Biomed Eng 2011; 39:1423-37. [DOI: 10.1007/s10439-010-0238-5] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2010] [Accepted: 12/18/2010] [Indexed: 02/04/2023]
|