槽道湍流涡粘张量模式的数值实验的研究.pdfVIP

槽道湍流涡粘张量模式的数值实验研究。 朱贞卫姜楠王振东 (天津大学力学系，天津3哪2) 摘要本文应用槽道湍流直接数值模拟的结果，对文献[1]提出的涡牯系数为四阶张量的涡粘张量模式进 行了数值计算．得出完全发展槽道湍流涡粘张量分量的分布规律，并对其作了简要分析． 关■词槽道湍流；湍流模式；涡粘张量 0前 言 湍流模式研究的目标是准确预测湍流的平均流动，其方向是寻找实用、准确的模式，比较简便地 预测复杂湍流的运动．文献[1]在分析工程中应用较广的原有涡粘模式中对湍流应力各向同性假设不 合理性的基础上，提出了能体现湍流应力各向异性的涡粘张量模式，并引用近代代数学中广义逆矩阵 理论，导出了涡粘张量分量的计算公式．本文在此基础上应用充分发展槽道湍流直接数值模拟的数值 结果，计算出了涡粘张量分量的分布规律，并对其进行了分析研究． 1涡粘张量模式简述 1877年，B0u鹕ir目q将湍流应力类比于与速度梯度成正比的粘性应力，引进了涡粘系数的概念，将 湍流应力表示为 一面=e。(差+鬻) ㈤ 式中涡粘系数￡。为标量．因式(1)收缩后与不可压缩流体的连续性方程矛盾．后来将其修正为 一“竹。8ml瓦+瓦J一了。庐一面：em(券+爰)一号啦 cz)¨’ 式中x是湍能． 文献[1]由张量识别定理断定，涡粘系数应是四阶张量[2： 一丽=s唰【号(鬻+鬻)】一号掣 c。， 由于湍流应力张量和变形速度张量的对称性，涡粘张量的81个分量中只有36个分量是独立的． 对于本文所作的槽道湍流，独立的涡粘张量分量只有4个不为零． 应用广义逆矩阵理论，将变形速度张量的9个分量作为一非零向量卢，可导出涡粘张量分量的计 算公式 (4) (q，，em，e蛐，s舛，，e弘，e拙，e叫，s僻，e脚)=(一面+号8乒)·一+ * 国家自然科学基金资助项目(1螂撇) ——百■■可酉——————一 l f au】 a如＼ l faulau2＼ au2 p iI瓦+瓦』’il瓦+瓦J’瓦’ 舭引 ¨亟嵋 叫一毗一‰ V、卜 帅“ 悒亟‰ + “黝坐岫引 上。坐嘞 一当‰厩一嵋 一塑‰一垫‰ 2槽道湍流直接数值模拟方法简介 本文使用清华大学工程力学系许春晓博士的槽道湍流直接数值模拟数据库，该槽道湍流是限制 在上下壁面之间，水平方向延伸至无穷的流动．在流向和展向采用周期边界条件，在横向采用固壁无 滑移条件．对于完全发展的槽道湍流，在流向和展向可认为是均匀的㈦． 物理空间的一个流动物理量q(x，y，z，t)，利用下式进行谱展开 ， V， = m Ⅱ ，箜 似 (5) P p 肛 一 岂一 黾gt ，∑脚 巧 式．配置点取为