An efficient 3D particle transport model for use in stratified flow

Random-walk models are a versatile tool for modelling dispersion of both passive and active tracers in turbulent flow. The physical and mathematical foundations of stochastic Lagrangian models of turbulent diffusion have become more and more solid over the years. An important aspect of these types of models that has not received much attention is the behaviour of the particles near boundaries. Often, a simple stochastic, numerical scheme is used. Because turbulent mixing in the vertical direction is much more complicated than in the two horizontal directions, it is in the vertical direction that a simple numerical scheme, such as the Euler scheme, may cause problems. In this paper our main goal is the development of an efficient 3D particle transport model that can be used in stratified flow. For this type of situation the vertical direction is of special interest. First, a closer look is taken at some considerations that should be regarded when choosing a numerical scheme. Specifically schemes are investigated that can be used in the vertical direction, where the diffusion coefficient is varying in that direction. Experiments are performed regarding the accuracy of different numerical schemes in various situations. The behaviour of the particles near an impermeable layer interface is investigated. The stochastic Heun and Runge–Kutta schemes turn out to be very attractive for this type of model. For the simulation of the transport of various physical quantities, such as salinity, heat, silt, oxygen, or bacteria, different types of models are available. In this case we will take a closer look at the modelling of the transport of pollutants from point sources (either instantaneous or continuous transport). For this purpose a 3D particle transport model has been developed that is especially suited for stratified situations such as can be found in estuaries. The main idea is to use a simple numerical scheme for the horizontal directions and a higher-order method for the vertical direction. The results play an important role in making specific choices for this type of particle transport model. Copyright © 2005 John Wiley & Sons, Ltd.

[1]  N. Dinar,et al.  A three-dimensional stochastic model for concentration fluctuation statistics inisotropic homogeneous turbulence , 1988 .

[2]  P. K. Yeung,et al.  Conditional Fluid-Particle Accelerations in Turbulence , 1998 .

[3]  W. Rüemelin Numerical Treatment of Stochastic Differential Equations , 1982 .

[4]  Paul A. Durbin,et al.  A stochastic model of two-particle dispersion and concentration fluctuations in homogeneous turbulence , 1980, Journal of Fluid Mechanics.

[5]  David J. Thomson,et al.  A stochastic model for the motion of particle pairs in isotropic high-Reynolds-number turbulence, and its application to the problem of concentration variance , 1990, Journal of Fluid Mechanics.

[6]  Recent developments in the Lagrangian stochastic theory of turbulent dispersion , 1993 .

[7]  Arnold W. Heemink,et al.  Stochastic modelling of dispersion in shallow water , 1990 .

[8]  D. Thomson Criteria for the selection of stochastic models of particle trajectories in turbulent flows , 1987, Journal of Fluid Mechanics.

[9]  Geoffrey Ingram Taylor,et al.  Diffusion by Continuous Movements , 1922 .

[10]  G. Stelling,et al.  On the construction of computational methods for shallow water flow problems , 1983 .

[11]  A. Jazwinski Stochastic Processes and Filtering Theory , 1970 .

[12]  P. Kloeden,et al.  Numerical Solution of Stochastic Differential Equations , 1992 .

[13]  S. Pope PDF methods for turbulent reactive flows , 1985 .

[14]  G. Batchelor Diffusion in a field of homogeneous turbulence , 1952, Mathematical Proceedings of the Cambridge Philosophical Society.

[15]  D. W. Dunsbergen Particle models for transport in three-dimensional shallow water flow , 1994 .

[16]  H. C. Rodean Stochastic Lagrangian Models of Turbulent Diffusion , 1996 .