Data-Driven Simulation of Fisher-Kolmogorov Tumor Growth Models Using Dynamic Mode Decomposition.

The computer simulation of organ-scale biomechanistic models of cancer personalized via routinely collected clinical and imaging data enables to obtain patient-specific predictions of tumor growth and treatment response over the anatomy of the patient's affected organ. These patient-specific computational forecasts have been regarded as a promising approach to personalize the clinical management of cancer and derive optimal treatment plans for individual patients, which constitute timely and critical needs in clinical oncology. However, the computer simulation of the underlying spatiotemporal models can entail a prohibitive computational cost, which constitutes a barrier to the successful development of clinically-actionable computational technologies for personalized tumor forecasting. To address this issue, here we propose to utilize dynamic-mode decomposition (DMD) to construct a low-dimensional representation of cancer models and accelerate their simulation. DMD is an unsupervised machine learning method based on the singular value decomposition that has proven useful in many applications as both a predictive and a diagnostic tool. We show that DMD may be applied to Fisher-Kolmogorov models, which constitute an established formulation to represent untreated solid tumor growth that can further accommodate other relevant cancer phenomena (e.g., therapeutic effects, mechanical deformation). Our results show that a DMD implementation of this model over a clinically relevant parameter space can yield promising predictions, with short to medium-term errors remaining under 1% and long-term errors remaining under 20%, despite very short training periods. In particular, we have found that, for moderate to high tumor cell diffusivity and low to moderate tumor cell proliferation rate, DMD reconstructions provide accurate, bounded-error reconstructions for all tested training periods. Additionally, we also show that the three-dimensional DMD reconstruction of the tumor field can be leveraged to accurately reconstruct the displacement fields of the tumor-induced deformation of the host tissue. Thus, we posit the proposed data-driven approach has the potential to greatly reduce the computational overhead of personalized simulations of cancer models, thereby facilitating tumor forecasting, parameter identification, uncertainty quantification, and treatment optimization.

Modeling nonlocal behavior in epidemics via a reaction-diffusion system incorporating population movement along a network.

The outbreak of COVID-19, beginning in 2019 and continuing through the time of writing, has led to renewed interest in the mathematical modeling of infectious disease. Recent works have focused on partial differential equation (PDE) models, particularly reaction-diffusion models, able to describe the progression of an epidemic in both space and time. These studies have shown generally promising results in describing and predicting COVID-19 progression. However, people often travel long distances in short periods of time, leading to nonlocal transmission of the disease. Such contagion dynamics are not well-represented by diffusion alone. In contrast, ordinary differential equation (ODE) models may easily account for this behavior by considering disparate regions as nodes in a network, with the edges defining nonlocal transmission. In this work, we attempt to combine these modeling paradigms via the introduction of a network structure within a reaction-diffusion PDE system. This is achieved through the definition of a population-transfer operator, which couples disjoint and potentially distant geographic regions, facilitating nonlocal population movement between them. We provide analytical results demonstrating that this operator does not disrupt the physical consistency or mathematical well-posedness of the system, and verify these results through numerical experiments. We then use this technique to simulate the COVID-19 epidemic in the Brazilian region of Rio de Janeiro, showcasing its ability to capture important nonlocal behaviors, while maintaining the advantages of a reaction-diffusion model for describing local dynamics.

Dynamic mode decomposition in adaptive mesh refinement and coarsening simulations.

Dynamic mode decomposition (DMD) is a powerful data-driven method used to extract spatio-temporal coherent structures that dictate a given dynamical system. The method consists of stacking collected temporal snapshots into a matrix and mapping the nonlinear dynamics using a linear operator. The classical procedure considers that snapshots possess the same dimensionality for all the observable data. However, this often does not occur in numerical simulations with adaptive mesh refinement/coarsening schemes (AMR/C). This paper proposes a strategy to enable DMD to extract features from observations with different mesh topologies and dimensions, such as those found in AMR/C simulations. For this purpose, the adaptive snapshots are projected onto the same reference function space, enabling the use of snapshot-based methods such as DMD. The present strategy is applied to challenging AMR/C simulations: a continuous diffusion-reaction epidemiological model for COVID-19, a density-driven gravity current simulation, and a bubble rising problem. We also evaluate the DMD efficiency to reconstruct the dynamics and some relevant quantities of interest. In particular, for the SEIRD model and the bubble rising problem, we evaluate DMD's ability to extrapolate in time (short-time future estimates).

Assessing the Spatio-temporal Spread of COVID-19 via Compartmental Models with Diffusion in Italy, USA, and Brazil.

The outbreak of COVID-19 in 2020 has led to a surge in interest in the mathematical modeling of infectious diseases. Such models are usually defined as compartmental models, in which the population under study is divided into compartments based on qualitative characteristics, with different assumptions about the nature and rate of transfer across compartments. Though most commonly formulated as ordinary differential equation models, in which the compartments depend only on time, recent works have also focused on partial differential equation (PDE) models, incorporating the variation of an epidemic in space. Such research on PDE models within a Susceptible, Infected, Exposed, Recovered, and Deceased framework has led to promising results in reproducing COVID-19 contagion dynamics. In this paper, we assess the robustness of this modeling framework by considering different geometries over more extended periods than in other similar studies. We first validate our code by reproducing previously shown results for Lombardy, Italy. We then focus on the U.S. state of Georgia and on the Brazilian state of Rio de Janeiro, one of the most impacted areas in the world. Our results show good agreement with real-world epidemiological data in both time and space for all regions across major areas and across three different continents, suggesting that the modeling approach is both valid and robust.