A new method for the calculation of steady periodic capillary-gravity waves on water of arbitrary uniform depth

This paper presents a method for the calculation of steady periodic capillary-gravity waves on water of arbitrary uniform depth. The method developed by Debiane and Kharif in 1997 for infinite depth is extended to finite depth. The water-wave problem is reduced to a system of nonlinear algebraic equations which is solved using Newton's method. For the resonant configurations, the method does not suffer from the Wilton's failures and is valid for all depths. In addition, it is shown that the method allows the computation of solitary waves and generalized solitary waves.     

水波动力学研究进展

水波动力学,历史悠久,内容丰富,此文仅进行简要的历史回顾,以阐明其进展的思想、智慧、途径以及这门学科对其它学科进展的贡献和影响。随后,引进一个新理论,能模拟三维、非定常、有完全非线性和频散(或色散)性的重力毛细波,在深度任变的水中传播和演化的现象。此外,对二维水波在海滩上之涨落问题,提供一个Lagrange-Euler观点相结合的精确计算方法,以供读者推广和应用,此外,对外力作用下的水波系统共振现象,作些补充论述。      

孤立波与淹没平板相互作用的三维波面和水动力实验研究

在波浪水池中进行了孤立波作用下有限长度和有限宽度淹没平板的三维模型水池实验. 首次应用多目视觉立体重构技术测量局部三维自由表面变形, 该系统的有效测量水平范围为1.7 m$\times $1.6 m. 用4个三分力测力传感器组成水下测力系统, 在不影响波面的情况下测量孤立波对平板的作用力和力矩. 针对波浪不破碎的情况, 选择0.4 m水深和0.16 m波高的来波条件, 平板淹没深度为0.1 m. 实验结果表明, 孤立波经过淹没平板时自由面有明显的三维变形, 导致孤立波波幅的时空变化. 波幅在平板尾缘中心线处达到最大值, 并沿展向逐渐减小. 利用多目视觉立体重构系统得到的波面变化过程与浪高仪给出定点波面时间序列相互印证, 表明建立标识码波面测量方法是有效的. 孤立波对淹没平板作用的水动力载荷变化分为6个典型阶段, 并与利用波面三维重构得到的波面测量标识码并讨论. 基于多目视觉立体重构技术得到了垂向力和俯仰力矩极值点出现时的三维波面形态. 建立的多目视觉立体重构系统将为海洋工程结构物的水池物理模型实验提供新的波面测量手段.      

浅水孤立波在三维浮体上的绕射

浅水域中非线性水波运动的控制方程通常是经过深度平均的Boussinesq方程。然而,这一方程在浮体近旁或水下障碍物附近不再适用,在这些区域,流动在水深方向的变化不容忽略,本文应用匹配渐近展开法和边缘层(edge layer)思想,建立了浅水弱非线性波与三维浮体相互作用的数学模型,作为算例,求解了浅水孤立波在垂直圆柱形浮体上的绕射.本方法可以推广到波在一般浮体上绕射的情况。     

求解浅水波方程的并行物理信息神经网络算法

双曲守恒律方程是一类比较特殊的偏微分方程,其数值求解方法的研究一直是一个热点问题,一个显著特性是即使初始条件是光滑的,其解也可能会发展成间断。浅水波方程作为非线性双曲守恒律方程,由于间断解的存在,其精确求解存在很大困难。针对浅水波方程数值求解问题,本文基于PINN(Physics informed neural networks)反问题网络结构构造新的网络,构造的网络结构包括两个并行的神经网络,其中一个网络与已知状态数据(熵稳定格式加密求出)相关,另一个网络与方程本身相关。利用已知速度数据结合浅水波方程本身求解未知水深,最终通过一些数值算例验证网络的可行性。结果表明,新的网络结构可用于浅水波方程求解,利用速度数据可以较为精确地推算出水深。     

一种用于计算板中缺陷兰姆波散射的混合边界元法

本文介绍了一种将传统边界元法和兰姆波的本征模式函数相结合的混合边界元法，用此方法解决兰姆波在板材中传播时遇到缺陷发生的散射问题，计算了兰姆波通过缺陷后的反射与透射系数。文中给出了若干兰姆波信号幅度与板中不同深度缺陷的相互关系。     

地震波倾斜入射下两相饱和土层响应数值模拟

