The effect of pathological shoulder rhythm on muscle and joint forces after reverse shoulder arthroplasty, a numerical analysis

BACKGROUND

Compromised abduction ability after reverse shoulder arthroplasty is primarily linked to limited glenohumeral range of motion while scapulothoracic mobility can typically be maintained. Glenohumeral joint forces strongly depend on the resulting scapulohumeral rhythm, however, an association between the acting muscle and joint forces and the subject-specific scapulohumeral rhythm after reverse shoulder arthroplasty has not been established.

METHODS

Eleven reverse shoulder arthroplasty patients were divided into groups of poor and excellent abduction ability. Subject-specific models were developed and scaled for each patient using existing motion capture data in AnyBody™. Shoulder muscle and joint forces were obtained using inverse dynamics calculations during shoulder abduction to 100° in the scapula plane. The scapulohumeral rhythm, the resting abduction angle and internal body forces between the outcome groups were compared using a Mann Whitney U test.

FINDINGS

The mean glenohumeral and scapulothoracic contribution to overall shoulder abduction for the excellent group was on average 9.7% higher and 21.4% lower, respectively, compared to the mean of the poor group. For shoulder abduction angles between 30° and 60°, the excellent group demonstrated on average 25% higher muscle forces in the anterior deltoid which was significantly higher compared to the poor outcome patients. Scapulothoracic muscle activity did not differ significantly between the two functional groups.

INTERPRETATION

Accordingly, rehabilitation strategies focusing on strengthening the anterior part of the deltoid in particular may improve clinical outcomes.

Hemodynamics and bio-mechanics of morphologically distinct saccular intracranial aneurysms at bifurcations: Idealised vs Patient-specific geometries

Understanding the factors that influence the rupture of aneurysms is of primary concern to the clinicians, who are grappled with patient management. It is important to know how the relation between morphological features of the cerebral aneurysm, and the mechanical stresses on the containing arterial walls are influenced by the hemodynamic forces. Present study investigates three different shapes, which have been identified correspondingly in patient-specific scenarios as well. The primary objective is to categorize the bifurcation aneurysms into standard shapes such as, spherical, beehive and pear-shaped, based on patient-specific clinical studies and further compare and contrast the model aneurysms with the patient specific configurations, for their hemodynamic factors as well as the attendant stresses on the wall. MethodsComputational fluid dynamic simulations are performed accounting for the fluid-structure interaction (FSI) effects between the flowing fluid and the containing vessel wall. Blood is assumed to be Newtonian, while the arterial walls are assumed to be linearly elastic. A commercial solver is used for performing detailed calculations. Hemodynamic and bio-mechanical rupture predictions are carried out for the three different shapes. Observations derived from the idealised simulations are compared and contrasted against their patient-specific counterparts. ResultsFrom detailed numerical simulations, it was observed that pear-shaped aneurysms exhibit large re-circulation bubble and flow stagnation zone, with higher residence time for the particles, which may lead to atherosclerotic lesions. Beehive shape allows for maximum flow into the aneurysmal sac with concentrated jet impinging on the dome, leading to high values of maximum WSS (MWSS) resulting in great propensity to form a secondary bleb. However, flow field inside a spherical aneurysm is found to be stable with fewer vortices, and nearly uniform distribution of wall stresses are observed though-out the sac, which perhaps signifies hemodynamically and bio-mechanically stable condition. ConclusionCategorizing patient-specific intracranial aneurysms into standard shapes viz, spherical, beehive and pear could generalize the process of prediction of hemodynamic and bio-mechanical rupture indicators. Comparative assessment of the flow field and stresses reported from the simulations on idealised models, with corresponding patient-specific simulations reveal that, these studies could aid in understanding the generalised shape dependence of hemodynamic and bio-mechanical behaviour of aneurysms.

MRI in CFD for chronic type B aortic dissection: Ready for prime time?

OBJECTIVES

Better tools are needed for risk assessment of Type B aortic dissection (TBAD) to determine optimal treatment for patients with uncomplicated disease. Magnetic resonance imaging (MRI) has the potential to inform computational fluid dynamics (CFD) simulations for TBAD by providing individualised quantification of haemodynamic parameters, for assessment of complication risks. This systematic review aims to present an overview of MRI applications for CFD studies of TBAD.

METHODS

Following PRISMA guidelines, a search in Medline, Embase, and the Scopus Library identified 136 potentially relevant articles. Studies were included if they used MRI to inform CFD simulation in TBAD.

RESULTS

There were 20 articles meeting the inclusion criteria. 19 studies used phase contrast MRI (PC-MRI) to provide data for CFD flow boundary conditions. In 12 studies, CFD haemodynamic parameter results were validated against PC-MRI. In eight studies, geometric models were developed from MR angiography. In three studies, aortic wall or intimal flap motion data were derived from PC/cine MRI.

CONCLUSIONS

MRI provides complementary patient-specific information in CFD haemodynamic studies for TBAD that can be used for personalised care. MRI provides structural, dynamic and flow data to inform CFD for pre-treatment planning, potentially advancing its integration into clinical decision-making. The use of MRI to inform CFD in TBAD surgical planning is promising, however further validation and larger cohort studies are required.

An image-based computational hemodynamics study of the Systolic Anterior Motion of the mitral valve

