二维潮流泥沙数值模拟及可视化研究

合集下载

泥沙局部冲淤二维数值模拟仿真

泥沙局部冲淤二维数值模拟仿真
( C o l l e g e o f E n g i n e e r i n g , O c e a n U n i v e r s i t y o f C h i n a , Q i n g d a o 2 6 6 1 0 0, C h i n a )
A b s t r a c t : I n o r d e r t o s t u d y t h e l o c a l e r o s i o n a n d d e p o s i t i o n i n a i r v e r o r a n o c e a n , a t w o d i m e n s i o n a l ( 2 一 D)n u m e r -
人结果基本符合 , 模 型的通用性较好. 该模型可 以应用 于其他形式 的局部 冲淤数值模拟 , 具有广阔的应用前景. 关键词 : 局部 冲淤 ; F L U E N T; 自定义 函数 ; 动网格 ;
d o i : 1 0 . 3 9 6 9 / j . i s s n . 1 0 0 6 - 7 0 0 3 . 2 0 1 2 0 3 0 5 5
i e l a mo d e l o n l o c a l e r o s i o n a nd d e p o s i t i o n i s e s t a b l i s h e d t h r o u g h t he r e de v e l o p me n t i n FL UENT c o d e .The lo f w f i e l d i s c a l c u l a t e d b y a p pl y i n g s t a n d a r d e q ua t i o ns .By u s i n g a u s e r d e ine f d f u n c t i o n,t h e s e d i me n t t r a n s p o r t mo d e l i s e m- b e d d e d i n t o FL UENT a n d b e s i d e s.t h e b e d p r o il f e i S s i mu l a t e d b y u t i l i z i n g a d y n a mi c me s h.Th e s t u d y u s e d t wo e x . p e ime r n t c a s e s t o v e if r y t h e r e l i a b i l i t y a n d u n i v e r s li a t y o f t h e s c o u in r g mo de 1 .Th e r e s e a r c h r e s u l t s s ho w t h a t t h e p r e d i c t i o n s f r o m t h e mo d e l a r e i n a g r e e me n t wi t h t he p io r r r e s ul t s.a nd t h e s c o u in r g mo d e l i S S O u n i v e r s a 1 t h a t i t ma y be u s e d t o s i mu l a t e l o c a l s c o ur i n d i f f e r e n t c a s e s .Th e r e f o r e.i t h a s a br o a d a p pl i c a t i o n p r o s p e c t . Ke y wo r d s: l o c a l e r o s i o n a nd de p o s i t i o n;FL UENT;u s e r — de in f e d f u n c t i o n;d y n a mi c me s h

二维潮流泥沙数学模型理论及工程应用的开题报告

二维潮流泥沙数学模型理论及工程应用的开题报告

二维潮流泥沙数学模型理论及工程应用的开题报告一、研究背景和意义潮流泥沙是河流流域中非常重要的物质组成部分,对于河流交通、水文环境工程建设、土地利用等方面的影响较大。

因此,建立潮流泥沙数学模型并研究其理论和工程应用,对于改善河流流域的环境、保护河流生态系统、开展水利工程建设等方面具有重要的现实意义。

二、研究目的和内容本文主要以二维潮流泥沙数学模型理论与工程应用为研究对象,旨在通过对潮流泥沙数学模型理论的深入剖析、建模、验证和应用实例剖析等方面的研究,进一步完善潮流泥沙数学模型在实际工程中的应用,提高潮流泥沙数学模型的预测精度,为实际工程应用提供数据和理论支持。

三、主要研究内容(1)二维潮流泥沙数学模型理论分析及建模;(2)二维潮流泥沙数学模型参数获取方法研究;(3)二维潮流泥沙数学模型验证方法和实验设计;(4)二维潮流泥沙数学模型在实际工程中的应用实例分析;(5)针对现有模型存在的问题对二维潮流泥沙数学模型优化的研究。

四、研究方法和流程本文主要采用数学模型与实验相结合的研究方法,通过对二维潮流泥沙数学模型理论的建模以及对模型进行实验验证,探究模型存在的不足之处并利用相关算法进行优化,力求探究二维潮流泥沙数学模型的科学性和可操作性,并为实际工程应用提供理论与相关数据。

五、预期研究成果(1)建立二维潮流泥沙数学模型并进行相关分析与验证;(2)提出针对模型存在的问题的解决方案;(3)应用二维潮流泥沙数学模型成功开展相关水文环境工程建设;(4)在环境保护、生态建设和水利工程等方面提供相关数据和理论支持。

六、论文框架本文计划包括绪论、潮流泥沙数学模型理论分析、模型应用实例分析、模型优化研究等章节,具体结构如下:1.绪论1.1 研究背景与意义1.2 国内外研究现状1.3 研究目的与意义1.4 研究方法与流程1.5 预期研究成果2.潮流泥沙数学模型理论分析2.1 潮流泥沙数学模型建立原理2.2 潮流泥沙数学模型的基本参数2.3 潮流泥沙数学模型的数学表达式分析3.模型应用实例分析3.1 模型应用实例介绍3.2 模型应用实例分析3.3 模型实例验证及评估结果4.模型优化研究4.1 模型存在的问题及原因分析4.2 模型优化解决方案4.3 优化后的模型应用实例研究5.结论与展望5.1 研究成果总结5.2 研究意义和价值5.3 研究展望七、研究进度安排1. 绪论部分2. 阅读相关文献并进行文献综述3. 完成潮流泥沙数学模型的建立和原理分析4. 采集相关数据并进行模型的参数获取5. 完成模型应用实例研究6. 完成模型优化研究7. 确认论文框架并完成论文初稿8. 完成论文修改9. 论文定稿及论文答辩10. 准备论文终稿提交和答辩。

平面二维水流泥沙数值模拟

平面二维水流泥沙数值模拟

平面二维水流泥沙数值模拟河流泥沙动力学研究。

666水科学进展第13卷22y +U+V=-g- + t+ x y(3)悬移质泥沙扩散方程+U+V= +-(S-S)s*H x yZb = (S-S*)22(4)(5)河床变形方程式中U、V分别为垂线平均流速在x、y方向上的分量;Zs、Zb和H分别为水位、河底高程和水深;g为重力加速度; 为水的密度; y方向上的分量:( x、y分别为底部切应力在x、x; y)=U+V(U;V);C为谢才系数; t为水流紊动粘性系数,由零方程紊流模型确定;S和S*CH 分别为垂线平均含沙量和挟沙力; s为泥沙紊动扩散系数; 为泥沙沉速; 为床沙干容重; 为悬移质泥沙恢复饱和系数。

水流挟沙力和分组挟沙力级配采用李义天方法进行计算。

数值计算定解条件为:进口给定流量、含沙量及级配,出口给定水位。

对固壁边界,取法向流速为0。

初始时刻给出地形、水位、流速和含沙量等物理量的初始值。

2 有限元方程为进行有限元计算,首先需对计算区域进行网格划分。

模型基于波前推进法思想,建立了河道三角形网格自动剖分系统。

即通过编制通用的程序,只需输入少量信息(边界线及地形),便可自动生成充满整个计算域的三角形网格,并配以自动绘图及屏幕显示。

该系统能自动获得网格的节点坐标、高程及单元关联信息,从而解决了有限元计算前处理工作量大的困难。

