基于Boussinesq方程的波浪分层数值模拟方法

对波浪进行精确的数值模拟一直是人们所关注的，采用在Boussinesq方程的基础上将垂向流速断面进行分层，在每一层中对流速断面进行积分，求得各层平均流速，用此平均流速来代替每一层的流速，并在各层间加上界面边界条件，以求出流速和自由表面位移，这样可以得出精确的波浪数学模型。本文例举了规则波经过淹没潜堤时的传播变形实验，验证了二层Boussinesq模型在实际运用中的精确性。     

基于Boussinesq方程的数值波浪水池

基于水深积分的非线性Boussinesq方程,采用高精度差分格式,开发了一个具有前后处理、可视化的模拟浅水波浪的通用数学模型系统NWB-BOUS。通过对一典型算例进行验算,表明该系统对波浪在港湾结构物附近的折射、绕射、反射及非线性浅化等多种现象都能进行很好的模拟。     

Boussinesq方程数学模型的域内造波源函数研究

研究Boussinesq方程模型的域内造波源函数,理论推导该数学模型的源函数,证明其中的波速应使用能量速度,即群速度。域内生成正弦波和椭圆余弦波的数值试验结果证实了该结论。     

基于改进型Boussinesq方程的二维波浪数值模型

基于改进型Boussinesq方程，在非交错网格下，建立了二维波浪数值模型。模型计算采用了有限差分法，时间格式上采用混合四阶Adams—Bashforth—Mouhon，空间格式上采用了Wei等（1995）给出的格式。数值计算中。采用了内部造渡技术。数值模拟针对3纽经典浅滩地形上波浪传播变形的实验进行，数值计算结果与实验结果吻合较好。验证了数值模型，该模型可期望用于实际港口渡浪预报。     

交错网格下四阶Boussinesq方程对潜堤上波浪演化的应用*

基于Madsen 和 Sch?ffer推导的四阶Boussinesq方程,在交错网格下建立一维数值模型。模型中的时间导数采用混合四阶Adams-Bashforth-Moulton 格式。为了避免二次反射,采用了Gobbi 和 Kirby的的内部源函数造波法。最后数值模拟了波浪在潜堤上的演化过程,并且与Luth的试验值进行对比,结果显示两者吻合较好,另外还给出了非交错网格下的模拟结果,通过对比可知,两个模型都能够很好地模拟波浪在潜堤上的演化过程。     

不同时空格式在求解Boussinesq水波模型中的应用

为了解不同时间空间差分格式在常用的二阶或四阶Boussinesq模型中的应用,针对4组近似到0（旷）阶完全非线性的二阶或四阶色散性的高阶Boussinesq水波方程,在非交错网格下,利用Crank-Nicolson格式、蛙跳格式、混合四阶Adams-Bashforth-Moulton格式,建立不同的数学模型。利用这些数值模型模拟波浪在潜堤上的传播变形,通过数值结果与试验结果的比较．考察时间格式及空间格式对模型的影响。结果表明：对同一方程,混合四阶Adams-Bashforth-Mouhon格式和Crank-Nicolson格式均能取得较好模拟效果,蛙跳格式的模拟效果最差;二阶Boussinesq模型采用追赶法求解已能满足要求;对四阶Boussineq模型,二阶空间导数色散项亦采用四阶精度,其数值效果会更好。     

二阶全非线性Boussinesq方程的二维数值模拟研究

在非交错网格下采用有限差分法首次对一组非线性精确至()Ομ2阶的全非线性Boussinesq方程数学模型进行了二维数值模拟分析。首先通过在方程的非线性项中引入缓坡假定,考察了其对模型数值精度的影响;其次,在模型中对二阶非线性项采用不同精度,考察了其对模型数值结果的影响。数值模拟结果表明,所建立的二阶完全非线性Boussinesq方程二维数值模型具有良好的适用性,模型非线性项中引入缓坡假定以及在二阶非线性项选用不同的精度对数值模拟结果影响不明显。     

双层Boussinesq模型非线性波浪模拟研究

