李家峡水电站两相非恒定流数值计算.pdfVIP

李家峡水电站两相非恒定流数值计算 郭仁东那秀艳 (沈阳大学建筑与环境工程学院110044) 摘要本文根据李家峡水电站备用机组的尾水两相非恒定流问题。利用持征线差分方程进行数值计算．对 上，下游水击压力及水轮机转速的稳定性进行了分析．给出了交界面波动的影响范围． 关键词《l鼯水电站-《朔非恒定流．水击压力 。’、 ＼ 1 引言 李家峡水电站的备用机组尾水隧漏通入施工导流洞出流，可以节箭＆多工程量．导流洞出口水 深为3m．属于明渠流，长约400 m，丽尾水隧洞中大部分为有压流．尾水隧洞与水轮机之间设一调 压井(臣1)。这样．当水轮机导叶启闭时，便会发生明满两相非恒定流现象。本文考虑了尾水波动对 发电水头的影响，对机组卸荷与增荷两种情况进行了计算。结果分析表明，增荷情况下两相流波动 比较咀显，而卸荷情况上、下游水击压力较为不利． 旦I兰．——————————一 图1车家峡水电站备用机组纵剖面图 2数学模型的描述 2．1基本方程 控制有压非恒定流的基本方程是 连续方程等+V夏aH十i一面aV=o (1) 运动方程g蓑+宝+v塞=g(s。一s。) (2) 式中H为水头；V为流速；t为时间；x为沿管轴线方向坐标；a为水击波波速；g为重力加速度；s。为 底坡；st为能坡。 与之对应的一维明渠流基本微分方程是： 连续方程Od【k．V夏3Y十．丁A磊3V—u_ (3) 运动方程g纂+宝+v塞=g(so--sO (4) 式中A为过水断面；Y为水深；T为水面宽。对于上游压力引水管的非恒定流计算，引用式(1)、(2) 即可。尾水隧洞则要同时引用(1)、(2)、(3)、(4)进行明满两相非恒定流计算。 2．2有压非恒定流的差分方程 根据国内外较为成熟的经验(33、o]、(j]对压力水管、尾水隧洞及导流洞明满非恒定流计算， 有压部分内部采用显式的中心差分格式。边界处利用特征线差分方程求解。方程(1)、(2)的显式差 分格式表达如下： H，=HM--i1赢AtL-V“(HR—H。)+詈(VR--VL)] (5) (6) V。一Vw一寺意[g(HR--HL)+V,a(VR--VO?+g(So--Sm)△t 式中脚标P是At时段后中点计算值，M是△t时段前中点已知值，L、R为At时段前左右相邻点已 知值。 2．3特征线差分方程 由于尾水隧洞中同时存在有压流与明渠流两种流态，采用简洁的特征线方法较为合理的将两 种流态的计算方法加以统一。以明渠流的特征线差分方程为例表示如下 C— VP=C．+C．Y， (7) C+ VP=c。+c．Y。 (8) 式中 C。=g／cw，Cu是重波波速。 外．由于明渠流Cu值是随时间变化的，所以明渠特征差分方程中的时间差分项应采用Friedrichs [5]格式。 对于参数Lit与At’，cw与¨之间有～定差别，由于“比C。的数值大很多，则4∥不得不比At 减少许多来满足库朗稳定条件。图2表示了有压流和明渠流交界处的差分网格划分系统。本文根据 库朗条件dc一0．3s，4∥=0．02s，则有压流计算次数为明渠流25倍之比。 3有压相与明渠相的界定 当水轮机处流入尾水隧洞的流量增加时，尾水隧洞与导流洞衔接处的流态将从明渠相变为有 压相．流量达到稳定后．再降为明渠相。这样，就存在一个相