垂向变化的波浪辐射应力影响下流场数值研究

采用浪流沙耦合模型COHEREN-SED模型进行黄河三角洲的海流三维运动数值模拟。文中将垂向平均的波浪辐射应力与垂向变化的波浪辐射应力引入至COHERENS-SED,通过该模型进行波浪辐射应力影响下的海流场演变数值研究。研究结果显示波浪辐射应力的存在显著改变了水深小于5m的流场。垂向变化的辐射应力则使得其表层波生沿岸流速大于传统的垂向均匀辐射应力引起的波生沿岸流速,底层则刚好相反。当波浪强度很大时,流场中的波生沿岸流对水深较浅的海岸区海流场影响较大。     

基于二阶时间步进法数值波浪水槽研究

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

基于Boussinesq数值模型的波浪速度垂向分布模拟研究

在Liu和Fang推导的双层Boussinesq方程基础上,将其简化为一层水波方程,并建立了基于混合4阶Adams-Bashforth-Moulton时间格式的立面二维数值模型;数值模拟了波浪在潜堤上的演化过程,并将数值计算结果与相关实验结果进行了对比,验证了该数值模型的正确性。进而在不同的入射波条件下,将沿着水深分布的水平速度和垂向速度的数值模拟结果与线性、二阶、三阶解析解解析结果进行综合对比。对比结果表明,在不同的无因次水深kh条件下,数值解与解析解的整体吻合程度较好。在不同的波陡H/L条件下,数值解展现了较好的非线性特征。在不同的波高水深比H/h条件下,数值解与解析解之间的整体差异较小。可以看到,该数值模型较好地模拟了波浪垂向速度场分布,体现了其优良的综合性能。     

辐射应力对台风风暴潮预报的影响和数值研究

台风过程期间,风暴潮和海浪是相伴相生的,相互作用的.波致辐射应力对于近岸风暴增、减水起着十分重要的作用,传统的海浪模式计算辐射应力耗时较多,不能满足业务化预报的要求.根据已有波浪辐射应力的理论表达式,经过严密的数学推导,适当的简化处理,提出了一个较为简单的波浪辐射应力表达式,并将其应用到业务化风暴潮数值预报模式中去,通...     

浮式防波堤性能理论计算

本文论述了在任意浪向微幅入射波作用下,浮式防波堤性能理论计算。在求解其绕射势和辐射势边值问题后,引入无穷远处波浪速度势渐近表达式,应用波浪能量守恒原理,分别计算出由浮式防波堤横荡、垂荡和横摇运动引起的堤后波和当其被固定时在入射波作用下在堤后所兴起的波浪。最后计算出表达浮式防波堤防波效果的,以堤后波幅值与入射波幅值之比所定义的堤后波系数G_T。 本文对二个计算实例进行了理论计算,结果与试验相比吻合良好。     

柱对称双浮子波能装置散射问题的解析方法

通过构造绕射势特解的一种新表达式,基于特征函数展开法提出了线性入射波作用下柱对称双振荡浮子波能装置散射问题的一种新的表达式,并在此基础上计算了垂向波浪激励力.从速度势的求解过程可知,绕射势和辐射势的计算量相当,但垂向波浪激励力的计算比采用Haskind关系方法简单得多,两种方法的结果一样.另外对不同工况下垂向波浪激励力的变化趋势进行了分析.     

波浪与水流共同作用下的流速场

本文讨论了波浪与稳定流相互作用的二元流条件下的流速场问题。分析表明由于波浪与水流相互作用的结果水流的垂线流速分布变得更为均匀,这一结果已经实验验证。结果还表明对于波流共同作用下的水平流速场可以应用波浪与水流二者水平流速值的叠加原理。数值计算的结果表明,对于波面高程和水平流速值,由非线性波理论所得的结果较线性波理论的结果为好。作者认为,在波流共同作用条件下,在工程实用上为了计算水平流速场,当相对水深 d/L 2大于0.1时,司采用斯托克斯三阶波或五阶波理论,当相对水深 d/L_2小于0.1时,司采用椭圆余弦波理论(d 为水深,L_2为静水中的波长)。     

