Search

Objective: We created a finite element model to predict the probability of dissection based on imaging-derived aortic stiffness and investigated the link between stiffness and wall tensile stress using our model., Methods: Transthoracic echocardiogram measurements were used to calculate aortic diameter change over the cardiac cycle. Aortic stiffness index was subsequently calculated based on diameter change and blood pressure. A series of logistic models were developed to predict the binary outcome of aortic dissection using 1 or more series of predictor parameters such as aortic stiffness index or patient characteristics. Finite element analysis was performed on a subset of diameter-matched patients exhibiting patient-specific material properties., Results: Transthoracic echocardiogram scans of patients with type A aortic dissection (n = 22) exhibited elevated baseline aortic stiffness index when compared with aneurysmal patients' scans with tricuspid aortic valve (n = 83, P < .001) and bicuspid aortic valve (n = 80, P < .001). Aortic stiffness index proved an excellent discriminator for a future dissection event (area under the curve, 0.9337, odds ratio, 2.896). From the parametric finite element study, we found a correlation between peak longitudinal wall tensile stress and stiffness index (ρ = .6268, P < .001, n = 28 pooled)., Conclusions: Noninvasive transthoracic echocardiogram-derived aortic stiffness measurements may serve as an impactful metric toward predicting aortic dissection or quantifying dissection risk. A correlation between longitudinal stress and stiffness establishes an evidence-based link between a noninvasive stiffness parameter and stress state of the aorta with clinically apparent dissection events., (Copyright © 2022 The American Association for Thoracic Surgery. Published by Elsevier Inc. All rights reserved.)

Heart growth in the pregnant patient helps maintain cardiovascular function while supporting the growing fetus. However, in some cases, the cardiovascular demand of pregnancy can trigger life-threatening conditions, including hypertensive disorders of pregnancy and peripartum cardiomyopathy. The mechanisms that control heart growth throughout pregnancy are unclear, and treating these diseases remains elusive. We previously developed a computational model that accounts for hormonal and hemodynamic interactions throughout pregnancy and demonstrated its ability to capture realistic cardiac growth in normal rat pregnancy. In this study, we evaluated whether this model could capture heart growth beyond normal pregnancy. After further validation of our normal pregnancy predictions, we tested our model predictions of three rat studies of hypertensive pregnancies. We next simulated the postpartum period and examined the impact of lactation on cardiac growth in rats. We demonstrate that our multiscale model can capture cardiac growth associated with new-onset hypertension during pregnancy and lactation status in the postpartum period. We conclude by elaborating on the potential clinical utility of our model in the future. NEW & NOTEWORTHY Our multiscale model predicts appropriate heart growth beyond normal pregnancy, including elevated heart weights in rats with induced hypertension during pregnancy and in lactating mice and decreased heart weight in nonlactating mice. Our model captures distinct mechanisms that result in similar organ-level growth, highlighting its potential to distinguish healthy from diseased pregnancy-induced growth.

Objective: This study aims to explore the potential of organic electrolytic photocapacitors (OEPCs), an innovative photovoltaic device, in mediating the activation of native voltage-gated Cav1.2 channels (I Ca,L ) in Guinea pig ventricular cardiomyocytes., Methods: Whole-cell patch-clamp recordings were employed to examine light-triggered OEPC mediated I Ca,L activation, integrating the channel's kinetic properties into a multicompartment cell model to take intracellular ion concentrations into account. A multidomain model was additionally incorporated to evaluate effects of OEPC-mediated stimulation. The final model combines external stimulation, multicompartmental cell simulation, and a patch-clamp amplifier equivalent circuit to assess the impact on achievable intracellular voltage changes., Results: Light pulses activated I Ca,L , with amplitudes similar to voltage-clamp activation and high sensitivity to the L-type Ca 2+ channel blocker, nifedipine. Light-triggered I Ca,L inactivation exhibited kinetic parameters comparable to voltage-induced inactivation., Conclusion: OEPC-mediated activation of I Ca,L demonstrates their potential for nongenetic optical modulation of cellular physiology potentially paving the way for the development of innovative therapies in cardiovascular health. The integrated model proves the light-mediated activation of I Ca,L and advances the understanding of the interplay between the patch-clamp amplifier and external stimulation devices., Significance: Treating cardiac conduction disorders by minimal-invasive means without genetic modifications could advance therapeutic approaches increasing patients' quality of life compared with conventional methods employing electronic devices.

Blood flow regulation within the microvasculature reflects a complex interaction of regulatory mechanisms and varies spatially and temporally according to conditions such as metabolism, growth, injury, and disease. Understanding the role of microvascular flow distributions across conditions is of interest to investigators spanning multiple disciplines; however, data collection within networks can be labor-intensive and challenging due to limited resolution. To overcome these experimental challenges, computational network models that can accurately simulate vascular behavior are highly beneficial. Constrained constructive optimization (CCO) is a commonly used algorithm for vascular simulation, particularly well known for its adaptability toward vascular modeling across tissues. The present work demonstrates an implementation of CCO aimed to simulate a branching arteriolar microvasculature in healthy skeletal muscle, validated against literature including comprehensive rat gluteus maximus vasculature datasets, and reviews a list of user-specified adjustable model parameters to understand how their variability affects the simulated networks. Network geometric properties, including mean element diameters, lengths, and numbers of bifurcations per order, Horton's law ratios, and fractal dimension, demonstrate good validation once model parameters are adjusted to experimental data. This model successfully demonstrates hemodynamic properties such as Murray's law and the network Fahraeus effect. Application of centrifugal and Strahler ordering schemes results in divergent descriptions of identical simulated networks. This work introduces a novel CCO-based model focused on generating branching skeletal muscle microvascular arteriolar networks based on adjustable model parameters, thus making it a valuable tool for investigations into skeletal muscle microvascular structure and tissue perfusion. NEW & NOTEWORTHY The present work introduces a CCO-based algorithm for generating branching arteriolar networks, with adjustable model parameters to enable modeling in varying skeletal muscle tissues. The geometric and hemodynamic parameters of the generated networks have been comprehensively validated using experimental data collected previously in-house and from literature. This is one of few validated CCO-based models to specialize in skeletal muscle microvasculature and acts as a beneficial tool for investigating the microvasculature for hypothesis testing and validation.

We sought to understand how leaflet forces change in response to annular dilation and leaflet tethering (LT) in single ventricle physiology. Explanted fetal bovine tricuspid valves were sutured onto image-derived annuli and ventricular mounts. Control valves (CON) were secured to a size-matched hypoplastic left heart syndrome (HLHS)-type annulus and compared to: (1) normal tricuspid valves secured to a size-matched saddle-shaped annulus, (2) HLHS-type annulus with LT, (3) HLHS-type annulus with annular dilation (dilation valves), or (4) a combined disease model with both dilation and tethering (disease valves). The specimens were tested in a systemic heart simulator at various single ventricle physiologies. Leaflet forces were measured using optical strain sensors sutured to each leaflet edge. Average force in the anterior leaflet was 43.2% lower in CON compared to normal tricuspid valves (P < 0.001). LT resulted in a 6.6% increase in average forces on the anterior leaflet (P = 0.04), 10.7% increase on the posterior leaflet (P = 0.03), and 14.1% increase on the septal leaflet (P < 0.001). In dilation valves, average septal leaflet forces increased relative to the CON by 42.2% (P = 0.01). In disease valves, average leaflet forces increased by 54.8% in the anterior leaflet (P < 0.001), 37.6% in the posterior leaflet (P = 0.03), and 79.9% in the septal leaflet (P < 0.001). The anterior leaflet experiences the highest forces in the normal tricuspid annulus under single ventricle physiology conditions. Annular dilation resulted in an increase in forces on the septal leaflet and LT resulted in an increase in forces across all 3 leaflets. Annular dilation and LT combined resulted in the largest increase in leaflet forces across all 3 leaflets., (Copyright © 2022 Elsevier Inc. All rights reserved.)

