水下平板与波浪相互作用完全非线性数值模拟(英文)

利用完全非线性数值波浪水槽技术研究水下平板与波浪的相互作用。假定水下平板厚度极薄、刚性,位于有限水深并且非常接近自由水面。应用四阶龙格库塔方法追踪每一时刻的波面形状,采用阻尼层来吸收反射波以保证算法的稳定性,同时引入平滑和重组的方法抑制自由表面控制点的较高梯度。通过对波浪与浮动圆柱相互作用的数值模拟证实了数值波浪水槽方法的有效性,计算结果与线性理论吻合良好。在波浪数值水槽方法中引入造波板模拟波浪产生并与水下平板发生相互作用,应用傅立叶解析方法对波面变形、波浪力作了分析。结果表明在板非常接近自由水面的情况下会表现出现很强的非线性,揭示了线性理论的局限性。     

Numerical simulation of sloshing waves in a 3D tank based on a finite element method

The sloshing waves in a three dimensional (3D) tank are analysed using a finite element method based on the fully non-linear wave potential theory. When the tank is undergoing two dimensional (2D) motion, the calculated results are found to be in very good agreement with other published data. Extensive calculation has been made for the tank in 3D motion. As in 2D motion, in addition to normal standing waves, travelling waves and bores are also observed. It is found that high pressures occur in various circumstances, which could have important implications for many engineering designs.     

Nonlinear response of membranes to ocean waves using boundary and finite elements

The behavior of a highly deformable membrane to ocean waves was studied by coupling a nonlinear boundary element model of the fluid domain to a nonlinear finite element model of the membrane. The hydrodynamic loadings induced by water waves are computed assuming large body hydrodynamics and ideal fluid flow and then solving the transient diffraction/radiation problem. Either linear waves or finite amplitude waves can be assumed in the model and thus the nonlinear kinematic and dynamic free surface boundary conditions are solved iteratively. The nonlinear nature of the boundary condition requires a time domain solution. To implicitly include time in the governing field equation, Volterra's method was used. The approach is the same as the typical boundary element method for a fluid domain where the governing field equation is the starting point. The difference is that in Volterra's method the time derivative of the governing field equation becomes the starting point.The boundary element model was then coupled through an iterative process to a finite element model of membrane structures. The coupled model predicts the nonlinear interaction of nonlinear water waves with highly deformable bodies. To verify the coupled model a large scale test was conducted in the OH Hinsdale wave Research Laboratory at Oregon State University on a 3-ft-diameter fabric cylinder submerged in the wave tank. The model data verified the numerical prediction of the structure displacements and of the changes in the wave field.The boundary element model is an ideal modeling technique for modeling the fluid domain when the governing field equations is the Laplace equation. In this case the nonlinear boundary element model was coupled with a finite element model of membrane structures, but the model could have been coupled with other finite element models of more rigid structures, such as a pontoon floating breakwater.     

Wavelet analysis for processing of ocean surface wave records

Wavelet analysis is a relatively new technique and in the recent years enormous interest in application of wavelets has been observed. This modern technique is particularly suitable for non-stationary processes as in contrast to the Fourier transform, (FT), the wavelet transform (WT) allows exceptional localization, both in time and frequency domains. The wavelet transform has been successfully implemented in signal and image processing, ordinary and partial differential equation theory, numerical analysis, communication theory and other fields. On the other hand, the application of the WT to ocean engineering and oceanography is rare. In this paper the WTs capability to give a full time–frequency representation of the wave signals is demonstrated. The processing of the time series of the non-stationary deep water waves, waves breaking at the tropical coral reefs and mechanically generated waves in the wave flume demonstrates the ability of the wavelet transform technique to detect a complex variability of these signals in the time–frequency domain. Various spectral representations resulting from the wavelet transform are discussed and their application for wave signals is shown.     

Measurement of ocean wave spectra using narrow-beam HE radar

A data interpretation algorithm is developed to extract ocean wave information from HF radar backscatter observed by a narrow-beam antenna system. The basis of this measurement is the inversion of the integral equation representing the second-order radar cross section of the ocean surface. This equation is numerically inverted by approximating it as a matrix equation and pseudoinverting the kernel matrix using a singular value decomposition. As a test of this algorithm, comparisons are made between wave spectrum estimates obtained from a WAVEC buoy and a pair of 25.4-MHz ground wave radars, using data collected during the 1986 Canadian Atlantic Storms Program (CASP). Overall, the results of this experiment have been positive and have demonstrated both the basic feasibility of the inversion algorithm and the wave sensing capability of HF radar. For example, significant wave height estimates deduced by two radars differed from the buoy, in an absolute value sense, by only 0.12 m on average. When using only one radar, the mean difference of this important parameter from the buoy was a reasonable 0.33 m     

A Calculation Method for Incidence and Reflection of Irregular Waves Along Seawalls

