CN113435136A - 耦合能量方程的输水管道气-液两相均质流的模拟方法 - Google Patents
耦合能量方程的输水管道气-液两相均质流的模拟方法 Download PDFInfo
- Publication number
- CN113435136A CN113435136A CN202110705741.XA CN202110705741A CN113435136A CN 113435136 A CN113435136 A CN 113435136A CN 202110705741 A CN202110705741 A CN 202110705741A CN 113435136 A CN113435136 A CN 113435136A
- Authority
- CN
- China
- Prior art keywords
- gas
- equation
- wave
- liquid
- mixed fluid
- 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
- 238000000034 method Methods 0.000 title claims abstract description 93
- 239000007788 liquid Substances 0.000 title claims abstract description 87
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 52
- 238000004088 simulation Methods 0.000 title claims abstract description 27
- 230000008878 coupling Effects 0.000 title claims abstract description 11
- 238000010168 coupling process Methods 0.000 title claims abstract description 11
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims abstract description 28
- 230000005514 two-phase flow Effects 0.000 claims abstract description 21
- 239000012530 fluid Substances 0.000 claims description 73
- 230000035939 shock Effects 0.000 claims description 32
- 239000000203 mixture Substances 0.000 claims description 18
- 230000008569 process Effects 0.000 claims description 14
- 230000004907 flux Effects 0.000 claims description 9
- 238000011144 upstream manufacturing Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000008859 change Effects 0.000 claims description 5
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 3
- 238000003379 elimination reaction Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 2
- 238000004590 computer program Methods 0.000 description 13
- 238000010586 diagram Methods 0.000 description 9
- 230000001052 transient effect Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000003068 static effect Effects 0.000 description 4
- 238000012546 transfer Methods 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 241000143252 Idaea infirmaria Species 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000009877 rendering Methods 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- 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)
- Fluid Mechanics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computing Systems (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于有限体积法耦合能量方程的输水管道气‑液两相均质流的模拟方法,包括:通过引入能量方程构建含自由气体均质两相流结合能量方程的双曲偏微分方程***模型;构建出***模型下各基础波的存在条件以及判定模型;构建***模型下断点处三种近似黎曼求解器;通过重建流动变量,使得***模型在空间上达到二阶精度;根据已有的边界条件,通过牛顿迭代法利用黎曼不变量迭代得到二阶的边界值;引入源项,求得格式离散方程;建立稳定性约束条件,更新初始变量进行下一时步计算。本发明基于有限体积法Godunov格式,通过引入了气液两相均质流能量方程的控制模型,提供了一种算法简单,易于实现的模拟方法,以在高计算效率的前提下得到高精度解。
Description
技术领域
本发明属于水电站水力学数值模拟计算技术领域领域,具体涉及一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法。
背景技术
在输水管道***中,阀门或机组的突然启闭可能导致管道内的压力突变,当管道内的压力变化超过管壁材料的承受能力时,可能会造成严重的管道***破坏,甚至会威胁人身安全。近些年,管道内的两相流问题一直是国内外研究的热门课题,不同于一般的单相流问题,两相流需要考虑不同相间的耦合作用。在两相流问题中,均质流模型被认为是应用最广的模型。
目前,针对输水管道内的气-液两相均质流问题,其模拟的方法主要为特征线法(MOC,Method ofCharacteristics)以及有限差分法(FDM,Finite Difference Method)。MOC由于其计算简单,可以较好的模拟管道内压力波动而被广泛应用。然而,由于均质流模型变波速、变密度的问题,MOC需要做插值运算而造成模型预测与实际相差较大。FDM则避免了线性插值问题,在计算精度上相较于MOC提升了不少。随着均质流模型的愈发成熟,寻求更高精度的模型成为了学术界和工业界的目标。最近,Guinot等人(2000)将该方案应用到水锤模拟中,提出了两种有效的近似求解方法。随后,Guinot等人(2001)建立了一阶Godunov方法描述的两相流模型,该格式大大的提高了模型的计算效率,并进一步指出数值扩散带来了严重的模型退化问题。
所以,需要一个新的技术方案来解决上述问题。
发明内容
发明目的:为了克服现有技术中存在的在模拟输水管道中气-液两相均质流时存在计算精度低等不足,提供一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,基于有限体积法Godunov格式,通过引入了气液两相均质流能量方程的控制模型,提供了一种算法简单,易于实现的模拟方法,以在高计算效率的前提下得到高精度解。
技术方案:为实现上述目的,本发明提供一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,包括如下步骤:
S1:通过引入能量方程构建含自由气体均质两相流结合能量方程的3×3双曲偏微分方程(PDE)***模型,根据模拟工况确定计算域、初始条件以及边界条件;
S2:构建出3×3双曲偏微分方程***模型下各基础波(接触波、激波、稀疏波)的存在条件以及判定模型;
S3:构建双曲偏微分方程***模型下断点处三种近似黎曼求解器,分别为一种迭代和两种非迭代近似黎曼求解器;
S4:通过重建流动变量,使得***模型在空间上达到二阶精度;
S5:根据已有的边界条件,通过牛顿迭代法利用黎曼不变量迭代得到二阶的边界值;
S6:引入源项,通过基于二阶Longge kuta离散格式的时间***法,求得时间梯度上二阶有限体积法Godunov格式离散方程,得到二阶精度;
S7:建立稳定性约束条件,更新初始变量进行下一时步计算,获取到最大允许时间步长。
进一步地,所述步骤S1中双曲偏微分方程(PDE)***模型的构建是通过在双曲偏微分方程(PDE)***体系下,构建含自由气体的两相均质流基本控制方程。
进一步地,所述步骤S1中需在水锤问题的基础上假定:(a)管道内自由气体的体积占比很小(<≈1%)且分散均匀,气-水两相流以等效单相流体处理;(b)瞬变过程时间尺度很小,整个过程中水气间的质量传递过程被忽略;(c)管道内流体为无粘性流动且整个过程发生在等温条件下不考虑热传递。在上述条件的基础上,含自由气体的两相均质流基本控制方程具体如下:
包含能量方程的水锤基本方程:
其中,矩阵代表气液混合流体的特征变量;矩阵代表通量;沿管线距离x与时间t是自变量;ρm(x,t)是平均截面气液混合流体的密度;V(x,t)是气液混合体平均截面速度;P(x,t)是气液混合流体的绝对压力;E(x,t)是气液混合流体单位体积总能量;
气-液两相流水锤波速与压力的关系方程:
气-液两相流压力与密度的关系方程:
进一步地,所述步骤S2中基础波包括接触波、激波、稀疏波,构建方法具体如下:
接触波:
对于两相流***中的接触波关系,可以利用欧拉方程中的特征结构分析断点处的初值,利用一般黎曼不变量,穿过波λ=V,得到,
这里Pr代表参考压力;u1、u2、u3见公式(2);ρmr为参考压力下气液混合流体的密度;cm为气液混合流体的波速;er为参考压力下气液混合体的内能;进一步简化后,黎曼不变量可以写做:
这里星号(*)代表中间状态,P、V、ρm分别代表混合流体的压力、速度、密度,L和R分别代表黎曼问题中接触波左右区间,在等熵条件下,穿过接触波,流体特征不发生变化;
稀疏波:
相较于激波,稀疏波表现的更为平滑,穿过波λ=V-cm,利用一般黎曼不变量,得到
进一步简化后,该微分方程可以写做
dy+cmd(lnpm)=0 (9)
对公式(9)积分后,得到
激波:
结合兰金-雨果关系(Rankine-Hugoniot conditions),进一步写做
将公式(13)带入(14)得到
引入热力学参数比焓h=e+Pv和比容v,得到方程
已知单位质量混合流体内能的表达式
这里ρmr,er分别代表参考压力Pr下气液混合体的密度和内能,结合公式(17)、(18),得到
通过公式(16)消去压力项P*、PL,得到左侧激波波速
同样地,穿过右激波得到波速
这里L和R分别对应波的左右侧。
进一步地,所述步骤S2中判定模型的构建方法为:
波型判定:
穿过左右侧非线性波fL(ρm*,UL)=VL-V*和fR(ρm*,UR)=V*-VR,得到有关混合流体中间状态处密度ρm*代数方程
f(ρm*,UL,UR)≡fL(ρm*,UL)+fR(ρm*,UR)+ΔV=0 (22)
接下来给出不同波型下fL(ρm*,UL)和fR(ρm*,UR)的表达式;
对于左侧波,当出现激波时,基于兰金-雨果关系(Rankine-Hugoniotconditions),公式(13)写做
同时公式(14)转化为
这里QmL为混合物质量通量,将公式(23)带入(24)得到
联合公式(24)(25),消除质量通量QmL,得到fL表达式
进一步地,得到
将公式(27)转化为更一般形式
当左侧出现稀疏波时,由公式(10),得到公式fL满足
同样地,将其简化为
同理,可以得到右侧波fR满足公式,汇总左右侧代数方程,可以得到有关黎曼问题中中间状态处气液混合流体速度和密度的表达式
或者
接下来,定义
结合方程(22),得到以下波型判定公式
对于双稀疏波和双激波类型,满足fmin·fmax>0,既
结合方程(32),得到
最终,给出以下判定公式。
相较于判定公式(34),判定条件(39)在编程计算更加的方便快捷。
进一步地,所述步骤S3中三种黎曼求解器的构建方法为:
迭代黎曼求解器:
基于密度方程(22),利用牛顿迭代法得到下式
这里给出初值ρm* 0的计算公式
两种非迭代黎曼求解器:
第一个非迭代求解器为两稀疏波求解器,在激波出现时会出现计算误差,这里构建第二个非迭代求解器
进一步地,所述步骤S4中通过含有Total Variation Diminishing(TVD)形式的MUSCL-Hancoke格式,得到重建的流动变量,具体如下:
S4-1:引入斜率限制器
S4-2:流体变量重构
S4-3:进一步变量重构
进一步地,所述步骤S5具体为:
以及
这里j代表迭代次数。
S6-1:对流项
S6-2:Δt/2时步更新
S6-3:Δt时步更新
源项S中梯度(S0)i定义为
进一步地,所述步骤S7具体为:
S7-1:由于采用了显式二阶龙格-库塔离散法将S引入求解,稳定性约束不仅要包括对流部分的Courant-Friedrichs-Lewy(CFL)准则,而且还要包括源项的约束。由CFL得到:
S7-2:显式二阶龙格-库塔离散化约束:
有益效果:本发明与现有技术相比,具备如下优点:
(1)本发明提供的基于耦合能量方程地有限体积法Godunov格式气-液两相均质流模型收敛到了正确结果并取得二阶精度,同时如MOC类方法一样简单且易于实现。
(2)本发明提出了一种精确地迭代型黎曼求解器和两种近似地非迭代黎曼求解器,相较于精确地黎曼求解器,两种非迭代求解器大幅提高了计算效率。
(3)与现有技术的二阶精度模型相比,本发明提供的模型在相同条件下能更为精确地拟合计算结果。
(4)能量方程的引入使得该均质流模型可以从能量角度分析管道中的瞬态过程,为以后的管道瞬变流研究提供了新的视角。
附图说明
图1为本发明的基本流程图;
图2为实施例对比实验装置示意图;
图3为实施例下,实验装置阀门末端(x=30.6)的压力曲线图;
图4为实施例下,网格数Nx=100,Crmax=0.95时,管道1号压力传感器处(x=8m),当前模型、León模型以及实验结果压力曲线对比图;
图5为实施例下,网格数Nx=100,Crmax=0.95时,实验管道2号压力传感器处(x=21.1m),当前模型、León模型以及实验结果压力曲线对比图。
具体实施方式
下面结合附图和具体实施例,进一步阐明本发明,应理解这些实施例仅用于说明本发明而不用于限制本发明的范围,在阅读了本发明之后,本领域技术人员对本发明的各种等价形式的修改均落于本申请所附权利要求所限定的范围。
本发明提供一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,如图1所示,其包括如下步骤:
S1:通过引入能量方程构建含自由气体均质两相流结合能量方程的3×3双曲偏微分方程(PDE)***模型,根据模拟工况确定计算域、初始条件以及边界条件:
通过在双曲偏微分方程(PDE)***体系下,构建含自由气体的两相均质流基本控制方程。
需在水锤问题的基础上假定:(a)管道内自由气体的体积占比很小(<≈1%)且分散均匀,气-水两相流以等效单相流体处理;(b)瞬变过程时间尺度很小,整个过程中水气间的质量传递过程被忽略;(c)管道内流体为无粘性流动且整个过程发生在等温条件下不考虑热传递。在上述条件的基础上,含自由气体的两相均质流基本控制方程具体如下:
包含能量方程的水锤基本方程:
其中,矩阵代表气液混合流体的特征变量;矩阵代表通量;沿管线距离x与时间t是自变量;ρm(x,t)是平均截面气液混合流体的密度;V(x,t)是气液混合体平均截面速度;P(x,t)是气液混合流体的绝对压力;E(x,t)是气液混合流体单位体积总能量;
气-液两相流水锤波速与压力的关系方程:
气-液两相流压力与密度的关系方程:
确定计算域、初始条件及边界条件:
计算域:从上游水库出水口至阀门之间的管道;
初始条件:球阀全开时,初始流速为设定值,各节点初始压力根据上游水库静压水头减去相应的稳态摩阻损失得到;
边界条件为:管道入口处,水库提供恒压边界,且等于上游水库静压水头;下游球阀快速关闭引发空穴流瞬变,此处压力变化由压力传感器收集得到。
S2:构建出3×3双曲偏微分方程***模型下各基础波(接触波、激波、稀疏波)的存在条件以及判定模型:
接触波:
对于两相流***中的接触波关系,可以利用欧拉方程中的特征结构分析断点处的初值,利用一般黎曼不变量,穿过波λ=V,得到,
这里Pr代表参考压力;u1、u2、u3见公式(2);ρmr为参考压力下气液混合流体的密度;cm为气液混合流体的波速;er为参考压力下气液混合体的内能;进一步简化后,黎曼不变量可以写做:
这里星号(*)代表中间状态,P、V、ρm分别代表混合流体的压力、速度、密度,L和R分别代表黎曼问题中接触波左右区间,在等熵条件下,穿过接触波,流体特征不发生变化;
稀疏波:
相较于激波,稀疏波表现的更为平滑,穿过波λ=V-cm,利用一般黎曼不变量,得到
进一步简化后,该微分方程可以写做
dy+cmd(lnρm)=0 (9)
对公式(9)积分后,得到
激波:
结合兰金-雨果关系(Rankine-Hugoniot conditions),进一步写做
将公式(13)带入(14)得到
引入热力学参数比焓h=e+Pv和比容v,得到方程
已知单位质量混合流体内能的表达式
这里ρmr,er分别代表参考压力Pr下气液混合体的密度和内能,结合公式(17)、(18),得到
通过公式(16)消去压力项P*、PL,得到左侧激波波速
同样地,穿过右激波得到波速
这里L和R分别对应波的左右侧。
判定模型的构建方法为:
波型判定:
穿过左右侧非线性波fL(ρm*,UL)=VL-V*和fR(ρm*,UR)=V*-VR,得到有关混合流体中间状态处密度ρm*代数方程
f(ρm*,UL,UR)≡fL(ρm*,UL)+fR(ρm*,UR)+ΔV=0 (22)接下来给出不同波型下fL(ρm*,UL)和fR(ρm*,UR)的表达式;
对于左侧波,当出现激波时,基于兰金-雨果关系(Rankine-Hugoniotconditions),公式(13)写做
同时公式(14)转化为
这里QmL为混合物质量通量,将公式(23)带入(24)得到
联合公式(24)(25),消除质量通量QmL,得到fL表达式
进一步地,得到
将公式(27)转化为更一般形式
当左侧出现稀疏波时,由公式(10),得到公式fL满足
同样地,将其简化为
同理,可以得到右侧波fR满足公式,汇总左右侧代数方程,可以得到有关黎曼问题中中间状态处气液混合流体速度和密度的表达式
或者
接下来,定义
结合方程(22),得到以下波型判定公式
对于双稀疏波和双激波类型,满足fmin·fmax>0,既
结合方程(32),得到
最终,给出以下判定公式。
S3:构建双曲偏微分方程***模型下断点处三种近似黎曼求解器,分别为一种迭代和两种非迭代近似黎曼求解器:
迭代黎曼求解器:
基于密度方程(22),利用牛顿迭代法得到下式
这里给出初值ρm* 0的计算公式
f(ρm*)=fL+fR+Δy= (41)
两种非迭代黎曼求解器:
第一个非迭代求解器为两稀疏波求解器,在激波出现时会出现计算误差,这里构建第二个非迭代求解器
S4:通过含有Total Variation Diminishing(TVD)形式的MUSCL-Hancoke格式,得到重建的流动变量,使其具有二阶精度:
S4-1:引入斜率限制器
S4-2:流体变量重构
S4-3:进一步变量重构
S5:根据已有的边界条件,通过牛顿迭代法利用黎曼不变量迭代得到二阶的边界值:
以及
这里j代表迭代次数。
S6:引入源项,通过基于二阶Longge kuta离散格式的时间***法,求得时间梯度上二阶有限体积法Godunov格式离散方程,得到二阶精度:
S6-1:对流项
S6-2:Δt/2时步更新
S6-3:Δt时步更新
源项S中梯度(S0)i定义为
S7:建立稳定性约束条件,更新初始变量进行下一时步计算,获取到最大允许时间步长:
S7-1:由于采用了显式二阶龙格-库塔离散法将S引入求解,稳定性约束不仅要包括对流部分的Courant-Friedrichs-Lewy(CFL)准则,而且还要包括源项的约束。由CFL得到:
S7-2:显式二阶龙格-库塔离散化约束:
本发明还提供一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟***,该***包括网络接口、存储器和处理器;其中,网络接口,用于在与其他外部网元之间进行收发信息过程中,实现信号的接收和发送;存储器,用于存储能够在所述处理器上运行的计算机程序指令;处理器,用于在运行计算机程序指令时,执行上述共识方法的步骤。
本发明还提供一种计算机存储介质,该计算机存储介质存储有计算机程序,在处理器执行所述计算机程序时可实现以上所描述的方法。所述计算机可读介质可以被认为是有形的且非暂时性的。非暂时性有形计算机可读介质的非限制性示例包括非易失性存储器电路(例如闪存电路、可擦除可编程只读存储器电路或掩膜只读存储器电路)、易失性存储器电路(例如静态随机存取存储器电路或动态随机存取存储器电路)、磁存储介质(例如模拟或数字磁带或硬盘驱动器)和光存储介质(例如CD、DVD或蓝光光盘)等。计算机程序包括存储在至少一个非暂时性有形计算机可读介质上的处理器可执行指令。计算机程序还可以包括或依赖于存储的数据。计算机程序可以包括与专用计算机的硬件交互的基本输入/输出***(BIOS)、与专用计算机的特定设备交互的设备驱动程序、一个或多个操作***、用户应用程序、后台服务、后台应用程序等。
本领域内的技术人员应明白,本申请的实施例可提供为方法、***、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
具体实施例方案:
基于上述方案,为了验证并分析本发明提供的基于有限体积法耦合能量方程的输水管道内水气两相均质流模型的模拟效果,本实施例选取Chaudry于1990年设计搭建的均质流实验装置***用于验证本发明方法的有效性。如图2所示,整个***由上游水箱,管道,下游球阀,下游水箱组成。管道总长30.6m,内径0.026mm,管道水平。达西韦斯巴赫系数为0.0195,纯水水锤波速为715m/s。瞬变过程由突然关闭下游球阀引起。实验工况参数为:初始流速2.940m/s,上游水库静压水头21.700m,稳定时气体质量流动率1.15×10-5kg/s,下游水体中气体体积含量0.0053。
通过以上方法编程计算后,将结合能量方程的Godunov格式模型计算结果与实验数据作对比。图3给出了实验装置阀门末端(x=30.6)的压力曲线图,图4和5分别给出了本实施例下网格数N=200,库朗特数Crmax=0.95时,压力传感器1、2处的压力曲线图。可以看出,基于本发明模型下的预测值能够很好的拟合实验检测值,本发明的有效性以及准确性得到了很好的验证。同时,针对输水管道中气-液两相均质流,相比MOC类方法以及现有文献中的守恒形式的Godunov格式方法,本发明提出的非守恒有限体积法Godunov方法,在模型预测方面更加的准确且更加的有效率。
分析总结:
本发明方法基于一个考虑能量效应的3×3双曲偏微分方程(PDE)***来描述均匀两相流的不稳定和可压缩行为;然后运用一种二阶精度的Godunov格式来实现偏微分方程的离散化,同时引入了一个精确的和两个非迭代的近似黎曼求解器来求解不连续处的黎曼问题,接着,利用数值测试对不同求解器进行比较,结果表明非迭代近似求解器可以更好地提高计算效率,两个非迭代近似求解器分别提高了64.28%和43.59%的计算效率。通过与以往模型和已发表文章中的实验测量数据进行比较来评估所提出模型的性能。与之前的保守型Godunov模型相比,本发明提出的方法在计算精度上表现更好,可以更准确地拟合实验曲线。同时能量方程的引入为研究两相流瞬态过程提供了一个新的视角,结果表明较大的柯朗数和较细的网格可以有效地减少***耗散。
Claims (10)
1.一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,其特征在于,包括如下步骤:
S1:通过引入能量方程构建含自由气体均质两相流结合能量方程的双曲偏微分方程***模型,根据模拟工况确定计算域、初始条件以及边界条件;
S2:构建出双曲偏微分方程***模型下各基础波的存在条件以及判定模型;
S3:构建双曲偏微分方程***模型下断点处三种近似黎曼求解器,分别为一种迭代和两种非迭代近似黎曼求解器;
S4:通过重建流动变量,使得***模型在空间上达到二阶精度;
S5:根据已有的边界条件,通过牛顿迭代法利用黎曼不变量迭代得到二阶的边界值;
S6:引入源项,通过基于二阶Longgekuta离散格式的时间***法,求得时间梯度上二阶有限体积法Godunov格式离散方程,得到二阶精度;
S7:建立稳定性约束条件,更新初始变量进行下一时步计算,获取到最大允许时间步长。
2.根据权利要求1所述的一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,其特征在于,所述步骤S1中双曲偏微分方程***模型的构建是通过在双曲偏微分方程***体系下,构建含自由气体的两相均质流基本控制方程。
3.根据权利要求2所述的一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,其特征在于,所述步骤S1中含自由气体的两相均质流基本控制方程具体如下:
包含能量方程的水锤基本方程:
其中,矩阵代表气液混合流体的特征变量;矩阵代表通量;沿管线距离x与时间t是自变量;ρm(x,t)是平均截面气液混合流体的密度;V(x,t)是气液混合体平均截面速度;P(x,t)是气液混合流体的绝对压力;E(x,t)是气液混合流体单位体积总能量;
气-液两相流水锤波速与压力的关系方程:
其中,cm为气液混合流体的波速;C为纯液体下水锤波速;Pref为参考压力;ρfref为参考压力下气液两相混合流体的密度;ψref为参考压力下混合流体内气体的初始体积分数;
气-液两相流压力与密度的关系方程:
其中,ρm为气液混合流体的密度,P为气液混合流体的压力;C为纯液体下水锤波速;Pref为参考压力;ρfref为参考压力下气液两相混合流体的密度;ψref为参考压力下混合流体内气体的初始体积分数。
4.根据权利要求1所述的一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,其特征在于,所述步骤S2中基础波包括接触波、激波、稀疏波,构建方法具体如下:
接触波:
对于两相流***中的接触波关系,利用欧拉方程中的特征结构分析断点处的初值,利用一般黎曼不变量,穿过波λ=V,得到,
这里Pr代表参考压力;u1、u2、u3见公式(2);ρmr为参考压力下气液混合流体的密度;cm为气液混合流体的波速;er为参考压力下气液混合体的内能;进一步简化后,黎曼不变量可以写做:
这里星号(*)代表中间状态,P、V、ρm分别代表混合流体的压力、速度、密度,L和R分别代表黎曼问题中接触波左右区间,在等熵条件下,穿过接触波,流体特征不发生变化;
稀疏波:
穿过波λ=V-cm,利用一般黎曼不变量,得到
进一步简化后,该微分方程可以写做
dV+cmd(lnρm)=0 (9)
对公式(9)积分后,得到
激波:
结合兰金-雨果关系,进一步写做
将公式(13)带入(14)得到
引入热力学参数比焓h=e+Pv和比容v,得到方程
已知单位质量混合流体内能的表达式
这里ρmr,er分别代表参考压力Pr下气液混合体的密度和内能,结合公式(17)、(18),得到
通过公式(16)消去压力项P*、PL,得到左侧激波波速
同样地,穿过右激波得到波速
这里L和R分别对应波的左右侧。
5.根据权利要求4所述的一种基于有限体积法耦合能量方程的输水管道气-液两相均质流的模拟方法,其特征在于,所述步骤S2中判定模型的构建方法为:
波型判定:
穿过左右侧非线性波fL(ρm*,UL)=VL-V*和fR(ρm*,UR)=V*-VR,得到有关混合流体中间状态处密度ρm*代数方程
f(ρm*,UL,UR)≡fL(ρm*,UL)+fR(ρm*,UR)+ΔV=0 (22)接下来给出不同波型下fL(ρm*,UL)和fR(ρm*,UR)的表达式;
对于左侧波,当出现激波时,基于兰金-雨果关系,公式(13)写做
同时公式(14)转化为
这里QmL为混合物质量通量,将公式(23)带入(24)得到
联合公式(24)(25),消除质量通量QmL,得到fL表达式
进一步地,得到
将公式(27)转化为更一般形式
当左侧出现稀疏波时,由公式(10),得到公式fL满足
同样地,将其简化为
同理,可以得到右侧波fR满足公式,汇总左右侧代数方程,可以得到有关黎曼问题中中间状态处气液混合流体速度和密度的表达式
或者
接下来,定义
结合方程(22),得到以下波型判定公式
对于双稀疏波和双激波类型,满足fmin·fmax>0,既
结合方程(32),得到
最终,给出以下判定公式
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110705741.XA CN113435136B (zh) | 2021-06-24 | 2021-06-24 | 耦合能量方程的输水管道气-液两相均质流的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110705741.XA CN113435136B (zh) | 2021-06-24 | 2021-06-24 | 耦合能量方程的输水管道气-液两相均质流的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113435136A true CN113435136A (zh) | 2021-09-24 |
CN113435136B CN113435136B (zh) | 2023-06-20 |
Family
ID=77754288
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110705741.XA Active CN113435136B (zh) | 2021-06-24 | 2021-06-24 | 耦合能量方程的输水管道气-液两相均质流的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113435136B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114398844A (zh) * | 2022-01-25 | 2022-04-26 | 南京航空航天大学 | 基于连续水膜流动的稳态防冰仿真方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105512363A (zh) * | 2015-11-25 | 2016-04-20 | 河海大学 | 基于Godunov格式的有压管道中水柱分离的模拟方法 |
CN106777770A (zh) * | 2017-01-09 | 2017-05-31 | 河海大学 | 基于有限体积法的输水管道中空穴流的模拟方法 |
CN109918787A (zh) * | 2019-03-08 | 2019-06-21 | 河海大学 | 基于有限体积法的输水管道内水气两相均质流的模拟方法 |
CN111339701A (zh) * | 2020-02-26 | 2020-06-26 | 河海大学 | 基于Brunone动态摩阻的管道泄漏特性Godunov模拟方法 |
CN111695307A (zh) * | 2020-05-21 | 2020-09-22 | 河海大学 | 显式考虑动态摩阻的水锤有限体积模拟方法 |
CN111985166A (zh) * | 2020-07-31 | 2020-11-24 | 河海大学 | 隐式考虑动态摩阻的管道水力瞬变模拟方法和存储介质 |
-
2021
- 2021-06-24 CN CN202110705741.XA patent/CN113435136B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105512363A (zh) * | 2015-11-25 | 2016-04-20 | 河海大学 | 基于Godunov格式的有压管道中水柱分离的模拟方法 |
CN106777770A (zh) * | 2017-01-09 | 2017-05-31 | 河海大学 | 基于有限体积法的输水管道中空穴流的模拟方法 |
CN109918787A (zh) * | 2019-03-08 | 2019-06-21 | 河海大学 | 基于有限体积法的输水管道内水气两相均质流的模拟方法 |
CN111339701A (zh) * | 2020-02-26 | 2020-06-26 | 河海大学 | 基于Brunone动态摩阻的管道泄漏特性Godunov模拟方法 |
CN111695307A (zh) * | 2020-05-21 | 2020-09-22 | 河海大学 | 显式考虑动态摩阻的水锤有限体积模拟方法 |
CN111985166A (zh) * | 2020-07-31 | 2020-11-24 | 河海大学 | 隐式考虑动态摩阻的管道水力瞬变模拟方法和存储介质 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114398844A (zh) * | 2022-01-25 | 2022-04-26 | 南京航空航天大学 | 基于连续水膜流动的稳态防冰仿真方法 |
CN114398844B (zh) * | 2022-01-25 | 2023-04-07 | 南京航空航天大学 | 基于连续水膜流动的稳态防冰仿真方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113435136B (zh) | 2023-06-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109918787B (zh) | 基于有限体积法的输水管道内气液两相均质流的模拟方法 | |
Younis et al. | Adaptively localized continuation-Newton method—nonlinear solvers that converge all the time | |
Zhou et al. | Investigation of hydraulic transients of two entrapped air pockets in a water pipeline | |
Page et al. | Pressure management of water distribution systems via the remote real-time control of variable speed pumps | |
Zhou et al. | Influence of entrapped air pockets on hydraulic transients in water pipelines | |
Zhou et al. | Phenomenon of white mist in pipelines rapidly filling with water with entrapped air pockets | |
Wu et al. | MOC-CFD coupled approach for the analysis of the fluid dynamic interaction between water hammer and pump | |
Zhang et al. | A damped iterative EnKF method to estimate relative permeability and capillary pressure for tight formations from displacement experiments | |
Voskov | An extended natural variable formulation for compositional simulation based on tie-line parameterization | |
Leon et al. | Efficient second-order accurate shock-capturing scheme for modeling one-and two-phase water hammer flows | |
Yang et al. | Dynamic analysis of the pump system based on MOC–CFD coupled method | |
Sarma et al. | Efficient and robust uncertainty quantification in reservoir simulation with polynomial chaos expansions and non-intrusive spectral projection | |
CN113435136A (zh) | 耦合能量方程的输水管道气-液两相均质流的模拟方法 | |
Pini et al. | Robust adjoint-based shape optimization of supersonic turbomachinery cascades | |
Vasconcelos et al. | Innovative simulation of unsteady low-pressure flows in water mains | |
Hu et al. | Application of implicit Roe-type scheme and Jacobian-Free Newton-Krylov method to two-phase flow problems | |
Todini et al. | Pressure flow–based algorithms for pressure-driven analysis of water distribution networks | |
Zhang | Simulation of the transient flow in a natural gas compression system using a high-order upwind scheme considering the real-gas behaviors | |
Moosavian | Multilinear method for hydraulic analysis of pipe networks | |
Zheng et al. | A multidimensional and multiscale model for pressure analysis in a reservoir-pipe-valve system | |
Nault et al. | Adaptive hybrid transient formulation for simulating incompressible pipe network hydraulics | |
Dong et al. | A new method for predicting CO2 minimum miscibility pressure MMP based on deep learning | |
Fuls | Enhancement to the traditional ellipse law for more accurate modeling of a turbine with a finite number of stages | |
Hanif Chaudhry et al. | Mathematical modeling of lake tap flows | |
Osigwe et al. | Performance Simulation to Understand the Effects of Multifluid Scaling of Gas Turbine Components for Generation IV Nuclear Power Plants |
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 |