1
|
Yi W, Otani T, Yoshida T, Endo S, Wada S. Computational study on hemodynamic effects of left superior pulmonary vein resection and associated physiological changes in the left atrium after left upper lobectomy. Comput Methods Biomech Biomed Engin 2024; 27:167-178. [PMID: 36790387 DOI: 10.1080/10255842.2023.2178258] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2022] [Accepted: 01/27/2023] [Indexed: 02/16/2023]
Abstract
Left upper lobectomy (LUL) with left superior pulmonary vein (LSPV) resection alters the left atrium (LA) physiological states and LA hemodynamics associated with thrombosis, although this underlying mechanism is poorly understood. Therefore, we investigated the effects of LSPV resection and associated LA physiological changes on LA hemodynamics using four-dimensional computed tomography image-based computational simulations. Three cases were considered: the LA before and after LUL extracted from computed tomography images and artificial LSPV resection without physiological changes. Comparisons among the three cases demonstrated that physiological changes associated with LSPV resection are the possible factors that affect the LA hemodynamics after LUL.
Collapse
Affiliation(s)
- Wentao Yi
- Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, Japan
| | - Tomohiro Otani
- Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, Japan
| | - Takuya Yoshida
- Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, Japan
| | - Shunsuke Endo
- Saitama Medical Center, Jichi Medical University, Omiya, Saitama, Japan
| | - Shiego Wada
- Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka, Japan
| |
Collapse
|
2
|
Rizzo M, Ait-Ali L, Federici D, Festa P, Piagneri V, Landra F, Cameli M, Montesi G. Surgical repair of partial anomalous pulmonary venous connection in adulthood: A 4-dimensional flow magnetic resonance imaging postoperative evaluation. JTCVS Tech 2023; 22:208-211. [PMID: 38152192 PMCID: PMC10750877 DOI: 10.1016/j.xjtc.2023.09.017] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2023] [Revised: 08/23/2023] [Accepted: 09/16/2023] [Indexed: 12/29/2023] Open
Affiliation(s)
- Martina Rizzo
- Cardiac Surgery Unit, Department of Thoracic and Cardiovascular Disease, Santa Maria alle Scotte Hospital, University of Siena, Siena, Italy
| | - Lamia Ait-Ali
- Institute of Clinical Physiology, CNR, Monasterio Foundation, Massa, Italy
| | - Duccio Federici
- Paediatric Cardiac Surgery and GUCH Unit, Heart Hospital, Monasterio Foundation, Massa, Italy
| | - Pierluigi Festa
- Paediatric Cardiology and GUCH Unit, Heart Hospital, Monasterio Foundation, Massa, Italy
| | - Valeria Piagneri
- Diagnostic Imaging Unit, Heart Hospital, Monasterio Foundation, Massa, Italy
| | - Federico Landra
- Division of Cardiology, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Matteo Cameli
- Division of Cardiology, Department of Medical Biotechnologies, University of Siena, Siena, Italy
| | - Gianfranco Montesi
- Cardiac Surgery Unit, Department of Thoracic and Cardiovascular Disease, Santa Maria alle Scotte Hospital, Siena, Italy
| |
Collapse
|
3
|
Carpenter HJ, Ghayesh MH, Zander AC, Psaltis PJ. On the nonlinear relationship between wall shear stress topology and multi-directionality in coronary atherosclerosis. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2023; 231:107418. [PMID: 36842347 DOI: 10.1016/j.cmpb.2023.107418] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Revised: 02/01/2023] [Accepted: 02/09/2023] [Indexed: 06/18/2023]
Abstract
BACKGROUND AND OBJECTIVE In this paper we investigate twelve multi-directional/topological wall shear stress (WSS) derived metrics and their relationships with the formation of coronary plaques in both computational fluid dynamics (CFD) and dynamic fluid-structure interaction (FSI) frameworks. While low WSS is one of the most established biomechanical markers associated with coronary atherosclerosis progression, alone it is limited. Multi-directional and topological WSS derived metrics have been shown to be important in atherosclerosis related mechanotransduction and near-wall transport processes. However, the relationships between these twelve WSS metrics and the influence of both FSI simulations and coronary dynamics is understudied. METHODS We first investigate the relationships between these twelve WSS derived metrics, stenosis percentage and lesion length through a parametric, transient CFD study. Secondly, we extend the parametric study to FSI, both with and without the addition of coronary dynamics, and assess their correlations. Finally, we present the case of a patient who underwent invasive coronary angiography and optical coherence tomography imaging at two time points 18 months apart. Associations between each of the twelve WSS derived metrics in CFD, static FSI and dynamic FSI simulations were assessed against areas of positive/negative vessel remodelling, and changes in plaque morphology. RESULTS 22-32% stenosis was the threshold beyond which adverse multi-directional/topological WSS results. Each metric produced a different relationship with changing stenoses and lesion length. Transient haemodynamics was impacted by coronary dynamics, with the topological shear variation index suppressed by up to 94%. These changes appear more critical at smaller stenosis levels, suggesting coronary dynamics could play a role in the earlier stages of atherosclerosis development. In the patient case, both dynamics and FSI vs CFD changes altered associations with measured changes in plaque morphology. An appendix of the linear fits between the various FSI- and CFD-based simulations is provided to assist in scaling CFD-based results to resemble the compliant walled characteristics of FSI more accurately. CONCLUSIONS These results highlight the potential for coronary dynamics to alter multi-directional/topological WSS metrics which could impact associations with changes in coronary atherosclerosis over time. These results warrant further investigation in a wider range of morphological settings and longitudinal cohort studies in the future.
Collapse
Affiliation(s)
- Harry J Carpenter
- School of Mechanical Engineering, University of Adelaide, Adelaide, South Australia 5005, Australia.
| | - Mergen H Ghayesh
- School of Mechanical Engineering, University of Adelaide, Adelaide, South Australia 5005, Australia.
| | - Anthony C Zander
- School of Mechanical Engineering, University of Adelaide, Adelaide, South Australia 5005, Australia
| | - Peter J Psaltis
- Vascular Research Centre, Lifelong Health Theme, South Australian Health and Medical Research Institute (SAHMRI), Adelaide, South Australia 5000, Australia; Adelaide Medical School, University of Adelaide, Adelaide, South Australia 5005, Australia; Department of Cardiology, Central Adelaide Local Health Network, Adelaide, South Australia 5000, Australia
| |
Collapse
|
4
|
Bass DI, Marsh LMM, Fillingham P, Lim D, Chivukula VK, Kim LJ, Aliseda A, Levitt MR. Modeling the Mechanical Microenvironment of Coiled Cerebral Aneurysms. J Biomech Eng 2023; 145:041005. [PMID: 36193892 PMCID: PMC9791668 DOI: 10.1115/1.4055857] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Revised: 09/09/2022] [Indexed: 12/30/2022]
Abstract
Successful occlusion of cerebral aneurysms using coil embolization is contingent upon stable thrombus formation, and the quality of the thrombus depends upon the biomechanical environment. The goal of this study was to investigate how coil embolization alters the mechanical micro-environment within the aneurysm dome. Inertialess particles were injected in three-dimensional, computational simulations of flow inside patient aneurysms using patient-specific boundary conditions. Coil embolization was simulated as a homogenous porous medium of known permeability and inertial constant. Lagrangian particle tracking was used to calculate the residence time and shear stress history for particles in the flow before and after treatment. The percentage of particles entering the aneurysm dome correlated with the neck surface area before and after treatment (pretreatment: R2 = 0.831, P < 0.001; post-treatment: R2 = 0.638, P < 0.001). There was an inverse relationship between the change in particles entering the dome and coil packing density (R2 = 0.600, P < 0.001). Following treatment, the particles with the longest residence times tended to remain within the dome even longer while accumulating lower shear stress. A significant correlation was observed between the treatment effect on residence time and the ratio of the neck surface area to porosity (R2 = 0.390, P = 0.007). The results of this study suggest that coil embolization triggers clot formation within the aneurysm dome via a low shear stress-mediated pathway. This hypothesis links independently observed findings from several benchtop and clinical studies, furthering our understanding of this treatment strategy.
Collapse
Affiliation(s)
- David I. Bass
- Department of Neurological Surgery, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104
| | - Laurel M. M. Marsh
- Department of Mechanical Engineering, University of Washington, 3900 East Stevens Way NE, Box 352600, Seattle, WA 98195
| | - Patrick Fillingham
- Department of Neurological Surgery, Stroke & Applied Neuroscience Center, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104
| | - Do Lim
- Department of Neurological Surgery, Stroke & Applied Neuroscience Center, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104
| | - V. Keshav Chivukula
- Department of Biomedical and Chemical Engineering and Sciences, Florida Institute of Technology, 150 West University Building, Melbourne, FL 32901
| | - Louis J. Kim
- Department of Neurological Surgery, Stroke & Applied Neuroscience Center, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104; Department of Radiology, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104
| | - Alberto Aliseda
- Department of Mechanical Engineering, Stroke & Applied Neuroscience Center, University of Washington, 3900 East Stevens Way NE, Box 352600, Seattle, WA 98195; Department of Neurological Surgery, Stroke & Applied Neuroscience Center, University of Washington, 3900 East Stevens Way NE, Box 352600, Seattle, WA 98195
| | - Michael R. Levitt
- Department of Neurological Surgery, Stroke & Applied Neuroscience Center, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104; Department of Radiology, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104; Department of Mechanical Engineering, University of Washington, 325 9th Avenue, Box 359924, Seattle, WA 98104
| |
Collapse
|
5
|
Flow residence time in intracranial aneurysms evaluated by in vitro 4D flow MRI. J Biomech 2022; 141:111211. [DOI: 10.1016/j.jbiomech.2022.111211] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2021] [Revised: 05/15/2022] [Accepted: 06/24/2022] [Indexed: 11/15/2022]
|
6
|
Abstract
Computational fluid dynamics (CFD) modeling of blood flow plays an important role in better understanding various medical conditions, designing more effective drug delivery systems, and developing novel diagnostic methods and treatments. However, despite significant advances in computational technology and resources, the expensive computational cost of these simulations still hinders their transformation from a research interest to a clinical tool. This bottleneck is even more severe for image-based, patient-specific CFD simulations with realistic boundary conditions and complex computational domains, which make such simulations excessively expensive. To address this issue, deep learning approaches have been recently explored to accelerate computational hemodynamics simulations. In this study, we review recent efforts to integrate deep learning with CFD and discuss the applications of this approach in solving hemodynamics problems, such as blood flow behavior in aorta and cerebral arteries. We also discuss potential future directions in the field. In this review, we suggest that incorporating physiologic understandings and underlying fluid mechanics laws in deep learning models will soon lead to a paradigm shift in the development novel non-invasive computational medical decisions.
Collapse
|
7
|
Mahmoudi M, Jennings C, Pereira K, Hall AF, Arzani A. Guiding the prostatic artery embolization procedure with computational uid dynamics. J Biomech Eng 2022; 144:1140781. [DOI: 10.1115/1.4054515] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2021] [Indexed: 11/08/2022]
Abstract
Abstract
Benign prostatic hyperplasia (BPH) is a common disease associated with lower urinary tract symptoms and the most frequent benign tumor in men. To reduce BPH therapy complications, prostatic artery embolization (PAE) was developed to replace the surgical options. PAE is a minimally invasive technique in which emboli are injected into the prostate arteries (PA), obstructing the blood flow in the hypervascular nodules. In this work, a personalized PAE treatment strategy was proposed using patient-specific computational fluid dynamics (CFD). First, the hemodynamics environment in the iliac arterial tree considering a large network of bifurcations was studied. The results showed complex blood flow patterns in the iliac arterial network. Subsequently, the transport of embolic particulates during PAE for the standard horizontal and a hypothetical vertical patient positioning was simulated using Lagrangian particle tracking. Emboli with different sizes were released at various locations across the iliac arterial tree. The emboli entering the PA were mapped back to their initial location to create emboli release maps (ERMs). The obtained ERMs during the standard patient positioning for smaller emboli at certain release locations showed distinct regions in which if the emboli were released within these regions, all of them would reach the PA without non-target embolization. During the hypothetical vertical patient positioning, the larger emboli formed a larger coherent region in the ERMs. Our patient-specific model can be used to find the best spatial location for emboli injection and perform the embolization procedure with minimal off-target delivery.
Collapse
Affiliation(s)
- Mostafa Mahmoudi
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, United States
| | - Chadrick Jennings
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, United States
| | - Keith Pereira
- Department of Radiology, Saint Louis University, St Louis, MO, United States
| | - Andrew F. Hall
- Department of Biomedical Engineering, Saint Louis University, St Louis, MO, United States
| | - Amirhossein Arzani
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, United States
| |
Collapse
|
8
|
Otani T, Yoshida T, Yi W, Endo S, Wada S. On the Impact of Left Upper Lobectomy on the Left Atrial Hemodynamics. Front Physiol 2022; 13:830436. [PMID: 35283800 PMCID: PMC8908205 DOI: 10.3389/fphys.2022.830436] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 01/27/2022] [Indexed: 11/13/2022] Open
Abstract
The left atrium (LA) functions to transport oxygenated blood from the pulmonary veins (PVs) to the left ventricle (LV). LA hemodynamics has received much attention because thrombosis in the LA in pathological states, such as atrial fibrillation, is a major factor leading to thromboembolic stroke. In the last 5 years, multiple cohort studies have revealed that left upper lobectomy (LUL) with PV resection risks thrombus formation in the PV stump even in the normal LA without a history of cardiac disease; the causal mechanism is, however, an open question. The present study investigated the potential effect of an LUL on LA hemodynamics associated with thrombus formation through computational simulation using four-dimensional computed tomography (4D-CT) images. Time series of patient-specific LA geometries before and after LUL were extracted from the 4D-CT images and these motions were estimated through non-rigid registration. Adopting the LA geometries and prescribed moving wall boundary conditions, the LA blood flow was determined using a Cartesian-grid computational fluid dynamics solver. The obtained results show that the LUL resulted in blood flow impingement from the left and right PV inflows into the LA upper region throughout most of the cardiac cycle. This characteristic alteration of the LA hemodynamics generated fine-scale vortices with viscous energy dissipations, enhancing the flow stasis associated with thrombus formation in the PV stump. These findings show that an LUL affects the hemodynamics not only in the PV stump but also throughout the LA region. They also highlight the importance of computational analysis of LA hemodynamics in understanding the underlying mechanism of LUL-induced thrombus formation.
Collapse
Affiliation(s)
- Tomohiro Otani
- Graduate School of Engineering Science, Osaka University, Osaka, Japan
- *Correspondence: Tomohiro Otani
| | - Takuya Yoshida
- Graduate School of Engineering Science, Osaka University, Osaka, Japan
| | - Wentao Yi
- Graduate School of Engineering Science, Osaka University, Osaka, Japan
| | - Shunsuke Endo
- Jichi Medical University Saitama Medical Center, Saitama, Japan
| | - Shigeo Wada
- Graduate School of Engineering Science, Osaka University, Osaka, Japan
| |
Collapse
|
9
|
Tercanlı MF, Olcay AB, Mutlu O, Bilgin C, Hakyemez B. Investigation of the effect of anticoagulant usage in the flow diverter stent treatment of the patient-specific cerebral aneurysm using the Lagrangian coherent structures. J Clin Neurosci 2021; 94:86-93. [PMID: 34863468 DOI: 10.1016/j.jocn.2021.10.007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2020] [Revised: 08/13/2021] [Accepted: 10/05/2021] [Indexed: 10/20/2022]
Abstract
Anticoagulants are prescribed to the flow diverter treated patients to diminish the risk of embolism in the arteries. In the present study, digital subtraction angiography images of a 49-year-old female patient with a left paraophthalmic aneurysm were used to build a numerical model to investigate the effect of an anticoagulant on hemodynamics at the aneurysm site. The Carreau-Yasuda viscosity model was utilized to define blood viscosity, and the coefficients of the viscosity model were updated based on the usage of warfarin. The five-cardiac cycle-long numerical simulations were performed, and Lagrangian coherent structures, hyperbolic time, and fluid particle analyses were also employed in the numerical models. These analyses allowed us to evaluate the formation of stagnated regions, recirculation zones, and the number of jailed particles inside the aneurysm sac following a flow diverter placement. It is realized that anticoagulant use caused blood to be less viscous, yielding a substantial amount of incoming blood flow to enter the aneurysm sac. Only 12% of the nearly 25,000 fluid particles seeded from the artery inlet have stayed inside the sac. Furthermore, the deviation between warfarin added blood and normal blood flow becomes more extensive, with every heartbeat undermining the effectiveness of patient-specific CFD models when the use of anticoagulants is overlooked in the viscosity models.
Collapse
Affiliation(s)
| | - Ali Bahadır Olcay
- Yeditepe University, Faculty of Engineering, Department of Mechanical Engineering, Kayisdagi Cad., 34755 Istanbul, Turkey.
| | - Onur Mutlu
- Biomedical Research Center, Qatar University, Doha P.O. Box 2713, Qatar
| | - Cem Bilgin
- University of Health Sciences, Bursa Yuksek Ihtisas Training and Research Hospital, Department of Radiology, Yildirim, Bursa 16310, Turkey
| | - Bahattin Hakyemez
- Uludag University School of Medicine, Department of Radiology, Gorukle, Bursa 16059, Turkey
| |
Collapse
|
10
|
Darwish A, Norouzi S, Kadem L. Spectral-Clustering of Lagrangian Trajectory Graphs: Application to Abdominal Aortic Aneurysms. Cardiovasc Eng Technol 2021; 13:504-513. [PMID: 34845627 DOI: 10.1007/s13239-021-00590-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 06/09/2021] [Accepted: 11/01/2021] [Indexed: 11/29/2022]
Abstract
PURPOSE Identification of coherent structures in cardiovascular flows is crucial to describe the transport and mixing of blood. Coherent structures can highlight locations where minimal blood mixing takes place, thus, potential thrombus formation can be expected thither. Graph-based approaches have recently been introduced in order to describe fluid transport and mixing between multiple Lagrangian trajectories, where each trajectory serves as a node that can be connected to another trajectory based on their relative distance during the course of time. METHODS In this study, we compute the Lagrangian trajectories from in vitro planar instantaneous velocity fields in two models of abdominal aortic aneurysms, (AAA) namely single bulge and bi-lobed. Then, we construct unweighted and undirected graphs based on the pairwise distance of Lagrangian trajectories. We report local measures of the graph namely the degree and the clustering coefficient. We also perform spectral clustering of the graph Laplacian to extract the flow coherent sets. RESULTS Local graph measures reveal fluid regions of high mixing such as vortex boundaries. Through spectral clustering, the fluid is partitioned into a reduced number of coherent sets where within each set, inner mixing of fluid is maximized while the fluid mixing between different coherent sets is minimized. The approach reveals multiple coherent sets adjacent to the AAA bulge that have sustained this adjacency to the wall through their coherent motion during one cardiac cycle. CONCLUSION Identifying coherent sets enables tracking their transport during the cardiac cycle and identify their role in the flow dynamics. Moreover, the size and the transport of the long residing coherent sets inside the AAA bulges can be deduced which may aid in predicting thrombus formation at such location.
Collapse
Affiliation(s)
- Ahmed Darwish
- Laboratory of Cardiovascular Fluid Dynamics, Concordia University, Montréal, QC, H3G 1M8, Canada. .,Mechanical Engineering Department, Assiut University, 71515, Assiut, Egypt.
| | - Shahrzad Norouzi
- Laboratory of Cardiovascular Fluid Dynamics, Concordia University, Montréal, QC, H3G 1M8, Canada
| | - Lyes Kadem
- Laboratory of Cardiovascular Fluid Dynamics, Concordia University, Montréal, QC, H3G 1M8, Canada
| |
Collapse
|
11
|
In Silico Hemodynamics and Filtering Evaluation of a Commercial Embolic Protection Device. Ann Biomed Eng 2021; 49:2659-2670. [PMID: 34405319 DOI: 10.1007/s10439-021-02846-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/08/2021] [Accepted: 08/01/2021] [Indexed: 10/20/2022]
Abstract
During the last years, several kinds of Embolic Protection Devices (EPD) have been developed, with the aim of minimizing complication caused by thrombi generated during Carotid Artery Stenting (CAS). These devices are capable of capturing small particles generated during the intervention, avoiding cerebral stroke and improving the outcomes of the surgery. However, they have associated complications, like the increase on flow resistance associated by their use or the lack of knowledge on their actual filtration efficiency for thrombi of low size. Current work proposes a validated computational methodology in order to predict the hemodynamic features and filtering efficiency of a commercial EPD. It will be observed how Computational Fluid Dynamics predicts pressure drop with fair agreement with the experimental measurements. Finally, this work analyzes the filtration efficiency and the influence of the distribution of injected particles on this parameter. The capabilities of the filter for retaining particles of diameter below the pore size is, additionally, discussed.
Collapse
|
12
|
Linden NJ, Tabuena DR, Steinmetz NA, Moody WJ, Brunton SL, Brunton BW. Go with the FLOW: visualizing spatiotemporal dynamics in optical widefield calcium imaging. J R Soc Interface 2021; 18:20210523. [PMID: 34428947 PMCID: PMC8385384 DOI: 10.1098/rsif.2021.0523] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Accepted: 08/02/2021] [Indexed: 12/11/2022] Open
Abstract
Widefield calcium imaging has recently emerged as a powerful experimental technique to record coordinated large-scale brain activity. These measurements present a unique opportunity to characterize spatiotemporally coherent structures that underlie neural activity across many regions of the brain. In this work, we leverage analytic techniques from fluid dynamics to develop a visualization framework that highlights features of flow across the cortex, mapping wavefronts that may be correlated with behavioural events. First, we transform the time series of widefield calcium images into time-varying vector fields using optic flow. Next, we extract concise diagrams summarizing the dynamics, which we refer to as FLOW (flow lines in optical widefield imaging) portraits. These FLOW portraits provide an intuitive map of dynamic calcium activity, including regions of initiation and termination, as well as the direction and extent of activity spread. To extract these structures, we use the finite-time Lyapunov exponent technique developed to analyse time-varying manifolds in unsteady fluids. Importantly, our approach captures coherent structures that are poorly represented by traditional modal decomposition techniques. We demonstrate the application of FLOW portraits on three simple synthetic datasets and two widefield calcium imaging datasets, including cortical waves in the developing mouse and spontaneous cortical activity in an adult mouse.
Collapse
Affiliation(s)
- Nathaniel J. Linden
- Department of Bioengineering, University of Washington, Seattle, WA 98195-0005, USA
- Department of Biology, University of Washington, Seattle, WA 98195-0005, USA
| | - Dennis R. Tabuena
- Department of Biology, University of Washington, Seattle, WA 98195-0005, USA
- Graduate Program in Neuroscience, University of Washington, Seattle, WA 98195-0005, USA
| | - Nicholas A. Steinmetz
- Department of Biological Structure, University of Washington, Seattle, WA 98195-0005, USA
| | - William J. Moody
- Department of Biology, University of Washington, Seattle, WA 98195-0005, USA
| | - Steven L. Brunton
- Department of Mechanical Engineering, University of Washington, Seattle, WA 98195-0005, USA
| | - Bingni W. Brunton
- Department of Biology, University of Washington, Seattle, WA 98195-0005, USA
| |
Collapse
|
13
|
Engineering Perspective on Cardiovascular Simulations of Fontan Hemodynamics: Where Do We Stand with a Look Towards Clinical Application. Cardiovasc Eng Technol 2021; 12:618-630. [PMID: 34114202 DOI: 10.1007/s13239-021-00541-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 09/20/2020] [Accepted: 04/30/2021] [Indexed: 01/02/2023]
Abstract
BACKGROUND Cardiovascular simulations for patients with single ventricles undergoing the Fontan procedure can assess patient-specific hemodynamics, explore surgical advances, and develop personalized strategies for surgery and patient care. These simulations have not yet been broadly accepted as a routine clinical tool owing to a number of limitations. Numerous approaches have been explored to seek innovative solutions for improving methodologies and eliminating these limitations. PURPOSE This article first reviews the current state of cardiovascular simulations of Fontan hemodynamics. Then, it will discuss the technical progress of Fontan simulations with the emphasis of its clinical impact, noting that substantial improvements have been made in the considerations of patient-specific anatomy, flow, and blood rheology. The article concludes with insights into potential future directions involving clinical validation, uncertainty quantification, and computational efficiency. The advancements in these aspects could promote the clinical usage of Fontan simulations, facilitating its integration into routine clinical practice.
Collapse
|
14
|
Fluid Flow Characteristics of Healthy and Calcified Aortic Valves Using Three-Dimensional Lagrangian Coherent Structures Analysis. FLUIDS 2021. [DOI: 10.3390/fluids6060203] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Aortic valve calcification is an important cardiovascular disorder that deteriorates the accurate functioning of the valve leaflets. The increasing stiffness due to the calcification prevents the complete closure of the valve and therefore leads to significant hemodynamic alterations. Computational fluid dynamics (CFD) modeling enables the investigation of the entire flow domain by processing medical images from aortic valve patients. In this study, we computationally modeled and simulated a 3D aortic valve using patient-specific dimensions of the aortic root and aortic sinus. Leaflet stiffness is deteriorated in aortic valve disease due to calcification. In order to investigate the influence of leaflet calcification on flow dynamics, three different leaflet-stiffness values were considered for healthy, mildly calcified, and severely calcified leaflets. Time-dependent CFD results were used for applying the Lagrangian coherent structures (LCS) technique by performing finite-time Lyapunov exponent (FTLE) computations along with Lagrangian particle residence time (PRT) analysis to identify unique vortex structures at the front and backside of the leaflets. Obtained results indicated that the peak flow velocity at the valve orifice increased with the calcification rate. For the healthy aortic valve, a low-pressure field was observed at the leaflet tips. This low-pressure field gradually expanded through the entire aortic sinus as the calcification level increased. FTLE field plots of the healthy and calcified valves showed a variety of differences in terms of flow structures. When the number of fluid particles in the healthy valve model was taken as reference, 1.59 and 1.74 times more particles accumulated in the mildly and severely calcified valves, respectively, indicating that the calcified valves were not sufficiently opened to allow normal mass flow rates.
Collapse
|
15
|
The association between hemodynamics and wall characteristics in human intracranial aneurysms: a review. Neurosurg Rev 2021; 45:49-61. [PMID: 33913050 DOI: 10.1007/s10143-021-01554-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2021] [Revised: 04/02/2021] [Accepted: 04/20/2021] [Indexed: 12/28/2022]
Abstract
Hemodynamics plays a key role in the natural history of intracranial aneurysms (IAs). However, studies exploring the association between aneurysmal hemodynamics and the biological and mechanical characteristics of the IA wall in humans are sparse. In this review, we survey the current body of literature, summarize the studies' methodologies and findings, and assess the degree of consensus among them. We used PubMed to perform a systematic review of studies that explored the association between hemodynamics and human IA wall features using different sources. We identified 28 publications characterizing aneurysmal flow and the IA wall: 4 using resected tissues, 17 using intraoperative images, and 7 using vessel wall magnetic resonance imaging (MRI). Based on correlation to IA tissue, higher flow conditions, such as high wall shear stress (WSS) with complex pattern and elevated pressure, were associated with degenerated walls and collagens with unphysiological orientation and faster synthesis. MRI studies strongly supported that low flow, characterized by low WSS and high blood residence time, was associated with thicker walls and post-contrast enhancement. While significant discrepancies were found among those utilized intraoperative images, they generally supported that thicker walls coexist at regions with prolonged residence time and that thinner regions are mainly exposed to higher pressure with complex WSS patterns. The current body of literature supports a theory of two general hemodynamic-biologic mechanisms for IA development. One, where low flow conditions are associated with thickening and atherosclerotic-like remodeling, and the other where high and impinging flow conditions are related to wall degeneration, thinning, and collagen remodeling.
Collapse
|
16
|
Wall Shear Stress Topological Skeleton Analysis in Cardiovascular Flows: Methods and Applications. MATHEMATICS 2021. [DOI: 10.3390/math9070720] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A marked interest has recently emerged regarding the analysis of the wall shear stress (WSS) vector field topological skeleton in cardiovascular flows. Based on dynamical system theory, the WSS topological skeleton is composed of fixed points, i.e., focal points where WSS locally vanishes, and unstable/stable manifolds, consisting of contraction/expansion regions linking fixed points. Such an interest arises from its ability to reflect the presence of near-wall hemodynamic features associated with the onset and progression of vascular diseases. Over the years, Lagrangian-based and Eulerian-based post-processing techniques have been proposed aiming at identifying the topological skeleton features of the WSS. Here, the theoretical and methodological bases supporting the Lagrangian- and Eulerian-based methods currently used in the literature are reported and discussed, highlighting their application to cardiovascular flows. The final aim is to promote the use of WSS topological skeleton analysis in hemodynamic applications and to encourage its application in future mechanobiology studies in order to increase the chance of elucidating the mechanistic links between blood flow disturbances, vascular disease, and clinical observations.
Collapse
|
17
|
Meschi SS, Farghadan A, Arzani A. Flow topology and targeted drug delivery in cardiovascular disease. J Biomech 2021; 119:110307. [PMID: 33676269 DOI: 10.1016/j.jbiomech.2021.110307] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Accepted: 02/03/2021] [Indexed: 12/13/2022]
Abstract
Targeted drug delivery is a promising technique to direct the drug to the specific diseased region. Nanoparticles have provided an attractive approach for this purpose. In practice, the major focus of targeted delivery has been on targeting cell receptors. However, the complex fluid mechanics in diseased biomedical flows questions if a sufficient number of nanoparticles can reach the desired region. In this paper, we propose that hidden topological structures in cardiovascular flows identified with Lagrangian coherent structures (LCS) control drug transport and provide valuable information for optimizing targeted drug delivery efficiency. We couple image-based computational fluid dynamics (CFD) with continuum transport models to study nanoparticle transport in coronary artery disease. We simulate nanoparticle transport as well as the recently proposed shear targeted drug delivery system that couples micro-carriers with nanoparticle drugs. The role of the LCS formed near the stenosed artery in controlling drug transport is discussed. Our results motivate the design of smart micro-needles guided by flow topology, which could achieve optimal drug delivery efficiency.
Collapse
Affiliation(s)
- Sara S Meschi
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, USA
| | - Ali Farghadan
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, USA; Department of Mechanical Engineering, University of Michigan, Ann Arbor, MI, USA
| | - Amirhossein Arzani
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, USA.
| |
Collapse
|
18
|
Plitman Mayo R, Yaakobovich H, Finkelstein A, Shadden SC, Marom G. Impact of BASILICA on the thrombogenicity potential of valve-in-valve implantations. J Biomech 2021; 118:110309. [PMID: 33601181 DOI: 10.1016/j.jbiomech.2021.110309] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Accepted: 02/03/2021] [Indexed: 01/17/2023]
Abstract
Subclinical leaflet thrombosis is becoming a major concern in valve-in-valve procedures, whereby a transcatheter aortic valve device is deployed inside a failed bioprosthetic surgical valve. Blood flow stagnation and prolonged residence times in the neo-sinuses have been suggested as possible explanations for leaflet thrombosis. The BASILICA technique, which was originally developed to treat coronary flow obstruction, has also been proposed as an alternative to reduce the risk of thrombus formation. The aim of this study is to understand the impact of BASILICA on the valve-in-valve thrombogenicity using computational fluid dynamics simulations. To this end, two Eulerian and two Lagrangian approaches were employed to estimate near-wall stagnation measures in eight valve-in-valve models. The models included an intact or lacerated Sorin Mitroflow surgical valve, and either a SAPIEN or Evolut transcatheter aortic valve device. The Lagrangian approaches predicted a high number of particles and vortices concentration in the proximal areas of the neo-sinuses, while the Eulerian approaches did so in the distal areas. As a consequence, this study demonstrated that Lagrangian approaches are better predictors of subclinical leaflet thrombosis, since they match experimental and clinical findings. Additionally, the SAPIEN valve possess a higher risk for developing leaflet thrombosis, and two lacerations are shown to provide the best results in terms of development of vortices and accumulation of particles within the neo-sinuses. This study highlights the potential of computational modeling in aiding clinicians in their decision-making in valve-in-valve and BASILICA procedures.
Collapse
Affiliation(s)
| | - Halit Yaakobovich
- Department of Biomedical Engineering, Tel Aviv University, Tel Aviv, Israel
| | | | - Shawn C Shadden
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | - Gil Marom
- School of Mechanical Engineering, Tel Aviv University, Tel Aviv, Israel.
| |
Collapse
|
19
|
Mirramezani M, Shadden SC. A Distributed Lumped Parameter Model of Blood Flow. Ann Biomed Eng 2020; 48:2870-2886. [PMID: 32613457 PMCID: PMC7725998 DOI: 10.1007/s10439-020-02545-6] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2020] [Accepted: 06/03/2020] [Indexed: 01/02/2023]
Abstract
We propose a distributed lumped parameter (DLP) modeling framework to efficiently compute blood flow and pressure in vascular domains. This is achieved by developing analytical expressions describing expected energy losses along vascular segments, including from viscous dissipation, unsteadiness, flow separation, vessel curvature and vessel bifurcations. We apply this methodology to solve for unsteady blood flow and pressure in a variety of complex 3D image-based vascular geometries, which are typically approached using computational fluid dynamics (CFD) simulations. The proposed DLP framework demonstrated consistent agreement with CFD simulations in terms of flow rate and pressure distribution, with mean errors less than 7% over a broad range of hemodynamic conditions and vascular geometries. The computational cost of the DLP framework is orders of magnitude lower than the computational cost of CFD, which opens new possibilities for hemodynamics modeling in timely decision support scenarios, and a multitude of applications of imaged-based modeling that require ensembles of numerical simulations.
Collapse
Affiliation(s)
- Mehran Mirramezani
- Mechanical Engineering, University of California, Berkeley, CA, 94720, USA
- Mathematics, University of California, Berkeley, CA, 94720, USA
| | - Shawn C Shadden
- Mechanical Engineering, University of California, Berkeley, CA, 94720, USA.
| |
Collapse
|
20
|
Plitman Mayo R, Yaakobovich H, Finkelstein A, Shadden SC, Marom G. Numerical models for assessing the risk of leaflet thrombosis post-transcatheter aortic valve-in-valve implantation. ROYAL SOCIETY OPEN SCIENCE 2020; 7:201838. [PMID: 33489295 PMCID: PMC7813235 DOI: 10.1098/rsos.201838] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2020] [Accepted: 11/20/2020] [Indexed: 05/04/2023]
Abstract
Leaflet thrombosis has been suggested as the reason for the reduced leaflet motion in cases of hypoattenuated leaflet thickening of bioprosthetic aortic valves. This work aimed to estimate the risk of leaflet thrombosis in two post-valve-in-valve (ViV) configurations, using five different numerical approaches. Realistic ViV configurations were calculated by modelling the deployments of the latest version of transcatheter aortic valve devices (Medtronic Evolut PRO, Edwards SAPIEN 3) in the surgical Sorin Mitroflow. Computational fluid dynamics simulations of blood flow followed the dry models. Lagrangian and Eulerian measures of near-wall stagnation were implemented by particle and concentration tracking, respectively, to estimate the thrombogenicity and to predict the risk locations. Most of the numerical approaches indicate a higher leaflet thrombosis risk in the Edwards SAPIEN 3 device because of its intra-annular implantation. The Eulerian approaches estimated high-risk locations in agreement with the wall sheer stress (WSS) separation points. On the other hand, the Lagrangian approaches predicted high-risk locations at the proximal regions of the leaflets matching the low WSS magnitude regions of both transcatheter aortic valve implantation models and reported clinical and experimental data. The proposed methods can help optimizing future designs of transcatheter aortic valves with minimal thrombotic risks.
Collapse
Affiliation(s)
- Romina Plitman Mayo
- School of Mechanical Engineering, Tel Aviv University, Tel Aviv, Israel
- Department of Biological Regulation, Weizmann Institute of Science, Herzl Street 234, Rehovot, Israel
| | - Halit Yaakobovich
- Department of Biomedical Engineering, Tel Aviv University, Tel Aviv, Israel
| | | | - Shawn C. Shadden
- Department of Mechanical Engineering, University of California, Berkeley, CA, USA
| | - Gil Marom
- School of Mechanical Engineering, Tel Aviv University, Tel Aviv, Israel
| |
Collapse
|
21
|
Fluid-Structure Interaction Simulation of an Intra-Atrial Fontan Connection. BIOLOGY 2020; 9:biology9120412. [PMID: 33255292 PMCID: PMC7760396 DOI: 10.3390/biology9120412] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 10/11/2020] [Revised: 11/19/2020] [Accepted: 11/20/2020] [Indexed: 12/15/2022]
Abstract
Simple Summary A fluid-structure interaction (FSI) simulation of an intra-atrial Fontan connection was performed. Power loss and pressure drop results fluctuated less during the FSI simulation than during the simulation run with rigid walls, but there were no observable differences in time-averaged pressure drop, connection power loss or hepatic flow distribution. These results suggested that employing a rigid wall is a reasonable assumption when evaluating time-averaged hemodynamic quantities of the Fontan connection under resting breath-held flow conditions. Abstract Total cavopulmonary connection (TCPC) hemodynamics has been hypothesized to be associated with long-term complications in single ventricle heart defect patients. Rigid wall assumption has been commonly used when evaluating TCPC hemodynamics using computational fluid dynamics (CFD) simulation. Previous study has evaluated impact of wall compliance on extra-cardiac TCPC hemodynamics using fluid-structure interaction (FSI) simulation. However, the impact of ignoring wall compliance on the presumably more compliant intra-atrial TCPC hemodynamics is not fully understood. To narrow this knowledge gap, this study aims to investigate impact of wall compliance on an intra-atrial TCPC hemodynamics. A patient-specific model of an intra-atrial TCPC is simulated with an FSI model. Patient-specific 3D TCPC anatomies were reconstructed from transverse cardiovascular magnetic resonance images. Patient-specific vessel flow rate from phase-contrast magnetic resonance imaging (MRI) at the Fontan pathway and the superior vena cava under resting condition were prescribed at the inlets. From the FSI simulation, the degree of wall deformation was compared with in vivo wall deformation from phase-contrast MRI data as validation of the FSI model. Then, TCPC flow structure, power loss and hepatic flow distribution (HFD) were compared between rigid wall and FSI simulation. There were differences in instantaneous pressure drop, power loss and HFD between rigid wall and FSI simulations, but no difference in the time-averaged quantities. The findings of this study support the use of a rigid wall assumption on evaluation of time-averaged intra-atrial TCPC hemodynamic metric under resting breath-held condition.
Collapse
|
22
|
Ameenuddin M, Anand M. A Mixture Theory Model for Blood Combined With Low-Density Lipoprotein Transport to Predict Early Atherosclerosis Regions in Idealized and Patient-Derived Abdominal Aorta. J Biomech Eng 2020; 142:101008. [PMID: 32507886 DOI: 10.1115/1.4047426] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2020] [Indexed: 07/25/2024]
Abstract
Genesis and onset of atherosclerosis are greatly influenced by hemodynamic forces. Two-phase transient computational fluid dynamic (CFD) simulations are performed using a mixture theory model for blood, and a transport equation for low-density lipoprotein (LDL), in idealized and patient-derived abdominal aorta to predict the sites at risk for atherosclerosis. Flow patterns at different time instants and relevant hemodynamic indicators-wall shear stress (WSS)-based (time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), and relative residence time (RRT)), and LDL concentration-are used concurrently to predict the susceptible sites of atherosclerosis. In the case of idealized geometry, flow recirculations are observed on the posterior wall opposite the superior mesenteric artery and below the renal bifurcations. Low TAWSS, high OSI, high RRT and high concentration of LDL are observed in these regions. This suggests that in idealized abdominal aorta, the posterior wall proximal to the renal artery junction is more prone to atherosclerosis. This matches qualitatively with the experimental and simulation data in the literature. In the case of patient-derived geometry, flow reversal, low TAWSS, high OSI and high RRT are observed infrarenal on the anterior wall. Further, high concentration of LDL is observed at the same location on the anterior wall suggesting anterior wall distal to the renal artery junction is more prone to atherosclerosis. These findings demonstrate the use of a novel method to predict the sites at risk for atherosclerosis in geometries where complexities like junctions and curvature play a major role.
Collapse
Affiliation(s)
- Mohammed Ameenuddin
- Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy, Telangana 502285, India
| | - Mohan Anand
- Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Sangareddy, Telangana 502285, India
| |
Collapse
|
23
|
Silva IM, Schelin AB, Viana RL, Caldas IL. Transport of blood particles: Chaotic advection even in a healthy scenario. CHAOS (WOODBURY, N.Y.) 2020; 30:093135. [PMID: 33003941 DOI: 10.1063/5.0013460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2020] [Accepted: 08/31/2020] [Indexed: 06/11/2023]
Abstract
We study the advection of blood particles in the carotid bifurcation, a site that is prone to plaque development. Previously, it has been shown that chaotic advection can take place in blood flows with diseases. Here, we show that even in a healthy scenario, chaotic advection can take place. To understand how the particle dynamics is affected by the emergence and growth of a plaque, we study the carotid bifurcation in three cases: a healthy bifurcation, a bifurcation with a mild stenosis, and the another with a severe stenosis. The result is non-intuitive: there is less chaos for the mild stenosis case even when compared to the healthy, non-stenosed, bifurcation. This happens because the partial obstruction of the mild stenosis generates a symmetry in the flow that does not exist for the healthy condition. For the severe stenosis, there is more irregular motion and more particle trapping as expected.
Collapse
Affiliation(s)
- I M Silva
- Instituto de Física, Universidade de Brasília, Campus Universitário Darcy Ribeiro, 70910-900 Brasília, Distrito Federal, Brazil
| | - A B Schelin
- Instituto de Física, Universidade de Brasília, Campus Universitário Darcy Ribeiro, 70910-900 Brasília, Distrito Federal, Brazil
| | - R L Viana
- Departamento de Física, Universidade Federal do Paraná, 81531-990 Curitiba, Paraná, Brazil
| | - I L Caldas
- Instituto de Física, Universidade de São Paulo, Rua do Matão, 1371, 05508-090 São Paulo, SP, Brazil
| |
Collapse
|
24
|
Mutlu O, Olcay AB, Bilgin C, Hakyemez B. Understanding the effect of effective metal surface area of flow diverter stent's on the patient-specific intracranial aneurysm numerical model using Lagrangian coherent structures. J Clin Neurosci 2020; 80:298-309. [PMID: 32712121 DOI: 10.1016/j.jocn.2020.04.111] [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] [Received: 03/06/2020] [Accepted: 04/19/2020] [Indexed: 11/27/2022]
Abstract
The effective metal surface area (EMSA) of flow diversions plays an essential role in the occlusion mechanism inside the aneurysm since the value of EMSA determines the amount of blood flow into the aneurysm sac. In the present study, three different models of a flow diverter stent, namely FRED 4017, FRED 4038, and FRED 4539, were virtually placed at the aneurysm neck of a 52-years-old female patient to identify the effect of EMSA on stagnation region formation inside the aneurysm sac. Lagrangian coherent structures (LCSs), hyperbolic time, and particle tracking analysis were employed to the velocity vectors obtained from computational fluid dynamics (CFD). It is noticed that use of FRED 4017 stent with 0.42 EMSA value caused nearly 40% of the weightless blood flow particles (more than FRED 4038 and FRED 4539) to stay inside the aneurysm while only 0.35% of the blood flow was remaining inside the aneurysm sac when no stent was placed into the aneurysm site. Furthermore, hyperbolic time computations illustrated the formation of stagnation fluid flow zones that can be associated with the residence time of the blood flow particles. Lastly, the results of hyperbolic time analysis are in good agreement with digital subtraction angiography (DSA) images taken in the clinic a few minutes after a FRED 4017 implantation.
Collapse
Affiliation(s)
- Onur Mutlu
- Yeditepe University, Faculty of Engineering, Department of Mechanical Engineering, Kayisdagi Cad., 34755 Istanbul, Turkey
| | - Ali Bahadır Olcay
- Yeditepe University, Faculty of Engineering, Department of Mechanical Engineering, Kayisdagi Cad., 34755 Istanbul, Turkey.
| | - Cem Bilgin
- Uludag University School of Medicine, Department of Radiology, Gorukle, Bursa 16059, Turkey
| | - Bahattin Hakyemez
- Uludag University School of Medicine, Department of Radiology, Gorukle, Bursa 16059, Turkey
| |
Collapse
|
25
|
Pushin DM, Salikhova TY, Zlobina KE, Guria GT. Platelet activation via dynamic conformational changes of von Willebrand factor under shear. PLoS One 2020; 15:e0234501. [PMID: 32525962 PMCID: PMC7289367 DOI: 10.1371/journal.pone.0234501] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2020] [Accepted: 05/26/2020] [Indexed: 01/10/2023] Open
Abstract
Shear-induced conformational changes of von Willebrand factor (VWF) play an important role in platelet activation. A novel approach describing VWF unfolding on the platelet surface under dynamic shear stress is proposed. Cumulative effect of dynamic shear on platelet activation via conformational changes of VWF is analysed. The critical condition of shear-induced platelet activation is formulated. The explicit expression for the threshold value of cumulative shear stress as a function of VWF multimer size is derived. The results open novel prospects for pharmacological regulation of shear-induced platelet activation through control of VWF multimers size distribution.
Collapse
Affiliation(s)
- Denis M. Pushin
- National Research Center for Hematology, Moscow, Russia
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | | | | | - Georgy Th. Guria
- National Research Center for Hematology, Moscow, Russia
- Moscow Institute of Physics and Technology, Dolgoprudny, Russia
- * E-mail:
| |
Collapse
|
26
|
Jayendiran R, Condemi F, Campisi S, Viallon M, Croisille P, Avril S. Computational prediction of hemodynamical and biomechanical alterations induced by aneurysm dilatation in patient-specific ascending thoracic aortas. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2020; 36:e3326. [PMID: 32087044 DOI: 10.1002/cnm.3326] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/27/2019] [Accepted: 02/15/2020] [Indexed: 06/10/2023]
Abstract
The aim of the present work is to propose a robust computational framework combining computational fluid dynamics (CFD) and 4D flow MRI to predict the progressive changes in hemodynamics and wall rupture index (RPI) induced by aortic morphological evolutions in patients harboring ascending thoracic aortic aneurysms (ATAAs). An analytical equation has been proposed to predict the aneurysm progression based on age, sex, and body surface area. Parameters such as helicity, wall shear stress (WSS), time-averaged WSS, oscillatory shear index, relative residence time, and viscosity were evaluated for two patients at different stages of aneurysm growth, and compared with age-sex-matched healthy subjects. The study shows that evolution of hemodynamics and RPI, despite being very slow in ATAAs, is strongly affected by morphological alterations and, in turn could impact biomechanical factors and aortic mechanobiology. An aspect of the current work is that the patient-specific 4D MRI data sets were obtained with a follow-up of 1 year and the measured time-averaged velocity maps and flow eccentricity were compared with the CFD simulation for validation. The computational framework presented here is capable of capturing the blood flow patterns and the hemodynamic descriptors during the various stages of aneurysm growth. Further investigations will be conducted in order to verify these results on a larger cohort of patients and with long follow-up times to finally elucidate the link between deranged hemodynamics, AA geometry, and wall mechanical properties in ATAAs.
Collapse
Affiliation(s)
- Raja Jayendiran
- Mines Saint-Etienne, Université de Lyon, INSERM, U1059, SAINBIOSE, Saint-Etienne F-42023, France
| | | | - Salvatore Campisi
- Department of Cardiovascular Surgery, University Hospital of Saint-Etienne, Saint-Etienne, France
| | - Magalie Viallon
- UJM-Saint-Etienne, INSA, CNRS UMR 5520, INSERM U1206, CREATIS, Université de Lyon, Saint-Etienne, France
- Department of Radiology, University Hospital of Saint-Etienne, Saint-Etienne, France
| | - Pierre Croisille
- UJM-Saint-Etienne, INSA, CNRS UMR 5520, INSERM U1206, CREATIS, Université de Lyon, Saint-Etienne, France
- Department of Radiology, University Hospital of Saint-Etienne, Saint-Etienne, France
| | - Stéphane Avril
- Mines Saint-Etienne, Université de Lyon, INSERM, U1059, SAINBIOSE, Saint-Etienne F-42023, France
| |
Collapse
|
27
|
Amili O, MacIver R, Coletti F. Magnetic Resonance Imaging Based Flow Field and Lagrangian Particle Tracking From a Left Ventricular Assist Device. J Biomech Eng 2020; 142:021007. [PMID: 31150534 DOI: 10.1115/1.4043939] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2019] [Indexed: 11/08/2022]
Abstract
This study explores the optimal left ventricular assist device (LVAD) cannula outflow configuration in a patient-specific replica of the aorta. The volumetric velocity field is measured using phase-contrast magnetic resonance imaging (PC-MRI) under a physiologically relevant steady flow. The effect of the LVAD outflow graft insertion site and anastomosis angle on the transport of embolic particles to cranial vessels is studied by solving the particle equation of motion for spheres in the range of 0.1-1.0 mm using the measured three-dimensional (3D) velocity field. Results show that for a given aorta anatomy, it is possible to design the cannula graft location and terminal curvature so that the probability of embolic transport to the cranial vessels is significantly minimized. This is particularly important since the complex flow pattern in each cannula case affects the embolic trajectories differently, and hence the common assumption that particles distribute by the volumetric flow division does not hold.
Collapse
Affiliation(s)
- Omid Amili
- Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, MN 55455
| | - Robroy MacIver
- Children's Heart Clinic, Children's Hospitals and Clinics of Minnesota, Minneapolis, MN 55404
| | - Filippo Coletti
- St. Anthony Falls Laboratory, Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, MN 55455
| |
Collapse
|
28
|
Varshney M, Haani Farooqi M, Usmani AY. Quantifying hemodynamics within an aneurysm exposed to prolonged exercise levels. COMPUTER METHODS AND PROGRAMS IN BIOMEDICINE 2020; 184:105124. [PMID: 31627149 DOI: 10.1016/j.cmpb.2019.105124] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2019] [Revised: 10/01/2019] [Accepted: 10/04/2019] [Indexed: 06/10/2023]
Abstract
BACKGROUND AND OBJECTIVE Non-invasive treatment of unruptured Abdominal Aortic Aneurysm involves subjecting the patients to certain physiological levels of the heart. Flow topology (Repeak = 200-1200, frequency: f = 1.18-2.41 Hz) within an aneurysm geometry (2-D) under resting and exercise (mild and moderate) conditions are explored in the present study. Blood is assumed to be Newtonian in nature. Spatio-temporal evolution of the flow patterns and vorticity are established. Hemodynamic indicators (TAWSS and OSI), movement of vortex cores and Particle Residence Index (PRI) are quantified to select an optimum exercise level in attenuating the disease. METHODS The finite volume method is employed for numerical solutions using ANSYS-FluentⓇ software. The SIMPLE scheme has been used for the pressure-velocity coupling. Least Square cell-based method is used for the spatial discretization of the gradients. Second order upwind scheme is considered for discretization of the pressure term. Third order upwind (QUICK) scheme is used to discretize the momentum equation. First order Implicit Scheme was used for the discretization of the temporal terms. Discrete Phase Material (DPM) technique is employed throughout, to visualize the signature of particle deposits within the aneurysm. RESULTS Vortex impingement induces a pressure peak within the aneurysm (moderate) while the peaks are anchored at the proximal and distal ends under resting and mild conditions. Along the averaged flow separation zone, exercise increases the maximum TAWSS from 1.21 N/m2 (mild) to 9.3 N/m2 (moderate). The distal site is exposed to oscillatory loading (OSI = 0.5) under mild activity whereas the loading becomes distributed almost over the entire wall, when subjected to moderate conditions. This in turn, reduces the time involved in 50 percent clearance of particles (PRI = 0.5) from 10.56 s (resting) to 3.98 s (mild) and 0.87 s (moderate), respectively. CONCLUSIONS Resting conditions manifests the aneurysmal wall to recirculating fluid for most of cycle time. Moderate exercise exhibits the least particle clearance time, but it exposes the aneurysmal wall and the distal end to high pressure, which otherwise has low intensity under mild activity. This in turn establishes that mild exercise for prolonged duration can be an optimum level for non-invasive aneurysmal treatment.
Collapse
Affiliation(s)
- Mehul Varshney
- Department of Mechanical Engineering, ZHCET, AMU, Aligarh 202002 India
| | - M Haani Farooqi
- Fluid Mechanics and Energetics Department, École Centrale de Nantes, Nantes 44300 France
| | - Abdullah Y Usmani
- Department of Mechanical Engineering, ZHCET, AMU, Aligarh 202002 India.
| |
Collapse
|
29
|
Evaluating the Effect of the Number of Wire of Flow Diverter Stents on the Nonstagnated Region Formation in an Aneurysm Sac Using Lagrangian Coherent Structure and Hyperbolic Time Analysis. World Neurosurg 2020; 133:e666-e682. [DOI: 10.1016/j.wneu.2019.09.116] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2019] [Revised: 09/21/2019] [Accepted: 09/23/2019] [Indexed: 01/18/2023]
|
30
|
Mazzi V, Gallo D, Calò K, Najafi M, Khan MO, De Nisco G, Steinman DA, Morbiducci U. A Eulerian method to analyze wall shear stress fixed points and manifolds in cardiovascular flows. Biomech Model Mechanobiol 2019; 19:1403-1423. [PMID: 31865482 DOI: 10.1007/s10237-019-01278-3] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2019] [Accepted: 12/08/2019] [Indexed: 12/12/2022]
Abstract
Based upon dynamical systems theory, a fixed point of a vector field such as the wall shear stress (WSS) at the luminal surface of a vessel is a point where the vector field vanishes. Unstable/stable manifolds identify contraction/expansion regions linking fixed points. The significance of such WSS topological features lies in their strong link with "disturbed" flow features like flow stagnation, separation and reversal, deemed responsible for vascular dysfunction initiation and progression. Here, we present a Eulerian method to analyze WSS topological skeleton through the identification and classification of WSS fixed points and manifolds in complex vascular geometries. The method rests on the volume contraction theory and analyzes the WSS topological skeleton through the WSS vector field divergence and Poincar[Formula: see text] index. The method is here applied to computational hemodynamics models of carotid bifurcation and intracranial aneurysm. An in-depth analysis of the time dependence of the WSS topological skeleton along the cardiac cycle is provided, enriching the information obtained from cycle-average WSS. Among the main findings, it emerges that on the carotid bifurcation, instantaneous WSS fixed points co-localize with cycle-average WSS fixed points for a fraction of the cardiac cycle ranging from 0 to [Formula: see text]; a persistent instantaneous WSS fixed point confined on the aneurysm dome does not co-localize with the cycle-average low-WSS region. In conclusion, the here presented approach shows the potential to speed up studies on the physiological significance of WSS topological skeleton in cardiovascular flows, ultimately increasing the chance of finding mechanistic explanations to clinical observations.
Collapse
Affiliation(s)
- Valentina Mazzi
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
- PolitoBIOMed Lab, Politecnico di Torino, Turin, Italy
| | - Diego Gallo
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
- PolitoBIOMed Lab, Politecnico di Torino, Turin, Italy
| | - Karol Calò
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
- PolitoBIOMed Lab, Politecnico di Torino, Turin, Italy
| | - Mehdi Najafi
- Biomedical Simulation Laboratory, Department of Mechanical and Industrial Engineering, University of Toronto, Toronto, ON, Canada
| | | | - Giuseppe De Nisco
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy
- PolitoBIOMed Lab, Politecnico di Torino, Turin, Italy
| | - David A Steinman
- Biomedical Simulation Laboratory, Department of Mechanical and Industrial Engineering, University of Toronto, Toronto, ON, Canada
| | - Umberto Morbiducci
- Department of Mechanical and Aerospace Engineering, Politecnico di Torino, Corso Duca degli Abruzzi 24, 10129, Turin, Italy.
- PolitoBIOMed Lab, Politecnico di Torino, Turin, Italy.
| |
Collapse
|
31
|
Farghadan A, Coletti F, Arzani A. Topological analysis of particle transport in lung airways: Predicting particle source and destination. Comput Biol Med 2019; 115:103497. [DOI: 10.1016/j.compbiomed.2019.103497] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2019] [Revised: 10/09/2019] [Accepted: 10/10/2019] [Indexed: 11/26/2022]
|
32
|
Malone F, McCarthy E, Delassus P, Buhk JH, Fiehler J, Morris L. Embolus Analog Trajectory Paths Under Physiological Flowrates Through Patient-Specific Aortic Arch Models. J Biomech Eng 2019; 141:2734765. [DOI: 10.1115/1.4043832] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2018] [Indexed: 01/10/2023]
Abstract
Atrial fibrillation (AF) is the most common irregular heartbeat among the world's population and is a major contributor to cardiogenic embolisms and acute ischemic stroke (AIS). However, the role AF flow plays in the trajectory paths of cardiogenic emboli has not been experimentally investigated. A physiological simulation system was designed to analyze the trajectory patterns of bovine embolus analogs (EAs) (n = 720) through four patient-specific models, under three flow conditions: steady flow, normal pulsatile flow, and AF pulsatile flow. It was seen that EA trajectory paths were proportional to the percentage flowrate split of 25–31% along the branching vessels. Overall, AF flow conditions increased trajectories through the left- (LCCA) and right (RCCA)-common carotid artery by 25% with respect to normal pulsatile flow. There was no statistical difference in the distribution of clot trajectories when the clot was released from the right, left, or anterior positions. Significantly, more EAs traveled through the brachiocephalic trunk (BCT) than through the LCCA or the left subclavian. Yet of the EAs that traveled through the common carotid arteries, there was a greater affiliation toward the LCCA compared to the RCCA (p < 0.05).
Collapse
Affiliation(s)
- F. Malone
- GMedTech, Department of Mechanical and Industrial Engineering, Galway-Mayo Institute of Technology, Galway H91 T8NW, Ireland e-mail:
| | - E. McCarthy
- GMedTech, Department of Mechanical and Industrial Engineering, Galway-Mayo Institute of Technology, Galway H91 T8NW, Ireland
| | - P. Delassus
- GMedTech, Department of Mechanical and Industrial Engineering, Galway-Mayo Institute of Technology, Galway H91 T8NW, Ireland
| | - J. H. Buhk
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Martinistr. 52, Hamburg 20246, Germany
| | - J. Fiehler
- Department of Diagnostic and Interventional Neuroradiology, University Medical Center Hamburg-Eppendorf, Martinistr. 52, Hamburg 20246, Germany
| | - L. Morris
- GMedTech, Department of Mechanical and Industrial Engineering, Galway-Mayo Institute of Technology, Galway H91 T8NW, Ireland e-mail:
| |
Collapse
|
33
|
Evaluating the Effectiveness of 2 Different Flow Diverter Stents Based on the Stagnation Region Formation in an Aneurysm Sac Using Lagrangian Coherent Structure. World Neurosurg 2019; 127:e727-e737. [DOI: 10.1016/j.wneu.2019.03.255] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2018] [Revised: 03/23/2019] [Accepted: 03/25/2019] [Indexed: 12/16/2022]
|
34
|
In Vitro Study of Particle Transport in Successively Bifurcating Vessels. Ann Biomed Eng 2019; 47:2271-2283. [PMID: 31165293 DOI: 10.1007/s10439-019-02293-2] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2019] [Accepted: 05/20/2019] [Indexed: 12/24/2022]
Abstract
To reach a predictive understanding of how particles travel through bifurcating vessels is of paramount importance in many biomedical settings, including embolization, thromboembolism, and drug delivery. Here we utilize an in vitro model in which solid particles are injected through a rigid vessel that symmetrically bifurcates in successive branching generations. The geometric proportion and fluid dynamics parameters are relevant to the liver embolization. The volumetric flow field is reconstructed via phase-contrast magnetic resonance imaging, from which the particle trajectories are calculated for a range of size and density using the particle equation of motion. The method is validated by directly tracking the injected particles via optical imaging. The results indicate that, opposite to the common assumption, the particles distribution is fundamentally different from the volumetric flow partition. In fact, the amount of delivered particles vary substantially between adjacent branches even when the flow is uniformly distributed. This is not due to the inertia of the particles, nor to gravity. The particle distribution is rather rooted in their different pathways, which in turn are linked to their release origin along the main vessel cross-section. Therefore, the tree geometry and the associated flow streamlines are the prime determinant of the particle fate, while local changes of volumetric flow rate to selected branches do not generally produce proportional changes of particle delivery.
Collapse
|
35
|
Salman HE, Ramazanli B, Yavuz MM, Yalcin HC. Biomechanical Investigation of Disturbed Hemodynamics-Induced Tissue Degeneration in Abdominal Aortic Aneurysms Using Computational and Experimental Techniques. Front Bioeng Biotechnol 2019; 7:111. [PMID: 31214581 PMCID: PMC6555197 DOI: 10.3389/fbioe.2019.00111] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2019] [Accepted: 05/02/2019] [Indexed: 11/13/2022] Open
Abstract
Abdominal aortic aneurysm (AAA) is the dilatation of the aorta beyond 50% of the normal vessel diameter. It is reported that 4-8% of men and 0.5-1% of women above 50 years of age bear an AAA and it accounts for ~15,000 deaths per year in the United States alone. If left untreated, AAA might gradually expand until rupture; the most catastrophic complication of the aneurysmal disease that is accompanied by a striking overall mortality of 80%. The precise mechanisms leading to AAA rupture remains unclear. Therefore, characterization of disturbed hemodynamics within AAAs will help to understand the mechanobiological development of the condition which will contribute to novel therapies for the condition. Due to geometrical complexities, it is challenging to directly quantify disturbed flows for AAAs clinically. Two other approaches for this investigation are computational modeling and experimental flow measurement. In computational modeling, the problem is first defined mathematically, and the solution is approximated with numerical techniques to get characteristics of flow. In experimental flow measurement, once the setup providing physiological flow pattern in a phantom geometry is constructed, velocity measurement system such as particle image velocimetry (PIV) enables characterization of the flow. We witness increasing number of applications of these complimentary approaches for AAA investigations in recent years. In this paper, we outline the details of computational modeling procedures and experimental settings and summarize important findings from recent studies, which will help researchers for AAA investigations and rupture mechanics.
Collapse
Affiliation(s)
| | - Burcu Ramazanli
- Department of Mechanical Engineering, Middle East Technical University, Ankara, Turkey
| | - Mehmet Metin Yavuz
- Department of Mechanical Engineering, Middle East Technical University, Ankara, Turkey
| | | |
Collapse
|
36
|
Sunderland K, Huang Q, Strother C, Jiang J. Two closely-spaced Aneurysms of the Supraclinoid Internal Carotid Artery: How Does One Influence the Other? J Biomech Eng 2019; 141:2735303. [PMID: 31141586 DOI: 10.1115/1.4043868] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2017] [Indexed: 11/08/2022]
Abstract
The objective of this study was to use image-based CFD simulation techniques to analyze the impact that multiple closely spaced IAs of the supra-clinioid segment of the ICA have on each other's hemodynamic characteristics. The vascular geometry of fifteen (15) subjects with 2 IAs were gathered using a 3D clinical system. Two groups of computer models were created for each subject's vascular geometry: both IAs present (Model A) and after removal of one IA (Model B). Models were separated into two groups based on IA separation: tandem (one proximal and one distal) and tandem (aneurysms directly opposite on a vessel). Simulations using a pulsatile velocity waveform were solved by a commercial CFD solver. Proximal IAs altered flow into distal IAs (5 of 7), increasing flow energy and spatial-temporally averaged wall shear stress (STA-WSS: 3-50\% comparing Model A to B) while decreasing flow stability within distal IAs. Thus, proximal IAs may ``protect" a distal aneurysm from destructive remodeling due to flow stagnation. Among adjacent IAs, the presence of both IAs decreased each other's flow characteristics, lowering WSS (Model A to B) and increasing flow stability: all changes statistically significant (t-test p < 0.05). A negative relationship exists between the mean percent change in flow stability in relation to adjacent IA volume and ostium area. Closely spaced IAs impact hemodynamic alterations onto each other concerning flow energy, stressors and stability. Understanding these alterations may improve clinical management of closely-spaced IAs.
Collapse
Affiliation(s)
- Kevin Sunderland
- Department of Biomedical Engineering, Michigan Technological University, Houghton, Michigan 49931
| | - Qinghai Huang
- Department of Neurosurgery, Chonghai Hospital, Second Military University, Shanghai, China
| | - Charlie Strother
- Department of Radiology, University of Wisconsin School of Medicine and Public Health, Madison, Wisconsin, 53705
| | - Jingfeng Jiang
- Department of Biomedical Engineering, Michigan Technological University, Houghton, Michigan 49931
| |
Collapse
|
37
|
Reza MMS, Arzani A. A critical comparison of different residence time measures in aneurysms. J Biomech 2019; 88:122-129. [DOI: 10.1016/j.jbiomech.2019.03.028] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2018] [Revised: 03/16/2019] [Accepted: 03/20/2019] [Indexed: 12/13/2022]
|
38
|
Leemans EL, Cornelissen BMW, Rosalini G, Verbaan D, Schneiders JJ, van den Berg R, Vandertop WP, van Bavel ET, Slump CH, Majoie CBLM, Marquering HA. Impact of Intracranial Aneurysm Morphology and Rupture Status on the Particle Residence Time. J Neuroimaging 2019; 29:487-492. [PMID: 31002750 PMCID: PMC6618041 DOI: 10.1111/jon.12618] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2019] [Revised: 03/29/2019] [Accepted: 04/02/2019] [Indexed: 11/30/2022] Open
Abstract
BACKGROUND AND PURPOSE Aneurysm hemodynamics play an important role in aneurysm growth and subsequent rupture. Within the available hemodynamic characteristics, particle residence time (PRT) is relatively unexplored. However, some studies have shown that PRT is related to thrombus formation and inflammation. The goal of this study is to evaluate the association between PRT and aneurysm rupture and morphology. METHODS We determined the PRT for 113 aneurysms (61 unruptured, 53 ruptured) based on computational fluid dynamic models. Virtual particles were injected into the parent vessel and followed during multiple cardiac cycles. PRT was defined as the time needed for 99% of the particles that entered an aneurysm to leave the aneurysm. Subsequently, we evaluated the association between PRT, rupture, and morphology (aneurysm type, presence of blebs, or multiple lobulations). RESULTS PRT showed no significant difference between unruptured (1.1 seconds interquartile range [IQR .39-2.0 seconds]) and ruptured aneurysms (1.2 seconds [IQR .47-2.3 seconds]). PRT was influenced by aneurysm morphology. Longer PRTs were seen in bifurcation aneurysms (1.3 seconds [IQR .54-2.4 seconds], P = .01) and aneurysms with blebs or multiple lobulations (1.92 seconds [IQR .94-2.8 seconds], P < .001). Four of five partially thrombosed aneurysms had a long residence time (>1.9 seconds). CONCLUSIONS Our study shows an influence of aneurysm morphology on PRT. Nevertheless, it suggests that PRT cannot be used to differentiate unruptured and ruptured aneurysms.
Collapse
Affiliation(s)
- E L Leemans
- Department of Biomedical Engineering and Physics, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands.,Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| | - B M W Cornelissen
- Department of Biomedical Engineering and Physics, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands.,Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands.,MIRA Institute for Biomedical Engineering and Technical Medicine, University of Twente, Enschede, the Netherlands
| | - G Rosalini
- Department of Biomedical Engineering and Physics, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands.,Department of Industrial Engineering and Information, University of Pavia, Pavia, Italy
| | - D Verbaan
- Neurosurgical Center Amsterdam, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| | - J J Schneiders
- Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| | - R van den Berg
- Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| | - W P Vandertop
- Neurosurgical Center Amsterdam, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| | - E T van Bavel
- Department of Biomedical Engineering and Physics, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| | - C H Slump
- MIRA Institute for Biomedical Engineering and Technical Medicine, University of Twente, Enschede, the Netherlands
| | - C B L M Majoie
- Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| | - H A Marquering
- Department of Biomedical Engineering and Physics, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands.,Department of Radiology and Nuclear Medicine, Amsterdam UMC, University of Amsterdam, Amsterdam, the Netherlands
| |
Collapse
|
39
|
Conti M, Vandenberghe S, Marconi S, Ferrari E, Romarowski RM, Morganti S, Auricchio F, Demertzis S. Reversed Auxiliary Flow to Reduce Embolism Risk During TAVI: A Computational Simulation and Experimental Study. Cardiovasc Eng Technol 2018; 10:124-135. [PMID: 30341729 DOI: 10.1007/s13239-018-00386-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/23/2018] [Accepted: 10/11/2018] [Indexed: 12/21/2022]
Abstract
INTRODUCTION Endovascular treatments, such as transcatheter aortic valve implantation (TAVI), carry a risk of embolization due to debris dislodgement during various procedural steps. Although embolic filters are already available and marketed, mechanisms underlying cerebral embolism still need to be elucidated in order to further reduce cerebrovascular events. METHODS We propose an experimental framework with an in silico duplicate allowing release of particles at the level of the aortic valve and their subsequent capture in the supra-aortic branches, simulating embolization under constant inflow and controlled hemodynamic conditions. The effect of a simple flow modulation, consisting of an auxiliary constant flow via the right subclavian artery (RSA), on the amount of particle entering the brachiocephalic trunk was investigated. Preliminary computational fluid dynamics (CFD) simulations were performed in order to assess the minimum retrograde flow-rate from RSA required to deviate particles. RESULTS Our results show that a constant reversed auxiliary flow of 0.5 L/min from the RSA under a constant inflow of 4 L/min from the ascending aorta is able to protect the brachiocephalic trunk from particle embolisms. Both computational and experimental results also demonstrate that the distribution of the bulk flow dictates the distribution of the particles along the aortic branches. This effect has also shown to be independent of release location and flow rate. CONCLUSIONS The present study confirms that the integration of in vitro experiments and in silico analyses allows designing and benchmarking novel solutions for cerebral embolic protection during TAVI such as the proposed embo-deviation technique based on an auxiliary retrograde flow from the right subclavian artery.
Collapse
Affiliation(s)
- Michele Conti
- Department of Civil Engineering and Architecture, University of Pavia, Via Ferrata 3, 27100, Pavia, Italy.
| | | | - Stefania Marconi
- Department of Civil Engineering and Architecture, University of Pavia, Via Ferrata 3, 27100, Pavia, Italy
| | - Enrico Ferrari
- Department of Cardiac Surgery, Cardiocentro Ticino, Lugano, Switzerland
| | - Rodrigo M Romarowski
- 3D and Computer Simulation Laboratory, IRCCS Policlinico San Donato, San Donato Milanese, Italy
| | - Simone Morganti
- Department of Electrical, Computer, and Biomedical Engineering, University of Pavia, Pavia, Italy
| | - Ferdinando Auricchio
- Department of Civil Engineering and Architecture, University of Pavia, Via Ferrata 3, 27100, Pavia, Italy
| | | |
Collapse
|
40
|
Arzani A. Accounting for residence-time in blood rheology models: do we really need non-Newtonian blood flow modelling in large arteries? J R Soc Interface 2018; 15:rsif.2018.0486. [PMID: 30257924 DOI: 10.1098/rsif.2018.0486] [Citation(s) in RCA: 63] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Accepted: 09/03/2018] [Indexed: 12/27/2022] Open
Abstract
Patient-specific computational fluid dynamics (CFD) is a promising tool that provides highly resolved haemodynamics information. The choice of blood rheology is an assumption in CFD models that has been subject to extensive debate. Blood is known to exhibit shear-thinning behaviour, and non-Newtonian modelling has been recommended for aneurysmal flows. Current non-Newtonian models ignore rouleaux formation, which is the key player in blood's shear-thinning behaviour. Experimental data suggest that red blood cell aggregation and rouleaux formation require notable red blood cell residence-time (RT) in a low shear rate regime. This study proposes a novel hybrid Newtonian and non-Newtonian rheology model where the shear-thinning behaviour is activated in high RT regions based on experimental data. Image-based abdominal aortic and cerebral aneurysm models are considered and highly resolved CFD simulations are performed using a minimally dissipative solver. Lagrangian particle tracking is used to define a backward particle RT measure and detect stagnant regions with increased rouleaux formation likelihood. Our novel RT-based non-Newtonian model shows a significant reduction in shear-thinning effects and provides haemodynamic results qualitatively identical and quantitatively close to the Newtonian model. Our results have important implications in patient-specific CFD modelling and suggest that non-Newtonian models should be revisited in large artery flows.
Collapse
Affiliation(s)
- Amirhossein Arzani
- Department of Mechanical Engineering, Northern Arizona University, Flagstaff, AZ, USA
| |
Collapse
|
41
|
Ghaffari M, Alaraj A, Du X, Zhou XJ, Charbel FT, Linninger AA. Quantification of near-wall hemodynamic risk factors in large-scale cerebral arterial trees. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2018; 34:e2987. [PMID: 29601146 PMCID: PMC6043404 DOI: 10.1002/cnm.2987] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/12/2017] [Revised: 03/21/2018] [Accepted: 03/23/2018] [Indexed: 05/18/2023]
Abstract
Detailed hemodynamic analysis of blood flow in pathological segments close to aneurysm and stenosis has provided physicians with invaluable information about the local flow patterns leading to vascular disease. However, these diseases have both local and global effects on the circulation of the blood within the cerebral tree. The aim of this paper is to demonstrate the importance of extending subject-specific hemodynamic simulations to the entire cerebral arterial tree with hundreds of bifurcations and vessels, as well as evaluate hemodynamic risk factors and waveform shape characteristics throughout the cerebral arterial trees. Angioarchitecture and in vivo blood flow measurement were acquired from healthy subjects and in cases with symptomatic intracranial aneurysm and stenosis. A global map of cerebral arterial blood flow distribution revealed regions of low to high hemodynamic risk that may significantly contribute to the development of intracranial aneurysms or atherosclerosis. Comparison of pre-intervention and post-intervention of pathological cases further shows large angular phase shift (~33.8°), and an augmentation of the peak-diastolic velocity. Hemodynamic indexes of waveform analysis revealed on average a 16.35% reduction in the pulsatility index after treatment from lesion site to downstream distal vessels. The lesion regions not only affect blood flow streamlines of the proximal sites but also generate pulse wave shift and disturbed flow in downstream vessels. This network effect necessitates the use of large-scale simulation to visualize both local and global effects of pathological lesions.
Collapse
Affiliation(s)
- Mahsa Ghaffari
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
| | - Ali Alaraj
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
- Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
| | - Xinjian Du
- Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
| | - Xiaohong Joe Zhou
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
- Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
- Department of Radiology, University of Illinois at Chicago, Chicago, IL, USA
- Center for MR Research, University of Illinois at Chicago, Chicago, IL, USA
| | - Fady T. Charbel
- Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
| | - Andreas A. Linninger
- Department of Bioengineering, University of Illinois at Chicago, Chicago, IL, USA
- Department of Neurosurgery, University of Illinois at Chicago, Chicago, IL, USA
| |
Collapse
|
42
|
Arzani A, Shadden SC. Wall shear stress fixed points in cardiovascular fluid mechanics. J Biomech 2018; 73:145-152. [DOI: 10.1016/j.jbiomech.2018.03.034] [Citation(s) in RCA: 38] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2018] [Revised: 03/09/2018] [Accepted: 03/21/2018] [Indexed: 12/13/2022]
|
43
|
Tan J, Sinno T, Diamond SL. A parallel fluid-solid coupling model using LAMMPS and Palabos based on the immersed boundary method. JOURNAL OF COMPUTATIONAL SCIENCE 2018; 25:89-100. [PMID: 30220942 PMCID: PMC6136258 DOI: 10.1016/j.jocs.2018.02.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/09/2023]
Abstract
The study of viscous fluid flow coupled with rigid or deformable solids has many applications in biological and engineering problems, e.g., blood cell transport, drug delivery, and particulate flow. We developed a partitioned approach to solve this coupled Multiphysics problem. The fluid motion was solved by Palabos (Parallel Lattice Boltzmann Solver), while the solid displacement and deformation was simulated by LAMMPS (Large-scale Atomic/Molecular Massively Parallel Simulator). The coupling was achieved through the immersed boundary method (IBM). The code modeled both rigid and deformable solids exposed to flow. The code was validated with the Jeffery orbits of an ellipsoid particle in shear flow, red blood cell stretching test, and effective blood viscosity flowing in tubes. It demonstrated essentially linear scaling from 512 to 8192 cores for both strong and weak scaling cases. The computing time for the coupling increased with the solid fraction. An example of the fluid-solid coupling was given for flexible filaments (drug carriers) transport in a flowing blood cell suspensions, highlighting the advantages and capabilities of the developed code.
Collapse
Affiliation(s)
- Jifu Tan
- Department of Mechanical Engineering, Northern Illinois University, DeKalb, IL 60115, USA
| | - Talid Sinno
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA
19104, USA
| | - Scott L Diamond
- Department of Chemical and Biomolecular Engineering, University of Pennsylvania, Philadelphia, PA
19104, USA
| |
Collapse
|
44
|
Hemodynamics in a giant intracranial aneurysm characterized by in vitro 4D flow MRI. PLoS One 2018; 13:e0188323. [PMID: 29300738 PMCID: PMC5754057 DOI: 10.1371/journal.pone.0188323] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/21/2017] [Accepted: 11/03/2017] [Indexed: 11/19/2022] Open
Abstract
Experimental and computational data suggest that hemodynamics play a critical role in the development, growth, and rupture of cerebral aneurysms. The flow structure, especially in aneurysms with a large sac, is highly complex and three-dimensional. Therefore, volumetric and time-resolved measurements of the flow properties are crucial to fully characterize the hemodynamics. In this study, phase-contrast Magnetic Resonance Imaging is used to assess the fluid dynamics inside a 3D-printed replica of a giant intracranial aneurysm, whose hemodynamics was previously simulated by multiple research groups. The physiological inflow waveform is imposed in a flow circuit with realistic cardiovascular impedance. Measurements are acquired with sub-millimeter spatial resolution for 16 time steps over a cardiac cycle, allowing for the detailed reconstruction of the flow evolution. Moreover, the three-dimensional and time-resolved pressure distribution is calculated from the velocity field by integrating the fluid dynamics equations, and is validated against differential pressure measurements using precision transducers. The flow structure is characterized by vortical motions that persist within the aneurysm sac for most of the cardiac cycle. All the main flow statistics including velocity, vorticity, pressure, and wall shear stress suggest that the flow pattern is dictated by the aneurysm morphology and is largely independent of the pulsatility of the inflow, at least for the flow regimes investigated here. Comparisons are carried out with previous computational simulations that used the same geometry and inflow conditions, both in terms of cycle-averaged and systolic quantities.
Collapse
|
45
|
Flow stagnation volume and abdominal aortic aneurysm growth: Insights from patient-specific computational flow dynamics of Lagrangian-coherent structures. Comput Biol Med 2018; 92:98-109. [DOI: 10.1016/j.compbiomed.2017.10.033] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2017] [Revised: 10/09/2017] [Accepted: 10/28/2017] [Indexed: 12/23/2022]
|
46
|
Kelsey LJ, Powell JT, Norman PE, Miller K, Doyle BJ. A comparison of hemodynamic metrics and intraluminal thrombus burden in a common iliac artery aneurysm. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2017; 33:e2821. [PMID: 27509188 DOI: 10.1002/cnm.2821] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/18/2016] [Revised: 06/06/2016] [Accepted: 07/02/2016] [Indexed: 06/06/2023]
Abstract
Aneurysms of the common iliac artery (CIAA) are typically found in association with an abdominal aortic aneurysm (AAA). Isolated CIAAs, in the absence of an AAA, are uncommon. Similar to AAAs, CIAA may develop intraluminal thrombus (ILT). As isolated CIAAs have a contralateral common iliac artery for comparison, they provide an opportunity to study the hemodynamic mechanisms behind ILT formation. In this study, we compared a large isolated CIAA and the contralateral iliac artery using computational fluid dynamics to determine if hemodynamic metrics correlate with the location of ILT. We performed a comprehensive computational fluid dynamics study and investigated the residence time of platelets and monocytes, velocity fields, time-averaged wall shear stress, oscillatory shear index, and endothelial cell activation potential. We then correlated these data to ILT burden determined with computed tomography. We found that high cell residence times, low time-averaged wall shear stress, high oscillatory shear index, and high endothelial cell activation potential all correlate with regions of ILT development. Our results show agreement with previous hypotheses of thrombus formation in AAA and provide insights into the computational hemodynamics of iliac artery aneurysms.
Collapse
Affiliation(s)
- Lachlan J Kelsey
- Vascular Engineering Laboratory, Harry Perkins Institute of Medical Research, Perth, WA, Australia
- Intelligent Systems for Medicine Laboratory, School of Mechanical and Chemical Engineering, The University of Western Australia, Crawley, WA, Australia
| | - Janet T Powell
- Vascular Surgery Research Group, Imperial College London, London, UK
| | - Paul E Norman
- Vascular Engineering Laboratory, Harry Perkins Institute of Medical Research, Perth, WA, Australia
- School of Surgery, The University of Western Australia, Crawley, WA, Australia
| | - Karol Miller
- Intelligent Systems for Medicine Laboratory, School of Mechanical and Chemical Engineering, The University of Western Australia, Crawley, WA, Australia
- Institute of Mechanics and Advanced Materials, Cardiff University, Cardiff, UK
| | - Barry J Doyle
- Vascular Engineering Laboratory, Harry Perkins Institute of Medical Research, Perth, WA, Australia
- School of Mechanical and Chemical Engineering, The University of Western Australia, Crawley, WA, Australia
- British Heart Foundation Centre for Cardiovascular Science, The University of Edinburgh, Edinburgh, UK
| |
Collapse
|
47
|
Inertial particle dynamics in large artery flows – Implications for modeling arterial embolisms. J Biomech 2017; 52:155-164. [DOI: 10.1016/j.jbiomech.2016.12.028] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2016] [Revised: 11/07/2016] [Accepted: 12/21/2016] [Indexed: 10/20/2022]
|
48
|
García Carrascal P, García García J, Sierra Pallares J, Castro Ruiz F, Manuel Martín FJ. Numerical Study of Blood Clots Influence on the Flow Pattern and Platelet Activation on a Stented Bifurcation Model. Ann Biomed Eng 2016; 45:1279-1291. [DOI: 10.1007/s10439-016-1782-4] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2016] [Accepted: 12/19/2016] [Indexed: 01/09/2023]
|
49
|
Arzani A, Gambaruto AM, Chen G, Shadden SC. Wall shear stress exposure time: a Lagrangian measure of near-wall stagnation and concentration in cardiovascular flows. Biomech Model Mechanobiol 2016; 16:787-803. [DOI: 10.1007/s10237-016-0853-7] [Citation(s) in RCA: 45] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2016] [Accepted: 11/08/2016] [Indexed: 12/18/2022]
|
50
|
Faghih MM, Keith Sharp M. Extending the Power-Law Hemolysis Model to Complex Flows. J Biomech Eng 2016; 138:2556264. [DOI: 10.1115/1.4034786] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2016] [Indexed: 11/08/2022]
Abstract
Hemolysis (damage to red blood cells) is a long-standing problem in blood contacting devices, and its prediction has been the goal of considerable research. The most popular model relating hemolysis to fluid stresses is the power-law model, which was developed from experiments in pure shear only. In the absence of better data, this model has been extended to more complex flows by replacing the shear stress in the power-law equation with a von Mises-like scalar stress. While the validity of the scalar stress also remains to be confirmed, inconsistencies exist in its application, in particular, two forms that vary by a factor of 2 have been used. This article will clarify the proper extension of the power law to complex flows in a way that maintains correct results in the limit of pure shear.
Collapse
Affiliation(s)
- Mohammad M. Faghih
- Biofluid Mechanics Laboratory, Department of Mechanical Engineering, University of Louisville, Louisville, KY 40292
| | - M. Keith Sharp
- Biofluid Mechanics Laboratory, Department of Mechanical Engineering, University of Louisville, Louisville, KY 40292
| |
Collapse
|