Code to the article: Patch-based multiscale algorithm for flow and reactive transport in fracture-microcrack systems in shales

A major challenge in modeling reactive transport in shales is the multiscale nature of the fractured rock system. Shales may contain a few main fractures and thousands of microcracks whose length and aperture are orders of magnitude smaller than the former. This renders fully resolved simulations too expensive, while traditional upscaling methods cannot accurately capture fracture clogging dynamics due to precipitation. We develop a fully coupled patch-based multiscale algorithm to address this issue: we solve reactive transport only in a few selected microcracks, while evaluating the others by interpolation and ensuring top-down bottom-up coupling between the main fracture and the microcracks.