This is a Boussinesq-type equation solver built within the OpenFOAM frame. - View it on GitHub