高精度模拟波浪在潜堤上的传播演变过程,要求数学模型具有优良的色散性、变浅性和非线性.为了考察双层Boussinesq方程模型模拟强非线性问题的精度,针对1:4潜堤上非线性较强的非破碎波传播变形开展了数值模拟研究,并设计了物理模型试验,确定了强非线性工况.数值模型求解利用预报-校正的有限差分法,时间差分格式采用混合4阶A...     

多孔介质中波浪传播的高阶Boussinesq方程*

准确模拟波浪在多孔介质中传播变形对于研究抛石防波堤等结构的消能作用是十分必要的。对Laplace方程、自由表面处的运动学方程和动力学方程以及海底运动学方程进行无因次化,且以自由表面处速度势为切点,进行幂级展开,最终给出4个不同的高阶Boussinesq水波方程。在常水深下对这些方程的一维问题进行了理论研究,并将无因次相速度和无因次虚波数与解析解结果进行对比,方程的相速度与解析解吻合程度较好,虚波数与解析解基本吻合,表明高阶Boussinesq方程可用于模拟波浪在多孔介质中的传播变形。     

基于Boussinesq方程的波浪破碎模型的研究综述

介绍了Boussinesq型波动方程的发展,并对基于Boussinesq方程的破碎模型和波生流计算相关的技术问题进行了系统的归纳总结和评述,以期对本研究领域的发展起到一定的促进作用。     

Boussinesq方程波浪数学模型计算马来西亚科美纳河河口的波浪场

Ｂｏｕｓｓｉｎｅｓｑ方程波浪数学模型计算马来西亚科美纳河河口单堤整治方案的波浪场,为分析整治建筑物对河口地区水动力条件的影响提供依据。     

Crack detection using a frequency response function in offshore platforms

Structural cracks can change the frequency response function (FRF) of an offshore platform. Thus, FRF shifts can be used to detect cracks. When a crack at a specific location and magnitude occurs in an offshore structure, changes in the FRF can be measured. In this way, shifts in FRF can be used to detect cracks. An experimental model was constructed to verify the FRF method. The relationship between FRF and cracks was found to be non-linear. The effect of multiple cracks on FRF was analyzed, and the shift due to multiple cracks was found to be much more than the summation of FRF shifts due to each of the cracks. Then the effects of noise and changes in the mass of the jacket on FRF were evaluated. The results show that significant damage to a beam can be detected by dramatic changes in the FRF, even when 10% random noise exists. FRF can also be used to approximately locate the breakage, but it can neither be efficiently used to predict the location of breakage nor the existence of small hairline cracks. The FRF shift caused by a 7% mass change is much less than the FRF shift caused by the breakage of any beam, but is larger than that caused by any early cracks.     

A B-spline Galerkin scheme for calculating the hydroelastic response of a very large floating structure in waves

This paper presents an effective scheme for computing the wave-induced hydroelastic response of a very large floating structure, and a validation of its usefulness. The calculation scheme developed is based on the pressure-distribution method of expressing the disturbance caused by a structure, and on the mode-expansion method for hydroelastic deflection with the superposition of orthogonal mode functions. The scheme uses bi-cubic B-spline functions to represent unknown pressures, and the Galerkin method to satisfy the body boundary conditions. Various numerical checks confirm that the computed results are extremely accurate, require relatively little computational time, and contain few unknowns, even in the region of very short wavelengths. Measurements of the vertical deflections in both head and oblique waves of relatively long wavelength are in good agreement with the computed results. Numerical examples using shorter wavelengths reveal that the hydroelastic deflection does not necessarily become negligible as the wavelength of incident waves decreases. The effects of finite water depth and incident wave angle are also discussed.     

基于一种固体区域迭代算法的圆柱涡激振动数值计算

利用Fluent平台的用户自定义程序（UDF）以及动网格模型,实现了圆柱运动方程的一种迭代求解算法,分别对层流、湍流状态下,弹性支承圆柱体在一定约化速度下的涡激响应进行了数值模拟,探讨了不同阻尼比对涡激响应的影响。结果表明：采用该迭代求解算法对弹性支承圆柱涡激振动的预测结果较为合理;随着阻尼比的逐渐增加,初始支振幅、升阻力系数时程曲线将由多频率拍振,最终变为单一频率主导的振动,且涡激振幅逐渐减小;除了质量-阻尼比联合参数m*ζ外,阻尼比ζ本身也应作为一个重要的涡激影响参数单独进行考量。     

