Heart failure is a long-term condition affecting 63 million people worldwide [33]. While heart trans-plantation is considered the gold standard for treating this kind of end-stage heart failure, the severes hortage of donor organs necessitates the exploration of alternative therapies such as total artificial hearts and ventricular assist devices. However, the current total artificial hearts are flawed; their rigid surfaces can induce thrombosis, haemolysis and infection [2][6][18][31]. A soft biomimetic approach could theoretically mitigate these flaws by providing a pumping mechanism that more closely replicates the physiology of the native ventricle [28][32]. Therefore, this thesis aims to manufacture and experimentally assess a soft biomimetic artificial ventricle.Osouli et al. developed a novel ventricular model, which describes the left ventricular myocardium as aseries of smoothly twisting surfaces [28], which can be used for this purpose. To mimic the ventricular sheet-like myocardial model, a flat pneumatic artificial muscle (FPAM) is used as an actuator. This typeof actuator is soft, strong and can contract under a low pressure of 50 kPa. Following systematic evaluation of materials and fabrication processes, nylon-coated TPU combined with heat pressing yields themost consistent results. Isotonic and isometric characterisation confirms state-of-the-art performance,with the FPAM achieving a mean contraction of 20.2 % and a peak force of 95.5 N.The FPAM forms the structural foundation of the artificial myocardium. The 3D model proposed by Osouli et al. [28] is transformed into a 2D surface to facilitate FPAM integration. Key design optimisations include the quantity and width of the zero-volume air chambers (ZACs). After which, this sheet is formed back into the myo architecture proposed by Osouli et al. Biomimetic motion evaluation revealsa ventricular twist of 9.03° (native: 17.0°), wall thickening of 12.0% (native: 37.5%), and longitudinal contraction of 1.37% (native: 11.9%). Functional assessment yields an ejection fraction of 38.4% anda maximum pressure of 82.9 mmHg.The prototype exhibits sub optimal performance relative to the native ventricle, a discrepancy primarily attributable to the limited contractile capacity of the artificial myocardium, which constrains both biomimetic motion and ventricular performance. Nevertheless, this thesis provides the first demonstration that a sheet-like actuator-based artificial ventricle can simultaneously induce bio mimicry andventricularlike performance, establishing this approach as a promising foundation for future devices. ... Read more

Post-myocardial infarction (MI) growth and remodeling (G&R), commonly referred to as fibrosis, involves both geometric deformation and progressive stiffening of infarcted tissue due to collagen accumulation. While zero-dimensional (0D) cardiac G&R models have successfully reproduced organ-level adaptations post-infarction, they often neglect evolving tissue properties associated with collagen turnover. In this study, we address this limitation by incorporating a time-dependent stiffening law into a strain-driven 0D framework, extending the original model by Witzenburg et al. (2018). Collagen turnover (CT) was modeled using a phenomenological exponential function, calibrated against experimental hydroxyproline data. The model was validated against independent canine datasets and benchmarked against both the original reference and a baseline No CT simulation. While full time-course verification was not achieved - due to inconsistencies in baseline reported parameters - control and acute states were accurately reproduced. Critically, the CT-enhanced model reduced the mean standardized z-score (MSZ) by 57.8%, with the most substantial improvements seen in ventricular volume and diastolic pressure predictions. These results confirm the added value of explicitly modeling tissue-level remodeling and highlight the importance of accurate initialization to ensure long-term prediction fidelity in reduced-order frameworks. ... Read more

Understanding how microstructural architecture governs macroscopic mechanical behavior is central to multiscale materials design, yet existing microstructure-informed workflows either rely on extensive experimental studies or costly microstructure-homogenization simulations. This thesis develops a unified, data-driven constitutive modeling framework that directly maps continuous two-phase microstructure to linear and nonlinear effective responses using statistical descriptors, bypassing reconstruction and generalizing across diverse morphologies. We first construct Micro3D, a statistically diverse synthetic dataset of binary microstructures using Gaussian random fields and multiple morphology generators, from which two-point statistics are extracted and compressed to serve as compact, physics-meaningful inputs. For the linear regime, a two-branch multilayer perceptron (MLP) is constructed with embedded symmetry and positive-definiteness constraints, using a Cholesky-based representation to predict the effective tensor. For the nonlinear regime, a hybrid framework combining a three-branch architecture, a hypernetwork, and an input-convex neural network (ICNN) is developed to capture complex material behaviors. Both models demonstrate strong generalization to unseen microstructures, with the nonlinear model accurately reproducing responses under previously unseen loading paths. Together, these components provide a practical route to microstructure-informed surrogate models that are interpretable, extensible, and suitable for downstream simulation. ... Read more

Engineered heart tissues (EHTs) provide a promising platform for modeling cardiac physiology, but their dense and heterogeneous fiber organization makes quantitative analysis highly challenging. This thesis presents an automated pipeline for fiber segmentation and structural analysis of confocal EHT images. The framework integrates frequency based preprocessing using FFT bandpass filtering, state of the art deep learning segmentation models (U-Net, Attention U-Net, and U-Net++), and post-processing refinement through a secondary U-Net. Evaluation was conducted on a synthetic labeled dataset and on real EHT slices with sparse annotations. The results highlight clear trade-offs between model architectures. U-Net produced the most complete and connected fibers but introduced substantial hallucinations. Attention U-Net generated clean outputs but with fragmented fibers, and U-Net++ achieved a balance by capturing directionality and coherence with reduced continuity. Refinement networks were effective at reducing thickness and noise in some cases, but they often removed true fibers and fragmented long structures, providing limited overall benefit. Fiber level metrics and human inspection confirmed these findings, showing that orientation is captured reliably across models, while continuity and connectivity remain major challenges. Overall, the pipeline demonstrates the feasibility of automated structural analysis of EHTs and establishes a foundation for future work with improved datasets, advanced refinement strategies, and broader use of pretrained or structurally informed models. ... Read more

