No Cover Image

Journal article 1340 views 832 downloads

A benchmark study of numerical schemes for one-dimensional arterial blood flow modelling

Etienne Boileau, Perumal Nithiarasu Orcid Logo, Pablo J Blanco, Lucas O. Müller, Fredrik Eikeland Fossan, Leif Rune Hellevik, Wouter P. Donders, Wouter Huberts, Marie Willemet, Jordi Alastruey

International Journal for Numerical Methods in Biomedical Engineering

Swansea University Author: Perumal Nithiarasu Orcid Logo

DOI (Published version): 10.1002/cnm.2732

Abstract

Hæmodynamical simulations using one-dimensional (1-D) computational models exhibit many of the features of the systemic circulation under normal and diseased conditions. Recent interest in verifying 1-D numerical schemes has led to the development of alternative experimental setups and the use of 3-...

Full description

Published in: International Journal for Numerical Methods in Biomedical Engineering
Published: 2015
URI: https://cronfa.swan.ac.uk/Record/cronfa22154
Tags: Add Tag
No Tags, Be the first to tag this record!
Abstract: Hæmodynamical simulations using one-dimensional (1-D) computational models exhibit many of the features of the systemic circulation under normal and diseased conditions. Recent interest in verifying 1-D numerical schemes has led to the development of alternative experimental setups and the use of 3-D numerical models to acquire data not easily measured in vivo. In most studies to date, only one particular 1-D scheme is tested. In this paper we present a systematic comparison of six commonly used numerical schemes for 1-D blood flow modelling: discontinuous Galerkin (DCG), locally conservative Galerkin (LCG), Galerkin least-squares finite element method (FEM), finite volume method (FVM), finite difference MacCormack method (McC), and a simplified trapezium rule method (STM). Comparisons are made in a series of six benchmark test cases with an increasing degree of complexity. The accuracy of the numerical schemes is assessed by comparison against theoretical results, 3-D numerical data in compatible domains with distensible walls, or experimental data in a network of silicone tubes. Results show a good agreement among all numerical schemes and their ability to capture the main features of pressure, flow and area waveforms in large arteries. All the information used in this study, including the input data for all benchmark cases, experimental data where available, and numerical solutions for each scheme, is made publicly available online, providing a comprehensive reference data set to support the development of 1-D models and numerical schemes. This article is protected by copyright. All rights reserved.
College: Faculty of Science and Engineering