首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 171 毫秒
1.
宁德志  苏晓杰  滕斌 《海洋学报》2015,37(3):126-133
针对波浪与带有窄缝多箱体结构作用产生的流体共振问题,建立了基于域内源造波技术的二维非线性时域数值波浪水槽模型,其中自由水面满足完全非线性运动学和动力学边界条件,窄缝内流体引入人工阻尼来等效由于涡旋运动和流动分离引起的黏性耗散,计算域边界采用高阶边界元进行离散。通过模拟三箱体间两窄缝内相对波高变化,并与已发表的数值与实验结果对比,验证了本模型的准确性。同时通过大量的数值计算,分析了箱体数量对窄缝内水体共振频率、共振波高以及对结构反射波高和透射波高的影响。  相似文献   

2.
柏威  滕斌 《海洋工程》2001,19(3):43-50
采用二阶时域理论对非线性波浪在任意三维物体周围的绕射问题进行了研究,对自由表面边界条件进行Taylor级数展开,应用摄动展开可以建立相应的边值问题,而且此边值问题的计算域不随时间变化,运用基于B-样条的边界元方法求解每一时刻的波浪场,二阶自由表面边界条件在时间上进行数值积分,在自由表面加了一个人工阻尼层以避免波浪的反射,速度势分解为已知的入射势和未知的散射势,初始条件采用二阶Stokes波浪场,通过加入物体表面边界条件,得到散射势在时间和空间上的发展,本文对圆柱所受规则波的二阶波浪力和波浪爬高进行了计算,数值结果表明此理论计算准确,效率高,数值稳定。  相似文献   

3.
基于势流理论和时域高阶边界元方法,建立了三维完全非线性数值波浪水槽模型.利用源造波法产生入射波浪,应用五阶斯托克斯波理论给定波浪速度;采用混合欧拉-拉格朗日方法追踪流体瞬时水面,将二阶泰勒级数展开法应用于更新下一时间步的波面和速度势;通过加速势的方法准确计算自由水面速度的法向导数和物面速度势的时间导数.对完全非线性波浪进行了模拟,得到了稳定的波形.当波浪非线性较小时,与四阶Runge- Kutta法(RK4)计算结果和五阶斯托克斯波理论解均吻合良好;随着波浪非线性的增大,计算结果误差逐渐增大.通过数值试验分析,在满足精度要求的基础上,本方法计算时间略大于四阶Runge- Kutta法的四分之一,大大减少了计算量.  相似文献   

4.
灾害性波浪是中国沿海地区最具破坏性的自然灾害之一。采用开源程序OpenFOAM中interFoam求解器,对低顶海堤(在风暴潮和海平面上升情况下所面临的不利工况)的孤立波越浪特性开展数值模拟研究。通过孤立波冲击海堤的基准算例,验证模型在模拟波浪爬升和越浪过程中大变形波面以及剧烈波浪力方面的精度。基于验证的数值模型,对孤立波在低顶海堤上的越浪特征以及防浪墙高度对越浪的影响开展参数化研究。结果表明堤顶超高减小导致更为剧烈的越浪。针对尚无低顶海堤孤立波越浪量经验公式的问题,提出新的适用于堤顶超高小或为0的孤立波越浪量经验公式。此外,研究发现增加防浪墙高度可有效减少越浪,但防浪墙所受的波浪力也增大。综合考虑防浪墙减少越浪以及自身所受波浪力,针对文中研究采用的海堤截面和波浪条件,建议无量纲防浪墙高度取为1.00。  相似文献   

5.
波浪方向对两矩形驳船间窄缝内流体共振影响的数值研究   总被引:2,自引:1,他引:1  
应用二阶非线性三维势流时域模型模拟两个并靠矩形驳船间的窄缝内流体共振现象。模型中,速度势分解为入射势和未知的散射势,其中散射势通过求解边界积分方程计算得到。采用四阶预报校正的方法对自由水面边界条件在时间上进行数值积分。本文研究了波浪方向对窄缝内一阶、二阶波高的影响。计算结果表明,入射波浪的角度不会影响窄缝内流体的共振频率,同时窄缝内流体发生共振时,最大的共振波高总是发生在窄缝中点处。然而,共振时对应的最大波高随着入射角度的变化而变化。同时,入射波浪的频率与窄缝内流体共振频率的相对大小导致窄缝内最大波高位置向上游或者下游移动。  相似文献   

6.
王鹏  孙大鹏  吴浩 《海洋工程》2011,29(4):97-102
基于线性波理论,采用VOF法、RNG k-ε的湍流模型,通过FLUENT软件平台的二次开发功能实现了质量源造波。利用多孔介质模型,在动量方程里添加多孔介质源项,对斜坡堤上铺设异型块体时的爬高与越浪进行了数值模拟,并通过模拟爬高值率定出FLUENT中多孔阻力系数与多种异形块体糙渗系数之间的关系,进而建立了可供实际应用的海堤上波浪爬高与越浪数值模型。  相似文献   

