Understanding tumor induced angiogenesis is a challenging problem with important consequences for diagnosis and treatment of cancer. Recently, strong evidences suggest the dual role of endothelial cells on the migrating tips and on the proliferating body of blood vessels, in consonance with further events behind lumen formation and vascular patterning. In this paper we present a multi-scale phase-field model that combines the benefits of continuum physics description and the capability of tracking individual cells. The model allows us to discuss the role of the endothelial cells' chemotactic response and proliferation rate as key factors that tailor the neovascular network. Importantly, we also test the predictions of our theoretical model against relevant experimental approaches in mice that displayed distinctive vascular patterns. The model reproduces the in vivo patterns of newly formed vascular networks, providing quantitative and qualitative results for branch density and vessel diameter on the order of the ones measured experimentally in mouse retinas. Our results highlight the ability of mathematical models to suggest relevant hypotheses with respect to the role of different parameters in this process, hence underlining the necessary collaboration between mathematical modeling, in vivo imaging and molecular biology techniques to improve current diagnostic and therapeutic tools.

We study the viscous fingering or Saffman–Taylor instability in two different dilute or semi-dilute polymer solutions. The different solutions exhibit only one non-Newtonian property, in the sense that other non-Newtonian effects can be neglected. The viscosity of solutions of stiff polymers has a strong shear rate dependence. Relative to Newtonian fluids, narrower fingers are found for rigid polymers. For solutions of flexible polymers, elastic effects such as normal stresses are dominant, whereas the shear viscosity is almost constant. Wider fingers are found in this case. We characterize the non-Newtonian flow properties of these polymer solutions completely, allowing for separate and quantitative investigation of the influence of the two most common non-Newtonian properties on the Saffman–Taylor instability. The effects of the non-Newtonian flow properties on the instability can in all cases be understood quantitatively by redefining the control parameter of the instability.

Haemodynamic simulations using one-dimensional (1-D) computational models exhibit many of the features of the systemic circulation under normal and diseased conditions. We propose a novel linear 1-D dynamical theory of blood flow in networks of flexible vessels that is based on a generalized Darcy's model and for which a full analytical solution exists in frequency domain. We assess the accuracy of this formulation in a series of benchmark test cases for which computational 1-D and 3-D solutions are available. Accordingly, we calculate blood flow and pressure waves, and velocity profiles in the human common carotid artery, upper thoracic aorta, aortic bifurcation, and a 20-artery model of the aorta and its larger branches. Our analytical solution is in good agreement with the available solutions and reproduces the main features of pulse waveforms in networks of large arteries under normal physiological conditions. Our model reduces computational time and provides a new approach for studying arterial pulse wave mechanics; e.g.,  the analyticity of our model allows for a direct identification of the role played by physical properties of the cardiovascular system on the pressure waves.

Rare events appear in a wide variety of phenomena such as rainfall, floods, earthquakes, and risk. We demonstrate that the stochastic behavior induced by the natural roughening present in standard microchannels is so important that the dynamics for the advancement of a water front displacing air has plenty of rare events. We observe that for low pressure differences the hydrophobic interactions of the water front with the walls of the microchannel put the front close to the pinning point. This causes a burstlike dynamics, characterized by series of pinning and avalanches, that leads to an extreme-value Gumbel distribution for the velocity fluctuations and a nonclassical time exponent for the advancement of the mean front position as low as 0.38.