Numerical modelling of multiphase multicomponent reactive transport in the Earth's interior

Beñat Oliveira, Juan Carlos Afonso*, Sergio Zlotnik, Pedro Diez

*Corresponding author for this work

    Research output: Contribution to journalArticlepeer-review

    14 Citations (Scopus)


    We present a conceptual and numerical approach to model processes in the Earth's interior that involve multiple phases that simultaneously interact thermally, mechanically and chemically. The approach is truly multiphase in the sense that each dynamic phase is explicitly modelled with an individual set of mass, momentum, energy and chemical mass balance equations coupled via interfacial interaction terms. It is also truly multicomponent in the sense that the compositions of the system and its constituent phases are expressed by a full set of fundamental chemical components (e.g. SiO2, Al2O3, MgO, etc.) rather than proxies. These chemical components evolve, react with and partition into different phases according to an internally consistent thermodynamic model. We combine concepts from Ensemble Averaging and Classical Irreversible Thermodynamics to obtain sets of macroscopic balance equations that describe the evolution of systems governed by multiphase multicomponent reactive transport (MPMCRT). Equilibrium mineral assemblages, their compositions and physical properties, and closure relations for the balance equations are obtained via a 'dynamic' Gibbs free-energy minimization procedure (i.e. minimizations are performed on-the-fly as needed by the simulation). Surface tension and surface energy contributions to the dynamics and energetics of the system are taken into account. We show how complex rheologies, that is, visco-elasto-plastic, and/or different interfacial models can be incorporated into ourMPMCRT ensemble-averaged formulation. The resulting model provides a reliable platform to study the dynamics and nonlinear feedbacks of MPMCRT systems of different nature and scales, as well as to make realistic comparisons with both geophysical and geochemical data sets. Several numerical examples are presented to illustrate the benefits and limitations of the model.

    Original languageEnglish
    Pages (from-to)345-388
    Number of pages44
    JournalGeophysical Journal International
    Issue number1
    Publication statusPublished - 1 Jan 2018


    • composition of the planets
    • dynamics of lithosphere and mantle
    • numerical modelling
    • magma genesis and partial melting
    • mechanics, theory, and modelling


    Dive into the research topics of 'Numerical modelling of multiphase multicomponent reactive transport in the Earth's interior'. Together they form a unique fingerprint.

    Cite this