基于Biot两相饱和多孔介质理论建立了频域内饱和土层的刚度矩阵,对从基岩倾斜入射地震波下饱和土层的波动响应进行了数值研究,入射波采取修正的Kanai-Tajimi谱。通过土层地表面处的水平以及竖直放大系数研究了在不同工况条件下土层的波动响应,研究结果表明,土层的波动响应在入射角接近临界入射角时达到峰值;土层的波动响应随着基岩埋深基以及土层阻尼比的增大显著减小;而土层与基岩的波速比对土层的波动响应的影响存在临界波速比,使土层响应达到极值。对于含有夹层的地质环境,研究了波动响应沿深度的分布,通过与均质土层下响应对比,显示非均质夹层能显著减小土层地表处的波动响应,尤其软夹层能较好地阻止入射波向上的传播。     

Oscillations of an inhomogeneous elastic layer

A method is proposed to analyze wave fields in an elastic layer with elastic properties varying arbitrarily with depth. The method is based on reducing the boundary-value problem to a system of Fredholm integral equations of the second kind, which is then analyzed numerically. Some features of the structure of dispersion sets are analyzed and, in particular, their asymptotes are constructed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 157–164, May–June, 2006.     

Nondestructive Evaluation of Crack Depth in Concrete Using PCA-compressed Wave Transmission Function and Neural Networks

Cracks in concrete are common defects that may enable rapid deterioration and failure of structures. Determination of a crack’s depth using surface wave transmission measurement and the cut-off frequency in the transmission function (TRF) is difficult, in part due to variability of the measurement data. In this study, use of complete TRF data as features for crack depth assessment is proposed. A principal component analysis (PCA) is employed to generate a basis for the measured TRFs for various simulated crack (notch) cases in concrete. The measured TRFs are represented by their projections onto the most significant PCs. Then neural networks (NN), using the PCA-compressed TRFs, are applied to estimate the crack depth. An experimental study is carried out for five different artificial crack (notch) cases to investigate the effectiveness of the proposed method. Results reveal that the proposed method can effectively estimate the artificial crack depth in concrete structures, even with incomplete NN training.     

On the accuracy of finite difference and modal methods for computing tidal and wind wave current profiles

This paper deals with the comparative accuracy of using finite difference grids or a modal representation through the vertical in modelling tidally or wind wave induced current profiles. A point model is used in the vertical, with a no-slip condition at the sea bed. In the finite difference approach the high-shear bottom layer is resolved using either a regular grid on a logarithmic or log-linear transformed co-ordinate or an irregular grid, varying in such a manner as to retain second-order accuracy. The accuracy of these various grid schemes is considered in detail. The relative merits of using either the Crank-Nicolson or Dufort-Frankel time integration methods are considered; in the case of a fine grid in a high-viscosity region, some numerical problems are found with the Dufort-Frankel method. An alternative approach to using a finite difference grid in the vertical, namely a modal (spectral) method, is described. The form of the modes is such that they can accurately resolve the high-shear bottom boundary layer. Calculations show that the thickness of the bottom boundary layer in relation to the total water depth is important in determining the choice of grid transform and rates of convergence of solutions using finite difference or modal methods. However, for the majority of problems the modal solution is numerically attractive owing to its computational efficiency and the ease with which solution algorithms based upon it can be coded in vectorizable form suitable for the new generation of vector computers. The influence of viscosity profile, its time variation and water depth upon tidally induced or wave induced currents is considered. Calculations suggest that near-bed measurements of tidal flow in shallow water together with associated modelling would enable appropriate formulations of eddy viscosity to be determined. Similar measurements, though using a laboratory flume, would be appropriate for wind wave problems.     

面元法求解有限水深船舶兴波及水底压力变化

应用势流理论中的格林函数方法计算了船舶定常运动的水动力参数,将有限水深Kelvin移动兴波源格林函数分解成三部分:简单Rankine源集合、局部扰动项和波函数项。在亚临界和超临界航速时,采用不同的积分顺序来消除被积函数的奇异性。利用面元法在船体表面上分布Kelvin源,计算了有限水深下船体表面的源强、压力分布及表面兴波,比较了有限与无限水深结果的区别和联系,进一步求解了船舶航行时引起的水底压力变化,计算结果与实验测量结果吻合良好。     

A Numerical Method for Wave Propagation in Viscoelastic Stratified Porous Media

This paper presents a numerical method, a transmission matrix method, for the wave propagation in viscoelastic stratified saturated porous media. The wave propagation in saturated media, based on Biot theory, is a coupled problem. In this stratified three-dimensional model we do the Laplace transform for the time variable and the Fourier transform for the horizontal space coordinate. The original problem is transformed into ordinary differential equations with six independent unknown variables, which are only the function of the coordinate of depth. Thus, we get a transmission matrix of the wave problem for each layer. In the process of solution we use numerical method to calculate the eigenvalues and the eigenvectors of the transmission matrices. In the first step of the solution process we can obtain the wave field in the transformed space. The fast Fourier transform (FFT) method is used to do the inverse Laplace and the inverse Fourier transforms to get the solution in the time space. The detailed formulae are derived and some numerical examples are given.     