Some researches have been made in this aspect.In the method by Walton Jr.(1992),incidentwaves are supposed to be the overlapping result of M component waves with different frequencies whichmay take different directions,the direction of incident waves should be available in advance,but in fact thedirection of incident waves is not available.In our study,incident waves are supposed to be composed ofM overlapping component waves with different frequencies,and different frequencies have different direc-tions.Based on the irregular wave reflection theory,the calculation formulas of wave direction,complexamplitude of incident waves,and complex amplitude of reflected waves in surface which are composed ofcomponent waves are derived by means of discrete Fourier transform.Then,the frequency spectra of inci-dent waves and reflected waves and the reflection coefficient of waves with corresponding frequencies areobtained.Verification of the method and the calculation results from in-situ measured data indicate thatt     

Efficient computation method for two-dimensional nonlinear waves

I~IOWNUmerics wave tank is a ~ tOOl by which variouS nonlinear wave - ac interactionpwhleln can be treated in time domain. As the face stage, we f~ our efforts on the devious of an effita saution tee~ of fully nofor waves. aam numtried of water - wavesconnected PwhlemS, es~ in a wave tank, is alwayS comas Of PartS: wave generation, P~ and a~. The effi~ Of a nUm~ wave tankdepends not Only on the qwhty Of the wave abotion tecboaue, which allowS the lOng times~ and keePS finite tank leng…     

EMD方法和Hilbert谱分析法的应用与探讨

利用EMD方法对海浪观测资料进行处理 ,通过在数据两端的“平衡位置”处分别附加平行直线段的方法进行端点抑制 ,分解出 1 0个内在模函数和 1个剩余趋势项 ,再对各内在模函数进行Hilbert变换 ,得到波浪的Hilbert谱。对所得结果的分析表明 ,各模态在Hilbert谱中的分布趋势和Fourier谱中谱线的变化趋势是一致的 ,第一模态的中心频率与Fourier谱的谱峰频率相对应 ;EMD方法是对非线性、非平稳过程数据进行距平化的好方法 ,距平化的过程和消除趋势项的处理是统一的。     

Diffraction of a directionally spread wave group by a cylinder

The problem of diffraction of a directionally spread focused wave group by a bottom-seated circular cylinder is considered from the viewpoint of second-order perturbation theory. After applying the time Fourier transform and separation of vertical variable the resulting two-dimensional non-homogeneous Helmholtz equations are solved numerically using finite differences. The detailed formulation of the second-order radiation condition is presented. Numerical solutions of the problem are obtained for JONSWAP amplitude spectra for the incoming wave group with various types of directional spreading. The results are compared with the corresponding results for a unidirectional wave group of the same amplitude spectrum. Finally we discuss the applicability of the averaged spreading angle concept for practical applications.     

基于HUST-Ship的船舶耐波性数值试验水池参数化建模方法研究

基于自研的HUST-Ship黏性流CFD求解器,提出了一种船舶耐波性数值试验水池参数化建模方法和网格生成技术。在自研的CFD软件平台中,根据波浪参数(规则波的波长、波高和不规则波的特征周期、有效波高)自动生成船舶耐波性数值试验水池。以结构化重叠网格技术为基础,分别对船舶与水池进行网格划分并进行网格的组装、重叠,实现船舶在耐波性数值试验水池中的数值“试验”,对船舶耐波性进行数值预报。根据波浪种类分为规则波数值试验水池与不规则波数值试验水池,可对各类波浪条件实现精确可靠的数值造波,模拟船舶在各种波浪环境中运动响应、耐波性事件与非线性运动特性。采用参数化建模方法生成的耐波性数值试验水池能够提高建模和网格生成效率,通过对国际标模进行数值预报并与船模试验结果对比,验证了参数化建模方法生成的耐波性数值试验水池对船舶在波浪中运动性能的计算精度。     

Solving the linearized forward-speed radiation problem using a high-order finite difference method on overlapping grids

The linearized potential flow approximation for the forward speed radiation problem is solved in the time domain using a high-order finite difference method. The finite-difference discretization is developed on overlapping, curvilinear body-fitted grids. To ensure numerical stability, the convective derivatives in the free-surface boundary conditions are treated using an upwind-biased stencil. Instead of solving for the radiation impulse response functions, a pseudo-impulsive Gaussian type displacement is employed in order to tailor the frequency-content to the discrete spatial resolution. Frequency-domain results are then obtained from a Fourier transform of the force and motion signals. In order to make a robust Fourier transform, and capture the response around the critical frequency, the tail of the force signal is asymptotically extrapolated assuming a linear decay rate. Fourth-order convergence of the calculations on simple geometries is demonstrated, along with a nearly linear scaling of the solution effort with increasing grid resolution. The code is validated by comparison with analytical and semi-analytical solutions using submerged and floating closed-form geometries. Calculations are also made for a modern bulk carrier, and good agreement is found with experimental measurements.     

Nonlinear effects of the distribution of amplitudes of sea waves

