Effects of atrial contraction, atrioventricular interaction and heart valve dynamics on human cardiovascular system response

Hemodynamic effects of pulsatile frequency of right ventricular assist device (RVAD) on pulmonary perfusion: a simulation study

Right ventricular assist devices (RVADs) have been extensively used to provide hemodynamic support for patients with end-stage right heart (RV) failure. However, conventional in-parallel RVADs can lead to an elevation of pulmonary artery (PA) pressure, consequently increasing the right ventricular (RV) afterload, which is unfavorable for the relaxation of cardiac muscles and reduction of valve complications. The aim of this study is to investigate the hemodynamic effects of the pulsatile frequency of the RVAD on pulmonary artery. Firstly, a mathematical model incorporating heart, systemic circulation, pulmonary circulation, and RVAD is developed to simulate the cardiovascular system. Subsequently, the frequency characteristics of the pulmonary circulation system are analyzed, and the calculated results demonstrate that the pulsatile frequency of the RVAD has a substantive impact on the pulmonary artery pressure. Finally, to verify the analysis results, the hemodynamic effects of the pulsatile frequency of the RVAD on pulmonary artery are compared under diffident support modes. It is found that the pulmonary artery pressure decreases by approximately 6% when the pulsatile frequency changes from 1 to 3 Hz. The increased pulsatile frequency of RA-PA support mode may facilitate the opening of the pulmonary valve, while the RV-PA support mode can more effectively reduce the load of RV. This work provides a useful method to decrease the pulmonary artery pressure during the RVAD supports and may be beneficial for improving myocardial function in patients with end-stage right heart failure, especially those with pulmonary hypertension.

Differences in Mechanical, Electrical and Calcium Transient Performance of the Isolated Right Atrial and Ventricular Myocardium of Guinea Pigs at Different Preloads (Lengths)

There are only a few studies devoted to the comparative and simultaneous study of the mechanisms of the length-dependent regulation of atrial and ventricular contractility. Therefore, an isometric force-length protocol was applied to isolated guinea pig right atrial (RA) strips and ventricular (RV) trabeculae, with a simultaneous measurement of force (Frank-Starling mechanism) and Ca2+ transients (CaT) or transmembrane action potentials (AP). Over the entire length-range studied, the duration of isometric contraction, CaT and AP, were shorter in the RA myocardium than in the RV myocardium. The RA myocardium was stiffer than the RV myocardium. With the increasing length of the RA and RV myocardium, the amplitude and duration of isometric contraction and CaT increased, as well as the amplitude and area of the "CaT difference curves" (shown for the first time). However, the rates of the tension development and relaxation decreased. No contribution of AP duration to the heterometric regulation of isometric tension was found in either the RA or RV myocardium of the guinea pig. Changes in the degree of overlap of the contractile proteins of the guinea pig RA and RV myocardium mainly affect CaT kinetics but not AP duration.

Contractile Behavior of Right Atrial Myocardium of Healthy Rats and Rats with the Experimental Model of Pulmonary Hypertension

There is a lack of data about the contractile behavior of the right atrial myocardium in chronic pulmonary heart disease. We thoroughly characterized the contractility and Ca transient of isolated right atrial strips of healthy rats (CONT) and rats with the experimental model of monocrotaline-induced pulmonary hypertension (MCT) in steady state at different preloads (isometric force-length), during slow force response to stretch (SFR), and during post-rest potentiation after a period of absence of electrical stimulation (PRP). The preload-dependent changes in the isometric twitch and Ca transient did not differ between CONT and MCT rats while the kinetics of the twitch and Ca transient were noticeably slowed down in the MCT rats. The magnitude of SFR was significantly elevated in the MCT right atrial strips and this was accompanied by the significantly higher elevation of the Ca transient relative amplitude at the end of SFR. The slow changes in the contractility and Ca transient in the PRP protocol did not differ between CONT and MCT. In conclusion, the alterations in the contractility and Ca transient of the right atrial myocardium of monocrotaline-treated rats with pulmonary hypertension mostly concern the elevation in SFR. We hypothesize that this positive inotropic effect in the atrial myocardium may (partly) compensate the systolic deficiency of the right ventricular failing myocardium.

Estimation of left ventricular stroke work based on a large cohort of healthy children