浅水浮式风机基础水动力特性及波浪荷载的CFD数值分析

采用计算流体动力学(CFD)方法对OC4-DeepCwind半潜式风机基础在30～100 m水深下的水动力特性和非线性波浪荷载进行了研究。基于验证的数值波浪水池,通过结构自由衰减运动及规则波下浮式基础的幅值响应算子(RAO)计算,分析了浮式风机基础的结构动力特性及水动力特性随水深的变化规律,研究了浅水条件下非线性波浪荷载特征。研究结果表明,随着水深的减小,风机基础结构垂荡固有周期略有增大,垂荡和纵摇的阻尼比均呈现逐渐增大的趋势;与黏性流理论计算结果相比,势流结果在浅水波浪工况下对垂荡和纵摇RAO存在显著低估;在浅水波浪作用下,风机基础所受波浪荷载的三倍波频成分占比可超过5%,体现了浅水条件下波浪与结构之间较强的非线性作用。     

规则波作用下四锚浮标系统动力分析

随着水深的增加,海洋浮标系泊缆索的长度也增加,对其浮标系泊系统进行动力分析是非常重要的.基于势流理论,应用边界元方法在频域内计算浮标的附加质量和阻尼、一阶波浪力和二阶平均漂移力,通过快速傅里叶变换将计算结果转换到时域,应用莫里森公式计算浮标及其skirt的粘性阻尼力,应用非线性有限元方法计算系泊缆索张力,最后在时域内应用四阶Runge- Kutta法计算浮标的运动响应和系泊缆索张力.比较分析了两种水深条件下,不同skirt直径和缆索长度情况下浮标的水动力特性、运动响应和系泊缆索张力.计算结果表明,浮标skirt的存在对纵荡方向的附加质量和辐射阻尼基本没有影响,而垂荡方向的附加质量增加,辐射阻尼减小,且能够减小二阶平均波浪力的峰值.     

浅水区波浪非线性效应对砂质海床动力响应的影响

以广义Biot动力固结理论为基础,运用一阶椭圆余弦波和二阶Stokes波等非线性波浪理论考虑浅水区波浪荷载的非线性效应,在时域上采用有限元方法对非线性波浪力作用下饱和砂质海床的动力响应进行了数值求解,并与线性波浪作用下海床动力响应特性进行了对比分析。结果表明,随着波长与水深之比L／d及无量纲参数T(g/d)^1/2的增大,非线性波浪对海床动力响应的影响增大。与线性波浪理论相比,孔隙水压力与有效应力幅值的增大效应非常显著。因此在近海海洋建筑物设计与工程场地评价中,波浪力的非线性特性必须引起注意。     

Nonlinear wave interaction with a vertical circular cylinder. Part II: Wave run-up

Theoretical results for second-order wave run-up around a large diameter vertical circular cylinder are compared to results of 22 laboratory experiments conducted in regular nonlinear waves. In general, the second-order theory explains a significant portion of the nonlinear wave run-up distribution measured at all angles around the cylinder. At the front of the cylinder, for example, measured maximum run-up exceeds linear theory by 44% on average but exceeds the nonlinear theory by only 11% on average. In some cases, both measured run-up and the second-order theory exceed the linear prediction by more than 50%. Similar results are found at the rear of the cylinder where the second-order theory predicts a large increase in wave amplitude for cases where the linear diffraction theory predicts little or no increase. Overall, the nonlinear diffraction theory is found to be valid for the same relative depth and wave steepness conditions applicable to Stokes second-order plane-wave theory. In the last section of the paper, design curves are presented for estimating the maximum second-order wave run-up for a wide range of conditions in terms of the relative depth, relative cylinder size, and wave steepness.     