7.
刘洪杰  刘勇  李玉成 《海洋学报》2009,31(4):159-166
基于线性势流理论,对斜向波与消浪室内带部分填料的开孔板式防波堤结构的相互作用进行了理论研究.利用匹配特征函数展开法给出了开孔防波堤结构反射系数和波浪力系数的理论计算方法,在极限情况下的计算结果与已有结果一致.利用数值算例分析了开孔结构反射系数和波浪力系数的主要影响因素.增加消浪室内填料厚度通常可以降低开孔墙结构的波浪力系数,加大结构的反射系数.随着波浪入射角度的增加,带填料防波堤结构的反射系数值先减小,达到最小值后又单调增加.  相似文献   

8.
三维完全非线性波浪水槽的数值模拟   总被引:7,自引:0,他引:7  
用有限元求解拉普拉斯方程,建立了三维完全非线性数值波浪水槽.跟踪流体自由表面的方法为满足完全非线性自由表面条件的半拉格朗日法,对离散单元采用20节点的六面体二次等参数单元.并把数值计算结果与水面初始升高产生箱体内流体运动解析解和二阶斯托克斯波理论解进行了对比,结果表明该模型是稳定的、守恒的,能精确模拟非线性波浪的产生和传播.  相似文献   

9.
在有、无礁冠两种典型岛礁剖面地形上开展孤立波冲击礁坪上直墙的二维水槽试验和数值模拟研究。通过OpenFOAM结合k-ε湍流模型进行RANS数值模拟,基于试验数据和验证后的RANS模拟结果,研究了岛礁剖面上孤立波传播的水动力特性及其对直墙的作用,并探讨了礁冠的存在对直墙所受冲击的影响。结果表明:RANS(k-ε)模型能够准确还原孤立波在岛礁地形上的主要传播特征,但在波浪破碎的情况下,RANS数值模拟不够精细,尤其对直墙所受瞬击荷载的捕捉能力欠缺;礁冠的存在,增强了波浪反射,使得波浪提前破碎,减小了直墙上的波浪爬升和最大动压。  相似文献   

10.
畸形波作用下JIP Spar平台波浪力分析   总被引:2,自引:1,他引:1  
利用时域高阶边界元法研究畸形波、系泊系统与浮式结构物耦合作用问题,运用四阶Runge-Kutta方法更新下一时间步的波面和速度势,采用非线性有限元法模拟系泊缆索的各种受力及缆索的弹性变形。以Longuet Hinggins模型为基础,调整组成波的初相位,在很短的波列中得到满足畸形波定义的波形。分析了聚焦位置、初始相位、谱峰周期、谱宽以及水流参数等对系泊JIP Spar平台缆绳力和波浪荷载的影响。将入射波浪的聚焦点分别设置在平台的迎浪侧和背浪侧,分析缆绳力和波浪荷载聚焦峰值,模拟结果发现其最大值产生在聚焦点为迎浪侧的情况。当聚焦点在迎浪侧时,缆绳力和波浪荷载聚焦峰值随着初相位的增加而减小,随着谱宽的增加而减小,而随着谱峰周期的增加而增加。  相似文献   

11.
波浪与带窄缝方箱作用共振现象的数值模拟   总被引:1,自引:1,他引:0  
By introducing a source term into the Laplace equation, a two-dimensional fully nonlinear time-domain numerical wave flume (NWF) is developed to investigate the resonance induced by the interaction bet...  相似文献   

12.
Gao  Jun-liang  Lyu  Jing  Wang  Jing-hua  Zhang  Jian  Liu  Qian  Zang  Jun  Zou  Tao 《中国海洋工程》2022,36(6):994-1006

In this paper, the transient fluid resonance phenomenon inside a narrow gap between two adjacent boxes excited by the incident focused waves with various spectral peak periods and focused wave amplitudes is simulated by utilizing the open-sourced computational fluid dynamics software, OpenFOAM. The weather-side box is allowed to heave freely under the action of waves, and the lee-side box keeps fixed. This paper mainly focuses on how both the spectral peak period and the focused wave amplitude affect the free-surface amplification inside the gap, the motion of the weather-side box, and the wave loads (including the vertical and the horizontal wave forces) acting on both boxes. For comparison, another two-box system with both boxes fixed is also considered as a control group. It is found that the motion of the weather-side box significantly changes the characteristics of the transient gap resonance, and it would cause that the fluid resonant period becomes 1.4–1.6 times that of the two-box system with both boxes fixed. All the concerned physical quantities (i.e., the free-surface amplification in the gap, the motion of the weather-side box, the wave loads) are found to closely depend on both the spectral peak period and the focused wave amplitude.

  相似文献   

