粘弹性减阻湍流模型研究.pdfVIP

中国工程热物理学会 传热传质 学术会议论文 编号：083026 粘弹性减阻湍流模型研究 1 1* 1 2 3 3 4 王艺 宇波 张劲军 魏进家 张红娜 李凤臣 川口靖夫 ( 1. 中国石油大学石油天然气工程学院，城市油气输配技术北京市重点实验室，北京 102249; 2. 西安交通大学能源与动力工程学院，陕西 西安 710049; 3. 哈尔滨工业大学能源科学与工程学院，黑龙江 哈尔滨 150001; 4. 东京理科大学机械工程系，日本 千叶县 2788510 ) *通信作者：电话：010 E-mail:yubobox@cup.edu.cn 摘要：基于 Giesekus 方程，建立了粘弹性减阻湍流的低 Re 数 k-ε模型和低 Re 数雷诺应力模型。与 DNS 数据比较，发现这两个模型对平均速度、平均变形和减阻率的预测均较准确，而对脉动速度强度的预 测精度有待提高。雷诺应力模型预测的雷诺应力、粘弹性应力和脉动强度等比 k-ε模型更接近 DNS ， 体现出了一定的优势。 关键词：k-ε模型，雷诺应力模型，湍流模型，减阻，粘弹性 0 前言 近年来直接数值模拟(DNS)研究和流体测量实验表明高分子和表面活性剂等添加剂 引起的湍流减阻与溶液的粘弹性有关。DNS 在探索添加剂引起的粘弹性减阻机理方面已 取得了丰硕的成果[1～4] ，大量的 DNS 数据库已建立。虽然 DNS 日益成为研究添加剂减 阻机理的一项基本手段，但是它需要非常大的计算机内存容量和机时耗费，目前只能对 中等以下雷诺数的流动进行研究，还远远不能满足工程需要。为切实可行地解决工程实 际问题，有必要吸收 DNS 和实验测量的最新研究成果，发展求解湍流时均值的粘弹性 湍流模型。众多湍流模型中 k-ε模型多年来得到了广泛应用，而在粘弹性减阻流体的 k-ε 模型方面有代表性的是Pinho 研究小组开发的低Re 数k-ε模型及其改进的模型[5~6] 。Pinho 的模型考虑了拉伸稠化这一减阻发生的重要因素，得到了与实验吻合良好的平均流速和 摩擦阻力系数，有力推动了减阻湍流模型的发展。另外，雷诺应力模型(RSM) 以其更好 的准确性和适应性，越来越受到人们的重视。Richard Leighton[7]等在使用 FENE-P 本构 方程模拟粘弹性减阻流动时，建立了低 Re 数雷诺应力模型，得到与 DNS 数据吻合较好 的平均速度和平均变形等。本文基于粘弹性本构关系中的Giesekus 方程，建立了减阻流 体的低 Re 数 k-ε模型和低 Re 数雷诺应力模型，并与 DNS 数据进行对比，以期找到更适 于粘弹性减阻湍流的模型构建方式。 *基金项目：教育部科学技术研究重点项目(106034)、高等学校博士学科点专项科研基金(20050425006)、 国家自然科学基金50521604) 。 1 控制方程及其时均化 [8] 采用 Giesekus 模型模拟粘弹性减阻槽道湍流的无量纲控制方程组 如下： ∂u+ i 0 (1) * ∂x i + + + 2 + + ∂ui + ∂ui ∂p β ∂ ui 1−β ∂cik +uk =− + +