porefoam1f

porefoam1f

Porefoam1f code solves for incompressible single-phase (1f) flow on 3D images of porous media using openfoam finite-volume library. The code computes and reports flow properties, such as the absolute permeability and formation factor, of the pore space.

---

Prerequisites

• A GNU Linux operating system (at least for compilation), tested on Ubuntu 18.04
• GNU make, cmake, and a c++11 compiler (available in most Linux distributions or can be installed through their package manager)
• MPI message passing interface. In Ubuntu (18.04) MPI another foamx4m prerequisites can be installed by typing in a terminal:
  sudo apt install mpi-default-dev flex libscotch-dev
• libvoxel (included) and its prerequisies:
  sudo apt install libjpeg-dev liblzma-dev
• A costomized foam-extend, foamx4m (included)
• libtiff and zlib (optional, both dependancies of libvoxel, included)
• script for convenient compilation (included)

Instructions

Download and compile

This repository is part of a larger group of partly independent packages. Please follow the instruction in README.md in the upper-most directory.

Running simulations

Please see the src/doc folder for installation and usage, and a sample input header file -- Image.mhd, which you can you in conjunction with segmented micro-CT image files in raw, raw.gz, tif, and amira (.am) files. For sample micro-CT images and their mhd header files, see Imperial College pore-scale modelling website.

In summary, you need to download a segmented micro-CT image (with a image.mhd header) and run in a bash terminal:

    source PATH/TO/src/script/bashrc # only once, change PATH/TO according to your porefoam installaion path
    #cd PATH/TO/IMAGES/
    # Set the number of processors (=$nProcX x $nProcY x $nProcZ) used to run the simulation,
    # based on the size of image and number of processors your machine:
    export nProcX=2;  export nProcY=2;  export nProcZ=2  
    AllRunImagePar "$(ls *.mhd)"  "X Y Z" # Run the simulations, 
    # This command runs flow simulations on all available mhd files in the current directory, in all 3 spatial directions, and will take a while.

Simulation results

The 3D simulation results are saved in OpenFOAM format, as well as converted into 3D images, both can be visualized using [Paraview]. The Input and output 3D images can also be viwed using [Fiji-is-ImageJ] software. A summary of the simulation results, including connected porosity, permeability, formation factor and velocity distributions, are saved in a text file named summary_...txt, alongside other log files. These log files should be monitoed in case of any anomoly in the produced results.

Contact and References

To report any problems, contact Sajjad Foroughi, email: [email protected]

For more contacts and references please visit:
https://www.imperial.ac.uk/earth-science/research/research-groups/pore-scale-modelling

License:
GPLv3