大尺度结构物的二阶波浪荷载

本文提出了一个在有限水深条件下大尺度物体上二阶波浪荷载的数值计算方法。用表面布源法求二阶散射势的特解,提出了有限水深满足非齐次表面边界条件的格林函数,给出了格林函数的积分和级数形式,研究了表面布源范围对解的影响。     

A special solution of wave dissipation by finite porous plates

The reflection and transmission of water waves caused by a small amplitude incident wave through finite fine porous plates with equal spacing and permeability in an infinitely long open channel of constant water depth and zero slope are studied. A special solution is obtained when the distance between the two neighbouring plates is an integral multiple of the half-wavelength of the incident wave. It is found, that when the dimensionless porous-effect parameter G_0 is equal to half the total plate number, the wave dissipation reaches a maximum, and only 50% of the incident wave energy remains in the reflected and transmitted waves. Meanwhile, the reflected and transmitted waves have the same amplitude.     

高填方路堤强夯效果的现场检测及三维有限元模拟

针对河南某岩石碎渣高填方路堤,提出了夯击能1000kN·m、虚铺厚度1.5m的强夯施工工艺,对试验段进行了现场检测：地质雷达波形较好地反映强夯后虚铺层厚度的减小幅度,并可判断强夯水平有效加固范围为4m,瑞利波速沿深度变化比较均匀,平均值达到了310m·s-1,重型动力触探结果表明试验段达到中密程度,灌水法测得的固体体积率都达到了83%以上,上述检测结果说明强夯的加固效果比较理想。采用ABAQUS有限元分析软件,将夯锤简化为刚性体,通过施加竖向初速度的方式模拟了强夯作用,结合雷达测试结果,认为夯沉量20mm可作为强夯的有效加固范围,最后模拟了存在涵洞结构物的路堤的强夯作用,指出涵洞上方路堤厚度必须达到8m以上,才能保证强夯安全作业。     

Head-on collision between two mKdV solitary waves in a two-layer fluid system

In this paper, based on the equations presented in [2], the head-on collision between two solitary waves described by the modified KdV equation (the mKdV equation, for short) is investigated by using the reductive perturbation method combined with the PLK method. These waves propagate at the interface of a two-fluid system, in which the density ratio of the two fluids equals the square of the depth ratio of the fluids. The second order perturbation solution is obtained. It is found that in the case of disregarding the nonuniform phase shift, the solitary waves preserve their original profiles after collision, which agrees with Fornberg and Whitham's numerical result of overtaking collision161 whereas after considering the nonuniform phase shift, the wave profiles may deform after collision.     

梯度介质半空间Rayleigh面波传播性能研究

对剪切弹性模量沿深度以指数函数变化的非均质半空间，本文用摄动法得到了Rayleigh面波的波函数解答及相速度方程。以不同金属与陶瓷复合而成的几种梯度材料为例，用数值方法求解了相速度方程，给出了相应的波的弥散曲线，结果表明，梯度介质半空间自由表面附近的Rayleigh波通常有两种不同的弥散形式，即正常弥散和非正常弥散。     

基于Boussinesq方程的波浪模型

从欧拉方程出发，提供了另一种推导完全非线性Boussinesq方程的方法，并对方程的 线性色散关系和线性变浅率进行了改进. 改进后方程的线性色散关系达到了一阶Stokes波 色散关系的Pad\'{e}[4,4]近似，在相对水深达1.0的强色散波浪时仍保持较高的准确性，并且方程的非线性和线性 变浅率都得到了不同程度的改善. 方程的水平一维形式用预估-校正的有限差分格式求解， 建立了一个适合较强非线性波浪的Boussinesq波浪数值模型. 作为验证，模拟了波浪在潜 堤上的传播变形，计算结果和实验数据的比较发现两者符合良好.     

The wave attenuation on the mud bed

In the paper the wave attenuation in a two layer fluid system is studied. The fluid in the top layer is ideal and that in the lower layer is the Voigt model of the viscoelastic medium. A dispersion relation is derived and the rate of the wave decay is computed. The approximate explicit expressions of the decay rate for different water depth are given, where the viscoelasticity is either very large or very small. Compared with the numerical results, our results are very accurate, which can be used by an engineer. The project supported by the National Natural Science Foundation of China and by the Lianyungang Port Office, China