Modelling two-dimensional complex flow using the lattice Boltzmann method

Lattice Boltzmann method is a computational technique that streams and collides "psuedo particles", or particle frequency distributions, on a lattice to approximate hydrodynamics.

• Julia language
• Python language

• PyCall
• PyPlot
• HDF5
• JLD
• ArgParse
• ProfileView
• Roots
• FastAnonymous

• PyYaml (python library)
• numpy (python library)
• matplotlib (python library)

Install the both the julia language and python. Note: the install script assumes an Ubuntu environment and that pip is installed.

git clone https://github.com/grasingerm/lbxflow.git
cd lbxflow
julia scripts/install_dependencies.jl
julia lbxflow.jl --help
julia lbxflow.jl -vf example_sims/poiseuille/poiseuille_velocity-profile.yaml

An introductory tutorial is available in the doc directory. Currently documentation is sparse, but the source files are all well documented and I am responsive to questions and feedback. Example input files are available in the example_sims directory.

Please report all bugs and issues using the github ticketing system.

• add capability for 3D simulations
• add animations and plots using gnuplot and Gadfly
• implement free surface flow
• break collision functions down into forcing method, constitutive eqn, etc.
• create alt. stream/collide methods using grids of obstacle flags
• write user documentation
• parallelize streaming and collision steps, multiscale mapping, and BCs
• wrap simulation code in a module
• profile simulation steps to identify bottlenecks
• write some tests?
• write mrt collision functions with "hard coded" default inverse transformation matrix
• verify Newtonian Poiseuille flow
• debug BCs at inlet/outlet (bounceback, periodic)
• implement MRT scheme for Bingham plastic flow
• clean up api
• make an autosave feature
• should be able to "catch" from exceptions and/or keyboard interrupt
• should be able to restart from interrupted model where left off
• clean up existing BCs / more general BCs
• develop post processing and animations
• implement HB-Bingham constitutive model used in Fluent