One of the fundamental problems in simulating the motion of sharp interfaces between immiscible fluids is a description of the transition that occurs when the interfaces merge and reconnect. It is well known that classical methods involving sharp interfaces fail to describe this type of phenomena. Following some previous work in this area, we suggest a physically motivated regularization of the Euler equations which allows topological transitions to occur smoothly. In this model, the sharp interface is replaced by a narrow transition layer across which the fluids may mix. The model describes a flow of a binary mixture, and the internal structure of the interface is determined by both diffusion and motion. An advantage of our regularization is that it automatically yields a continuous description of surface tension, which can play an important role in topological transitions. An additional scalar field is introduced to describe the concentration of one of the fluid components and the resulting system of equations couples the Euler (or Navier-Stokes) and the Cahn-Hilliard equations. The model takes into account weak non-locality (dispersion) associated with an internal length scale and localized dissipation due to mixing. The non-locality introduces a dimensional surface energy; dissipation is added to handle the loss of regularity of solutions to the sharp interface equations and to provide a mechanism for topological changes. In particular, we study a nontrivial limit when both components are incompressible, the pressure is kinematic but the velocity field is non-solenoidal (quasi-incompressibility). To demonstrate the effects of quasi-incompressibility, we analyse the linear stage of spinodal decomposition in one dimension. We show that when the densities of the fluids are not perfectly matched, the evolution of the concentration field causes fluid motion even if the fluids are inviscid. In the limit of infinitely thin and well-separated interfacial layers, an appropriately scaled quasi-incompressible Euler-Cahn-Hilliard system converges to the classical sharp interface model. In order to investigate the behaviour of the model outside the range of parameters where the sharp interface approximation is sufficient, we consider a simple example of a change of topology and show that the model permits the transition to occur without an associated singularity.

A new formulation and new methods are presented for computing the motion of fluid interfaces with surface tension in two-dimensional, irrotational, and incompressible fluids. Through the Laplace-Young condition at the interface, surface tension introduces high-order terms, both nonlinear and nonlocal, into the dynamics. This leads to severe stability constraints for explicit time integration methods and makes the application of implicit methods difficult. This new formulation has all the nice properties for time integration methods that are associated with having a linear, constant coefficient, highest order term. That is, using this formulation, we give implicit time integration methods that have no high order time step stability constraint associated with surface tension and are explicit in Fourier space. The approach is based on a boundary integral formulation and applies more generally, even to problems beyond the fluid mechanical context. Here they are applied to computing with high resolution the motion of interfaces in Hele-Shaw flows and the motion of free surfaces in inviscid flows governed by the Euler equations. One Hele-Shaw computation shows the behavior of an expanding gas bubble over long-time as the interface undergoes successive tip-splittings and finger competition. A second computation shows the formation of a very ramified interface through the interaction of surface tension with an unstable density stratification. In Euler flows, the computation of a vortex sheet shows its roll-up through the Kelvin-Helmholtz instability. This motion culminates in the late time self-intersection of the interface, creating trapped bubbles of fluid. This is, we believe, a type of singularity formation previously unobserved for such flows in 2D. Finally, computations of falling plumes in an unstably stratified Boussinesq fluid show a very similar behavior.

We study solid tumor ( carcinoma) growth in the nonlinear regime using boundary-integral simulations. The tumor core is nonnecrotic and no inhibitor chemical species are present. A new formulation of the classical models [18,24,8,3] is developed and it is demonstrated that tumor evolution is described by a reduced set of two dimensionless parameters and is qualitatively unaffected by the number of spatial dimensions. One parameter describes the relative rate of mitosis to the relaxation mechanisms (cell mobility and cell-to-cell adhesion). The other describes the balance between apoptosis (programmed cell-death) and mitosis. Both parameters also include the effect of vascularization. Our analysis and nonlinear simulations reveal that the two new dimensionless groups uniquely subdivide tumor growth into three regimes associated with increasing degrees of vascularization: low (diffusion dominated, e.g., in vitro), moderate and high vascularization, that correspond to the regimes observed in vivo. We demonstrate that critical conditions exist for which the tumor evolves to nontrivial dormant states or grows self-similarly (i.e., shape invariant) in the first two regimes. This leads to the possibility of shape control and of controlling the release of tumor angiogenic factors by restricting the tumor volume-to-surface-area ratio. Away from these critical conditions, evolution may be unstable leading to invasive fingering into the external tissues and to topological transitions such as tumor breakup and reconnection. Interestingly we find that for highly vascularized tumors, while they grow unbounded, their shape always stays compact and invasive fingering does not occur. This is in agreement with recent experimental observations [30] of in vivo tumor growth, and suggests that the invasive growth of highly-vascularized tumors is associated to vascular and elastic anisotropies, which are not included in the model studied here.

