The computational studies employ various numerical techniques, combining mesoscopic modelling using Lattice Boltzmann Method (LBM), multi-physics modelling and conventional CFD to investigate the two-phase flow physics, and fluid-fluid and fluid-rock interactions at sub-pore levels:

i) An efficient two-phase flow LBM model for application of modelling CO2 migration in brine has been developed based on evaluating, comparing and improving the various existing models, which will be optimised for the particular fluid properties and thermodynamic conditions.

ii) a finite-volume solver for the fundamental equations governing the basic flow phenomena is being developed based on first principles to describe flows at sub-pore levels, but with simplifications taking advantage of the low flow rate conditions to reduce computing resources. Both of the above solvers will then be used to study the physical problems and generate further detailed information which is not available from experiments.

iii) Exercises using CFD simulations with commercial software are also being undertaken which will produce complementary data to compare with our new methods.

  • LBM for Brine/CO2 Flow in Rocks (Sheffield & Leeds)



Numerical rock (in green) (left) and phase distributions of brine and CO2 (in blue) (middle and right).


Relative permeabilities of brine and CO2 obtained by the dynamic saturation of fluids in the rock with 0.5 mm and resolution 2.5 µm (left) and in the rock with 1.0 mm and resolutions 5.0 µm, 3.3 µm and 2.5 µm.

  • CFD of CO2 Migration and Fluid-rock Reaction (Tsinghua, Leeds, Aberdeen & Sheffield)


Pore scale modelling of the effect of inlet velocity on water saturation.

  • CFD Based on Stoke Flow (Aberdeen & Sheffield)


Predicted pressure, velocity and permeability tensor in a 5 µm resolution image.