k-ε双方程水库水温模型浮力模拟分析 - 中国水利水电科学研究院学报.pdf

第15卷 第1期 中国水利水电科学研究院学报 Vol.15 No.1 2017年2月 Journal of China Institute of Water Resources and Hydropower Research February，2017 文章编号：1672-3031（2017）01-0037-07 k-ε双方程水库水温模型浮力模拟分析 林 玲，梁瑞峰，李克锋，李 嘉 （四川大学 水力学与山区河流开发保护国家重点实验室，四川 成都 610065） 摘要：在比较国内外常用水库水温紊流模型的基础上，针对采用k-ε双方程紊流模型的宽度平均立面二维水库水温 模型，从理论分析和数值模拟两方面讨论了k-ε紊流模型在模拟温差异重流时紊动动能生成项和浮力项的作用机 理。分析了紊动动能生成项和浮力项对紊动动能和耗散率各自的影响程度，对水温垂向梯度与紊动动能的相互制约 关系进行了探讨，研究了紊动动能浮力项和生成项在水库垂向的分布规律。W2模型根据实测数据率定的涡黏系数 或混合长度在通用性上存在一定不足，但在使用者经验丰富的条件下仍能得到足够精确的预测结果。k-ε模型紊动 动能生成项和浮力项引起的耗散率改变对紊动动能的影响作为小量可予忽略；k-ε模型可通过紊动动能生成项和浮 力项计算与垂向流速和温度分布相适应的涡黏系数，将温差浮力影响反映在垂向动量方程中；浮力项主要受垂向温 度梯度控制，紊动动能与垂向温度梯度的相互制约关系是k-ε紊流模型模拟温差异重流的关键因素。 关键词：紊流模型；紊动动能；浮力项；温差异重流 中图分类号：TV131.3 文献标识码：A doi ：10.13244/ki.jiwhr.2017.01.005 水库水温分层对水库水质、水生动植物及下游水生环境都具有重大影响。水温分层环境下的水 库一般表层温度与气温变化同步，而具有较稳定的库底水温、一定的表层与底层温差和温跃层，温 跃层范围内的水温在垂直方向上存在剧烈变化。这种环境下的水体流动属于湍浮力流，表现在垂向 上由于温度差产生浮力，出现分层流动等异于等密度流的特殊流动形式。伴随这种流动，水库垂向 出现水温分层形成、发展和消失的过程。水库的流场分布与水库的温度分布具有强烈的相互作用， 因此准确模拟水库温度场的关键在于精确模拟浮力流流场。 k-ε紊流双方程模型是在单方程模型的基础上，增加确定紊动特征长度的微分输运方程而得到 的，因而这种模型仍然采用了各向同性的紊动黏性系数的概念。在双方程模型中，标准k-ε模型在工 程实际中得到了最广泛的应用和验证，它考虑了紊动速度比尺和紊动长度比尺的输运，对于大多数 ［1-3］ 水流问题，标准k-ε双方程模型或其改进模型一般能得到较好的结果 。 标准k-ε模型也有其不足，如模型中经验常数的通用性尚不十分令人满意，另外，对于有些各向 异性明显的水流或流动区域，有必要精确地描述应力各分量的输运，各向同性涡黏性系数的概念和 据此建立起来的标准k-ε模型的适用性受到怀疑。有学者尝试采用各向异性的紊动黏性系数，如 ［4］ Young 在其二维水库模型中对垂向上的动量和热量的紊动扩散采用了不同经验公式。 ［5］ 也有学者直接采用各向同性的k-ε模型计算各向异性紊流流动。Farrell等 将k-ε模型应用于水 库密度流的模拟，试验性地计算了一个100m长的水库的下潜流过程，结果认为k-ε模型能够模拟出 ［6-7］ 水库密度流的下潜、垂向旋涡和温度分层的特征现象。余真真 采用k-ε模型计算了温度分层环境 ［8］ ［9］ 对潮成内波和悬浮颗粒垂向输运的影响。邓云等 、刘兰芬等 也对水库温差异重流进行了模拟并得 ［8-9］ 到了实测资料的验证 。 综上可知，k-ε模型能够模拟某些