No Arabic abstract
The human heart is enclosed in the pericardial cavity. The pericardium consists of a layered thin sac and is separated from the myocardium by a thin film of fluid. It provides a fixture in space and frictionless sliding of the myocardium. The influence of the pericardium is essential for predictive mechanical simulations of the heart. However, there is no consensus on physiologically correct and computationally tractable pericardial boundary conditions. Here we propose to model the pericardial influence as a parallel spring and dashpot acting in normal direction to the epicardium. Using a four-chamber geometry, we compare a model with pericardial boundary conditions to a model with fixated apex. The influence of pericardial stiffness is demonstrated in a parametric study. Comparing simulation results to measurements from cine magnetic resonance imaging reveals that adding pericardial boundary conditions yields a better approximation with respect to atrioventricular plane displacement, atrial filling, and overall spatial approximation error. We demonstrate that this simple model of pericardial-myocardial interaction can correctly predict the pumping mechanisms of the heart as previously assessed in clinical studies. Utilizing a pericardial model can not only provide much more realistic cardiac mechanics simulations but also allows new insights into pericardial-myocardial interaction which cannot be assessed in clinical measurements yet.
Patients resuscitated from cardiac arrest (CA) face a high risk of neurological disability and death, however pragmatic methods are lacking for accurate and reliable prognostication. The aim of this study was to build computational models to predict post-CA outcome by leveraging high-dimensional patient data available early after admission to the intensive care unit (ICU). We hypothesized that model performance could be enhanced by integrating physiological time series (PTS) data and by training machine learning (ML) classifiers. We compared three models integrating features extracted from the electronic health records (EHR) alone, features derived from PTS collected in the first 24hrs after ICU admission (PTS24), and models integrating PTS24 and EHR. Outcomes of interest were survival and neurological outcome at ICU discharge. Combined EHR-PTS24 models had higher discrimination (area under the receiver operating characteristic curve [AUC]) than models which used either EHR or PTS24 alone, for the prediction of survival (AUC 0.85, 0.80 and 0.68 respectively) and neurological outcome (0.87, 0.83 and 0.78). The best ML classifier achieved higher discrimination than the reference logistic regression model (APACHE III) for survival (AUC 0.85 vs 0.70) and neurological outcome prediction (AUC 0.87 vs 0.75). Feature analysis revealed previously unknown factors to be associated with post-CA recovery. Results attest to the effectiveness of ML models for post-CA predictive modeling and suggest that PTS recorded in very early phase after resuscitation encode short-term outcome probabilities.
Predictive high-fidelity finite element simulations of human cardiac mechanics co-mmon-ly require a large number of structural degrees of freedom. Additionally, these models are often coupled with lumped-parameter models of hemodynamics. High computational demands, however, slow down model calibration and therefore limit the use of cardiac simulations in clinical practice. As cardiac models rely on several patient-specific parameters, just one solution corresponding to one specific parameter set does not at all meet clinical demands. Moreover, while solving the nonlinear problem, 90% of the computation time is spent solving linear systems of equations. We propose a novel approach to reduce only the structural dimension of the monolithically coupled structure-windkessel system by projection onto a lower-dimensional subspace. We obtain a good approximation of the displacement field as well as of key scalar cardiac outputs even with very few reduced degrees of freedom while achieving considerable speedups. For subspace generation, we use proper orthogonal decomposition of displacement snapshots. To incorporate changes in the parameter set into our reduced order model, we provide a comparison of subspace interpolation methods. We further show how projection-based model order reduction can be easily integrated into a gradient-based optimization and demonstrate its performance in a real-world multivariate inverse analysis scenario. Using the presented projection-based model order reduction approach can significantly speed up model personalization and could be used for many-query tasks in a clinical setting.
An electrocardiogram (EKG) is a common, non-invasive test that measures the electrical activity of a patients heart. EKGs contain useful diagnostic information about patient health that may be absent from other electronic health record (EHR) data. As multi-dimensional waveforms, they could be modeled using generic machine learning tools, such as a linear factor model or a variational autoencoder. We take a different approach:~we specify a model that directly represents the underlying electrophysiology of the heart and the EKG measurement process. We apply our model to two datasets, including a sample of emergency department EKG reports with missing data. We show that our model can more accurately reconstruct missing data (measured by test reconstruction error) than a standard baseline when there is significant missing data. More broadly, this physiological representation of heart function may be useful in a variety of settings, including prediction, causal analysis, and discovery.
When they are damaged or injured, soft biological tissues are able to self-repair and heal. Mechanics is critical during the healing process, as the damaged extracellular matrix (ECM) tends to be replaced with a new undamaged ECM supporting homeostatic stresses. Computational modeling has been commonly used to simulate the healing process. However, there is a pressing need to have a unified thermodynamics theory for healing. From the viewpoint of continuum damage mechanics, some key parameters related to healing processes, for instance, the volume fraction of newly grown soft tissue and the growth deformation, can be regarded as internal variables and have related evolution equations. This paper is aiming to establish this unified framework inspired by thermodynamics for continuum damage models for the healing of soft biological tissues. The significant advantage of the proposed model is that no textit{ad hoc} equations are required for describing the healing process. Therefore, this new model is more concise and offers a universal approach to simulate the healing process. Three numerical examples are provided to demonstrate the effectiveness of the proposed model, which is in good agreement with the existing works, including an application for balloon angioplasty in an arteriosclerotic artery with a fiber cap.
Frequently during its lifetime a human organism is subjected to the acoustical and similar to them vibrating impacts. Under the certain conditions such influence may cause physiological changes in the organs functioning. Thus the study of the oscillatory mechanical impacts to the organism is very important task of the numerical physiology. It allows to investigate the endurance limits of the organism and to develop protective measures in order to extend them. The noise nuisances affects to the most parts of the organism disrupting their functions. The vibrating disturbances caused to the lung function as one of the most sensitive to the acoustical impacts is considered in this work. The model proposed to describe the air motion in trachea-bronchial tree is based on the one dimensional no-linear theory including mass and momentum conservation for the air flow in flexible tubes.