{"title":"An automated and time-efficient framework for simulation of coronary blood flow under steady and pulsatile conditions","authors":"","doi":"10.1016/j.cmpb.2024.108415","DOIUrl":null,"url":null,"abstract":"<div><h3>Background and objective</h3><p>Invasive fractional flow reserve (FFR) measurement is the gold standard method for coronary artery disease (CAD) diagnosis. FFR-CT exploits computational fluid dynamics (CFD) for non-invasive evaluation of FFR, simulating coronary flow in virtual geometries reconstructed from computed tomography (CT), but suffers from cost-intensive computing process and uncertainties in the definition of patient specific boundary conditions (BCs). In this work, we investigated the use of time-averaged steady BCs, compared to pulsatile to reduce the computational time and deployed a self-adjusting method for the tuning of BCs to patient-specific clinical data.</p></div><div><h3>Methods</h3><p>133 coronary arteries were reconstructed form CT images of patients suffering from CAD. For each vessel, invasive FFR was measured. After segmentation, the geometries were prepared for CFD simulation by clipping the outlets and discretizing into tetrahedral mesh. Steady BCs were defined in two steps: (<em>i</em>) rest BCs were extrapolated from clinical and image-derived data; (<em>ii</em>) hyperemic BCs were computed from resting conditions. Flow rate was iteratively adjusted during the simulation, until patient's aortic pressure was matched. Pulsatile BCs were defined exploiting the convergence values of steady BCs. After CFD simulation, lesion-specific hemodynamic indexes were computed and compared between group of patients for which surgery was indicated and not. The whole pipeline was implemented as a straightforward process, in which each single step is performed automatically.</p></div><div><h3>Results</h3><p>Steady and pulsatile FFR-CT yielded a strong correlation (<em>r</em> = 0.988, <em>p</em> < 0.001) and correlated with invasive FFR (<em>r</em> = 0.797, <em>p</em> < 0.001). The per-point difference between the pressure and FFR-CT field predicted by the two methods was below 1 % and 2 %, respectively. Both approaches exhibited a good diagnostic performance: accuracy was 0.860 and 0.864, the AUC was 0.923 and 0.912, for steady and pulsatile case, respectively. The computational time required by steady BCs CFD was approximatively 30-folds lower than pulsatile case.</p></div><div><h3>Conclusions</h3><p>This work shows the feasibility of using steady BCs CFD for computing the FFR-CT in coronary arteries, as well as its computational and diagnostic performance within a fully automated pipeline.</p></div>","PeriodicalId":10624,"journal":{"name":"Computer methods and programs in biomedicine","volume":null,"pages":null},"PeriodicalIF":4.9000,"publicationDate":"2024-09-06","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://www.sciencedirect.com/science/article/pii/S0169260724004085/pdfft?md5=041392929b671b0f9c1c0a15b7e05d66&pid=1-s2.0-S0169260724004085-main.pdf","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Computer methods and programs in biomedicine","FirstCategoryId":"5","ListUrlMain":"https://www.sciencedirect.com/science/article/pii/S0169260724004085","RegionNum":2,"RegionCategory":"医学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"COMPUTER SCIENCE, INTERDISCIPLINARY APPLICATIONS","Score":null,"Total":0}
引用次数: 0
Abstract
Background and objective
Invasive fractional flow reserve (FFR) measurement is the gold standard method for coronary artery disease (CAD) diagnosis. FFR-CT exploits computational fluid dynamics (CFD) for non-invasive evaluation of FFR, simulating coronary flow in virtual geometries reconstructed from computed tomography (CT), but suffers from cost-intensive computing process and uncertainties in the definition of patient specific boundary conditions (BCs). In this work, we investigated the use of time-averaged steady BCs, compared to pulsatile to reduce the computational time and deployed a self-adjusting method for the tuning of BCs to patient-specific clinical data.
Methods
133 coronary arteries were reconstructed form CT images of patients suffering from CAD. For each vessel, invasive FFR was measured. After segmentation, the geometries were prepared for CFD simulation by clipping the outlets and discretizing into tetrahedral mesh. Steady BCs were defined in two steps: (i) rest BCs were extrapolated from clinical and image-derived data; (ii) hyperemic BCs were computed from resting conditions. Flow rate was iteratively adjusted during the simulation, until patient's aortic pressure was matched. Pulsatile BCs were defined exploiting the convergence values of steady BCs. After CFD simulation, lesion-specific hemodynamic indexes were computed and compared between group of patients for which surgery was indicated and not. The whole pipeline was implemented as a straightforward process, in which each single step is performed automatically.
Results
Steady and pulsatile FFR-CT yielded a strong correlation (r = 0.988, p < 0.001) and correlated with invasive FFR (r = 0.797, p < 0.001). The per-point difference between the pressure and FFR-CT field predicted by the two methods was below 1 % and 2 %, respectively. Both approaches exhibited a good diagnostic performance: accuracy was 0.860 and 0.864, the AUC was 0.923 and 0.912, for steady and pulsatile case, respectively. The computational time required by steady BCs CFD was approximatively 30-folds lower than pulsatile case.
Conclusions
This work shows the feasibility of using steady BCs CFD for computing the FFR-CT in coronary arteries, as well as its computational and diagnostic performance within a fully automated pipeline.
期刊介绍:
To encourage the development of formal computing methods, and their application in biomedical research and medical practice, by illustration of fundamental principles in biomedical informatics research; to stimulate basic research into application software design; to report the state of research of biomedical information processing projects; to report new computer methodologies applied in biomedical areas; the eventual distribution of demonstrable software to avoid duplication of effort; to provide a forum for discussion and improvement of existing software; to optimize contact between national organizations and regional user groups by promoting an international exchange of information on formal methods, standards and software in biomedicine.
Computer Methods and Programs in Biomedicine covers computing methodology and software systems derived from computing science for implementation in all aspects of biomedical research and medical practice. It is designed to serve: biochemists; biologists; geneticists; immunologists; neuroscientists; pharmacologists; toxicologists; clinicians; epidemiologists; psychiatrists; psychologists; cardiologists; chemists; (radio)physicists; computer scientists; programmers and systems analysts; biomedical, clinical, electrical and other engineers; teachers of medical informatics and users of educational software.