Table of Contents Author Guidelines Submit a Manuscript
Modelling and Simulation in Engineering
Volume 2008, Article ID 408075, 6 pages
Research Article

Simulation of Snow Drift and the Effects of Snow Particles on Wind

Key Laboratory of Mechanics on Western Disaster and Environment, The Ministry of Education of China, Department of Mechanics, Lanzhou University, Lanzhou 730000, Gansu, China

Received 14 June 2008; Accepted 18 December 2008

Academic Editor: Andrzej Dzielinski

Copyright © 2008 Jie Zhang and Ning Huang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Coupled equations between wind and saltating particles are presented for a stable wind blowing over an infinite plane bed and the equations are solved for a simplified particle-bed impact process. The calculated results show that the saltating snow particles strongly affect the velocity distribution of the wind, causing a deviation from a logarithmically distributed wind velocity profile. The average height and length of saltating snow particle trajectories exponentially increase as the friction velocity increases; the ejected snow number flux and the streamwise snow transport rate also increase as the friction velocity increases.

1. Introduction

Blowing snow occurs in many parts of the world and can originate from a variety of sources, including ice cover, mountain glaciers, and avalanches [1]. It also contributes to natural disasters such as ice-snow floods, avalanches, and mudslides [2] that negatively impact local environments and economies and sometimes cause the loss of human life.

In China, blowing snow is common throughout the provinces of Qinghai, Tibet, Gansu, Sichuan, and Inner Mongolia. Snow frequently accumulates on highways in these regions, causing considerable traffic delays, and the diminished visibility due to blowing snow results in hazardous driving conditions [3].

Wind-snow interactions have been studied as far back as the early twentieth century, when the characteristics of blowing snow were investigated in Antarctica [4]. Later, many other studies examined the movement and trajectory of snow particles via wind tunnel experiments, high-speed photography, and theoretical and numerical analyses (e.g., [58]). In China, extensive wind tunnel simulations of wind-snow interactions have been conducted with the aim of developing a qualitative approach to minimize the destructive effects of blowing snow [3].

Several recent theoretical analyses of saltating snow particles (such as [2, 911]) do not address how such particles affect the wind velocity distribution, thus their conclusions may have limited practical use, given that saltating particles have long been known to appreciably change the speed of the wind that drives them [6, 7]. This paper will present coupled equations between wind and saltating particles for a stable wind blowing over an infinite plane bed, and the equations will be solved for a simplified particle-bed impact process.

2. Force Analysis of a Snow Particle in Wind

High-speed photographs of snow particles in a wind tunnel [12] revealed several distinct types of motion that arise at different wind speeds. As wind speed becomes progressively higher, the following types of particle motion appear: vibrating in place, rolling along the snow surface, slipping along the snow surface, saltating (bouncing) along the snow surface, and suspension in air while moving along with the wind. When saltation is present, it is the predominant mode of particle movement [12] and so it is our main concern.

Surface grains with a dendritic shape tend to be broken by impacting grains, thus grains become small and more elastic in time and the bulk density increases [13], so people usually assumed the drifting snow grains to be spherical in their research [1416]). Figure 1 depicts the forces (not to scale) exerted on a spherical snow particle in a region, where wind flows in the positive direction and where is the particle height above the bed. and are the snow particle diameter and density, respectively.

Figure 1: Forces acting on a saltating particle.

and are the air density and gravitational acceleration, respectively; and forces acting on the particle are (1) weight ; (2) aerodynamic drag ; (3) aerodynamic lift ; (4) Magnus lift (due to particle spin) ; (5) electrical force .

The weight may be expressed as

The wind speed, , is greater than the speed of the saltating particle with respect to the snow surface, so the drag force on the particle acts in the positive direction and can be expressed as [17, 18] where is the drag coefficient and is the velocity of the particle relative to the air. is given by

For a spherical particle, the drag coefficient is given by the empirical formula [19] where Re is the Reynolds number, which may be expressed as where is the kinematic viscosity of air.

