大气湍流扩散随机游走模拟理论探讨.pdfVIP

大气湍流扩散的随机游走模拟理论探讨 林官明 北京大学环境学院(100871) 环境模拟与污染控制国家重点实验室 email: gmlin@ 摘 要：探讨了湍流扩散随机游走模拟的理论基础，从湍流脉动动量方程出发导出均匀湍流 的 Langevin 方程，进一步得到非均匀湍流马尔可夫链形式的随机游走方程，非均匀湍流项 可以自然地出现在方程中。分析指出，可作为随机项加以模拟的是脉动压力梯度项和粘性项。 关键词：湍流扩散 随机游走 Langevin 方程 马尔可夫链 1．引言 随机游走方法最先应用于计算中子链式反应，后来在大气污染扩散研究方面也得到广 泛应用。Smith (1968)首先建议用 马尔可夫随机游走方法模拟大气扩 散 问题[1] ， Thomson (1971)首先将此法用于边界层扩散[2]，之后在大气扩散研究领域有大量文献涌现 [3][4][5]。 随机游走模拟的基本思想是用随机过程的实现来模拟湍流脉动。在大气扩散研究中，风 速的平均运动部分由其它模式提供，如应用最小变分原理的诊断模式或湍流模式等。在构造 的湍流场中释放示踪粒子，粒子在其中如流体质点一样被动运动，统计粒子的运动轨迹即得 到湍流扩散的规律。本文直接从湍流脉动方程出发，通过分析方程中各项的特征，得到用随 机过程模拟湍流流动的理论基础。 2．脉动方程量阶分析 在直角坐标系中，湍流运动的脉动动量方程为 ∂u / ∂t +(u +u )u +u u =−p / ρ+νu +u u , （1） ′ ′ ′ ′ ′ ′ ′ ′ i j j i ,j j i ,j ,i i ,jj i j j 或 / =− − / + + , du dt u u p ρ νu u u （2） ′ ′ ′ ′ ′ ′ i j i i i jj i j j ,j , , 用雷诺平均法则考虑流体质点的脉动情况时，实际上已经假设了时间和空间上两种尺度 的存在：空间上的两个特征尺度是流体团的直径D 和平均运算的几何尺度D ；时间上的特征 p a 尺度也有两个，分别为流体团的特征时间尺度T 和平均运算的特征时间T 。这样做对应于雷 p a 诺平均时人为地对速度进行平均运算，即只要进行了雷诺平均，必定以相应的时间和空间特 - 1 - 征尺度为参照，这是雷诺平均运算的物理意义。同样，由于雷诺平均决定了特征时间和空间 的存在，在湍流研究中，凡应用到雷诺平均的地方都应当考虑这一特征时间和空间的存在， 对它们的忽视可能会失去物理上的意义。无论是物理模拟还是数值模拟，其分析与运