Left ventricular stroke work is an important prognostic marker to analyze cardiac function. Standard values for children are, however, missing. For clinicians, standards can help to improve the treatment decision of heart failures. For engineers, they can help to optimize medical devices. In this study, we estimated the left ventricular stroke work for children based on modeled pressure-volume loops. A lumped parameter model was fitted to clinical data of 340 healthy children. Reference curves for standard values were created over age, weight, and height. Left ventricular volume was measured with 3D echocardiography, while maximal ventricular pressure was approximated with a regression model from the literature. For validation of this method, we used 18 measurements acquired by a conductance catheter in 11 patients. The method demonstrated a low absolute mean difference of 0.033 J (SD: 0.031 J) for stroke work between measurement and estimation, while the percentage error was 21.66 %. According to the resulting reference curves, left ventricular stroke work of newborns has a median of 0.06 J and increases to 1.15 J at the age of 18 years. Stroke work increases over weight and height in a similar trend. The percentile curves depict the distribution. We demonstrate how reference curves can be used for quantification of differences and comparison in patients.

Numerical Simulation of the Effect of Pulmonary Vascular Resistance on the Hemodynamics of Reoperation After Failure of One and a Half Ventricle Repair

OBJECTIVE

The one and a half ventricle repair (1.5VR) is a common clinical choice for patients with right heart dysfunction. Considering the influence of blood circulation failure and reoperation in urgent need, this essay aims to explore the hemodynamic effects of different pulmonary vascular resistance (PVR) values on reoperation after 1.5VR failure.

METHODS

