潜堤上波流传播的完全非线性数值模拟

利用时域高阶边界元方法建立模拟波流混合作用的完全非线性数值水槽模 型, 其中自由水面满足完全非线性自由水面条件. 采用混合欧拉-拉格朗日方法追踪流体瞬时 水面, 运用4阶Runge-Kutta方法更新下一时间步的波面和速度势. 为了减少计算域, 提高 计算速度, 采用同时消除底面和侧面的镜像格林函数; 在每一时间步内, 对网格进行重新划分 以避免由于网格运动变形而引起的数值不稳定问题. 对水流中淹没潜堤上的波浪变形在水槽 中开展了物理模型试验, 并把试验结果和数值结果进行了对比, 吻合得很好. 进一步研究了 水流及潜堤的存在对高阶谐波产生的影响.     

波流与结构物相互作用的数值模拟

采用高阶边界元法建立波流与任意形状结构物相互作用的时域数学模型。在小流速假设下,将速度势进行摄动展开,边界条件分解为双物体假设下的零阶波陡稳定问题和一阶波陡下的不稳定波浪问题。物面速度势及自由水面速度势的法向导数等未知量通过求解高阶边界元积分方程得到,而积分方程的求解则通过一个数值程序来实现。在每个时间步上采用四阶Runge-Kutta法更新下一时刻自由水面波面和速度势,自由水面上采用人工阻尼层消除散射波。通过波流与直立圆柱相互作用的数值计算研究了一阶激振力和二阶慢漂力随波数的变化关系以及圆柱周围波幅的分布曲线,并与已有频域结果和时域结果对比验证了所建模型的准确性,进而应用本文模型研究了波流与实际工程结构物相互作用的问题。     

波浪作用下的泥沙沉降速度

本文认为泥沙在波动水流中的沉降与在静水中的沉降不同,即泥沙还受波浪紊动的举力作用而使沉降速度减慢.波浪运动由水面向下逐渐衰减,所以,波浪对泥沙沉降速度的减缓作用也由水面向下逐渐衰减.由此出发推得了相应的理论公式,计算结果与实验数据的比较在规律性上是一致的,但计算值的精度尚嫌不足.     

基于预修正快速傅里叶变换高阶边界元方法的多体水动力分析

高阶边界元方法在求解波浪对海上建筑物作用问题中具有诸多优势，但由于它所形成的矩阵是一个满阵，计算量和存储量均为未知量的平方量级，很难满足如多体水动力分析等大尺度多未知量问题的计算需要。本文采用预修正快速傅里叶变换高阶边界元方法（pFFT-HOBEM），将计算量和存储量均降低到未知量的线性量级。通过对不同未知量时该方法与传统边界元法的计算量与存储量的对比，以及该方法自身各步骤计算时间的对比，研究了不同pFFT网格方案对计算量和存储量的影响，并提出了基于计算时间最小化原则的pFFT网格优化方法。采用本文方法研究了四柱结构在不同频率波浪作用下的作用力及波面分布，并对近场干涉发生时的物理现象进行了分析和讨论。     

台阶式溢洪道滑行水流水面线和消能效果的试验研究

通过模型试验研究了台阶式溢洪道滑行水流的流态、水面线和消能效果。试验得出,台阶式溢洪道的水流流态分为台阶内部的旋滚流态和虚拟底板以上主流区的流态。台阶式溢洪道水面线的变化与来流量、台阶尺寸和掺气浓度有关。当水流进入台阶段时,水面由于失重有一定的降低,然后逐渐抬高至某一高度后再略有降低;在掺气发生点以前,水面线沿程降低,掺气发生点以后,水面线沿程升高;当水流掺气饱和时,水深沿程不再变化,为明渠均匀流。台阶式溢洪道的消能率是逐级台阶的累积效应。水流通过与台阶之间的碰撞、台阶上的水流漩滚以及水流内部的紊动剪切强化了溢洪道的消能效果。在试验范围内,台阶式溢洪道与光滑溢洪道比较,消能率可以提高40% ~70%。通过试验,得出了台阶式溢洪道滑行水流水面线和消能率的计算方法。     

非周期波浪与直墙作用的非线性数值研究

