CFD-DEM modeling of snowdrifts on stepped flat roofs

Snowdrift formation on roofs should be considered in snowy and windy areas to ensure the safety of buildings. Presently, the prediction of snowdrifts on roofs relies heavily on field measurements, wind tunnel tests and numerical simulations. In this paper, a new snowdrift modeling method by using CFD (Computational Fluid Dynamics) coupled with DEM (Discrete Element Method) is presented, including material parameters and particle size, collision parameters, particle numbers and input modes, boundary conditions of CFD, simulation time and inlet velocity, and coupling calculation process. Not only is the two-way coupling between wind and snow particles which includes the transient changes in snow surface topography, but also the cohesion and collision between snow particles are taken into account. The numerical method is applied to simulate the snowdrift on a typical stepped flat roof. The feasibility of using coupled CFD with DEM to study snowdrift is verified by comparing the simulation results with field measurement results on the snow depth distribution of the lower roof.