Transport-Controlled Metamorphic Reaction Localization and Completion: A Scaling Framework Coupling Deformation, Diffusion, and Drainage

Full description

Research context (extract from the abstract of the associated publication)

Metamorphic reactions commonly localize into narrow fronts and remain incomplete, producing spatially heterogeneous textures under broadly similar pressure–temperature conditions. These observations indicate that reaction affinity alone does not determine reaction progress; the redistribution and retention of reactive components are also required. This study develops a compact scaling framework based on two measures: a reaction–transport length scale (λ), defined by the balance between effective diffusivity and effective reaction rate, and a reaction–drainage stability number (Λ), which compares reaction advance with advective removal of components.

What this record contains

Four standalone Python scripts implementing a hierarchy of one-dimensional numerical models of coupled transport and reaction in metamorphic rocks, together with the code that generates the published figures:

• Model I — transient diffusion–reaction; recovers the reaction localization thickness λ. Generates Figure 3.
• Model II — steady advection–diffusion–reaction; tests the collapse of reaction progress onto the reaction–drainage stability number Λ. Generates Figure 4.
• Model III — porosity-based deformation-dependent transport; porosity is integrated from the volumetric strain rate, with permeability and diffusivity following porosity power laws. Generates Figure 5.
• Model IV — reaction-generated porosity and transformation-dependent transport; produces spatially variable reaction–drainage fields. Generates Figure 6.

How to use it

Each script is standalone and runs under Python 3 with NumPy, SciPy and Matplotlib. Running a script reproduces its manuscript figure in PNG, PDF and SVG, and prints the diagnostic values quoted in the text. Every parameter is declared at the top of each file with its units and its literature source. No input data files are required.

The full methodology is available in the publication linked under Related publications below