Error analysis and correction for Lattice Boltzmann simulated flow conductance in capillaries of different shapes and alignments

  • Authors:
  • Ankush Sengupta;Paul S. Hammond;Daan Frenkel;Edo S. Boek

  • Affiliations:
  • Fluids Department, Schlumberger Cambridge Research Centre, High Cross, Madingley Road, Cambridge CB3 0EL, United Kingdom and Department of Chemistry, University of Cambridge, Lensfield Road, Cambr ...;Fluids Department, Schlumberger Cambridge Research Centre, High Cross, Madingley Road, Cambridge CB3 0EL, United Kingdom;Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom;Department of Chemical Engineering, Imperial College London, South Kensington Campus, London SW7 2AZ, United Kingdom

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2012

Quantified Score

Hi-index 31.45

Visualization

Abstract

We study the relative error in conductance calculations, for simulated flow of a single component single phase fluid through a capillary in three dimensions, by the Lattice Boltzmann (LB) method with bounce-back boundary conditions. The relative error with respect to analytical results for capillary cross-sections of circular, triangular and square shapes are calculated as a function of the cross-section diameter, a, and for different alignment of the cross-section relative to the underlying lattice grid. It is shown, when the shapes are not aligned perfectly to the lattice, that the relative error decreases systematically with the size, a, as ~1/a when a is evaluated by mapping the computed cross-sectional area, in terms of the enclosed number of grid points, to the respective geometrical shapes concerned. For perfectly aligned geometries, viz. the square capillary aligned to the LB lattice grid or rotated with its side along the diagonal of the LB grid, the relative error decreases as ~1/a^2. A simple method is suggested to locate the boundary wall depending on its orientation relative to the grid, such that the exact conductance of the new shape matches the LB computed conductance.