Fluid-Solid Interaction Analysis for Developing In-Situ Strain and Flow Sensors for Prosthetic Valve Monitoring.

Transcatheter aortic valve implantation (TAVI) was initially developed for adult patients, but there is a growing interest to expand this procedure to younger individuals with longer life expectancies. However, the gradual degradation of biological valve leaflets in transcatheter heart valves (THV) presents significant challenges for this extension. This study aimed to establish a multiphysics computational framework to analyze structural and flow measurements of TAVI and evaluate the integration of optical fiber and photoplethysmography (PPG) sensors for monitoring valve function. A two-way fluid-solid interaction (FSI) analysis was performed on an idealized aortic vessel before and after the virtual deployment of the SAPIEN 3 Ultra (S3) THV. Subsequently, an analytical analysis was conducted to estimate the PPG signal using computational flow predictions and to analyze the effect of different pressure gradients and distances between PPG sensors. Circumferential strain estimates from the embedded optical fiber in the FSI model were highest in the sinus of Valsalva; however, the optimal fiber positioning was found to be distal to the sino-tubular junction to minimize bending effects. The findings also demonstrated that positioning PPG sensors both upstream and downstream of the bioprosthesis can be used to effectively assess the pressure gradient across the valve. We concluded that computational modeling allows sensor design to quantify vessel wall strain and pressure gradients across valve leaflets, with the ultimate goal of developing low-cost monitoring systems for detecting valve deterioration.

A Study on the Effect of Contact Pressure during Physical Activity on Photoplethysmographic Heart Rate Measurements.

Heart rate (HR) as an important physiological indicator could properly describe global subject's physical status. Photoplethysmographic (PPG) sensors are catching on in field of wearable sensors, combining the advantages in costs, weight and size. Nevertheless, accuracy in HR readings is unreliable specifically during physical activity. Among several identified sources that affect PPG recording, contact pressure (CP) between the PPG sensor and skin greatly influences the signals. METHODS: In this study, the accuracy of HR measurements of a PPG sensor at different CP was investigated when compared with a commercial ECG-based chest strap used as a test control, with the aim of determining the optimal CP to produce a reliable signal during physical activity. Seventeen subjects were enrolled for the study to perform a physical activity at three different rates repeated at three different contact pressures of the PPG-based wristband. RESULTS: The results show that the CP of 54 mmHg provides the most accurate outcome with a Pearson correlation coefficient ranging from 0.81 to 0.95 and a mean average percentage error ranging from 3.8% to 2.4%, based on the physical activity rate. CONCLUSION: Authors found that changes in the CP have greater effects on PPG-HR signal quality than those deriving from the intensity of the physical activity and specifically, the individual best CP for each subject provided reliable HR measurements even for a high intensity of physical exercise with a Bland-Altman plot within ±11 bpm. Although future studies on a larger cohort of subjects are still needed, this study could contribute a profitable indication to enhance accuracy of PPG-based wearable devices.

Computational modeling of bicuspid aortopathy: Towards personalized risk strategies.

This paper describes current advances on the application of in-silico for the understanding of bicuspid aortopathy and future perspectives of this technology on routine clinical care. This includes the impact that artificial intelligence can provide to develop computer-based clinical decision support system and that wearable sensors can offer to remotely monitor high-risk bicuspid aortic valve (BAV) patients. First, we discussed the benefit of computational modeling by providing tangible examples of in-silico software products based on computational fluid-dynamic (CFD) and finite-element method (FEM) that are currently transforming the way we diagnose and treat cardiovascular diseases. Then, we presented recent findings on computational hemodynamic and structural mechanics of BAV to highlight the potentiality of patient-specific metrics (not-based on aortic size) to support the clinical-decision making process of BAV-associated aneurysms. Examples of BAV-related personalized healthcare solutions are illustrated.

Modeling Right Ventricle Failure After Continuous Flow Left Ventricular Assist Device: A Biventricular Finite-Element and Lumped-Parameter Analysis.