Experiments have shown that the forces , , and are typically two orders of magnitude smaller than and [17], so we will not consider them.

3. Equations of Motion and Boundary Conditions

To analyze the saltation of snow particles, we must consider (1) the acceleration of the particles by the wind, (2) the force on the wind due to saltating particles, (3) the collision of particles with the snow surface. For a stable wind blowing over an infinite plane bed, Ungar and Haff [17], Werner [20], and so on applied a modified Navier-Stokes equation: where is the horizontal force per unit volume on the wind due to the presence of saltating particles, is von Karman's constant, y is the height above the snow surface, and u is the wind speed in the X-direction.

The two dominant forces acting on a snow particle can be decomposed into x and y components to yield the equations of motion for its two-dimensional trajectory: where x and y are the horizontal and vertical coordinates, respectively, g is the acceleration due to gravity, and are the snow particle mass and diameter, and , are the and components of the net force experienced by the particle.

During the process of saltation, each particle will cross the same height y twice, once upward and once downward. Let and be equal to the number density of upward and downward particles, respectively, as a function of height. Similarly, let and be equal to the horizontal force on a single upward and downward particle. By Newton's third law, per unit volume, the horizontal force exerted on the wind must be equal and opposite to the horizontal force exerted on the particles: This equation can be rewritten as where s is the number of particles ejected per unit area of bed per unit time (assuming all particles are ejected at the same initial speed), and are the vertical speed of a single upward and downward particle, respectively.

When a saltating grain approaches the snow bed with a velocity , the number of ejected grains with is given by [17] where is the magnitude of , is the unit vector in the positive y direction, is the lift-off velocity of a vertically ejected grain, A and B are dimensionless constants. Obviously, (10) means that a saltating grain impacts the bed with a speed , then the number of ejected grains with a same speed is .

For (6), the effective bed roughness is , giving us the boundary condition , .

Assuming all grains are ejected at the same speed, (7) are solved using the boundary condition

4. Computational Method

According to Bagnold [21], the transport of grains is limited, for a given wind speed, by an “equilibrium mechanism” whereby the number of grains per unit time ejected into the saltation layer is equal to that falling onto the bed. Thus, for each grain that impacts the bed, another grain enters the saltation layer. So in (10), the value of is unity.

Based on observations at air temperatures less than −6°C, Zhonglong and Yuan [1] classified surface snow particles into three types: new snow (diameter less than 1.0 mm, average density 0.06 g/cm), fine snow (diameter less than 0.5 mm, average density less than 0.18 g/cm), and old fine snow (diameter from 0.5 mm to 1.0 mm, average density 0.23 g/cm). The experimental results with which we compare our theoretical results are limited to fine snow.

Zhonglong [3] identified a threshold wind speed for fine snow particles, at a height of 1.85 m, of 3.65 m/s; below this speed, the number of saltating grains is too small to cause a measurable decrease in wind speed. From this result, the particle lift-off speed can be calculated and then used to solve (7) to find . Then, by solving (6) and (7), we can obtain the wind profile and the snow particle trajectory.

To reach a numerical solution, (1)input values for s and for at and then solve (6),(2)substitute the result from (6) into (7) to calculate the particle trajectory and find the value of ,(3)substitute the results from (7) into (9) to calculate and then substitute into (6),(4)repeat the process 1–3 until the wind profile calculated in step (3) adequately approaches the value calculated in step (1),(5)if the precision condition ( is a positive number with a sufficient small value) is not satisfied, then the input values must be adjusted accordingly, and (6), (7), and (9) solved iteratively until the condition is satisfied.

5. Analysis of Results

For snow-wind coupling, the particle trajectory heights and lengths will increase with the friction velocity (Figure 2(a)).

Figure 2: (a) Saltation trajectories at different friction velocities. (b) (i) Photograph of particle trajectories. (ii) Trace of particle motion on the internal of each spot is 3.7 m/s.

Figures 3(a) and 3(b) compare the coupled and uncoupled particle trajectories at two different friction velocities.

Figure 3: (a) A comparison of coupled and uncoupled particle trajectories for  m/s. (b) A comparison of coupled and uncoupled particle trajectories for  m/s.

The particle trajectory is strongly influenced by coupling. When  m/s, the coupled trajectory height exceeds the uncoupled height by 48%; when  m/s, the coupled trajectory height exceeds the uncoupled height by 20%.

Kikuchi [6], Maeno et al. [7], and Zhonglong [3] observed that snow-wind coupling causes the vertical wind velocity profile to deviate from a logarithmic distribution. In Figure 4, our calculated results illustrate this deviation at different friction velocities; the effect tends to be more pronounced at higher wind speeds when the ejected particle rate s is greater.

Figure 4: Calculated wind profiles with drifting fine snow.

Figures 5(a) and 5(b) show our calculated results for coupled and uncoupled wind profiles based on the experimental data of Zhonglong [3]. Our coupled profile is qualitatively consistent with Wang's result.

Figure 5: (a) A comparison of the coupled and uncoupled wind velocity profiles. (b) The experimentally determined wind profile of Zhonglong [3].

Figure 6 compares our calculations of trajectory height versus friction velocity with Wang's experimental result (grain diameter = 0.5 mm, mean density = 0.9 g/c).

Figure 6: Trajectory height versus friction velocity.

According to Zhonglong [3], the average height of saltating snow particles increases exponentially with friction velocity. Figure 6 shows a numerical fit of our calculated values, which is given by where

Figures 7(a) and 7(b) compare our calculated values of particle trajectory length versus friction velocity with the experimental result of Maeno et al. [8], showing qualitative agreement. A numerical fit for our values is given by where

Figure 7: (a) Our results for trajectory length versus friction velocity. (b) The experimental result of Maeno et al. [8].

We have used our numerical model to calculate s, the number of ejected particles per unit area of bed per unit time, based on the experimental results of Zhonglong [3]. Figure 8 shows these values as a function of friction velocity, an exponential relationship that is given by where , , and .

Figure 8: Ejected grain number versus friction velocity.

Zhonglong [3] proposed, based on experimental results, that the saltation transport rate varies as the cube of friction velocity. Figure 9 shows our calculated values for Q, the transport rate per unit width of bed, which is linear to the cube of friction velocity for greater than 0.2 . Our result is consistent with that of Dziji [5], validating our theoretical model in this respect.

Figure 9: Saltation transport rate versus cube of friction velocity.

6. Conclusion

In this paper, coupled equations between wind and saltating snow particles are presented for a stable wind blowing over an infinite plane bed, and the equations are solved for a simplified particle-bed impact process. (1)Snow-wind coupling results in considerably different particle trajectories when compared to a model in which no coupling exists.(2)The average height and length of saltating snow particles increase exponentially with friction velocity.(3)Our model provides the number of particles ejected from a unit area of bed per unit time, given experimental data for particle trajectory height and friction velocity.(4)The transport rate of saltating particles is linear to the cube of friction velocity, in agreement with Kobayashi [5].


