动波浪壁圆柱绕流三维数值模拟

利用计算流体力学软件Fluent开展了三维动波浪壁圆柱绕流的数值模拟,建立了三维运动波浪壁圆柱模型,通过C语言自编程序实现波浪壁面的运动控制,并保证壁面变形时网格的高质量。在来流速度u=0.125 m/s、雷诺数Re=12 500的情况下,开展了动波浪壁波动速度w=0、0.062 5、0.125、0.187 5 m/s四个工况的计算分析,并比较了不同波动速度对流场结构、升力、阻力特性的影响。结果表明:动波浪壁圆柱能有效抑制流动的分离,消除交替脱落的尾涡,从而消除周期振荡的升力;在消除卡门涡街的同时,圆柱后驻点处的涡量值随波动速度增加而增加,其原因在于波形移动加大了壁面流体的速度,从而减小了圆柱前后的压力差,减小了阻力;随着波动速度的增大,平均阻力系数呈明显下降趋势,当波动速度为来流速度的1.5倍时,平均阻力系数相对于光滑圆柱下降了53.76%。     

寒潮影响下江苏沿海风浪场数值模拟研究

基于第三代浅水波浪数值预报模型SWAN,建立自西北太平洋嵌套至东中国海、江苏沿海的三重嵌套模型,对2010年12月12日至15日江苏沿海寒潮大风引起的风浪过程进行了数值模拟研究。利用西北太平洋和江苏沿海实测数据对模型进行了验证,结果表明SWAN嵌套模型能较好地模拟江苏沿海寒潮风浪场的时空分布。通过响水站实测数据对江苏沿海底摩擦系数进行了率定,研究表明选取Collins拖曳理论中摩擦因数C_f=0.001时,有效波高模拟误差相对较小。寒潮风浪场的特征分析表明,有效波高分布与风场分布基本一致,寒潮风浪在江苏沿海北部影响较为显著,辐射沙洲附近由于其特殊地形影响相对较小。     

Extreme wave crest distribution by direct numerical simulations of long-crested nonlinear wave fields

The wave crest height qualification checks are required during the wave calibration before the model test in wave basin. However, the reliable criteria of nonlinear wave crest probability distribution in 3-h duration (full-scale) has not been well established yet. We investigate wave crest-height statistics of long-crested nonlinear wave fields using high-order spectral (HOS) method, which can take the effects of both second-order bound waves and third-order free waves into account. The energy dissipation effects due to wave breaking were included by employing an eddy viscosity model. Sensitivity analyses to the wave breaking onset criterion have been performed. Validation is provided by comparing the obtained numerical results with the available calibration test data. Based on extensive and direct numerical simulations, semi-empirical single realization distributions for wave calibration have been developed through 3-parameter Weibull fitting and systematic regression analyses. Particular attention has been paid to the tail of upper bound of wave crest distributions. The effects of wave steepness and water depth on the maximum wave crest height in 3-h duration have been examined. It is found that with the increase of wave steepness, the extreme wave crest height increases until it reaches a critical value. In addition, for the scale water depth k_ph < 1.36, the maximum crest height decreases as the water depth increases, while in the opposite case the maximum crest height increases as the water depth increases. Moreover, it is confirmed that that the fourth-order nonlinearity does not have significant effects on the distribution of the wave crest height.     

槽道湍流的展向振荡电磁力壁面减阻

采用直接数值模拟方法,对槽道湍流的展向振荡电磁力的减阻效果和减阻机理进行了研究,讨论了电磁力强度和振荡频率对湍流猝发事件以及壁面减阻率的影响.结果表明,电磁力强度或振荡频率变化时,湍流猝发频率和猝发强度的变化趋势是相反的,所以存在最优参数使得减阻效果最好.等价壁面展向速度可以很好地描述电磁力强度和振荡频率的变化对减阻效果的综合效应.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

悬链线管体三维流场数值模拟研究

