Quantification of total uncertainty in the physics-informed reconstruction of CVSim-6 physiology

When predicting physical phenomena through simulation, quantification of the total uncertainty due to multiple sources is as crucial as making sure the underlying numerical model is accurate. Possible sources include irreducible aleatoric uncertainty due to noise in the data, epistemic uncertainty induced by insufficient data or inadequate parameterization and model-form uncertainty related to the use of misspecified model equations. In addition, recently proposed approaches provide flexible ways to combine information from data with full or partial satisfaction of equations that typically encode physical principles. Physics-based regularization interacts in non-trivial ways with aleatoric, epistemic and model-form uncertainty and their combination, and a better understanding of this interaction is needed to improve the predictive performance of physics-informed digital twins that operate under real conditions. To better understand this interaction, with a specific focus on biological and physiological models, this study investigates the decomposition of total uncertainty in the estimation of states and parameters of a differential system simulated with MC X-TFC, a new physics-informed approach for uncertainty quantification based on random projections and Monte Carlo sampling. After an introductory comparison between approaches for physics-informed estimation, MC X-TFC is applied to a six-compartment stiff ODE system, the CVSim-6 model, developed in the context of human physiology. The system is first analysed by progressively removing data while estimating an increasing number of parameters, and subsequently by investigating total uncertainty under model-form misspecification of nonlinear resistance in the pulmonary compartment. In particular, we focus on the interaction between the formulation of the discrepancy term and quantification of model-form uncertainty, and show how additional physics can help in the estimation process. The method demonstrates robustness and efficiency in estimating unknown states and parameters, even with limited, sparse and noisy data. It also offers great flexibility in integrating data with physics for improved estimation, even in cases of model misspecification.This article is part of the theme issue 'Uncertainty quantification for healthcare and biological systems (Part 1)'.

Computational modelling of valvular heart disease: haemodynamic insights and clinical implications

An aging population and an increasing incidence of cardiovascular risk factors form the basis for a global rising prevalence of valvular heart disease (VHD). Research to further our understanding of the pathophysiology of VHD is often confined to the clinical setting. However, in recent years, sophisticated computational models of the cardiovascular system have been increasingly used to investigate a variety of VHD states. Computational modelling provides new opportunities to gain insights into pathophysiological processes that may otherwise be difficult, or even impossible, to attain in human or animal studies. Simulations of co-existing cardiac pathologies, such as heart failure, atrial fibrillation, and mixed valvular disease, have unveiled new insights that can inform clinical research and practice. More recently, advancements have been made in using models for making patient-specific diagnostic predictions. This review showcases valuable insights gained from computational studies on VHD and their clinical implications.

An In Silico Model for Predicting the Efficacy of Edge-to-Edge Repair for Mitral Regurgitation

In recent years, transcatheter edge-to-edge repair (TEER) has been widely adopted as an effective treatment for mitral regurgitation (MR). The aim of this study is to develop a personalized in silico model to predict the effect of edge-to-edge repair in advance to the procedure for each individual patient. For this purpose, we propose a combination of a valve deformation model for computing the mitral valve (MV) orifice area (MVOA) and a lumped parameter model for the hemodynamics, specifically mitral regurgitation volume (RVol). Although we cannot obtain detailed information on the three-dimensional flow field near the mitral valve, we can rapidly simulate the important medical parameters for the clinical decision support. In the present method, we construct the patient-specific pre-operative models by using the parameter optimization and then simulate the postoperative state by applying the additional clipping condition. The computed preclip MVOAs show good agreement with the clinical measurements, and the correlation coefficient takes 0.998. In addition, the MR grade in terms of RVol also has good correlation with the grade by ground truth MVOA. Finally, we try to investigate the applicability for the predicting the postclip state. The simulated valve shapes clearly show the well-known double orifice and the improvement of the MVOA, compared with the preclip state. Similarly, we confirmed the improved reverse flow and MR grade in terms of RVol. A total computational time is approximately 8 h by using general-purpose PC. These results obviously indicate that the present in silico model has good capability for the assessment of edge-to-edge repair.

Effects of Choice of Medical Imaging Modalities on a Non-invasive Diagnostic and Monitoring Computational Framework for Patients With Complex Valvular, Vascular, and Ventricular Diseases Who Undergo Transcatheter Aortic Valve Replacement

Due to the high individual differences in the anatomy and pathophysiology of patients, planning individualized treatment requires patient-specific diagnosis. Indeed, hemodynamic quantification can be immensely valuable for accurate diagnosis, however, we still lack precise diagnostic methods for numerous cardiovascular diseases including complex (and mixed) valvular, vascular, and ventricular interactions (C3VI) which is a complicated situation made even more challenging in the face of other cardiovascular pathologies. Transcatheter aortic valve replacement (TAVR) is a new less invasive intervention and is a growing alternative for patients with aortic stenosis. In a recent paper, we developed a non-invasive and Doppler-based diagnostic and monitoring computational mechanics framework for C3VI, called C3VI-DE that uses input parameters measured reliably using Doppler echocardiography. In the present work, we have developed another computational-mechanics framework for C3VI (called C3VI-CT). C3VI-CT uses the same lumped-parameter model core as C3VI-DE but its input parameters are measured using computed tomography and a sphygmomanometer. Both frameworks can quantify: (1) global hemodynamics (metrics of cardiac function); (2) local hemodynamics (metrics of circulatory function). We compared accuracy of the results obtained using C3VI-DE and C3VI-CT against catheterization data (gold standard) using a C3VI dataset (N = 49) for patients with C3VI who undergo TAVR in both pre and post-TAVR with a high variability. Because of the dataset variability and the broad range of diseases that it covers, it enables determining which framework can yield the most accurate results. In contrast with C3VI-CT, C3VI-DE tracks both the cardiac and vascular status and is in great agreement with cardiac catheter data.

Varying speed modulation of continuous-flow left ventricular assist device based on cardiovascular coupling numerical model

Continuous-flow left ventricular assist devices (CFLVADs) routinely operate at a constant speed for the support of a failing heart, which decreases the pulsatility in the arteries. Some late complications could be related to a long-term lack of pulsatility. Modulating the CFLVAD speed is a solution to enhance the pulsatility. The purpose of this study is to modulate multiple varying speed patterns and investigate their effects on the ventricle and vascular system. A cardiovascular coupling numerical model is developed to provide a simulation platform for testing the varying speed patterns. The varying speed patterns are modulated by combining the shape, amplitude, frequency, phase shift, and pulsatile duty cycle of the speed profile. The influence of varying speed support is examined by analyzing the indexes of pulsatility, indexes of ventricular unloading, and hemodynamic variables. The results show that the synchronous counterpulsation pattern can effectively reduce the ventricular unloading indexes, whereas the low-frequency asynchronous pattern can effectively increase the vascular pulsatility indexes. Also, the hemodynamics with synchronous varying speed support is more physiological than that with asynchronous varying speed support. This study provides valuable insight for further optimization of varying speed modulation by weighing vascular pulsatility, ventricular unloading, and hemodynamics.