This research was partially supported by a grant from the National Natural Science Foundation of China (10772073, 40571018), Joint Project of NSFC and ANR (10811130470), the New Century Outstanding Talent of the Ministry of Education of China, and the Science Fund of the Ministry of Education of China for Ph.D. Program.


  1. W. Zhonglong and C. Yuan, “Research on prevention of snow-drifts by blower fences,” Journal of Glaciology, vol. 26, no. 94, pp. 435–445, 1980. View at Google Scholar
  2. J. L. Michaux, F. Naaim-Bouvet, K. Kosugi, A. Sato, and T. Sato, “Study in a climatic wind tunnel (Cryospheric Environment Simulator) of the influence of the type of grain of snow and of the flow on the formation of a snow-drif,” La Houille Blanche—Revue Internationale de l'Eau, no. 6-7, pp. 79–83, 2002. View at Google Scholar
  3. W. Zhonglong, Research on Snow Drift and Its Hazard Control Engineering in China, Lanzhou University, Lanzhou, China, 2001.
  4. G. C. Simpson, British Antarctic Expedition, 1910–1913. Meteorology. Volume 1, Thacker, Spink, Calcutta, India, 1919.
  5. K. Dziji, “Studies of snow transport in low-level drifting snow,” Low Temperature Science, vol. 12, no. 8, pp. 73–79, 1972. View at Google Scholar
  6. T. Kikuchi, “A wind tunnel study of the aerodynamic roughness associated with drifting snow,” Cold Regions Science and Technology, vol. 5, no. 2, pp. 107–118, 1981. View at Google Scholar
  7. N. Maeno, K. Araoka, K. Nishimura, and Y. Kaneda, “Physical aspects of the wind-snow interaction in blowing snow,” Journal of the Faculty of Science, Hokkaido University. Series 7, vol. 6, no. 1, pp. 127–141, 1980. View at Google Scholar
  8. N. Maeno, R. Naruse, K. Nishimura et al., “Wind-tunnel experiments on blowing snow,” Annals of Glaciology, vol. 6, pp. 63–67, 1985. View at Google Scholar
  9. R. Bintanja, “Modification of the wind speed profile caused by snowdrift: results from observations,” Quarterly Journal of the Royal Meteorological Society, vol. 127, no. 577, pp. 2417–2434, 2001. View at Publisher · View at Google Scholar
  10. R. Bintanja, “Buoyancy effects induced by drifting snow particles,” Annals of Glaciology, vol. 32, pp. 147–152, 2001. View at Google Scholar
  11. S. Alhajraf, “Computational fluid dynamic modeling of drifting particles at porous fences,” Environmental Modelling & Software, vol. 19, no. 2, pp. 163–170, 2004. View at Publisher · View at Google Scholar
  12. W. Zhonglong and X. Weirong, Study on the Evaporation of Grassland in the East Part of Qilian Mountain, vol. 3 of Memories of Glaciology and Geocryology, Institute of CAS, Science Press, Beijing, China, 1982.
  13. J. J. J. Doorschot and M. Lehning, “Equilibrium saltation: mass fluxes, aerodynamic entrainment, and dependence on grain properties,” Boundary-Layer Meteorology, vol. 104, no. 1, pp. 111–130, 2002. View at Publisher · View at Google Scholar
  14. R. A. Schmidt, “Threshold wind-speeds and elastic impact in snow transport,” Journal of Glaciology, vol. 26, no. 94, pp. 453–467, 1980. View at Google Scholar
  15. R. Bintanja, “Snowdrift suspension and atmospheric turbulence—part I: theoretical background and model description,” Boundary-Layer Meteorology, vol. 95, no. 3, pp. 343–368, 2000. View at Publisher · View at Google Scholar
  16. M. Nemoto and K. Nishimura, “Numerical simulation of snow saltation and suspension in a turbulent boundary layer,” Journal of Geophysical Research D, vol. 109, no. 18, Article ID D18206, 14 pages, 2004. View at Publisher · View at Google Scholar
  17. J. E. Ungar and P. K. Haff, “Steady state saltation in air,” Sedimentology, vol. 34, no. 2, pp. 289–299, 1987. View at Google Scholar
  18. R. S. Anderson and P. K. Haff, “Wind modification and bed response during saltation of sand in air,” Acta Mechanica, vol. 1, supplement, pp. 21–51, 1991. View at Google Scholar
  19. F. Write, Viscous Fluid Flow, Mcgraw-Hill, New York, NY, USA, 1974.
  20. B. T. Wener, “A steady-state modrl of wind-blown sand transport,” Journal of Geology, vol. 98, pp. 1–17, 1990. View at Google Scholar
  21. R. A. Bagnold, Physics of Blown Sand and Desert Dunes, Mathuen, London, UK, 1941.