CN102496168B - 一种用于河道水文数值模拟的复杂河道网格化方法 - Google Patents
一种用于河道水文数值模拟的复杂河道网格化方法 Download PDFInfo
- Publication number
- CN102496168B CN102496168B CN201110372235.XA CN201110372235A CN102496168B CN 102496168 B CN102496168 B CN 102496168B CN 201110372235 A CN201110372235 A CN 201110372235A CN 102496168 B CN102496168 B CN 102496168B
- Authority
- CN
- China
- Prior art keywords
- grid
- river course
- river
- river channel
- mesh
- 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 44
- 238000004088 simulation Methods 0.000 title claims abstract description 17
- 238000011144 upstream manufacturing Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000013459 approach Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明用于河道水文数值模拟的复杂河道网格化方法,首先将带有支流出入汇、江心洲等各种复杂河道区域自由分割为简单区域集合,根据河道数值模拟需要设置区域控制曲线网格间距,生成二维网格,归并简单区域网格得到复杂河道剖分后的初始网格。对于河道岸线特定区域,可根据需求灵活增加区域网点密度实现该区域内网格加密划分。通过求解拉普拉斯方程,设定收敛精度,本发明能够实现初始网格正交化,最终得到河道二维正交网格。本发明构思巧妙,适应性强,可对含有江心洲的复杂河道进行网格化处理,且可根据需要实际进行指定网格加密。
Description
技术领域
本发明涉及一种用于河道水文数值模拟的复杂河道网格化方法,属于流体力学中的网格生成技术及地理信息***(GIS)中的空间数据内插技术领域。
背景技术
河道数值模拟综合了计算机科学、流体力学、偏微分方程数学理论、计算几何、数值分析等学科,它的发展除依赖于这些学科的发展外,更直接表现在对河道平面模拟的二维网格生成技术的依赖。在应用诸如有限差分法、有限体积法等方法对河道数值模拟进行数值求解时,需要生成相应的计算网格,生成网格质量的优劣直接影响流畅计算的精度和效率。
网格的生成问题可以看成是一个坐标变换的边值问题,建立一种物理域的边界点同计算域的边界点的映射关系,从而求出物理域与计算域内部结点间的对应关系,将不规则的物理域转化成规则的计算域,这种方法可归结为在给定的第一边界条件下求解偏微分方程。这种转换关系如图1所示,左侧为物理域,右侧为计算域。
目前常用的河道网格剖分方法有河道Thompson法、椭圆型网格生成法、代数法等。河道Thompson法主要是通过物理平面(天然河道平面)与变换平面(数模计算平面)之间Poisson方程边值问题数值解实现二维正交网格的生成;椭圆型网格生成法,求解椭圆形偏微分方程具有曲线坐标系的光滑性和可微性,并且所生成的网格在解决复杂外形方面具有生成原理直观、易于编制程序等优点。代数法是指通过一些代数关系式把物理平面上的不规则区域转换成计算平面上矩形区域的方法。这种网格生成方法要找出合适的插值函数,一般采用Langrange插值为基础构造双边界法、无限插值法等。由于无限插值法是把2个方向4个边界的信息同时向内部传递,插值的点数是无限的,用它来生成初始网格,迭代以后得到的是均分网格。这些方法虽然能实现对河道区域的网格剖分,但自身也存在一些问题。
河道Thompson法存在的主要问题有:(1)复杂洲、滩及岸线河道岸线、非恒定流动岸和数值求解引起的动边界等情况下,网格与河道岸线之间的拟合同样会出现实际偏离;(2)对于常见的宽、窄相间的河道平面形态,二维网格不均匀间距可能导致的数值计算精度问题;椭圆型网格生成法虽然保证了网格线与边界的正交性,但不能很好地控制网格的疏密;代数法未对不规则区域网格进行疏密调整。
这些问题导致上述方法在复杂河道网格化过程中不能有效处理存在众多江心洲的多汊河道,并且缺乏根据河道岸线地势变化等相关因素合理控制网格疏密的能力。在河道网格化过程中网格设计应综合考虑地形变化、物理量的变化梯度、动力特性、计算费用、模型的目的等因素。国内外已有的网格自动剖分算法和软件缺乏根据复杂河道地势变化控制网格大小与节点密度的灵活性。
发明内容
本发明要解决技术问题是:克服现有技术的上述不足,提供一种用于河道水文数值模拟的复杂河道网格化方法,其能为复杂河道数值模拟提供分布均匀、数量可根据用户需求灵活设定的正交网格及离散节点。
为了解决以上技术问题,本发明提供的一种用于河道水文数值模拟的复杂河道网格化方法,包括以下步骤:
第一步、数据读取——读取待网格化的二维河道矢量数据;
第二步、河道分割——沿河道方向将河道分割为若干多边形集合,对于含有江心洲的河道区域,对江心洲两侧的河道分别进行分割,所述分割后的多边形沿河道纵向方向分布;
第三步、网格间距计算——构建多边形边界曲线参数方程,计算多边形边界曲线长度,根据网格化所需行列数目计算划分网格单元时边界曲线网格间距;
第四步、内插网点——根据多边形边界曲线参数方程及网格化所需网格间距,将网点内插于多边形边界曲线上,并计算内插网点的二维坐标;
第五步、网格剖分——根据多边形边界曲线内插网点的数量和二维坐标,计算多边形网格化后网格节点的坐标,并生成网格节点,将生成的网格节点用线段依次连接形成多边形网格,完成多边形的网格剖分;
第六步、网格归并——归并各多边形网格,获得河道网格;
第七步、网格正交化——求解拉普拉斯方程,调整河道网格节点位置,完成河道网格正交化,得到河道平面二维正交网格。
进一步的,本发明用于河道水文数值模拟的复杂河道网格化方法的第二步中,对江心洲两侧河道分别进行分割的方法是:在江心洲上下游的弧段上分别选取一至两个端点,在河岸上选取至所述端点最近的点,连接端点和距端点最近的河岸上的点,若该连接线未与其他江心洲相交,则该连线作为划分河道的边界线,否则在与该端点最近的其他江心洲上选取与之最近的点,两者的连接线作为划分河道的边界线。
本发明所述第五步完成后,针对需要进行网格加密的河道局部区域,选定所需加密的网格、设定加密的行列数目,并针对各网格进行加密,网格加密方法参见所述第三步至第五步。
第七步完成之后,根据河道平面二维正交网格各网格节点得二维坐标及其连接顺序,输出河道最终网格化后的节点坐标。
本发明方法首先将带有支流出入汇、江心洲等各种复杂河道区域自由分割为简单区域集合,根据河道数值模拟需要设置区域控制曲线网格间距,生成二维网格,归并简单区域网格得到复杂河道剖分后的初始网格。对于河道岸线特定区域,可根据需求灵活增加区域网点密度实现该区域内网格加密划分。通过求解拉普拉斯方程,设定收敛精度,本发明能够实现初始网格正交化,最终得到河道二维正交网格。
本发明将河道进行分段(分割),并且针对有江心洲的复杂区域进行巧妙分割为互相独立的多边形,使分割的多边形内均不包含有江心洲,因此本发明方法可处理带有支流出入汇、江心洲等各种复杂河道区域,相对于传统河道网格化方法本发明具有更强的适应性和通用性;本发明方法可针对所需加密区域进行网格加密划分,使河道网格化与复杂河道地势变化相适应,为水文相关运算提供良好的数据基础。
附图说明
下面结合附图对本发明作进一步的说明。
图1为坐标变化示意图。
图2为本发明方法流程图。
图3为待网格化的河道示意图。
图4为划分图3所示河道形成简单多边形集合示意图。
图5为图4的局部放大图。
图6为图3所示河道网格化后生成的网格示意图。
具体实施方式
下面根据附图详细说明本发明,本发明的目的和效果将变得更加明显。
如图2所示,为本发明一种用于河道水文数值模拟的复杂河道网格化方法实施例的流程图,包括以下步骤:
第一步、数据读取——读取待网格化的二维河道矢量数据;如图3所示,河道由非自交弧段组成。
第二步、河道分割——沿河道方向将河道分割为若干多边形集合,对于含有江心洲的河道区域,对江心洲两侧的河道分别进行分割,所述分割后的多边形沿河道纵向方向分布;
本步骤中,对江心洲两侧河道分别进行分割的方法是,在江心洲上下游的弧段上分别选取一至两个端点,在河岸上选取至所述端点最近的点,连接端点和距端点最近的河岸上的点,若该连接线未与其他江心洲相交,则该连线作为划分河道的边界线,否则在与该端点最近的其他江心洲上选取与之最近的点,并进行连接作为划分河道的边界线。
如图5所示,江心洲的上下游弧段分别选取两个端点,在河岸上选取了与其最近的点,并进行了连接,将河道进行了分割。图4中,最右侧的三角洲的上游(左侧)弧段上则选取了一个端点与河岸连接。可根据根据实际情况确定在三角洲上下游弧段选取端点的个数,一般情况选取2个端点的情况较为通用。
计算机具体执行过程如下:输入网格控制弧段,将河道区域划分成简单的连通空间区域。计算弧段与河道岸线交点,利用该点打断河道岸线,并将其加入到边界弧段中。对于每个简单区域,应用GIS拓扑重建功能,将简单区域边界弧段首尾依次相连,构建简单区域多边形。最终将图3所示的复杂河道组织为图4所示的多边形集合。
第三步、网格间距计算——构建多边形边界曲线参数方程,计算多边形边界曲线长度,根据网格化所需行列数目计算划分网格单元时边界曲线网格间距;
根据需求设定多边形网格化所需行列数目即位于边界弧段上的网格节点数目,计算划分多边形平面时横向曲线与纵向曲线网格间距。设定网格行列数目必须保证每个多边形对应横向曲线网格列数相同、对应纵向曲线网格行数相等。对于任意两个拥有相同横向曲线或纵向曲线的多边形也要保证合并后仍满足此条件。
图5为复杂河道划分为简单区域后的局部放大图,本发明以图5中多边形为例描述设定边界网点数目时应满足的条件。图中A、B、C、D表示简单区域多边形,Si(i=1,2,3,...)为多边形边界弧段,Ni表示Si中节点数目。设定多边形A边界网点数目时应满足条件:N2=N3,N1=N4+N5+N6。按此条件设定多边形B、C、D边界网点数目,此外还需满足条件N1=N11+N12+N13,N2+N10=N3+N7,以保证多边形A、B、C、D组成的连通区域对应同向控制曲线网点数目相同。
第四步、内插网点——根据多边形边界曲线参数方程及网格化所需网格间距,将网点内插于多边形边界曲线上,并计算内插网点的二维坐标;
第五步、网格剖分——根据多边形边界曲线内插网点的数量和二维坐标,计算多边形网格化后网格节点的坐标,并生成网格节点,将生成的网格节点用线段依次连接形成多边形网格,完成多边形的网格剖分;
第六步、网格归并——归并各多边形网格,获得河道网格;
第七步、网格正交化——求解拉普拉斯方程,调整河道网格节点位置,完成河道网格正交化,得到河道平面二维正交网格。
第七步完成之后,根据河道平面二维正交网格各网格节点得二维坐标及其连接顺序,输出河道最终网格化后的节点坐标。
针对需要进行网格加密的河道局部区域,在第五步完成后,选定所需加密的网格、设定加密的行列数目,并针对各网格进行加密,网格加密方法参见所述第三步至第五步。
除上述实施例外,本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案,均落在本发明要求的保护范围。
Claims (3)
1.一种用于河道水文数值模拟的复杂河道网格化方法,包括以下步骤:
第一步、数据读取——读取待网格化的二维河道矢量数据;
第二步、河道分割——沿河道方向将河道分割为若干多边形集合,对于含有江心洲的河道区域,对江心洲两侧的河道分别进行分割,所述分割后的多边形沿河道纵向方向分布;
第三步、网格间距计算——构建多边形边界曲线参数方程,计算多边形边界曲线长度,根据网格化所需行列数目计算划分网格单元时边界曲线网格间距;
第四步、内插网点——根据多边形边界曲线参数方程及网格化所需网格间距,将网点内插于多边形边界曲线上,并计算内插网点的二维坐标;
第五步、网格剖分——根据多边形边界曲线内插网点的数量和二维坐标,计算多边形网格化后网格节点的坐标,并生成网格节点,将生成的网格节点用线段依次连接形成多边形网格,完成多边形的网格剖分;
第六步、网格归并——归并各多边形网格,获得河道网格;
第七步、网格正交化——求解拉普拉斯方程,调整河道网格节点位置,完成河道网格正交化,得到河道平面二维正交网格;
所述第二步中,对江心洲两侧河道分别进行分割的方法是,在江心洲上下游的弧段上分别选取一至两个端点,在河岸上选取至所述端点最近的点,连接端点和距端点最近的河岸上的点得到连接线,若该连接线未与其他江心洲相交,则该连接线作为划分河道的边界线,否则在与该端点最近的其他江心洲上选取与之最近的点,并进行连接作为划分河道的边界线。
2.根据权利要求1所述的用于河道水文数值模拟的复杂河道网格化方法,其特征在于:所述第五步完成后,针对需要进行网格加密的河道局部区域,选定所需加密的网格、设定加密的行列数目,并针对各网格进行加密,网格加密方法参见所述第三步至第五步。
3.根据权利要求1或2所述的用于河道水文数值模拟的复杂河道网格化方法,其特征在于:第七步完成之后,根据河道平面二维正交网格各网格节点得二维坐标及其连接顺序,输出河道最终网格化后的节点坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110372235.XA CN102496168B (zh) | 2011-11-22 | 2011-11-22 | 一种用于河道水文数值模拟的复杂河道网格化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110372235.XA CN102496168B (zh) | 2011-11-22 | 2011-11-22 | 一种用于河道水文数值模拟的复杂河道网格化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102496168A CN102496168A (zh) | 2012-06-13 |
CN102496168B true CN102496168B (zh) | 2014-04-16 |
Family
ID=46187990
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110372235.XA Expired - Fee Related CN102496168B (zh) | 2011-11-22 | 2011-11-22 | 一种用于河道水文数值模拟的复杂河道网格化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102496168B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427330A (zh) * | 2015-12-21 | 2016-03-23 | 中科宇图科技股份有限公司 | 一种矢量多边形的任意多比例分割方法 |
CN108021780B (zh) * | 2018-01-24 | 2021-08-13 | 国电南瑞科技股份有限公司 | 一种基于无规则非结构网格模型的山洪动态仿真方法 |
CN108986222B (zh) * | 2018-08-15 | 2022-10-14 | 四川大学 | 无汊河道数字地形生成方法 |
CN109101732B (zh) * | 2018-08-15 | 2021-05-04 | 四川大学 | 基于地形特征界线的无汊河道二维结构网格剖分方法 |
CN110135080B (zh) * | 2019-05-17 | 2020-04-21 | 中科三清科技有限公司 | 水污染的模拟方法及装置 |
CN110189618A (zh) * | 2019-05-28 | 2019-08-30 | 南京大学 | 一种顾及密度差异的河渠线状水系要素自动制图综合方法 |
CN111311734A (zh) * | 2020-02-14 | 2020-06-19 | 中核第四研究设计工程有限公司 | 铀尾矿库排洪***水力特性三维数值模拟方法 |
CN111444561A (zh) * | 2020-03-04 | 2020-07-24 | 中电建生态环境集团有限公司 | 一种河道地形模型建立方法及装置 |
CN115222913B (zh) * | 2022-03-29 | 2024-03-15 | 广州汽车集团股份有限公司 | 网格生成、结构仿真分析方法、装置、设备及存储介质 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009544101A (ja) * | 2006-07-17 | 2009-12-10 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | 医用画像セグメンテーション用の多角形メッシュを用いた効率的ユーザインタラクション |
CN102129715B (zh) * | 2011-03-24 | 2015-05-20 | 山东大学 | 具有任意内部特征约束的几何模型的四边形网格生成方法 |
-
2011
- 2011-11-22 CN CN201110372235.XA patent/CN102496168B/zh not_active Expired - Fee Related
Non-Patent Citations (4)
Title |
---|
TK-2DC软件前处理模块-贴体正交网格生成软件开发;张明进等;《水道港口》;20080415;第29卷(第2期);第3节、图2-12 * |
二维网格快速生成的一种方法;翁培奋;《上海大学学报(自然科学版)》;19961030;第2卷(第5期);摘要、第1节、图3 * |
张明进等.TK-2DC软件前处理模块-贴体正交网格生成软件开发.《水道港口》.2008,第29卷(第2期), |
翁培奋.二维网格快速生成的一种方法.《上海大学学报(自然科学版)》.1996,第2卷(第5期), |
Also Published As
Publication number | Publication date |
---|---|
CN102496168A (zh) | 2012-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102496168B (zh) | 一种用于河道水文数值模拟的复杂河道网格化方法 | |
CN104460343B (zh) | 一种基于水位监测数据的河道洪水淹没模拟方法 | |
CN106682262A (zh) | 一种获取飞行器流场的数值模拟方法 | |
CN105118091B (zh) | 一种构建多精度非均匀地质网格曲面模型的方法和*** | |
CN101287889A (zh) | 利用代数串联级线性解算器改进储层模拟的装置、方法和*** | |
CN102521863A (zh) | 一种基于粒子***的三维流体标矢量统一动态表示方法 | |
CN108986222A (zh) | 无汊河道数字地形生成方法 | |
CN107038308B (zh) | 一种基于线性内插的规则格网地形建模方法 | |
CN107657075B (zh) | 模拟地下水介质交界面处达西速度的区域分解有限元法 | |
CN106844963B (zh) | 模拟开挖至运行全过程的拱坝三维网格模型自动剖分方法 | |
CN104635262A (zh) | 一种基于增强型矩形网格的正逆断层等值线自动生成方法 | |
CN108399290A (zh) | 基于裂隙网络连通性的裂隙流的modflow模拟方法 | |
Evans et al. | Table cartogram | |
Song et al. | Effective simulation of flow in a moderately curved bend with a single short branch to support the design optimization of river-branch–plant configurations | |
CN101105864B (zh) | 一种基于序列剖面拓扑推理的三维地质体动态重构方法 | |
CN107871048B (zh) | 大尺度流域网格分块方法 | |
Chang | Study of a solitary wave interacting with a surface piercing square cylinder using a three-dimensional fully nonlinear model with grid-refinement technique on surface layers | |
Michálková et al. | Imposing angle boundary conditions on B-spline/NURBS surfaces | |
Yan et al. | Analysis free surface of nonlinear seepage using the MQRBF method | |
CN107871041B (zh) | 流域复杂地形的网格处理方法 | |
Holzbecher et al. | The doublet system simulator | |
KR20090069604A (ko) | Cdg 기법을 적용한 유체 흐름의 2차원 유한요소수치해석 방법 | |
Sun | An automated spatial flow layout algorithm using triangulation, approximate Steiner tree, and path smoothing | |
Rajani et al. | Numerical Simulation of Flow past a Circular Cylinder with Varying Tunnel Height to Cylinder Diameter at Re 40 | |
Ko | A study on subdivision scheme-draft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140416 Termination date: 20171122 |
|
CF01 | Termination of patent right due to non-payment of annual fee |