The random walk technique is commonly used to model diffusion in the environment. For a constant diffusivityK and model time-step @dt, the random step should be chosen from a distribution with variance [email protected] However, ifK varies spatially, this choice of step leads to the accumulation of particles in regions of low diffusivity. This problem may be overcome either by the incorporation of an apparent advection velocity, or by transforming to a coordinate system in which the diffusivity is constant. The latter technique requires no immediate approximations, is applicable to any reasonable diffusivity field and is therefore the preferred approach. In this case, as with constantK, the random step should be chosen from a normal distribution, for reasons of both theoretical accuracy and computational efficiency. Three important aspects of model design are discussed: the selection of the random number generator, the time step and the total number of particles.
[1]
D. Ramsden,et al.
Timestepping Lagrangian particles in two dimensional Eulerian flow fields
,
1991
.
[2]
William H. Press,et al.
Numerical recipes
,
1990
.
[3]
Ae Hill,et al.
Vertical migration in tidal currents
,
1991
.
[4]
G. Csanady.
Turbulent Diffusion in the Environment
,
1973
.
[5]
A. Elliott,et al.
Shear diffusion and the spread of oil in the surface layers of the North Sea
,
1986
.
[6]
Frederick Mosteller,et al.
Probability With Statistical Applications
,
1962
.
[7]
J. Church,et al.
Modelling the advection of vertically migrating shrimp larvae
,
1983
.
[8]
E. Maier‐Reimer,et al.
On the spreading of power plant cooling water in a tidal river applied to the river Elbe
,
1978
.