A Patient-Specific Anisotropic Diffusion Model for Brain Tumour Spread.

Gliomas are primary brain tumours arising from the glial cells of the nervous system. The diffuse nature of spread, coupled with proximity to critical brain structures, makes treatment a challenge. Pathological analysis confirms that the extent of glioma spread exceeds the extent of the grossly visible mass, seen on conventional magnetic resonance imaging (MRI) scans. Gliomas show faster spread along white matter tracts than in grey matter, leading to irregular patterns of spread. We propose a mathematical model based on Diffusion Tensor Imaging, a new MRI imaging technique that offers a methodology to delineate the major white matter tracts in the brain. We apply the anisotropic diffusion model of Painter and Hillen (J Thoer Biol 323:25-39, 2013) to data from 10 patients with gliomas. Moreover, we compare the anisotropic model to the state-of-the-art Proliferation-Infiltration (PI) model of Swanson et al. (Cell Prolif 33:317-329, 2000). We find that the anisotropic model offers a slight improvement over the standard PI model. For tumours with low anisotropy, the predictions of the two models are virtually identical, but for patients whose tumours show higher anisotropy, the results differ. We also suggest using the data from the contralateral hemisphere to further improve the model fit. Finally, we discuss the potential use of this model in clinical treatment planning.

Robustness of the inverse cascade in two-dimensional turbulence.

We study quasisteady inverse cascades in unbounded and bounded two-dimensional turbulence driven by time-independent injection and dissipated by molecular viscosity. It is shown that an inverse cascade that carries only a fraction r of the energy input to the largest scales requires the enstrophy-range energy spectrum to be steeper than k(-5) (ruling out a direct cascade) unless 1-r<<1. A direct cascade requires the presence of an inverse cascade that carries virtually all energy to the largest scales (1-r<<1). These facts underlie the robustness of the Kolmogorov-Kraichnan k(-5/3) inverse cascade, which is readily observable in numerical simulations without an accompanying direct enstrophy cascade. We numerically demonstrate an instance where the k(-5/3) inverse-cascading range is realizable with 79% of the energy injection dissipated within the energy range and virtually all of the enstrophy dissipated in the vicinity of the forcing region. As equilibrium is approached, the respective logarithmic slopes -alpha and -beta of the ranges of wave numbers lower and higher than the forcing wave number satisfy alpha+beta approximately 8. These results are consistent with recent theoretical predictions.

Energy budgets in Charney-Hasegawa-Mima and surface quasigeostrophic turbulence.

We study energy transfer in unbounded Charney-Hasegawa-Mima and surface quasigeostrophic turbulence. The possible inverse-cascading quantities in these systems are, respectively, I identical with integral ( infinity )(0)k(-2)E(k) dk and J identical with integral ( infinity )(0)k(-1)E(k) dk, where E(k) is the kinetic energy spectrum. The supposed direct-cascading quantities for both surface quasigeostrophic and Navier-Stokes turbulence are shown to be bounded. We derive a constraint on E(k) for the surface quasigeostrophic system.