The lumped parameter model (LPM) was used to simulate the reoperation, including the return biventricular repair (2VR), ligation of azygos vein (1.5VR') and return single ventricular repair (1.0VR). Firstly, the debugging parameters were used to simulate the hemodynamics of 2VR. Secondly, the value of PVR was changed from one to four times while the other parameters remained unchanged. Finally, 15 cardiac cycles were simulated and the 15th result was obtained. In this work, the left and right ventricular stroke work and their sum (Plv, Prv, Ptotal), the left and right ventricular ejection fraction (LVEF, RVEF), the mean Cardiac Output (mCO) and the mean pressure and flow-rate ratio of superior and inferior vena cava (mPsvc\mPivc and mQsvc\mQivc), respectively, were used to describe the hemodynamics of reoperation.

RESULTS

With the change of PVR from one to four times, the values of Plv, Prv, Ptotal, LVEF, and RVEF gradually decreased. The change rate of Plv, Ptotal and LVEF of 1.0VR were the largest in the three kinds of reoperation. The change rate of Prv of 1.5VR' was larger than that of 2VR, but it was the opposite for their EF change rate. The mCO of 2VR, 1.5VR', and 1.0VR decreased by 18.53%, 37.58%, and 48.07%, respectively. The mPsvc\mPivc of 1.5VR' increased from 3.76 to 6.77 and the mQsvc\mQivc decreased from 0.55 to 0.36, while the mPsvc\mPivc and mQsvc\mQivc of 2VR and 1.0VR remained 1 and 0.67, respectively. The peak value of the tricuspid flow-rate (Qti) waveform of 2VR and 1.5VR' changed from "E peak" to "A peak."

CONCLUSION

The numerical results demonstrate the highly reoperation-dependent hemodynamic consequences and their responses to variations in PVR. Comprehensive analysis of EF, mCO and ventricular stroke work indicates that PVR has a greater impact on 1.5VR' and 1.0VR. Therefore, we suggest that the selection strategy of reoperation should focus on PVR.

Mechano-electric finite element model of the left atrium

Mechanical stretch plays a major role in modulating atrial function, being responsible for beat-by-beat responses to changes in chamber preload, enabling a prompt regulation of cardiac function. Mechano-electric coupling (MEC) operates through many mechanisms and has many targets, making it experimentally difficult to isolate causes and effects especially under sinus conditions where effects are more transient and subtle. Therefore, modelling is a powerful tool to help understand the role of MEC with respect to the atrial electromechanical interaction. We propose a cellular-based computational model of the left atrium that includes a strongly coupled MEC component and mitral flow component to account for correct pressure generation in the atrial chamber as a consequence of blood volume and contraction. The method was applied to a healthy porcine left atrium. Results of the strongly coupled simulation show that strains are higher in the areas adjacent to the mitral annulus, the rim of the appendage, around the pulmonary venous trunks and at the location of the Bachmann's bundle, approximately between the mitral annulus and the region where the venous tissue transitions into atrial. These are regions where arrhythmias are likely to originate. The role of stretch-activated channels was very small for sinus rhythm for the single cardiac beat simulation, although tension development was very sensitive to stretch. The method could be applied to investigate potential therapeutic interventions acting on the mechano-electrical properties of the left atrium.

Numerical modeling of heart valves using resistive Eulerian surfaces

The goal of this work is the development and numerical implementation of a mathematical model describing the functioning of heart valves. To couple the pulsatile blood flow with a highly deformable thin structure (the valve's leaflets), a resistive Eulerian surfaces framework is adopted. A lumped-parameter model helps to couple the movement of the leaflets with the blood dynamics. A reduced circulation model describes the systemic hemodynamics and provides a physiological pressure profile at the downstream boundary of the valve. The resulting model is relatively simple to describe for a healthy valve and pathological heart valve functioning while featuring an affordable computational burden. Efficient time and spatial discretizations are considered and implemented. We address in detail the main features of the proposed method, and we report several numerical experiments for both two-dimensional and three-dimensional cases with the aim of illustrating its accuracy. Copyright © 2015 John Wiley & Sons, Ltd.

In-Series Versus In-Parallel Mechanical Circulatory Support for the Right Heart: A Simulation Study

Right heart failure (RHF) is a serious health issue with increasing incidence and high mortality. Right ventricular assist devices (RVADs) have been used to support the end-stage failing right ventricle (RV). Current RVADs operate in parallel with native RV, which alter blood flow pattern and increase RV afterload, associated with high tension in cardiac muscles and long-term valve complications. We are developing an in-series RVAD for better RV unloading. This article presents a mathematical model to compare the effects of RV unloading and hemodynamic restoration on an overloaded or failing RV. The model was used to simulate both in-series (sRVAD) and in-parallel (pRVAD) (right atrium-pulmonary artery cannulation) support for severe RHF. The results demonstrated that sRVAD more effectively unloads the RV and restores the balance between RV oxygen supply and demand in RHF patients. In comparison to simulated pRVAD and published clinical and in silico studies, the sRVAD was able to provide comparable restoration of key hemodynamic parameters and demonstrated superior afterload and volume reduction. This study concluded that in-series support was able to produce effective afterload reduction and preserve the valve functionality and native blood flow pattern, eliminating complications associated with in-parallel support.

In vitro comparison of two different mechanical circulatory support devices installed in series and in parallel

This study investigates the novel approach of placing a ventricular assist pump in the descending aorta in series configuration with the heart and compares it with the two traditional approaches of left-ventricle-to-ascending-aorta (LV-AA) and left-ventricle-to-descending-aorta (LV-DA) placement in parallel with the heart. Experiments were conducted by using the in-house simulator of the cardiovascular blood-flow loop (SCVL). The results indicate that the use of the LV-AA in-parallel configuration leads to a significant improvement in the systemic and pulmonic flow as the level of continuous flow is increased; however, this approach is considered highly invasive. The use of the LV-DA in-parallel configuration leads to an improvement in the systemic and pulmonic flow at lower levels of continuous flow but at higher levels of pump support leads to retrograde flow. In both in-parallel configurations, increasing the level of pump continuous flow leads to a decrease in pulsatility to a certain extent. The results of placing the pump in the descending aorta in series configuration show that the pressure drop upstream of the pump facilitates cardiac output as a result of afterload reduction. In addition, the pressure rise downstream of the pump may assist with renal perfusion. However, at the same time, the pressure drop generated at the proximal part of the descending aorta induces a slight drop in carotid perfusion, which would be autoregulated by the brain in a native cardiovascular system. The pulse wave analysis shows that placing the pump in the descending aorta leads to improved pulsatility in comparison with the traditional in-parallel configurations.

In vitro cardiovascular system emulator (bioreactor) for the simulation of normal and diseased conditions with and without mechanical circulatory support

This article presents a new device designed to simulate in vitro flow rates, pressures, and other parameters representing normal and diseased conditions of the human cardiovascular system. Such devices are sometimes called bioreactors or "mock" simulator of cardiovascular loops (SCVLs) in literature. Most SCVLs simulate the systemic circulation only and have inherent limitations in studying the interaction of left and right sides of circulation. Those SCVLs that include both left and right sides of the circulation utilize header reservoirs simulating cycles with constant atrial pressures. The SCVL described in this article includes models for all four chambers of the heart, and the systemic and pulmonary circulation loops. Each heart chamber is accurately activated by a separate linear motor to simulate the suction and ejection stages, thus capturing important features in the perfusion waveforms. Four mechanical heart valves corresponding to mitral, pulmonary, tricuspid, and aortic are used to control the desired unidirectional flow. This SCVL can emulate different physiological and pathological conditions of the human cardiovascular system by controlling the different parameters of blood circulation through the vascular tree (mainly the resistance, compliance, and elastance of the heart chambers). In this study, four cases were simulated: healthy, congestive heart failure, left ventricular diastolic dysfunction conditions, and left ventricular dysfunction with the addition of a mechanical circulatory support (MCS) device. Hemodynamic parameters including resistance, pressure, and flow have been investigated at aortic sinus, carotid artery, and pulmonary artery, respectively. The addition of an MCS device resulted in a significant reduction in mean blood pressure and re-establishment of cardiac output. In all cases, the experimental results are compared with human physiology and numerical simulations. The results show the capability of the SCVL to replicate various physiological and pathological conditions with and without MCS.

Review of zero-D and 1-D models of blood flow in the cardiovascular system

Background

Zero-dimensional (lumped parameter) and one dimensional models, based on simplified representations of the components of the cardiovascular system, can contribute strongly to our understanding of circulatory physiology. Zero-D models provide a concise way to evaluate the haemodynamic interactions among the cardiovascular organs, whilst one-D (distributed parameter) models add the facility to represent efficiently the effects of pulse wave transmission in the arterial network at greatly reduced computational expense compared to higher dimensional computational fluid dynamics studies. There is extensive literature on both types of models.

Method and Results

The purpose of this review article is to summarise published 0D and 1D models of the cardiovascular system, to explore their limitations and range of application, and to provide an indication of the physiological phenomena that can be included in these representations. The review on 0D models collects together in one place a description of the range of models that have been used to describe the various characteristics of cardiovascular response, together with the factors that influence it. Such models generally feature the major components of the system, such as the heart, the heart valves and the vasculature. The models are categorised in terms of the features of the system that they are able to represent, their complexity and range of application: representations of effects including pressure-dependent vessel properties, interaction between the heart chambers, neuro-regulation and auto-regulation are explored. The examination on 1D models covers various methods for the assembly, discretisation and solution of the governing equations, in conjunction with a report of the definition and treatment of boundary conditions. Increasingly, 0D and 1D models are used in multi-scale models, in which their primary role is to provide boundary conditions for sophisticate, and often patient-specific, 2D and 3D models, and this application is also addressed. As an example of 0D cardiovascular modelling, a small selection of simple models have been represented in the CellML mark-up language and uploaded to the CellML model repository http://models.cellml.org/. They are freely available to the research and education communities.

Conclusion

Each published cardiovascular model has merit for particular applications. This review categorises 0D and 1D models, highlights their advantages and disadvantages, and thus provides guidance on the selection of models to assist various cardiovascular modelling studies. It also identifies directions for further development, as well as current challenges in the wider use of these models including service to represent boundary conditions for local 3D models and translation to clinical application.

The role of the autonomic nervous system in hypertension: a bond graph model study

A bond graph model of the cardiovascular system with embedded autonomic nervous regulation was developed for a better understanding of the role of the autonomic nervous system (ANS) in hypertension. The model is described by a pump model of the heart and a detailed representation of the head and neck, pulmonary, coronary, abdomen and extremity circulation. It responds to sympathetic and parasympathetic activities by modifying systemic peripheral vascular resistance, heart rate, ventricular end-systolic elastance and venous unstressed volumes. The impairment of ANS is represented by an elevation of the baroreflex set point. The simulation results show that, compared with normotensive, in hypertension the systolic and diastolic blood pressure (SBP/DBP) rose from 112/77 mmHg to 144/94 mmHg and the left ventricular wall thickness (LVWT) increased from 10 mm to 12.74 mm. In the case that ANS regulation was absent, both the SBP and DBP further increased by 8 mmHg and the LVWT increased to 13.22 mm. The results also demonstrate that when ANS regulation is not severely damaged, e.g. the baroreflex set point is 97 mmHg, it still has an effect in preventing the rapid rise of blood pressure in hypertension; however, with the worsening of ANS regulation, its protective role weakens. The results agree with human physiological and pathological features in hemodynamic parameters and carotid baroreflex function curves, and indicate the role of ANS in blood pressure regulation and heart protection. In conclusion, the present model may provide a valid approach to study the pathophysiological conditions of the cardiovascular system and the mechanism of ANS regulation.