Hydrology solvers for OpenFOAM

The OpenFOAM user guide can be found on https://www.openfoam.com/documentation/user-guide/

• Author: Laurent Orgogozo
• Contact: Laurent.Orgogozo (@) get (.) omp (.) eu

Introduction

Source code

Licence: GPL-3.0-or-later

Installation

0. You need an installation of OpenFOAM on your computer and compiler/make tools.

   • Tested with OpenFOAM-v1912 and OpenFOAM-v2106.

1. Clone the repository or download a tar-archive to a user-writeable location

2. Run the top-level ./Allwmake script, which will compile the hydrology library and the solvers.

The hydrology library primarily supplies boundary conditions such as rainFlux and noRainFlux as well as some additional functions for seasonal time variation of boundary conditions:(seasonal, thawing, cosine1, sine1, square1).

Testing

To test the solver permaFoam (see Orgogozo et al., Permafrost and Periglacial Processes 2019 for details) and the solver RichardsFoam3 (see Orgogozo et al., 2014 and Orgogozo 2015 in Computer Physics Communications for details) there are a few demonstration cases. You can copy any of these cases to the "run/" directory of your OpenFOAM workspace (create it if it does not already exist), or simply run inplace. The usual ./Allrun and ./Allclean commands are used to run the case (create mesh, run solver, postProcess).

The post-processing output files will contain averaged pressure heads, temperature, total water content, liquid water content, and ice content (in postProcessing/internal/0/volFieldValue.dat), water fluxes at the top boundary (in postProcessing/top/0/surfaceFieldValue.dat). Besides, of course, all the output fields at each time of writing of the results will be contained in the time directories (named along the associated time coordinate).

Note: The directory expected_postProcessing contains the expected results of the postprocessing operations (permaFoam -postProcess), so that they can be compared to the ones that the execution of the commandsRun procedure provides, for validating your installation.

The full fields can be visualized in paraview, with the regular ParaView/VTK OpenFOAM reader and/or with the so-called paraFoam reader module.

tutorials/permaFoam/demoCase

This case is a simple setup of simulation of permafrost thawing.

tutorials/permaFoam/demoCase_sinusoidalClimatForcings

This case illustrates the handling of seasonal variabilities in the boundary conditions.

tutorials/RichardsFoam/demoCase

This case is a simple setup of simulation of rain water infiltration.

tutorials/RichardsFoam/demoCase_sinusoidalClimatForcings

This case illustrates the handling of seasonal variabilities in the boundary conditions.

Further References

The system of equations, the input fields and the post-processing information are further described in the permaFoam wiki and the RichardsFoam wiki.

---