In this paper we describe the OpenGeoSys (OGS) project, which is a scientific open source initiative for numerical simulation of thermo-hydro-mechanical-chemical (THMC) processes in porous media. The basic concept is to provide a flexible numerical framework (using primarily the Finite Element Method (FEM)) for solving multi-field problems in porous and fractured media for applications in geoscience and hydrology. To this purpose OGS is based on an object-oriented FEM concept including a broad spectrum of interfaces for pre-and post-processing. The OGS idea has been in development since the mid eighties. We provide a short historical note about the continuous process of concept and software development having evolved through Fortran, C, and C++ implementations. The idea behind OGS is to provide an open platform to the community, outfitted with professional software engineering tools such as platform-independent compiling and automated benchmarking. A comprehensive benchmarking book has been prepared for publication. Benchmarking has been proven to be a valuable tool for cooperation between different developer teams, e.g. for code comparison and validation purposes (DEVOVALEX and CO2 BENCH projects). On one hand, object-orientation (OO) provides a suitable framework for distributed code development; however the parallelization of OO codes still lacks efficiency. High-performance-computin (HPC) efficiency of OO codes is subject to future research.

No abstract

No abstract

Abstract. Most large-scale hydrologic models fall short in reproducing groundwater head dynamics and simulating transport process due to their oversimplified representation of groundwater flow. In this study, we aim to extend the applicability of the mesoscale Hydrologic Model (mHM v5.7) to subsurface hydrology by coupling it with the porous media simulator OpenGeoSys (OGS). The two models are one-way coupled through model interfaces GIS2FEM and RIV2FEM, by which the grid-based fluxes of groundwater recharge and the river–groundwater exchange generated by mHM are converted to fixed-flux boundary conditions of the groundwater model OGS. Specifically, the grid-based vertical reservoirs in mHM are completely preserved for the estimation of land-surface fluxes, while OGS acts as a plug-in to the original mHM modeling framework for groundwater flow and transport modeling. The applicability of the coupled model (mHM–OGS v1.0) is evaluated by a case study in the central European mesoscale river basin – Nägelstedt. Different time steps, i.e., daily in mHM and monthly in OGS, are used to account for fast surface flow and slow groundwater flow. Model calibration is conducted following a two-step procedure using discharge for mHM and long-term mean of groundwater head measurements for OGS. Based on the model summary statistics, namely the Nash–Sutcliffe model efficiency (NSE), the mean absolute error (MAE), and the interquartile range error (QRE), the coupled model is able to satisfactorily represent the dynamics of discharge and groundwater heads at several locations across the study basin. Our exemplary calculations show that the one-way coupled model can take advantage of the spatially explicit modeling capabilities of surface and groundwater hydrologic models and provide an adequate representation of the spatiotemporal behaviors of groundwater storage and heads, thus making it a valuable tool for addressing water resources and management problems.

Scientific visualization is an integral part of the modeling workflow, enabling researchers to understand complex or large data sets and simulation results. A highresolution stereoscopic virtual reality (VR) environment further enhances the possibilities of visualization. Such an environment also allows collaboration in work groups including people of different backgrounds and to present results of research projects to stakeholders or the public. The requirements for the computing equipment driving the VR environment demand specialized software applications which can be run in a parallel fashion on a set of interconnected machines. Another challenge is to devise a useful data workflow from source data sets onto the display system. Therefore, we develop software applications like the OpenGeoSys Data Explorer, custom data conversion tools for established visualization packages such as ParaView and Visualization Toolkit as well as presentation and interaction techniques for 3D applications like Unity. We demonstrate our workflow by presenting visualization results for case studies from a broad range of applications. An outlook on how visualization techniques can be deeply integrated into the simulation process is given and future technical improvements such as a simplified hardware setup are outlined.