Cardiovascular diseases (CVDs) are a group of disorders of the heart and blood vessels.
CVDs are the leading cause of death worldwide. To diagnose and treat CVDs, clinicians and cardiologists use multiple noninvasive imaging techniques. These scans are used to segment certain structures of the heart. Deep learning-based cardiac segmentation on short-axis cardiac magnetic resonance images (CMRI) has gained popularity over the past few years because of its generalisability and accuracy. This has exponentially reduced contouring times for clinicians. The development of such deep learning techniques has seen a common trend. In order to accommodate learning for larger cardiac datasets, the depth and size of segmentation networks have been increased. Unfortunately, the environmental impact of exploding such networks is not taken into account. One solution to mitigate having computationally expensive networks is to incorporate anatomical knowledge in the form of shape priors. The Gridnet and UNet with a shape prior are computationally efficient networks that are used to evaluate segmentation performance on a large and varied cardiac dataset (Combination of the Automated Cardiac Diagnosis Challenge - ACDC and Multi-Centre, Multi-Vendor and Multi-Disease Cardiac Segmentation challenge - M&M datasets). On average, these networks segment CMRIs with an average dice score of 0.87 and a Hausdorff distance of 11.7mm. In parallel, one of the major issues in cardiac technology is the under-representation of women in cardiac datasets. Purposefully curated cardiac datasets such as ACDC and M&M try and maintain equal representation. In real-world scenarios, this might not always be the case. Clinical trials to collect such data often report female representation as low as 25%. Evaluation of segmentation performance between a balanced and skewed dataset is conducted. This is to address if bias in such cardiac training datasets affects the performance of segmentation networks between male and female test patients.

... Read more

This thesis investigates the ability of Bayesian EUCLID to retrieve a predictive approximate material model for the myocardium in the presence of heterogeneous deformation fields due to simulated biaxial stretch tests. The Holzapfel-Ogden material model is used as the ground-truth material model of the simulation, since it is capable of describing the orthotropic and hyperelastic behavior of the myocardium. In total, four material models were retrieved from Bayesian-EUCLID by considering two fiber orientations within the sample and two displacement measurement techniques. The two fiber directions included one in which the sheet-like structure of the sample lies within the test plane and one in which it lies perpendicular to the test plane. The displacement field measurements considered were a full displacement field measurement and an approximation thereof based on a stereo digital image correlation method. The results show that when the sheet-like structure of the myocardium lies within the test plane of the biaxial stretch test and the full displacement field is available, Bayesian EUCLID is capable of finding orthotropic and hyperelastic material models that correspond well to the ground truth. This shows that a biaxial stretch test can adequately characterize the material behavior of the myocardium. ... Read more

A congenital heart defect (CHD) is an anomaly in the structure of the heart that is present at birth. In the last 15 years, a CHD is present in 9 per 1,000 live births, making it the most prevalent birth defect (Linde et al., 2011).
CHD’s prevalence coupled with its inherent complexity culminates into situations that are both complicated and unique. This makes the treatment of patients with CHDs highly patient-specific. Patient-specific physics-based computational models can aid in selecting the optimal course of treatment for the patient and serve as a predictive tool for different surgical plans (Capelli et al., 2011). The creation of patient-specific geometric heart models of children with CHDs will be the first crucial step towards a patient-specific approach in ameliorating this challenge. Approaches have to be found to face challenges including: a lack of high quality imaging data of the target population; slice misalignment due to breathing in a breath-hold procedure; and a lack of protocols in image segmentation.

To accomplish the objective, seven biventricular geometric models and meshes are created for healthy hearts, tetralogy of Fallot hearts and Fontan hearts. Slice misalignment is corrected using contours of the papillary muscle and the epicardium as a reference. Lastly, ground truth geometries are constructed by stacking disks that originate directly from the MRI segmentation.
To validate the geometric models, global and local approaches are used. The global validation includes ventricular mass and volume analyses with a relative error (RE) ranging from 0.00036 to 0.79. The Dice similarity coefficient (DSC) for the local validation is between 0.711 +/- 0.138 and 0.968 +/- 0.018. The slice misalignment correction is validated using qualitative and quantitative measures. A quantitative measure is the 4-chamber long-axis local validation which has a DSC of 0.877 for the epicardium and 0.916 for the papillary muscle correction method. Global validation in the ground truth exercise has a mean RE of 0.010 +/- 0.0066. Results of the local validation show a minimum DSC of 0.960 +/- 0.008.

Large differences in RE can be explained by the small number of MRI slices available, segmentation variability, the complex geometry of the right ventricle, and the methods of volume and mass computations. The DSC computed in the local validation of the models shows good to excellent agreement for all models. According to the long-axis local validation, the papillary muscle is the best reference method proposed for the slice misalignment correction. The ground truth exercise reveals close correlation to the Medis Suite measurements, possibly leading to underestimation of the computed volumes.

Patient-specific biventricular geometric models are created and validated in this thesis. The models are representative of a range of patients including healthy hearts and hearts with complex CHDs. A patient-specific model of the heart would be beneficial for clinicians to understand complex pathological cases, as the ones described in this work. Herein, a patient-specific biventricular geometric model is presented. By including a healthy population as well as patients affected by complex congenital heart defects, this work can be used as a starting point for the development of more patient-specific computational models of the heart.
... Read more