基于时域高阶边界元方法,建立了完全非线性二维数值波浪水槽,对非周期波浪与直墙的相互作用问题进行了模拟和研究.自由表面满足完全非线性自由水面运动学和动力学边界条件,采用混合欧拉-拉格朗日方法追踪瞬时自由面流体质点,采用四阶Runge-Kutta法对下一时间步的波面和自由面速度势进行更新.采用加速度式法求解直墙表面速度势的时间导数,对瞬时物体湿表面上的水动力压强积分,得到作用在物体上的瞬时波浪力.首先,将全非线性与Serre-Green-Naghdi(SGN)模型的结果进行了对比分析,发现对于大幅值双入射波问题,仅满足弱色散关系的SGN模型大大低估了最大波浪爬高;其次,研究了双入射波与直墙的非线性作用问题,发现线性预报对波浪最大爬高有较大低估,而波浪的非线性成分不只导致了自由面爬高的异常增大,也引起了局部自由面的高频振荡,该物理过程中,直墙所受的波浪载荷,也展示出了与波浪爬高相似的非线性特性;最后,对波浪爬升和波浪力的时间历程进行了频谱分析,发现入射主频波的部分能量传递给了更高频的波浪成分,反映出该问题具有典型的非线性特性.     

非均匀水流中非线性波传播的数值模拟

以一种考虑波流相互作用的新型{Boussinesq}型方程为控制方程组， 采用五阶{Runge}-{Kutta}-{England}格式离散时间积分，采用七点 差分格式离散空间导数，并通过采用恰当的出流边界条件，从而建立了非均匀水流中非线性 波传播的数值模拟模型. 通过对均匀水流与水深水域内和潜堤地形上存在弱流或强流时波浪 传播的数值模拟，说明模型能有效地反映水流对波浪传播的影响.     

模拟畸形波的聚焦波浪模型

利用改进的高阶谱方法建立了模拟极限波的二维聚焦模型，通过与Baldock(1996) 的实验结果和理论 值的比较，验证了模型的正确性，并分析了波浪非线性的相互作用对聚焦结果的影响. 通过 改进Longuet-Higgins海浪模型，给出了4种实验室聚焦模拟畸形波的波浪模型：极限波 聚焦模型+随机波模型；极限波聚焦模型+规则波模型；相位角分布范围调制聚焦模型； 相同相位角组成波个数调制聚焦模型. 基于上述完全非线性数值波浪模型，采用不同的能 量分配方式，在有限模拟长度和时间内得到了具有不同$H_{\max}/H_{s}$值的畸形波.     

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡,了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提.基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽,模型采用S-A IDDES湍流模型进行湍流封闭,并采用修正的VOF方法捕捉自由液面,数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程,重点分析了结构物附近的水动力和掺气水流运动特性.结果表明,建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程,较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程;波浪与直立结构物相互作用产生强湍动掺气水流,其运动过程十分复杂;掺气流体输运过程中水气界面周围一直伴随着涡的存在,其中,气泡的分裂与周围正负涡量剪切作用密切相关,且其输运轨迹主要受周围流场的影响;研究揭示了结构物附近湍动能与掺气特性的关系,发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.     

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡, 了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提. 基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽, 模型采用S-A IDDES湍流模型进行湍流封闭, 并采用修正的VOF 方法捕捉自由液面, 数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程, 重点分析了结构物附近的水动力和掺气水流运动特性. 结果表明, 建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程, 较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程; 波浪与直立结构物相互作用产生强湍动掺气水流, 其运动过程十分复杂; 掺气流体输运过程中水气界面周围一直伴随着涡的存在, 其中, 气泡的分裂与周围正负涡量剪切作用密切相关, 且其输运轨迹主要受周围流场的影响; 研究揭示了结构物附近湍动能与掺气特性的关系, 发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.      

Phase velocity effects of the wave interaction with exponentially sheared current

The nonlinear interaction between the unidirectional bichromatic wave-train and exponentially sheared current in water of an infinite depth is investigated. The model is based on the vorticity transport equation and the exact free surface conditions, without any assumptions for the existence of small physical parameters. Earlier works of the wave–current interaction were mainly restricted to either current acted on the monochromatic wave or irregular waves limited to irrotational current. Different from these previous works, no constraint is made in our model for amplitudes of the primary wave, and the current owns an exponential type profile along the vertical line. To ensure that the effect of vorticity on the phase velocity is consistent with earlier derivation, the case of a small amplitude wave traveling on the exponentially sheared current is examined firstly. Then the effect of nonlinearity on the phase velocity of primary waves in a bichromatic wave-train is considered. Accurate high-order approximations of the phase velocity are obtained under consideration of both the nonlinear wave self–self and mutual interactions. Finally, the combined effect of vorticity and nonlinearity on the phase velocity is investigated through the case of a bichromatic wave-train propagating on an exponentially sheared current. It is found that the characteristic current slope determines the effect of vorticity on the phase velocity caused by nonlinear wave self–self and mutual interactions, and the surface current strength may amplify/reduce this effect.     

