Pore Morphology Method

This code calculates the pressure-saturation relations using the Pore Morphology Method with trapped mechanism.

A computationally efficient method to calculate the capillary pressure-saturation relations of im-miscible multiphase flow on two-dimensional pore morphologies is presented in this code The method is an extension of the porous morphology method that includes the wetting angle and trapped mechanism of the displaced flui. Because the method uses only morphological binary operations, it is more efficient than well-established high-resolution voxel dynamics methods such as Lattice-Boltzmann methods and level-set computational fluid dynamics. Please run the "RunME" code and enjoy.