基于计算流体力学开源代码Open FOAM开展不同雷诺数(100≤Re≤3 900)下弯曲管体绕流流场三维数值模拟研究。对比分析不同雷诺数条件下弯曲管体的瞬时和时均流场特征,研究了弯曲管体不同位置处的时均压力系数和旋涡脱落频率。结果表明:当Re=100时,尾流中仅有展向涡存在;当300≤Re≤3 900时,流向涡出现且强度逐渐增加,展向涡减弱,随着局部曲率增大,沿弯管展向涡脱强度减小。随着雷诺数的增大,回流区沿流向的延伸尺度减小,但其沿展向的分布范围增大。沿弯管展向随着深度减小,前后驻点时均压力系数的绝对值递减,圆周时均压力系数变化幅度减小。在300≤Re≤3 900范围内,管体曲率会对展向涡脱频率产生一定影响。     

西沙湾工程波浪推算及波浪场数值模拟

本研究根据崇武水文站1976—2007年共31 a的波浪观测资料,使用皮尔逊Ⅲ型频率曲线进行拟合分析,推算出崇武在不同重现期条件下的特征波浪要素,再将其作为外海波浪要素输入,通过基于椭圆型缓坡方程的CGWAVE近岸波浪数值模型,模拟西沙湾的重现期波浪场,得到其设计波浪要素,得到较为合理的计算结果。西沙湾的海底地形强烈地影响近岸波浪,导致了波高和波向的不均匀分布特征,岛屿和礁石都是天然屏障,东侧水域受龟屿及闽台码头的掩护,波浪反射、绕射明显,对避风坞起到了有效的保护作用。此外,对于不同地质的边界,反射系数不同使得近岸波高分布有较大差异。     

Three-dimensional numerical simulation of oblique wave action on vertical walls

1 Introduction A vertical wall is one of the typical maritime sheltering structures in the coastal region, which is mainly subjected to wave forces. It has been a com- mon engineering assumption adopted for design pur- poses that normally incident wave fo…     

Three-dimensional numerical simulation of obfique wave action on vertical walls

A three-dimensional numerical model based on the potential theory was developed to study the oblique wave action on vertical walls. A source term inside the domain was used to generate incident waves and outgoing waves were dissipated by sponge layers and transmitted by radiation boundaries. The finite difference method was used to solve the governing equations and boundary conditions in the regular transformed domain in σ-coordinate. Satisfactory agreements between the numerical predictions and experimental results of wave force were obtained. It is concluded that the maximum wave force acting on the vertical walls is induced by the obliquely incident waves rather than the normally incident waves.     

孤立波爬坡的二维数值模拟

波浪在滩地上以及遇海岸工程后传播发生变形、爬坡等现象,对其进行数值模拟具有广阔的工程应用背景.应用基于Boltzmann方程的KFVS(kinetic flux vector splitting)格式求解二维浅水方程,同时采用干底Riemann解模拟动边界问题.模型模拟了孤立波在滩地上传播变形、爬坡的过程,以及孤立波在滩地上遇圆柱后绕射、变形和爬坡的过程.计算结果与实验结果非常吻合,表明模型具有较大的推广应用价值.     

网衣对波浪传播影响的三维数值模拟

为研究网衣对波浪传播的影响,采用多孔介质模型模拟网衣,建立了模拟网衣与波浪相互作用的三维数值波浪水槽模型。基于该数值模型研究了波浪经过网衣作用后波高衰减的变化规律,并与物理模型试验结果进行比较。通过对比数值模拟和物理模型试验的结果,证明了该三维数值模型模拟网衣对波浪传播影响的可行性。     

基于紊流随机理论的航槽三维流动数学模型

根据窦国仁的紊流随机理论，建立了一种模拟河口海岸水域中航槽三维流动的数学模型。采用控制体积法导出三维偏微分方程的离散格式；将水压力分解为动水压力和静水压力，用Patankar和Spalding提出的压力校正法求解动水压力，通过求解水位控制方程来得到自由表面；紊流模型采用安国仁提出的紊流随机理论，克服了k—ε模型中采用各向同性紊动粘滞系数的不足，而k—ε模型可作为紊流随机理论的一个特例。该模型计算了各种不同挖深比，各种航槽与水流交角的航槽流速分布，利用该模型计算得到的挖槽中的流速分布与水槽及水池试验资料相吻合。利用该模型可为开敞水域中开挖航槽的选择提供依据。     

聚焦波浪在浅堤上传播变形高精度数值模拟研究