In computational fluid dynamic studies related to blood flow, investigating the behavior of blood particles is crucial, especially red blood cells as they constitute a significant proportion of blood particles. Additionally, studying red blood cell movements is necessary, especially in stenotic artery geometries. A new multiphase scheme was utilized to demonstrate the effect of red blood cells on hemodynamics in complex coronary arteries and investigate the consequence of their motion. To investigate the effect of red blood cell movement on flow, the dense discrete phase model (DDPM) was used. This simulation was performed in 3D coronary arteries with different degrees of stenosis, utilizing blood pressure as inlet and outlet boundary conditions while assuming the arterial wall to be rigid. The model prediction shows good agreement with experimental data. Velocity values were comparable in both single-phase and two-phase flow simulations, but the shear stress in two-phase modeling had higher values. In the two-phase DDPM modeling, the recirculation areas indicated a higher probability of atherosclerosis plaque re-formation in the pre-stenosis area compared to the stenosis and post-stenosis areas. The DDPM model was found to be more effective in obtaining shear stress values in the artery. Additionally, this model provides good results compared to the single-phase model in investigating the movement of particles along the artery as well as recirculation areas that lead to the deposition of particles., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships and funding that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Ltd.)

Age-associated cardiovascular diseases (CVDs), predominantly resulting from artery-related disorders such as atherosclerosis, stand as a leading cause of morbidity and mortality among the elderly population. Consequently, there is a growing interest in the development of clinically relevant bioengineered models of CVDs. Recent developments in bioengineering and material sciences have paved the way for the creation of intricate models that closely mimic the structure and surroundings of native cardiac tissues and blood vessels. These models can be utilized for basic research purposes and for identifying pharmaceutical interventions and facilitating drug discovery. The advancement of vessel-on-a-chip technologies and the development of bioengineered and humanized in vitro models of the cardiovascular system have the potential to revolutionize CVD disease modelling. These technologies offer pathophysiologically relevant models at a fraction of the cost and time required for traditional experimentation required in vivo. This progress signifies a significant advancement in the field, transitioning from conventional 2D cell culture models to advanced 3D organoid and vessel-on-a-chip models. These innovative models are specifically designed to explore the complexities of vascular aging and stiffening, crucial factors in the development of cardiovascular diseases. This review summarizes the recent progress of various bioengineered in vitro platforms developed for investigating the pathophysiology of human cardiovascular system with more focus on advanced 3D vascular platforms., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier B.V. All rights reserved.)

Aortic valve replacement has become an increasing concern due to the rising prevalence of aortic stenosis in an ageing population. Existing replacement options have limitations, necessitating the development of improved prosthetic aortic valves. In this study, flow characteristics during systole in a stenotic aortic valve case are compared with those downstream of two newly designed surgical bioprosthetic aortic valves (BioAVs). To do so, advanced three-dimensional fluid-structure interaction simulations are conducted and dedicated analysis methods to investigate jet flow configuration and vortex dynamics are developed. Our findings reveal that the stenotic case maintains a high jet flow eccentricity due to a fixed orifice geometry, resulting in flow separation and increased vortex stretching and tilting in the commissural low-flow regions. One BioAV design introduces non-axisymmetric leaflet motion, which reduces the maximum jet velocity and forms more vortical structures. The other BioAV design produces a fixed symmetric triangular jet shape due to non-moving leaflets and exhibits favourable vorticity attenuation, revealed by negative temporally and spatially averaged projected vortex stretching values, and significantly reduced drag. Therefore, this study highlights the benefits of custom-designed aortic valves in the context of their replacement through comprehensive and novel flow analyses. The results emphasise the importance of analysing jet flow, vortical structures, momentum balance and vorticity transport for thoroughly evaluating aortic valve performance., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Objective: Iatrogenic coronary artery dissection is a complication of coronary intimal injury and dissection due to improper catheter manipulation. The impact of tear direction on the prognosis of coronary artery dissection (CAD) remains unclear. This study examines the hemodynamic effects of different tear directions (transverse and longitudinal) of CAD and evaluates the risk of thrombosis, rupture and further dilatation of CAD., Methods: Two types of CAD models (Type I: transverse tear, Type II: longitudinal tear) were reconstructed from the aorto-coronary CTA dataset of 8 healthy cases. Four WSS-based indicators were analyzed, including time-averaged wall shear stress (TAWSS), oscillatory shear index (OSI), relative residence time (RRT), and cross flow index (CFI). A thrombus growth model was also introduced to predict the trend of thrombus growth in CAD with two different tear directions., Results: For most of the WSS-based indicators, including TAWSS, RRT, and CFI, no statistically significant differences were observed across the CAD models with varying tear directions, except for OSI, where a significant difference was noted (p < 0.05). Meanwhile, in terms of thrombus growth, the thrombus growing at the tear of the Type I (transverse tear) CAD model extended into the true lumen earlier than that of the Type II (longitudinal tear) model., Conclusions: Numerical simulations suggest that: (1) The CAD with transverse tear have a high risk of further tearing of the dissection at the distal end of the tear. (2) The CAD with longitudinal tear create a hemodynamic environment characterized by low TAWSS and high OSI in the false lumen, which may additionally increase the risk of vessel wall injury. (3) The CAD with transverse tear may have a higher risk of thrombosis and coronary obstruction and myocardial ischemia in the early phase of the dissection., Competing Interests: Declaration of competing interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2024. Published by Elsevier B.V.)

Single ventricle patients, including those with hypoplastic left heart syndrome (HLHS), typically undergo three palliative heart surgeries culminating in the Fontan procedure. HLHS is associated with high rates of morbidity and mortality, and many patients develop arrhythmias, electrical dyssynchrony, and eventually ventricular failure. However, the correlation between ventricular enlargement and electrical dysfunction in HLHS physiology remains poorly understood. Here we characterize the relationship between growth and electrophysiology in HLHS using computational modeling. We integrate a personalized finite element model, a volumetric growth model, and a personalized electrophysiology model to perform controlled in silico experiments. We show that right ventricle enlargement negatively affects QRS duration and interventricular dyssynchrony. Conversely, left ventricle enlargement can partially compensate for this dyssynchrony. These findings have potential implications on our understanding of the origins of electrical dyssynchrony and, ultimately, the treatment of HLHS patients.

Despite the key role of fibrosis in atrial fibrillation (AF), the effects of different spatial distributions and textures of fibrosis on wave propagation mechanisms in AF are not fully understood. To clarify these aspects, we performed a systematic computational study to assess fibrosis effects on the characteristics and stability of re-entrant waves in electrically-remodelled atrial tissues. A stochastic algorithm, which generated fibrotic distributions with controlled overall amount, average size, and orientation of fibrosis elements, was implemented on a monolayer spheric atrial model. 245 simulations were run at changing fibrosis parameters. The emerging propagation patterns were quantified in terms of rate, regularity, and coupling by frequency-domain analysis of correspondent synthetic bipolar electrograms. At the increase of fibrosis amount, the rate of reentrant waves significantly decreased and higher levels of regularity and coupling were observed (p < 0.0001). Higher spatial variability and pattern stochasticity over repetitions was observed for larger amount of fibrosis, especially in the presence of patchy and compact fibrosis. Overall, propagation slowing and organization led to higher stability of re-entrant waves. These results strengthen the evidence that the amount and spatial distribution of fibrosis concur in dictating re-entry dynamics in remodeled tissue and represent key factors in AF maintenance., (© 2024. The Author(s).)

Uncertainty quantification is becoming a key tool to ensure that numerical models can be sufficiently trusted to be used in domains such as medical device design. Demonstration of how input parameters impact the quantities of interest generated by any numerical model is essential to understanding the limits of its reliability. With the lattice Boltzmann method now a widely used approach for computational fluid dynamics, building greater understanding of its numerical uncertainty characteristics will support its further use in science and industry. In this study we apply an in-depth uncertainty quantification study of the lattice Boltzmann method in a canonical bifurcating geometry that is representative of the vascular junctions present in arterial and venous domains. These campaigns examine how quantities of interest-pressure and velocity along the central axes of the bifurcation-are influenced by the algorithmic parameters of the lattice Boltzmann method and the parameters controlling the values imposed at inlet velocity and outlet pressure boundary conditions. We also conduct a similar campaign on a set of personalised vessels to further illustrate the application of these techniques. Our work provides insights into how input parameters and boundary conditions impact the velocity and pressure distributions calculated in a simulation and can guide the choices of such values when applied to vascular studies of patient specific geometries. We observe that, from an algorithmic perspective, the number of time steps and the size of the grid spacing are the most influential parameters. When considering the influence of boundary conditions, we note that the magnitude of the inlet velocity and the mean pressure applied within sinusoidal pressure outlets have the greatest impact on output quantities of interest. We also observe that, when comparing the magnitude of variation imposed in the input parameters with that observed in the output quantities, this variability is particularly magnified when the input velocity is altered. This study also demonstrates how open-source toolkits for validation, verification and uncertainty quantification can be applied to numerical models deployed on high-performance computers without the need for modifying the simulation code itself. Such an ability is key to the more widespread adoption of the analysis of uncertainty in numerical models by significantly reducing the complexity of their execution and analysis., (© 2024. The Author(s).)

