Search

Braided stents are self-expandable devices widely used in many different clinical applications. In-silico methods could be a useful tool to improve the design stage and preoperative planning; however, numerical modeling of braided structures is not trivial. The geometries are often challenging, and a parametric representation is not always easily achieved. Moreover, in the literature, different options have been proposed to handle the contact among the wires, but an extensive comparison of these modeling techniques is missing.In this work, both the geometry and contact issues are discussed. Firstly, an effective strategy based on parametric equations to draw complex braided geometries is illustrated and exploited to build three beam meshes resembling commercial devices. Secondly, three finite element simulations (bending, crimping and confined release) were carried out to compare simplified contact techniques involving connector elements with the more realistic but computationally expensive option based on the general contact algorithm, which has already been validated in the literature through comparisons with experimental results. Both local (stress distribution) and global quantities (forces/displacements) were analyzed.The results obtained using the connectors are significantly affected by wire interpenetrations and over-constraint.The percentage errors reached considerably high values, exceeding 100% in the confined release test and 50% in the remaining cases study. Moreover, the errors do not show uniform trends but vary according to the stent geometry, boundary conditions, connector type and investigated entity, suggesting that it is not possible to replace the use of the general contact algorithm with simplified approaches., 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 © 2020 Elsevier Ltd. All rights reserved.)

Initial palliation for univentricular hearts can be achieved via a systemic-to-pulmonary shunt (SPS). SPS configurations differ depending on the proximal anastomosis location, which might lead to dissimilar coronary and upper body perfusions. Mathematical modeling can be used to explore the local and global hemodynamic effects of the SPSs. In literature there are few patient-specific models of SPS that specifically address the influence of both the local and peripheral vasculature. In this study, multi-domain models of univentricular circulations were developed to investigate local hemodynamics and flow distribution in the presence of two shunt configurations. We also analyzed the relative impact of local and peripheral vascular resistances on coronary perfusion and flows through the upper aortic branches. A two-step approach was followed. First, two patient-specific models were based on clinical data collected from univentricular patients having different shunts and peripheral vasculatures. Each model coupled a three-dimensional representation of SPS, aortic arch (AA) and pulmonary arteries, with a lumped parameter model (LPM) of peripheral vasculature closing the circulatory loop. Then, two additional models of hypothetical subjects were created by coupling each customized LPM with the other patient's three-dimensional anatomy. Flow rates and pressures predicted by the patient-specific models revealed overall agreement with clinical data. Differences in the local hemodynamics were seen during diastole between the two models. Varying the three-dimensional models, while keeping an identical LPM, led to comparable flow distribution through the AA, suggesting that peripheral vasculatures have a dominant effect on local hemodynamics with respect to the shunt configuration., (Copyright © 2018 Elsevier Ltd. All rights reserved.)

Subclavian arteries are a possible alternate location for left ventricular assist device (LVAD) outflow grafts due to easier surgical access and application in high risk patients. As vascular blood flow mechanics strongly influence the clinical outcome, insights into the hemodynamics during LVAD support can be used to evaluate different grafting locations. In this study, the feasibility of left and right subclavian artery (SA) grafting was investigated for the HeartWare HVAD with a numerical multiscale model. A 3-D CFD model of the aortic arch was coupled to a lumped parameter model of the cardiovascular system under LVAD support. Grafts in the left and right SA were placed at three different anastomoses angles (90°, 60° and 30°). Additionally, standard grafting of the ascending and descending aorta was modelled. Full support LVAD (5l/min) and partial support LVAD (3l/min) in co-pulsation and counter-pulsation mode were analysed. The grafting positions were investigated regarding coronary and cerebral perfusion. Furthermore, the influence of the anastomosis angle on wall shear stress (WSS) was evaluated. Grafting of left or right subclavian arteries has similar hemodynamic performance in comparison to standard cannula positions. Angularity change of the graft anastomosis from 90° to 30° slightly increases the coronary and cerebral blood flow by 6-9% while significantly reduces the WSS by 35%. Cannulation of the SA is a feasible anastomosis location for the HVAD in the investigated vessel geometry., (Copyright © 2016 Elsevier Ltd. All rights reserved.)

