CN118467881B 一种海岸沙丘平面二维风沙输沙率的模拟计算方法 （河海大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN118467881B(45)授权公告日2025.07.04

(21)申请号202410093235.3

(22)申请日2024.01.23

(65)同一申请的已公布的文献号申请公布号CN118467881A

(43)申请公布日2024.08.09

(73)专利权人河海大学

地址211100江苏省南京市江宁区佛城西

路8号

(72)发明人李元李沁倚张弛何岩雨刘建辉郑金海

(74)专利代理机构南京经纬专利商标代理有限公司32200

专利代理师徐尔东

GO6F30/20(2020.01)

GO6F119/14(2020.01)

(56)对比文件

CN109582996A,2019.04.05

邹学勇等.输沙率对土壤颗粒特性和气流湍流脉动的响应.《地球科学进展》.2019,第34卷(第08期),787-800.

审查员毛小丽

(51)Int.CI.

GO6F17/10(2006.01)权利要求书3页说明书8页附图4页

(54)发明名称

一种海岸沙丘平面二维风沙输沙率的模拟

计算方法

(57)摘要

CN118467881B本发明涉及一种海岸沙丘平面二维风沙输沙率的模拟计算方法，读取海岸地形、边界、平面三角形网格以及关键字文件，获取海岸地形数据、泥沙代表粒径、边界处的风速数据；继续获取临界起动风速、推求高程梯度并进行插值计算；计算风沙输沙率；本发明实现在以一定角度，一定高度入射海滩的风在经过滩面不同高程地形

CN118467881B

二维结四风》

二维结四风》

读取地形、边界、网格、关键一

字文件

地面高程

推求高程梯度并插值

如果该点临界启动风速小于该点风速，带入输沙公式计算

血率*0

输出结果文件

如果该点临界启动风速大于该点风速

该点输沙率为计算值

泥沙代表粒径

临界启动风速

各点风速

边界风速

CN118467881B权利要求书1/3页

2

1.一种海岸沙丘平面二维风沙输沙率的模拟计算方法，其特征在于：具体包括以下步

骤：

步骤S1:读取海岸地形、边界、平面三角形网格以及关键字文件，获取海岸地形数据、泥沙代表粒径、边界处的风速数据；

步骤S2:设定x方向代表向离岸方向，且离岸方向为正方向；y方向代表沿岸方向，且下游方向为正方向；z方向为地形高程方向，且竖直向上为正；假设位于边界处的风以预设角度入射海岸，为入射方向与x轴正方向形成的夹角；

其中表征风剪切应力在地形剖面上分布的公式为：

(1)

(2)

公式(1)和公式(2)中，表示x方向上的表面剪切应力扰动值，表示y方向上的表面

剪切应力扰动值，表示地面高程在x方向上的偏导数，表示地面高程在y方向上的偏导数；

步骤S3:采用平面三角形网格计算地形高程梯度与

步骤S4:基于步骤S3获取的地形高程梯度，在拟计算风速的平面三角形网格节点(x,y)

上沿着x轴以及y轴分别做剖面，在剖面上取等间距的插值节点，在插值节点上对分别插值，以获得平面三角形网格任意位置处的地形高程梯度；

步骤S5:将剖面上的任意插值节点记为(x?,y?),选取离(x?,y?)最近的四个网格节点(x,y),并计算其平均值作为插值节点处的函数值；

步骤S6:基于步骤S5获得的插值节点上的函数值，对公式(1)和(2)在x剖面与y剖面上

进行积分算出待求点(x,y)某一高度处的风速分量：

(9)

此处的待求点(x,y)某一高度处指的是与实际海岸上测量风速的风速仪等高的高度；步骤S7:根据步骤S1获得的泥沙代表粒径计算近底临界起动风速：

CN118467881B权利要求书2/3页

3

(10)

公式(10)中，u*为近底临界起动风速，A为经验系数，取0.08,p,[kg/m3]为泥沙颗粒密度，pa[kg/m3]为空气密度，d[m]为泥沙代表粒径；

步骤S8:通过壁面函数将近底临界起动风速转化为一定高度处的临界起动风速：

u.=but(11)

(12)

公式(11)和公式(12)中，z为近底粗糙高度，即理论