UCL Discovery
UCL home » Library Services » Electronic resources » UCL Discovery

Choice of boundary condition for lattice-Boltzmann simulation of moderate-Reynolds-number flow in complex domains

Nash, RW; Carver, HB; Bernabeu, MO; Hetherington, J; Groen, D; Krueger, T; Coveney, PV; (2014) Choice of boundary condition for lattice-Boltzmann simulation of moderate-Reynolds-number flow in complex domains. PHYSICAL REVIEW E , 89 (2) , Article 023303. 10.1103/PhysRevE.89.023303. Green open access

[thumbnail of PhysRevE.89.023303.pdf] PDF
PhysRevE.89.023303.pdf

Download (3MB)

Abstract

Modeling blood flow in larger vessels using lattice-Boltzmann methods comes with a challenging set of constraints: a complex geometry with walls and inlets and outlets at arbitrary orientations with respect to the lattice, intermediate Reynolds (Re) number, and unsteady flow. Simple bounce-back is one of the most commonly used, simplest, and most computationally efficient boundary conditions, but many others have been proposed. We implement three other methods applicable to complex geometries [Guo, Zheng, and Shi, Phys. Fluids 14, 2007(2002); Bouzidi, Firdaouss, and Lallemand, Phys. Fluids 13, 3452 (2001); Junk and Yang, Phys.Rev.E 72,066701(2005)] in our open-source application HEMELB. We use these to simulate Poiseuille and Womersley flows in a cylindrical pipe with an arbitrary orientation at physiologically relevant Re number (1–300) and Womersley (4–12) numbers and steady flow in a curved pipe at relevant Dean number (100–200) and compare the accuracy to analytical solutions. We find that both the Bouzidi-Firdaouss-Lallemand (BFL) and Guo-Zheng-Shi (GZS) methods give second-order convergence in space while simple bounce-back degrades to first order. The BFL method appears to perform better than GZS in unsteady flows and is significantly less computationally expensive. The Junk-Yang method shows poor stability at larger Re number and so cannot be recommended here. The choice of collision operator (lattice Bhatnagar-Gross-Krook vs multiple relaxation time) and velocity set (D3Q15 vs D3Q19 vs D3Q27) does not significantly affect the accuracy in the problems studied.

Type: Article
Title: Choice of boundary condition for lattice-Boltzmann simulation of moderate-Reynolds-number flow in complex domains
Open access status: An open access version is available from UCL Discovery
DOI: 10.1103/PhysRevE.89.023303
Publisher version: http://dx.doi.org/10.1103/PhysRevE.89.023303
Language: English
Additional information: Published by the American Physical Society under the terms of the Creative Commons Attribution 3.0 License. Further distribution of this work must maintain attribution to the author(s) and the published article's title, journal citation, and DOI.
UCL classification: UCL
UCL > Provost and Vice Provost Offices > UCL BEAMS
UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences
UCL > Provost and Vice Provost Offices > UCL BEAMS > Faculty of Maths and Physical Sciences > Dept of Chemistry
URI: https://discovery.ucl.ac.uk/id/eprint/1375762
Downloads since deposit
196Downloads
Download activity - last month
Download activity - last 12 months
Downloads by country - last 12 months

Archive Staff Only

View Item View Item