The Finite Element Method Applied to a Problem of Blood Flow in Vessels

We use the finite element method to solve a convection-diffusion equation when convection is dominating, a problem which describes the behavior of the concentration of a solute in a blood vessel. A new technique for computing the discrete problem is used

One-Dimensional Navier-Stokes Finite Element Flow Model

This technical report documents the theoretical, computational, and practical aspects of the one-dimensional Navier-Stokes finite element flow model. The document is particularly useful to those who are interested in implementing, validating and utilizing this relatively-simple and widely-used model.Comment: 46 pages, 1 tabl

Application of the penalty coupling method for the analysis of blood vessels

Due to the significant health and economic impact of blood vessel diseases on modern society, its analysis is becoming of increasing importance for the medical sciences. The complexity of the vascular system, its dynamics and material characteristics all make it an ideal candidate for analysis through fluid structure interaction (FSI) simulations. FSI is a relatively new approach in numerical analysis and enables the multi-physical analysis of problems, yielding a higher accuracy of results than could be possible when using a single physics code to analyse the same category of problems. This paper introduces the concepts behind the Arbitrary Lagrangian Eulerian (ALE) formulation using the penalty coupling method. It moves on to present a validation case and compares it to available simulation results from the literature using a different FSI method. Results were found to correspond well to the comparison case as well as basic theory

Non-Newtonian Rheology in Blood Circulation

Blood is a complex suspension that demonstrates several non-Newtonian rheological characteristics such as deformation-rate dependency, viscoelasticity and yield stress. In this paper we outline some issues related to the non-Newtonian effects in blood circulation system and present modeling approaches based mostly on the past work in this field.Comment: 26 pages, 5 figures, 2 table

Multiphase modelling of vascular tumour growth in two spatial dimensions

In this paper we present a continuum mathematical model of vascular tumour growth which is based on a multiphase framework in which the tissue is decomposed into four distinct phases and the principles of conservation of mass and momentum are applied to the normal/healthy cells, tumour cells, blood vessels and extracellular material. The inclusion of a diffusible nutrient, supplied by the blood vessels, allows the vasculature to have a nonlocal influence on the other phases. Two-dimensional computational simulations are carried out on unstructured, triangular meshes to allow a natural treatment of irregular geometries, and the tumour boundary is captured as a diffuse interface on this mesh, thereby obviating the need to explicitly track the (potentially highly irregular and ill-defined) tumour boundary. A hybrid finite volume/finite element algorithm is used to discretise the continuum model: the application of a conservative, upwind, finite volume scheme to the hyperbolic mass balance equations and a finite element scheme with a stable element pair to the generalised Stokes equations derived from momentum balance, leads to a robust algorithm which does not use any form of artificial stabilisation. The use of a matrix-free Newton iteration with a finite element scheme for the nutrient reaction-diffusion equations allows full nonlinearity in the source terms of the mathematical model. Numerical simulations reveal that this four-phase model reproduces the characteristic pattern of tumour growth in which a necrotic core forms behind an expanding rim of well-vascularised proliferating tumour cells. The simulations consistently predict linear tumour growth rates. The dependence of both the speed with which the tumour grows and the irregularity of the invading tumour front on the model parameters are investigated