双材料中矩形裂纹问题的超奇异积分方程方法

使用超奇异积分方程方法,对双材料空间中垂直于界面的矩形裂纹Ⅰ型问题进行了研究.首先根据双材料空间的弹性力学基本解,使用边界积分方程方法,在有限部积分的意义下导出了以裂纹面位移间断为未知函数的超奇异积分方程.根据裂纹面上位移函数的分布特性,通过将位移间断函数表示为特征函数和一组多项式乘积的形式,为其建立了数值方法.数值结果表明,该方法不仅具有较好的收敛性和较高的数值计算精度,而且能够精确满足裂纹面上的边界条件.在此基础上,对不同材料组合界面对裂纹前沿应力强度因子的影响进行了分析,取得了较好的数值结果.     

全垫升气垫船内部流场的数值模拟与改进

应用计算流体力学(CFD)方法对全垫升气垫船气囊内部的气动特性进行数值模拟计算,获得了气囊内的三维流场分布.其中的关键问题是气囊的排气孔壁面的处理,文中采用了多孔介质模型进行近似,针对气囊前后端压力不平衡的问题,分别采用在进气道内设置导流叶片和在气囊通道中设置隔断两种方法使压力趋于平衡,并通过数值试验找出其最佳位置.此外,文章还对实船相似设计的合理性进行了讨论.     

Frequency-Domain 3D Computer Program for Predicting Motions and Loads on a Ship in Regular Waves

The development of an in-house computer program for determining the motions and loads of advancing ships through sea waves in the frequency domain,is described in this paper.The code is based on the potential flow formulation and originates from a double-body code enhanced with the regular part of the velocity potential computed using the pulsing source Green function.The code is fully developed in C++ language with extensive use of the object-oriented paradigm.The code is capable of estimating the excitation and inertial radiation loads or arbitrary incoming wave frequencies and incidence angles.The hydrodynamic responses such as hydrodynamic coefficients,ship motions,the vertical shear force and the vertical bending moment are estimated.A benchmark container ship and an LNG carrier are selected for testing and validating the computer code.The obtained results are compared with the available experimental data which demonstrate the acceptable compliance for the zero speed whereas there are some discrepancies over the range of frequencies for the advancing ship in different heading angles.     

Estimation of TKE dissipation rates in dense bottom plumes using a Pulse Coherent Acoustic Doppler Profiler (PC-ADP) — Structure function approach

During a hydrographic survey in January 2006 the spreading of inflowing saline water was observed in the Arkona Basin (Western Baltic Sea). Two bottom mounted ‘pulse coherent’ acoustic Doppler profilers (PC-ADP) were used to measure the near-bottom current field of the dense plume with a high temporal (1 s) and spatial resolution (5 cm). In order to estimate the dissipation rate of turbulent kinetic energy () a structure function approach was applied to the beam velocity data. Simultaneous measurements with a microstructure shear profiler (MSS) and an acoustic Doppler velocimeter (ADV) supplied independent data for the verification of the structure function method. Additional measurements with standard CTD, near-bottom towed and vessel mounted acoustic Doppler current profilers (ADCP) completed the data set.The estimated dissipation rates from the structure function approach fit well with the values derived from the ADV and the MSS probe. It is shown that the structure function approach is a reliable and easily applicable method to derive estimates of TKE dissipation rates from PC-ADP beam velocities. The observed dissipation rates ranged between 5 · 10^− 6 and 1 · 10^− 8 W kg^− 1 depending on the hydrographic conditions. Inside the plume the dissipation rates exceeded that of the overlaying brackish water by two orders of magnitude. Since the noise level of velocity data in pulse coherent mode is considerably lower than in the Doppler mode the PC-ADP can also be used for estimates in marine environments with low turbulence level. Reynolds stresses estimated from the PC-ADP and the ADV agreed well at the same depth level. TKE production derived from PC-ADP measurements compared reasonably well with the dissipation rate of TKE in a varying environment.