snowpackBuoyantPimpleFoam is a two-phase solver implemented to model convection of water vapor with phase change in snowpacks. This new solver is based on the standard solver of buoyantPimpleFoam in OpenFOAM 5.
This python code uses the Finite Element library FENICS (via docker) to solve the one dimensional partial differential equations for heat and mass transfer in snow. The results are written in vtk format from which the paper figure is reproduced.
StreamFlow is an extension to the spatially-distributed snow model Alpine3D which allows the user to perform hydrological simulations. On top of discharge and water height, StreamFlow can also compute stream temperature at any point along the stream.