3D free‐surface flow computation using a RANSE/Fourier–Kochin coupling

A coupling method for numerical calculations of steady free‐surface flows around a body is presented. The fluid domain in the neighbourhood of the hull is divided into two overlapping zones. Viscous effects are taken in account near the hull using Reynolds‐averaged Navier–Stokes equations (RANSE), whereas potential flow provides the flow away from the hull. In the internal domain, RANSE are solved by a fully coupled velocity, pressure and free‐surface elevation method. In the external domain, potential‐flow theory with linearized free‐surface condition is used to provide boundary conditions to the RANSE solver. The Fourier–Kochin method based on the Fourier–Kochin formulation, which defines the velocity field in a potential‐flow region in terms of the velocity distribution at a boundary surface, is used for that purpose. Moreover, the free‐surface Green function satisfying this linearized free‐surface condition is used. Calculations have been successfully performed for steady ship‐waves past a serie 60 and then have demonstrated abilities of the present coupling algorithm. Copyright © 2003 John Wiley & Sons, Ltd.     

A BOUNDARY ELEMENT METHOD FOR STEADY SHIP WAVE-MAKING POTENTIAL PROBLEMS

A boundary element method for three-dimensional steady ship wave-making potentialproblems is established with the Rankine souree function as its fundamental solution.In the treatmentof the linearized free surface condition,one sided.upstream finite difference operator(FDO)is usedto suppress the upstream waves,and the equation of the disturbance velocity is established so that thefirst order FDO can be used in place of the second order FDO.Compared with the method with the sec-ond order FDO,the current method gives better precision and stability.Numerical examples are pre-sented for verification.     

Semi‐coupled air/water immersed boundary approach for curvilinear dynamic overset grids with application to ship hydrodynamics

For many problems in ship hydrodynamics, the effects of air flow on the water flow are negligible (the frequently called free surface conditions), but the air flow around the ship is still of interest. A method is presented where the water flow is decoupled from the air solution, but the air flow uses the unsteady water flow as a boundary condition. The authors call this a semi‐coupled air/water flow approach. The method can be divided into two steps. At each time step the free surface water flow is computed first with a single‐phase method assuming constant pressure and zero stress on the interface. The second step is to compute the air flow assuming the free surface as a moving immersed boundary (IB). The IB method developed for Cartesian grids (Annu. Rev. Fluid Mech. 2005; 37 :239–261) is extended to curvilinear grids, where no‐slip and continuity conditions are used to enforce velocity and pressure boundary conditions for the air flow. The forcing points close to the IB can be computed and corrected under a sharp interface condition, which makes the computation very stable. The overset implementation is similar to that of the single‐phase solver (Comput. Fluids 2007; 36 :1415–1433), with the difference that points in water are set as IB points even if they are fringe points. Pressure–velocity coupling through pressure implicit with splitting of operators or projection methods is used for water computations, and a projection method is used for the air. The method on each fluid is a single‐phase method, thus avoiding ill‐conditioned numerical systems caused by large differences of fluid properties between air and water. The computation is only slightly slower than the single‐phase version, with complete absence of spurious velocity oscillations near the free surface, frequently present in fully coupled approaches. Validations are performed for laminar Couette flow over a wavy boundary by comparing with the analytical solution, and for the surface combatant model David Taylor Model Basin (DTMB) 5512 by comparing with Experimental Fluid Dynamics (EFD) and the results of two‐phase level set computations. Complex flow computations are demonstrated for the ONR Tumblehome DTMB 5613 with superstructure subject to waves and wind, including 6DOF motions and broaching in SS7 irregular waves and wind. Copyright © 2008 John Wiley & Sons, Ltd.     

Unsteady RANS method for surface ship boundary layer and wake and wave field

Results are reported of an unsteady Reynolds‐averaged Navier–Stokes (RANS) method for simulation of the boundary layer and wake and wave field for a surface ship advancing in regular head waves, but restrained from body motions. Second‐order finite differences are used for both spatial and temporal discretization and a Poisson equation projection method is used for velocity–pressure coupling. The exact kinematic free‐surface boundary condition is solved for the free‐surface elevation using a body‐fitted/free‐surface conforming grid updated in each time step. The simulations are for the model problem of a Wigley hull advancing in calm water and in regular head waves. Verification and validation procedures are followed, which include careful consideration of both simulation and experimental uncertainties. The steady flow results are comparable to other steady RANS methods in predicting resistance, boundary layer and wake, and free‐surface effects. The unsteady flow results cover a wide range of Froude number, wavelength, and amplitude for which first harmonic amplitude and phase force and moment experimental data are available for validation along with frequency domain, linear potential flow results for comparisons. The present results, which include the effects of turbulent flow and non‐linear interactions, are in good agreement with the data and overall show better capability than the potential flow results. The physics of the unsteady boundary layer and wake and wave field response are explained with regard to frequency of encounter and seakeeping theory. The results of the present study suggest applicability for additional complexities such as practical ship geometry, ship motion, and maneuvering in arbitrary ambient waves. Copyright © 2001 John Wiley & Sons, Ltd.     