整个计算域共划分为NE个单元,NP个节点,节点整体编号为i,i=1,2, ,NP。

对水沙方程(1)~(4),用Calerkin加权余量法进行积分,将单元形函数经整理得到如下有限元方程:AijAijdZsjdt=-D1ij(HU)j-D2ij(HV)jdUj xj=-BijUj-gD1ijZsj-Aij- t(Cij+Eij)Ujdt(6)(7)(8)(9)代入积分方程,dVj yjAijdt=-BijVj-gD2ijZsj-Aij- t(Cij+Eij)VjAij其中Aij=B2ijkdSj=-BijSj- Aij[(S-S*)/H]js(Cij+Eij)Sj- dt = dxdy , jkij idxdy, Bij=B1ijkUk+B2ijkVk, B1ijk=Cij= k i jdxdy, i j i jdxdy,+第6期张细兵、殷瑞兰:平面二维水流泥沙数值模拟667D1ij= i jdxdy, D2ij= i jdxdy, Eij= j j idx- idy 为形函数;i、j、k为节点整体编号。

二维非均匀水沙问题的数值模拟的开题报告

二维非均匀水沙问题的数值模拟的开题报告

二维非均匀水沙问题的数值模拟的开题报告1. 题目二维非均匀水沙问题的数值模拟2. 研究目的水沙运动是地表水环境研究中的重要内容,对于水资源的管理和保护具有重要意义。

本研究旨在通过数值模拟的方法,研究二维非均匀水沙问题的运动规律和演化过程,为河流污染的预测和治理提供科学依据。

3. 研究内容(1)建立数学模型,包括流体力学方程、输沙方程和边界条件等;(2)采用有限差分或有限元等数值方法,对模型进行数值求解;(3)模拟并分析不同水文条件下,河流污染物携带和输移规律;(4)研究不同底质条件下,河流泥沙的运动规律和治理措施。

4. 研究方法(1)建立非均匀水沙数值模型,包括流体力学方程、输沙方程等;(2)采用有限差分或有限元等数值方法,对模型进行数值求解;(3)运用可视化技术,对数据进行可视化分析,并与实测数据进行对比;(4)对模型进行敏感性分析,确定模型的可靠性和适用范围。

5. 研究意义(1)为河流水污染的预测和治理提供科学依据;(2)为河流工程设计提供重要参考;(3)完善水环境数值模拟方法,提高对水环境的管理和保护能力。

6. 参考文献[1] Liu Y, Torfs P J J F, Uijttewaal W S J. A new two-dimensional non-hydrostatic model of waves, currents and sediment transport in coastal waters [J]. Coastal Engineering, 2007, 54(1): 85-103.[2] Wang X, Leng G, Tang C, et al. A two-dimensional hydro-morphological model for simulating sediment spills in shallow rivers [J]. Environmental Modelling & Software, 2017, 96: 322-341.[3] Yang S, Zhang X, Shi G. Numerical simulation of water and sediment transport in upper reaches of Yellow River based on 2D mathematical model [J]. Journal of China University of Petroleum, 2012, 36(1): 133-140.。

二维垂向淹没射流冲刷粗砂砂床的数值模拟

二维垂向淹没射流冲刷粗砂砂床的数值模拟

吉 ( R 。 o y ) + G + - b y 一 R S O R (
O W +
) ( 3 )
警 ) =
u a O W+
r A
w R O


a y

向: 疏 浚设备 与技术 。E - m a i l : j i a n j u n w a n g h h u c @1 2 6 . t o m。
中图法分类号
U 4 4 2 . 3 2 ;
文献标 志码
A 泥船 在工 作 时是航 行 的 , 向 压 水 射 流 在 栗 处 冲 刷 的
耙 吸式挖 泥船 , 是水 力式 挖泥 船 中大型 自航 、 自
载式挖泥船 , 因其具有广泛的适应性和经济性 , 是目 前世界上使用较广泛 的一种挖泥船 , 它适用于开挖 淤泥、 黏土 、 砂壤土及各种砂土。现有耙 吸挖泥船在 挖 掘 泥土时 , 一般 配备 高压 水流 破土 , 工作 时高 压水 流从喷嘴射出 , 将板结砂土冲碎松散 , 再通过耙齿 的 切削 , 冲切相互作用使泥砂容 易悬浮 , 便 于吸入泥 管, 提 高挖 泥效 率 。高 压 水 射 流 冲刷 技 术 在耙 吸式 挖泥 船上 的应 用 有 效 地 提 高 了耙 吸 式 挖 泥 船 的产 量, 对 其进 行研 究是很 有 必要 的 , 本 文研究 二维 垂 向 淹没 射流 对无 黏性粗 砂砂 床 的冲刷 效果 。 对 于射 流 局 部 冲刷 问 题 , 目前 国 内外 学 者 已开 展了广泛的研究。李文学根据实验结果得出射流冲 刷坑 的几何相似性并建立了射流冲刷坑特征长度的 经验公示 I 2 ; 槐文信在计算和实验研究 的基 础上 给 出 了射 流 冲刷平 衡 时冲刷 坑特 征长度 的半 经验 表 达 式 J 。近 年 来 , 随着 计算 流体力 学 ( C F D) 的 发

二维潮流数学模型可视化平台建设

二维潮流数学模型可视化平台建设
( 天 津大学 河流海 岸工 程泥 沙研究所 , 天津 3 0 0 0 7 2 ) 摘 要 : 采用 V B与 F O R T R A N 混 合 编 程 技 术 为 天 津 大 学 泥 沙 实 验 室 自主 开发 的 I S T U 数 学 模 型 制 作 可
视 化操 作 平 台, 并 对 AU T O CA D 制 图软 件 进 行 二 次 开 发 实 现 二 维 潮 流 流 场 图 的绘 制 。 本 平 台 集 成 度 高 、 界 面友 好 、 操 作 简便 、 计 算 快 速 可 广 泛 适 用 于二 维 潮 流 的数 值 模 拟 。
关键词 : 混 合 编 程 ;二 次 开 发 ;可 视 化
中图分 类号 : G 4 3 2 ; T P 6 3 5
文献 标志 码 : A
文章编 号 : 1 0 0 6— 7 1 6 7 ( 2 0 l 3 ) 0 3— 0 0 9 2— 0 3
Co n s t r u c t i o n o f a T wo . d i me n s i o n a I Mo d e l f o r Ti d a I Vi s u a l i z a t i o n P I a t f O r m
的数值 模拟 软件 , 具有 潮流 、 波浪 、 温度 、 泥沙 等独立 模
块, 已成功 应用 于数 项大 型工 程 , 采用 F O R T R A N语 言
编 写 。
制数 学模 型操 作界 面 , 实 现参数 输入 数组 赋值等 功能 ,
采用 F O R T R A N进 行 循 环 计算 , 则 可 以 同 时 利 朋 两 种
BAI 一 c hu a n. Y I NG Ch a o
( I n s t i t u t e o f S e d i me n t o n Ri v e r a n d C o a s t E n g i n e e r i n g ,T i a n j i n Un i v e r s i t y ,T i a n j i n 3 0 0 0 7 2,Ch i n a )

盘锦新港区二维潮流泥沙数学模型研究

盘锦新港区二维潮流泥沙数学模型研究