Background: Integration of a patient's non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019-10-07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013-11-12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient ( ICC ) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error ( χ 2 ) of LV myocardial strain, strain rate, and cavity volume., Results: A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients ( χ 2  < 1.6), but minimum parameter reproducibility was poor ( ICC min  = 0.01). Iterative reduction yielded a reproducible ( ICC min  = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs ( χ 2  < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05)., Conclusions: By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient's underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making., (© 2024. The Author(s).)

The efficacy of flow diverters is influenced by the strut configuration changes resulting from size discrepancies between the stent and the parent artery. This study aimed to quantitatively analyze the impact of size discrepancies between flow diverters and parent arteries on the flow diversion effects, using computational fluid dynamics. Four silicone models with varying parent artery sizes were developed. Real flow diverters were deployed in these models to assess stent configurations at the aneurysm neck. Virtual stents were generated based on these configurations for computational fluid dynamics analysis. The changes in the reduction rate of the hemodynamic parameters were quantified to evaluate the flow diversion effect. Implanting 4.0 mm flow diverters in aneurysm models with parent artery diameters of 3.0-4.5 mm, in 0.5 mm increments, revealed that a shift from oversized to undersized flow diverters led to an increase in the reduction rates of hemodynamic parameter, accompanied by enhanced metal coverage rate and pore density. However, the flow diversion effect observed transitioning from oversizing to matching was less pronounced when moving from matching to undersizing. This emphasizes the importance of proper sizing of flow diverters, considering the benefits of undersizing and not to exceed the threshold of advantages., (© 2024. The Author(s).)

Solitary wave solutions are of great interest to bio-mathematicians and other scientists because they provide a basic description of nonlinear phenomena with many practical applications. They provide a strong foundation for the development of novel biological and medical models and therapies because of their remarkable behavior and persistence. They have the potential to improve our comprehension of intricate biological systems and help us create novel therapeutic approaches, which is something that researchers are actively investigating. In this study, solitary wave solutions of the nonlinear Murray equation will be discovered using a modified extended direct algebraic method. These solutions represent a uniform variation in blood vessel shape and diameter that can be used to stimulate blood flow in patients with cardiovascular disease. These solutions are newly in the literature, and give researchers an important tool for grasping complex biological systems. To see how the solitary wave solutions behave, graphs are displayed using Matlab., (© 2024. The Author(s).)

Competing Interests: Conflict of interest: none declared.

The coronary artery constitutes a vital vascular system that sustains cardiac function, with its primary role being the conveyance of indispensable nutrients to the myocardial tissue. When coronary artery disease occurs, it will affect the blood supply of the heart and induce myocardial ischemia. Therefore, it is of great significance to numerically simulate the coronary artery and evaluate its blood supply capacity. In this article, the coronary artery lumped parameter model was derived based on the relationship between circuit system parameters and cardiovascular system parameters, and the blood supply capacity of the coronary artery in healthy and stenosis states was studied. The aortic root pressure calculated by the aortic valve fluid-structure interaction (AV FSI) simulator was employed as the inlet boundary condition. To emulate the physiological phenomenon of sudden pressure drops resulting from an abrupt reduction in blood vessel radius, a head loss model was connected at the coronary artery's entrance. For each coronary artery outlet, the symmetric structured tree model was appended to simulate the terminal impedance of the missing downstream coronary arteries. The particle swarm optimization (PSO) algorithm was used to optimize the blood flow viscous resistance, blood flow inertia, and vascular compliance of the coronary artery model. In the stenosis states, the relative flow and fractional flow reserve (FFR) calculated by numerical simulation corresponded to the published literature data. It was anticipated that the proposed model can be readily adapted for clinical application, serving as a valuable reference for diagnosing and treating patients.

Objective. Significant aortic regurgitation is a common complication following left ventricular assist device (LVAD) intervention, and existing studies have not attempted to monitor regurgitation signals and undertake preventive measures during full support. Regurgitation is an adverse event that can lead to inadequate left ventricular unloading, insufficient peripheral perfusion, and repeated episodes of heart failure. Moreover, regurgitation occurring during full support due to pump position offset cannot be directly controlled through control algorithms. Therefore, accurate estimation of regurgitation during percutaneous left ventricular assist device (PLVAD) full support is critical for clinical management and patient safety. Approach. An estimation system based on the regurgitation model is built in this paper, and the unscented Kalman filter estimator (UKF) is introduced as an estimation approach. Three offset degrees and three heart failure states are considered in the investigation. Using the mock circulatory loop experimental platform, compare the regurgitation estimated by the UKF algorithm with the actual measured regurgitation; the errors are analyzed using standard confidence intervals of ±2 SDs, and the effectiveness of the mentioned algorithms is thus assessed. The generalization ability of the proposed algorithm is verified by setting different heart failure conditions and different rotational speeds. The root mean square error and correlation coefficient between the estimated and actual values are quantified and the statistical significance of accuracy differences in estimation is illustrated using one-way analysis of variance (One-Way ANOVA), which in turn assessed the accuracy and stability of the UKF algorithm. Main results. The research findings demonstrate that the regurgitation estimation system based on the regurgitation model and UKF can relatively accurately estimate the regurgitation status of patients during PLVAD full support, but the effect of myocardial contractility on the estimation accuracy still needs to be taken into account. Significance. The proposed estimation method in this study provides essential reference information for clinical practitioners, enabling them to promptly manage potential complications arising from regurgitation. By sensitively detecting LVAD adverse events, valuable insights into the performance and reliability of the LVAD device can be obtained, offering crucial feedback and data support for device improvement and optimization., (© 2024 Institute of Physics and Engineering in Medicine.)

Aims: Understanding of the tissue cooling properties of cryoballoon ablation during pulmonary vein (PV) isolation is lacking. The purpose of this study was to delineate the depth of the tissue cooling effect during cryoballoon freezing at the pulmonary venous ostium., Methods and Results: A left atrial-PV model was constructed using a three-dimensional printer with data from a patient to which porcine thigh muscle of various thicknesses could be affixed. The model was placed in a 37°C water tank with a PV water flow at a rate that mimicked biological blood flow. Cryofreezing at the PV ostium was performed five times each for sliced porcine thigh muscle of 2, 4, and 6 mm thickness, and sliced muscle cooling on the side opposite the balloon was monitored. The cooling effect was assessed using the average temperature of 12 evenly distributed thermocouples covering the roof region of the left superior PV. Tissue cooling effects were in the order of the 2, 4, and 6 mm thicknesses, with an average temperature of -41.4 ± 4.2°C for 2 mm, -33.0 ± 4.0°C for 4 mm, and 8.0 ± 8.7°C for 6 mm at 180 s (P for trend <0.0001). In addition, tissue temperature drops were steeper in thin muscle (maximum temperature drop per 5 s: 5.2 ± 0.9°C, 3.9 ± 0.7°C, and 1.3 ± 0.7°C, P for trend <0.0001)., Conclusion: The cooling effect of cryoballoon freezing is weaker in the deeper layers. Cryoballoon ablation should be performed with consideration to myocardial thickness., (© The Author(s) 2024. Published by Oxford University Press on behalf of the European Society of Cardiology.)

Left atrial appendage occlusion (LAAO) is a percutaneous procedure to prevent thromboembolism in patients affected by atrial fibrillation. Despite its demonstrated efficacy, the LAA morphological complexity hinders the procedure, resulting in postprocedural drawbacks (device-related thrombus and peri-device leakage). Local anatomical features may cause difficulties in the device's positioning and affect the effectiveness of the device's implant. The current work proposes a detailed FE model of the LAAO useful to investigate implant scenarios and derive clinical indications. A high-fidelity model of the Watchman FLX device and simplified parametric conduits mimicking the zone of the LAA where the device is deployed were developed. Device-conduit interactions were evaluated by looking at clinical indicators such as device-wall gap, possible cause of leakage, and device protrusion. As expected, the positioning of the crimped device before the deployment was found to significantly affect the implant outcomes: clinician's choices can be improved if FE models are used to optimize the pre-operative planning. Remarkably, also the wall mechanical stiffness plays an important role. However, this parameter value is unknown for a specific LAA, a crucial point that must be correctly defined for developing an accurate FE model. Finally, numerical simulations outlined how the device's configuration on which the clinician relies to assess the implant success (i.e., the deployed configuration with the device still attached to the catheter) may differ from the actual final device's configuration, relevant for achieving a safe intervention., (© 2024 John Wiley & Sons Ltd.)