Study on the head-on collisions of two-dimensional trains of solitons

The head on collisions of trains of solitons induced by a two-dimensional submerged elliptical cylinder at critical speed in shallow water are studied based on velocity potential theory. The boundary value problems are solved through boundary element method (BEM). The nonlinear free surface boundary conditions are satisfied. The mixed Euler–Lagrangian method is adopted to track the free surface through a time stepping scheme. The effects of thickness and velocity of the elliptical cylinder on the evolution of solitary waves have been investigated. Two sets of solitons are truncated from these trains of solitary waves. The head-on collisions of these solitons have been simulated. The wave profiles and velocity fields during collision have been analysed. The propagation of solitary waves is the transmissions of kinetic energy and the collision processes are the results of the dynamic balance of potential energy and kinematic energy.     

An explicit formulation for the evolution of nonlinear surface waves interacting with a submerged body

An explicit formulation to study nonlinear waves interacting with a submerged body in an ideal fluid of infinite depth is presented. The formulation allows one to decompose the nonlinear wave–body interaction problem into body and free‐surface problems. After the decomposition, the body problem satisfies a modified body boundary condition in an unbounded fluid domain, while the free‐surface problem satisfies modified nonlinear free‐surface boundary conditions. It is then shown that the nonlinear free‐surface problem can be further reduced to a closed system of two nonlinear evolution equations expanded in infinite series for the free‐surface elevation and the velocity potential at the free surface. For numerical experiments, the body problem is solved using a distribution of singularities along the body surface and the system of evolution equations, truncated at third order in wave steepness, is then solved using a pseudo‐spectral method based on the fast Fourier transform. A circular cylinder translating steadily near the free surface is considered and it is found that our numerical solutions show excellent agreement with the fully nonlinear solution using a boundary integral method. We further validate our solutions for a submerged circular cylinder oscillating vertically or fixed under incoming nonlinear waves with other analytical and numerical results. Copyright © 2007 John Wiley & Sons, Ltd.     

A two‐dimensional vertical non‐hydrostatic σ model with an implicit method for free‐surface flows

An implicit finite difference model in the σ co‐ordinate system is developed for non‐hydrostatic, two‐dimensional vertical plane free‐surface flows. To accurately simulate interaction of free‐surface flows with uneven bottoms, the unsteady Navier–Stokes equations and the free‐surface boundary condition are solved simultaneously in a regular transformed σ domain using a fully implicit method in two steps. First, the vertical velocity and pressure are expressed as functions of horizontal velocity. Second, substituting these relationship into the horizontal momentum equation provides a block tri‐diagonal matrix system with the unknown of horizontal velocity, which can be solved by a direct matrix solver without iteration. A new treatment of non‐hydrostatic pressure condition at the top‐layer cell is developed and found to be important for resolving the phase of wave propagation. Additional terms introduced by the σ co‐ordinate transformation are discretized appropriately in order to obtain accurate and stable numerical results. The developed model has been validated by several tests involving free‐surface flows with strong vertical accelerations and non‐linear waves interacting with uneven bottoms. Comparisons among numerical results, analytical solutions and experimental data show the capability of the model to simulate free‐surface flow problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Non-linear free surface flows and an application of the Orlanski boundary condition

The focus of this paper is the analysis of spatially two-dimensional non-linear free surface problems. The critical aspects of the problem concern the treatment of the non-linear free surface, the body boundary condition for large motions and the imposition of suitable radiation conditions. To address such complexities, time domain simulation was chosen as the method of analysis. With the use of a finite domain for simulation, a major concern is with the radiation condition to be applied at the open or truncation boundary. For the two-dimensional problem at hand, no theoretical radiation conditions are known to exist. An extension of the Orlanski open boundary condition, based on phase velocity determination at the free surface, is proposed. Three categories of problems were analysed using numerical simulation-namely, freely moving steep waves, waves over a submerged body and forced body motion. Simulation results have been compared with linear theory and experiments.