The risk of right ventricle (RV) failure remains a major contraindication for continuous-flow left ventricular assist device (CF-LVAD) implantation in patients with heart failure. It is therefore critical to identify the patients who will benefit from early intervention to avoid adverse outcomes. We sought to advance the computational modeling description of the mechanisms underlying RV failure in LVAD-supported patients. RV failure was studied by computational modeling of hemodynamic and cardiac mechanics using lumped-parameter and biventricular finite element (FE) analysis. Findings were validated by comparison of bi-dimensional speckle-tracking echocardiographic strain assessment of the RV free wall vs. patient-specific computational strain estimations, and by non-invasive lumped-based hemodynamic predictions vs. invasive right heart catheterization data. Correlation analysis revealed that lumped-derived RV cardiac output (R = 0.94) and RV stroke work index (R = 0.85) were in good agreement with catheterization data collected from 7 patients with CF-LVAD. Biventricular FE analysis showed abnormal motion of the interventricular septum towards the left ventricular free wall, suggesting impaired right heart mechanics. Good agreement between computationally predicted and echocardiographic measured longitudinal strains was found at basal (- 19.1 ± 3.0% for ECHO, and - 16.4 ± 3.2% for FEM), apical (- 20.0 ± 3.7% for ECHO, and - 17.4 ± 2.7% for FEM), and mid-level of the RV free wall (- 20.1 ± 5.9% for echo, and - 18.0 ± 5.4% for FEM). Simulation approach here presented could serve as a tool for less invasive and early diagnosis of the severity of RV failure in patients with LVAD, although future studies are needed to validate our findings against clinical outcomes.

Biomechanical Determinants of Right Ventricular Failure in Pulmonary Hypertension.

Pulmonary hypertension (PH) is a disease characterized by progressive adverse remodeling of the distal pulmonary arteries, resulting in elevated pulmonary vascular resistance and load pressure on the right ventricle (RV), ultimately leading to RV failure. Invasive hemodynamic testing is the gold standard for diagnosing PH and guiding patient therapy. We hypothesized that lumped-parameter and biventricular finite-element (FE) modeling may lead to noninvasive predictions of both PH-related hemodynamic and biomechanical parameters that induce PH. We created patient-specific biventricular FE models that characterize the biomechanical response of the heart and coupled them with a lumped-parameter model that represents the systemic and pulmonic circulation. Simulations were calibrated by adjusting the pulmonary vascular resistance and myocardial contractility parameters through matching imaging data of ventricular chambers. Linear regression analysis demonstrated that the lumped-derived RV cardiac index (CI) was in good agreement with catheterization measurements collected from 10 patients with PH (R = 0.82; p < 0.001). Biventricular FE analysis revealed a paradoxical leftward shift of the interventricular septum, and this correlated with invasive measurements of pulmonary vascular resistances (R = 0.70; p = 0.048) as found by Pearson's coefficient. A significant difference was noted for RV myocardial fiber stress in healthy control patients (4.5 ± 0.7 kPa) compared with that of patients with PH at either rest (30.1 ± 12.1 kPa; p = 0.005) or simulated exercise conditions (69.6 ± 24.8 kPa; p < 0.001), thus suggesting adverse RV remodeling. This approach may become a useful and versatile tool for noninvasively assessing RV impairment induced by PH and realistically predicting ventricular mechanics and interactions for an improved management of patients with PH.

Particle image velocimetry study of the celiac trunk hemodynamic induced by continuous-flow left ventricular assist device.

Whereas left ventricular assist device (LVAD) is the gold-standard therapy for patients with heart failure, gastrointestinal bleeding is one of the most common complications. LVAD implantation may remarkably impact aortic hemodynamics so that experimental and computational flow analyses can be used to study the disease mechanisms. Here we present an experimentally-calibrated computational model of the celiac trunk hemodynamic of a LVAD-supported patient who experienced bleeding after device implantation. Specifically, both particle image velocimetry (PIV) and echocardiography were used to measure and compare flow distributions in each branch of a phantom model of the patient abdominal aorta. Then, the distribution of wall shear stress (WSS) was estimated by computational flow analysis. At a cardiac output of 5 L/min, the highest flow division was found in the mesenteric artery (13.6% for PIV and 14.6% for echocardiography), while the left renal artery exhibited the lowest amount in the celiac trunk model (2.6% for PIV and 2.4% for echocardiography). Bland-Altman analysis demonstrated a high agreement between echocardiographic and PIV-related flow measurements, while computational flow analysis revealed that WSS was high in the LVAD graft anastomosis site and just after the ostia of both the celiac trunk and mesenteric artery. This altered shear stress distribution in the celiac trunk may lead to a flow-mediated mechanism of adverse remodeling of the von Willebrand factor and ultimately to gastrointestinal bleeding as seen clinically in this patient.