Systolic Anterior Motion (SAM) of the mitral valve - often associated with Hypertrophic Obstructive Cardiomyopathy (HOCM) - is a cardiac pathology in which a functional subaortic stenosis is induced during systole by the mitral leaflets partially obstructing the outflow tract of the left ventricle. Its assessment by diagnostic tests is often difficult, possibly underestimating its severity and thus increasing the risk of heart failure. In this paper, we propose a new computational pipeline, based on cardiac cine Magnetic Resonance Imaging (cine-MRI) data, for the assessment of SAM. The pipeline encompasses image processing of the left ventricle and the mitral valve, and numerical investigation of cardiac hemodynamics by means of Computational Fluid Dynamics (CFD) in a moving domain with image-based prescribed displacement. Patient-specific geometry and motion of the left ventricle are considered in view of an Arbitrary Lagrangian-Eulerian approach for CFD, while the reconstructed mitral valve is immersed in the computational domain by means of a resistive method. We assess clinically relevant flow and pressure indicators in a parametric study for different degrees of SAM severity, in order to provide a better quantitative evaluation of the pathological condition. Moreover, we provide specific indications for its possible surgical treatment, i.e. septal myectomy.

A rule-based method for predicting the electrical activation of the heart with cardiac resynchronization therapy from non-invasive clinical data

Background

Cardiac Resynchronization Therapy (CRT) is one of the few effective treatments for heart failure patients with ventricular dyssynchrony. The pacing location of the left ventricle is indicated as a determinant of CRT outcome.

Objective

Patient specific computational models allow the activation pattern following CRT implant to be predicted and this may be used to optimize CRT lead placement.

Methods

In this study, the effects of heterogeneous cardiac substrate (scar, fast endocardial conduction, slow septal conduction, functional block) on accurately predicting the electrical activation of the LV epicardium were tested to determine the minimal detail required to create a rule based model of cardiac electrophysiology. Non-invasive clinical data (CT or CMR images and 12 lead ECG) from eighteen patients from two centers were used to investigate the models.

Results

Validation with invasive electro-anatomical mapping data identified that computer models with fast endocardial conduction were able to predict the electrical activation with a mean distance errors of 9.2 ± 0.5 mm (CMR data) or (CT data) 7.5 ± 0.7 mm.

Conclusion

This study identified a simple rule-based fast endocardial conduction model, built using non-invasive clinical data that can be used to rapidly and robustly predict the electrical activation of the heart. Pre-procedural prediction of the latest electrically activating region to identify the optimal LV pacing site could potentially be a useful clinical planning tool for CRT procedures.

The effect of chordae tendineae on systolic flow

When using Computational Fluid Dynamics to simulate ventricular blood flow in the heart, it has been common practice to neglect the effect of the sub-valvular apparatus and the trabeculae on the flow conditions. In this study, we analyze the effect of neglecting the chordae tendineae on the fluid flow and pressure drop. To test the assumption we use a previously developed dynamic 3D model of the left ventricle, aorta and valves that is based on 3D echocardiographic recordings. To this model we add the chordae tendineae as a sub-grid model. The previously developed 3D model for the left ventricle during systole is based on real-time three-dimensional echocardiography (RT3DE) recordings of a 30 years old female volunteer. The segmented ventricular wall does not include details of the aorta and the mitral valve, so these were reconstructed. The subgrid model for the flow across the chordae tendineae is based on the Actuator Line Method, which means that they are represented by drag coefficients. The analysis shows that the effect of the chordae tendineae on the pressure drop and work efficiency of the normal heart during systole is minor, and it seems that for simulating ventricular fluid flow and pressure drop during systole, one can follow the current practice and ignore the chordae. However, there can be local effects such as small vortices behind the chordae. Whether such effects are important for a particular application must be evaluated for the given case.

Prediction of stenting related adverse events through patient-specific finite element modelling

Right ventricular outflow tract (RVOT) calcific obstruction is frequent after homograft conduit implantation to treat congenital heart disease. Stenting and percutaneous pulmonary valve implantation (PPVI) can relieve the obstruction and prolong the conduit lifespan, but require accurate pre-procedural evaluation to minimize the risk of coronary artery (CA) compression, stent fracture, conduit injury or arterial distortion. Herein, we test patient-specific finite element (FE) modeling as a tool to assess stenting feasibility and investigate clinically relevant risks associated to the percutaneous intervention. Three patients undergoing attempted PPVI due to calcific RVOT conduit failure were enrolled; the calcific RVOT, the aortic root and the proximal CA were segmented on CT scans for each patient. We numerically reproduced RVOT balloon angioplasty to test procedure feasibility and the subsequent RVOT pre-stenting expanding the stent through a balloon-in-balloon delivery system. Our FE framework predicted the occurrence of CA compression in the patient excluded from the real procedure. In the two patients undergoing RVOT stenting, numerical results were consistent with intraprocedural in-vivo fluoroscopic evidences. Furthermore, it quantified the stresses on the stent and on the relevant native structures, highlighting their marked dependence on the extent, shape and location of the calcific deposits. Stent deployment induced displacement and mechanical loading of the calcific deposits, also impacting on the adjacent anatomical structures. This novel workflow has the potential to tackle the analysis of complex RVOT clinical scenarios, pinpointing the procedure impact on the dysfunctional anatomy and elucidating potential periprocedural complications.

An optimization algorithm for individualized biomechanical analysis and simulation of tibia fractures

An algorithmic strategy to determine the minimal fusion area of a tibia pseudarthrosis to achieve mechanical stability is presented. For this purpose, a workflow capable for implementation into clinical routine workup of tibia pseudarthrosis was developed using visual computing algorithms for image segmentation, that is a coarsening protocol to reduce computational effort resulting in an individualized volume-mesh based on computed tomography data. An algorithm detecting the minimal amount of fracture union necessary to allow physiological loading without subjecting the implant to stresses and strains that might result in implant failure is developed. The feasibility of the algorithm in terms of computational effort is demonstrated. Numerical finite element simulations show that the minimal fusion area of a tibia pseudarthrosis can be less than 90% of the full circumferential area given a defined maximal von Mises stress in the implant of 80% of the total stress arising in a complete pseudarthrosis of the tibia.