Direct two phase flow solver
The latest experimental version of the code based on the conservative curvature force formulation by Shams et al (2018) and the filtering of surface forces by Raeini et al. (2012) [doi] [PDF] is available from github:
See this Fishare link for the 2015 version of porefoam code which is based on the algorithm published in Raeini et al. (2012) [doi] [PDF] with modifications to improve the stability of the code for unstructured meshes. It includes the required pre- and post-processing tools, sample input files and a short documentation. It has been used to study pore-scale events such as snap-off and piston-like displacement (Raeini et al. 2014a) [DOI] [PDF], to upscale pore-scale forces in two-phase flow (Raeini et al. 2014b) [DOI] [PDF] and to study capillary trapping on micro-CT images of porous media (Raeini et al. 2015) [DOI].
An older version (2011-2012) of the direct two-phase flow solver used for modelling two-phase flow at the pore/micron scale can be downloaded here. It is closer to the algorithm published in Raeini et al. (2012) [DOI] [PDF] than the links above. 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.