Establishment of Numerical Wave Flume Based on the Second-Order Wave-Maker Theory

With growing computational power, the first-order wave-maker theory has become well established and is widely used for numerical wave flumes. However, existing numerical models based on the first-order wave-maker theory lose accuracy as nonlinear effects become prominent. Because spurious harmonic waves and primary waves have different propagation velocities, waves simulated by using the first-order wave-maker theory have an unstable wave profile. In this paper, a numerical wave flume with a piston-type wave-maker based on the second-order wave-maker theory has been established. Dynamic mesh technique was developed. The boundary treatment for irregular wave simulation was specially dealt with. Comparisons of the free-surface elevations using the first-order and second-order wave-maker theory prove that second-order wave-maker theory can generate stable wave profiles in both the spatial and time domains. Harmonic analysis and spectral analysis were used to prove the superiority of the second-order wave-maker theory from other two aspects. To simulate irregular waves, the numerical flume was improved to solve the problem of the water depth variation due to low-frequency motion of the wave board. In summary, the new numerical flume using the second-order wave-maker theory can guarantee the accuracy of waves by adding an extra motion of the wave board. The boundary treatment method can provide a reference for the improvement of nonlinear numerical flume.     

Second-order theory for coupling 2D numerical and physical wave tanks: Derivation,evaluation and experimental validation

A full second-order theory for coupling numerical and physical wave tanks is presented. The ad hoc unified wave generation approach developed by Zhang et al. [Zhang, H., Schäffer, H.A., Jakobsen, K.P., 2007. Deterministic combination of numerical and physical coastal wave models. Coast. Eng. 54, 171–186] is extended to include the second-order dispersive correction. The new formulation is presented in a unified form that includes both progressive and evanescent modes and covers wavemaker configurations of the piston- and flap-type. The second order paddle stroke correction allows for improved nonlinear wave generation in the physical wave tank based on target numerical solutions. The performance and efficiency of the new model is first evaluated theoretically based on second order Stokes waves. Due to the complexity of the problem, the proposed method has been truncated at 2D and the treatment of regular waves, and the re-reflection control on the wave paddle is also not included. In order to validate the solution methodology further, a series of nonlinear, periodic waves based on stream function theory are generated in a physical wave tank using a piston-type wavemaker. These experiments show that the new second-order coupling theory provides an improvement in the quality of nonlinear wave generation when compared to existing techniques.     

Applicability of numerical models to nonlinear dispersive waves

The applicability of three different wave-propagation models in nonlinear dispersive wave fields has been investigated. The numerical models tested here are based on three different wave theories: a fully nonlinear potential theory, a Stokes second-order theory, and a Boussinesq-type theory with an improved dispersion relation. Physical experiments and computations were conducted for wave evolutions during passage over a submerged shelf under various wave conditions. As expected, the fully nonlinear solutions agree better with the measurements than do the other solutions. Although the second-order solution has sufficient accuracy for smaller-amplitude wave cases, the truncation after the third harmonics causes significant discrepancies in wave form for larger waves. In addition, the second-order model markedly overestimates the first- and second-harmonic amplitudes in transmitted waves. The Boussinesq model provides excellent predictions of wave profile over the shelf even in larger wave cases. However, this model also overestimates the magnitudes of several higher harmonics in transmitted waves. These facts may indicate that energy transfer from bound components into free waves in these higher harmonics cannot be accurately evaluated by the Boussinesq-type equations.     

Nonlinear wave interaction with a vertical circular cylinder: wave forces

Second-order wave forces on a large diameter vertical circular cylinder, computed according to a semi-analytic nonlinear diffraction theory, are compared to results of 22 laboratory experiments with regular waves. In general, predicted forces agree quite well with measured forces. In most tests, both measured and predicted maximum forces exceeded linear theory by 5 to 15%. In a few cases, however, the measured forces were less than those predicted by linear theory, in contrast to the second-order predictions. It is shown that these results are related to the phasing of various linear and nonlinear wave force components, and are consistent with those obtained by other investigators.     