A closed-loop lumped parameter model of blood circulation is considered for single-ventricle shunt physiology. Its parameters are estimated by an inverse problem based on patient-specific haemodynamics measurements. As opposed to a black-box approach, maximizing the number of parameters that are related to physically measurable quantities motivates the present model. Heart chambers are described by a single-fibre mechanics model, and valve function is modelled with smooth opening and closure. A model for valve prolapse leading to valve regurgitation is proposed. The method of data assimilation, in particular the unscented Kalman filter, is used to estimate the model parameters from time-varying clinical measurements. This method takes into account both the uncertainty in prior knowledge related to the parameters and the uncertainty associated with the clinical measurements. Two patient-specific cases - one without regurgitation and one with atrioventricular valve regurgitation - are presented. Pulmonary and systemic circulation parameters are successfully estimated, without assumptions on their relationships. Parameters governing the behaviour of heart chambers and valves are either fixed based on biomechanics, or estimated. Results of the inverse problem are validated qualitatively through clinical measurements or clinical estimates that were not included in the parameter estimation procedure. The model and the estimation method are shown to successfully capture patient-specific clinical observations, even with regurgitation, such as the double peaked nature of valvular flows and anomalies in electrocardiogram readings. Lastly, biomechanical implications of the results are discussed., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Coronary stent fracture is still an unresolved issue in the field of minimally invasive cardiovascular interventions due to its high rate of incidence and uncertain clinical consequences. Recent studies, based on clinical data, proved that there are several factors which can be identified as independently responsible of coronary stent fracture. Among these, calcifications, which increase the local stiffness and heterogeneity of atherosclerotic plaques, seem to play a major role. From a mechanical point of view, stent fracture in coronary arteries is triggered by the cyclic loading of pulsatile blood pressure combined with the movement of cardiac wall. In this context, this study aims at simulating the stent expansion in a model of epicardial atherosclerotic coronary artery and correlating the effects of cyclic blood pressure and cardiac wall movement on the stent fatigue resistance. Two ideal cases of atherosclerotic plaques were modelled: the first one included a localised plaque calcification; the latter one did not include such calcification. Results of stress/strain and fatigue analyses confirmed the influence of the plaque calcification on potential fracture of the devices. In addition, the effects of cardiac wall movement were quantified as more dangerous causes of the stent fatigue fracture with respect to the internal blood pressure oscillations. In conclusion, this study demonstrates the increased risk of coronary stent fracture associated to the presence of localised plaque calcifications. This work also suggests the necessity of more realistic biomechanical models which takes into account the heterogeneity of atherosclerotic plaques in order to assess the mechanical performances of coronary stents., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

Single ventricle hearts are congenital cardiovascular defects in which the heart has only one functional pumping chamber. The treatment for these conditions typically requires a three-staged operative process where Stage 1 is typically achieved by a shunt between the systemic and pulmonary arteries, and Stage 2 by connecting the superior venous return to the pulmonary circulation. Surgically, the Stage 2 circulation can be achieved through a procedure called the Hemi-Fontan, which reconstructs the right atrium and pulmonary artery to allow for an enlarged confluence with the superior vena cava. Based on pre-operative data obtained from two patients prior to Stage 2 surgery, we developed two patient-specific multi-scale computational models, each including the 3D geometrical model of the surgical junction constructed from magnetic resonance imaging, and a closed-loop systemic lumped-parameter network derived from clinical measurements. "Virtual" Hemi-Fontan surgery was performed on the 3D model with guidance from clinical surgeons, and a corresponding multi-scale simulation predicts the patient's post-operative hemodynamic and physiologic conditions. For each patient, a post-operative active scenario with an increase in the heart rate (HR) and a decrease in the pulmonary and systemic vascular resistance (PVR and SVR) was also performed. Results between the baseline and this "active" state were compared to evaluate the hemodynamic and physiologic implications of changing conditions. Simulation results revealed a characteristic swirling vortex in the Hemi-Fontan in both patients, with flow hugging the wall along the SVC to Hemi-Fontan confluence. One patient model had higher levels of swirling, recirculation, and flow stagnation. However, in both models, the power loss within the surgical junction was less than 13% of the total power loss in the pulmonary circulation, and less than 2% of the total ventricular power. This implies little impact of the surgical junction geometry on the SVC pressure, cardiac output, and other systemic parameters. In contrast, varying HR, PVR, and SVR led to significant changes in theses clinically relevant global parameters. Adopting a work-flow of customized virtual planning of the Hemi-Fontan procedure with patient-specific data, this study demonstrates the ability of multi-scale modeling to reproduce patient specific flow conditions under differing physiological states. Results demonstrate that the same operation performed in two different patients can lead to different hemodynamic characteristics, and that modeling can be used to uncover physiologic changes associated with different clinical conditions., (Copyright © 2012 Elsevier Ltd. All rights reserved.)

