Numerical simulation of three-dimensional, transient snow drifting around a cube

This paper presents the results of the numerical simulation of drifting snow surrounding a simple 2 m cubicle structure. These results are compared and verified against data from snow drifting experimental tests carried out at SANAE TV research station, Antarctica, during the summer research period of January 2002. In addition to the snow drifting field observations, wind profile data were also obtained from cup-type anemometers mounted on a 6 m wind mast. These data were employed to derive the characteristic surface roughness, shear velocity and approaching wind profile functions. The present work numerically simulates a transient three-dimensional turbulent viscous flow in an Eulerian coordinate frame including snow advection. A modified turbulent wall law is employed that accounts for the effects of snow saltation on the effective aerodynamic surface roughness. The numerical simulation employed the commercial CFD code, FLOW-3D, with additional user Fortran coding added to model the snow entrainment, subsequent accumulation or erosion of snow as well as temporal snow surface changes. The snow accumulation and erosion model is presented and discussed. The snow accumulation predicted by this numerical simulation compares favourably with the experimental results obtained from the Antarctica field testing. (C) 2004 Elsevier Ltd. All rights reserved.