Spectral analysis of Stokes waves

The spectral properties of Stokes waves are shown in this paper by theoretical and numerical methods. This is done by expressing wave profiles and velocities of water particles as nonlinear combinations of the first order component of wave profiles. Under the assumption of the first order wave profiles being zero mean Gaussian processes, the relationship between autocorrelation functions of wave profiles and velocities of water particles and the first order component of wave profiles is established using the nonlinear spectral analysis. The spectral densities of nonlinear random waves, the velocities and accelerations of water particles are then obtained. Numerical computations are carried out to analyze the effect of fundamental parameters of waves. The results indicate that wave height is the most sensible parameter to the root mean squares related and wave depth is the least sensible one of all.     

Spectral Analysis of Nonlinear Drag Forces

The spectral properties of nonlinear drag forces of random waves on vertical circular cylinders are analyzed in this paper by means of nonlinear spectral analysis. The analysis provides basic parameters for estimation of the characteristic drag forces. Numerical computation is also performed for the investigation of the effects of nonlinearity of the drag forces.The results indicate that the wave drag forces calculated by linear wave theory are larger than those calculated by the third order Stokes wave theory for given waves. The difference between them increases with wave height. The wave drag forces calculated by use of hnear approximation are about 5% smaller than their actual values when measured in the peak values of spectral densities. This will result in a safety problem for the design of offshore structures. Therefore, the nonlinear effect of wave drag forces should be taken into comidemtion in design and application of important offshore structures.     

A Three-Dimensional Nearshore Hydrodynamic Model with Depth-Dependent Radiation Stresses

For the simulation of the three-dimensional (3D) nearshore circulation,a 3D hydrodynamic model is developed by taking into account the depth-dependent radiation stresses.Expressions for depth-depeedent radiation stresses in the Cartsian coordinates are introduced on the basis of the linear wave theory,and then vertical variations of depth-dependent radiation stresses are discussed.The 3D hydrodynamic model of EICIRC (Eulerian-Lagrangian CIRCulation) is extended by adding the terms of the depth-dependent or depth-averaged radiation stresses in the momentum equations.The wave set-up,set-down and undertow are simulated by the extended ELCIRC model based on the wave fields provided by the experiment or the REF/DIF wave model.The simulated results with the depth-dependent and depth-averaged radiation stresses both show good agreement with the experimental data for wave set-up and set-down.The undertow profiles predicted by the model with the depth-dependent radiation stresses are also consistent with the experimental results,while the model with the depth-averaged radiation stresses can not reflect the vertical distribution of undertow.     

Hydroelasticity of a floating plate in multidirectional waves

The membrane forces are included in the hydroelastic analysis of a floating plate undergoing large vertical deflections in regular monochromatic multidirectional waves. The first-order vertical displacements induced by the linear wave exciting forces are calculated by the mode expansion method in the frequency domain. The second-order vertical displacements induced by the membrane forces are calculated by the von Karman plate theory. The results show that the membrane contribution both in terms of the axial stresses and the effect on the bending stresses can be important.     

Nonlinear wave profiles of wave–wave interaction in a finite water depth by fixed point approach

In this paper, a superposition of two periodic wave profiles in a finite water depth was investigated. This paper is focused on the improvement of a wave profile on the linear superposition of two waves. This improvement was realized by introducing an iterative method, which was based on a fixed point approach. Application of the fixed point approach to the wave superposition made it possible to obtain a wave profile of wave–wave interaction. The improved result of the wave profile was in good agreement with that of the nonlinear perturbation solution of the second order. It was interesting that the improved result revealed the higher-order nonlinear frequencies for two interacting Stokes waves while Dalzell's solution by a perturbation method could not predict them.