The effect of nonlinearities, such as wave-breaking and vertical asymmetry associated with sea waves, on the distribution of wave amplitudes is explored. Semiclosed theoretical expressions are derived to describe the distributions of breaking-limited crest and trough amplitudes for Stokes-type nonlinear sea waves. These are compared with the conventional Rayleigh distribution appropriate to linear wave amplitudes. The construction of nonlinear wave envelopes with the fast Fourier transform technique is described. The technique can be utilized to enlarge the data base in empirical analyses of field records which typically contain limited information on amplitude characteristics. The theoretical distributions and the proposed data enlargement technique are demonstrated with the analysis of a nonlinear wave record.     

Application of Wavelet Algorithm to Spectral Analysis of Oceanic Waves and Offshore Structure Responses

Fourier transform (FT) is a commonly used method in spectral analysis of ocean wave and offshore structure responses,but it is not suitable for records of short length.In this paper another method,wavelet transform (WT),is applied to analyze the data of short length.The Morlet wavelet is employed to calculate the spectral density functions for wave records and simulated Floating Production Storage and Offloading (FPSO) vessels' responses.Computed wave data include simulated wave data based on JONSWAP spectr...     

The limiting wave height in wind-induced wave trains

The vertical acceleration threshold concept has been applied to evaluate the limiting wave height in the train of wind-induced waves propagated over a horizontal bottom. This concept yields very simple computation of the probability of breaking for stochastic sea in deep and finite water depths. The computations confirmed the available field and laboratory observations that the limiting wave steepness in the deep water is lower than the steepness predicted by Stokes. For shallow water depth, the limiting wave height is smaller than 0.55h. This conclusion is consistent with field as well as wave tank observations.     

Hartley变换在海洋工程中的应用研究

Hartley变换是一种基于谐函数的实变换,具有Fourier变换的特性,在海洋工程中使用Fourier变换的大部分应用都可以通过Hartley变换实现,并且快速Hartley变换的计算效率是快速Fourier变换的2倍,并易于实现,所以Hartley变换可以在许多应用中代替Fourier变换。详细介绍了Hartley变换在海洋工程领域的应用,包括:波谱分析、波浪分离、PIV(粒子图像测速法)技术和微分方程的Hartley域变换,举例说明了Hartley变换在节约计算空间和提高计算效率上具备的优越性。     

波浪水槽中随机波的模拟

根据线性系统理论,首先由给定的靶谱模拟海浪信号,再把此信号转换为造波控制输入信号,最后在水槽中指定的截面处产生具有给定靶谱的随机波。这种模拟方法已在青岛海洋大学物理海洋实验室的大型水槽中实现,并取得良好的效果。     

Measurement of ocean wave spectra using a ship-mounted HF radar

The work describes an inversion algorithm for HF radar measurement of nondirectional wave spectra using an omnidirectional receive/transmit antenna. Such a radar would be suitable for deployment on a stationary ship or drill rig. In this approach, wave information is extracted from the radar observations by numerically inverting the integral equation representing the backscatter return from the ocean. Test results of this technique applied to data collected using a 25.4-MHz radar installed on a ship have been very positive. For the two measurements collected, there is a high degree of correlation between the radar wave estimates and those of a WAVE-TRACK buoy     

Wave grouping characteristics in nearshore Great Lakes

The recently advanced approach of wavelet transform is applied to the analysis of wave data measured in the nearshore areas of the Great Lakes. The conventional spectrum analysis of wave time series in the frequency domain can be readily generalized to the frequency and time domain using the wavelet transform. The traditional Fourier transform approach has not been able to directly assess the time localized nature of wave groups. With the application of wavelet transformation, the relatively unexplored wave grouping characteristics come to light as the predominant feature of wave processes.     

The application of Hartley transform to ocean engineering

The Hartley transform is a real integral transform based on harmonic functions and hassome characteristics similar to the Fourier transform. Most applications in ocean engineering requiring the Fourier transform can also be performed by the Hartley transform. The fast Hartley transform is twice faster and more convenient to handle than the corresponding fast Fourier transform, so it is a real valued alternative to the complex Fourier transform in many applications. The use of the Hartley transform in ocean engineering is presented in detail in this paper, including wave spectral analysis, separation of waves, cross-correlation in PIV technique and expression of equation in the Hartley domain. The examples in the paper show deeply the advantage and efficiency of the Hartley transform over the Fourier transform.     

圆形数值波浪水池斜向和多向入射波浪模拟研究

为了研究斜向入射波浪,基于三维不可压缩两相流模型,开发了一套圆形数值波浪水池数值模型。在圆形波浪水池中,通过源项造波法成功生成了任意入射方向的波浪,并且利用人工摩擦项模拟阻尼区以数值耗散反射波浪。模型基于嵌入式多块网格体系,采用FVM法(finite volume method)离散Navier-Stokes方程,VOF法(volume of fluid)追踪自由水面。试验结果表明,斜向入射波浪的模拟结果与理论值基本一致,圆形波浪水池在模拟斜向入射波浪时,有效区域的面积较传统波浪水池显著增大,而且有效区域受波浪入射角度的影响也较小。同时,通过叠加多列斜向入射波浪,模拟出了多向交叉波列,并通过与理论结果对比,发现其具有较高的精度。