A less-invasive procedure that combines interventional stent placement in the ductus arteriosus and surgical banding of the branch pulmonary arteries has been recently introduced in the treatment of the hypoplastic left heart syndrome (HLHS). The hemodynamic behaviour of this hybrid approach has not been examined before in a mathematical model. In this study, a mathematical model of the hybrid procedure for HLHS is described, applying a multiscale approach that couples 3D models of the area of the surgical operation and lumped parameter models of the remaining circulation. The effects of various degrees of pulmonary banding and different stent sizes inserted in the ductus arteriosus on pulmonary-systemic flow ratio, cardiac output and oxygen delivery were assessed. Computational results suggest that balanced systemic and pulmonary blood flow and optimal systemic oxygen delivery are sensitive to the degree of pulmonary arterial banding and not to the size of the ductal stent., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Experimental tests, such as the confined and unconfined compression and the indentation tests, are traditionally used to determine the poroelastic properties of hydrated soft tissues (HSTs). The purpose of this study was to quantitatively evaluate the reliability of H(A) and K values as identified from experimental confined test data, estimating the errors that could occur in several situations with more realistic sample geometry and boundary conditions. Finite element models of the step-wise stress-relaxation confined compression tests on HSTs were developed including geometrical imperfections of the sample and the presence of a gap between the piston and the confining chamber. The errors occurring when H(A) and K were estimated by means of the analytical solution of the 1-D confined compression problem were assessed. Results of the analysis indicate that errors in the parameter estimation due to geometrical inaccuracies of the sample can be eliminated by applying a 5% strain pre-compression to the sample. Gap errors are negligible for H(A), can reach 20% for K, and cannot be eliminated by a pre-compression of the sample.

Multiscale computing is a challenging area even in biomechanics. Application of such a methodology to quantitatively compare postoperative hemodynamics in congenital heart diseases is very promising. In the treatment of hypoplastic left heart syndrome, which is a congenital heart disease where the left ventricle is missing or very small, the necessity to feed the pulmonary and systemic circulations is obtained with an interposition shunt. Two main options are available and differ from the sites of anastomoses: (i) the systemic-to-pulmonary conduit (Blalock-Taussig shunt known as the Norwood Operation (NO)) connecting the innominate artery (NO-BT) or the aorta (NO-CS) to the right pulmonary artery and (ii) the right ventricle to pulmonary artery shunt (known as Sano operation (SO)). The proposition that the SO is superior to the NO remains controversial. 3-D computer models of the NO (NO-BT and NO-CS) and SO were developed and investigated using the finite volume method. Conduits of 3, 3.5 and 4 mm were used in the NO models, whereas conduits of 4, 5 and 6 mm were used in the SO model. The hydraulic nets (lumped resistances, compliances, inertances and elastances) which represent the systemic, coronary and pulmonary circulations and the heart were identical in the two models. A multiscale approach was adopted to couple the 3-D models with the circulation net. Computer simulation results were compared with post-operative catheterization data. Results showed that (i) there is a good correlation between predicted and observed data: higher aortic diastolic pressure, decreased pulmonary arterial pressure, lower pulmonary-to-systemic flow ratio and higher coronary perfusion pressure in SO; (ii) there is a minimal regurgitant flow in the SO conduit. The close correlation between predicted and observed clinical data supports the use of mathematical modelling, with a mandatory multiscale approach, in the design and assessment of surgical procedures.

The objective of this study is to compare the coronary and pulmonary blood flow dynamics resulting from two configurations of systemic-to-pulmonary artery shunts currently utilized during the Norwood procedure: the central (CS) and modified Blalock Taussig (MBTS) shunts. A lumped parameter model of the neonatal cardiovascular circulation and detailed 3-D models of the shunt based on the finite volume method were constructed. Shunt sizes of 3, 3.5 and 4 mm were considered. A multiscale approach was adopted to prescribe appropriate and realistic boundary conditions for the 3-D models of the Norwood circulation. Results showed that the average shunt flow rate is higher for the CS option than for the MBTS and that pulmonary flow increases with shunt size for both options. Cardiac output is higher for the CS option for all shunt sizes. Flow distribution between the left and the right pulmonary arteries is not completely balanced, although for the CS option the discrepancy is low (50-51% of the pulmonary flow to the right lung) while for the MBTS it is more pronounced with larger shunt sizes (51-54% to the left lung). The CS option favors perfusion to the right lung while the MBTS favors the left. In the CS option, a smaller percentage of aortic flow is distributed to the coronary circulation, while that percentage rises for the MBTS. These findings may have important implications for coronary blood flow and ventricular function.