基于开源程序REEF3D,通过建立高精度二维数值波浪水槽,系统研究了聚焦波浪在浅堤上传播变形的规律,着重分析了聚焦波浪通过浅堤的水动力过程及能量变化规律,讨论了不同波浪要素对聚焦波浪传播特性的影响。除此之外,还考虑了双浅堤布置对聚焦波浪传播变形的影响。研究结果表明:极端波浪通过浅堤时,堤顶水深越小,波浪主频能量衰减越显著。在给定堤顶水深条件下,聚焦点与浅堤的相对位置对聚焦波浪能量的衰减影响较小。在双浅堤布置条件下,随着浅堤间距的增加,上下游浅堤的相互影响逐渐减弱,高频段的波浪能量也随之减小。     

基于FLUENT的波浪溢流水动力数值模拟

波浪溢流现象使得海堤受到了越浪和溢流的联合作用,复杂的水动力过程会引起海堤后坡产生严重的侵蚀破坏。基于FLUENT软件建立了二维数值波浪溢流水槽模型,该模型运用UDF速度边界造波法分析在不同超高条件下海堤后坡流量和水流厚度的水力学特性。结果表明数值模拟结果与前人物理模型试验结果吻合,该模型可以真实地模拟出海堤波浪溢流现象。在此基础上进一步研究了波浪溢流中越浪和溢流在不同相对超高条件下的主导性作用,而后建立了十分准确的波浪溢流海堤后坡稳定水流厚度计算公式。     

Linear irregular wave generation in a numerical wave tank

In the design of any floating or fixed marine structure, it is vital to test models in order to understand the fluid/structure interaction involved. A relatively inexpensive method, compared to physical model testing, of achieving this is to numerically model the structure and the wave conditions in a numerical wave tank. In this paper, a methodology for accurately replicating measured ocean waves in a numerical model at full scale is detailed. A Fourier analysis of the measured record allows the wave to be defined as a summation of linear waves and, therefore, Airy's linear wave theory may be used to input the wave elevation and associated water particle velocities. Furthermore, a structure is introduced into the model to display the ability of the model to accurately predict wave–structure interaction. A case study of three individual measured waves, which are recorded at the Atlantic marine energy test site, off the west coast of Ireland, is also presented. The accuracy of the model to replicate the measured waves and perform wave–structure interaction is found to be very high. Additionally, the absolute water particle velocity profile below the wave from the numerical model is compared to a filtered analytical approximation of the measured wave at a number of time-steps and is in very good agreement.     

Effect of wave action of a two-layer flow on underwater obstacles

A force action caused by the generation of internal waves in the vicinity of a density jump layer in a marine medium and exerted on underwater obstacles streamlined by a two-layer flow has been studied as part of a theoretical model. A horizontal element of an engineering structure, for example, a transport pipeline, which is modeled by a point dipole, is selected as an underwater obstacle. The integral representations for the wave resistance and the lift force are obtained. The characteristic features of variation in the hydrodynamic reactions to the streamlined obstacle, as well as the conditions of their significant intensification are revealed.     

均匀流中直立圆柱体绕流三维数值模拟

研究直立贯底圆柱体的三维粘性绕流问题。以不可压缩Navier-Stokes方程为控制方程,采用有限体积法和SIMPLE算法,建立了数值模拟方法。考察在不同水平和垂直断面上,圆柱体绕流产生的尾涡和流动速度场的分布特性,成功地数值模拟了直立贯底圆柱体绕流场的三维特性。结果表明,在考虑重力影响的情况下,直立圆柱体周围的流动具有明显的三维特性,而且沿圆柱体轴向不同断面上的尾涡分布是不相同的。     

Determination of wave energy dissipation factor and numerical simulation of wave height in the surf zone

Based on the time-dependent mild slope equation including the effect of wave energy dissipation, an expression for the energy dissipation factor is derived in conjunction with the wave energy balance equation. The wave height of regular and irregular waves is numerically simulated by use of the parabolic mild slope equation considering the energy dissipation due to wave breaking. Comparison of numerical results with experimental data shows that the expression for the energy dissipation factor is reasonable. The effects of the wave breaking coefficient on the breaking point and the distribution of wave height after breaking are discussed through the study of a specific experimental topography.     

湍流数值模拟中封闭模式应用的局限性

通过对琼州海峡的潮流场特征进行数值模拟，指出了选择不同的特征混合长度表达式对数值模拟结果的影响，表明了基于特征混合长度理论的流封闭模式在近海湍流数值模拟中应用的局限性。