Publication Date:
2019
Description:
〈span〉〈div〉Abstract〈/div〉The Lattice Boltzmann Method (LBM) is a method to simulate fluid dynamics based on modelling distributions of particles moving and colliding on a lattice. The 〈span〉Python〈/span〉 scripting language provides a clean programming paradigm to develop codes based on the Lattice Boltzmann Method, however in order to reach performance comparable to compiled languages, it needs to be carefully implemented, maximizing its vectorized tools, mostly integrated in the NumPy module. We present here the details of a Python implementation of a concise LBM code, with the purpose of offering a pedagogical tool for students and professionals in the geosciences who are approaching this technique for the first time. The first half of the paper focuses on how to vectorize a 2D LBM code and show how if carefully done, this allows performance close to a compiled code. In the second part of the paper, we use the vectorization described earlier to naturally write a parallel implementation using MPI and test both weak and hard scaling up to 1280 cores. One benchmark, Poiseuille flow, and two applications, one on sound wave propagation and another on fluid-flow through a simplified model of a rock matrix are finally shown.〈/span〉
Print ISSN:
2051-1965
Electronic ISSN:
1365-246X
Topics:
Geosciences
Published by
Oxford University Press
on behalf of
The Deutsche Geophysikalische Gesellschaft (DGG) and the Royal Astronomical Society (RAS).
