rtmf6 - Reactive Transport Modeling with MODFLOW 6 and PHREEQC 3¶
Status of rtmf6¶
rtmf6 is in alpha stage. It is working for test models. It is currently tested on larger cases. You can try it out by following the instructions in the documentation for rtmf6.
Please contact us for more information.
Coupling flow, transport, and chemical reactions¶
rtmf6 couples the models MODFLOW 6 and PHREEQC into a reactive groundwater transport model. It supports parallel computing both in transport and the reaction part.
Coupling the transport part - MODFLOW 6¶
A running MODFLOW model is accessed with pymf6. It uses Modflowapi and xmipy that, in turn, uses ctypes to access bmi-python, which implements the Basic Model Interface (BMI).
Coupling the reaction part - PHREEQC¶
The PHREEQC model is accessed with PhreeqPy that builds on PhreeqcRM.
Architecture¶
This schematics shows how all these parts work together.
rtmf6 is Open Source.
Workflow¶
Preprocessing¶
After create a MODFLOW 6 transport model, using a GUI tool or flopy, a preprocessing steps has to be applied. This steps clones the model by coping its input files into different sub-directories. This creates one model for each reaction component. Typically, the components O and H are explicitly model. Therefore, each will be represented by its own sub-directory. This figures shows the how ths preprocessing steps work.
Running¶
The model runs in parallel. Each transportmodel runs its own process, speeding up the execution time considerably. In addition, PhreeqcRM provides parallel execution of the chemical reaction computations.
