Measurements of wall shear stress with the lattice Boltzmann method and staircase approximation of boundaries

• Published in: Computers & fluids Vol. 39; no. 9; pp. 1625 - 1633
• Main Authors: Stahl, B., Chopard, B., Latt, J.
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Ltd 01-10-2010; Elsevier
• Subjects: Accuracy; Approximation; Biological and medical sciences; Blood and lymphatic vessels; Boundaries; Cardiology. Vascular system; CFD; Computational methods in fluid dynamics; Computer simulation; Diseases of the peripheral vessels. Diseases of the vena cava. Miscellaneous; Exact sciences and technology; Flows in ducts, channels, nozzles, and conduits; Fluid dynamics; Fundamental areas of phenomenology (including applications); Lattice Boltzmann; Lattices; Mathematical analysis; Medical sciences; Physics; Voxel geometry; Wall shear stress; Wall shear stresses; Walls; CFD; Lattice Boltzmann; Wall shear stress; Voxel geometry; Boltzmann equation; Pipe flow; Computational fluid dynamics; Digital simulation; Aneurysm; Cardiovascular disease; Vascular disease; Blood circulation; Voxel; Lattice model; Curved pipe; Shear stress; Modelling; Inclined pipe
• DOI: 10.1016/j.compfluid.2010.05.015
• EISSN: 1879-0747
• ISSN: 0045-7930

We analyze the accuracy of wall shear stress measurements in lattice Boltzmann simulations that are based on a voxel representation of the geometry and staircase approximation of boundaries. Such approximations are commonly used in the context of lattice Boltzmann simulations, because they favor the use of simple and highly efficient data structures. We show on several two- and three-dimensional simulations that this low-order approximation of the boundary affects the accuracy of wall shear stress measurements in areas directly adjacent to the wall. A few lattice nodes apart from the wall, the accuracy is however largely improved, and can be considered to be compatible with the overall accuracy of a simulation at a given coarseness level of the grid. This result is interpreted as a justification for the use of walls with staircase shape, even in simulations with high expectations regarding the level of accuracy. Furthermore, we propose a novel method for establishing the direction of the wall normal, a quantity which is required for the computation of the wall shear stress. With this method, the wall normal is computed from local data that is extracted from the results of the fluid flow simulation. Owing to the nature of the flow dynamics, which tends to smooth out the asperities of the wall, the information on the wall orientation obtained in this way is observed to be of high quality.