Background: This study investigates the hypothesis that presence of atrial fibrillation (AF) in LVAD patients increases thrombogenicity in the left ventricle (LV) and exacerbates stroke risk., Methods: Using an anatomical LV model implanted with an LVAD inflow cannula, we analyze thrombogenic risk and blood flow patterns in either AF or sinus rhythm (SR) using unsteady computational fluid dynamics (CFD). To analyze platelet activation and thrombogenesis in the LV, hundreds of thousands of platelets are individually tracked to quantify platelet residence time (RT) and shear stress accumulation history (SH)., Results: The irregular and chaotic mitral inflow associated with AF results in markedly different intraventricular flow patterns, with profoundly negative impact on blood flow-induced stimuli experienced by platelets as they traverse the LV. Twice as many platelets accumulated very high SH in the LVAD + AF case, resulting in a 36% increase in thrombogenic potential score, relative to the LVAD + SR case., Conclusions: This supports the hypothesis that AF results in unfavorable blood flow patterns in the LV adding to an increased stroke risk for LVAD + AF patients. Quantification of thrombogenic risk associated with AF for LVAD patients may help guide clinical decision-making on interventions to mitigate the increased risk of thromboembolic events., Competing Interests: Declaration of conflicting interestsThe author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: CM is an investigator/consultant for Abbott, Abiomed, and Endotronix. VKC, JB, SL, NA, and AA have no disclosures.

This study aims to extend earlier Krogh Cylinder Models of an oxygen profile by considering axial diffusion and analytically solving Fick's Law Partial Differential Equation with novel boundary conditions via the separation of variables. We next prospectively collected a total of 20 animals, which were randomly assigned to receive either fresh or two-week-old stored red blood cell (RBC) transfusions and PQM oxygen data were measured acutely (90 min) or chronically (24 h). Transfusion effects were evaluated in vivo using intravital microscopy of the dorsal skinfold window chamber in Golden Syrian Hamsters. Hamsters were initially hemorrhaged by 50% of total blood volume and resuscitated 1-h post hemorrhage. PQM data were subsequently collected and fit the derived 2D Krogh cylinder model. Systemic hemodynamics (mean arterial pressure, heart rate) were similar in both pre and post-transfusion with either stored or fresh cells. Transfusion with stored cells was found to impair axial and radial oxygen gradients as quantified by our model and consistent with previous studies. Specifically, we observed a statistically significant decrease in the arteriolar tissue radial oxygen gradient after transfusion with stored RBCs at 24 h compared with fresh RBCs (0.33 ± 0.17 mmHg μ m-1 vs, 0.14 ± 0.12 mmHg μ m-1; p = 0.0280). We also observed a deficit in the arteriolar tissue oxygen gradient (0.03 ± 0.01 mmHg μ m-1 fresh vs. 0.018 ± 0.007 mmHg μ m-1 stored; p = 0.0185). We successfully derived and validated an analytical 2D Krogh cylinder model in an animal model of microhemodynamic oxygen diffusion aberration secondary to storage lesions., Competing Interests: Declaration of competing interest None Declared., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Mechanical circulatory support (MCS) devices can bridge the gap to transplant whilst awaiting a viable donor heart. The Realheart Total Artificial Heart is a novel positive-displacement MCS that generates pulsatile flow via bileaflet mechanical valves. This study developed a combined computational fluid dynamics and fluid-structure interaction (FSI) methodology for simulating positive displacement bileaflet valves. Overset meshing discretised the fluid domain, and a blended weak-strong coupling FSI algorithm was combined with variable time-stepping. Four operating conditions of relevant stroke lengths and rates were assessed. The results demonstrated this modelling strategy is stable and efficient for modelling positive-displacement artificial hearts.

Background: Peripheral intravenous catheters (PIVCs) are the most commonly used invasive medical device, yet despite best efforts by end-users, PIVCs experience unacceptably high early failure rates. We aimed to design a new PIVC that reduces the early failure rate of in-dwelling PIVCs and we conducted preliminary tests to assess its efficacy and safety in a porcine model of intravenous access., Methods: We used computer-aided design and simulation to create a PIVC with a ramped tip geometry, which directs the infused fluid away from the vein wall; we called the design the FloRamp™. We created FloRamp prototypes (test device) and tested them against a market-leading device (BD Insyte™; control device) in a highly-controlled setting with five insertion sites per device in four pigs. We measured resistance to infusion and visual infusion phlebitis (VIP) every 6 h and terminated the experiment at 48 h. Veins were harvested for histology and seven pathological markers were assessed., Results: Computer simulations showed that the optimum FloRamp tip reduced maximum endothelial shear stress by 60%, from 12.7 Pa to 5.1 Pa, compared to a typical PIVC tip and improved the infusion dynamics of saline in the blood stream. In the animal study, we found that 2/5 of the control devices were occluded after 24 h, whereas all test devices remained patent and functional. The FloRamp created less resistance to infusion (0.73 ± 0.81 vs 0.47 ± 0.50, p  = 0.06) and lower VIP scores (0.60 ± 0.93 vs 0.31 ± 0.70, p  = 0.09) than the control device, although neither findings were significantly different. Histopathology revealed that 5/7 of the assessed markers were lower in veins with the FloRamp., Conclusions: Herein we report preliminary assessment of a novel PIVC design, which could be advantageous in clinical settings through decreased device occlusion and reduced early failure rates., Competing Interests: Declaration of conflicting interestsThe author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: BJD, LJK and CS are named inventors on a patent describing the ramped tip design (WO/2020/237286). The study was funded by Flomatrix Pty Ltd of which CS is an employee and in which both BJD and CS hold equity. SK reports monies received by her employer QUT on her behalf from for educational consultancies with BD Medical and ITL Biomedical unrelated to this study. The remaining authors have no conflicts to report.

Background: Nutcracker syndrome is a disease characterized by complex symptoms, making its diagnosis challenging and often delayed, often resulting in a painful experience for the patients., Objective: This study aimed to investigate the pathogenesis of nutcracker syndrome through the perspective of hemodynamics by simulating blood flow with varying compression degrees of the left renal vein., Methods: 3D patient-specific vascular models of the abdominal aorta, superior mesenteric artery and left renal vein were constructed based on CT images of patients suspected of having nutcracker syndrome. A hemodynamic simulation was then conducted using computational fluid dynamics to identify the correlation between alterations in hemodynamic parameters and varying degrees of compression., Results: The study indicated the presence of an evident gradient in velocity distribution over the left renal vein with relatively high degrees of stenosis (α ≤ 50°), with maximum velocity in the central region of the stenosis. Additionally, when the compression degree of the left renal vein increases, the pressure distribution of the left renal vein presents an increasing number of gradient layers. Furthermore, the wall shear stress shows a correlation with the variation of blood flow velocity, i.e., the increase of wall shear stress correlates with the acceleration of the blood flow velocity., Conclusions: Using computational fluid dynamics as a non-invasive instrument to obtain the hemodynamic characteristics of nutcracker syndrome is feasible and could provide insights into the pathological mechanisms of the nutcracker syndrome supporting clinicians in diagnosis., (© 2024. The Author(s) under exclusive licence to Italian Society of Nephrology.)

Congenital arterial stenosis such as supravalvar aortic stenosis (SVAS) are highly prevalent in Williams syndrome (WS) and other arteriopathies pose a substantial health risk. Conventional tools for severity assessment, including clinical findings and pressure gradient estimations, often fall short due to their susceptibility to transient physiological changes and disease stage influences. Moreover, in the pediatric population, the severity of these and other congenital heart defects (CHDs) often restricts the applicability of invasive techniques for obtaining crucial physiological data. Conversely, evaluating CHDs and their progression requires a comprehensive understanding of intracardiac blood flow. Current imaging modalities, such as blood speckle imaging (BSI) and four-dimensional magnetic resonance imaging (4D MRI) face limitations in resolving flow data, especially in cases of elevated flow velocities. To address these challenges, we devised a computational framework employing zero-dimensional (0D) lumped parameter models coupled with patient-specific reconstructed geometries pre- and post-surgical intervention to execute computational fluid dynamic (CFD) simulations. This framework facilitates the analysis and visualization of intricate blood flow patterns, offering insights into geometry and flow dynamics alterations impacting cardiac function. In this study, we aim to assess the efficacy of surgical intervention in correcting an extreme aortic defect in a patient with WS, leading to reductions in wall shear stress (WSS), maximum velocity magnitude, pressure drop, and ultimately a decrease in cardiac workload., Competing Interests: Declaration of competing interest The authors declare no conflicts of interest., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

