A two-dimensional numerical model of the saltation process was developed on
a parallel computer in order to investigate the temporal behaviour of tran
sport rate as well as its downwind distribution. Results show that the effe
cts of unsteady flow on the transportation of particulates (sediment) have
to be considered in two spatial dimensions (x, y).
Transport rate Q(x, t) appears in the transport equation for mass M(x, t).
1 partial derivative M/A partial derivative t = - partial derivative Q/part
ial derivative x + S
where A = Delta xW denotes unit area composed of unit streamwise length Del
ta x and width W. S(x, t) (units kg m(-2) s(-1)) stands for the balance ove
r the splash process. A transport equation for transport rate itself
partial derivative Q/partial derivative t = Uc partial derivative Q/partial
derivative x - Q partial derivative U-c.partial derivative x + partial der
ivative/partial derivative t (Delta xS)
is suggested with U-c (x, t) a mean particle velocity at location x as the
characteristic velocity of the grain cloud.
For a steadily blowing wind over a 50 m long sediment bed it was found that
downwind changes in Q cease after roughly 10-40 m, depending on the streng
th of the wind. The onset of stationarity (partial derivative/partial deriv
ative t = 0) was found to be a function of the friction velocity and locati
on. The local equilibrium between transport rate and wind was obtained at d
ifferent times for different downstream locations. Two time scales were fou
nd. One fast response (in the order of 1) to incipient wind and a longer ti
me for equilibrium to be reached throughout the simulation length. Transpor
t rate also has different equilibrium Values at different locations.
A series of numerical experiments was conducted to determine a propagation
speed of the grain cloud. It was found that this velocity relates linearly
to friction velocity. Copyright (C) 2000 John Wiley & Sons, Ltd.