CN114662425A - 一种水轮机启停工况流场仿真预测方法及系统 - Google Patents
一种水轮机启停工况流场仿真预测方法及系统 Download PDFInfo
- Publication number
- CN114662425A CN114662425A CN202210572725.2A CN202210572725A CN114662425A CN 114662425 A CN114662425 A CN 114662425A CN 202210572725 A CN202210572725 A CN 202210572725A CN 114662425 A CN114662425 A CN 114662425A
- Authority
- CN
- China
- Prior art keywords
- grid
- water turbine
- equation
- flow field
- starting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 103
- 238000000034 method Methods 0.000 title claims abstract description 74
- 238000004088 simulation Methods 0.000 title claims abstract description 51
- 239000012530 fluid Substances 0.000 claims abstract description 86
- 230000033001 locomotion Effects 0.000 claims abstract description 67
- 230000008859 change Effects 0.000 claims abstract description 56
- 238000004364 calculation method Methods 0.000 claims abstract description 54
- 230000010349 pulsation Effects 0.000 claims abstract description 25
- 238000012423 maintenance Methods 0.000 claims abstract description 22
- 230000008569 process Effects 0.000 claims abstract description 18
- 239000006185 dispersion Substances 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 4
- 238000010276 construction Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 238000013277 forecasting method Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 230000002123 temporal effect Effects 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 6
- 238000010168 coupling process Methods 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008034 disappearance Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/12—Simultaneous equations, e.g. systems of linear equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Control Of Water Turbines (AREA)
Abstract
本发明公开了一种水轮机启停工况流场仿真预测方法及系统,属于水力机械技术领域。现有技术无法科学、准确地预测水轮机在启停工况下的流场状态,进而将无法保证水轮机的安全运行。本发明构建动网格模型,从而可以根据水轮机启停工况下活动导叶运动过程,确定网格运动方式,在计算过程中,流体域网格几何形状可以自动变化,以匹配真实的物理状况。同时利用计算流体力学CFD模型,得到水轮机启停工况下状态变化,实现水轮机启停工况下内流场状态和压力脉动情况的仿真以及预测,进而可指导水电站运维。进而本发明能够准确实现启停工况下的数值模拟,有效提升水轮机的数值模拟精度,方案科学、合理,切实可行,利于推广,便于实施。
Description
技术领域
本发明涉及一种水轮机启停工况流场仿真预测方法及系统,属于水力机械技术领域。
背景技术
水电站凭借着启停迅速的优点,在电力系统中还发挥着系统调峰、调频、调相及事故备用的作用。
中国专利(公布号CN103853884A)公开了一种水轮机活动导叶振动特性预测方法,其利用设计工况流量作为流体计算的进口边界条件,求得导叶初始设计构型下的稳态流场作为流场初值条件;在一个时间步内交替调用结构计算模块和流体计算模块,满足收敛条件后,流体和结构计算整体同步向前推进;通过界面信息交换模块传递流固边界信息;输出时间历程上的结构振动位移。上述发明实现导叶与流场的耦合计算,本发明准确性较以往单纯的结构动力学方法和单向耦合方法有了显著的提高,而且能够观察整个导叶振动发展过程,有利于更好地指导结构设计,使流固耦合从理论研究走上实际工程应用。
然而水轮机启停的过程中,水轮机会经历从空载到满载的过程。由于处在非设计工况,尾水管进口处的水流具有较大的切向速度,使水流产生涡带,这会导致尾水管产生强烈的压力脉动现象。
这种突变的运行工况,由于负载的变化,会导致振动频发,这对整个机组的稳定性乃至电站安全都会有较大的影响。同时,这种振动对设备的冲击较大,特别是会对过流部件造成疲劳损伤,缩短设备运行寿命。因此,科学、准确地预测水轮机在启停工况下的流场状态对整个水电站的运维具有较大的指导意义。
但上述发明方案没有公开如何科学、合理地预测启停工况下的流场状态以及压力脉动信息,从而影响水轮机的安全运行。
同时,由于水轮机内部流场复杂,且无法安装监测系统来获得流体流动情况,想要获得准确的速度和压力脉动信息存在较大难度。为了配合电网调度,电站水轮机会频繁地经历启停工况,如果不能科学地预测启停工况下的压力脉动信息,将无法保证水轮机的安全运行。
进一步,随着计算流体力学(CFD)的不断发展,数值模拟的计算速度和结果准确性的大大提升,研究人员提出了通过数值仿真的方式来模拟水轮机内流场流动情况,获得尾水管压力脉动信息。
目前,针对各类水轮机各种定工况下的数值模拟研究已获得了不错的成果,但针对启停这一特殊工况下的研究目前尚未得到充分的研究。
网格作为CFD计算的输入文件,不仅包含了模拟区域的几何信息,网格质量也直接影响计算的速度和精度。水轮机通过改变导叶开度来调节工况,一次启动,活动导叶都经历从全闭到全开的过程,这就会造成水轮机流体域几何结构的变化,而且这个变化是基本覆盖水轮机整个启动过程。基于以上原因,加上水轮机内流场本身复杂程度较高,这就导致了启停工况数值模拟工作难以展开。
发明内容
针对现有技术的缺陷,本发明的目的一在于提供一种构建网格模型以及动网格模型,从而可以根据水轮机启停工况下活动导叶运动过程,确定网格运动方式,在计算过程中,流体域网格几何形状可以自动变化,以匹配真实的物理状况,实现流体动、静状态的准确描述;同时利用计算流体力学CFD模型,得到水轮机启停工况下状态变化,实现水轮机启停工况下内流场状态和压力脉动情况的仿真以及预测,进而可指导水电站运维的水轮机启停工况流场仿真预测方法。
本发明的目的二在于提供一种设置几何模块、网格模块、动网格模块对启停工况下,导叶开度变化导致的流体域几何模型变化进行准确描述,进而能够实现启停工况下的数值模拟仿真,有效提升水轮机的数值模拟精度;并利用流体力学计算模块、水电站运维模块,得到水轮机启停工况下状态变化,实现水轮机启停工况下内流场状态和压力脉动情况的仿真以及预测,进而可指导水电站运维,方案科学、合理,切实可行,利于推广,便于实施的水轮机启停工况流场预测系统。
为实现上述目的之一,本发明的第一技术方案为:
一种水轮机启停工况流场仿真预测方法,
包括如下步骤:
步骤1,根据水轮机流体域的特征,构建几何模型;
步骤2,对步骤1中的几何模型进行分块,生成网格模型,用以描述流体域静止状态下的几何形状;
步骤3,根据步骤2中的网格模型,确定边界面及边界条件;
步骤4,根据水轮机启停工况过程中,导叶开度变化情况,以及步骤3中的导叶边界面,确定流体域网格移动速度;
步骤5,根据步骤4中的流体域网格移动速度,构建动网格模型,用以生成变化的流体域网格,描述变化的几何形状,以匹配水轮机在实际启停工况下的物理状况;
步骤6,将步骤5中的动网格模型,装载到计算流体力学CFD模型中,求解得到水轮机启停工况下状态变化,实现水轮机启停工况流场的仿真;
所述状态变化包括内流场的状态变化、压力脉动状态信息;
步骤7,根据步骤6中的内流场的状态变化、压力脉动状态信息,预测水轮机启停工况下内流场状态和压力脉动情况,辅助水电站运维。
本发明经过不断探索以及试验,构建网格模型以及动网格模型,从而可以根据水轮机启停工况下活动导叶运动过程,确定网格运动方式,在计算过程中,流体域网格几何形状可以自动变化,以匹配真实的物理状况。同时利用计算流体力学CFD模型,得到水轮机启停工况下状态变化,实现水轮机启停工况下内流场状态和压力脉动情况的仿真以及预测,进而可指导水电站运维。
本发明充分考虑预测水轮机启停工况流场状态的需求,以及数值模拟存在的几个难题,设置网格模型,用于描述静止状态的流体域网格几何形状;并设置动网格模型对启停工况下,导叶开度变化导致的流体域几何模型变化进行准确描述,进而根据流体动、静状态的准确描述,能够实现启停工况下的数值模拟仿真,有效提升水轮机的数值模拟精度,能够有效确保水轮机的安全运行,方案科学、合理,切实可行,利于推广,便于实施。
作为优选技术措施:
所述步骤1中,几何模型利用水轮机流体域全流道的CAD模型进行构建。
作为优选技术措施:
所述步骤2中,所述网格模型根据湍流模型进行构建,并设置壁面网格厚度。
作为优选技术措施:
所述步骤4中,流体域网格移动速度的确定方法如下:
步骤41,获取每片导叶旋转轴的旋转轴坐标,确定每片导叶旋转中心;
步骤42,遍历边界面上的网格单元的坐标信息,获取每个边界面上的中心坐标,并使用绕轴旋转公式,确定叶片绕轴旋转的方向;
步骤43,根据步骤41中的旋转中心以及步骤42中的中心坐标,得到边界面相对于旋转中心的坐标;
步骤44,根据步骤43中的边界面相对于旋转中心的坐标,得到边界面距旋转中心的距离和相对于坐标轴的夹角;
步骤45,根据步骤44中的边界面距旋转中心的距离和相对于坐标轴的夹角,以及启停工况下的导叶转动速率变化,计算导叶随时间变化的转速;
步骤46,根据步骤45中的转速,计算每个边界面的网格运动速度。
作为优选技术措施:
边界面相对于旋转中心坐标的计算公式如下:
x0, y0为导叶旋转中心的坐标;
xn, yn为边界面的中心坐标;
边界面距旋转中心距离的计算公式如下:
边界面相对于坐标轴夹角的计算公式如下:
边界面网格运动速度的计算公式如下:
其中, w为导叶转动速度;
作为优选技术措施:
所述步骤5中,动网格模型通过使用任意拉格朗日-欧拉方法,求解网格运动的控制方程组,得到全局网格节点的运动信息,其生成方法如下:
步骤51,获取网格运动的速度场,其计算公式如下:
其中, v ( x , t)表示网格运动的速度场;
x ( x 0 , t)表示初始位于 x 0 的点在t时刻的坐标;
步骤52,根据步骤51中的速度场,通过莱布尼茨积分定则,计算上一时刻网格的某单元内张量场 v 在Ωc(t)中随时间的变化,其计算公式如下:
其中,Ω(t)为随时间变化的流体域,dS为面微元;
Ωc(t)为流体域Ω(t)的子域,其对应于上一时刻网格的一个单元;
d/dt表示积分量的总导数,积分随Ωc(t)变化而变化;
步骤53,根据步骤52中的张量场 v ,并利用莱布尼茨积分定则,得到网格单元的体积运动方程,其计算公式如下:
步骤54,对步骤53中的体积运动方程、网格质量方程、动量守恒方程作一个时间步上的积分,得到网格运动时间相关的控制方程组:
其中,所有的n上标表示第n时间步时的物理量,、、分别表示第n时间步
时,该网格单元的体积、质量和动量;tn为第n时间步;为Ωc的所有面单元,f表示中的
一个面,、、分别表示面f的面积、速度和应力张量。
步骤55,根据控制方程组,得到网格移动的全场解,根据全场解改变流体域网格几何形状,以匹配水轮机在实际启停工况下的物理状况。
在确定网格运动时,只需要确定边界面处的网格移动,通过使用任意拉格朗日-欧拉方法,求解网格运动守恒方程,得到全局网格节点的运动信息,方案简单实用,切实可行。
作为优选技术措施:
所述网格运动的速度场通过求解流体域Ω(t)上的泊松方程得到,所述泊松方程的计算公式如下:
其中:∂Ωimp(t)表示具有网格移动速度或位移的边界;
∂Ω\∂Ωimp(t)是固定边界;
Imposed velocity为定义的网格运动速度;
λ为单位矩阵。
水轮机启停工况下,由活动导叶开度变化造成流体域几何结构变化,在此动网格模型的基础上,只需得到活动导叶边界面上网格移动的速度,就能得到网格移动的全场解,实现网格变化的定量计算,方案切实可行。
作为优选技术措施:
所述步骤6中,计算流体力学CFD模型根据流体力学控制方程进行构建,其具体的构建方法如下:
步骤61,构建质量守恒方程和动量守恒方程,其计算公式分别如下所示:
步骤62,对步骤61中的质量守恒方程和动量守恒方程分别进行时间和空间上的离散,其具体包括以下内容:
所述质量守恒方程进行时间上离散的方程如下:
其中Δt为时间步的步长;
所述动量守恒方程进行时间上离散的方程如下:
对质量守恒方程和动量守恒方程,利用有限体积方法,进行空间上的离散;
步骤63,步骤62中的离散完成后,使用显式压力场进行速度预测,再进行一个校正计算,并利用连续性方程用于计算压力的变化;
所述校正的方法为通过调整质量通量以确保质量守恒,然后更新速度场。
作为优选技术措施:
所述单元体积的计算公式如下:
所述单元质量的计算公式如下:
所述单元动量的计算公式如下:
对Δt,第n步到第n+1步,一个时间步的积分的计算公式如下:
(.)n表示第n时间步的物理量。
为实现上述目的之一,本发明的第二技术方案为:
一种水轮机启停工况流场预测系统,
应用上述的一种水轮机启停工况流场仿真预测方法;
其包括几何模块、网格模块、动网格模块、流体力学计算模块、水电站运维模块;
所述几何模块,用于描述水轮机流体域的几何特征;
所述网格模块,用于描述初始状态的流体域网格几何形状;
所述动网格模块,用于描述变化状态的流体域网格几何形状;
所述流体力学计算模块,用于计算水轮机启停工况下内流场的状态变化、压力脉动状态信息;
所述水电站运维模块,用于指导水电站运维。
本发明充分考虑预测水轮机启停工况流场状态的需求,以及数值模拟存在的几个难题,设置几何模块、网格模块、动网格模块对启停工况下,导叶开度变化导致的流体域几何模型变化进行准确描述,进而能够实现启停工况下的数值模拟仿真,有效提升水轮机的数值模拟精度;并利用流体力学计算模块、水电站运维模块,得到水轮机启停工况下状态变化,实现水轮机启停工况下内流场状态和压力脉动情况的仿真以及预测,进而可指导水电站运维,方案科学、合理,切实可行,利于推广,便于实施。
与现有技术相比,本发明具有以下有益效果:
本发明经过不断探索以及试验,构建网格模型以及动网格模型,从而可以根据水轮机启停工况下活动导叶运动过程,确定网格运动方式,在计算过程中,流体域网格几何形状可以自动变化,以匹配真实的物理状况。同时利用计算流体力学CFD模型,得到水轮机启停工况下状态变化,实现水轮机启停工况下内流场状态和压力脉动情况的仿真以及预测,进而可指导水电站运维。
进一步,本发明充分考虑预测水轮机启停工况流场状态的需求,以及数值模拟存在的几个难题,设置几何模块、网格模块、动网格模块对启停工况下,导叶开度变化导致的流体域几何模型变化进行准确描述,进而根据流体动、静状态的准确描述,能够实现启停工况下的数值模拟仿真,有效提升水轮机的数值模拟精度;并利用流体力学计算模块、水电站运维模块,得到水轮机启停工况下状态变化,实现水轮机启停工况下内流场状态和压力脉动情况的仿真以及预测,进而可指导水电站运维,方案科学、合理,切实可行,利于推广,便于实施。
附图说明
图1为本发明水轮机启停工况下预测方法流程图;
图2为本发明导叶旋转导致开度变化状态图;
图3为本发明仿真实施例水轮机结构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
相反,本发明涵盖任何由权利要求定义的在本发明的精髓和范围上做的替代、修改、等效方法以及方案。进一步,为了使公众对本发明有更好的了解,在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。对本领域技术人员来说没有这些细节部分的描述也可以完全理解本发明。
如图1所示,本发明水轮机启停工况流场仿真预测方法的一种具体实施例:
一种水轮机启停工况流场仿真预测方法,
包括如下步骤:
步骤1,根据水轮机流体域的特征,构建几何模型;
步骤2,对步骤1中的几何模型进行分块,生成网格模型,用以描述流体域静止状态下的几何形状;
步骤3,根据步骤2中的网格模型,确定边界面及边界条件;
步骤4,根据水轮机启停工况过程中,导叶开度变化情况,可参见图2,以及步骤3中的导叶边界面,确定流体域网格移动速度;
步骤5,根据步骤4中的流体域网格移动速度,构建动网格模型,用以生成变化的流体域网格,描述变化的几何形状,以匹配水轮机在实际启停工况下的物理状况;
步骤6,将步骤5中的动网格模型,装载到计算流体力学CFD模型中,求解得到水轮机启停工况下状态变化,实现水轮机启停工况流场的仿真;
所述状态变化包括内流场的状态变化、压力脉动状态信息;
步骤7,根据步骤6中的内流场的状态变化、压力脉动状态信息,预测水轮机启停工况下内流场状态和压力脉动情况,辅助水电站运维。
本发明水轮机启停工况流场仿真预测方法的一种最佳具体实施例:
一种水轮机启停工况流场仿真预测方法,包括以下内容:
首先,给出在随时间变化的流体域Ω(t)内的守恒控制方程。用v(x, t)表示网格运动的速度场。其中,点x(x 0, t)表示初始位于x 0的点在t时刻的坐标。因此网格运动的速度场的计算公式如下:
现在考虑Ω(t)的子域Ωc(t),该域对应于上一时刻网格的一个单元。
使用莱布尼茨积分定则(或Reynolds输运定理),可以得到任意张量场 v 在Ωc(t)中随时间的变化(这里v为一阶张量,即向量,莱布尼茨积分定则适用于任意阶数的张量),其具体的计算公式如下:
这里d/dt表示积分量的总导数,积分随Ωc(t)变化而变化。
对单位标量场使用莱布尼茨积分定则,得到网格单元体积运动的形式。
对体积、质量、动量三个方程作一个时间步上的积分,得到网格运动时间相关的守恒方程:
这样就得到了网格运动的控制方程组,任意的网格运动速度场v(x, t)均需满足这三个控制方程。
针对具体的网格移动速度场则需要通过求解Ω(t)上的泊松方程得到:
其中∂Ωimp(t)表示具有网格移动速度或位移的边界,∂Ω\∂Ωimp(t)是固定边界。也可以在域内定义网格位移或者固定网格。
张量λ通常被定义为单位矩阵,也可以定义各向异性的张量λ以使网格在特定方向和给定区域中具有更强的刚性。但是,具有各向异性的λ会导致网格速度v的分量之间存在强耦合。
水轮机启停工况下,由活动导叶开度变化造成流体域几何结构变化,在此动网格算法的基础上,只需确定活动导叶边界面上网格移动的速度,就能得到网格移动的全场解,实现网格变化。
一般情况下,一台水轮机具有20片活动导叶,以圆周均匀间隔的方式分布在导水区。每片导叶都有固定的旋转轴,首先定位旋转轴的坐标 (x0, y0, 0) ,再获取边界面上的网格位置信息(xn, yn, zn)。一般情况下,活动导叶的旋转为竖直方向即z轴,故保持z不变考虑平面转动即可。进行坐标变化,将坐标原点变为旋转中心,其计算公式如下:
距旋转中心的距离r和与x轴的夹角θ为:
在转动角速度为ω的情况下,点(xn, yn, zn)转动速度为:
在计算流体力学CFD模型求解中,通过使用UDF(用户自定义函数)的方法,确定网格运动的边界条件,确定每个网格节点的运动速度,其具体包括以下内容:
首先使用获取网格坐标的函数,获得边界面上网格位置信息;
再对边界面上的网格节点做循环,使用上述方法依次计算距旋转中心的距离r和与x轴的夹角θ,再根据转动速度公式,计算节点的运动速度,最终确定边界面网格的运动速度。
本发明充分考虑现有技术中针对预测水轮机启停工况流场状态的需求,以及数值模拟存在的难题:如何准确描述以及计算启停工况下,导叶开度变化导致的流体域几何模型变化。
进而,本发明构建若干模型,可以根据水轮机启停工况下活动导叶运动过程,确定网格运动方式,在计算过程中,流体域网格几何形状自动变化,以匹配真实的物理状况。在确定网格运动时,只需要确定边界面处的网格移动,通过使用任意拉格朗日-欧拉方法,求解网格运动守恒方程,得到全局网格节点的运动信息。
本发明计算流体力学CFD模型的一种具体实施例:
计算流体力学CFD模型根据流体力学控制方程进行构建。
流体力学控制方程包括质量守恒和动量守恒方程,其具体的计算公式如下:
接下来对流体力学控制方程进行时间和空间上的离散。
该计算流体力学CFD模型使用有限体积方法进行计算,通过对网格单元进行积分,确定外延量,其计算公式如下:
对Δt(第n步到第n+1步,一个时间步)的积分计算公式如下:
其中,(.)n表示第n时间步的物理量。
对质量守恒方程进行时间上的离散,其计算公式如下:
对动量守恒方程进行时间上的离散,其计算公式如下:
对空间上的离散,选择CFD常用的有限体积方法离散方程,首先使用显式压力场进行速度预测,再进行一个校正计算,其中连续性方程用于计算压力的变化。
校正步骤调整质量通量以确保质量守恒,然后更新速度场。
本发明网格运动求解方法的一种具体实施例:
网格运动求解方法包括以下内容:
首先,在计算流体力学CFD模型中添加动网格模型,该动网格模型添加在速度场求解之后。
然后求解Ωn域上的预测速度并得到修正后的速度场u c n+1,k后(n为当前时间步,n+1为下一时间步,k为当前迭代步)。
再根据用户输入的网格位移边界条件求解网格运动的三个控制方程和泊松方程,得到网格运动的全场解v c|n n+1,k。
最后,根据网格运动的全场解,重构出每个网格节点的速度v f|n n+1,k,由此得到每个节点的位移,构建出网格运动后的求解域Ωn+1。
本发明应用在如图3所示的Francis 99混流式水轮机上的一种具体实施例:
一种应用在Francis 99混流式水轮机的启停工况流场预测方法,包括以下内容:
步骤1:构建水轮机流体域全流道的CAD模型。由于水轮机结构复杂,需要根据详细的工程图,构建出准确的CAD模型。在CAD模型的基础上,使用网格生成技术,对几何体分块,生成结构化的网格模型,并根据将要选用的湍流模型对第一层网格Y+要求,设置壁面网格厚度。
步骤2:将网格模型导入计算流体力学CFD模型,根据启停工况参数确定边界条件。该步骤需要获取启停工况下,入口流量和出口压力随时间变化的值。同时,使用旋转机械模块,用以模拟转子-定子间的相互作用。在该模块下,确定转轮区域的旋转轴以及启停工况下,随时间变化的转轮转速。
步骤3:根据导叶在启停工况下的运动规律,确定导叶边界面上的网格移动速度,并在每片导叶的边界面上确定一个边界条件,其具体包括内容:
首先,确定每片导叶旋转轴的坐标,在遍历边界面上的网格单元的坐标信息,使用绕轴旋转公式,确定叶片绕轴旋转的方向。并根据启停工况下的导叶转动速率变化确定随时间变化的转速。该部分可以使用UDF(用户自定义函数)完成。
根据本发明描述的叶片旋转计算公式,首先确定一片叶片的旋转中心(x0, y0,0),再遍历该叶片上的所有边界面单元,获取每个边界面单元的中心坐标(xn, yn, zn)。
假设导叶转动速度为ω,则确定每个边界面单元的网格运动速度为:
以此类推,根据每个叶片的旋转中心,确定所有叶片边界面网格的旋转速度。
步骤4:最后完成数值计算相关参数设置,湍流模型选择K-ω SST、时间步长设置为Δt=1*10-5s、离散方式选择二阶中心差分,求解精度设置为1*10-8。确定完所有数值参数后,使用计算流体力学CFD模型运算。
在计算结束后,获得整个启停工况的全场解。
本发明运维过程的一种具体实施例:
首先,确定尾水管测点的位置,使用后处理软件,获得测点位置的压力变化信息。
再以速度云图形式显示结果,使用后处理软件的播放功能,观察在整个启停工况下,内流场的状态变化,以及尾水管涡带的形成、增强和消失。
最后结合压力脉动信息,分析获取压力脉动最为剧烈的工况区间,指导运维人员尽量避免该工况。
应用本发明方法的一种装置实施例:
一种计算机设备,其包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序;
当所述一个或多个程序被所述一个或多个处理器执行时,使得所述一个或多个处理器实现上述的一种水轮机启停工况流场仿真预测方法。
应用本发明方法的一种计算机介质实施例:
一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现上述的一种水轮机启停工况流场仿真预测方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。
Claims (10)
1.一种水轮机启停工况流场仿真预测方法,其特征在于,
包括如下步骤:
步骤1,根据水轮机流体域的特征,构建几何模型;
步骤2,对步骤1中的几何模型进行分块,生成网格模型,用以描述流体域静止状态下的几何形状;
步骤3,根据步骤2中的网格模型,确定边界面及边界条件;
步骤4,根据水轮机启停工况过程中,导叶开度变化情况,以及步骤3中的导叶边界面,确定流体域网格移动速度;
步骤5,根据步骤4中的流体域网格移动速度,构建动网格模型,用以生成变化的流体域网格,描述变化的几何形状,以匹配水轮机在实际启停工况下的物理状况;
步骤6,将步骤5中的动网格模型,装载到计算流体力学CFD模型中,求解得到水轮机启停工况下状态变化,实现水轮机启停工况流场的仿真;
所述状态变化包括内流场的状态变化、压力脉动状态信息;
步骤7,根据步骤6中的内流场的状态变化、压力脉动状态信息,预测水轮机启停工况下内流场状态和压力脉动情况,辅助水电站运维。
2.如权利要求1所述的一种水轮机启停工况流场仿真预测方法,其特征在于,
所述步骤1中,几何模型利用水轮机流体域全流道的CAD模型进行构建。
3.如权利要求1所述的一种水轮机启停工况流场仿真预测方法,其特征在于,
所述步骤2中,所述网格模型根据湍流模型进行构建,并设置壁面网格厚度。
4.如权利要求1所述的一种水轮机启停工况流场仿真预测方法,其特征在于,
所述步骤4中,流体域网格移动速度的确定方法如下:
步骤41,获取每片导叶旋转轴的旋转轴坐标,确定每片导叶旋转中心;
步骤42,遍历边界面上的网格单元的坐标信息,获取每个边界面上的中心坐标,并使用绕轴旋转公式,确定叶片绕轴旋转的方向;
步骤43,根据步骤41中的旋转中心以及步骤42中的中心坐标,得到边界面相对于旋转中心的坐标;
步骤44,根据步骤43中的边界面相对于旋转中心的坐标,得到边界面距旋转中心的距离和相对于坐标轴的夹角;
步骤45,根据步骤44中的边界面距旋转中心的距离和相对于坐标轴的夹角,以及启停工况下的导叶转动速率变化,计算导叶随时间变化的转速;
步骤46,根据步骤45中的转速,计算每个边界面的网格运动速度。
6.如权利要求1所述的一种水轮机启停工况流场仿真预测方法,其特征在于,
所述步骤5中,动网格模型通过使用任意拉格朗日-欧拉方法,求解网格运动的控制方程组,得到全局网格节点的运动信息,其生成方法如下:
步骤51,获取网格运动的速度场,其计算公式如下:
其中, v ( x , t)表示网格运动的速度场;
x ( x 0 , t)表示初始位于 x 0 的点在t时刻的坐标;
步骤52,根据步骤51中的速度场,通过莱布尼茨积分定则,计算上一时刻网格的某单元内张量场 v 在Ωc(t)中随时间的变化,其计算公式如下:
其中,Ω(t)为随时间变化的流体域;dS为面微元;
Ωc(t)为流体域Ω(t)的子域,其对应于上一时刻网格的一个单元;
d/dt表示积分量的总导数,积分随Ωc(t)变化而变化;
步骤53,根据步骤52中的张量场 v ,并利用莱布尼茨积分定则,得到网格单元的体积运动方程,其计算公式如下:
步骤54,对步骤53中的体积运动方程、网格质量方程、动量守恒方程作一个时间步上的积分,得到网格运动时间相关的控制方程组:
步骤55,根据控制方程组,得到网格移动的全场解,根据全场解改变流体域网格几何形状,以匹配水轮机在实际启停工况下的物理状况。
8.如权利要求7所述的一种水轮机启停工况流场仿真预测方法,其特征在于,
所述步骤6中,计算流体力学CFD模型根据流体力学控制方程进行构建,其具体的构建方法如下:
步骤61,构建质量守恒方程和动量守恒方程,其计算公式分别如下所示:
步骤62,对步骤61中的质量守恒方程和动量守恒方程分别进行时间和空间上的离散,其具体包括以下内容:
所述质量守恒方程进行时间上离散的方程如下:
其中Δt为时间步的步长;
所述动量守恒方程进行时间上离散的方程如下:
对质量守恒方程和动量守恒方程,利用有限体积方法,进行空间上的离散;
步骤63,步骤62中的离散完成后,使用显式压力场进行速度预测,再进行一个校正计算,并利用连续性方程用于计算压力的变化;
所述校正的方法为通过调整质量通量以确保质量守恒,然后更新速度场。
10.一种水轮机启停工况流场预测系统,其特征在于,
应用如权利要求1-9任一所述的一种水轮机启停工况流场仿真预测方法;
其包括几何模块、网格模块、动网格模块、流体力学计算模块、水电站运维模块;
所述几何模块,用于描述水轮机流体域的几何特征;
所述网格模块,用于描述初始状态的流体域网格几何形状;
所述动网格模块,用于描述变化状态的流体域网格几何形状;
所述流体力学计算模块,用于计算水轮机启停工况下内流场的状态变化、压力脉动状态信息;
所述水电站运维模块,用于指导水电站运维。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210572725.2A CN114662425B (zh) | 2022-05-25 | 2022-05-25 | 一种水轮机启停工况流场仿真预测方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210572725.2A CN114662425B (zh) | 2022-05-25 | 2022-05-25 | 一种水轮机启停工况流场仿真预测方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114662425A true CN114662425A (zh) | 2022-06-24 |
CN114662425B CN114662425B (zh) | 2022-09-20 |
Family
ID=82038197
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210572725.2A Active CN114662425B (zh) | 2022-05-25 | 2022-05-25 | 一种水轮机启停工况流场仿真预测方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114662425B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115906718A (zh) * | 2023-03-09 | 2023-04-04 | 陕西空天信息技术有限公司 | 一种旋转机械cfd系统 |
CN116090138A (zh) * | 2023-04-03 | 2023-05-09 | 浙江远算科技有限公司 | 一种基于数据监测的水轮机转轮疲劳仿真计算方法和系统 |
CN116992747A (zh) * | 2023-09-28 | 2023-11-03 | 深圳十沣科技有限公司 | 一种基于sph流固耦合的冲击式水轮机动力学分析方法 |
CN117436322A (zh) * | 2023-12-21 | 2024-01-23 | 浙江远算科技有限公司 | 基于叶素理论的风力机叶片气动弹性仿真方法和介质 |
CN117473906A (zh) * | 2023-12-26 | 2024-01-30 | 浙江远算科技有限公司 | 一种基于流体动力学仿真的风电机舱后处理方法和介质 |
CN118070456A (zh) * | 2024-04-19 | 2024-05-24 | 上海东方泵业集团南通有限公司 | 一种大型双吸泵意外停机瞬态评估方法及水锤防护装置 |
CN118313227A (zh) * | 2024-06-07 | 2024-07-09 | 中国地质科学院地质力学研究所 | 一种高位滑坡堰塞坝数值仿真方法及装置 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102799730A (zh) * | 2012-07-13 | 2012-11-28 | 北京航空航天大学 | 一种燃气轮机风扇叶片反扭过程的预估方法 |
US8577648B1 (en) * | 2011-04-11 | 2013-11-05 | The United States Of America As Represented By The Secretary Of The Navy | Simulating fluid flow at a moving boundary |
CN103853884A (zh) * | 2014-02-24 | 2014-06-11 | 昆明理工大学 | 一种水轮机活动导叶振动特性预测方法 |
WO2014102552A1 (ru) * | 2012-12-28 | 2014-07-03 | Abramyan Vitaly | Способ придания движения рабочему колесу и рабочее колесо гидротурбины |
CN107480324A (zh) * | 2017-06-30 | 2017-12-15 | 河海大学 | 一种风电机组叶片状态监测系统构建方法 |
CN109035387A (zh) * | 2018-06-26 | 2018-12-18 | 国家电网有限公司 | 一种基于水锤效应和动网格理论的抽水蓄能电站过渡过程三维模拟方法 |
CN109977345A (zh) * | 2019-01-29 | 2019-07-05 | 河海大学 | 一种轴流泵叶顶间隙泄漏涡空化的数值模拟方法 |
CN110727996A (zh) * | 2019-09-17 | 2020-01-24 | 北京理工大学 | 适用于动边界绕流的湍流模型修正方法 |
CN113723030A (zh) * | 2021-10-18 | 2021-11-30 | 山东大学 | 基于计算流体力学的实际气体物性仿真方法及系统 |
CN113947003A (zh) * | 2021-10-15 | 2022-01-18 | 西安交通大学 | 一种面向热流耦合场景的粒子型无网格仿真系统 |
-
2022
- 2022-05-25 CN CN202210572725.2A patent/CN114662425B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8577648B1 (en) * | 2011-04-11 | 2013-11-05 | The United States Of America As Represented By The Secretary Of The Navy | Simulating fluid flow at a moving boundary |
CN102799730A (zh) * | 2012-07-13 | 2012-11-28 | 北京航空航天大学 | 一种燃气轮机风扇叶片反扭过程的预估方法 |
WO2014102552A1 (ru) * | 2012-12-28 | 2014-07-03 | Abramyan Vitaly | Способ придания движения рабочему колесу и рабочее колесо гидротурбины |
CN103853884A (zh) * | 2014-02-24 | 2014-06-11 | 昆明理工大学 | 一种水轮机活动导叶振动特性预测方法 |
CN107480324A (zh) * | 2017-06-30 | 2017-12-15 | 河海大学 | 一种风电机组叶片状态监测系统构建方法 |
CN109035387A (zh) * | 2018-06-26 | 2018-12-18 | 国家电网有限公司 | 一种基于水锤效应和动网格理论的抽水蓄能电站过渡过程三维模拟方法 |
CN109977345A (zh) * | 2019-01-29 | 2019-07-05 | 河海大学 | 一种轴流泵叶顶间隙泄漏涡空化的数值模拟方法 |
CN110727996A (zh) * | 2019-09-17 | 2020-01-24 | 北京理工大学 | 适用于动边界绕流的湍流模型修正方法 |
CN113947003A (zh) * | 2021-10-15 | 2022-01-18 | 西安交通大学 | 一种面向热流耦合场景的粒子型无网格仿真系统 |
CN113723030A (zh) * | 2021-10-18 | 2021-11-30 | 山东大学 | 基于计算流体力学的实际气体物性仿真方法及系统 |
Non-Patent Citations (9)
Title |
---|
HASAN AKIN等: "A CFD aided hydraulic turbine design methodology applied to Francis turbines", 《IEEE》 * |
刘宜等: "水泵水轮机同步启动过程内部流场分析", 《甘肃科学学报》 * |
徐连奎等: "基于CFD的水轮机模型建立与数值仿真", 《软件导刊》 * |
王洪杰等: "基于动网格水泵水轮机泵工况的压力脉动分析", 《大电机技术》 * |
程志远等: "移动网格技术在计算流体动力学数值仿真中的应用", 《重庆大学学报》 * |
闫云雷: "贯流式水轮机多相流数值模拟与整机振动特性分析", 《中国优秀硕士学位论文全文数据库(电子期刊)》 * |
霍天浪等: "导叶关闭速度对尾水管流场影响的数值模拟", 《陕西水利》 * |
黄剑峰等: "基于动网格的活动导叶流道内湍流场数值模拟", 《排灌机械工程学报》 * |
黄剑峰等: "混流式水轮机全流道三维非定常湍流场的动态大涡模拟", 《中国农村水利水电》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115906718A (zh) * | 2023-03-09 | 2023-04-04 | 陕西空天信息技术有限公司 | 一种旋转机械cfd系统 |
CN116090138A (zh) * | 2023-04-03 | 2023-05-09 | 浙江远算科技有限公司 | 一种基于数据监测的水轮机转轮疲劳仿真计算方法和系统 |
CN116992747A (zh) * | 2023-09-28 | 2023-11-03 | 深圳十沣科技有限公司 | 一种基于sph流固耦合的冲击式水轮机动力学分析方法 |
CN116992747B (zh) * | 2023-09-28 | 2024-03-22 | 深圳十沣科技有限公司 | 一种基于sph流固耦合的冲击式水轮机动力学分析方法 |
CN117436322A (zh) * | 2023-12-21 | 2024-01-23 | 浙江远算科技有限公司 | 基于叶素理论的风力机叶片气动弹性仿真方法和介质 |
CN117436322B (zh) * | 2023-12-21 | 2024-04-19 | 浙江远算科技有限公司 | 基于叶素理论的风力机叶片气动弹性仿真方法和介质 |
CN117473906A (zh) * | 2023-12-26 | 2024-01-30 | 浙江远算科技有限公司 | 一种基于流体动力学仿真的风电机舱后处理方法和介质 |
CN117473906B (zh) * | 2023-12-26 | 2024-04-19 | 浙江远算科技有限公司 | 一种基于流体动力学仿真的风电机舱后处理方法和介质 |
CN118070456A (zh) * | 2024-04-19 | 2024-05-24 | 上海东方泵业集团南通有限公司 | 一种大型双吸泵意外停机瞬态评估方法及水锤防护装置 |
CN118313227A (zh) * | 2024-06-07 | 2024-07-09 | 中国地质科学院地质力学研究所 | 一种高位滑坡堰塞坝数值仿真方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN114662425B (zh) | 2022-09-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114662425B (zh) | 一种水轮机启停工况流场仿真预测方法及系统 | |
Wu et al. | Numerical prediction and similarity study of pressure fluctuation in a prototype Kaplan turbine and the model turbine | |
CN102913464B (zh) | 离心泵转子瞬态流固耦合特性预测方法 | |
CN109753716B (zh) | 基于流场仿真的核/火电汽轮机组流体激励数值计算方法及系统 | |
Wang et al. | Simulation of multistage compressor at off-design conditions | |
Ekici et al. | Harmonic balance analysis of blade row interactions in a transonic compressor | |
Kersken et al. | Time-Linearized and Time-Accurate 3D RANS Methods for Aeroelastic Analysis in Turbomachinery | |
CN115906718B (zh) | 一种旋转机械cfd系统 | |
Guo et al. | A CFD/CSD model for aeroelastic calculations of large-scale wind turbines | |
Vahdati et al. | Unsteady flow and aeroelasticity behavior of aeroengine core compressors during rotating stall and surge | |
Wang et al. | Towards massively parallel large eddy simulation of turbine stages | |
Anciger et al. | Prediction of rotating stall and cavitation inception in pump turbines | |
CN104166752A (zh) | 液力变矩器全流道瞬态数值模拟计算方法 | |
CN103631994A (zh) | 一种风力机噪声声辐射规律数值预测的方法 | |
Weiss et al. | Simulation of unsteady turbomachinery flows using an implicitly coupled nonlinear harmonic balance method | |
CN111985166A (zh) | 隐式考虑动态摩阻的管道水力瞬变模拟方法和存储介质 | |
Gourdain et al. | High-performance computing to simulate large-scale industrial flows in multistage compressors | |
CN105989205A (zh) | 飞行器表面脉动压力的确定方法 | |
Chwalowski et al. | Fun3d analyses in support of the second aeroelastic prediction workshop | |
Yang et al. | Gas kinetic flux solver based high-order finite-volume method for simulation of two-dimensional compressible flows | |
Brandsen | Prediction of axial compressor blade vibration by modelling fluid-structure interaction | |
CN106599395A (zh) | 一种油浸变压器的噪声数值仿真计算方法 | |
Camara | Validation of time domain flutter predictiontool with experimental results | |
Spiering | Coupling of TAU and TRACE for parallel accurate flow simulations | |
CN104753441B (zh) | 一种伺服电机的基于k‑观测器的滑模预测控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |