用分块耦合技术模拟三峡永久船闸反向弧形阀周围的湍流场.pdfVIP

用分块耦合技术模拟三峡永久船闸 反向弧形阀周围的湍流场① 陈文学 谢省宗 (中国水利水电科学研究院水力学所) 藏光清 高旭东 (四川联合大学高速水力学国家重点实验室) 摘要 本文讨论了如何用分块耦合技术生成网格疏密状况和正交性能较好的网格，并用标准K—e模 式模拟了三峡永久船闸反向弧形阀周围的紊流场。激光测试结果表明，计算结果与激光测速值吻合得 。 较好。 关键词分块耦合，数值模拟，湍流流场 1、前言 三峡双线连续五级船闸是三峡水利枢纽三大主体工程之～，也是世界上最大的内河 通航建筑物，其总水头为113m，其中中间级船闸的作用水头达45．2m。为保证三峡船闸 的安全运行，如防止空化、空蚀和“声振”现象等，国内许多研究单位对三峡船闸的结 构和水力学特性进行了深入的分析和研究，得出了许多有价值的结论11】【2】【3114】，但阀门段 湍流场的定量研究分析尚不多见。 · 从流动特征来看，船闸阀门段的流动属动边界非恒定流动，研究难度较大。此外， 问题的复杂性还体现在：(1)阀后有不同类型的突扩，如顶扩、底扩、侧扩等，流体的 运动边界比较复杂，当计及阀门门体的厚度时．底缘区出现奇点，计算难度增加；(2)阀 门和门楣之间有一极小的缝隙，缝隙射流使阀后出流的流态更加复杂，若将此窄缝考虑 到全场数值模拟中，则给计算网格的设计(尤其是网格比的调整)带来了一定的困难。 本文用分块耦合技术将计算区域分成若干子区，子区之间或相互搭接、或由纯连接区连 2D方法生成，网格的疏密易于 接起来，各子区的网格由TM方法【5】【研和Hilgenstockf7l[SJ 控制，网格的正交性也比较理想。然后用混合有限分析法【9】成功地模拟了反弧形阀门周 围的定常湍流场。 ①国家自然科学基金委员会和中国长江三峡工程开发总公司联合资助的国家自然科学基金重大项目资 助(批准文号；5舛936002--2) 382 2、网格生成技术 工程中的流动边界往往比较复杂，为了保证阏格质量和提高计算精度，最好采用分 块耦合方法，即根据流动边界的特点，将流场空间分成若干块，分别生成计算网格／／o／／“】。 本文将船闸反弧形阀门周围的流场空间分为6个子区(见图1)，根据子区的结构特点， 可将子区分为三种类型：(1)单独生成的子区，该子区的边界一般比较复杂，如子区2 和4；(2)与别的区域有搭接部分的区域，该子区的边界比较简单，其网格的生成与相 邻的子区有关，如子区1和3：(3)纯连接块，它将左右或上下子区连接起来，如子区5 将子区2和3连接起来，子区6将子区2和4连接起来。各子区独立编号，分别求解， 搭接区和连接区起传递块与块之间信息的作用。在计算时应注意连接区、搭接区与相邻 块间的压强传递，以保证压强场的连续。 各子区的网格生成中，应注意两个问题，一是网格的疏密控制，其二是网格的正交 ● 方法在控制网格疏密时需要较强的技巧，给应用带了诸多不便：T＆M方法可以利用边界 网格节点的信息自动调整计算区域内的网格分布，但网格的正交状况不太理想： 修正源项，从而生成疏密情况和正交性比较理想的网格。 对于二维问题，物理平面上的网格节点坐标由以下控制方程确定“1 a(k+螟)一2羼。+，(％+惦)=0 【1) · F=(x，Y) 口=Rf2， ∥=弓‘， y=睡12 (2)