Parametric uncertainty analysis of pulse wave propagation in a model of a human arterial network

Accepted versio

Multiscale modelling of vascular tumour growth in 3D: the roles of domain size & boundary condition

We investigate a three-dimensional multiscale model of vascular tumour growth, which couples blood flow, angiogenesis, vascular remodelling, nutrient/growth factor transport, movement of, and interactions between, normal and tumour cells, and nutrient-dependent cell cycle dynamics within each cell. In particular, we determine how the domain size, aspect ratio and initial vascular network influence the tumour's growth dynamics and its long-time composition. We establish whether it is possible to extrapolate simulation results obtained for small domains to larger ones, by constructing a large simulation domain from a number of identical subdomains, each subsystem initially comprising two parallel parent vessels, with associated cells and diffusible substances. We find that the subsystem is not representative of the full domain and conclude that, for this initial vessel geometry, interactions between adjacent subsystems contribute to the overall growth dynamics. We then show that extrapolation of results from a small subdomain to a larger domain can only be made if the subdomain is sufficiently large and is initialised with a sufficiently complex vascular network. Motivated by these results, we perform simulations to investigate the tumour's response to therapy and show that the probability of tumour elimination in a larger domain can be extrapolated from simulation results on a smaller domain. Finally, we demonstrate how our model may be combined with experimental data, to predict the spatio-temporal evolution of a vascular tumour

Blood flow in microvascular networks

This paper was presented at the 2nd Micro and Nano Flows Conference (MNF2009), which was held at Brunel University, West London, UK. The conference was organised by Brunel University and supported by the Institution of Mechanical Engineers, IPEM, the Italian Union of Thermofluid dynamics, the Process Intensification Network, HEXAG - the Heat Exchange Action Group and the Institute of Mathematics and its Applications.Simulation of blood presents a very complex haemodynamics problem especially in relation to the understanding of atherogenesis. In many simulations, blood has been treated as a single-phase homogeneous fluid, a classical approach that does not account for the presence of red blood cells (RBCs). Although this approach provides satisfactory tools to describe certain aspects of blood flow in large arteries, it fails to give an adequate representation of the flow field in the vessels of smaller diameter where the size of the RBC becomes significant relative to vessel diameter. So, this article is concerned with the study of non-Newtonian blood flow in microvascular networks with the intention of developing a new cell depletion layer model to represent the behaviour of RBCs through bifurcating networks. The model is tested in a microvascular network constructed possessing realistic bifurcation features, with controlled dimensions and angles. The RBC depletion model treats blood as two continuum layers, with a central, non-Newtonian core region of concentrated red cell suspension that is surrounded by a layer of plasma (Newtonian fluid) adjacent to the vessel wall. In the central core region, blood is described by Quemada's non-Newtonian rheological model. Geometry differences are shown to significantly affect flow rates, haematocrit distributions and the corresponding cell depletion layers

Blood Vessel Tortuosity Selects against Evolution of Agressive Tumor Cells in Confined Tissue Environments: a Modeling Approach

Cancer is a disease of cellular regulation, often initiated by genetic mutation within cells, and leading to a heterogeneous cell population within tissues. In the competition for nutrients and growth space within the tumors the phenotype of each cell determines its success. Selection in this process is imposed by both the microenvironment (neighboring cells, extracellular matrix, and diffusing substances), and the whole of the organism through for example the blood supply. In this view, the development of tumor cells is in close interaction with their increasingly changing environment: the more cells can change, the more their environment will change. Furthermore, instabilities are also introduced on the organism level: blood supply can be blocked by increased tissue pressure or the tortuosity of the tumor-neovascular vessels. This coupling between cell, microenvironment, and organism results in behavior that is hard to predict. Here we introduce a cell-based computational model to study the effect of blood flow obstruction on the micro-evolution of cells within a cancerous tissue. We demonstrate that stages of tumor development emerge naturally, without the need for sequential mutation of specific genes. Secondly, we show that instabilities in blood supply can impact the overall development of tumors and lead to the extinction of the dominant aggressive phenotype, showing a clear distinction between the fitness at the cell level and survival of the population. This provides new insights into potential side effects of recent tumor vasculature renormalization approaches

A multiple scale model for tumor growth

We present a physiologically structured lattice model for vascular tumor growth which accounts for blood flow and structural adaptation of the vasculature, transport of oxygen, interaction between cancerous and normal tissue, cell division, apoptosis, vascular endothelial growth factor release, and the coupling between these processes. Simulations of the model are used to investigate the effects of nutrient heterogeneity, growth and invasion of cancerous tissue, and emergent growth laws

Molecular Dynamics Simulation of Vascular Network Formation

Endothelial cells are responsible for the formation of the capillary blood vessel network. We describe a system of endothelial cells by means of two-dimensional molecular dynamics simulations of point-like particles. Cells' motion is governed by the gradient of the concentration of a chemical substance that they produce (chemotaxis). The typical time of degradation of the chemical substance introduces a characteristic length in the system. We show that point-like model cells form network resembling structures tuned by this characteristic length, before collapsing altogether. Successively, we improve the non-realistic point-like model cells by introducing an isotropic strong repulsive force between them and a velocity dependent force mimicking the observed peculiarity of endothelial cells to preserve the direction of their motion (persistence). This more realistic model does not show a clear network formation. We ascribe this partial fault in reproducing the experiments to the static geometry of our model cells that, in reality, change their shapes by elongating toward neighboring cells.Comment: 10 pages, 3 figures, 2 of which composite with 8 pictures each. Accepted on J.Stat.Mech. (2009). Appeared at the poster session of StatPhys23, Genoa, Italy, July 13 (2007