b i r e l f y i n t r o d u c e s he t n a t u r l a c o n it d i o n nd a l f u v i a l p r o c e s s o f t h e s e a b e d ,a n d a l s o r e s e a r c h e s t h e h y d I o d y I l c s s u r r o u n d i n g i n i t i a l p r o j e c t nd a p l a n n i  ̄p r o j e c t .T h e r e s l u t s s h o w ha t t h y a  ̄ o d y n a ie m i f e l d ft o h e s u r r o u n i d n g a r e a fp o r o j e c t c h ng a e s  ̄ i g h l f y ,a n d t h e
GAO Xi a n g — y u, D OU Xi - p i n g, P AN Yu n
( K e y L a b o f P o r t ,Wa t e r w a y a n d S e d i m e n t a t i o n E n g i n e e r i n g f o M i n i s t r y o f
高 祥 宇 ,窦 希萍 ,潘 昀
( 南京 水利科 学研 究院港 口航 道 泥 沙工程 交通部 重 点 实验 室 ,江苏 南京
2 1 0 0 2 9பைடு நூலகம்)

要 :针对盘锦港新 港区规划建设方案 ,在分析该 海域 自然 条件和岸滩 演变 的基础 上 ,建立 了工程 海域

波—流共同作用下长江口二维悬沙数值模拟

波—流共同作用下长江口二维悬沙数值模拟

波—流共同作用下长江口二维悬沙数值模拟波—流共同作用下长江口二维悬沙数值模拟摘要:本文通过数值模拟方法,研究了波—流共同作用下长江口内二维悬沙运移的情况。

通过建立数学模型和计算方法,模拟了长江口内悬沙的输运过程,并对波浪和流场的相互影响进行了探讨。

结果表明,波浪和流场的共同作用对悬沙输运过程具有显著影响,进一步深化了对长江口悬沙运动规律的认识。

一、引言长江口是我国最重要的河口之一,其水动力和悬沙运动规律的研究对于河口水文、水沙动力学的理解具有重要的意义。

波浪和流场作为长江口内两个主要的力学驱动因素,其相互作用对悬沙的输运过程有着重要影响。

因此,基于数值模拟方法,研究波流共同作用下长江口内二维悬沙运移的规律具有重要的理论和实际意义。

二、模型建立在本文研究中,建立了二维悬沙运移数值模型,考虑了波浪、流场以及悬沙之间的相互作用。

模型采用了数值解的方法,通过离散化长江口内区域,并运用流动方程、输运方程和释放方程来模拟波流共同作用下的悬沙运动过程。

模型的输入数据包括入口流速、入口浓度和入口波浪参数等,模拟的输出结果为悬沙在长江口内的输运情况和浓度分布。

三、数值模拟结果通过数值模拟得出的结果显示,长江口内的悬沙输运过程受到波浪和流场的共同作用。

在波浪作用下,悬沙会随着波动而上升和下降,这会影响悬沙的输运速度和方向。

而流场则会决定悬沙的平均输运速度和方向。

模拟结果还显示,在大波浪和强流场共同作用的情况下,悬沙会有较大的浓度差异和输运速度。

四、悬沙输运规律分析根据数值模拟结果,我们进一步分析了悬沙输运的规律。

首先,当波浪和流场方向相一致时,悬沙的输运速度较快;当方向相反时,输运速度较慢。

其次,波浪的幅度和周期对悬沙的上升和下降速度有较大影响,较大的波浪幅度和周期会引起更剧烈的悬沙运动。

此外,在大流量的情况下,悬沙的输运速度较快,浓度较高。

五、结论通过本文的研究,我们深入了解了长江口内波流共同作用下的悬沙运移规律。

射阳河口二维潮流、泥沙数模研究的开题报告

射阳河口二维潮流、泥沙数模研究的开题报告

射阳河口二维潮流、泥沙数模研究的开题报告
一、研究背景
射阳河口是长江下游区域中的一个典型河口,是江苏省南通市的一
条重要的出海航线,也是运输航道,是南通港口的重要组成部分。

射阳
河口水流、泥沙输移特征是研究该河口水文及海洋动力过程的重要方面。

因此,对射阳河口的潮流和泥沙数值模拟具有重要的科学意义和实际意义。

二、研究目的
通过对射阳河口的潮流和泥沙数值模拟研究,对河口水文和海洋动
力特征进行深入分析和探讨,对河口周边的生态和经济环境提供科学的
参考和决策支持。

三、研究内容和方法
1.研究内容:
(1)射阳河口潮流特征的数值模拟
(2)射阳河口泥沙输移特征的数值模拟
(3)对河口水动力和泥沙输移的影响因素进行分析
2.研究方法:
(1)收集射阳河口相关实测数据
(2)建立射阳河口二维潮流、泥沙数值模拟模型
(3)对模型进行校验和验证
(4)模拟并分析河口水动力和泥沙输移的影响因素
四、研究意义
1. 对射阳河口的水动力和泥沙输移特征进行深入分析和研究,为该区的生态环境和经济建设提供科学支持和决策参考。

2. 丰富河口水文和海洋动力研究领域,扩大模拟模型的应用范围和科学价值。

3. 对于其他地区的类似研究具有借鉴和参考意义。

基于GIS泥石流二维数值模拟_王纯祥

基于GIS泥石流二维数值模拟_王纯祥

− ⎟ ⎟ ⎠
⎛ M in, j +1/ 2 + M in−1, j +1/ 2 ⎜ 1 ⎜ 2hin−+ 1/ 2, j −1/ 2 ⎝
networks)。其中最常用的是栅格网。在这里采用有
限差分法求解泥石流的连续方程和运动方程,栅格 网作为有限差分法的网格。 在数字高程模型中,每一栅格有 8 个可能的流 用角度表示分别为 0°(东), 动方向, 如图 1(a)所示,
1 前 言
在泥石流的灾害和风险性评价中,关键是预测 其沿着复杂三维地形的传播运动轨迹、流动距离和 泛滥范围,并预测可能受到灾害的范围、房屋和路 段。数值模拟方法和 GIS 的集成是一重要的分析工 具 。 本文首先给出了一个泥石流的二维数值模拟公 式。该方程是在假定泥石和雨水组成的固液混合物 遵循均匀、连续、不可压缩的非定常的牛顿流体的 运动规律, 基于质量守恒方程和 Naiver-stokes 方程, 采用深度积分方法推导出的[2]。由于 GIS 中的栅格
GIS-based two-dimensional numerical simulation of debris flow
WANG Chun-xiang1, BAI Shi-wei2, ESAKI Tetsuro1, MITANI Yasuhiro1
(1. Institute of Environment Systems, Kyushu University, Fukuoka, 812-8581, Japan; 2. Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China)
Abstract: To predict the run out distance and extent of the hazard area of rainfall-induced debris flows, we proposed a depth-averaged two-dimensional numerical model, in which the debris and water mixture is assumed to be a uniform continuous, incompressible, unsteady Newtonian turbulence fluid. As raster grid networks of digital elevation model in GIS can be used as the finite difference mesh, the continuity and momentum equations are solved numerically using the finite difference method. All the input data and output are processed in GIS. The model is applied to simulate the rainfall-induced real debris flow of July 20, 2003, in Minamata City of southern Kyushu, Japan. The simulation reproduces the propagation of the real debris flow; and the results are in good agreement with the field investigation. The two-dimensional numerical simulation using GIS can be used to estimate the potentially inundated areas, the potentially dangerous homes and roads for hazard and risk analysis. Key words: debris flow; two-dimensional numerical simulation; geographical information systems (GIS); raster grid networks; finite difference method

