snowBedFoam: an OpenFOAM Eulerian-Lagrangian solver for modelling snow transport

snowBedFoam 1.0. is a snow transport solver implemented in the computational fluid dynamics software OpenFOAM. It is adapted from the standard multi-phase flow solver DPMFoam for application in snow-influenced environments. To simulate aeolian snow transport, snowBedFoam 1.0. handles coupled Eulerian–Lagrangian phases, which involve a finite number of particles (snow) spread in a continuous phase (air). The snow erosion and deposition are modelled through physics-based equations similar to the ones employed in the well-established LES-Lagrangian Stochastic Model (Comola and Lehning, 2017 ; Sharma et al., 2018 ; Melo et al., 2022). This modelling approach is computationally intensive and thus adapted to simulate snow movement and distribution on small scale terrain. First, snowBedFoam 1.0. was applied to topographical data collected on Arctic sea ice during the MOSAiC expedition (Clemens-Sewall, 2021). Together with atmospheric data from the MOSAiC Met City (Shupe et al., 2021) used for the fluid forcing, the model was able to accurately simulate the zones of erosion and deposition of snow along a complex ice ridge structure (Hames et al., 2022). Second, snowBedFoam 1.0. was used to simulate the snow distribution around the German Antarctic research station Neumayer Station III. The effect of snow properties, fluid forcing and aerodynamic structures on the snow accumulation were assessed. snowBedFoam 1.0 was implemented in 2 different OpenFOAM versions, namely OpenFOAM-2.3.0 and OpenFOAM-5.0. The latter offers more options for turbulence models and boundary conditions. The fundamental model equations were not changed from one implementation to the other, thus both still correspond to snowBedFoam 1.0. The two branches are called snowBedFoam-v1-2.3.0 (OpenFOAM-2.3.0) and snowBedFoam-v1-5.0 (OpenFOAM-5.0). The core codes of snowBedFoam 1.0. are directly accessible on the WSL/SLF GitLab repository (more details in the Resources section).