A two-dimensional, finite-element calculation was performed to simulate the development of a snowdrift. The underlying two-dimensional topography was taken as the lower boundary, and this changed as the drift developed. The two-dimensional wind field was calculated from the Navier–Stokes equations. The calculation was performed until a quasi-steady state was reached, and the pattern of erosion and deposition was estimated by computing the divergence of snowdrift transport. Considering a snowdrift as a new surface, the wind field was then recalculated. Simulated results were compared with outdoor observations of the wind field and snowdrifts. The results were in fairly good agreement with each other, but more detailed research is regarded as being necessary.