用能耗极值原理模拟水力机械中含沙水流的数学模型.pdfVIP

中国工程热物理学会 987037 工程热力学与能源利用 学术会议 用能耗极值原理模拟水力机械中含沙水流的数学模型 李仁年 张建华 (甘肃工业大学 兰州 930050) 摘 要 本文从连续介质流体的能量方程出发，推导了含沙水流中固液两相的能量方程及 能量耗散项，并将能量耗散极值原理引入水力机械中含沙水流的研究中，从中得NT～ 些新的认识，提出了～种全新的数值模型。 关键词 水力机械 过流部件含沙水流 能量耗散极值原理 数值模拟 1． 引言 水力机械中含沙水流流动现象在我国普遍存在，其中含沙水流对过流部件的磨蚀 破坏作用己成为人们尤为关心的问题，许多学者在这一领域进行了卓有成效的工作，得 出了一些有益的结论。但由于含沙水流流动现象的复杂性，既包括了水力机械内部含沙 水流两相流动分析问题，又包括了沙粒对水力机械过流部件表面的破坏过程，使得这些 结论都有很大的局限性。近年来，在多泥沙水流情况下导叶表面形状、叶片的型线如何 适合沙粒群引起的流场改变，在研究水轮机的抗磨蚀过程中，越来越引起人们的重视， 而解决这～问题的前提是弄清过流部件中含沙水流流场情况。随着计算机的发展，用数 值模拟的方法来研究水力机械内部的含沙水流流动现象取得了很大的进展，众多学者提 出了～些不同的数学模型”】，主要有：单颗粒动力学模型、两流体模型和颗粒轨道模型， 在不同的范围内计算得各自的结果，但是人们对这些基本两相流模型还存在争议，主要 是如何确定两相间的相互作用以及变形剪切应力r，这些模型基本方程都是连续性方 程和运动方程，这样做的结果使得最后的模型方程较复杂。本文从能量变化的角度出 发，运用能量耗散极值原理，提出了一种新的数值模拟模型。 2． 能量耗散极值原理 1978年，黄万里在他的论文中首次提出能量耗散极值原理口j，即：凡属固体、液 体或其它临界物体连续介质独立系统的整体，以给定的初始和边界条件，在外力作用下 运动时，其所有质点将随时按整体机械能产生最大的热能转化率，即消散率，形成其应 力场和应变场或流速场、压力场和密度场。把这个原理应用于含沙水流的研究中，则可 芝 进一步阐述为：在同时存在浓度分布和流速分布的含沙水流中，水和沙粒的速度是在系 统的实际能量消散率较虚拟的而言总是每时每刻达到极值的条件下分布的。应用此原 毒 理，倪晋仁p1计算了一维明槽流的两相流速分布，得到了和实验结果较一致的结果。其 后至今，未发现此原理在艇决两相流流场中的进一步应用。本文的目的就是把能量耗散 极值原理引入水力机械中的含沙水流流场的研究中，利用雷诺输运定理，从能量方程着 《 手，推导出水、沙两相的能量耗散项，得出一种新的数值模型，为数值模拟含沙水流中 的水力机械过流部件区域的流场探索一条新的方法。 { 3．同液两相流的能量方程及数值模型 首先，为了能较方便地研究水沙两相流而又符合实际情况，作如下假设： I一23 1)固、液两相间无质量交换； 2)将固相看成“拟流体”，即沙粒群是充满整个研究空间的无间隙的“连续介质”； 3)所有颗粒均为尺寸相同的球粒； 4)两相温度相同，且通过水力机械时温度不变： 对于水力机械中的含沙水流来说，大量沙粒粒径约为O．03mm左右，极限容积中 包含，矿个沙粒的立方体边长为lcm，远小于实际水力机械的特征尺度，因此将沙粒群 近似为拟流体是合理的。 基于上述假设，应用欧拉法和拉格朗目法相结合的雷诺输运定理，对固、液两相 分别推导，得出各自的连续睦方程和动量方程如下： 液 续 方程 警 出～掣 固相连续性方程： 劐+生趋：。