Code based on the OpenFOAM(R) toolbox