Computational Fluid-Dynamic Analysis after Carotid Endarterectomy: Patch Graft versus Direct Suture Closure.

BACKGROUND: Closure technique after carotid endarterectomy (CEA) still remains an issue of debate. Routine use of patch graft (PG) has been advocated to reduce restenosis, stroke, and death, but its protective effect, particularly from late restenosis, is less evident and recent studies call into question this thesis. This study aims to compare PG and direct suture (DS) by means of computational fluid dynamics (CFD). To identify carotid regions with flow recirculation more prone to restenosis development, we analyzed time-averaged oscillatory shear index (OSI) and relative residence time (RRT), that are well-known indices correlated with plaque formation. METHODS: CFD was performed in 12 patients (13 carotids) who underwent surgery for stenosis >70%, 9 with PG, and 4 with DS. Flow conditions were modeled using patient-specific boundary conditions derived from Doppler ultrasound and geometries from magnetic resonance angiography. RESULTS: Mean value of the spatial averaged OSI resulted 0.07 for PG group and 0.03 for DS group, the percentage of area with OSI above a threshold of 0.2 resulted 10.1% and 3.7%, respectively. The mean of averaged-in-space RRT values was 4.4 1/Pa for PG group and 1.6 1/Pa for DS group, the percentage of area with RRT values above a threshold of 4 1/Pa resulted 22.5% and 6.5%, respectively. CONCLUSIONS: Both OSI and RRT values resulted higher when PG was preferred to DS and also areas with disturbed flow resulted wider. The absolute higher values computed by means of CFD were observed when PG was used indiscriminately regardless of carotid diameters. DS does not seem to create negative hemodynamic conditions with potential adverse effects on long-term outcomes, in particular when CEA is performed at the common carotid artery and/or the bulb or when ICA diameter is greater than 5.0 mm.

Shear stress alterations in the celiac trunk of patients with a continuous-flow left ventricular assist device as shown by in-silico and in-vitro flow analyses.

BACKGROUND: The use of left ventricular assist devices (LVADs) to treat advanced cardiac heart failure is constantly increasing, although this device leads to high risk for gastrointestinal bleeding. METHODS: Using in-silico flow analysis, we quantified hemodynamic alterations due to continuous-flow LVAD (HeartWare, Inc., Framingham, MA) in the celiac trunk and major branches of the abdominal aorta, and then explored the relationship between wall shear stress (WSS) and celiac trunk orientation. To assess outflow from the aortic branch, a 3-dimensional-printed patient-specific model of the celiac trunk reconstructed from an LVAD-supported patient was used to estimate echocardiographic outflow velocities under continuous-flow conditions, and then to calibrate computational simulations. Moreover, flow pattern and resulting WSS values were computed for 5 patients with LVAD implantation. RESULTS: Peak WSS values were estimated on the 3 branches of the celiac trunk and the LVAD cannula. The mean WSSs demonstrated that the left gastric artery underwent the highest WSS of 9.08 ± 5.45 Pa, with an average flow velocity of 0.57 ± 0.25 m/s compared with that of other vessel districts. The common hepatic artery had a less critical WSS of 4.58 ± 1.77 Pa. A positive correlation was found between the celiac trunk angulation and the WSS stress just distal to the ostium of the celiac trunk (R = 0.9), which may increase vulnerability of this vessel to bleeding. CONCLUSIONS: Although further studies are needed to confirm these findings in a larger patient cohort, computational flow simulations may enhance the information of clinical image data and may have an application in clinical investigations of hemodynamic changes in LVAD-supported patients.

Three-dimensional parametric modeling of bicuspid aortopathy and comparison with computational flow predictions.

