Propagation of uncertainty in the mechanical and biological response of growing tissues using multi-fidelity Gaussian process regression

Abstract

A key feature of living tissues is their capacity to remodel and grow in response to environmental cues. Within continuum mechanics, this process can be captured with the multiplicative split of the deformation gradient into growth and elastic contributions. The mechanical and biological response during tissue adaptation is characterized by inherent variability. Accounting for this uncertainty is critical to better understand tissue mechanobiology, and, moreover, it is of practical importance if we aim to develop predictive models for clinical use. However, the current gold standard in computational models of growth and remodeling remains the use of deterministic finite element (FE) simulations. Here we focus on tissue expansion, a popular technique in which skin is stretched by a balloon-like device inducing its growth. We construct FE models of tissue expansion with various levels of detail, and show that a sufficiently broad set of FE simulations from these models can be used to train an accurate and efficient multi-fidelity Gaussian process (GP) surrogate. The approach is not limited to simulation data, rather, it can fuse different kinds of data, including from experiments. The main appeal of the framework relies on the common experience that highly detailed models (or experiments) are more accurate but also more costly, while simpler models (or experiments) can be easily evaluated but are bound to have some error. In these situations, doing uncertainty analysis tasks with the high fidelity models alone is not feasible and, conversely, relying solely on low fidelity approximations is also undesirable. We show that a multi-fidelity GP outperforms the high fidelity GP and low fidelity GP when tested against the most detailed FE model. In turn, having trained the multi-fidelity GP model, we showcase the propagation of uncertainty from the mechanical and biological response parameters to the spatio-temporal growth outcomes. We expect that the methods and applications in this paper will enable future research in parameter calibration under uncertainty and uncertainty propagation in real clinical scenarios involving tissue growth and remodeling.Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

Multiscale Modeling Meets Machine Learning: What Can We Learn?

“…We can also use neural networks to directly map pressure and topology as input onto the deformation field as output, for example to accelerate decision making in electrosurgery [42]. Finally, natural questions that machine learning can help us answer focus on sensitivity analysis [62] and uncertainty quantification [61,146]. We have structured this introduction to the special issue around four methodological areas, ordinary and partial differential equations, and data and theory driven machine learning [3].…”

Gaussian Process Surrogates for Modeling Uncertainties in a Use Case of Forging Superalloys

“…where κ 1 = 5 3 mm 2 /kg and κ 2 = 5 3 mm 2 /kgs are constants. We utilize a designed forging press specific GP with data generated by using (15) to model the mean and input dependent deviations in respect to the manufacturing process characteristic Z (i.e., Z represents the velocity profile of the forging die v die ). Z is defined by a distribution with mean m(Z) and aleatoric standard deviation Σ al (Z).…”

“…We utilize (15) and different input parameter combinations to generate training data for the forging press GP, see Table 1. In terms of time step size t, we assume that each forging stroke lasts one second, and we model each stroke with a resolution of 100 time steps.…”

“…Venkatraman et al [14] use GP regression as a surrogate model of texture in micro-springs. GP regression can also be used on data with multiple fidelity levels, where Lee et al [15] investigate GP regression surrogate modeling with uncertain material properties of soft tissues and multi-fidelity data. Brevault et al [16] provide an overview of multi-fidelity GP regression techniques in the field of aerospace systems.…”

Sensitivity analysis of non‐local damage in soft biological tissues

“…The large computational cost is a main challenge in uncertainty analysis when a probabilistic method is used 19,20 . Although the Monte Carlo method is the most straightforward way to conduct probability analysis, 21 deterministic problems must usually be solved many times to guarantee the accuracy of the Monte Carlo method.…”