圆形数值波浪水池斜向和多向入射波浪模拟研究
水波实验报告
水波实验报告水波实验报告引言:水波实验是一种经典的物理实验,通过观察水波的传播和相互干涉现象,可以帮助我们更好地理解波动现象的特性和规律。
本次实验旨在通过搭建水波实验装置,观察和分析水波的传播、干涉和衍射现象,并探讨实验结果与理论预期的一致性。
实验装置:我们使用了一个长方形水槽作为实验装置,水槽内部填满了适量的水。
在水槽的一端,我们固定了一个发生器,通过激发水面上的振动源来产生水波。
在水槽的另一端,我们安装了一个接收器,用于观察和记录水波的传播情况。
为了更好地观察水波的细节,我们在水槽上方安装了一台高速摄像机,以便捕捉水波的运动轨迹。
实验过程:1. 单一波源实验:我们首先将发生器设置为单一波源模式,即只有一个振动源激发水面上的波动。
在这种情况下,我们观察到波浪以同心圆的形式从波源处向四周传播。
通过调节振动源的频率和振幅,我们发现波浪的传播速度和振幅大小之间存在一定的关系。
2. 双波源干涉实验:接下来,我们将实验装置调整为双波源模式,即在水槽的两端同时设置两个振动源。
在这种情况下,我们观察到两组波浪从波源处同时向四周传播,并在某些区域发生干涉现象。
通过调节两个振动源的频率和相位差,我们可以观察到干涉现象的变化,包括增强干涉和相消干涉。
3. 衍射实验:最后,我们将实验装置调整为衍射模式,即在水槽的一端设置一个振动源,而在另一端设置一个障碍物。
在这种情况下,我们观察到波浪从振动源处向障碍物传播,并在障碍物后形成衍射现象。
通过调节振动源的频率和障碍物的形状,我们可以观察到衍射现象的变化,包括衍射角度和衍射图案的形态。
实验结果与讨论:通过以上实验,我们得出了以下结论:1. 水波的传播速度与波长和频率有关,符合波动方程的理论预期。
2. 在双波源干涉实验中,当两个波源的频率相同且相位差为整数倍时,会出现增强干涉现象;当相位差为奇数倍时,会出现相消干涉现象。
3. 在衍射实验中,障碍物的形状和波源的频率对衍射现象的角度和形态有显著影响。
斜向和多向不规则波在斜坡堤上的单波越浪量的实验研究
块体 护面 , 前者 共有 三 个 坡 比( 1:3 1: , , 2 1:1 5 .,
图 l)后者 只有 一个坡 比( 15 图 l ) a目 : 等 学 校博 士 学 科 点 专 项 科 研 基 金 资 助 项 目(0 3 1 10 ) 高 2 00 4 0 2 .
小 , 以工 程设计 中通 常 要选 择 某 个 累积 概 率很 小 所 ( 0 0 1 的单 波 越 浪量 来 作 为 设 计 值. 如 .0) 由于 越 浪
收 稿 日期 :0 60 — 7 修 订 日期 : 0 6 1 —O 20 —71 ; 2 0 — 03 .
实验共 设置 了一 个 堤 顶 宽度 ( . 4m) 两 个 堤 0 1 、 顶高 度 ( . 8 . 6m) 三 个 水 深 ( . 0 . 6 , 0 5 ,0 5 、 0 4 ,0 4 6 0 4 5m) . 8 以及 两种 护 面 形式 : 凝 土 护 面 和扭 工 字 混
作 者 简 介 : 晓 亮 (9 9 )男 , 宁 省 抚 顺 市 人 , 士 生 , 究 方 向 为 不规 则 波 越 浪 量 的 研 究 . — i ll su e tdu.d .n 李 1 7一 , 辽 博 研 E ma : x@ td n. lteu c l
维普资讯
l 引 言
出于经 济 和安全 等 因素 的考 虑现 在 国 内外对海 岸建 筑物 多按 允许 发 生 越 浪 的标 准 来 进 行设 计 , 越
的机 理 十分复 杂 , 现有 的估 算 越 浪 量 的 方法 几 乎 都 是基 于物 理模 型实 验 的结果 . S u等 提 出不 规 则 波 单 波 越 浪 量 的概 率 分 布 符合 Web l分 布 , iu 1 之后 该 结 论 又 被 不 同 的 实验 研 究l 2 所证 实. Web l分 布 函数 为基 础通 过进 行 。 以 iu1
数值波浪水池中
图3 Wigley III 船模垂荡运动时的垂荡力时历 (T=2.0s)
采用同样的计算和分析方法,船体分别作单一的垂荡、纵摇和纵荡运动,计算分析可得垂荡、纵摇和纵荡的无因次的附加质量和阻尼系数。如图4~图6所示,图中离散点值为本文CFD理论计算结果,连续线为笔者根据文献[1]的原理编程计算所得的线性势流结果。基于CFD方法计算各模态,垂荡运动位移幅值为0.02m,纵摇运动幅值为0.02rad,纵荡运动幅值为0.04m。
近年来,随着计算机技术和计算技术的飞速发展,计算流体力学(CFD)也有了长足的进步。基于CFD理论的船舶水动力学方面的数值模拟,因为具有费用低、无触点流场测量、无比尺效应、能消除物模中由传感器尺寸及模型变形等因素对流场的影响、可获得较为详细的流场信息等优点而广受关注,且应用的范围越来越广[8-11]。
表1 Wigley III 模型主尺度
0.1875 0.6667 0.4622 4.2 船体水动力系数的计算分析和验证
本文就Wigley III 船模分别所作的"纯垂荡"、"纯纵摇"及"纯纵荡"运动下的流场进行了数值模拟, 以验证本文方法的有效性和准确性。当船体在数值波浪水池中以固定的频率和幅值作单一模态的强迫振荡运动时,通过流场的模拟计算,可以获得相应模态下的力和力矩,以及静力和静力矩。如图3所示,当船模以周期为2.0s,位移幅值为20mm作"纯垂荡"运动,即以的位移规律运动时,通过模拟计算获得整个船模的垂向静力和垂荡力的时历曲线。计算中采用了固定的时间步长约为振荡运动周期的六百分之一。图3中的横坐标为时间,纵坐标为计算所得的垂向静力及垂荡力。对力和力矩的时历曲线进行相位分解。由于CFD模拟的计算结果都总存在噪声,用最小二乘原理以正弦曲线对平稳的时历曲线拟合,进而获得船体的附加质量和阻尼系数。对应于图3的时历,附加质量与阻尼系数分别为:=113.784,=140.795。计算发现垂荡和纵摇之间耦合的附加质量和阻尼的量级比垂荡的小6~7个量级,这与Wigley数学船型关于横中剖面对称是一致的。
数值波浪水池中船舶顶浪运动模拟研究
对船 舶运动 的数值 预报也 相 当重要 ; 最后 就是 波浪 中船 舶运 动的数 值模 拟 。
文献(】 数值 波 浪 水池 的定 义及 其 发展 历史 和研 究 现状 作 了较 为详 细 的介 绍 , 1 对 并且 进行 了规则 波 和不规则 波的模 拟 , 获得 了满足 进一 步数 值计 算要求 的 波浪环 境 , 同时 还对 规则 波 中顶浪前 进 的拘 束船模 进行 了数 值模 拟 , 得到 的船 模在 波浪 中的水动力 预报 结果 与 D T的试验 数据 吻合 良好 。所有 U
Ab ta t N meia i lt n o a e id c d s i t n n rg lrh a ig wa e n ie a 3 i— sr c : u rc lsmu ai fw v -n u e hp moi si e ua e dn v si s D vs o o d
第l 2卷第 5期
20 0 8年 1 0月
文 章编 号 :1 0 — 2 4(0 8) 5 0 9 — 5 0 7 7 9 20 0 — 6 2 0
船 舶力 学
J u n lo hp Me h nc o r a fS i c a is
V0 .2 1】 No5 .
Oc . 0 t 2 08
的数 值计算都是 基 于 N S方程 的 。 —
文献 [ 则基 于 N S方程 , 航行船 舶辐 射 问题 进行 了数值模 拟研 究 。文 中计 算 了不 同航 速 、 同 2 ] - 对 不
频 率 下船模 垂荡 和 纵摇 的力 、 矩 、 力 附加质 量及 阻 尼 系数 , 与 D T的 试验 数 据 进行 了 比较 , 者符 并 U 二 合 良好 。
s i fn me ia i l t n a e c mp r d w t h x e i n a a a f m l v r i f e h o o y u t o u r ls s c mu ai r o a e i t e e p rme tld t r Def Un e st o c n lg , o h o t i y T a d s o ar o d a r e n t h m. S b s d NWT sc n t ce y t e r s a c r . n h w f i y g o ge me t h t e An N- a e l wi r i o sr t d b h e e r h wok u Ke r s y wo d :N- q a in ; e d n a e ; v - n u e t n S e u t s h a i g w v s wa e i d c d moi s o o
工程流体力学中的波浪与海浪模拟探讨
工程流体力学中的波浪与海浪模拟探讨工程流体力学是研究流体在各种工程领域中的力学行为的学科。
其中,波浪和海浪作为海洋工程中的重要现象,对于海岸工程、沿海建筑物设计等领域具有重要意义。
如何准确模拟波浪和海浪对于这些工程的设计、建设和维护至关重要。
本文将探讨工程流体力学中的波浪和海浪模拟方法。
波浪模拟是模拟频率较低的周期性波浪,如在湖泊、港口和与海洋接触的近海水域中观测到的波浪。
常用的波浪模拟方法有线性波浪理论和非线性波浪模型。
线性波浪理论是基于近似假设的波浪模拟方法。
它假设波浪幅度较小,且海洋介质为无限深度和均匀密度。
线性波浪理论可以通过解线性波浪方程获得波浪的水平位移、速度和加速度等参数。
然而,这种方法无法准确模拟高振幅和非线性波浪,因此在实际工程中应用受到限制。
非线性波浪模拟是一种更加准确的方法,它将波浪的非线性效应考虑在内。
基于非线性波浪理论的模型可以模拟各种频率和振幅的波浪,因此在复杂海况下具有广泛应用。
其中最常用的非线性波浪模拟方法是第三代数值模型,如SWAN、MIKE21等。
这些模型基于数值方法求解非线性波浪方程组,可以模拟包括海浪传播、折射、反射、干涉等现象,以及海浪对工程结构的冲击力。
海浪模拟是模拟海洋波浪的过程。
海浪是复杂的现象,受到风力、地球自转、地形等多种因素的影响。
海浪模拟的目的是预测和理解海浪的性质和行为,以便进行海洋工程的设计和规划。
在海浪模拟中,常用的方法包括统计方法、天气预报模型和数值模型。
统计方法是通过对海浪历史数据的统计分析来预测未来海浪的情况。
这些方法通常基于长期观测得到的海浪数据,并使用统计学方法进行分析和建模。
然而,由于统计方法无法考虑到具体的海洋环境和地形,其预测结果存在一定的不确定性。
天气预报模型是基于大气环流和海洋动力学的模型。
这些模型可以通过对大气条件和海洋环流的模拟,预测海浪的生成和演化过程。
天气预报模型的优点是可以通过实时数据输入进行预测,但其精度受到模型参数和边界条件的限制。
波浪波浪模拟实验
3.16波浪波浪模拟实验海洋、湖泊、水库等宽广的水面下可能发生较大的波浪,波浪将影响船舶进航和停泊的安全;对堤防闸坝以及其他岸边建筑物具有强大的冲击作用;也会引起近岸浅水地带水底泥沙运动;淘刷岸坡和护岸建筑的基础,使航道和港口受淤、岸坡崩溃。
因此波浪成为水力学的重要课题之一。
一、实验目的:1、掌握波浪水槽和造波机的基本结构和原理2、掌握描述波浪基本要素的定义和测量3、了解波浪水槽模型实验的基本方法二、实验装置:1、水槽:1、电脑2、液压系统和造波板3、浪高仪4、消波器2、造波机:造波机安装在波浪水槽一端,造波机后侧直立式消能网,水槽的另一端设有消能坡消除波浪反射影响。
整个造波系统由造波板、液压伺服作动器、液压泵站、伺服放大器、AD/DA接口及计算机与外设等部分组成。
实验中要模拟一个波谱时,首先根据目标谱利用傅利叶变换将其展开成一个电压时间序列值控制信号,经D/A接口转换成不规则的电压信号,送给伺服控制放大器,驱动造波机推板作相应的推挽运动,推动水体而产生波列,位移传感器实时测出推板的运动轨迹,实时反馈到控制放大器,修正机械惯性带来的误差,以确保推波板能准确地跟踪计算机给定的信号运行。
造波的同时,浪高仪将波浪物理量转换成电量信号送A/D转换器进行数据采集,并暂存于内存中,供谱分析。
一般情况下,每次谱模拟不得少于120个波。
由于传递函数拟合时产生的误差及机械系统的影响,很难一次模拟成功,必须按以下公式修正。
[])(ϖα−=ωϖ+S*ϖ(DS(S))(S)式中:---------------修正后的控制谱S*ϖ()-----------------实测模拟谱)(S ϖ α -------------------修正系数)(DS ϖ --------------目标谱按重新计算出电压时间序列值,再一次控制造波机造波,分析比较,直至得到理想的模拟谱为止。
一般情况下经过二到五次修正就可基本成功。
)(S *ϖ三、实验原理:1、波浪要素的定义:波浪现象的特征是水的自由界面出现有规律的起伏形态,水的质点则作有规律的振荡运动,同时波形以一定速度向前传播。
基于CFD的规则波顺浪数值水池模拟
基于CFD的规则波顺浪数值水池模拟近年来,众多研究人员越来越多地将计算流体力学(CFD)应用于波浪水池的数值模拟中。
这种方法基于纳维-斯托克斯方程组和运动比例分析,使用分析方法来描述流体的运动、静力学和流动物理量。
利用CFD技术可以建立非常准确的数值模型,从而对各种复杂的物理现象进行研究,并为设计新型波浪水池提供基础实验数据。
规则波顺浪数值水池,是一种基于CFD技术的数值水池模拟,用于研究顺浪波的反射、折射以及波浪相互作用等问题。
该模型通常由多个二维波浪水池组成,可以模拟出多种不同波形及频率的波浪。
这种模型可以极大地提高波浪水池实验的效率和精确度,同时可以避免一些实验中不可控因素,如风、池面污染等。
数值模拟的关键在于模型的建立。
在建模的过程中,需要将波浪水池外的空间分为不同的层,分别表示空气、水、沙子等物质。
通过计算流体力学的理论公式,可以获得层间的物理数量关系,并建立起有效的计算模型。
在数值模拟中,波浪水池的边界需要被精确建模,以确保水池中的波浪表现与实验所见一致。
同时,水池中液体的物理参数,如形态、粘性和密度等,也需要被精确建模。
这些因素都对波浪水池的模拟结果产生很大的影响。
模型建立完成后,便可以进行不同波形波频的波浪模拟。
在模拟中,波浪水池主要分为两个区域,分别为波浪产生的区域和波浪衰减的区域。
在波浪水池的产生区域,可以通过输入不同波形或不同频率的激励,启动不同种类的波浪。
在波浪衰减区域,可以模拟不同波形波频下的波浪衰减以及波浪与边界的相互作用。
在模拟过程中,所得到的波浪高度、波速以及波长等参数,均可以用于分析和研究目的,如波浪对海岸线的破坏、波浪对结构物的影响以及海洋工程设计等。
总的来说,规则波顺浪数值水池模拟是一种基于CFD技术的新型水池数值模拟方法,可以有效地模拟不同波形波频下的波浪,具有高效、精准和可控等特点,适合于海洋工程、造浪研究等领域的应用。
在规则波顺浪数值水池模拟中,以下数据是最为重要的参数:1. 波高:波高是指波浪峰和波谷之间的垂直距离。
物理实验技术中的波浪现象实验指导
物理实验技术中的波浪现象实验指导波浪现象是自然界中广泛存在的一种物理现象,也是物理学中一个重要的研究领域。
在物理实验技术中,我们经常使用波浪现象来进行实验和研究。
本文将为您介绍一些常见的波浪现象实验指导。
一、水波实验在物理实验中,水波实验是常见且容易进行的一种波浪实验。
通过在水面上制造波浪,我们可以观察和研究波浪的特性和传播规律。
1. 实验器材准备进行水波实验时,我们需要准备以下器材:水槽、振荡器、振荡器控制器、水波测量仪等。
2. 实验步骤(1)将水槽放置在平稳的实验台上,并保证水面平整。
(2)将振荡器安装在水槽底部,并连接振荡器控制器。
(3)根据实验需求,设置振荡器的频率和振幅。
(4)观察水波的传播和现象,可以使用水波测量仪等工具进行测量和记录。
3. 实验注意事项在进行水波实验时,需要注意以下事项:(1)确保实验环境的稳定,避免外界干扰。
(2)水槽和器材的清洁,以保证实验的准确性和可靠性。
(3)根据实验要求调整振荡器的频率和振幅,避免过大或过小,以免造成实验结果的偏差。
二、光波实验光波实验是物理实验技术中另一个重要的波浪实验。
通过对光的干涉、衍射等现象的研究,我们可以了解光的波动特性和传播规律。
1. 干涉实验干涉实验是通过两束光的叠加来观察光波的干涉现象。
常见的干涉实验有双缝干涉、薄膜干涉等。
(1)双缝干涉实验准备:准备双缝实验装置,包括光源、双缝装置、光屏、测量工具等。
(2)双缝干涉实验步骤:①将光源放置在一定距离处,光源越亮度越好。
②在光源后面放置双缝装置,调整缝隙的大小和距离。
③将光屏放置在双缝后方,调整距离和观察角度。
④观察光屏上的干涉条纹并记录。
2. 衍射实验衍射实验是通过光波通过障碍物或过孔时产生的衍射现象来观察光的波动性。
(1)单缝衍射实验准备:准备单缝实验装置,包括光源、单缝装置、光屏、测量工具等。
(2)单缝衍射实验步骤:①将光源放置在一定距离处,光源越亮度越好。
②在光源后面放置单缝装置,调整缝隙的大小和距离。
模拟水波反射_折射和衍射现象的解析模型的研究进展
刘焕文
( 广西民族学院 计算机与信息科学学院, 广西 南宁 530006)
摘 要 : 本文对模拟水波反射、 折射和衍射现 象的解析模型的研究进展给出了一个简单的综述 . 关键词 : 水波; 折射; 衍射; 解析模型 中图分类号 : O1751 11 文献标识码: A 文章编号: 1007- 0311( 2004) 02- 0073- 06
图1
矩形障碍物连接直线斜坡对水波的反射
当 s 2 = 1, x 0 v ] 时, 经代数推导, 公式( 1) 将退 化成 KR = 11+ h 1/ h 0 h 1/ h 0 ( 2)
这正是 Lamb( 1932) 关于长波被无限长大陆架 作用 的结果 . 当 s 1 = s 2 , 即 h 1 = h 2 时, 公式 ( 1) 又将退化成 v0 L 1 + v1 L 0 - 2 v 0 v 1 - 8 s1/ P / A 0 v0 L 1 + v1 L 0 - 2 v 0 v 1 + 8 s1/ P / A 0 ( 3) 其中, v 0 = J 0( A 0) + J 1( A 0) , v 1 = J 0( A 1) + J 1( A 1) , L0 = Y 0 ( A 0) + Y 1( A 0) , L1 = Y 0 ( A 1) + Y 1( A 1) , v 0 = J 0 ( A0 ) Y 0 ( A 0) + J 1( A 0) Y 1( A 0) , v 1 = J 0 ( A1 ) Y 0 ( A 1) + J 1( A 1) Y 1( A 1) . 公式( 3) 正是 Dean( 1964) 的当大陆架与直线斜坡连 接时长波方程的解析解( 见 p19- 20) . 当 s 2 = 1, x 0 v x 1 时, 公式( 1) 还将退化成 KR = ( s 1 - 1) | sin k 1 a | ( s 1 + 1) - ( s 1 - 1) cos k 1 a
水波演示实验报告
水波演示实验报告篇一:大学物理演示实验报告水波盘【实验目的】利用水波的投影显示波的形成、传播、反射、干涉和衍射等的形象。
【实验器材】水波盘演示仪,如图20-1所示。
有水槽、振动源、光源、各种振子(包括单振子、双振子、平面波振子)及挡板2块水槽及壳体水槽是用底部装有密封、透明玻璃的不锈钢盆制成。
壳体用金属材料制成,上面放有水槽,正面竖直安装毛玻璃,作为水波投影的屏幕。
框架内部倾斜45°装有平面镜,用来反射水面的影象到屏幕上,底部装有变压电源,后面装有一立杆。
立杆上端安装光源盒,中部安装振动源盒,在立杆的中部开有长槽孔,用来调节振动源盒的高度。
振动源振动源采用电磁、激励式。
它是由电磁铁、电位器、振杆、振子、主板等组成。
振频调节是一个与电磁铁线圈串联的可调电阻,控制其电流以改变振动的频率。
调节振幅螺丝,可使投影波形的清晰度达到最佳。
振动源盒后面有一插孔,使用时与光源盒插头相接。
光源光源为盒式机械遮挡频闪光源,灯泡为12V100w幻灯机溴钨灯,频闪器是由直流微型电动机驱动一个可旋转的遮挡叶片,使其遮光次数为50?60次/秒。
盒的顶部开有散热窗,当电机旋转时,带动遮光叶片,驱动盒内热空气排除盒外,使之降温。
实验原理】两列同频率、振动方向相同、相差恒定的波在某一区域相遇,会产生干涉现象。
有的地方振动始终干涉加强,有的地方干涉减弱。
理论计算表明,干涉加强与干涉减弱由下式确定:(20-1)其中,【实验操作与现象】1.准备工作,为波长。
(a)在水槽内注3?8mm深的清水,充分湿润水槽四周及实验用的附件。
(b)将振荡波所需的振子固定在承接块上,调节振动源盒的高度,使振子插入水面1-2mm。
(c)先把光源电源插头、直流电源插头及振源插头插入到相对应的插座中,再插上交流220伏输入电源插头。
2.实验操作(a)圆形波将单振子固定在承接块上,打开电源开关,溴钨灯亮,遮光叶轮转动,振杆弹动,屏幕上即显示图象,为圆形波,如图20-2所示。
波浪水池消波装置选型的试验研究
波浪水池消波装置选型的试验研究兰 波, 缪泉明, 姚木林,胡定健,向 旭(中国船舶科学研究中心,江苏 无锡 214082 )摘要:主要考虑到不同消波模型消波效果的不同,对直壁墙、铁丝网阵模型、箱式模型以及圆弧形模型的消波效果进行试验比较研究。
通过试验得到不同消波装置模型的反射系数,并对其进行比较,得出最理想的消波装置。
计算程序采用了五种不同的理论计算方法,并将计算结果进行了比较。
本研究为水池消波装置的选型提供依据,并为该类型试验的试验方法提供依据。
关键词:消波装置;消波效果;反射系数;选型消波装置的作用在于减少水池波浪受反射波的干扰影响,所以各国水池对消波装置的研究都非常的重视。
对于不同的消波装置其消波效果是不同的,而消波效果的好坏会直接影响到波浪水池试验的精度。
目前有关港工方面的消波装置的试验研究资料很丰富,但是由于它们对消波装置的要求比较低,其所要求的反射系数不能满足船舶耐波性/海洋工程水池的要求,所以有必要对消波装置进行系统的研究,特别是对各种不同消波装置进行试验比较,为水池消波装置的选型提供依据。
1 试验原理消波装置的消波效果可以用反射系数的大小进行评估。
对于规则波,其反射系数R K 可以表示为反射波波幅R A 与入射波波幅I A 的比值,即/R R I K A A =。
而对于不规则波,其反射系数()R K f 可以用反射波的谱密度()R S f 与入射波的谱密度()I S f 来表示:()R K f 的时间历时,并对其进行分离计算来得到反射系数。
2 模型试验试验是在二维波浪水槽中进行,水槽长40m ,宽0.8m ,深1.4m ,试验水深为0.9m 。
试验一端装有推板式造波机,可以产生规则波、不规则波。
在水槽中装有5个电容式浪高仪,用来测量消波装置前的波浪波高。
布置如图1所示。
试验主要有三种试验模型:铁丝网阵模型、箱式模型、圆弧形试验模型。
铁丝网阵消波装置由5块透空率为49%的铁丝网(参见图 2)根据一定间距布置在远离水槽造波机的一端。
数值波浪水槽的造波及消波方法
数值波浪水槽的造波及消波方法李晖;何宏舟;杨绍辉【摘要】In order to improve the efficiency of solving marine engineering related problems, and to pro⁃vide reference for the design of the physical wave generator, the FLUENT software was used to simulate the process of wave generation and absorption in a wave tank. In the simulation, the motion conditions of the swinging⁃plate wave generator were defined by the user defined function UDF and the fluid motion was realized by employing dynamic mesh technology. The reflection of the wave was eliminated by the method of setting porous media area and adding momentum source term of viscous drag force. The VOF method was utilized to simulate the free surface. The simulation results show that the swinging-plate wave generator proposed in this paper can be used to produce a stable linear wave in the numerical flume. The periodic error is small enough to meet the requirements of the simulation experiments. The wave absorption method of momentum source term with a viscous drag force in a porous media area is effective enough to eliminate the reflected wave from the end wall of the tank.%为了提高海洋工程相关问题的研究效率,并为物理造波机的设计提供参考,应用FLUENT软件模拟了数值波浪水槽的造波及消波过程。
波浪在斜坡堤上传播的数值模拟
波浪在斜坡堤上传播的数值模拟安蒙华;蒋勤;张长宽【摘要】基于VOF自由表面追踪技术,通过直接求解Navier-Stokes方程及κ-ε方程,建立了求解波浪与建筑物相互作用的数值波浪水槽,并将其用于研究波浪在不透水斜坡堤上传播的物理过程.首先,通过对水平床上非线性波传播过程的模拟分析验证了该数值波浪水槽的模拟精度;在此基础上,对规则波在斜波堤上传播的水动力过程,包括波浪爬坡、越浪及越流等进行了数值模拟,并通过物理模型试验结果对堤前水位及越浪量变化的模型模拟结果进行验证.【期刊名称】《水运工程》【年(卷),期】2014(000)006【总页数】6页(P25-29,35)【关键词】波浪;斜坡堤;爬坡;越浪;数值模拟【作者】安蒙华;蒋勤;张长宽【作者单位】河海大学,江苏南京210029;河海大学,江苏南京210029;河海大学,江苏南京210029【正文语种】中文【中图分类】TV139.2+6波浪传播到斜坡堤时,受斜坡堤底坡的影响,波浪剖面变陡,波峰质点速度增加,致使波浪发生破碎,并沿斜坡爬升。
当波浪的上爬高度超过坡顶时,即发生越浪越流。
正确评估斜坡堤上波浪的爬坡、越浪和越流与堤顶的设计及其安全稳定性密切相关。
因此,研究波浪在斜坡堤上的传播过程,特别是评估波浪的爬高和越浪量,是国内外学者研究的传统热点课题。
早在1955年,Saville T等[1]就针对规则波在斜坡堤上的爬坡以及越浪现象进行了物理模型试验研究。
此后,Owen[2]和Van Der Meer等[3]先后对斜坡堤的越浪量问题做了大量的研究工作。
20世纪90年代以后,Holger Schüttrumpf[4]及Van-Gent[5]通过物理模型试验研究给出了斜坡堤越浪量、越浪流的流速及流深的理论公式。
20世纪70年代后期,随着计算机技术的发展,数值波浪水槽的研究和应用受到人们的关注。
关于波浪在斜坡堤上传播的数值模拟,Hu等[6]基于非线性浅水方程,使用AMAZON模型对斜坡堤上波浪爬高及越浪过程进行了模拟,并与Saville T的试验结果进行了对比。
数值模拟推板式波浪水槽的关键参数取值探讨
数值模拟推板式波浪水槽的关键参数取值探讨关超;吴静萍;詹成胜;刘博宇【摘要】Numerical simulation is quite time-consuming for water wave problems due to unsteady flow. In order to solve this problem ,in addition to using the high-performance computer and parallel tech-nique ,it is necessary to optimize the mesh size and time step in the numerical simulation which play key roles in improving calculation accuracy and saving calculation time.In this paper ,assuming that the flow waslaminar ,the two-dimensional wave propagation of piston-type wave flume was numeri-cally simulated with variation of 7 kinds of mesh size in wave direction ,4 kinds of mesh size in wave height direction and 3 kinds of time steps.The accuracy verification was analyzed by the numerical and analytical solutions in the wave height error and degree of wave height attenuation in unit distance. Comprehensively considering the calculation accuracy and time consumption , the optimal mesh size and time step are recommended.The optimal mesh size is related to the aspect ratio of element ,and the aspect ratio is presented in this paper.Then compared history wave elevation of simulation results with experiment results ,they are agree in a satisfactory extent.%对于水波问题,属于非定常流动,数值模拟相当耗时.为了解决这个难题,除了采用高性能计算机和并行技术外,需要优选数值模拟中网格尺度和时间步长,它们在提高计算精度和节约计算耗时中起着关键作用.文中假设流动为黏性层流,改变了7种波长方向网格尺度,改变了四种自由面附近波高方向的网格尺度和改变了三种时间步长,对二维推板式波浪水槽中波浪传播展开了数值模拟.数值模拟的精度用数值计算和解析解的波高误差及单位距离波高相对衰减量来衡量.综合考虑计算精度和计算耗时,推荐了较优的网格尺度和时间步长的取值.较优的网格尺度还与单元网格长宽比有关,文中给出了较优网格尺度对应的单元网格长宽比.将推荐取值下数值模拟的时历波形与试验波形进行比较,二者吻合效果良好.【期刊名称】《武汉理工大学学报(交通科学与工程版)》【年(卷),期】2018(042)003【总页数】5页(P515-519)【关键词】数值模拟;二维推板式波浪水槽;网格尺度;时间步长【作者】关超;吴静萍;詹成胜;刘博宇【作者单位】武汉理工大学交通学院武汉 430063;武汉理工大学交通学院武汉430063;武汉理工大学交通学院武汉 430063;中交机电工程局有限公司船舶及贸易事业部北京 100088【正文语种】中文【中图分类】U656.20 引言随着计算机技术和计算流体动力学(computational fluid dynamics,CFD)的发展,从20世纪60年代开始至今,计算流体动力学取得了丰硕的成果,出现了众多的CFD商业软件,如PHOENICS,CFX,FLUENT,STAR CCM+,FLOW-3D,FIDIP等,使得数值模拟方法以其成本低、周期短、精度高等优势在流体力学学科的科学研究和工程项目中已经得到认可和被推广使用[1-2].数值模拟方法在波浪运动、波浪与结构物相互作用问题中的应用,因自由面捕捉技术、造波技术、消波技术、动网格技术等关键技术起步相对较晚,与物理波浪水槽并称为数值波浪水槽,在海岸工程、港口与航道工程、船舶与海洋工程专业的波浪问题研究中起着重要作用.数值波浪水槽造波有3种方法:仿物理(推板式或摇板式)造波法[3-4]、边界造波法[5]和源造波法[6-7].对于水波问题,属于非定常流动.推板式造波适用于浅水波及中等水深波,而摇板式造波则适用于深水波.数值模拟非定常的水波问题,计算相当耗时,除了必须采用高性能计算机和并行技术之外,其网格尺度和时间步长对减少计算时间和保证计算精度起着关键作用.网格尺度、时间步长的取值问题一直受到研究者重视.郭晓宇[8]对源造波的数值波浪水槽进行验证时采用的波参数为波长L=2 m,波高H=0.05 m,水平方向网格尺度取L/20和L/40两种,垂向波高范围内网格尺度取H/8,结果显示两者与解析波形均吻合良好;此外又通过一系列具有不同非线性程度的波浪输入条件进行验证,波长分别为L=30,60和90 m,波高分别为H=1,3和4 m,水平方向网格尺度为L/120,L/180,L/240,发现随着波浪非线性程度的增大,单位波长需要较多网格节点才可以满足足够的精度.查晶晶等[9]用OpenFOAM实现数值波浪水槽造波,采用推板式造波法模拟波高H=0.02 m,周期T=1 s的线性波,所选网格尺度在30 m水池长度范围内的单元数为600,800,1 000,1 500,在垂向波高范围内的单元数为10,20,30,时间步长选取0.008,0.005,0.002,0.001 s进行探讨.通过与解析解比较,最终确定纵向网格、垂向网格和时间步长为1 500,20,0.002和800,10,0.002为最优参数,并得出结论,计算结果的精度主要受波面区域矩形单元的长宽比影响,网格单元越细长,引入的数值黏性越大,由此造出的波与目标波波幅相差太大.Finnegan等[10]采用CFX构建二维摇板式数值波浪水槽,模拟线性深水波与有限水深波,讨论网格依赖性,推荐了最小网格尺度计算式;推荐的时间步长为T/50.管陈等[11]采用FLUENT软件实现摇板式造波,模拟波长L=2 m,波高分别为H=0.04,0.06和0.1 m的规则波,选取RNG 雷诺时均模型,时间步长约为0.001 s,自由液面附近网格最小尺度为波高的1/10~1/20,水平方向网格尺度沿池长方向逐渐加大,并将数值结果与解析解进行比较,发现同一波长下,波高越小,其误差越小,模拟精度越高.此外,虽然一些学者研究数值波浪水槽时未对相关数值造波技术参数进行依赖性的讨论,但在文献中给出了其取值.严汝建等[12]使用FLUENT软件基于摇板式造波,数值模拟波参数为L=8 m,H=0.2 m的规则波,水平方向和自由液面附近垂向网格尺度分别为L/20和H/10,并选取RNGk-ε湍流模型.李宏伟[13]使用FLUENT建立二维数值水槽,采用摇板式造波法模拟参数为L=8 m,H=0.4 m,T=2.2678 s的规则波,水槽长度方向上网格尺度为L/100,自由液面附近进行加密,时间步长为0.001 s,采用RNGk-ε湍流模型.黄华等[14]利用推板式造波法模拟二维二阶Stokes波,波浪参数为L=5.783 2 m,H=0.129 m,T=2 s,给出的网格尺度水平和竖直方向分别约为L/111和H/21,时间步长为T/400.廉静静等[15]采用FLUENT软件,通过摇板式造波模拟参数为L=8 m,T=2.284 8 s的规则波时,时间步长选为0.001 s,采用k-ε湍流模型.本文以某流体力学实验室的小型推板式波浪水槽为原型,基于商业软件FLUENT建立推板式造波的二维数值波浪水槽.选取波浪水槽有代表性的波浪参数造波.考虑水波理论多数基于势流理论,而真实的流体具有黏性,所以假设本文模拟的波浪运动为层流流动.首先讨论了沿波长方向不同网格尺度、自由面附近沿波高方向上不同网格尺度以及不同时间步长对计算精度和计算耗时的影响,通过将数值结果与解析解的比较,综合考虑波高误差、单位距离波高相对衰减量和计算耗时,推荐综合性能较优的网格尺度、时间步长;然后将基于推荐的网格尺度、时间步长的数值结果与试验结果进行比较,两种吻合程度相当令人满意.1 数值方法的基本控制方程本文中二维波浪流动属于黏性、不可压缩和非定常类型,其波浪传播方向为x轴正向,波浪高度方向向上为z轴正向,则流体连续性方程为(1)动量方程为(2)(3)式中:u为x方向速度;w为z方向速度;t为时间;ρ为流体密度;p为压强;ν为运动黏度;g为重力加速度.采用VOF方法捕捉自由液面,控制方程为(4)式中:F为体积分数的函数,F=1为单元体全为水,F=0为单元体全为空气,0<F<1为单元体处于空气与水的交界面上.2 数值计算方法本文建立二维数值波浪水槽.水槽长为19 m(后4 m为消波段),高为0.8 m,造波水深为0.5 m.生成目标波形的参数为L=1.5 m,T=0.995 1 s,H=0.06 m(为波峰致波谷之间的高度差),属于二阶Stokes波.采用建模软件Gambit进行数值波浪水槽的模型建立及网格划分,网格均为结构化网格.水槽底部设置边界层网格,自由液面附近进行网格加密,加密总高度为波高的2倍.数值波浪水槽与物理波浪水槽一样,需采要取一定的措施消除水槽末端池壁反射波的影响.数值消波技术主要采用人工阻尼区[16]、辐射边界条件[17]和主动吸收[18],或者在数值波浪水槽末端使用稀疏网格进行数值耗散[19].本文在消波段采用渐变的稀疏网格耗散波能,从而消除水槽末端反射波的影响.数值波浪水槽上方设为压力出口边界条件;两端及底部设为壁面边界条件,其中左边模拟物理推板造波,赋予推板一个UDF函数,采用layer的动网格方法,使其按照给定的方式进行运动.数值波浪水槽模型和边界条件见图1,网格划分见图2.图1 数值波浪水槽模型和边界条件图2 数值波浪水槽网格划分3 数值波浪水槽模拟计算3.1 时间步长变化时间步长的选取会影响数值计算结果的准确度和计算耗时.本文仅考虑波浪传播方向及自由液面附近网格尺度,首先取计算网格尺度Δx=L/70,Δz=H/10,考虑黏性,采用层流模型,取Δt=T/2 000,T/1 000,T/500三种时间步长进行数值模拟.图3为t=19 s时沿池长方向的波形空间分布,z为波面抬高.表1为x=5 m与x=8 m处时间步长Δt变化时数值波高与解析波高比较的相对误差,以及单位距离波高相对衰减量.图3 时间步长变化对数值模拟波浪传播的影响表1 时间步长变化对数值波浪传播的影响Δt网格单元长宽比H1/m(x=5 m)波高相对误差|H1-H|H×100%H2/m(x=8 m)波高相对误差|H2-H|H×100%单位距离波高相对衰减量/%T/2 0003.570.05488.670.053 710.500.61T/1 0003.570.056 46.000.05626.330.11T/5003.570.057 14.830.058 03.330.50注:H1为x=5 m处数值模拟的波高,H2为x=8 m处数值模拟的波高.单位距离波高相对衰减量采用表达式计算.由表1可知,随着时间步长的增大,数值结果与解析解的误差反而更小.由图3的波形及单位距离波高相对衰减量来看,当Δt=T/1 000时,波形的稳定性相对较好.综合考虑计算精度和波形稳定性,Δt=T/1 000是合适的时间步长选择.3.2 网格尺度变化网格尺度会对计算结果的精度产生影响,因此需对计算网格进行依赖性测试.首先讨论波长方向上网格尺度对数值造波的影响.取自由面附近波浪高度方向上网格尺度为Δz=H/10不变,波浪传播方向(水平方向)上网格尺度分别取Δx=L/50,L/60,L/70,L/80,L/90,L/100,L/110.图4为在t=19 s时波形的空间分布.表2为x=5 m与x=8 m处网格尺度Δx变化时数值波高与解析波高比较的相对误差,以及单位距离波高相对衰减量.图4 水平网格尺度变化对数值模拟波浪传播的影响表2 水平网格尺度变化对数值模拟波浪传播的影响Δx网格单元长宽比H1/m(x=5 m)波高相对误差|H1-H|H×100%H2/m(x=8 m)波高相对误差|H2-H|H×100%单位距离波高相对衰减量/%L/505.000.05449.330.053 910.170.28L/604.170.055 87.000.055 57.500.17L/703.570.056 46.000.056 26.330.11L/803.130.057 05.000.056 95.170.06L/902.780.057 24.670.057 14.830.05L/1002.500.057 44.330.057 34.500.06L/1102.270.057 54.170.057 44.330.05由图4与表2可知,当波浪传播方向上网格尺度Δx=L/80时,数值结果与解析解的误差已比较小.保持Δx=L/80不变,取波浪高度方向上、在自由面附近加密的网格尺度分别为Δz=H/10, H/20, H/30, H/40.图5为在t=19 s时波形的空间分布.表3为x=5 m 与x=8 m处网格尺度Δz变化时数值波高与解析波高比较的相对误差,以及单位距离波高相对衰减量.图5 竖直网格尺度变化对数值模拟波浪传播的影响表3 竖直网格尺度变化对数值模拟波浪传播的影响Δz网格单元长宽比H1/m(x=5 m)波高相对误差|H1-H|H×100%H2/m(x=8 m)波高相对误差|H2-H|H×100%单位距离波高相对衰减量/%H/103.130.05705.000.056 95.170.06H/206.250.056 85.330.05646.000.22H/309.380.056 85.330.056 46.000.22H/4012.500.056 95.170.056 46.000.28由图5和表3可知,波浪高度方向上网格尺度为Δz=H/10时数值计算结果基本趋近解析解.一味减小Δz并不能增加计算精度,要考虑合适的网格单元长宽比.综合上述工作,本文波浪数值模拟的计算网格尺度为Δx=L/80,Δz=H/10已经能够达到较好的精度.4 与物理试验结果比较依据数值模拟推板式波浪水槽波浪传播的关键参数取值分析,所推荐的网格尺度Δx=L/80,Δz=H/10和时间步长Δt=T/1 000,选取层流模型,数值模拟波长为1.5 m、波高为0.06 m的波,将此数值计算结果与物理试验比较,进一步验证该数值算法的可靠性.物理试验水槽为武汉理工大学流体力学实验室小型推板式波浪水槽,长为18 m、宽为0.6 m、高为0.8 m、水深为0.5 m.试验采用伺服电机控制推板造波,在水槽x=5 m与x=10 m处各布置一个浪高仪,用以测量波高随时间的变化.在数值计算中,同样在x=5 m与10 m处设置虚拟浪高仪,得出随时间变化的波形图.图6为数值模拟结果与试验结果的比较.表4为x=5 m与10 m处数值平均波高与试验平均波高的比较及相对误差(图6方框内平均波高).图6 数值结果与试验结果的比较表4 数值结果与试验结果的比较H1/m(x=5 m)H2/m(x=10 m)试验结果 0.057 20.055 9数值结果 0.056 50.055 3相对误差/%1.221.07由图6与表4可知,数值模拟的波形与试验波形吻合程度很好,其变化趋势表现出良好的一致性,验证了该数值算法的有效性.5 结论1) 本文数值波浪水槽中推荐采用:波浪传播方向网格尺度取Δx=L/80,自由液面附近波高方向网格加密尺度取Δz=H/10,时间步长取Δt=T/1 000;网格尺度确定还要满足一定的单元长宽比.2) 数值模拟的结果与物理试验结果吻合程度很好,验证了本文数值算法的有效性. 在展开新型浮式防波堤衰减长周期波性能研究时遭遇物理波浪水槽长度尺度限制的问题.本文工作为发挥数值波浪水槽的优势,利用数值模拟手段辅助物理试验进行研究奠定了基础.也为读者应用数值波浪水槽进行研究提供有益参考.参考文献[1] 王福军.计算流体动力学分析:CFD软件原理与应用[M].北京:清华大学出版社,2004.[2] 韩占忠.FLUENT:流体工程仿真计算实例与应用[M].北京:北京理工大学出版社,2010.[3] HEMMING A,SCHFFER. Second-order wavemaker theory for irregular waves [J]. Ocean Engineering, 1996,23(1):47-88.[4] PRASAD D D, AHMED M R, LEE Y H, et al. Validation of a piston type wave-maker using Numerical Wave Tank [J]. Ocean Engineering,2017,131:57-67.[5] 刘霞,谭国焕,王大国.基于边界造波法的二阶Stokes波的数值生成[J].辽宁工程技术大学学报(自然科学版),2010,29(1):107-111.[6] LIN P, LIU P, L F. Internal wave-maker for navier-stokes equations models[J]. Journal of Waterway, Port, Coastal and Ocean Engineering, 1999,125(4):207-217.[7] WEI G, KIRBY J T, SINHA A. Generation of waves in boussinesq models using a source function method[J]. Coastal Engineering, 1999,36(4):271-299.[8] 郭晓宇.数值波浪水槽及其应用研究[D].上海:上海交通大学,2011.[9] 查晶晶,万德成.用OpenFOAM实现数值水池造波和消波[J].海洋工程,2011,29(3):1-12.[10] FINNEGAN W, GOGGINS J. Numerical simulation of linear water waves and wave-structure interaction [J]. Ocean Engineering, 2012,43(4):23-31.[11] 管陈,董国祥,金允龙.三维数值波浪水池造波技术研究[J].上海船舶运输科学研究所学报,2013,36(2):11-15.[12] 严汝建,庞永杰,李宏伟,等.深水池造波系统数值造波仿真研究[J].哈尔滨工程大学学报,2010,31(1):32-41.[13] 李宏伟.数值水池造波方法研究[D].哈尔滨:哈尔滨工程大学,2011.[14] 黄华,邓冰,陈昱松,等.数值波浪水槽构建与二阶Stokes波仿真[J].系统仿真学报,2012,24(1):227-231.[15] 廉静静,尹勇,杨晓,等.基于黏性流船舶数值波浪水池造波和消波方法研究[J].船舶力学,2013(1):56-62.[16] 赵西增,孙昭晨,梁书秀.非线性波浪水槽内波浪的产生和传播[J].船舶力学,2010,14(3):203-216.[17] ORLANSKI I. A simple boundary condition for unbounded hyperbolic flows[J]. J Comput Phys, 1976,21(3):251-269.[18] DEEPAK D P,MOHAMMED R A,YOUNG H L,et al. Validation of a piston type wave-maker using Numerical Wave Tank [J]. Ocean Engineering, 2017,131:57-67.[19] PARK J C, UNO Y, SATO T, et al. Numerical reproduction of fully nonlinear multi-directional waves by a viscous 3D numerical wave tank[J]. Ocean Engineering, 2004,31(11):1549-1565.。
斜向和多向不规则波在斜坡堤上越浪量的研究的开题报告
斜向和多向不规则波在斜坡堤上越浪量的研究的开题报告题目:斜向和多向不规则波在斜坡堤上越浪量的研究摘要:斜向和多向不规则波是海洋中常见的波浪形式,其在海岸工程和船舶工程中的作用十分显著。
本文将通过实验的方法,研究斜向和多向不规则波在斜坡堤上的波浪量,以期为海洋工程领域的相关设计提供有效的数据支持。
关键词:斜向波、多向波、不规则波、波浪量、斜坡堤1.研究背景海岸工程和船舶工程中,波浪的特性和波浪量是设计和工程施工过程中必须考虑的重要因素。
而不规则波则是海洋中常见的波浪形式,其能给海岸线和海上构筑物带来很大的冲击力和摩擦力,极大地影响海岸工程和船舶运行的安全性。
因此,对于不规则波的研究成为了海洋工程领域的热点之一。
而斜坡堤作为一种常用的海岸工程结构,其与海浪的相互作用也备受关注。
斜坡堤上不同波浪形态的波浪量对于斜坡堤的稳定性和工程设计都有重要影响。
因此,探究斜向和多向不规则波对斜坡堤的波浪量的影响,对于提高海岸工程和航运的安全性和可靠性有重要的意义。
2.研究目的本文旨在通过实验方法,研究斜向和多向不规则波在斜坡堤上的波浪量及其变化规律。
通过在实验过程中改变波浪的入射方向、波高等参数,获取不同条件下的波浪量数据,并进行分析和处理,得出斜向和多向不规则波对斜坡堤波浪量的影响规律,为海洋工程的相关设计提供参考数据。
3.研究内容本文的研究内容主要包括以下几部分:(1)斜向和多向不规则波的基本概念和特性。
(2)斜坡堤结构的基本概念和特性,以及其与波浪的相互作用。
(3)实验设计:确定实验参数,包括斜坡角度、波高、波浪入射方向等,设计实验方案。
(4)实验过程记录和数据分析:在实验过程中,记录不同条件下的波浪量数据,并进行分析和处理。
(5)结果展示和讨论:根据实验结果,分析斜向和多向不规则波对斜坡堤波浪量的影响规律,并探究其原因。
(6)结论:总结研究结果,得出结论,并展望今后可能的研究方向。
4.研究方法本研究采用实验方法,通过在水槽中放置斜坡堤,测量不同角度、不同波高、不同入射方向的斜向和多向不规则波在斜坡堤上的波浪量,并对数据进行分析和处理。
基于粘性流船舶数值波浪水池造波和消波方法研究
1 引
言
物理实验波浪水池因造价成本高 , 试验周期长 , 维修麻烦以及物理实验波浪水池很难观测到水池 中流体在船舶周围细微的变化等原因 , 越来越多的水动力研究者将兴趣转移到数值波浪水池 , 并取得 了不少成果 [2-4]。 船舶数值波浪水池是通过计算机仿真 ,尽可能逼真地模拟真实的物理实验波浪水池的 各种功能 ,并最终能够代替物理实验波浪水池完成相应的科学研究及工程设计等任务 [1]。 不断地完善 收稿日期 : 2012-06-26
2 VOF 方法
流体体积法 VOF (Volume of Fluid ) 是由美国学者 Hirt 和 Nichols 等 (1981 )[5]提出的 , 它是一种可以 处理任意自由面的方法 , 其基本原理通过网格单元的流体体积分数 Cq 来确定运动界面 , 追踪流体变 化 , 该方法可以处理自由面重构等强非线性现象 , 所需要计算时间更短 、 存储量更少 。 定义为单元内第
k0 lsinh k0 h -cosh k0 h +cosh k0 h-k0 l k0 h+sinh k0 h cos k0 h kn lsinh kn h -cosh kn h +cosh kn h-kn l kn h+sinh kn h cos kn h
(7) (8)
(6 ) 式中 , 等式右边第一部分为行进波 , 第二部分为造波板产生的驻波 , 驻波在离开造波板一定距离以 后可以忽略 ,则可知造波板前的波幅为
第 17 卷第 1-2 期 2013 年 2 月
文章编号 : 1007-7294 (2013 )01-0056-07
船舶力学
Journal of Ship Mechanics
Vol.17 No.1-2 Feb. 2013
数值水池短峰不规则波模拟研究_冯光
第14卷第4期船舶力学Vol.14No.4 2010年4月Journal of Ship Mechanics Apr.2010文章编号:1007-7294(2010)04-0347-08数值水池短峰不规则波模拟研究冯光1,吴乘胜1,郑文涛1,顾民1,朱德祥2(1中国船舶科学研究中心,江苏无锡214082;2上海市计算技术研究所,上海200040)摘要:基于粘性数值波浪水池技术,对短峰不规则波进行数值模拟。
文中首先改进了长峰不规则波的数值模拟方法,并使用改进后的方法进行了不同海况下长峰不规则波数值模拟,模拟效果明显改善,数值模拟结果与目标值/谱吻合相当好。
之后进行了短峰不规则波的数值模拟,模拟结果与特征值/目标谱也相当接近。
关键词:短峰不规则波;长峰不规则波;数值模拟;RANS方程;VOF方法中图分类号:U661.2文献标识码:AReproduction of short-crested irregular waves bya viscous numerical wave tankFENG Guang1,WU Cheng-sheng1,ZHENG Wen-tao1,GU Min1,ZHU De-xiang2(1China Ship Scientific Research Center,Wuxi214082,China;2Shanghai Insitute of Computing Technology,Shanghai200040,China)Abstract:The short-crested irregular waves are reproduced by using a viscous numerical wave tank simu-lation technique.First,the simulation method for long-crested irregular waves is improved in this paper.The improved method is used for the simulating of long-crested irregular waves for different sea states.Finally, the results of numerical simulation are compared with the targets and show quite good agreement with them.Numerical simulation of short-crested irregular waves are carried out.The results also show good a-greement with the targets.Key words:short-crested irregular waves;long-crested irregular waves;numerical simulation;RANSE;VOF method1引言近年来,波浪及船舶运动相关研究领域正呈现出三个重要发展趋势[1]:从频域转向时域、从线性转向非线性、从势流转向粘性。
波浪水槽综合实验
波浪水槽综合实验一、实验目的:1、了解波浪水槽实验的基本原理和理论基础:包括基本造波方法、波浪理论、相似理论和近岸波浪传播现象2、了解造波机、浪高仪的基本构成和测量原理。
3、 通过实验采集一组波浪信号,分析波浪频谱特征4、 观测海堤附近波浪现象和越浪形态。
二、实验原理:1.造波方法和基本波浪理论自由表面重力波是船舶工程、海洋工程和海岸工程领域十分普遍的现象,配备造波机的波浪水槽是模拟波浪与二维结构物相互作用的常用实验设备。
通过给定造波信号由液压泵或步进电机控制推板运动,在波浪水槽中产生特定波列。
距离造波板2-3个波长外可以略去局部非传播模态的影响,可认为水槽中为行进波。
在水槽中通过浪高仪可以测量水槽中不同位置的波面时间过程线。
水槽中常用测力天平和压力传感器测量水动力载荷。
水槽末端设置多孔介质构成的消波区,消除反射波。
图1 波浪水槽示意图2.相似原理自由表面重力波的恢复力是重力,进行以重力为主要作用的流动实验通常采用重力相似准则或傅汝德数相似,其定义为/Fr v =,其中为流速,L 为特征长度,为重力加速度。
v g 波浪断面模型实验一般按重力相似准则设计。
若取几何比尺/2L p m L L 0λ==,有关物理模型比尺如下:时间比尺: 4.47t λ==速度比尺:4.47v λ==重度比尺: 38000WL λλ==单宽流量比尺:89.44Q λ==式中为工程原型长度,为模型长度。
pL mL 风速模拟通常按重力相似,风速测点位于测量断面上方中心。
3 近岸波浪现象3.1 线性波浪理论在平底均匀水深域中,根据势流理论波浪呈周期性分布。
单色行波波浪参数包括波浪周期T ,波长L ,波高H 和水深h ,如图2所示。
周期、波长和水深满足色散关系,对于线性波浪其表达式为,,其中波浪圆频率2tanh gk kh ω=2/T ωπ=,波数2/k L π=。
波高水深比为小量的波浪称为小振幅波,可用线性波浪理论描述,见图3。
第三章 波浪与波浪载荷
φ ( z , x, t ) =
三、线性波浪水质点运动特性
1 水质点速度
u=
∂φ kgH chk (d + z ) = cos(kx − ωt ) chkd ∂x 2ϖ ∂φ kgH shk (d + z ) w= = sin( kx − ωt ) chkd ∂z 2ϖ
2 加速度
ax ≈
∂u kgH chk (d + z ) = sin( kx − ωt ) chkd 2 ∂t kgH shk (d + z ) ∂w az ≈ =− cos(kx − ωt ) chkd ∂t 2
4 ∂t πD ∂u f I 2 ( x, z , t ) = ρ Cm 4 ∂t CM = 1 + Cm
s 2
f I 1 ( x, z , t ) = ∫∫ pds =ρ
πD 2 ∂u
Cm为附加质量系数(Added mass coefficient) 2 阻力 根据粘性流体力学得到:C D 为粘滞阻力系数。 园柱总载荷为:
H max ≈ 2.0 H 1 3 周期为10.7s。试求波浪为H max 时,总的水平力F H
(1)计算深水波长
L0 = L= gT 2 thkd 2π
gT 2 = 178.6(m) 2π
(2)利用微幅波理论查相关图表或利用公式 计算浅水波长L=174.29(m)
(3)计算桩柱相对直径 D/L=4.88/174.29=0.03<0.2,可以使用莫里森(Morrison)公式 (4)计算 FIm ax 及FD max , 求最大FH 出现的时刻 因为波高相对于水深较小,积分范围可以从海底到静水面。
式中:C M --惯性力系数 (Inertia Coefficient) C d --阻力系数(Drag Coefficient) u(x,z,t)为水质点速度。 ρ为水的密度,D为园柱直径。 讨论 1 惯性力部分: 惯性力被分成两部分, 一部分为假设园柱体不存在, 占据该体积的流体加速运动所需推力; 第二部分为跟随园柱运动的部分流体加速运动所需推力。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第36卷第4期海洋工程Vol. 36No.4 2018 年7 月THE OCEAN ENGINEERING Jul. 2018文章编号:1005-9865 (2018) 04-0078-10圆形数值波浪水池斜向和多向入射波浪模拟研究任兴月、陶军1,彭伟2(1.海南大学海洋工程与技术系,海南海口570228; 2.河海大学港口海岸与近海工程系,江苏南京210098)摘要:为了研究斜向人射波浪,基于三维不可压缩两相流模型,开发了一套圆形数值波浪水池数值模型。
在圆形波浪水池 中,通过源项造波法成功生成了任意人射方向的波浪,并且利用人工摩擦项模拟阻尼区以数值耗散反射波浪。
模型基于嵌人 式多块网格体系,米用F V M法(finite volume method)离散Navier-Stokes方程,V O F法(volume of flu id)追踪自由水面。
试验结 果表明,斜向人射波浪的模拟结果与理论值基本一致,圆形波浪水池在模拟斜向人射波浪时,有效区域的面积较传统波浪水 池显著增大,而且有效区域受波浪人射角度的影响也较小。
同时,通过叠加多列斜向人射波浪,模拟出了多向交叉波列,并通过与理论结果对比,发现其具有较髙的精度。
关键词:斜向人射波浪;圆形波浪水池;柱坐标系;V O F法;F V M法中图分类号:T V139.2文献标志码:A D O I:10.16483/j.issn.1005-9865.2018.04.009Study on oblique and multi-directional waves based on numerical circular wave tankR E N X in y u e1,T A O J u n1,P E N G W e i2(1.School of Civil Engineering and Architecture,Hainan University,Haikou570228, China; 2.College of Port,Coastal and Offshore Engineering,Hohai University,Nanjing210098, China)Abstract :To study the directional waves,a numerical circular wave tank was developed based on the three-dimensional incompressible two-phase flow model.W ithin the circular wave tank,waves with arbitrary incident direction could be generated by adopting a wave source term to the continuity equation,and the reflection waves could be numerically dissipated using a damping zone by adding an artificial friction ter^n to the momentum equation.The numerical computation was solved based on an embedded grid system,and the control equations were discreted by using the finite volume method (FVM).To track the free surface,the volume of fluid (VOF)was applied.The numerical results with the oblique waves agree well with the theoretical data.The results show that the circular wave tank has an advantage in simulating direction waves,and the effective area in the circular wave tank for oblique wave simulation is significantly larger than that in a traditional rectangular wave tank.And the incident angle of the oblique wave w ill not affect the effective area in the circular wave tank,while it is mainly influenced by the wave incident direction,wave maker in a rectangular wave tank.Moreover,multi-directional cross waves were also simulated by a superimposition of the oblique waves in this work with high accuracy.Keywords :oblique wave;circular numerical wave tank;cylindrical coordinate;VOF;FVM海洋工程结构物长期受到风浪流等海洋环境条件的影响,其中波浪是影响海洋工程结构物最重要的环 境因素。
波向是海洋波浪的另一个重要特性,其在海洋和海洋结构水动力荷载的应用中扮演很重要的角色。
收稿日期:2017-12-19基金项目:国家自然科学基金(51609056)作者简介:任兴月(1983-),男,博士,讲师,主要从事计算流体动力学研究。
通信作者:陶军,男,主要从事港口、海岸近海工程研究。
E-m ail:reiixiiy@第4期任兴月,等:圆形数值波浪水池斜向和多向入射波浪模拟研究79为了研究多向不规则波浪,20世纪50年代,法国的S o g re a h 实验室研制了分段式造波机,并将之安装于一个 矩形的平面水池中,利用各单块板之间运动的相位差成功制造了斜向入射波浪。
20世纪70年代,英国爱丁 堡大学利用分段式造波机进一步产生了多向不规则波浪。
李玉成等[1]利用海洋工程水池研究了斜向不规 则波浪与直立墙的作用;俞聿修等[2]研究了单向与多向不规则波浪作用下斜坡式建筑物护面块体的稳 定性。
由于分段式造波机产生的波浪是由各单块板产生波浪叠加的结果,每块板产生的波浪向外耗散传播,具 有三维性,同时相邻板之间缝隙间的水体流动,也会造成产生波浪的横向大小不均匀。
另外,制造多向波浪 或斜向入射波浪时,波浪入射方向与造波机和水池边壁都呈一定的夹角,生成的波浪遇到水池边壁会发生反 射,从而影响试验效果。
柳淑学[3]分析了分段式造波机的性能参数,发现造波机板总宽与单元板宽是影响 波浪的重要因素。
李木国等[4]也研究了分段式造波机的总宽度对波浪均匀性的影响,并证实了波浪在水池 两侧的反射等误差源会增加波浪能量的不均匀性。
俞聿修和柳学淑[5]进一步给出了多向不规则波浪水池 试验存在的问题,认为模型所能产生的波浪的方向范围受到分段式造波机性能的限制,而且波浪的有效区域 只是分段式造波机一侧一个较小的范围。
W illia m s 和C m ll [6]探讨了斜向入射波浪有效区的位置和形状,指 出其为由波浪的入射方向、分段式造波机与波浪外边界构成的一个三角形区域,并且研究了不同入射方向斜 向入射波的有效区,发现随着波向角的增加,水池边壁的反射波对实验的影响也逐步增加,实验有效区的面 积会逐步减小。
可见,在多向造波系统中,造波机板和波浪入射方向的夹角是影响波浪入射质量的重要因 素。
所以,多向造波系统的性能对波浪入射角的依赖性较大,波向角越大,多段式造波机的造波性能越差,水 池边壁反射波的影响区域越大。
为了改善矩形波浪水池的多向造波系统,研究者进行了大量的研究与尝试。
柳淑学[7]在产生斜向波和 多向不规则波时,考虑了试验水池边壁的反射波;李俊等[8—9]、冯光等[10]在数值波浪水池中,采用双边多段式 造波机,成功地增大了实验有效区的范围。
N a ito [11^l j 建造了一个新型的圆形波浪水池,并且给出了斜向入 射波浪的生成方法。
不同于传统的矩形波浪水池,其分段式造波系统沿水池周向呈圆形布置,基于可吸收反 射波浪的造波机,可以产生任意方向的波浪。
在这个圆形水池中,圆形多段式造波机沿径向运动产生柱面 波,通过叠加多列柱面波来生成平面波,相较于传统的矩形波浪水池,其对方向波浪具有优越的适应性。
然而针对圆形造波水池的研究只是刚刚起步,造波理论还没有健全,一直没有解决多向不规则波浪的生成问题。
N a tio [11]的研究基于柱面波叠加合成平面波,采用贝塞尔函数控制多段式造波机生成斜向入射波浪,如果基于该方法生成多向不规则波浪,造波机的控制函数会非常复杂。
T a n a k a 等[12]在水池内部采用弧形多段式造波机模拟斜向入射波浪,可以通过叠加多个弧形造波机的控制函数来生成多向不规则波浪,这为圆形波浪水池多向不规则波浪的生成理论奠定了基础。
为了研究圆形水池,Ren等[13-14]在T a n a k a 等[12]理论的基础上,建立了圆形数值波浪水池的数值计算模型(如图1所示),并成功制作了任意方向的斜向入射波浪。
在R e n 等[13-14]模型的基础上,通过修改控制方 程,在计算区域的内部设置了一个造波区,并施加了一个质量源项(B ro rs e n 和L a rse n [15];L i n 和L iu [16];K a w a sa ki [17])来生成任意方向的波浪,试验建筑物的反射波浪可以透过造波机传出计算区域而不会被再次 反射。
并且,采用人工黏度法(董志等[18];韩朋等[19])消除外围波浪,进一步完善了圆形数值波浪水池数值 模型。
模型采用嵌入式多块网格技术产生计算网格,并在计算网格上,采用有限体积法离散N a v ie r-S to k e s 方 程,通过V 0F 法追踪自由水面。
模拟试验表明圆形波浪水池中的有效区域面积较传统波浪水池显著增大, 而且有效区域受波浪入射角度的影响极小。
\造波区图1圆形波浪水槽示意 Fig . 1 Schematic view of circular wave tank80海 洋工程第36卷1控制方程数值波浪水池是一种利用C F D (计算流体力学)方法,在有限计算区域中,施加造波与消波条件,求解控 制方程,时域的数值模拟自由表面重力波浪的方法。