Bicuspid aortic valve (BAV)-associated ascending aneurysmal aortopathy (namely "bicuspid aortopathy") is a heterogeneous disease making surgeon predictions particularly challenging. Computational flow analysis can be used to evaluate the BAV-related hemodynamic disturbances, which likely lead to aneurysm enlargement and progression. However, the anatomic reconstruction process is time consuming so that predicting hemodynamic and structural evolution by computational modeling is unfeasible in routine clinical practice. The aim of the study was to design and develop a parametric program for three-dimensional (3D) representations of aneurysmal aorta and different BAV phenotypes starting from several measures derived by computed-tomography angiography (CTA). Assuming that wall shear stress (WSS) has an important implication on bicuspid aortopathy, computational flow analyses were then performed to estimate how different would such an important parameter be, if a parametric aortic geometry was used as compared to standard geometric reconstructions obtained by CTA scans. Morphologic parameters here documented can be used to rapidly model the aorta and any phenotypes of BAV. t-test and Bland-Altman plot demonstrated that WSS obtained by flow analysis of parametric aortic geometries was in good agreement with that obtained from the flow analysis of CTA-related geometries. The proposed program offers a rapid and automated tool for 3D anatomic representations of bicuspid aortopathy with promising application in routine clinical practice by reducing the amount of time for anatomic reconstructions.

Early distal remodeling after elephant trunk repair of thoraco-abdominal aortic aneurysms.

Hemodynamic alterations occur when the elephant trunk (ET) technique is adopted to treat extensive aortic aneurysms. In planning the 2nd stage operation to complete ET repair, surgeons must weigh an adequate recovery time after initial surgery against the risk of postoperative ET-related complications. The purpose of this study was to understand the mechanistic link between the flow alteration caused by the ET graft and the development of premature aortic rupture before the 2nd stage operation. Specifically, fluid-structure interaction (FSI) analysis was performed using the CT imaging data of aorta at different stages of ET repair, and then computational variables were compared to those observed in patients who underwent a prophylactic 2nd stage operation to complete aortic repair. Results show that intramural stress exerted near the distal ET anastomosis (IMS=37.5kPa) at the time of urgent intervention was comparable to that of the extensive aortic aneurysm (IMS=47.4kPa) at initial in-hospital admission, but was considerably higher than that occurring after the 1st stage procedure (IMS=3.5kPa). Pressure index suggested higher peri-graft pressurization than aortic lumen pressure during diastole, imparting an apparent risk of aortic dilatation. These critical hemodynamic and structural parameters are related to the impending rupture of descending aorta observed clinically and can thus guide prophylactic intervention and optimal timing for the 2nd stage operation of a ET technique.

An In Vitro Phantom Study on the Role of the Bird-Beak Configuration in Endograft Infolding in the Aortic Arch.

PURPOSE: To assess endograft infolding for excessive bird-beak configurations in the aortic arch in relation to hemodynamic variables by quantifying device displacement and rotation of oversized stent-grafts deployed in a phantom model. METHODS: A patient-specific, compliant, phantom pulsatile flow model was reconstructed from a patient who presented with collapse of a Gore TAG thoracic endoprosthesis. Device infolding was measured under different flow and pressure conditions for 3 protrusion extensions (13, 19, and 24 mm) of the bird-beak configuration resulting from 2 TAG endografts with oversizing of 11% and 45%, respectively. RESULTS: The bird-beak configuration with the greatest protrusion extension exhibited the maximum TAG device displacement (1.66 mm), while the lowest protrusion extension configuration led to the minimum amount of both displacement and rotation parameters (0.25 mm and 0.6°, respectively). A positive relationship was found between the infolding parameters and the flow circulating in the aorta and left subclavian artery. Similarly, TAG device displacement was positively and significantly (p<0.05) correlated with the pulse pressure for all bird-beak configurations and device sizes. However, no collapse was observed under chronic perfusion testing maintained for 30 days and pulse pressure of 100 mm Hg. CONCLUSION: These findings suggest that endograft infolding depends primarily on the amount of aortic pulsatility and flow rate and that physiological flows do not necessarily engender hemodynamic loads on the proximal bird-beak segment sufficient to cause TAG collapse. Hemodynamic variables may allow for identification of patients at high risk of endograft infolding and help guide preventive intervention to avert its occurrence.

