CN106844964B - 一种压力管道输水系统非恒定流模型的优化方法 - Google Patents
一种压力管道输水系统非恒定流模型的优化方法 Download PDFInfo
- Publication number
- CN106844964B CN106844964B CN201710048311.9A CN201710048311A CN106844964B CN 106844964 B CN106844964 B CN 106844964B CN 201710048311 A CN201710048311 A CN 201710048311A CN 106844964 B CN106844964 B CN 106844964B
- Authority
- CN
- China
- Prior art keywords
- pressure
- water
- pipe section
- wave
- pressure pipe
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000005457 optimization Methods 0.000 title claims abstract description 9
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 99
- 239000004071 soot Substances 0.000 claims abstract description 23
- 230000011218 segmentation Effects 0.000 claims abstract description 19
- 238000004364 calculation method Methods 0.000 claims description 9
- 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
- 238000007796 conventional method Methods 0.000 abstract description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000011438 discrete method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pipeline Systems (AREA)
Abstract
本发明公开了一种压力管道输水系统非恒定流模型的优化方法,其特征在于,涉及水利工程技术领域。所述方法包括:获取压力管道输水系统中的有压管段,对水锤波传播时间最短有压管段的水锤波速进行调整,获取最大时间步长Δtmax,计算获取最大时间步长Δtmax所对应的最短有压管段的最小分段数N0min;通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优;通过最小波速调整幅度fimin,调整水锤波速,保证调整前、后水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同。本发明使压力管道输水系统的波速调整幅度与比传统方法相比显著减小。
Description
技术领域
本发明涉及水利工程技术领域,尤其涉及一种压力管道输水系统非恒定流模型的优化方法。
背景技术
长距离输水工程的特点为输水管线长、地形起伏多变和水流条件复杂等,其过渡过程模拟计算的影响因素比较多,除了输水系统自身影响因素外,自身影响因素如摩阻、压力、温度等,计算模型因素包括时间步长选取、采用的计算方法等也至关重要。长距离输水管道瞬变流计算时,时间步长的取值大小是影响水力过渡过程计算结果的重要因素之一,时间步长取值太大,无法满足计算精度要求,方程的解难以达到稳定;取值太小,若输水管线长,导致计算工作量较大。因此合理确定计算管道分段数,既保证计算的精度,又减少计算的工作量,是长距离管道输水泵站工程的水力过渡过程数值模拟必须考虑的一个重要因素。
现有通常调整波速法计算时间步长,调整波速法是通过在一定范围内对波速进行调整,使压力管道输水系统内各管段的分段数是整数且时间步长相等。该法是目前最为常用的压力管道(复杂管系)分段方法。
调整波速法存在的问题如下:
(1)由于调整波速法是依据压力管道输水系统内波传播时间最短的一条管段(通常为管长最短的管段)作为计算基础,在计算时间步长时,不对短管的波速并没有调整,使得压力管道输水系统的管道分段受到限制,有可能使时间步长Δt的取值过小。
(2)波速调整后,水锤波从压力管道输水系统起点传播到终点的总时间T′与未调整波速时的总时间T很难一致,存在一个差值ΔT=T′-T引起误差。
(3)现有调整波速法没有体现优化思想,分段时只要满足公式a就算满足要求。
尽管调整波速分段法的提出主要是基于水锤波速影响因素多,不能精确得知时间步长,得到的时间步长与真实值间存在误差。
发明内容
本发明的目的在于提供一种压力管道输水系统非恒定流模型的优化方法,从而解决现有技术中存在的前述问题。
为了实现上述目的,本发明所述压力管道输水系统非恒定流模型的优化方法,所述方法包括:
首先,获取压力管道输水系统中的有压管段,对水锤波传播时间最短有压管段的水锤波速进行调整,获取最大时间步长Δtmax,计算获取最大时间步长Δtmax所对应的最短有压管段的最小分段数N0min;
然后,通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优;
最后,通过最小波速调整幅度fimin,调整水锤波速,保证调整前、后水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同。
优选地,对水锤波传播时间最短有压管段的水锤波速进行调整,获取最大时间步长Δtmax,计算获取最大时间步长Δtmax所对应的最短有压管段最小分段数,具体按照下述步骤实现:
计算压力管道输水系统内水锤波传播时间最短有压管段的管长L0、波速a0和分段数N0,其中,管长L0、波速a0存在关系式(1):
压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段分段数Ni满足公式(2),所述Ni为正整数:
其中,i表示压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的编号,i为大于等于1的正整数,m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量,α表示压力管道输水系统中除水锤波传播时间最短有压管段外其他有压管段的水锤波波速;fi表示编号为i的有压管段的波速调整系数;Li表示编号为i的有压管段的长度;△t表示水锤波传播时间最短管段的关系时间步长;
所述△t的取值范围符合公式(3):
fmax表示|fi|的最大取值,在计算水锤波传播时间最短有压管段的波速调整的关系时间步长时,|fi|=0.15;T表示未调整波速时,水锤波从压力管道输水系统起点传播到终点的实际总时间;
由于水锤波传播时间最短有压管段的波速调整不会影响对压力管道输水系统的计算结果,故,在获取△t的取值范围时,设|fi|=0.15,则fmax=-0.15或fmax=0.15,因此,△t的取值范围符合公式(4):
通过公式(4),得到最大关系时间步长Δtmax的计算公式(5);
在公式(5)的基础上,根据公式(6)计算最大的关系时间步长Δtmax所对应的最小分段数N0min为:
优选地,通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优,具体按照下述步骤实现:
S1,采用公式(7)计算任意一个有压管段i的关系时间步长Δtj;
△tj=L0/(N0mina0(1+f0)),-0.15≤f0≤0.15 (7);
L0、a0分别表示压力管道输水系统内水锤波传播时间最短有压管段的管长和波速;f0表示压力管系内波传播时间最短管段的波速调整系数;当f0在[-0.15,0.15]范围内以一微小增量j改变时,由公式(7)计算当f0=f0j时,与f0j对应的Δtj值,所述微小增量j的取值增量为0.001;
S2,根据关系时间步长Δtj,计算压力管道输水系统内编号为i的有压管段的分段数Ni,i=1,2,…,m,m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量;判断分段数Ni是否为式(8)所示的范围内的正整数;
如果是,进入S3,计算得到编号为i的有压管道的最小波速调整幅度fimin:
如果否,表示f0=f0j求得的关系时间步长Δtj不符合要求,则返回S1,应用公式(7)计算当f0=f0(j+1)时,与f0(j+1)所对应的△t(j+1),继续S2,判断通过f0=[-0.15,0.15]中的所有取值是否能够计算得到符合式(8)的正整数分段数;如果是,则进入S3;如果否,则返回S1,并将公式(7)中的N0min修改为N0min+1,继续计算,直至计算得到编号为i的管段至少存在一个符合式(8)的正整数分段数为止;
S3,将符合式(8)的编号为i的有压管道分段数按下式(9)进行计算,得到编号为i的有压管道的最小波速调整幅度fimin;
fimin所对应的分段数即为有压管段i的最优分段数Ni最优,Li表示有压管段i的管长,ai表示管段i的波速。
优选地,通过最小波速调整幅度fimin,调整水锤波速,保证调整前、后,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同,具体按照下述步骤实现:
水锤波速未调整,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的实际总时间为T,水锤波速调整后,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间为T′,T与T′的关系用公式(10)表示:
T'ε=T (10);
其中,系数ε符合公式(11);
则,最优波速调整系数为f′i:f′i=(1+fimin)/ε-1 (12);
最优时间步长为Δt′:△t′=△tjε (13);
最后,判断压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的最优波速调整系数为fi′是否小于0.10,如果是,则该有压管段的最优时间步长为Δt′,如果否,则将水锤波传播时间最短有压管段的分段数加1后,重新计算最大的时间步长Δtmax及最短有压管段的最小分段数N0min,然后进入S2,直至计算得到的所有结果均满足条件为止。
优选地,调整水锤波速调整前、后的总时间相同后,还包括以下步骤:
计算与有压管道相邻的无压渠道的时间步长;设无压管道的时间步长为△t0,管道流的时间步长为△tc,令△t0与△tc符合公式(14):
△tc=N△t0,N为整数 (14);
设明渠和管道流的空间步长分别为△x0与△xc,为保证明渠与管道流的衔接,将明渠在连接断面附近的空间步长△x0设置为△x0/N,然后按库朗条件将明渠在连接断面的△t0缩小1/N×△t0,从而与管道的时间步长△tc相匹配。
本发明是一种通过改进调整波速法压力管道确定时间步长的压力管道输水系统非恒定流模型的优化技术,采用改进的调整波速法进行分段,首先对波传播时间最短管段的波速进行调整,选取尽可能大的时间步长Δt,最小化整个管系波速调整幅度,再通过各管段最优波速调整系数,使波速调整前后水锤波从压力管道起点传播到终点的总时间保持一致。
本发明的有益效果是:本发明所述方法使压力管道输水系统的波速调整幅度与比传统方法相比显著减小;通过ε的调整,使时差ΔT=0s,保证了波速调整前后水锤波从压力管系起点传播到终点的总时间完全一致。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不用于限定本发明。
实施例
本实施例所述压力管道输水系统非恒定流模型的优化方法,所述方法包括:
首先,获取压力管道输水系统中的有压管段,对水锤波传播时间最短有压管段的水锤波速进行调整,获取最大时间步长Δtmax,计算获取最大时间步长Δtmax所对应的最短有压管段的最小分段数N0min;
然后,通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优;
最后,通过最小波速调整幅度fimin,调整水锤波速,保证调整前、后水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同。
更详细的解释说明:
(一)对水锤波传播时间最短有压管段的水锤波速进行调整,获取最大时间步长Δtmax,计算获取最大时间步长Δtmax所对应的最短有压管段最小分段数,具体按照下述步骤实现:
计算压力管道输水系统内水锤波传播时间最短有压管段的管长L0、波速a0和分段数N0,其中,管长L0、波速a0存在关系式(1):
压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段分段数Ni满足公式(2),所述Ni为正整数:
其中,i表示压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的编号,i为大于等于1的正整数,m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量,α表示压力管道输水系统中除水锤波传播时间最短有压管段外其他有压管段的水锤波波速;fi表示编号为i的有压管段的波速调整系数;Li表示编号为i的有压管段的长度;△t表示水锤波传播时间最短管段的关系时间步长;
所述△t的取值范围符合公式(3):
fmax表示|fi|的最大取值,在计算水锤波传播时间最短有压管段的波速调整的关系时间步长时,|fi|=0.15;T表示未调整波速时,水锤波从压力管道输水系统起点传播到终点的实际总时间;
由于水锤波传播时间最短有压管段的波速调整不会影响对压力管道输水系统的计算结果,故,在获取△t的取值范围时,设|fi|=0.15,则fmax=-0.15或fmax=0.15,因此,△t的取值范围符合公式(4):
通过公式(4),得到最大关系时间步长Δtmax的计算公式(5);
在公式(5)的基础上,根据公式(6)计算最大的关系时间步长Δtmax所对应的最小分段数N0min为:
(二)通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优,具体按照下述步骤实现:
S1,采用公式(7)计算任意一个有压管段i的关系时间步长Δtj;
△tj=L0/(N0mina0(1+f0)),-0.15≤f0≤0.15 (7);
L0、a0分别表示压力管道输水系统内水锤波传播时间最短有压管段的管长和波速;f0表示压力管系内波传播时间最短管段的波速调整系数;当f0在[-0.15,0.15]范围内以一微小增量j改变时,由公式(7)计算当f0=f0j时,与f0j对应的Δtj值,所述微小增量j的取值增量为0.001;
S2,根据关系时间步长Δtj,计算压力管道输水系统内编号为i的有压管段的分段数Ni,i=1,2,…,m,m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量;判断分段数Ni是否为式(8)所示的范围内的正整数;
如果是,进入S3,计算得到编号为i的有压管道的最小波速调整幅度fimin:
如果否,表示f0=f0j求得的关系时间步长Δtj不符合要求,则返回S1,应用公式(7)计算当f0=f0(j+1)时,与f0(j+1)所对应的△t(j+1),继续S2,判断通过f0=[-0.15,0.15]中的所有取值是否能够计算得到符合式(8)的正整数分段数;如果是,则进入S3;如果否,则返回S1,并将公式(7)中的N0min修改为N0min+1,继续计算,直至计算得到编号为i的管段至少存在一个符合式(8)的正整数分段数为止;
S3,将符合式(8)的编号为i的有压管道分段数按下式(9)进行计算,得到编号为i的有压管道的最小波速调整幅度fimin;
fimin所对应的分段数即为有压管段i的最优分段数Ni最优,Li表示有压管段i的管长,ai表示管段i的波速。
(三)通过最小波速调整幅度fimin,调整水锤波速,保证调整前、后,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同,具体按照下述步骤实现:
水锤波速未调整,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的实际总时间为T,水锤波速调整后,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间为T′,T与T′的关系用公式(10)表示:
T'ε=T (10);
其中,系数ε符合公式(11);
则,最优波速调整系数为fi′:fi'=(1+fimin)/ε-1 (12);
最优时间步长为Δt′:△t′=△tjε (13);
最后,判断压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的最优波速调整系数为fi′是否小于0.10,如果是,则该有压管段的最优时间步长为Δt′,如果否,则将水锤波传播时间最短有压管段的分段数加1后,重新计算最大的时间步长Δtmax及最短有压管段的最小分段数N0min,然后进入S2,直至计算得到的所有结果均满足条件为止。
(四)调整水锤波速调整前、后的总时间相同后,还包括以下步骤:
计算与有压管道相邻的无压渠道的时间步长;设无压管道的时间步长为△t0,管道流的时间步长为△tc,令△t0与△tc符合公式(14):
△tc=N△t0,N为整数 (14);
设明渠和管道流的空间步长分别为△x0与△xc,为保证明渠与管道流的衔接,将明渠在连接断面附近的空间步长△x0设置为△x0/N,然后按库朗条件将明渠在连接断面的△t0缩小1/N×△t0,从而与管道的时间步长△tc相匹配。
通过采用本发明公开的上述技术方案,得到了如下有益的效果:本发明所述方法使压力管道输水系统的波速调整幅度与比传统方法相比显著减小;通过ε的调整,使时差ΔT=0s,保证了波速调整前后水锤波从压力管系起点传播到终点的总时间完全一致。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (4)
1.一种压力管道输水系统非恒定流模型的优化方法,其特征在于,所述方法包括:
首先,获取压力管道输水系统中的有压管段,对水锤波传播时间最短有压管段的水锤波速进行调整,获取最大时间步长Δtmax,计算获取最大时间步长Δtmax所对应的最短有压管段的最小分段数N0min;
然后,通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优;
最后,通过最小波速调整幅度fimin,调整水锤波速,保证调整前、后水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同;
对水锤波传播时间最短有压管段的水锤波速进行调整,获取最大时间步长Δtmax,计算获取最大时间步长Δtmax所对应的最短有压管段最小分段数,具体按照下述步骤实现:
计算压力管道输水系统内水锤波传播时间最短有压管段的管长L0、波速a0和分段数N0,其中,管长L0、波速a0存在关系式(1):
压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段分段数Ni满足公式(2),所述Ni为正整数:
其中,i表示压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的编号,i为大于等于1的正整数,m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量,α表示压力管道输水系统中除水锤波传播时间最短有压管段外其他有压管段的水锤波波速;fi表示编号为i的有压管段的波速调整系数;Li表示编号为i的有压管段的长度;△t表示水锤波传播时间最短管段的关系时间步长;
所述△t的取值范围符合公式(3):
fmax表示|fi|的最大取值,在计算水锤波传播时间最短有压管段的波速调整的关系时间步长时,|fi|=0.15;T表示未调整波速时,水锤波从压力管道输水系统起点传播到终点的实际总时间;
由于水锤波传播时间最短有压管段的波速调整不会影响对压力管道输水系统的计算结果,故,在获取△t的取值范围时,设|fi|=0.15,则fmax=-0.15或fmax=0.15,因此,△t的取值范围符合公式(4):
通过公式(4),得到最大关系时间步长Δtmax的计算公式(5);
在公式(5)的基础上,根据公式(6)计算最大的关系时间步长Δtmax所对应的最小分段数N0min为:
2.根据权利要求1所述方法,其特征在于,通过最小分段数N0min计算其他任意一个有压管段i的最小波速调整幅度fimin及最小波速调整幅度fimin所对应的最优分段数Ni最优,具体按照下述步骤实现:
S1,采用公式(7)计算任意一个有压管段i的关系时间步长Δtj;
△tj=L0/(N0mina0(1+f0)),-0.15≤f0≤0.15 (7);
L0、a0分别表示压力管道输水系统内水锤波传播时间最短有压管段的管长和波速;f0表示压力管系内波传播时间最短管段的波速调整系数;当f0在[-0.15,0.15]范围内以一微小增量j改变时,由公式(7)计算当f0=f0j时,与f0j对应的Δtj值,所述微小增量j的取值增量为0.001;
S2,根据关系时间步长Δtj,计算压力管道输水系统内编号为i的有压管段的分段数Ni,i=1,2,…,m,m表示压力管道输水系统中除水锤波传播时间最短有压管段外有压管段的总数量;判断分段数Ni是否为式(8)所示的范围内的正整数;
如果是,进入S3,计算得到编号为i的有压管道的最小波速调整幅度fimin:
如果否,表示f0=f0j求得的关系时间步长Δtj不符合要求,则返回S1,应用公式(7)计算当f0=f0(j+1)时,与f0(j+1)所对应的△t(j+1),继续S2,判断通过f0=[-0.15,0.15]中的所有取值是否能够计算得到符合式(8)的正整数分段数;如果是,则进入S3;如果否,则返回S1,并将公式(7)中的N0min修改为N0min+1,继续计算,直至计算得到编号为i的管段至少存在一个符合式(8)的正整数分段数为止;
S3,将符合式(8)的编号为i的有压管道分段数按下式(9)进行计算,得到编号为i的有压管道的最小波速调整幅度fimin;
fimin所对应的分段数即为有压管段i的最优分段数Ni最优,Li表示有压管段i的管长,ai表示管段i的波速。
3.根据权利要求2所述方法,其特征在于,通过最小波速调整幅度fimin,调整水锤波速,保证调整前、后,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间相同,具体按照下述步骤实现:
水锤波速未调整,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的实际总时间为T,水锤波速调整后,水锤波从压力管道输水系统起点传播到压力管道输水系统终点的总时间为T′,T与T′的关系用公式(10)表示:
T'ε=T (10);
其中,系数ε符合公式(11);
则,最优波速调整系数为fi′:fi'=(1+fimin)/ε-1 (12);
最优时间步长为Δt′:△t′=△tjε (13);
最后,判断压力管道输水系统中除水锤波传播时间最短有压管段外任意一个有压管段的最优波速调整系数为fi′是否小于0.10,如果是,则该有压管段的最优时间步长为Δt′,如果否,则将水锤波传播时间最短有压管段的分段数加1后,重新计算最大的时间步长Δtmax及最短有压管段的最小分段数N0min,然后进入S2,直至计算得到的所有结果均满足条件为止。
4.根据权利要求1所述方法,其特征在于,调整水锤波速调整前、后的总时间相同后,还包括以下步骤:
计算与有压管道相邻的无压管道的时间步长;设无压管道的时间步长为△t0,管道流的时间步长为△tc,令△t0与△tc符合公式(14):
△tc=N△t0,N为整数 (14);
设明渠和管道流的空间步长分别为△x0与△xc,为保证明渠与管道流的衔接,将明渠在连接断面附近的空间步长△x0设置为△x0/N,然后按库朗条件将明渠在连接断面的△t0缩小1/N×△t0,从而与管道的时间步长△tc相匹配。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710048311.9A CN106844964B (zh) | 2017-01-20 | 2017-01-20 | 一种压力管道输水系统非恒定流模型的优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710048311.9A CN106844964B (zh) | 2017-01-20 | 2017-01-20 | 一种压力管道输水系统非恒定流模型的优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106844964A CN106844964A (zh) | 2017-06-13 |
CN106844964B true CN106844964B (zh) | 2018-12-18 |
Family
ID=59121117
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710048311.9A Expired - Fee Related CN106844964B (zh) | 2017-01-20 | 2017-01-20 | 一种压力管道输水系统非恒定流模型的优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106844964B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108319745B (zh) * | 2017-12-18 | 2020-12-08 | 中国水利水电科学研究院 | 渠道非恒定流计算方法及装置 |
CN113705120B (zh) * | 2021-08-25 | 2023-09-22 | 山东省调水工程运行维护中心 | 一种调水工程泵后阀最优调控方案制定方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102337752A (zh) * | 2010-07-29 | 2012-02-01 | 河南省电力勘测设计院 | 输水管道等径准调压井 |
-
2017
- 2017-01-20 CN CN201710048311.9A patent/CN106844964B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102337752A (zh) * | 2010-07-29 | 2012-02-01 | 河南省电力勘测设计院 | 输水管道等径准调压井 |
Non-Patent Citations (1)
Title |
---|
长距离输水系统的非恒定流特性研究;万五一;《中国博士学位论文全文数据库工程科技II辑》;20050715;第36-56页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106844964A (zh) | 2017-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106844964B (zh) | 一种压力管道输水系统非恒定流模型的优化方法 | |
CN107702933B (zh) | 测定平原河网灌溉泵站水量-电量转换系数的装置及方法 | |
CN102183895B (zh) | 一种新型的排水管网水力学系统建模与控制方法 | |
CN110443407B (zh) | 一种提高供水串联管道瞬态水力模拟精度的节点流量优化分配方法 | |
CN106869918A (zh) | 海上油井产能测试实时调整方法 | |
CN106704163A (zh) | 一种水泵变频调速控制方法、装置及系统 | |
CN105115550A (zh) | 一种双压差气液流量在线测量装置及在线测量方法 | |
CN108182319A (zh) | 一种超声速一体化喷管设计方法 | |
CN102590557A (zh) | 一种可变径负压式微风速标定装置 | |
CN211015169U (zh) | 一种管道汇流液体温度的控制系统 | |
CN105436213B (zh) | 一种轧后冷却装置集管流量前馈设定方法 | |
CN110908414A (zh) | 一种管道汇流液体温度的控制系统及方法 | |
CN110705109A (zh) | 套筒调节阀等百分比节流窗口型线确定方法 | |
CN115790758A (zh) | 基于温度补偿的伽马传感器计数矫正方法 | |
CN106092178B (zh) | 提高测量精度的数据修正方法 | |
CN105091967B (zh) | 一种利用管道本身或管道内障碍物的流量测量方法 | |
CN206146471U (zh) | 宽量程线性可变面积环形孔板流量装置 | |
CN110260953A (zh) | 一种修正音速喷嘴的流出系数的方法 | |
CN105181049B (zh) | 管道内流体流量的测量方法、测量系统及辅助测量系统 | |
WO2021031217A1 (zh) | 一种基于科氏力质量流量计的湿气流量计量方法及其装置 | |
CN206919944U (zh) | 一种节能型喉部取压长径喷嘴装置 | |
CN117109676B (zh) | 一种压力损失小、抗流体扰动性强的超声水表流场设计 | |
CN105257518A (zh) | 一种蠕动泵以及蠕动泵的精确定量校准方法 | |
CN218955825U (zh) | 一种带自动标定功能的矩阵流量计 | |
AU2003236775A1 (en) | Method for measuring the flow rate in channel shafts |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181218 |