二维非恒定流流场与温度场的数值模拟.pdfVIP

第四届全国海事技术研讨会吏集 二维非恒定流流场和温度场的数值模拟 陆敬荣 (中国水利水电科学研究院，北京100044) 摘要本文采用二维非恒定流数学模型模拟预报火力发电厂的流场和温度场，其计算步骤 主要是根据连续方程、动量方程和能量方程，采用显示计算格式和正方形均匀网格法编制相应 的计算机程序，分别计算对应水位、流速和温度，并通过具体实例进行阐述和论证。成果表明 计算值与现场实测值十分吻合．这种方法，曾用在许多拟建火力发电厂工程中模拟预报各种不 同温排水域(海域、湖泊、水库，河流)的流场和温度场，为设计单位提供了可靠依据。 关键词 均匀闷格，流场，温度场．显示格式 0前言 二维非恒定流数学模型计算方法，曾用来解决许多工程实际问题．但过去在计算中蔫 要输入大量的相关数据，耗费巨多人力，而且一旦出现差错，极难奄找，工作异常繁重。为 此，作者根据多年来工程实践经验，反复思考，总结研制出一套适用性强、使用方便的温排 水通用数值模拟计算程序，只须输入边界值和少量有关数据，就可准确预报拟建火电厂各 种类型浅水域湓排水的流场和温度场，历经考验，屡试不爽。本文现以“马来西亚吉打洲 YAN镶火电厂温排水数值预报”课题中一种工况加以论证。 ：搬N镇电厂拟建于沿海Pulau 机组，冷却水量115．2ma／S，凝汽器温升7．13℃，以海水作为冷却水的直流供水系统水源， 采用上述计算程序预报电厂温排水随潮流运动的规律及取水温升，满足了设计需要。 J 数学模型 连续方程： 鲁+掣+掣=。 ㈤ 动量方程： X方向 争芸+V考=歉》恕勃一g象嘻华㈤ Y方向 cs， 能量方程： 百枷瓦椰石。否卜刮+万r刮一 署+社罢+V雾=丢(E为+号(E訇一 式中： 曰一水深，日=嚏+刁，魄一基准面以下水深， 卵一相对基准面潮位(参见右图)； 甜，V一分别为x，y方向流速； t一时间变量； c一谢才系数(采用曼宁公式计算)，c：三Ⅳ：，n一糙率系数； 玎 z乙，D，，￡k，zk一分别为各方向的广义粘性系数： p一水面超温，护=r一气，，一水面温度，7卜自然水温：E—熟扩散系数 △≯；肺口，肛-水面综合散热系数； p一水密度；g一重力加速度；Cp一水的比热。 2计算方法 深为5～15m，电厂附近平均水深约8m。由于水浅，涨落潮过程中同一垂线表底流向基 本一致，故选用二维非恒定流数学模型。 将以上计算水域划分为7200个正方形网格单元。取边长为250w，单元间水边界称为 通道。方程(1)～(4)为非线性偏微分方程组，目前尚未有解析解，因此将上述方程组转化 离散成代数方程组，然后求解。实际计算时，在单元形心处计算水位、水温。在通道上计 算流速、流量。取x，y同向流速为正(参见图1)。 L．． ．j‘ 计算《等吒．亡l喂屹 口一水位，木匠计算点 -卅l 卜-一直建计算点 什再尝‘·弓‘D嚣虬 图1计算格式示意图 ?15· 根据方程(1)，采用控制容积法计算出 各单元的水位。根据方程(2)～(3)，采用