13.
基于高阶边界元的三维数值波浪港池--波浪破碎的模拟   总被引:5,自引:1,他引:4  
在势流理论的框架内,采用高阶边界元方法和混合欧拉-拉格朗日法,实现了对三维波浪破碎过程的数值模拟.数值模型使用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解自由表面的演化过程.在所使用的边界元方法中,采用16节点三次滑移四边形单元来表示,这种单元在单元内具有高阶的精度同时在单元之间具有良好的连续性.给出了孤立波的传播和周期性非线性波浪沿缓坡传播的计算结果,表明数值模型具有良好的稳定性.  相似文献   

14.
X.T. Zhang  B.C. Khoo  J. Lou 《Ocean Engineering》2006,33(17-18):2310-2331
The problem of wave propagation in a fully nonlinear numerical wave tank is studied using desingularized boundary integral equation method coupled with mixed Eulerian–Lagrangian formulation. The present method is employed to solve the potential flow boundary value problem at each time step. The fourth-order predictor–corrector Adams–Bashforth–Moulton scheme is used for the time-stepping integration of the free surface boundary conditions. A damping layer near the end-wall of wave tank is added to absorb the outgoing waves with as little wave reflection back into the wave tank as possible. The saw-tooth instability is overcome via a five-point Chebyshev smoothing scheme. The model is applied to several wave propagations including solitary, irregular and random incident waves.  相似文献   

15.
The three-dimensional numerical model with σ-coordinate transformation in the vertical direction is applied to the simulation of surface water waves and wave-induced laminar boundary layers. Unlike most of the previous investigations that solved the simplified one-dimensional boundary layer equation of motion and neglected the interaction between boundary layer and outside flow, the present model solves the full Navier–Stokes equations (NSE) in the entire domain from bottom to free surface. A non-uniform mesh system is used in the vertical direction to resolve the thin boundary layer. Linear wave, Stokes wave, cnoidal wave and solitary wave are considered. The numerical results are compared to analytical solutions and available experimental data. The numerical results agree favorably to all of the experimental data. It is found that the analytical solutions are accurate for both linear wave and Stokes wave but inadequate for cnoidal wave or solitary wave. The possible reason is that the existing analytical solutions for cnoidal and solitary waves adopt the first-order approximation for free stream velocity and thus overestimate the near bottom velocity. Besides velocity, the present model also provides accurate results for wave-induced bed shear stress.  相似文献   

16.
极限波浪运动特性的非线性数值模拟   总被引:1,自引:0,他引:1       下载免费PDF全文
宁德志  滕斌  姜立明  臧军 《海洋学报》2008,30(3):126-132
利用时域高阶边界元方法建立了模拟极限波浪运动的完全非线性数值模型,其中自由水面满足完全非线性自由水面条件.采用半混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶Runge-Kutta方法更新下一时间步的波面和速度势,同时应用镜像格林函数消除水槽两个侧面和底面上的积分.研究中利用波浪聚焦的方法产生极限波浪,并且在水槽中开展了物理模型实验,将测点试验数据与数值结果进行了对比,两者吻合得很好.对极限波浪运动的非线性和流域内速度分布进行了研究.  相似文献   

17.
A fully nonlinear numerical model based on a time-domain higher-order boundary element method (HOBEM) is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here.  相似文献   

18.
A new coupling model of wave interaction with porous medium is established in which the wave field solver is based on the two dimensional Reynolds Averaged Navier-Stokes (RANS) equations with a closure. Incident waves, which could be linear waves, cnoidal waves or solitary waves, are produced by a piston-type wave maker in the computational domain and the free surface is traced through the Piecewise Linear Interface Construction-Volume of Fluid (PLIC-VOF) method. Nonlinear Forchheimer equations are adopted to calculate the flow field within the porous media. By introducing a velocity–pressure correction equation, the wave field and the porous flow field are highly and efficiently coupled. The two fields are solved simultaneously and no boundary condition is needed at the interface of the internal porous flow and the external wave. The newly developed numerical model is used to simulate wave interaction with porous seabed and the numerical results agree well with the experimental data. The additional numerical tests are also conducted to study the effects of seabed thickness, porosity and permeability coefficient on wave damping and the pore water pressure responses.  相似文献   

19.
The unsteady, two-dimensional Navier–Stokes equations and the exact free surface boundary conditions were solved to study the interaction of a solitary wave and a submerged dike. A piston-type wavemaker was set up in the computational domain to produce the incident solitary waves. The incident wave and the associated boundary layer flow in a wave tank with a flat bed were compared with the analytical solutions to verify the accuracy of this numerical scheme. Effects of the incident wave height and the size of the dike on the wave transformation, the flow fields, and the drag forces on the dike were discussed. Our numerical results showed that even though the induced local shear stress on the top surface of the dike is large at some particular locations, the resultant pressure drag is much larger than the friction drag. The primary vortex generated at the lee side of the dike and the secondary vortex at the right toe of the dike may scour the bottom and cause a severe problem for the dike.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号