Click here for a recent version of the code from the figshare website.

This code can be used to perform two-phase flow simulations on micro-CT images of porous media. It is based on the algorithm published in Raeini et al. (2012) [doi] [PDF download], with some modifications for improving the stability of the code for unstructured meshes.

It includes the required pre- and post-processing tools, sample input files and a short documentation. For examples of the applications of this code, please see the publications by A.Q. Raeini, B. Bijeljic and M.J. Blunt, for studying pore-scale events such as snap-off and piston-like displacement (2014a) [DOI] [PDF download], upscaling two-phase flow properties (2014b) [DOI] [PDF download] and studying capillary trapping on micro-CT images (2015) [DOI].

A previous version (2011-2012) of the direct two-phase flow solver used for modelling two-phase flow at the pore/micron scale can be downloaded.

The algorithm implemented in this code is closer to the algorithm published in Raeini et al. (2012) [doi] [PDF download].

The code should be linked to OpenFOAM; it has been tested to work with OpenFOAM version 1.6-ext, please see OpenFOAM-extend website or its SourceForge repository.

For any enquiry, please contact Ali Raeini: a.qaseminejad-raeini09@imperial.ac.uk.