While a large and growing literature exists on mathematical and computational models of tumor growth, to date tumor growth models are largely qualitative in nature, and fall far short of being able to provide predictive results important in life-and-death decisions. This is largely due to the enormous complexity of evolving biological and chemical processes in living tissue and the complex interactions of many cellular and vascular constituents in living organisms. Several new technologies have emerged, however, which could lead to significant progress in this important area: (i) the development of so-called phase-field, or diffuse-interface models, which can be developed using continuum mixture theory, and which provide a general framework for modeling the action of multiple interacting constituents. These are based on generalizations of the Cahn–Hilliard models for spinodal decomposition, and have been used recently in certain tumor growth theories; (ii) the emergence of predictive computational methods based on the use of statistical methods for calibration, model validation, and uncertainty quantification; (iii) advances in imaging, experimental cell biology, and other medical observational methodologies; and (iv) the advent of petascale computing that makes possible the resolution of features at scales and at speeds that were unattainable only a short time in the past. Here we develop a general phenomenological thermomechanical theory of mixtures that employs phase-field or diffuse interface models of surface energies and reactions and which provides a framework for generalizing existing theories of the types that are in use in tumor growth modeling. In principle, the framework provides for the effects of M solid constituents, which may undergo large deformations, and for the effect of N - M fluid constituents, which could include highly nonlinear, non-Newtonian fluids. We then describe several special cases which have the potential of providing acceptable models of tumor growth. We then describe the beginning steps of the development of Bayesian methods for statistical calibration, model validation, and uncertainty quantification, which, with further work, could produce a truly predictive tool for studying tumor growth. In particular, we outline the processes of statistical calibration and validation that can be employed to determine if tumor growth models, drawn from the broad class of models developed here, are valid for prediction of key quantities of interest critical to making decisions on various medical protocols. We also describe how uncertainties in such key quantities of interest can be quantified in ways that can be used to establish confidence in predicted outcomes.

Laser surgery, or laser-induced thermal therapy, is a minimally invasive alternative or adjuvant to surgical resection in treating tumors embedded in vital organs with poorly defined boundaries. Its use, however, is limited due to the lack of precise control of heating and slow rate of thermal diffusion in the tissue. Nanoparticles, such as nanoshells, can act as intense heat absorbers when they are injected into tumors. These nanoshells can enhance thermal energy deposition into target regions to improve the ability for destroying larger cancerous tissue volumes with lower thermal doses. The goal of this paper is to present an integrated computer model using a so-called nestedblock optimization algorithm to simulate laser surgery and provide transient temperature field predictions. In particular, this algorithm aims to capture changes in optical and thermal properties due to nanoshell inclusion and tissue property variation during laser surgery. Numerical results show that this model is able to characterize variation of tissue properties for laser surgical procedures and predict transient temperature fields comparable to those measured by in vivo NIH Public Access

Predicting the outcome of thermotherapies in cancer treatment requires an accurate characterization of the bioheat transfer processes in soft tissues. Due to the biological and structural complexity of tumor (soft tissue) composition and vasculature, it is often very difficult to obtain reliable tissue properties that is one of the key factors for the accurate treatment outcome prediction. Efficient algorithms employing in vivo thermal measurements to determine heterogeneous thermal tissues properties in conjunction with a detailed sensitivity analysis can produce essential information for model development and optimal control. The goals of this paper are to present a general formulation of the bioheat transfer equation for heterogeneous soft tissues, review models and algorithms developed for cell damage, heat shock proteins, and soft tissues with nanoparticle inclusion, and demonstrate an overall computational strategy for developing a laser treatment framework with the ability to perform real-time robust calibrations and optimal control. This computational strategy can be applied to other thermotherapies using the heat source such as radio frequency or high intensity focused ultrasound.

For cancerous tumors in vital internal organs, minimally invasive laser surgery may be a desirable choice for cancer treatment due to its precise control and compatibility with most of the imaging modalities such as MRI (magnetic resonance imaging). However, the complexity of tumor composition and tissue response to a thermal dose demands real time optimization and control. In the previous work, we have developed a quite general computational framework that is capable of processing MRI anatomical data, providing pretreatment surgical protocol, and controlling tissue damage based on in vivo MRTI (magnetic resonance thermal imaging) data. In this paper, we describe computational techniques that are involved in real time optimization and control for laser surgical protocols of cancer treatment.