Evaluation of ventricular wall stress and cardiac function in patients with dilated cardiomyopathy.

Dilated cardiomyopathy is a heart disease characterized by both left ventricular dilatation and left ventricular systolic dysfunction, leading to cardiac remodeling and ultimately heart failure. We aimed to investigate the effect of dilated cardiomyopathy on the pump performance and myocardial wall mechanics using patient-specific finite element analysis. Results evinced pronounced end-systolic wall stress on left ventricular wall of patients with dilated cardiomyopathy as compared to that of normal hearts. In dilated cardiomyopathy, both end-diastolic and end-systolic pressure-volume relationships of left ventricle and right ventricle were shifted to the right compared to controls, suggesting reduced myocardial contractility. We hereby propose that finite element analysis represents a useful tool to assess the myocardial wall stress and cardiac work, which are responsible for progressive left ventricular deterioration and poor clinical course.

Biomechanical implications of excessive endograft protrusion into the aortic arch after thoracic endovascular repair.

Endografts placed in the aorta for thoracic endovascular aortic repair (TEVAR) may determine malappositioning to the lesser curvature of the aortic wall, thus resulting in a devastating complication known as endograft collapse. This premature device failure commonly occurs in young individuals after TEVAR for traumatic aortic injuries as a result of applications outside the physical conditions for which the endograft was designed. In this study, an experimentally-calibrated fluid-structure interaction (FSI) model was developed to assess the hemodynamic and stress/strain distributions acting on the excessive protrusion extension (PE) of endografts deployed in four young patients underwent TEVAR. Endograft infolding was experimentally measured for different hemodynamic scenarios by perfusion testing and then used to numerically calibrate the mechanical behavior of endograft PE. Results evinced that the extent of endograft can severely alter the hemodynamic and structural loads exerted on the endograft PE. Specifically, PE determined a physiological aortic coarctation into the aortic arch characterized by a helical flow in the distal descending aorta. High device displacement and transmural pressure across the stent-graft wall were found for a PE longer than 21 mm. Finally, marked intramural stress and principal strain distributions on the protruded segment of the endograft wall may suggest failure due to material fatigue. These critical parameters may contribute to the endograft collapse observed clinically and can be used to design new devices more suitable for young individuals to be treated with an endoprosthesis for TEVAR of blunt traumatic aortic injuries.

Mechanics of pericardial effusion: a simulation study.

Pericardial effusion is a pathological accumulation of fluid within pericardial cavity, which may compress heart chambers with hemodynamic impairment. We sought to determine the mechanics underlying the physiology of the hemodynamic impairment due to pericardial effusion using patient-specific computational modeling. Computational models of left ventricle and right ventricle were based on magnetic resonance images obtained from patients with pericardial effusion and controls. Myocardial material parameters were adjusted, so that volumes of ventricular chambers and pericardial effusion agreed with magnetic resonance imaging data. End-diastolic and end-systolic pressure-volume relationships as well as stroke volume were determined to evaluate impaired cardiac function of biventricular model. Distributions of myocardial fiber stresses and their regional variation along left ventricular wall were compared between patient groups. Both end-diastolic and end-systolic pressure-volume relationships shifted to the left for patients with pericardial effusion, with right ventricle diastolic filling particularly restricted. Left ventricle function as estimated by Starling curve was reduced by pericardial effusion. End-systolic fiber stress of left ventricle was significantly reduced as compared to that found for healthy patients. Myocardial stress was found increased at interventricular septum when compared to that exerted at lateral wall of left ventricle. Right ventricular myocardial stress was reduced as a consequence of the pressure equalization between right ventricle and pericardial effusion. Diastolic right ventricle collapse in patients with pericardial effusion is related to higher myocardial fiber stress on interventricular septum and to an extensible pericardium reducing motion of ventricular chambers, with right ventricle particularly restrained. These findings likely portend progression of pericardial effusion to cardiac tamponade.