各向异性密度流模型及其验证 An anisotropic buoyant turbulence model and its validation.pdfVIP

VoL 第2l卷第1期 水科学进展 21．No．1 ADVANCESIN SCIENCE Jan．．2010 2010年1月 WATER 各向异性密度流模型及其验证 陆俊卿1，张小峰2，崔占峰3 (1．环境保护部华南环境科学研究所，广东广州510655；2．武汉大学水资源与水电工程科学国家重点实验室， 湖北武汉430072；3．长江科学院河流研究所，湖北武汉430010) 摘要：针对密度流流动各向异性的特性，采用非线性五嚼模式求取雷诺应力，结合RNG模型及浮力修正建立了三 维密度流模型。应用该三维模型计算了水库温差异重流，分析了温差异重流演进过程，并就流速分布与水温下泄过 程与实测资料进行了比较。结果表明，所建模型由于引入非线性RNG模式，能有效反映密度流的各向异性特性， 较好地模拟了流场和温度场强烈耦合的水库温差异重流，能准确预测温度分层过程与异重流潜行时间，具有较高 精度。 关键词：数值模拟；密度流；各向异性；紊流模型 中图分类号：X143 文献标志码：A 文章编号：1001-6791(2010)01-0095436 由于水体密度、压力、温度和浓度的不均匀分布，天然环境水流常表现为多相、多过程、极为复杂的受 浮力影响的密度流流动。虽然对密度流流动国内外已进行了大量研究，但直到今天，用紊流模型来精细模拟 密度流动仍面临相当大的挑战。 k-e模型是目前为止得到最广泛应用和检验的紊流模型，是最适合于一般目的紊流模拟的模型…。浮力 修正的k-e模型在模拟分层流流动中也有了一定的应用旧引，但由于湍浮力流动的各向异性强烈，线性k-e 模型基于各向同性假设，Rodi虽对紊动动能方程和耗散方程作了浮力修正，但未改变线性k-8模型应力的各 向同性【3 k-e湍流模式以来，许多学者对其进行了修正和发展，并将其应用 J。自Yakhot和Orszag提出RNG k-8模 型与非线性模式相结合，分别对湍流冲击射流问题以及三维直方管内充分发展的湍流流动进行了预测，得到 的结果与标准k-8模型相比具有显著的提高…]。本文对非线性RNGk-e模型进行浮力修正建立了各向异性 的湍浮力流模型并采用实验室资料对模型进行了验证。 1 模型方程 在RNG而喏模型中，通过在大尺度运动和修正后的粘性项体现小尺度的影响，而使这些小尺度运动系 k-6模型的最大特点在于在占方程产生项的系数C。。计算中引入了主流的时均应变率S。，这样在RNGk-E模 k-8模型方 型中的C，。值不仅与流动情况相关，而且在同一问题中也还是空间坐标的函数。浮力修正的RNG 程及其参数详见参考文献[10]。本文在浮力修正的RNGk-8模型的基础上，引入非线性的雷诺应力模式建 立非线性RNGk-s模型。模型中非线性雷诺应力模式如下： (1) 一五丐=一寺觞口+_|}1尼国F+cD尸(D拥西呵一÷D。西。岛)+cez2(西F一{面。6#) 式中瓦为雷诺应力；k为紊动动能；占为紊动动能耗散率；z为混合长度尺度，z=2v。／k∽；玩为平均应变 收稿日期：2008．1l-06 基金项目：国家自然科学基金资助项目 作者简介：陆俊卿(1979一)，男，江苏南通人，工程师，博士，主要从事水环境数值模拟研究。 E—mail：lujunqing@scies．org 万方数据 水科