CN111695304B - 一种水位梯度的加权平均计算方法 - Google Patents
一种水位梯度的加权平均计算方法 Download PDFInfo
- Publication number
- CN111695304B CN111695304B CN202010383921.6A CN202010383921A CN111695304B CN 111695304 B CN111695304 B CN 111695304B CN 202010383921 A CN202010383921 A CN 202010383921A CN 111695304 B CN111695304 B CN 111695304B
- Authority
- CN
- China
- Prior art keywords
- grid
- downstream
- interface
- upstream
- water level
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 58
- 238000004364 calculation method Methods 0.000 title claims abstract description 40
- 238000011144 upstream manufacturing Methods 0.000 claims description 63
- 238000000034 method Methods 0.000 claims description 33
- 238000004088 simulation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000002301 combined effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000002123 temporal effect 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
- 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)
- 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)
- Measuring Volume Flow (AREA)
Abstract
本发明提供一种水位梯度的加权平均计算方法,包括以下具体步骤:确定当前要计算的网格结点编号i及其左、右的网格结点编号i‑1、i+1;获取网格i‑1、i、i+1的水位zs、水深h、流速u、过水断面面积A数据;计算网格i的左、右界面的水位梯度;计算网格i的左、右界面水位梯度的加权系数wL、wR;通过加权计算得到计算网格i的水位梯度;本发明可以根据不同需求而使用相应的水位梯度计算方法,具有稳定性好、精度高的优点。
Description
技术领域
本发明涉及水流运动数值模拟领域,具体的说是一种水位梯度的加权平均计算方法。
背景技术
在进行水流运动数值模拟时,数值计算方法的误差控制能力,直接决定一个数值模型优劣甚至成败。好的数值格式,能够兼顾计算效率与误差控制,并至少在某一个方面有其优势。事实上,数值计算方法的不断进步,是数值模拟技术发展的主要驱动力之一。水位梯度项,是一、二维水动力方程中的一项,它表征的是水流单元各表面正压力的综合效应。为了避免分拆造成的“和谐性”问题,许多水动力模型将水位梯度项作为源项进行处理。作为源项处理时,这一项的格式涉及到时间、空间上的构造。在隐式模型中,目前一般采用时间加权的顺风格式。而在显式模型中,目前常用的有顺风格式、中心差分以及修正类方法。中心差分无法抑制非物理震荡,顺风格式精度相对较低。已有的修正类方法,相对于顺风或者中心格式,在许多情况下有一定的改善,但仍有提高的空间,且其在某些情况如水跃流条件下,修正效果不佳。
发明内容
针对现有技术的不足,本发明提出了提供一种水位梯度的加权平均计算方法,可以根据不同需求而使用,具有稳定性好、精度高的优点。
本发明的一种水位梯度的加权平均计算方法,包括以下具体步骤:
S1、确定当前要计算的网格结点编号为i,其左、右网格结点编号为i-1、i+1;
S2、获取网格i-1、i、i+1的水位zs、水深h、流速u、过水断面面积A数据;
S3、计算网格i的左、右界面的水位梯度:
其中,x为空间坐标,下标L、R分别代表左、右;
S4、计算网格i的左、右界面水位梯度的加权系数wL、wR;
S5、通过加权计算得到计算网格i的水位梯度:
所述步骤S4中推求加权系数wL、wR有3种方法,分别为稳定性更好的方法A;精度更高的方法B;所需参数较少的方法C。
所述方法A的具体步骤为,
A1、通过网格i-1、i、i+1的流速判断流向,如果网格i、i+1的流速均大于零,则网格i+1为下游、网格i-1为上游;如果网格i、i-1的流速均小于零,则网格i-1为下游、网格i+1为上游;其它情况则不区分上下游;
A2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
A3、如果区分上下游,则先计算网格i左界面的移动速度,其计算方法为:
上式中,下标L和R是表示相对于所求的界面而不是网格的左、右;
所述方法B的具体步骤为,
B1、通过网格i-1、i、i+1的流速判断流向,如果网格i、i+1的流速均大于零,则网格i+1为下游、网格i-1为上游;如果网格i、i-1的流速均小于零,则网格i-1为下游、网格i+1为上游;其它情况则不区分上下游;
B2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
B3、如果区分上下游,则先计算网格左界面的移动速度,其计算方法为:
所述方法C的具体步骤为,
C1、通过网格i的流速判断流向,如果网格i的流速大于零,则网格i的右侧为下游、左侧为上游;如果网格i的流速小于零,则网格i的左侧为下游、右侧为上游;其它情况则不区分上下游;
C2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
C3、如果区分上下游,则先计算网格i的Fr数,其计算方法为:
C4、然后计算水位差Δzs,如果网格i的右侧为下游,则Δzs=zs,i-zs,i+1,如果网格i的左侧为下游,则Δzs=zs,i-zs,i-1;
C5、先计算上游界面的权重因子wup,如果水位差Δzs>0则wup=min(0.5,Fr2),否则wup=0;
C6、如果左界面位于网格i的上游侧,则左界面的权重系数wL=wup,如果左界面位于网格i的下游侧,则wL=1-wup;
C7、右界面的权重系数wR=1-wL。
本发明的技术效果:本发明可以根据不同需求而使用的水位梯度计算方法,具有稳定性好、精度高的优点。
附图说明
图1为本发明的方法流程示意图。
图2为本发明的计算网格示意图。
图3为本发明实施案例的计算初始水深分布图。
图4为本发明实施案例的计算结果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1至图4所示,一种水位梯度的加权平均计算方法,包括以下具体步骤:
S1、确定当前要计算的网格结点编号为i,其左、右网格结点编号为i-1、i+1;
S2、获取网格结点i-1、i、i+1的水位zs、水深h、流速u、过水断面面积A数据;
S3、计算网格i的左、右界面的水位梯度:
其中,x为空间坐标,下标L、R分别代表左、右;
S4、计算网格i的左、右界面水位梯度的加权系数wL、wR;
S5、通过加权计算得到计算网格i的水位梯度:
所述步骤S4中推求加权系数wL、wR有3种方法,分别为稳定性更好的方法A;精度更高的方法B;所需参数较少的方法C。
所述方法A的具体步骤为,
A1、通过网格i-1、i、i+1的流速判断流向,如果网格i、i+1的流速均大于零,则网格i+1为下游、网格i-1为上游;如果网格i、i-1的流速均小于零,则网格i-1为下游、网格i+1为上游;其它情况则不区分上下游;
A2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
A3、如果区分上下游,则先计算网格i左界面的移动速度,其计算方法为:
上式中,下标L和R是表示相对于所求的界面而不是网格的左、右;
所述方法B的具体步骤为,
B1、通过网格i-1、i、i+1的流速判断流向,如果网格i、i+1的流速均大于零,则网格i+1为下游、网格i-1为上游;如果网格i、i-1的流速均小于零,则网格i-1为下游、网格i+1为上游;其它情况则不区分上下游;
B2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
B3、如果区分上下游,则先计算网格i的左界面的移动速度,其计算方法为:
所述方法C的具体步骤为,
C1、通过网格i的流速判断流向,如果网格i的流速大于零,则网格i的右侧为下游、左侧为上游;如果网格i的流速小于零,则网格i的左侧为下游、右侧为上游;其它情况则不区分上下游;
C2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
C3、如果区分上下游,则先计算网格i的Fr数,其计算方法为:
C4、然后计算水位差Δzs,如果网格i的右侧为下游,则Δzs=zs,i-zs,i+1,如果网格i的左侧为下游,则Δzs=zs,i-zs,i-1;
C5、先计算上游界面的权重因子wup,如果水位差Δzs>0则wup=min(0.5,Fr2),否则wup=0;
C6、如果左界面在网格i的上游侧,则左界面的权重系数wL=wup,如果左界面在网格i的下游侧,则wL=1-wup;
C7、右界面的权重系数wR=1-wL。
下面通过一个案例,并结合附图,对本发明的具体实施方式,做具体说明。
案例:光滑、平底湿河床薄壁坝瞬时全溃
案例简介:平底、光滑、湿河床薄壁瞬时溃坝。河道总长L=1000m,计算时长T=10s。网格步长Δx=5m,时间步长Δt=0.01s。
利用本发明提出的水位梯度计算方法对流场进行计算的步骤为:
根据网格步长,确定网格结点的总数量N及各网格结点的空间位置xi(i=1,...,N);
根据初始条件,给定各网格点的初始流场单宽流量Qi、水深Hi、流速Ui,其中i=1,...,N;
选择合适的物理方程,构建数学物理模型,本案例采用一维圣维南方程组作为数学模型,其形式如下:
式中,t、x为时间、一维空间坐标,zs为水位,g为重力加速度,n为曼宁糙率系数,其它如前所述;
对一维圣维南方程组进行离散,其离散结果可整理为如下形式:
其中,下标i-1、i、i+1为结点编号,上标n为流场情况已知的当前时刻时间步,上标n+1为待求的下一个时间步。
应用本发明的水位梯度计算方法,计算方程(4)中右边的水位梯度项;
通过方程(3)和(4)的反复利用,实现计算在时间域的步进,从而完整流场求解。
计算结果如图4,为便于对比,图4还展示了采用顺风格式方法计算水位梯度项的水深模拟结果。经统计,顺风格式、方法A、方法B、方法C四种方法计算结果与理论解的相关系数分别为0.99921、0.99904、0.99926、0.99956。可见,在本案例中,本发明所提出的3种方法中有2种方法对于水深的计算精度更高。因方法的性质不同,不同案例中,各方法表现略有差异,这是正常的。
此外,应当理解,虽然本说明书按照实施方式加以描述,但并非每个实施方式仅包含一个独立的技术方案,说明书的这种叙述方式仅仅是为清楚起见,本领域技术人员应当将说明书作为一个整体,各实施例中的技术方案也可以经适当组合,形成本领域技术人员可以理解的其他实施方式。
Claims (5)
2.根据权利要求1所述的一种水位梯度的加权平均计算方法,其特征在于:所述步骤S4中推求加权系数wL、wR有3种方法,分别为稳定性更好的方法A;精度更高的方法B;所需参数较少的方法C。
3.根据权利要求2所述的一种水位梯度的加权平均计算方法,其特征在于:所述方法A的具体步骤为,
A1、通过网格i-1、i、i+1的流速判断流向,如果网格i、i+1的流速均大于零,则网格i+1为下游、网格i-1为上游;如果网格i、i-1的流速均小于零,则网格i-1为下游、网格i+1为上游;其它情况则不区分上下游;
A2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
A3、如果区分上下游,则先计算网格i的左界面的移动速度,其计算方法为:
4.根据权利要求2所述的一种水位梯度的加权平均计算方法,其特征在于:所述方法B的具体步骤为,
B1、通过网格i-1、i、i+1的流速判断流向,如果网格i、i+1的流速均大于零,则网格i+1为下游、网格i-1为上游;如果网格i、i-1的流速均小于零,则网格i-1为下游、网格i+1为上游;其它情况则不区分上下游;
B2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
B3、如果区分上下游,则先计算网格i的左界面的移动速度,其计算方法为:
5.根据权利要求2所述的一种水位梯度的加权平均计算方法,其特征在于:所述方法C的具体步骤为,
C1、通过网格i的流速判断流向,如果网格i的流速大于零,则网格i的右侧为下游、左侧为上游;如果网格i的流速小于零,则网格i的左侧为下游、右侧为上游;其它情况则不区分上下游;
C2、如果不区分上下游,则wL、wR的值均为0.5,无需后续步骤;
C3、如果区分上下游,则先计算网格i的Fr数,其计算方法为:
C4、然后计算水位差Δzs,如果网格i的右侧为下游,则Δzs=zs,i-zs,i+1,如果网格i的左侧为下游,则Δzs=zs,i-zs,i-1;
C5、先计算网格i的上游侧界面的权重因子wup,如果水位差Δzs>0则wup=min(0.5,Fr2),否则wup=0;
C6、如果网格i的左界面为上游侧,则左界面的权重系数wL=wup,如果左界面在下游侧,则wL=1-wup;
C7、右界面的权重系数wR=1-wL。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383921.6A CN111695304B (zh) | 2020-05-08 | 2020-05-08 | 一种水位梯度的加权平均计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383921.6A CN111695304B (zh) | 2020-05-08 | 2020-05-08 | 一种水位梯度的加权平均计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111695304A CN111695304A (zh) | 2020-09-22 |
CN111695304B true CN111695304B (zh) | 2022-08-02 |
Family
ID=72477406
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010383921.6A Active CN111695304B (zh) | 2020-05-08 | 2020-05-08 | 一种水位梯度的加权平均计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111695304B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111046605A (zh) * | 2019-12-09 | 2020-04-21 | 天津大学 | 基于人工神经网络的河道-泛区洪水演进异维融合模拟原理与计算方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2933743B1 (en) * | 2014-04-18 | 2019-05-08 | Université de Rennes 1 | Method of characterizing molecular diffusion within a body from a set of diffusion-weighted magnetic resonance signals and apparatus for carrying out such a method |
CN108629135B (zh) * | 2018-05-11 | 2021-07-20 | 中国水利水电科学研究院 | 非统一高精度曲面网格水流水质模拟及可视化方法和系统 |
-
2020
- 2020-05-08 CN CN202010383921.6A patent/CN111695304B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111046605A (zh) * | 2019-12-09 | 2020-04-21 | 天津大学 | 基于人工神经网络的河道-泛区洪水演进异维融合模拟原理与计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111695304A (zh) | 2020-09-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109669049B (zh) | 一种基于卷积神经网络的粒子图像测速方法 | |
CN114091163B (zh) | 一种基于有限体积法的河道构筑物过流数值模拟方法 | |
CN104268943B (zh) | 一种基于欧拉‑拉格朗日耦合方法的流体仿真方法 | |
CN107451372B (zh) | 一种运动波与动力波相结合的山区洪水过程数值模拟方法 | |
CN112270115B (zh) | 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法 | |
CN106372320B (zh) | 一种采用亚滤波尺度模型对公路隧道湍流进行大涡模拟的方法 | |
KR101319996B1 (ko) | 유체 시뮬레이션 방법 | |
Thürey et al. | Interactive free surface fluids with the lattice Boltzmann method | |
CN105022928B (zh) | 一种飞行器燃油系统重心位置的数字化实时确定方法 | |
CN105913457A (zh) | 基于加强粒子群优化的粒子滤波目标跟踪方法 | |
CN113779671A (zh) | 一种基于时空步长自适应技术的明渠调水工程水动力实时计算方法 | |
CN107038308B (zh) | 一种基于线性内插的规则格网地形建模方法 | |
KR101106548B1 (ko) | 텐서형 와점성계수를 가진 2차원 하천흐름모형을 이용하여 천수흐름을 해석하는 방법 | |
CN110147646B (zh) | 一种数值模拟框架下线性挡水构筑物的过流处理方法 | |
CN110059399A (zh) | 一种数值模拟跌流水舌的方法 | |
CN116108965A (zh) | 一种城市洪涝预测模型的计算方法及装置 | |
CN111695304B (zh) | 一种水位梯度的加权平均计算方法 | |
CN110175377A (zh) | 一种关于水翼云状空化的数值预报方法 | |
CN106844889A (zh) | 基于非等分区块化原则的四边形结构网格划分方法 | |
CN102592304A (zh) | 一种分类应变限制的动态服装仿真方法 | |
CN108509741A (zh) | 一种泥石流方程的有限元数值求解方法 | |
CN116595745B (zh) | 基于粒子水平集的强冲击物质界面高精度捕捉方法 | |
CN117669408A (zh) | 基于水动力模拟的河网安全评估方法和系统 | |
Dagenais et al. | Real-Time Virtual Pipes Simulation and Modeling for Small-Scale Shallow Water. | |
KR101135307B1 (ko) | 유한체적법에 의한 보존변수 구성 방법 및 이를 이용한 2차원 유체 흐름 해석 방법 |
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 |