This is the first paper in a two-part series in which we develop, analyze, and simulate a diffuse interface continuum model of multispecies tumor growth and tumor-induced angiogenesis in two and three dimensions. Three-dimensional simulations of nonlinear tumor growth and neovascularization using this diffuse interface model were recently presented in Frieboes et al. [2007. Computer simulation of glioma growth and morphology. NeuroImage S59-S70], but that paper did not describe the details of the model or the numerical algorithm. This is done here. In this diffuse interface approach, sharp interfaces are replaced by narrow transition layers that arise due to differential adhesive forces among the cell species. Accordingly, a continuum model of adhesion is introduced. The model is thermodynamically consistent, is related to recently developed mixture models, and thus is capable of providing a detailed description of tumor progression. The model is well-posed and consists of fourth-order nonlinear advection-reaction-diffusion equations (of Cahn-Hilliard-type) for the cell species coupled with reaction-diffusion equations for the substrate components. We demonstrate analytically and numerically that when the diffuse interface thickness tends to zero, the system reduces to a classical sharp interface model. Using a new fully adaptive and nonlinear multigrid/finite difference method, the system is simulated efficiently. In this first paper, we present simulations of unstable avascular tumor growth in two and three dimensions and demonstrate that our techniques now make large-scale three-dimensional simulations of tumors with complex morphologies computationally feasible. In part II of this study, we will investigate multispecies tumor invasion, tumor-induced angiogenesis, and focus on the morphological instabilities that may underlie invasive phenotypes.

Despite major scientific, medical and technological advances over the last few decades, a cure for cancer remains elusive. The disease initiation is complex, and including initiation and avascular growth, onset of hypoxia and acidosis due to accumulation of cells beyond normal physiological conditions, inducement of angiogenesis from the surrounding vasculature, tumour vascularization and further growth, and invasion of surrounding tissue and metastasis. Although the focus historically has been to study these events through experimental and clinical observations, mathematical modelling and simulation that enable analysis at multiple time and spatial scales have also complemented these efforts. Here, we provide an overview of this multiscale modelling focusing on the growth phase of tumours and bypassing the initial stage of tumourigenesis. While we briefly review discrete modelling, our focus is on the continuum approach. We limit the scope further by considering models of tumour progression that do not distinguish tumour cells by their age. We also do not consider immune system interactions nor do we describe models of therapy. We do discuss hybrid-modelling frameworks, where the tumour tissue is modelled using both discrete (cell-scale) and continuum (tumour-scale) elements, thus connecting the micrometre to the centimetre tumour scale. We review recent examples that incorporate experimental data into model parameters. We show that recent mathematical modelling predicts that transport limitations of cell nutrients, oxygen and growth factors may result in cell death that leads to morphological instability, providing a mechanism for invasion via tumour fingering and fragmentation. These conditions induce selection pressure for cell survivability, and may lead to additional genetic mutations. Mathematical modelling further shows that parameters that control the tumour mass shape also control its ability to invade. Thus, tumour morphology may serve as a predictor of invasiveness and treatment prognosis.

In this article, we present a new multiscale mathematical model for solid tumour growth which couples an improved model of tumour invasion with a model of tumour-induced angiogenesis. We perform nonlinear simulations of the multi-scale model that demonstrate the importance of the coupling between the development and remodeling of the vascular network, the blood flow through the network and the tumour progression. Consistent with clinical observations, the hydrostatic stress generated by tumour cell proliferation shuts down large portions of the vascular network dramatically affecting the flow, the subsequent network remodeling, the delivery of nutrients to the tumour and the subsequent tumour progression. In addition, extracellular matrix degradation by tumour cells is seen to have a dramatic affect on both the development of the vascular network and the growth response of the tumour. In particular, the newly developing vessels tend to encapsulate, rather than penetrate, the tumour and are thus less effective in delivering nutrients.

We present an unconditionally energy stable finite-difference scheme for the phase field crystal equation. The method is based on a convex splitting of a discrete energy and is semiimplicit. The equation at the implicit time level is nonlinear but represents the gradient of a strictly convex function and is thus uniquely solvable, regardless of time step size. We present local-in-time error estimates that ensure the convergence of the scheme. While this paper is primarily concerned with the phase field crystal equation, most of the theoretical results hold for the related Swift-Hohenberg equation as well.

We develop a thermodynamically consistent mixture model for avascular solid tumor growth which takes into account the effects of cell-to-cell adhesion, and taxis inducing chemical and molecular species. The mixture model is well-posed and the governing equations are of Cahn–Hilliard type. When there are only two phases, our asymptotic analysis shows that earlier single-phase models may be recovered as limiting cases of a two-phase model. To solve the governing equations, we develop a numerical algorithm based on an adaptive Cartesian block-structured mesh refinement scheme. A centered-difference approximation is used for the space discretization so that the scheme is second order accurate in space. An implicit discretization in time is used which results in nonlinear equations at implicit time levels. We further employ a gradient stable discretization scheme so that the nonlinear equations are solvable for very large time steps. To solve those equations we use a nonlinear multilevel/multigrid method which is of an optimal order O (N) where N is the number of grid points. Spherically symmetric and fully two dimensional nonlinear numerical simulations are performed. We investigate tumor evolution in nutrient-rich and nutrient-poor tissues. A number of important results have been uncovered. For example, we demonstrate that the tumor may suffer from taxis-driven fingering instabilities which are most dramatic when cell proliferation is low, as predicted by linear stability theory. This is also observed in experiments. This work shows that taxis may play a role in tumor invasion and that when nutrient plays the role of a chemoattractant, the diffusional instability is exacerbated by nutrient gradients. Accordingly, we believe this model is capable of describing complex invasive patterns observed in experiments.