In the subsurface, fluids play a critical role by transporting dissolved minerals, colloids, and contaminants (sometimes over long distances); by mediating dissolution and precipitation processes; and by enabling chemical transformations in solution and at mineral surfaces. Although the complex geometries of fracture apertures, fracture networks, and pore spaces may make it difficult to accurately predict fluid flow in saturated (single‐phase) subsurface systems, well‐developed methods are available. The simulation of multiphase fluid flow in the subsurface is much more challenging because of the large density and/or viscosity ratios found in important applications (water/air in the vadose zone; water/oil, water/gas, gas/oil, and water/oil/gas in hydrocarbon reservoirs; water/air/nonaqueous phase liquids (nonaqueous phase liquids/dense nonaqueous phase liquids) in contaminated vadose zone systems; and gas/molten rock in volcanic systems, for example). In addition, the complex behavior of fluid‐fluid‐solid contact lines and their impact on dynamic contact angles must also be taken into account and coupled with the fluid flow. Here we review the methods that are currently being used to simulate pore‐scale multiphase fluid flow and reactive transport in fractured and porous media. After the introduction, the review begins with an overview of the fundamental physics of multiphase fluids flow followed by a more detailed discussion of the complex dynamic behavior of contact lines and contact angles, an important barrier to accurate pore‐scale modeling and simulation. The main part of the review focuses on five different approaches: pore network models, lattice gas and lattice Boltzmann methods, Monte Carlo methods, particle methods (molecular dynamics, dissipative particle dynamics, and smoothed particle hydrodynamics), and traditional grid‐based computational fluid dynamics coupled with interface tracking and a contact angle model. Finally, the review closes with a discussion of future trends and challenges.