A profound investigation of the interaction mechanics between blood vessels and guidewires is necessary to achieve safe intervention. An interactive force model between guidewires and blood vessels is established based on cardiovascular fluid dynamics theory and contact mechanics, considering two intervention phases (straight intervention and contact intervention at a corner named "J-vessel"). The contributing factors of the force model, including intervention conditions, guidewire characteristics, and intravascular environment, are analyzed. A series of experiments were performed to validate the availability of the interactive force model and explore the effects of influential factors on intervention force. The intervention force data were collected using a 2-DOF mechanical testing system instrumented with a force sensor. The guidewire diameter and material were found to significantly impact the intervention force. Additionally, the intervention force was influenced by factors such as blood viscosity, blood vessel wall thickness, blood flow velocity, as well as the interventional velocity and interventional mode. The experiment of the intervention in a coronary artery physical vascular model confirms the practicality validation of the predicted force model and can provide an optimized interventional strategy for vascular interventional surgery. The enhanced intervention strategy has resulted in a considerable reduction of approximately 21.97 % in the force exerted on blood vessels, effectively minimizing the potential for complications associated with the interventional surgery., Competing Interests: Declaration of competing interest None declared., (Copyright © 2024. Published by Elsevier Ltd.)

Repaired Tetralogy of Fallot (rTOF) patients suffer from pulmonary regurgitation and may require pulmonary valve replacement (PVR). Cardiac magnetic resonance imaging (cMRI) guides therapy, but conventional measurements do not quantify the intracardiac flow effects from pulmonary regurgitation or PVR. This study investigates intracardiac flow parameters of the right ventricle (RV) of rTOF by computational fluid dynamics (CFD). cMRI of rTOF patients and controls were retrospectively included. Feature-tracking captured RV endocardial contours from long-axis/short-axis cine. Ventricular motion was reconstructed via diffeomorphic mapping, serving as domain boundary for CFD simulations. Vorticity (1/s), viscous energy loss (ELoss, mJ/L) and turbulent kinetic energy (TKE, mJ/L) were quantified in RV outflow tract (RVOT) and RV inflow. These parameters were normalized against total RV kinetic energy (KE) and RV inflow vorticity to derive dimensionless metrics. Vorticity contours by Q-criterion were qualitatively compared. rTOF patients (n = 15) had mean regurgitant fraction 38 ± 12% and RV size 162 ± 35 mL/m 2 . Compared to controls (n = 12), rTOF had increased RVOT vorticity (142.6 ± 75.6/s vs. 40.4 ± 11.8/s, p < 0.0001), Eloss (55.6 ± 42.5 vs. 5.2 ± 4.4 mJ/L, p = 0.0004), and TKE (5.7 ± 5.9 vs. 0.84 ± 0.46 mJ/L, p = 0.0003). After PVR, there was decrease in normalized RVOT Eloss/TKE (p = 0.0009, p = 0.029) and increase in normalized tricuspid inflow vorticity/KE (p = 0.0136, p = 0.043), corresponding to reorganization of the "donut"-shaped tricuspid ring-vortex. The intracardiac flow in rTOF patients can be simulated to determine the impact of PVR and improve the clinical indications guided by cardiac imaging., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

The progression of cardiovascular disease is intricately influenced by a complex interplay between physiological pathways, biochemical processes, and physical mechanisms. This study aimed to develop an in-silico physics-based approach to comprehensively model the multifaceted vascular pathophysiological adaptations. This approach focused on capturing the progression of proximal pulmonary arterial hypertension, which is significantly associated with the irreversible degradation of arterial walls and compensatory stress-induced growth and remodeling. This study incorporated critical characteristics related to the distinct time scales for the deformation, thus reflecting the impact of mean pressure on artery growth and tissue damage. The in-silico simulation of the progression of pulmonary hypertension was realized based on computational code combined with the finite element method (FEM) for the simulation of disease progression. The parametric studies further explored the consequences of these irreversible processes. This computational modeling approach may advance our understanding of pulmonary hypertension and its progression., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Background: Percutaneous pulmonary valve implantation (PPVI) is a non-surgical treatment for right ventricular outflow tract (RVOT) dysfunction. During PPVI, a stented valve, delivered via catheter, replaces the dysfunctional pulmonary valve. Stent oversizing allows valve anchoring within the RVOT, but overexpansion can intrude on the surrounding structures. Potentially dangerous outcomes include aortic valve insufficiency (AVI) from aortic root (AR) distortion and myocardial ischemia from coronary artery (CA) compression. Currently, risks are evaluated via balloon angioplasty/sizing before stent deployment. Patient-specific finite element (FE) analysis frameworks can improve pre-procedural risk assessment, but current methods require hundreds of hours of high-performance computation., Methods: We created a simplified method to simulate the procedure using patient-specific FE models for accurate, efficient pre-procedural PPVI (using balloon expandable valves) risk assessment. The methodology was tested by retrospectively evaluating the clinical outcome of 12 PPVI candidates., Results: Of 12 patients (median age 14.5 years) with dysfunctional RVOT, 7 had native RVOT and 5 had RV-PA conduits. Seven patients had undergone successful RVOT stent/valve placement, three had significant AVI on balloon testing, one had left CA compression, and one had both AVI and left CA compression. A model-calculated change of more than 20% in lumen diameter of the AR or coronary arteries correctly predicted aortic valve sufficiency and/or CA compression in all the patients., Conclusion: Agreement between FE results and clinical outcomes is excellent. Additionally, these models run in 2-6 min on a desktop computer, demonstrating potential use of FE analysis for pre-procedural risk assessment of PPVI in a clinically relevant timeframe., (© 2024 The Authors. Catheterization and Cardiovascular Interventions published by Wiley Periodicals LLC.)

Blood vessels are subjected to complex biomechanical loads, primarily from pressure-driven blood flow. Abnormal loading associated with vascular grafts, arising from altered hemodynamics or wall mechanics, can cause acute and progressive vascular failure and end-organ dysfunction. Perturbations to mechanobiological stimuli experienced by vascular cells contribute to remodeling of the vascular wall via activation of mechanosensitive signaling pathways and subsequent changes in gene expression and associated turnover of cells and extracellular matrix. In this review, we outline experimental and computational tools used to quantify metrics of biomechanical loading in vascular grafts and highlight those that show potential in predicting graft failure for diverse disease contexts. We include metrics derived from both fluid and solid mechanics that drive feedback loops between mechanobiological processes and changes in the biomechanical state that govern the natural history of vascular grafts. As illustrative examples, we consider application-specific coronary artery bypass grafts, peripheral vascular grafts, and tissue-engineered vascular grafts for congenital heart surgery as each of these involves unique circulatory environments, loading magnitudes, and graft materials., Competing Interests: Disclosures None.

The substantial computational cost of high-fidelity models in numerical hemodynamics has, so far, relegated their use mainly to offline treatment planning. New breakthroughs in data-driven architectures and optimization techniques for fast surrogate modeling provide an exciting opportunity to overcome these limitations, enabling the use of such technology for time-critical decisions. We discuss an application to the repair of multiple stenosis in peripheral pulmonary artery disease through either transcatheter pulmonary artery rehabilitation or surgery, where it is of interest to achieve desired pressures and flows at specific locations in the pulmonary artery tree, while minimizing the risk for the patient. Since different degrees of success can be achieved in practice during treatment, we formulate the problem in probability, and solve it through a sample-based approach. We propose a new offline-online pipeline for probabilistic real-time treatment planning which combines offline assimilation of boundary conditions, model reduction, and training dataset generation with online estimation of marginal probabilities, possibly conditioned on the degree of augmentation observed in already repaired lesions. Moreover, we propose a new approach for the parametrization of arbitrarily shaped vascular repairs through iterative corrections of a zero-dimensional approximant. We demonstrate this pipeline for a diseased model of the pulmonary artery tree available through the Vascular Model Repository., (© 2024 John Wiley & Sons Ltd.)