宁夏沙坡头河段二维水沙运动数值模拟研究

宁夏沙坡头河段二维水沙运动数值模拟研究
水 力 发 电
第3 9 卷第 5 期
2 0 1 3 年5月
宁 夏 坡ห้องสมุดไป่ตู้头 河 段 二 维 水 沙 运 动
数 值 模 拟 研 究
杨 程 1 , 2 , 3 , 4 ,李春 光 , 一 ,景何 仿 2 , 吕岁 菊
( 1 . 宁 夏大 学土 木 与水 利 工 程 学 院 ,宁 夏 银 川 7 5 0 0 2 1 ;2 . 北 方 民族 大 学数 值 计 算 与工 程 应 用 研 究所 ,宁 夏 银 川 7 5 0 0 2 1 ;3 . 宁夏 节 水 灌 溉 与 水 资 源调 控 工程 技 术 研 究 中心 , 宁夏
Y A N G C h e n g , L I C h u n g u a n g , J I N G He f a n g  ̄ , L U S u i j u
( 1 . C o l l e g e o f C i v i l a n d Hy d r a u l i c En g i n e e r i n g , Ni n g x i a Un i v e r s i t y . Yi n c h u a n 7 5 0 0 2 1 , Ni n g x i a , C h i n a ; 2 . Re s e a r c h I n s t i t u t e o f
Nu me i r c a l C o mp u t a t i o n a n d E n g i n e e r i n g A p p l i c a t i o n , B e i f a n g Un i v e r s i t y f o r Na t i o n a l i t i e s , Yi n c h u a n 7 5 0 0 2 1 , Ni n g x i a , Ch i n a ; 3 . Ni n g x i a R e s e a r c h C e n t e r o f T e c h n o l o g y o n Wa t e r - s a v i n g I r r i g a t i o n a n d Wa t e r Re s o u r c e s Re g u l a t i o n , Yi n c h u a n 7 5 0 0 21 , N i n g x i a , C h i n a ; 4 . En g i n e e r i n g Re s e a r c h C e n t e r or f E ic f i e n t Ut i l i z a t i o n o f Wa t e r Re s o u r c e s i n Mo d e m Ag i r c u l t u r e i n A r i d

水中泥沙颗粒聚集的二维数值模拟——单体凝聚模型

水中泥沙颗粒聚集的二维数值模拟——单体凝聚模型

t e r . eef c fsikn r b b l is t kn o iin,mo in t  ̄e t r n a tce h o y Th fe to tc i g p o a i t ,si i g p st ie c o to r co ya d p ril
n um b r c nc nt a i n h v e n ob a n d. e o e r to a e b e t i e Num e i a i u a i n c n be us d t t y s di rc lsm l to a e o s ud e —
第 7卷 第 4期 21 0 0年 1 2月
长 沙 理 工 大 学 学 报 (自 然 科 学 版 )
J u n l f h n s aU ie i f c n ea dT c n lg ( au a S i c ) o r a o a gh nv r t o i c n eh oo y N tr l c n e C sy S e e
TAN a — h n,LI W ncu AO a — i IJeh a Xio x ,L i— u ,LIHa — a g, oy n
W A NG n b Yu - o,SUN — u n,W U n —on Shiq a Fa g t g
( c o lo d a l g n e i g, a g h i e st fS in e a d Te h o o y Ch n s a 4 0 0 Ch n ) S h o f Hy r u i En i e rn Ch n s a Un v r i o ce c n c n l g , a g h 1 0 4, i a c y
me tp r i l g r g t n s c e s u l . e f c o s wh c a o e s u i d i x e i e t n a tc e a g e a i u c s f l Th a t r i h c n n t b t d e n e p rm n o y a c r t l a e c n r l d Th s i n e f c i e s p l m e t a d n w y f r r s a c f c u a ey c n b o t o l . e i s a f e tv u p e n n e wa o e e r h o

河流水流泥沙运动数值模拟技术研究

河流水流泥沙运动数值模拟技术研究

图 1-11 单连通域网格的生成 1.2.3 网格拼接方法 网格拼接是把分散的区域网格拼接起来,形成整体,以便于完整网格节点坐标的输 出。 下图为网格拼接后图形。
A D B
1
C A A
A
D B
2-1
D C C C B
2-2
D B

图 1-12 2.2.4 网格节点数据的导出
网格拼接实例
网格节点信息的导出形式有三种: (1).dxf 形式; (2).grd 形式; (3).shp 形式, 这三种形式网格信息可以分别用 CAD、Delft3D 和 GIS 软件直接打开。 2.2.5 地形插值计算 网格数据导出前可以对地形进行插值计算,从而得到网格节点处地形数据“地形插 值计算”界面如图 1-13 所示。
图 2-2
二维水沙数值模拟软件系统
二维水沙数值模拟软件系统数据流图如图 2-3 所示。
现场数据
软件生成 用户输入
输入文件
程序调用
前处理
求解器
软件生成 计算生成
后处理
用户操作
计算结果
图 2-3 软件数据流示意图
输出文件
交通部天津水运工程科学研究所
10
河流水流泥沙运动数值模拟技术研究 简本
2.1.1 TK-2DC 的构成 TK-2DC 包括前处理模块、中心计算模块(求解器) 、后处理模块。其中,中心计 算模块包括水流计算模块和泥沙计算模块,均采用 FORTRAN90 语言进行编制,其源程 序可以在 Fortran 编译系统中编译成 *.exe 可执行文件,该文件可以脱离 Fortran 系统单 独执行。 2.1.1.1 前处理模块 (1) 网格划分 网格划分是前处理中的重要部分,通过用户给定计算域散点地形高程,网格划分模 块程序还可以将地形插值到网格节点上,完成网格地形文件的生成。详见子题一报告。

泥沙局部冲淤二维数值模拟仿真

泥沙局部冲淤二维数值模拟仿真

泥沙局部冲淤二维数值模拟仿真张芝永;拾兵【摘要】为了更好地对河流、海洋的局部冲淤问题进行研究,在FLUENT中通过二次开发建立了泥沙局部冲淤二维数值模型.模型采用标准k-ε紊流模型来计算流场,通过自定义函数功能将泥沙输运模型嵌入到FLUENT中,床面地形的变化采用动网格技术来进行实现.通过2个算例对该模型的可靠性及通用性进行了验证.结果表明,该模型模拟结果与前人结果基本符合,模型的通用性较好.该模型可以应用于其他形式的局部冲淤数值模拟,具有广阔的应用前景.%In order to study the local erosion and deposition in a river or an ocean, a two dimensional (2-D) numerical model on local erosion and deposition is established through the redevelopment in FLUENT code. The flow field is calculated by applying standard equations. By using a user defined function, the sediment transport model is embedded into FLUENT and besides, the bed profile is simulated by utilizing a dynamic mesh. The study used two experiment cases to verify the reliability and universality of the scouring model. The research results show that the predictions from the model are in agreement with the prior results, and the scouring model is so universal that it may be used to simulate local scour in different cases. Therefore, it has a broad application prospect.【期刊名称】《哈尔滨工程大学学报》【年(卷),期】2013(034)002【总页数】6页(P145-150)【关键词】局部冲淤;FLUENT;自定义函数;动网格【作者】张芝永;拾兵【作者单位】中国海洋大学工程学院,山东青岛266100【正文语种】中文【中图分类】P751泥沙冲淤是河流、海洋中比较普遍的一种自然现象,在近床面位置存在结构物情况下,床面泥沙运动剧烈,局部地形变化强烈,局部冲淤问题严重.对于此类问题,国内外许多学者主要是通过模型实验手段来进行研究分析,其研究对象主要包括桥墩[1]、丁坝[2-3]、海底管道[4-5]等结构物周围的局部冲淤问题.随着计算机技术的发展,局部冲淤的数值模拟方法越来越引起研究者的重视.Ouillon等[6]采用标准的三维k-ε紊流模型模拟了丁坝周围的流场,并对平整床面时的床面切应力沿河床的分布与最终冲刷坑的床面形态进行了定性比较.彭静等[7]将线性与非线性三维k-ε模型应用于丁坝绕流模拟,并比较了各种模型的模拟效果.Brφrs[8]提出了新的二维冲淤模型采用标准的k-ε紊流模型和同时考虑悬移质和推移质输运,对海底管道局部冲刷过程进行了模拟.Liang和Cheng [9-10]基于有限差分法发展了一个可以精确地预测海底管道冲刷时间历程的二维模型.Liu[11]利用VOF和动网格技术对水平射流情况下的底床冲淤进行了数值模拟,该模型考虑了水面的变化情况,冲坑变化过程模拟值与试验值吻合较好.韦燕机[12]利用OpenFOAM和动网格技术对桩周局部冲刷进行了三维数值模拟.Zhi-wen Zhu[13]利用网格技术对桥墩冲刷进行了数值模拟,验证了网格较好的地形边界拟合效果.综上所述,以上大部分模型均是通过动网格技术来模拟泥沙的局部冲刷,这说明动网格技术要比其他如两相流技术来模拟泥沙更为准确.但上述模型也存在一些不足,如大部分局部冲淤模拟均是通过非通用计算程序来完成的,这限制了其数值模拟方法的推广.FLUENT是一款目前较为流行的商用CFD软件,它在航空航天,能源、水利、环境等领域得到了广泛的应用,是目前全球功能最强大的计算流体力学商用软件.该软件在水动力计算模块有着广泛的通用性,但该软件并没有利用动网格来计算泥沙冲淤的模块.有鉴于此,作者充分利用FLUENT的二次开发接口,即自定义函数功能,通过C++编写了泥沙输运模块,并将其嵌入到FLUENT中,实现了泥沙局部冲淤的二维数值模拟.1 水动力模型1.1 控制方程水动力模块中,其控制方程为连续性方程和动量方程:式中:ui为流体速度在i方向上的分量;ρ是流体的密度;μ是流体粘度;p是作用在流体微元上的压力;v是流体的运动粘滞系数,Sij是平均应变速率张量.ui'是i方向速度脉动值为雷诺应力张量,vT紊流运动粘滞系数,δij是克罗内克尔符号.1.2 湍流模型对于不可压缩流体且不考虑用户自定义的源项时,标准k-ε模型方程为k方程: ε方程:式中:Gk是由于平均速度梯度引起的湍动能k的产生项;C1ε和 C2ε为经验常数,取 C1ε =1.44,C2ε =1.92;σk和σε分别是与湍动能k和耗散率ε对应的 Prandtl 数,取σk=1.0,σε=1.3.对于床面近底层,采用壁面函数法将壁面上的物理量与湍流核心区求得物理量联系起来.2 泥沙输运模型2.1 推移质输沙模型床面泥沙的起动可以通过临界希尔兹数来确定,其希尔兹数的表达式为式中:τ为床面切应力,ρ为水的密度,g为重力加速度,s为泥沙比重,d50为泥沙中值粒径.水平床面上泥沙临界起动希尔兹数为其中,D*是无量纲泥沙颗粒尺寸,定义为在有坡度的床面上泥沙临界起动希尔兹数按式(10)进行修正[14]:式中:α是砂床面和水平面的坡角,规定上坡为正角,下坡为负角;φ是泥沙的休止角.平面单宽体积输沙率q0由式(11)求得推移质单宽体积输沙率可用下式表示:式中:q0为平面单宽体积输沙率,τ为床面剪应力,τx为床面剪应力x方向分量,h为床面高程,经验系数 C 取值范围为 1.5 ~2.3[8].2.2 床面变形模型根据输沙平衡,床面高程h的变化可以用冲淤方程表示为式中:p0为床沙孔隙率,qbx为x方向推移质单宽体积输沙率.3 动边界处理方法3.1 动边界处理方法动网格模型可以用来模拟流场形状由于边界运动而随时间改变的问题.其广泛应用于活塞、阀门及柔性体等有运动边界工况的模拟.床面的冲淤变形形式类似于柔性体变形,需要对各个节点的运动情况进行描述.因此本文采用FLUENT的DEFINE_GRID_MOTION宏命令来给出边界节点的运动方式.在床面节点位置更新后,需要对区域内部网格进行调整.FLUENT提供了3种网格更新方式:1)光顺网格法,2)动态层网格法,3)局部网格重画法.非结构化网格比结构化网格更加适用于不规则边界的变形问题,因此本模型采用非结构化网格,其网格更新方法为光顺网格法和局部网格重画法相结合的方式,这2种方法的网格更新效果可见图1.图1为下文算例中水流下坡面处局部网格变化情况.从图中可以看出,在底部边界发生变化后,区域内网格进行重新划分,各区域网格尺寸尺度基本保持不变,靠近下部变形边界的网格仍然比较密集,这对于确保近壁面流场计算的准确性至关重要. 图1 动网格更新效果Fig.1 Mesh deformation3.2 沙滑模型由于泥沙传输的复杂性和与流场相互作用的非线性,在海床高程更新过程中会发生反射,导致网格畸变,造成不存在的冲淤形态和数值的不稳定.为解决此问题,采用Liang提出的沙滑模型[10]:当局部冲淤斜坡角度超过泥沙休止角时,对相应的海床节点进行调节,使冲淤斜坡角等于休止角.4 数值解法及计算流程水动力控制方程及紊流方程使用基于单元中心的有限体积法离散,利用非稳态求解器进行求解,瞬态项采用二阶隐格式,对流项采用二阶迎风差分格式,扩散项采用中心差分格式;压力与速度的耦合使用SIMPLE算法.床面变形方程式(13)通过有限差分法进行离散,其离散后形式为式中:p代表床面节点序号;n为当前时间步;n+1为下一时间步;Δtb为床面更新时间步长.为节省计算时间,流场计算的时间步长和床面更新的时间步长采用不同值,床面更新时间步长Δtb要远大于流场计算时间步长Δtflow.其数值模拟过程可概括为以下几个过程:1)计算多个步长的速度、压力、湍流数及切应力等;2)利用最后一步得到的床面切应力数据,调用泥沙输运模块,计算推移质输沙率和床面高程变化情况;3)采用动网格技术,改变床面节点位置,同时重新调整计算区域内部网格;4)返回步骤1),重新上面的计算,直至到达指定时间.5 模型验证为验证本文所建模型的可靠性,通过2个算例对其进行了验证.5.1 水渠冲淤Van Rijn[15]对有坡度水渠的水流情况进行了试验研究,李昌良[16]在其基础上又对其冲淤情况进行了进一步的研究.其试验布置方案如图2所示.左边进口平均流速0.51 m/s.水深0.39 m,沙床泥沙平均粒径d50=0.16 mm.在数值模拟中,其计算区域与试验布置相同,左边进口设置为速度入口,右边为自由出流边界,水面边界设置为对称边界.底面边界设置为wall边界,其粗糙度为2.5d50.整个区域采用非结构化网格进行离散.图2 试验布置Fig.2 Layout of experiment图3列出了在初始地形情况下4个垂直断面的流速分布情况,从图中可以看出数值模拟结果与试验结果是较为一致的,这说明了数值模拟中水面采用刚盖假设的合理性,同时为进一步的泥沙冲淤过程研究奠定了基础.图4为不同时段的网格情况及床面地形变化情况,可以发现在床面地形发生变化后,采用光顺网格法和局部网格重画法更新后的非结构化网格很好地拟合了床面地形的变化,在变化过程中,网格的疏密程度也得到了较好的保持.图5给出了3 h后床面冲淤情况.从图中可以看出,在沟内由于水深的增加,水流流速减小,进而导致床面切应力减小,泥沙在此淤积.在水流下坡区域,床面切应力减小,泥沙淤积较为严重,而在水流上坡区域,下游流速加大,床面切应力增大,从而导致该区域冲刷严重.数值结果与前人结果比较一致,准确地反映了沟内床面的变化规律.图3 各断面流速分布情况Fig.3 Velocity distributions at different crossvsections图4 不同时刻网格更新情况Fig.4 Mesh deformation in different time图5 3 h后床面地形Fig.5 Bed profile at t=3.0 h5.2 海底管道局部冲刷本算例对近床面的海底管道局部冲刷的冲刷历程进行详细的分析.Liang在文献[17]中对有间隙海底管道的局部冲刷进行了详细的数值模拟研究.在这里采用其中一计算算例进行了数值计算,算例中管道直径D为10 cm,管道与床面之间垂直距离G为0.5D,来流流速为0.5 m/s,泥沙中值粒径为 0.36 mm,水深 3.5D.选择距管道中心上下游各20D距离的区域为计算区域.左边界为速度进口边界,右边界为自由出流边界,上部边界采用对称边界,下部边界为wall边界(如图6所示).整个区域采用非结构化网格进行离散,在近壁面处和管道周围进行加密.图6 计算区域Fig.6 Computational domain图7 不同时刻x方向流速等值线(单位:m/s)Fig.7 Contours of x veolicity indifferent time(unit:m/s)由于动网格计算耗时较长,本文只模拟了80 min内的床面冲刷过程.图7列出了冲刷过程中,不同时刻的x方向流速分布图,该图同样直观地描述了床面的冲刷变形情况.在t=5 min时,由于初始冲刷阶段间隙内流速较大,管道局部冲刷程度较为严重,管道正下方出现沙坑,下游区域出现沙丘.在t=15 min时,上游冲坑深度加大,下游区沙丘逐渐向下游移动并并且沙丘高度逐渐减小.在t=80 min时,下游沙丘消失,整个过程与实际的冲刷过程规律是较为一致的.从图中还可以发现,管道下游并未出现涡旋脱落现象,这可能是由网格尺度的限制以及启动动网格模型而导致的计算精度下降等原因造成的.而在实际试验中,管后区域是存在涡旋脱落现象的.Liang[17]的研究结果表明,有涡旋脱落时,管道后尾流区范围比较大而且变化比较剧烈,大范围的尾流区使得更多水流通过管道下方孔道流向下游,进而导致近床面流速加大,管下和管后床面切应力也相应增大,同时,涡旋的不断脱落又使得管道下方及后方床面的切应力波动变化,波动幅值有可能会是平均值的数倍.因此,有涡旋释脱落时的冲刷深度要大于无涡旋脱落情况.图8列出了t=80 min时,Liang的模拟结果与本文模拟结果的对比情况.从图中可以看出,管道上游的床面冲刷模拟结果较为接近,管道下游区域床面冲刷模拟结果与Liang的无涡旋脱落时的冲刷结果比较接近,而与有涡旋脱落时的模拟结果则有较大偏差.其原因在于本文模拟过程中并未出现涡旋脱落现象,从而造成管道下方及后方床面切应力计算值偏小,冲刷程度较小.图8 t=80 min时床面地形Fig.8 Bed profile at t=80 min6 结论通过对商用CFD软件FLUENT的二次开发,建立了局部冲淤二维数值模型.该模型的泥沙输运模型通过FLUENT的DEFINE_GRID_MOTION宏命令嵌入,床面边界的冲淤变化利用动网格技术来实现.数值计算的结果表明,泥沙冲淤模型可以比较准确的预测不同条件下泥沙冲淤过程,可用于模拟以推移质输沙为主的二维泥沙冲淤问题,具有广阔的应用前景.但该模型同样存在一些问题,如对于有涡旋脱落时的泥沙冲淤模拟偏差较大以及计算时间过长等问题.这都有待于进一步优化和研究.参考文献:【相关文献】[1]MELVILLE B W.Pier and abutment scour:integrated Approach[J].Journal of Hydraulic Engineering,1997,123(2):125-136.[2]彭静,河原能久.丁坝群近体流动结构的可视化实验研究[J].水利学报,2000,31(3):44-47.PENG Jing,YOSHIHISA Kawahara.Visualization of flow structure around submerged spur dikes[J].Journal of Hydraulic Engineering,2000,31(3):44-47.[3]ROGER A,KUHNLE C V,ALONSO F,et al.Local scour associated with angled spur dikes[J].Journal of Hydraulic Engineering,2002,128(12):1087-1093.[4]SUMER B M,JENSEN H R,FREDSØE J.Effect of leewake on scour below pipelines in current[J].J Waterway,Port,Coastal and Ocean Engineering,1988,114(5):599-614. [5]CHIEW Y M.Prediction of maximum scour depth at submarine pipelines[J].Journal of Hydraulic Engineering,1991,117(4):452-466.[6]OUILLON S,DARTUS D.Three-dimensional computation of flow around groyne [J].Journal of Hydraulic Engineering,1997,123(11):962-970.[7]彭静,河源能久.线性与非线性紊流模型及其在丁坝绕流中的应用[J].水动力学研究与进展:A辑.2003,18(5):589-594.PENG Jing ,KAWAHARA Yoshihisa.Application of linear and non-linear turbulent models in spur dike flow[J].Chinese Journal of Hydrodynamics,2003,18(5):589-594.[8]BRØRS B.Numerical modeling of flow and scour at pipelines[J].Journal of Hydraulic Engineering,1999,125(5):511-523.[9]LIANG D F,CHENG L,LI F.Numerical modelling of scour below a pipeline in currents.part I:flow simulation[J].Coastal Engineering,2004,52(1):25-42.[10]LIANG D F,CHENG L,LI F.Numerical modelling of scour below a pipeline incurrents.partⅡ:Scou r simulation[J].Coastal Engineering,2004,52(1):43-62.[11]LIU Xiaofeng,GARCÍA M H.Three-dimensional numerical model with free water surface and mesh deformation for local sediment scour[J].J Waterway,Port,Coastal and Ocean Engineering,2008,134(4):203-217.[12]韦雁机,叶银灿.床面上短圆柱体局部冲刷三维数值模拟[J].水动力学研究与进展:A辑,2008,23(6):655-661.WEI Yanji,YE Yincan.3D numerical modeling of flow and scour around short cylinder[J].Chinese Journal of Hydrodynamics,2008,23(6):655-661. [13]ZHU Zhiwen,LIU Zhenqing.CFD prediction of local scour hole around bridge piers [J].Journal of Central South University of Technology,2012,19(1):273-281.[14]ALLEN J R L.Simple models for the shape and symmetry of tidal sandwaves:statically stable equilibrium forms[J].Marine Geology,1982,48(12):31-49. [15]VAN RIJN L C.Principles of sediment transport in rivers,estuaries and coastal seas [M].Amsterdam:Aqua Publications,1993:1245-1247.[16]李昌良.泥沙运动与底床变形的数值模拟[D].青岛:中国海洋大学,2008:58-61.LI Changliang.The numerical simulation of sediment transport and bed evolution[D].Qingdao:Ocean University of China,2008:58-61.[17]LIANG Dongfang,CHENG Liang,YEOW K.Numerical study of the Reynolds-number dependence of two-dimensional scour beneath offshore pipelines in steady currents[J].Ocean Engineering,2005,32(4):1590-1607.。

关于三峡库区二维泥沙冲淤计算可视化的思考

关于三峡库区二维泥沙冲淤计算可视化的思考

关于三峡库区二维泥沙冲淤计算可视化的思考摘要自1994年以来,我国不断加大了三峡水利工程的建设力度,以起到改善水土流失、防洪抗汛等作用,有助于保护长江沿岸人民的生命财产安全,拉动经济增长。

但在实际发展中,三峡库区受多方面因素的影响,库区内泥沙淤积量显著提高,对库区的蓄水性能、经济效益、安全性造成了严重的影响,需要对三峡库区的冲淤量进行计算,便于调整。

而随着科学技术的发展,对冲淤计算方法有了更高的要求,应向现代化、信息化、可视化方向发展,发挥其综合效益价值。

关键词三峡库区;泥沙冲淤计算;可视化1 浅析三峡库区概况三峡库区属于亚热带季风气候,且周围多山地、丘陵,导致该地区降水较多,能满足三峡蓄水需要。

同时,该地区江河纵横分布,加上降水的影响,河流汛期较长,水量充沛。

自从2003年三峡库区正是蓄水以来,库区泥沙冲淤量不断增加,使得河流特性、河床的稳定性等受到了影响。

近些年来,根据各个控制站的观测数据可知,截止2013年底,三峡库区入库泥沙总量达到了20.3亿吨,包括悬移质泥沙、泥沙淤积等,存在严重的安全隐患。

同时,当长江排汛时,会将库区中的部分泥沙量流入到下游库区当中,进而整体抬高洪水位,降低库容。

因此,需要工作人员做好泥沙冲淤计算工作,及时了解三峡库区的泥沙淤积量,合理制定清理计划,确保三峡库区正常运行。

但由于传统的计算方法存在很多限制性问题,使得计算结果的真实性、准确性、实时性等受到了影响,需要基于计算机网络技术,开发专业的泥沙冲淤计算系统,确保计算过程、结果的可视化,保障三峡库区的控制效果[1]。

2 探究三峡库区二维泥沙冲淤计算可视化系统开发策略2.1 开发程序和语言VB工具是Windows系统中的重要内容,该工具的开发语言面向整个系统,通过进行语言编程,有利于形成不同的可视化程序窗口,便于工作人员进行操作。

通常情况下,利用VB来开发系统,最基础的可视化内容为界面外观。

同时,还需要利用对象程序设计方法来将程序与数据连接在一起,形成可视化对象或数据。

基于Vb建立的二维潮流泥沙数值模拟界面

基于Vb建立的二维潮流泥沙数值模拟界面
2 1 0 0 0 9) ( 南 京 市 河 海 大 学 港 口海 岸 与近 海工 程 学 院 ,江 苏 南 京

鉴 于 计 算机 技 术 在 工 程 中 的广 泛 渗 入 ,文 中应 用 v b 为主 生成 界 要 :数 值 模 拟 技 术 广 泛 应 用 于 工 程 中各 方 面 ,
( VB S c r i p t ) 是广泛使 用的脚本语言 ,它是 Vi s u a l B a s i c 语言 的子 集 。从开发个人使用或小组使 用的工具 ,到大 型企
业 应 用 系 统 ,甚 至 通 过 I n t e r n e t的 遍 及 全 球 分 布 式应 用程 序, 都 可 在 Vi s u a l Ba s i c 提 供 的工 具 中各 取 所 需 。 Vb作 为 软件 编 程 中界 面 显 示 的先 行 者 ,有 着 非常 大 的 优势 。制 作界 面使 用 VB 的 优 点 :
作 以 及显 示 上 。
本文 采 用 v i s u a l b a s i c 进 行 软 件 界 面 的 开 发 。 因 为
v i s u a l b a s i c在 界 面 显示 方 面拥 有 很 大 的优 势 。
“ V i s u a l ”指 的 是开 发 图形 用户 界 面
行 近 似 求 解 以显 示 某 种 过程 ,与 物 理 模 型 相 比 ,数 值 模 拟 具
语言 。 Vi s u a l Ba s i c在原 有 B AS I C 语言 的基础上进一步发
展 ,至今 包含 了数 百条语 句、 函数及 关键 词 ,其 中很 多和 Wi n d o ws GUI有直接 关系 。 专业人员可 以用 Vi s u a l Ba s i c 实现其它任何 Wi n d o ws 编程语言 的功能 , 而初学者只要掌

黄河河口海岸二维非恒定水流泥沙数学模型

黄河河口海岸二维非恒定水流泥沙数学模型

黄河河口海岸二维非恒定水流泥沙数学模型曹文洪,何少苓,方春明(中国水利水电科学研究院泥沙研究所)摘要:针对黄河河口海岸岸线变化剧烈和含沙量变幅大的特点,开发和建立了适合黄河河口海岸应用的平面二维动边界非恒定水流泥沙数学模型。

验证表明,本模型可以较好地模拟黄河河口海岸泥沙输移和冲淤变化,为研究和解决多沙河口海岸的泥沙问题提供技术手段。

关键词:黄河口;挟沙能力;窄缝法;非恒定流;数学模型收稿日期:2000-01-06基金项目:国家重点基础研究发展规划项目(G1*******)资助作者简介:曹文洪(1963-),男,(满族),黑龙江省人,中国水利水电科学研究院教授级高工,博士。

自本世纪七十年代以来,由于计算机技术的迅猛发展,国内外相继出现了众多的河口海岸泥沙数学模型[1-7],有力地促进了河口海岸的泥沙研究的发展。

然而,已有的河口海岸数学模型大多是模拟含沙量较低的河口海岸的泥沙运动,而能够模拟多沙和岸线延伸剧烈的河口海岸泥沙数学模型还极为少见。

近年来,已有个别学者尝试用泥沙数学模型模拟黄河河口海岸的泥沙运动,如张世奇开发了一套黄河口平面二维泥沙冲淤数学模型,得到了较好的效果[18~20]。

为了全面系统地反映黄河三角洲海陆动态交互影响机理和泥沙运动与湿地演替关系,本文开发和建立了径流、潮流和波浪作用下的黄河河口海岸平面二维动边界非恒定流非均匀沙不平衡输沙数学模型。

1 模型结构1.1 水流运动基本方程(1)(2)(3)为谢才式中:U、V分别为潮流速在x及y方向的垂线平均值分量;Z为潮位;Cf系数;F为柯氏系数,F=2ωsinφ,式中ω为自转角速度,φ为地理纬度;h为水深;Z为海底起始高程。

b1.2 潮流和波浪共同作用下泥沙运动方程在河口海岸地区,潮流和波浪是泥沙运动的最主要动力。

“波浪掀沙和潮流输沙”使河口海岸地区的泥沙运动极为活跃,但也更为复杂。

因此,在计算和预报河口海岸地区的泥沙运动和冲淤变化时,仅仅考虑潮流的作用是不全面的,还应考虑波浪的作用。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
及 深入 地 系统研 究 二维 水 沙运动 基本 规律 . 是 建
象. 并 利 用 数 值 计 算 进 行 近 似 求 解 以复 演 某 种 自 然 或人 工 过程 这些 在严 世 强. 熊德 琪 I l l 对潮 流 数 值 模拟 研究 进展 中均有所 提及 。与物 理模 型相 比。

l 0一
吉林 水利
二 维 潮流 泥沙 数值 模 拟及 可视 化研 究
朱嘉 玲等 2 0 1 3年 1 O月

O h+ 生 + 生 : 0
O t O x a v
计 算 区域 内 的任 意位 置 及相 应 模拟 时 间段 以 内 的 任意 时刻进 行 潮 流及 泥 沙 运动 复演 .同时 还 可对 港 口、 航道 、 围垦 、 码 头 等水 利 及 海 洋 工 程 实 施 之 后 的潮 流及 泥 沙 运 动进 行模 拟 研究 。 主界 面分 为
第 1 0期 ( 总第 3 7 7期 )
[ 文章编号]1 0 0 9 — 2 8 4 6( 2 0 1 3 )1 0 — 0 0 1 0 — 0 6
吉 林 水 利
2 0 1 3年 l 0月
二维 潮流 泥 沙数 值 模 拟 及 可 视 化 研 究
朱嘉 玲 .石 然 2 1 0 0 0 9 ) ( 河海 大学 港 口海岸 与近 海工程 学院 ,江 苏 南京
是 二 维 泥沙 数 学模 型 问题 至 今 尚未 能很 好 地 解 决 有 两 点 重要 原 因 .其 一是 由于 泥 沙运 动 规 律极 为 复 杂 .对 冲积 河道 的水 流 和泥 沙运 动 规 律 认 识 不
足, 很 多 问题 , 诸 如输 沙 率 、 阻 力 等 问题 , 均 没 有 较 好 地 解 决 现 有 的二 维 泥沙 数 值 模 拟大 都 是 直 接 引用 一 维计 算 公 式 .对 于 这样 做 是 否适 合 很 少 涉
[ 摘 要 ]本 研 究在 以往 研 究 成 果 基础 上 , 以 水 动 力 学基 本 理 论 与 软 件 工 程 理 论 相 结 合 ,研 究 开发 了以 V i s u a l B a .
s i c编 程 语 言 为依 托 的 二 维 潮 流 泥 沙 数 值 模 拟 及 可视 化研 究模 型 , 共 分 为三 个相 互 连 接 的界 面 , 实现 了潮 流 泥
立 二 维 泥沙 数 学模 型 的主要 任 务 之 一 其 二 是 由 于河 道 边 界极 不 规 则 .在 计算 中要 求 划分 的 网格 较密 . 所需 计 算量 及 贮 存量 太 大 . 目前 常 用 的有 限 差 分法 及 有 限元 法 难 以较 好 地 解决 这 一 问 题 寻 找 计算 量 及 计算 贮 存 量较 小 且适 用 于不 规 则 河道 边 界 的计 算 方法 .是 二维 泥 沙 数学 模 拟 急 待解 决
的问题 。
数值模拟具有经济 、 快速 、 修 改方 便 等 优 点 , 并 且
不 受 比尺 的限制 . 已广 泛应 用 于解 决 各种 海 岸 、 河 口动力 问题 的研 究. 如 河 口整 治工 程 、 环境 保 护工 程 等 此外 . 还可 以用 来揭 示 许 多海 岸 、 河 口物理 现象, 例如 : 通 过潮 流 、 泥 沙数 值模 拟 , 可 了解 水质 点在 任 意 时刻 的流 速 、 流 向和 水 量 的输送 、 湍 流 的 强 度 以及 涡旋 、环 流状 况 和泥 沙 输 移 以及 底 床 变 化等。
环境 模 型界 面层 次 清晰 . 同时对各 方 案 的计算 结 果 进行 直 观 、 动态 显示
四个 因素 : f 1 )人 类 生产 活动 和 发展 的需 求 : ( 2 ) 计 算 数 学 的发 展 , 即提供 控 制方 程 更 多 、 更方 便 的离
散 求 解方 法 : ( 3 ) 计算 机 c p u 、 内存 、 运 算 速度 、 图像 制 作 技 术 等 的发展 ; ( 4 ) 河流 海 洋 资料 的测 量 精准
潮 流 与泥 沙数 值 模 拟 的不 断发 展取 决 于 以下
模 拟 开发依 托 Wi n d o w s 技术 平 台 . 采用 V i s u a l
B a s i c编程语 言 进行 可视 化 程序设 计 .是 一个 集 输
入、 输出、 计算、 结 果 显 示 等 功 能 为 一 体 的 可 化效 果 图显 示 ,操 作 简便 ,可 运 用 性 强 。 [ 关键 词 ]二 维 ;泥 沙 ;数 值 模 拟 ;可视 化 ;v b
[ 中 图分 类 号 ], I T v 1 4 9
[ 文 献 标识 码 ]A
与 进展 及 河 口水 动 力 泥 沙数 值模 拟 进 行 总 结 但
[ 作 者 简 介 ] 朱 嘉玲 ( 1 9 9 2 一 ) , 女, 江 苏 宜 兴人 , 本科 专业 : 港 口航 道 与 海 岸 工 程 。
[ 基 金 项 目 ] 国 家创 新 训 练 项 2 0 1 2 0 5 x e x 0 2 9 ]项 目名 称 : 二 维 潮 流泥 沙数 值 模 拟 及 可视 化 研 究 。
度。
2 基 本 原 理
平 面 二维 潮 流泥 沙 数 值模 拟 所 依 据 的基 本 方 程 可归结 为 :
水 流连 续方 程 式
现今 . 国 内外很 多学 者对 此 展 开 研究 . 李 孟 国
l 2 1

曹 祖德【 , 】 等 分别 海 岸河 口潮 流 数值 模拟 的研 究
[ 收 稿 日期 ]2 0 1 3 — 0 5 — 2 1
1 概 述
潮流与泥沙数值模拟始于 2 0世 纪 6 0年 代 . 是一 门新兴 的综 合 性模 拟 技 术. 它 以流 体 力学 、 泥 沙 动 力 学 以 及计 算 数 学 等 学 科 为 基 础 . 结 合 具 体
工程技术. 通 过 采 用 数 学 模 型 来 模 拟 某 种 物 理 现
相关文档
最新文档