Premium
Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones
Author(s) -
Li Bing,
Liu Xiaofeng,
Kaufman Matthew H.,
Turetcaia Anna,
Chen Xingyuan,
Cardenas M. Bayani
Publication year - 2020
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2019wr026528
Subject(s) - multiphysics , modular design , coupling (piping) , flow (mathematics) , computer science , flux (metallurgy) , interface (matter) , flexibility (engineering) , boundary value problem , computational science , mechanics , finite element method , engineering , chemistry , mechanical engineering , mathematics , physics , mathematical analysis , statistics , structural engineering , organic chemistry , bubble , maximum bubble pressure method , parallel computing , operating system
Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one‐way coupling between the surface and subsurface domains. Such one‐way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam , was developed using the open‐source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod‐type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three‐dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.