The reduced-order lumped parameter model (LPM) has great computational efficiency in real-time numerical simulations of haemodynamics but is limited by the accuracy of patient-specific computation. This study proposed a method to achieve the individual LPM modeling with high accuracy to improve the practical clinical applicability of LPM. Clinical data was collected from two medical centres comprising haemodynamic indicators from 323 individuals, including brachial artery pressure waveforms, cardiac output data, and internal carotid artery flow waveforms. The data were expanded to 5000 synthesised cases that all fell within the physiological range of each indicator. LPM of the human blood circulation system was established. A double-path neural network (DPNN) was designed to input the waveforms of each haemodynamic indicator and their key features and then output the individual parameters of the LPM, which was labelled using a conventional optimization algorithm. Clinically collected data from the other 100 cases were used as the test set to verify the accuracy of the individual LPM parameters predicted by DPNN. The results show that DPNN provided good convergence in the training process. In the test set, compared with clinical measurements, the mean differences between each haemodynamic indicator and the estimate calculated by the individual LPM based on the DPNN were about 10 %. Furthermore, DPNN prediction only takes 4 s for 100 cases. The DPNN proposed in this study permits real-time and accurate individualization of LPM's. When facing medical issues involving haemodynamics, it lays the foundation for patient-specific numerical simulation, which may be beneficial for potential clinical application., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 Elsevier Ltd. All rights reserved.)

Arteriosclerosis is a major risk factor for cardiovascular disease and results in arterial vessel stiffening. Velocity estimation of the pulse wave sent by the heart and propagating into the arteries is a widely accepted biomarker. This symmetrical pulse wave propagates at a speed which is related to the Young's modulus through the Moens Korteweg (MK) equation. Recently, an antisymmetric flexural wave has been observed in vivo. Unlike the symmetrical wave, it is highly dispersive. This property offers promising applications for monitoring arterial stiffness and early detection of atheromatous plaque. However, as far as it is known, no equivalent of the MK equation exists for flexural pulse waves. To bridge this gap, a beam based theory was developed, and approximate analytical solutions were reached. An experiment in soft polymer artery phantoms was built to observe the dispersion of flexural waves. A good agreement was found between the analytical expression derived from beam theory and experiments. Moreover, numerical simulations validated wave speed dependence on the elastic and geometric parameters at low frequencies. Clinical applications, such as arterial age estimation and arterial pressure measurement, are foreseen., (© 2024 Acoustical Society of America.)

Iatrogenic damage to the cardiac conduction system (CCS) remains a significant risk during congenital heart surgery. Current surgical best practice involves using superficial anatomical landmarks to locate and avoid damaging the CCS. Prior work indicates inherent variability in the anatomy of the CCS and supporting tissues. This study introduces high-resolution, 3D models of the CCS in normal pediatric human hearts to evaluate variability in the nodes and surrounding structures. Human pediatric hearts were obtained with an average donor age of 2.7 days. A pipeline was developed to excise, section, stain, and image atrioventricular (AVN) and sinus nodal (SN) tissue regions. A convolutional neural network was trained to enable precise multi-class segmentation of whole-slide images, which were subsequently used to generate high- resolution 3D tissue models. Nodal tissue region models were created. All models (10 AVN, 8 SN) contain tissue composition of neural tissue, vasculature, and nodal tissues at micrometer resolution. We describe novel nodal anatomical variations. We found that the depth of the His bundle in females was on average 304 μm shallower than those of male patients. These models provide surgeons with insight into the heterogeneity of the nodal regions and the intricate relationships between the CCS and surrounding structures., Competing Interests: Declaration of competing interest None., (Copyright © 2024 Elsevier Inc. All rights reserved.)

The dialysis catheter indwelling in human bodies has a high risk of inducing thrombus and stenosis. Biomechanical research showed that such physiological complications are triggered by the wall shear stress of the vascular vessel. This study aimed to assess the impact of CVC implantation on central venous haemodynamics and the potential alterations in the haemodynamic environment related to thrombus development. The SVC structure was built from the images from computed tomography. The blood flow was calculated using the Carreau model, and the fluid domain was determined by CFD. The vascular wall and the CVC were computed using FEA. The elastic interaction between the vessel wall and the flow field was considered using FSI simulation. With consideration of the effect of coupling, it was shown that the catheter vibrated in the vascular systems due to the periodic variation of blood pressure, with an amplitude of up to 10% of the vessel width. Spiral flow was observed along the catheter after CVC indwelling, and recirculation flow appeared near the catheter tip. High OSI and WSS regions occurred at the catheter tip and the vascular junction. The arterial lumen tip had a larger effect on the WSS and OSI values on the vascular wall. Considering FSI simulation, the movement of the catheter inside the blood flow was simulated in the deformable vessel. After CVC indwelling, spiral flow and recirculation flow were observed near the regions with high WSS and OSI values., (© 2024 John Wiley & Sons Ltd.)

Background: Pulsed Field Ablation (PFA) has recently been proposed as a non-thermal energy to treat atrial fibrillation by selective ablation of ganglionated plexi (GP) embedded in epicardial fat. While some of PFA-technologies use an endocardial approach, others use epicardial access with promising pre-clinical results. However, as each technology uses a different and sometimes proprietary pulse application protocol, the comparation between endocardial vs. epicardial approach is almost impossible in experimental terms. For this reason, our study, based on a computational model, allows a direct comparison of electric field distribution and thermal-side effects of both approaches under equal conditions in terms of electrode design, pulse protocol and anatomical characteristics of the tissues involved., Methods: 2D computational models with axial symmetry were built for endocardial and epicardial approaches. Atrial (1.5-2.5 mm) and fat (1-5 mm) thicknesses were varied to simulate a representative sample of what happens during PFA ablation for different applied voltage values (1000, 1500 and 2000 V) and number of pulses (30 and 50)., Results: The epicardial approach was superior for capturing greater volumes of fat when the applied voltage was increased: 231 mm 3 /kV with the epicardial approach vs. 182 mm 3 /kV with the endocardial approach. In relation to collateral damage to the myocardium, the epicardial approach considerably spares the myocardium, unlike what happens with the endocardial approach. Although the epicardial approach caused much more thermal damage in the fat, there is not a significant difference between the approaches in terms of size of thermal damage in the myocardium., Conclusions: Our results suggest that epicardial PFA ablation of GPs is more effective than an endocardial approach. The proximity and directionality of the electric field deposited using an epicardial approach are key to ensuring that higher electric field strengths and increased temperatures are obtained within the epicardial fat, thus contributing to selective ablation of the GPs with minimal myocardial damage., Competing Interests: Declaration of competing interest Barry O'Brien is employee of AtriAN Medical, Galway, Ireland. The rest of the authors declare no conflict of interest., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Objectives: To investigate the risk factors and biomechanical mechanisms of in-stent restenosis (ISR) in central venous occlusive disease (CVOD)., Patients and Methods: This retrospective study consecutively included 77 maintenance hemodialysis (MHD) patients with CVOD who received the first percutaneous transluminal angioplasty with stenting (PTS) due to symptomatic CVOD in a tertiary hospital. The mean age was 59.7 ± 14.0 years, and 51.9% of patients were male. The clinical characteristics, occurrence of ISR and patency rates were recorded. Finite element method was applied to assess the biomechanical properties of stents., Results: Among 77 patients with a mean CVS score of 8.0 ± 2.8, 20.8%, 62.3%, and 16.9% of patients had the main vein of CVOD in the subclavian vein, brachiocephalic vein, and superior vena cava, respectively. A total of 72 (93.5%) patients received successful PTS treatment, for which the stents implanted were mainly Fluency covered stent (48.6%) and SMART bare stent (31.9%). During 15 (10-24)-months of follow-up, ISR occurred in 36.1% of the 72 patients. The primary and assisted primary patency rates at 6, 12, and 18 months were 78%, 56%, 42% and 95%, 90%, 87%, respectively. A prolonged dialysis vintage was an independent risk factor for ISR, yet the stent type or the main vein location was not correlated with ISR. Among three laser-engraving stents, the SMART stent was the best in terms of flexibility, stress, and strain on stents but worst in stress or strain on vessels. The Luminexx stent was the best in radial force and worst in stress or strain on stents. The Vici stent was the best in stress and strain on vessels and worst in radial force and flexibility., Conclusions: An unsatisfactory comprehensive biomechanical performance from configurations rooted in existing stents may account for the high incidence of ISR in CVOD., Competing Interests: Declaration of conflicting interestsThe author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.

