CN107436977A - 复杂河网分流的数值模拟方法 - Google Patents
复杂河网分流的数值模拟方法 Download PDFInfo
- Publication number
- CN107436977A CN107436977A CN201710604786.1A CN201710604786A CN107436977A CN 107436977 A CN107436977 A CN 107436977A CN 201710604786 A CN201710604786 A CN 201710604786A CN 107436977 A CN107436977 A CN 107436977A
- Authority
- CN
- China
- Prior art keywords
- mrow
- mfrac
- msub
- mover
- river
- 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.)
- Pending
Links
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)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及复杂河网流场数值模拟的技术领域,公开了复杂河网分流的数值模拟方法,利用二维非恒定浅水方程组,搭建四边形网格与三角形网格混交的网格体系并用有限体积法进行离散,从而建立复杂河网的平面二维水流数学模型,模拟河网中实际的水流运动特征;利用实测资料率定模拟结果,确定模型的相关参数;随后引入分流不均匀系数衡量流量分配特征,计算河网中多个分汊节点的分流不均匀系数。这样,不仅避免了一维模型在空间上的局限性和三维模型计算时间成本过大问题,而且克服了多维耦合模型在衔接中出现质量及动量局部不守恒现象的出现。
Description
技术领域
本发明涉及复杂河网流场数值模拟的技术领域,尤其是涉及具有多汊道的复杂河网流量分配的数值模拟。
背景技术
世界范围内许多大型河流,比如尼罗河、湄公河、密西西比河和珠江等河流进入下游河口区往往都会出现河道分汊现象,从而发育成形态复杂的河网。这些河网由众多的分汊河道和汊口组成,水流在某一支汊的运动,都对河网的整体有“牵一发而动全身”的影响。由于河网多位于河口地区,它不仅受到上游径流的作用,而且还受到外海潮汐的影响,使得河网中的动力作用十分复杂。
分汊河道的流量分配一直是河流动力学研究中的经典问题。它关系到泥沙、营养物及污染物等入海物质通量的分配,从而影响下游河道的生态系统和相关的航道、桥梁工程。而复杂河网中流量分配则涉及到更多的影响因素和更复杂的作用机制,十分值得探讨。
本发明拟针对具有复杂河网的珠江河口,建立整体的平面二维数值模型,模拟在径流和潮汐作用下,河网内的流量分配过程。区别以往的分流比系数,本发明引入分流不均匀系数用来衡量流量分配的不均匀性。并选取河网中多个主要河道的分汊口,计算一系列的流量不均匀系数,定量研究河网分流特征。克服了对复杂河网分流研究的难点,以期对河网中物质通量分配和河口区洪水、咸潮入侵等灾害研究提供参考。
发明内容
本发明提供的复杂河网分流的数值模拟方法,旨在建立一种在径流和潮汐共同作用下河网的平面二维数值模型,引入分流不均匀系数描述河网流量分配的不均匀性,克服现今对复杂河网分流研究的不足。
本发明是这样实现的,利用二维非恒定浅水方程组,搭建四边形网格与三角形网格混交的网格体系并用有限体积法进行离散,从而建立复杂河网的平面二维水流数学模型,模拟河网中实际的水流运动特征;利用实测资料率定模拟结果,确定模型的相关参数;随后引入分流不均匀系数衡量流量分配特征,计算河网中多个分汊节点的分流不均匀系数。
进一步地,所述二维非恒定浅水方程组为:
式中:η为水位;h为静水深;t为时间;x、y、z为右手笛卡尔坐标系;u、v为流速在x、y方向上的分量;Pa为当地气压;ρ为水密度,ρ0为参考水密度;为科氏力;和为地球自转引起的加速度;Sxx、Sxy、Syx、Syy为辐射应力在各方向的分量;τsx、τsy、τbx、τby为水面与河床边界水流切应力;Txx、Txy、Tyx、Tyy为水平黏滞应力;S为源汇项;us、vs为源汇项水流流速。
进一步地,所述四边形网格设置在河道中,所述三角形网格设置在河道相交处及河口外的海滨区。
进一步地,所述分流不均匀系数为ψ,其计算方程式为:
式中:ψ为分流不均匀系数;<Q>为潮周期平均流量,1、2下标分别代表两个支汊河道。
与现有技术相比,本发明提供的复杂河网分流的数值模拟方法,通过建立复杂河网的平面二维水流数学模型,模拟河网中实际的水流运动特征,并利用实测资料率定模拟结果,确定模型的相关参数,随后引入分流不均匀系数衡量流量分配特征,计算河网中多个分汊节点的分流不均匀系数,这样,不仅避免了一维模型在空间上的局限性和三维模型计算时间成本过大问题,而且克服了多维耦合模型在衔接中出现质量及动量局部不守恒现象的出现。
本发明具有的有益效果:
1、充分考虑复杂河网中各个汊道之间的相互影响,精确地计算复杂河网中多汊口的分流情况,从而研究河网中的分流特征。
2、适用于在径流和潮汐共同作用下的三角洲河网地区,引入分流不均匀系数有利于分离潮汐作用对分流的影响。
3、适用范围广,对不同地形和不同动力条件下的河网普遍适用。
附图说明
图1是本发明实施例提供的复杂河网分流的数值模拟方法的平面二维水流数学模型范围和河网区网格的示意图;
图2是本发明实施例提供的复杂河网分流的数值模拟方法的平面二维水流数学验证模型的水位站点位置图及汊口位置示意图;
图3是本发明实施例提供的复杂河网分流的数值模拟方法的部分站点的分潮振幅验证图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
以下结合具体实施例对本发明的实现进行详细的描述。
本实施例的附图中相同或相似的标号对应相同或相似的部件;在本发明的描述中,需要理解的是,若有术语“上”、“下”、“左”、“右”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制,对于本领域的普通技术人员而言,可以根据具体情况理解上述术语的具体含义。
参照图1~3所示,为本发明提供的较佳实施例。
本发明提供的复杂河网分流的数值模拟方法,用于研究复杂河网的流量分配。其中,珠江河口位于中国南部,其三大支流西江、北江和东江在河口区交汇形成复杂河网,因此选取珠江河网为研究区域,模型范围见图1,上边界为高要、石角、博罗、老鸦岗和石咀,下边界到外海-20m左右水深。下面以珠江河网为例,给出本发明提出的复杂河网分流的数值模拟方法的具体步骤和实施方案:
(1)建立二维非恒定浅水方程组,具体如下:
式中:η为水位;h为静水深;t为时间;x、y、z为右手笛卡尔坐标系;u、v为流速在x、y方向上的分量;Pa为当地气压;ρ为水密度,ρ0为参考水密度;为科氏力;和为地球自转引起的加速度;Sxx、Sxy、Syx、Syy为辐射应力在各方向的分量;τsx、τsy、τbx、τby为水面与河床边界水流切应力;Txx、Txy、Tyx、Tyy为水平黏滞应力;S为源汇项;us、vs为源汇项水流流速。
(2)搭建四边形网格和三角形网格混交的网格系统
在河道中采用四边形网格,该网格空间步长根据实际地形取50m-400m;在河道相交处采用三角形网格过渡该网格空间步长介于50-2000m;在河口外海滨区同样采用三角形网格,该网格空间步长介于300-2000m。整体来看,网格的空间步长介于50m-2000m,总数共计56911个。这些网格采用有限体积法进行离散,从而建立复杂河网的平面二维水流数学模型,模拟河网中实际的水流运动特征。
(3)模型参数设置
平面二维水流数学模型的上边界为高要、石角、博罗、老鸭岗和石咀站实测逐时流量,模型下边界为中国近海潮波模型预测的逐时水位。平面二维水流数学模型设定的水流涡黏系数为0.28,河床糙率则根据实际地形率定之后取0.015-0.030。
(4)模型验证
平面二维水流数学模型验证选择珠江河网区的分布的26个水位站点(图2),将站点实测的逐时水位资料利用潮汐调和分析分离出主要的分潮,M2,S2,N2,K2,K1,O1,P1,Q1,分别验证各个分潮的振幅。图3是其中八个水位站点的模型验证结果。图中实测水位和模拟水位大致相等,模型验证良好。
(4)计算分流不均匀系数
在模型验证良好的基础上,选取珠江河网中八个主要的汊口(图2),分别是天河/南华,登州头左/登州头右,黄沙左/黄沙右,左一/右一,左二/右二,麻涌/大盛,百顷/大敖,竹洲/竹银,石龙一/石龙二,计算这八个分汊口的分流不均匀系数定量研究分流,从而研究分析复杂河网的分流特征。
分流不均匀系数的计算公式如下所示:
式中,ψ为分流不均匀系数;<Q>为潮周期平均流量,1、2下标分别代表两个支汊河道。分流不均匀系数ψ的取值范围为[-1,1],如果分流不均匀系数为0,则表明在汊口处的潮周期平均流量均匀分配;如果分流不均匀系数大于0,表明潮周期平均流量通过支汊河道1输运的比支汊河道2的多;特别地,当分流不均匀系数为1时表明潮周期平均流量完全通过支汊河道1输运;反之当分流不均匀系数小于0时,支汊河道1通过的潮周期平均流量则比支汊河道2少。
其中,2001年和2005年枯季平均的汊口分流不均匀系数计算结果如表1所示。
表1分流不均匀系数计算结果
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.复杂河网分流的数值模拟方法,其特征在于,利用二维非恒定浅水方程组,搭建四边形网格与三角形网格混交的网格体系并用有限体积法进行离散,从而建立复杂河网的平面二维水流数学模型,模拟河网中实际的水流运动特征;利用实测资料率定模拟结果,确定模型的相关参数;随后引入分流不均匀系数衡量流量分配特征,计算河网中多个分汊节点的分流不均匀系数。
2.如权利要求1所述的复杂河网分流的数值模拟方法,其特征在于,所述二维非恒定浅水方程组为:
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>v</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>h</mi>
<mi>S</mi>
</mrow>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<msup>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<msup>
<mover>
<mrow>
<mi>&upsi;</mi>
<mi>u</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>f</mi>
<mover>
<mi>&upsi;</mi>
<mo>&OverBar;</mo>
</mover>
<mi>h</mi>
<mo>-</mo>
<mi>g</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&eta;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mi>h</mi>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>gh</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&rho;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>s</mi>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>b</mi>
<mi>x</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>&rho;</mi>
</mfrac>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>x</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>x</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>hu</mi>
<mi>s</mi>
</msub>
<mi>S</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mi>&upsi;</mi>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<mover>
<mrow>
<mi>u</mi>
<mi>&upsi;</mi>
</mrow>
<mo>&OverBar;</mo>
</mover>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>h</mi>
<msup>
<mover>
<mi>&upsi;</mi>
<mo>&OverBar;</mo>
</mover>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<mi>f</mi>
<mover>
<mi>u</mi>
<mo>&OverBar;</mo>
</mover>
<mi>h</mi>
<mo>-</mo>
<mi>g</mi>
<mi>h</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&eta;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mi>h</mi>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>p</mi>
<mi>a</mi>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<mi>gh</mi>
<mn>2</mn>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mrow>
</mfrac>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>&rho;</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>s</mi>
<mi>y</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>&tau;</mi>
<mrow>
<mi>b</mi>
<mi>y</mi>
</mrow>
</msub>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>&rho;</mi>
<mn>0</mn>
</msub>
</mfrac>
<mo>(</mo>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msub>
<mi>S</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>x</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>y</mi>
<mi>x</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mo>&part;</mo>
<mrow>
<mo>&part;</mo>
<mi>y</mi>
</mrow>
</mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>hT</mi>
<mrow>
<mi>y</mi>
<mi>y</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>hv</mi>
<mi>s</mi>
</msub>
<mi>S</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
式中:η为水位;h为静水深;t为时间;x、y、z为右手笛卡尔坐标系;u、v为流速在x、y方向上的分量;Pa为当地气压;ρ为水密度,ρ0为参考水密度;为科氏力;和为地球自转引起的加速度;Sxx、Sxy、Syx、Syy为辐射应力在各方向的分量;τsx、τsy、τbx、τby为水面与河床边界水流切应力;Txx、Txy、Tyx、Tyy为水平黏滞应力;S为源汇项;us、vs为源汇项水流流速。
3.如权利要求1所述的复杂河网分流的数值模拟方法,其特征在于,所述四边形网格设置在河道中,所述三角形网格设置在汊道相交处及河口外的海滨区。
4.如权利要求1所述的复杂河网分流的数值模拟方法,其特征在于,所述分流不均匀系数为ψ,其计算方程式为:
<mrow>
<mi>&psi;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mo><</mo>
<mi>Q</mi>
<msub>
<mo>></mo>
<mn>1</mn>
</msub>
<mo>-</mo>
<mo><</mo>
<mi>Q</mi>
<msub>
<mo>></mo>
<mn>2</mn>
</msub>
</mrow>
<mrow>
<mo><</mo>
<mi>Q</mi>
<msub>
<mo>></mo>
<mn>1</mn>
</msub>
<mo>+</mo>
<mo><</mo>
<mi>Q</mi>
<msub>
<mo>></mo>
<mn>2</mn>
</msub>
</mrow>
</mfrac>
</mrow>
式中:ψ为分流不均匀系数;<Q>为潮周期平均流量,1、2下标分别代表两个支汊河道。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710604786.1A CN107436977A (zh) | 2017-07-24 | 2017-07-24 | 复杂河网分流的数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710604786.1A CN107436977A (zh) | 2017-07-24 | 2017-07-24 | 复杂河网分流的数值模拟方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107436977A true CN107436977A (zh) | 2017-12-05 |
Family
ID=60459799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710604786.1A Pending CN107436977A (zh) | 2017-07-24 | 2017-07-24 | 复杂河网分流的数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107436977A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108090030A (zh) * | 2017-12-22 | 2018-05-29 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种圆形单桩局部流场的处理方法 |
CN108446502A (zh) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
CN109033725A (zh) * | 2018-09-14 | 2018-12-18 | 中国水利水电科学研究院 | 一种定床河工模型试验大面积床面切应力的估算方法 |
CN110442997A (zh) * | 2019-08-14 | 2019-11-12 | 中交(天津)生态环保设计研究院有限公司 | 一种用于指导河口生境改善工程的河口水域范围划定方法 |
CN115470965A (zh) * | 2022-08-08 | 2022-12-13 | 长沙理工大学 | 基于径潮对抗模式的潮汐分汊河口汊道落潮分流比快速确定方法及预测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091065A (zh) * | 2014-07-03 | 2014-10-08 | 南京信息工程大学 | 一种求解浅水问题模拟间断水流数值的方法 |
CN106372375A (zh) * | 2016-10-26 | 2017-02-01 | 武汉大学 | 水库下游两级分汊河道汊河分流比变化的预测方法 |
-
2017
- 2017-07-24 CN CN201710604786.1A patent/CN107436977A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104091065A (zh) * | 2014-07-03 | 2014-10-08 | 南京信息工程大学 | 一种求解浅水问题模拟间断水流数值的方法 |
CN106372375A (zh) * | 2016-10-26 | 2017-02-01 | 武汉大学 | 水库下游两级分汊河道汊河分流比变化的预测方法 |
Non-Patent Citations (3)
Title |
---|
孙卫红等: "基于不均匀系数的水环境容量计算方法探讨", 《水资源保护》 * |
张细兵: "江湖河网水沙运动数值模拟技术研究及应用", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
张蔚等: "长江口分流过程对周期性潮波运动的响应机制", 《第十七届中国海洋(岸)工程学术讨论会论文集》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108090030A (zh) * | 2017-12-22 | 2018-05-29 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种圆形单桩局部流场的处理方法 |
CN108090030B (zh) * | 2017-12-22 | 2021-02-19 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种圆形单桩局部流场的处理方法 |
CN108446502A (zh) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
CN108446502B (zh) * | 2018-03-22 | 2019-01-29 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
CN109033725A (zh) * | 2018-09-14 | 2018-12-18 | 中国水利水电科学研究院 | 一种定床河工模型试验大面积床面切应力的估算方法 |
CN109033725B (zh) * | 2018-09-14 | 2020-05-05 | 中国水利水电科学研究院 | 一种定床河工模型试验大面积床面切应力的估算方法 |
CN110442997A (zh) * | 2019-08-14 | 2019-11-12 | 中交(天津)生态环保设计研究院有限公司 | 一种用于指导河口生境改善工程的河口水域范围划定方法 |
CN115470965A (zh) * | 2022-08-08 | 2022-12-13 | 长沙理工大学 | 基于径潮对抗模式的潮汐分汊河口汊道落潮分流比快速确定方法及预测方法 |
CN115470965B (zh) * | 2022-08-08 | 2023-04-25 | 长沙理工大学 | 基于径潮对抗模式的潮汐分汊河口汊道落潮分流比快速确定方法及预测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107436977A (zh) | 复杂河网分流的数值模拟方法 | |
Huang et al. | A physical similarity model of an impulsive wave generated by Gongjiafang landslide in Three Gorges Reservoir, China | |
Bosa et al. | Shallow water numerical model of the wave generated by the Vajont landslide | |
Lai et al. | Impact of multichannel river network on the plume dynamics in the P earl R iver estuary | |
CN108021780B (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
CN104091065A (zh) | 一种求解浅水问题模拟间断水流数值的方法 | |
CN105512417A (zh) | 基于粒子追踪的孔隙地下水污染物三维运移模拟方法 | |
Tran et al. | Evaluate the influence of groynes system on the hydraulic regime in the ha thanh river, binh dinh Province, vietnam | |
Leake et al. | Assignment of boundary conditions in embedded ground water flow models | |
CN106055774B (zh) | 地面灌溉中地表水流运动的模拟方法 | |
Dey et al. | Hydrodynamics of flow over a gradually varied bed roughness | |
Hasan et al. | Flow simulation for lake Harapan using CCHE2D—a case study | |
Nigussie et al. | Performance evaluation of synthetic unit hydrograph methods in mediterranean climate. A case study at guvenc micro-watershed, Turkey | |
Herrero Casas | Experimental and theoretical analysis of flow and sediment transport in 90-degree fluvial diversions. | |
Pak et al. | Sensitivity analysis for sediment transport in the hydrologic modeling system (HEC-HMS) | |
Khuzaiea et al. | A Hydraulic Model for Identification of Surface Friction Coefficient for Euphrates River within Al Muthanna Governorate, Iraq | |
Halaj et al. | Longitudinal dispersion coefficient impact assessment on HEC-RAS water quality model outputs | |
CN103530462A (zh) | 面向山洪演进数值模拟的计算网格流出率的修正方法 | |
Ashraf et al. | River meandering prediction: Case studies for four rivers in Texas | |
He et al. | Modeling bedload transport trajectories along a Sine-generated channel | |
Jia et al. | Modeling River Morphodynamic Process Using a Depth-Averaged Computational Model and an Application to a Mountain River | |
Scott | Qualitative and Quantitative Flow Visualization Studies on a Distorted Hydraulic Physical Model | |
CN108154006A (zh) | 泥石流流动路径的搜索方法 | |
Šikšnys et al. | Numerical and field investigations of local bridge abutment scour and unsteady downstream river flow from a nearby hydropower plant | |
Đorđević | On the application of the double-averaging methodology in modelling vegetated open-channel flows |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20171205 |