Computational fluid dynamics (CFD) is an emerging technology applied in the field of cardiovascular medicine, which can obtain hemodynamic data by simulating the blood flow in the patient's heart for cardiac function assessment and disease diagnosis. Left ventricular function plays a key role in the occurrence and development of cardiomyopathies and coronary disease. CFD can reconstruct the left ventricular anatomic structures of patients to clarify pathophysiologic mechanisms and analyze hemodynamic parameters to evaluate left ventricular function, verify surgical efficacy, and guide surgical strategy, which has a positive effect on achieving early diagnosis and reducing mortality from cardiomyopathies and coronary disease. At present, there are still technical limitations in the large-scale clinical application of CFD, and various solutions are being developed and tested, and further improvement and refinement are needed.

Computational fluid dynamics (CFD) is a valuable asset for patient-specific cardiovascular-disease diagnosis and prognosis, but its high computational demands hamper its adoption in practice. Machine-learning methods that estimate blood flow in individual patients could accelerate or replace CFD simulation to overcome these limitations. In this work, we consider the estimation of vector-valued quantities on the wall of three-dimensional geometric artery models. We employ group-equivariant graph convolution in an end-to-end SE(3)-equivariant neural network that operates directly on triangular surface meshes and makes efficient use of training data. We run experiments on a large dataset of synthetic coronary arteries and find that our method estimates directional wall shear stress (WSS) with an approximation error of 7.6% and normalised mean absolute error (NMAE) of 0.4% while up to two orders of magnitude faster than CFD. Furthermore, we show that our method is powerful enough to accurately predict transient, vector-valued WSS over the cardiac cycle while conditioned on a range of different inflow boundary conditions. These results demonstrate the potential of our proposed method as a plugin replacement for CFD in the personalised prediction of hemodynamic vector and scalar fields., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Jelmer M. Wolterink reports financial support was provided by Dutch Research Council. Pim de Haan reports a relationship with Qualcomm Technologies, Inc. that includes: employment., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Cardiac-coronary interaction is fundamental to the function of the heart. As one of the highest metabolic organs in the body, the cardiac oxygen demand is met by blood perfusion through the coronary vasculature. The coronary vasculature is largely embedded within the myocardial tissue which is continually contracting and hence squeezing the blood vessels. The myocardium-coronary vessel interaction is two-ways and complex. Here, we review the different types of cardiac-coronary interactions with a focus on insights gained from mathematical models. Specifically, we will consider the following: (1) myocardial-vessel mechanical interaction; (2) metabolic-flow interaction and regulation; (3) perfusion-contraction matching, and (4) chronic interactions between the myocardium and coronary vasculature. We also provide a discussion of the relevant experimental and clinical studies of different types of cardiac-coronary interactions. Finally, we highlight knowledge gaps, key challenges, and limitations of existing mathematical models along with future research directions to understand the unique myocardium-coronary coupling in the heart. This article is categorized under: Cardiovascular Diseases > Computational Models Cardiovascular Diseases > Biomedical Engineering Cardiovascular Diseases > Molecular and Cellular Physiology., (© 2024 The Authors. WIREs Mechanisms of Disease published by Wiley Periodicals LLC.)

Objective: Arteriovenous malformation is a disease of the vascular system that occurs mainly in the cerebral arteries and spine. Numerical simulation as a powerful method is used to investigate the Cerebral Arteriovenous Malformation hemodynamic after occlusion of the abnormality step by step by embolization., Methods: The computed tomography (CT) Angiographic imaging data of 2 patients are used and a geometric model is extracted by the Mimics software. Numerical simulation of blood flow is performed in both Newtonian and non-Newtonian models. The Navier-Stokes and continuity governing equations are solved by a finite element method using the COMSOL Multiphysics software (the commercial computational fluid dynamics (CFD) simulation software). To validate the numerical results, the real data on blood flow rate in the feeding artery and draining veins are used, as well as angiographic images at different times., Results: Regarding the comparison of pressure contours for different occlusions of 0, 30, 50, and 90%, by increasing the amount of occlusion in the nidus, there is an increase in the blood pressure. Regarding the comparison of the blood flow velocity in the feeding artery, draining veins, and inside the AVM nidus for Newtonian and non-Newtonian models, there is a significant difference between these 2 simulations in vessels with smaller dimensions (such as vessels inside the nidus)., Conclusions: By increasing the amount of nidus occlusion, the blood pressure is increased, so the blood supply process is better. According to a significant difference between the Newtonian and non-Newtonian simulations in vessels with smaller dimensions (such as vessels inside the nidus), therefore, non-Newtonian simulation should be done for different occlusions of 30, 50, and 90%., (Copyright © 2024. Published by Elsevier Inc.)

Importance: Accurate risk prediction of morbidity and mortality in patients with heart failure with preserved ejection fraction (HFpEF) may help clinicians risk stratify and inform care decisions., Objective: To develop and validate a novel prediction model for clinical outcomes in patients with HFpEF using routinely collected variables and to compare it with a biomarker-driven approach., Design, Setting, and Participants: Data were used from the Dapagliflozin Evaluation to Improve the Lives of Patients With Preserved Ejection Fraction Heart Failure (DELIVER) trial to derive the prediction model, and data from the Angiotensin Receptor Neprilysin Inhibition in Heart Failure With Preserved Ejection Fraction (PARAGON-HF) and the Irbesartan in Heart Failure With Preserved Ejection Fraction Study (I-PRESERVE) trials were used to validate it. The outcomes were the composite of HF hospitalization (HFH) or cardiovascular death, cardiovascular death, and all-cause death. A total of 30 baseline candidate variables were selected in a stepwise fashion using multivariable analyses to create the models. Data were analyzed from January 2023 to June 2023., Exposures: Models to estimate the 1-year and 2-year risk of cardiovascular death or hospitalization for heart failure, cardiovascular death, and all-cause death., Results: Data from 6263 individuals in the DELIVER trial were used to derive the prediction model and data from 4796 individuals in the PARAGON-HF trial and 4128 individuals in the I-PRESERVE trial were used to validate it. The final prediction model for the composite outcome included 11 variables: N-terminal pro-brain natriuretic peptide (NT-proBNP) level, HFH within the past 6 months, creatinine level, diabetes, geographic region, HF duration, treatment with a sodium-glucose cotransporter 2 inhibitor, chronic obstructive pulmonary disease, transient ischemic attack/stroke, any previous HFH, and heart rate. This model showed good discrimination (C statistic at 1 year, 0.73; 95% CI, 0.71-0.75) in both validation cohorts (C statistic at 1 year, 0.71; 95% CI, 0.69-0.74 in PARAGON-HF and 0.75; 95% CI, 0.73-0.78 in I-PRESERVE) and calibration. The model showed similar discrimination to a biomarker-driven model including high-sensitivity cardiac troponin T and significantly better discrimination than the Meta-Analysis Global Group in Chronic (MAGGIC) risk score (C statistic at 1 year, 0.60; 95% CI, 0.58-0.63; delta C statistic, 0.13; 95% CI, 0.10-0.15; P < .001) and NT-proBNP level alone (C statistic at 1 year, 0.66; 95% CI, 0.64-0.68; delta C statistic, 0.07; 95% CI, 0.05-0.08; P < .001). Models derived for the prediction of all-cause and cardiovascular death also performed well. An online calculator was created to allow calculation of an individual's risk., Conclusions and Relevance: In this prognostic study, a robust prediction model for clinical outcomes in HFpEF was developed and validated using routinely collected variables. The model performed better than NT-proBNP level alone. The model may help clinicians to identify high-risk patients and guide treatment decisions in HFpEF.

The study aimed to investigate the mechanical factors for distal stent graft-induced new entry (dSINE) in aortic dissection patients and discussed these factors in conjunction with aortic morphology. Two patients (one dSINE and one non-dSINE), with the same age, gender, and type of implanted stent, were selected, then aortic morphological parameters were calculated. In addition, the stent material parameters used by the patients were also fitted. Simulations were performed based on the patient's aortic model and the stent graft used. The true lumen segment at the distal stent graft was designated as the "dSINE risk zone," and mechanical parameters (maximum principal strain, maximum principal stress) were computed. When approaching the area with higher mechanical parameters in the dSINE risk zone, dSINE patient exhibited higher values and growth rates in mechanical parameters compared to non-dSINE patient. Furthermore, dSINE patient also presented larger aortic taper ratio, stent oversizing ratio, and expansion mismatch ratio of the distal true lumen (EMRDTR). The larger mechanical parameters and growth rates in dSINE patient corresponded to a greater aortic taper ratio, stent oversizing ratio, and EMRDTR. The failure of dSINE prediction by the stent tortuosity index indicated that mechanical parameters were the fundamental reasons for dSINE development., (© 2024 John Wiley & Sons Ltd.)

The healthy adult aorta is a remarkably resilient structure, able to resist relentless cardiac-induced and hemodynamic loads under normal conditions. Fundamental to such mechanical homeostasis is the mechano-sensitive cell signaling that controls gene products and thus the structural integrity of the wall. Mouse models have shown that smooth muscle cell-specific disruption of transforming growth factor-beta (TGFβ) signaling during postnatal development compromises this resiliency, rendering the aortic wall susceptible to aneurysm and dissection under normal mechanical loading. By contrast, disruption of such signaling in the adult aorta appears to introduce a vulnerability that remains hidden under normal loading, but manifests under increased loading as experienced during hypertension. We present a multiscale (transcript to tissue) computational model to examine possible reasons for compromised mechanical homeostasis in the adult aorta following reduced TGFβ signaling in smooth muscle cells., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Ltd.)

The coronary sinus reducer (CSR) is an emerging medical device for treating patients with refractory angina, often associated with myocardial ischemia. Patients implanted with CSR have shown positive outcomes, but the underlying mechanisms are unclear. This study sought to understand the mechanisms of CSR by investigating its effects on coronary microcirculation hemodynamics that may help explain the therapy's efficacy. We applied a validated computer model of the coronary microcirculation to investigate how CSR affects hemodynamics under different degrees of coronary artery stenosis. With moderate coronary stenosis, an increase in capillary transit time (CTT) [up to 69% with near-complete coronary sinus (CS) occlusion] is the key change associated with CSR. Because capillaries in the microcirculation can still receive oxygenated blood from the upstream artery with moderate stenosis, the increase in CTT allows more time for the exchange of gases and nutrients, aiding tissue oxygenation. With severe coronary stenosis; however, the redistribution of blood draining from the nonischemic region to the ischemic region (up to 96% with near-complete CS occlusion) and the reduction in capillary flow heterogeneity are the key changes associated with CSR. Because blood draining from the nonischemic region is not completely devoid of O 2 , the redistribution of blood to the capillaries in the ischemic region by CSR is beneficial especially when little or no oxygenated blood reaches these capillaries. This simulation study provides insights into the mechanisms of CSR in improving clinical symptoms. The mechanisms differ with the severity of the upstream stenosis. NEW & NOTEWORTHY Emerging coronary venous retroperfusion treatments, particularly coronary sinus reducer (CSR) for refractory angina linked to myocardial ischemia, show promise; however, their mechanisms of action are not well understood. We find that CSR's effectiveness varies with the severity of coronary stenosis. In moderate stenosis, CSR improves tissue oxygenation by increasing capillary transit time, whereas in severe stenosis, it redistributes blood from nonischemic to ischemic regions and reduces capillary flow heterogeneity.

Dynamic cerebral autoregulation (dCA) is the mechanism that describes how the brain maintains cerebral blood flow approximately constant in response to short-term changes in arterial blood pressure. This is known to be impaired in many different pathological conditions, including ischaemic and haemorrhagic stroke, dementia and traumatic brain injury. Many different approaches have thus been used both to analyse and to quantify this mechanism in a range of healthy and diseased subjects, including data-driven models (in both the time and the frequency domain) and biophysical models. However, despite the substantial body of work on both biophysical models and data-driven models of dCA, there remains little work that links the two together. One of the reasons for this is proposed to be the discrepancies between the time constants that govern dCA in models and in experimental data. In this study, the processes that govern dCA are examined and it is proposed that the application of biophysical models remains limited due to a lack of understanding about the physical processes that are being modelled, partly due to the specific model formulation that has been most widely used (the equivalent electrical circuit). Based on the analysis presented here, it is proposed that the two most important time constants are arterial transit time and feedback time constant. It is therefore time to revisit equivalent electrical circuit models of dCA and to develop a more physiologically realistic alternative, one that can more easily be related to experimental data. KEY POINTS: Dynamic cerebral autoregulation is governed by two time constants. The first time constant is the arterial transit time, rather than the traditional 'RC' time constant widely used in previous models. This arterial transit time is approximately 1 s in the brain. The second time constant is the feedback time constant, which is less accurately known, although it is somewhat larger than the arterial transit time. The equivalent electrical circuit model of dynamic cerebral autoregulation should be replaced with a more physiologically representative model., (© 2024 The Authors. The Journal of Physiology © 2024 The Physiological Society.)

Purpose: In arterial disease, the presence of two or more serial stenotic lesions is common. For mild lesions, it is difficult to predict whether their combined effect is hemodynamically significant. This study assessed the hemodynamic significance of idealized serial stenotic lesions by simulating their hemodynamic interaction in a computational flow model., Materials and Methods: Flow was simulated with SimVascular software in 34 serial lesions, using moderate (15 mL/s) and high (30 mL/s) flow rates. Combinations of one concentric and two eccentric lesions, all 50% area reduction, were designed with variations in interstenotic distance and in relative direction of eccentricity. Fluid and fluid-structure simulations were performed to quantify the combined pressure gradient., Results: At a moderate flow rate, the combined pressure gradient of two lesions ranged from 3.8 to 7.7 mmHg, which increased to a range of 12.5-24.3 mmHg for a high flow rate. Eccentricity caused an up to two-fold increase in pressure gradient relative to concentric lesions. At a high flow rate, the combined pressure gradient for serial eccentric lesions often exceeded the sum of the individual lesions. The relative direction of eccentricity altered the pressure gradient by 15-25%. The impact of flow pulsatility and wall deformability was minor., Conclusion: This flow simulation study revealed that lesion eccentricity is an adverse factor in the hemodynamic significance of isolated stenotic lesions and in serial stenotic lesions. Two 50% lesions that are individually non-significant can combine more often than thought to hemodynamic significance in hyperemic conditions., (© 2024. The Author(s).)

Transposition of the great arteries (TGA) is a cyanotic congenital heart disease characterized by ventriculoarterial discordance and atrioventricular concordance with the great arteries in a parallel relationship. Prenatal diagnosis of TGA has implications for postnatal outcomes, allowing for planned delivery and perinatal management. Three-dimensional virtual or physical models of fetal TGA allow better understanding of fetal cardiac anomalies by parents and interactive discussion among the multidisciplinary team (obstetricians, pediatricians, maternal-fetal specialists, pediatric cardiologists, and cardiovascular surgeons), as well as continuing medical education., (© 2024. The Author(s), under exclusive licence to Springer Nature B.V.)

Mutations in cardiac myosin-binding protein C (cMyBP-C) or titin may respectively lead to hypertrophic (HCM) or dilated (DCM) cardiomyopathies. The mechanisms leading to these phenotypes remain unclear because of the challenge of translating cellular abnormalities to whole-heart and system function. We developed and validated a novel computer model of calcium-contraction coupling incorporating the role of cMyBP-C and titin based on the key assumptions: 1) tension in the thick filament promotes cross-bridge attachment mechanochemically, 2) with increasing titin tension, more myosin heads are unlocked for attachment, and 3) cMyBP-C suppresses cross-bridge attachment. Simulated stationary calcium-tension curves, isotonic and isometric contractions, and quick release agreed with experimental data. The model predicted that a loss of cMyBP-C function decreases the steepness of the calcium-tension curve, and that more compliant titin decreases the level of passive and active tension and its dependency on sarcomere length. Integrating this cellular model in the CircAdapt model of the human heart and circulation showed that a loss of cMyBP-C function resulted in HCM-like hemodynamics with higher left ventricular end-diastolic pressures and smaller volumes. More compliant titin led to higher diastolic pressures and ventricular dilation, suggesting DCM-like hemodynamics. The novel model of calcium-contraction coupling incorporates the role of cMyBP-C and titin. Its coupling to whole-heart mechanics translates changes in cellular calcium-contraction coupling to changes in cardiac pump and circulatory function and identifies potential mechanisms by which cMyBP-C and titin abnormalities may develop into HCM and DCM phenotypes. This modeling platform may help identify distinct mechanisms underlying clinical phenotypes in cardiac diseases., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)