非线性波浪波面追踪的一种新模式

基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。     

Simulation of Fully Nonlinear 3-D Numerical Wave Tank

A fully nonlinear numerical wave tank (NWT) has been simulated by use of a three-dimensional higher order bouodary element method (HOBEM) in the time domain. Within the frame of potential flow and the adoption of simply Rankine source, the resulting boundary integral equation is repeatedly solved at each time step and the fully nonlinear free surface boundary conditions are integrated with time to update its position and boundary values. A smooth technique is also adopted in order to eliminate the possible saw-tooth numerical instabilities. The incident wave at the uptank is given as theoretical wave in this paper. The outgoing waves are absorbed inside a damping zone by spatially varying artificial damping on the free surface at the wave tank end. The numerical results show that the NWT developed by these approaches has a high accuracy and good numerical stability.     

Evaluation of the Diurnal Warming of Sea Surface Temperature Using Satellite-Derived Marine Meteorological Data

In order to produce a high-quality sea surface temperature (SST) data set, the daily amplitude of SST (ΔSST) should be accurately known. The purpose of this study was to evaluate the diurnal variation of sea surface temperature in a simple manner. The authors first simulated ΔSST with a one-dimensional numerical model using buoy-observed meteorological data and satellite-derived solar radiation data. When insolation is strong, the model-simulated 1-m-depth ΔSST becomes much smaller than the in situ value as wind speed decreases. By forcibly mixing the sea surface layer, the model ΔSST becomes closer to the in situ value. It can be considered that part of this difference is due to the turbulence induced by the buoy hull. Then, on the assumption that the model results were reliable, the authors derived a regression equation to evaluate ΔSST at the skin and 1-m depth from daily mean wind speed (U) and daily peak solar radiation (PS). ΔSST is approximately proportional to In(U) and (PS)², and the skin ΔSST estimated by the equation is not inconsistent with in situ observation results reported in past studies. The authors prepared maps of PS and U using only satellite data, and demonstrated the ΔSST evaluation over a wide area. The result showed that some wide patchy areas where the skin ΔSST exceeds 3.0 K can appear in the tropics and the mid-latitudes in summer. This revised version was published online in July 2006 with corrections to the Cover Date.     

Free-surface evolution and wave kinematics for nonlinear uni-directional focused wave groups

This paper concerns the propagation of transient wave groups, focused at a point in time and space to produce locally large waves having a range of steepness. The experimental study was carried out in a wave flume at Dalian University of Technology. The numerical simulations were based on a nonlinear boundary integral equation solved by a higher-order boundary element method (HOBEM). Rather than simulate the whole experimental tank, local surface elevation measurements were used to drive the numerical solution from a point less than two wavelengths upstream of the focus position, leading to significant savings in computational time. Excellent agreement is achieved between the water surface elevations and the water particle kinematics measured in the experiments and those predicted numerically at wave group focus, even for near-breaking waves up to a steepness of kA=0.405 for which even locally matched 2nd-order theory is inadequate. Results based on the linear and 2nd-order theory are also presented in the comparisons. When compared with the first- and 2nd-order solutions, the fully nonlinear wave–wave interactions produce a steeper wave envelope in which the central wave crest is higher and narrower, while the adjacent wave troughs are broader and less deep.     

波浪与带窄缝方箱作用共振现象的数值模拟

By introducing a source term into the Laplace equation, a two-dimensional fully nonlinear time-domain numerical wave flume (NWF) is developed to investigate the resonance induced by the interaction bet...     

Wave propagation in a fully nonlinear numerical wave tank: A desingularized method

The problem of wave propagation in a fully nonlinear numerical wave tank is studied using desingularized boundary integral equation method coupled with mixed Eulerian–Lagrangian formulation. The present method is employed to solve the potential flow boundary value problem at each time step. The fourth-order predictor–corrector Adams–Bashforth–Moulton scheme is used for the time-stepping integration of the free surface boundary conditions. A damping layer near the end-wall of wave tank is added to absorb the outgoing waves with as little wave reflection back into the wave tank as possible. The saw-tooth instability is overcome via a five-point Chebyshev smoothing scheme. The model is applied to several wave propagations including solitary, irregular and random incident waves.     

Nonlinear numerical simulation on extreme-wave kinematics

A fully nonlinear numerical model based on a time-domain higher-order boundary element method (HOBEM) is founded to simulate the kinematics of extreme waves. In the model, the fully nonlinear free surface boundary conditions are satisfied and a semi-mixed Euler-Lagrange method is used to track free surface; a fourth-order Runga-Kutta technique is adopted to refresh the wave elevation and velocity potential on the free surface at each time step; an image Green function is used in the numerical wave tank so that the integrations on the lateral surfaces and bottom are excluded. The extreme waves are generated by the method of wave focusing. The physical experiments are carried out in a wave flume. On the horizontal velocity of the measured point, numerical solutions agree well with experimental results. The characteristics of the nonlinear extreme-wave kinematics and the velocity distribution are studied here.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depth

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropriate wave theory. At the downstream boundary, an artificial damping zone is used to prevent wave reflection back into the computational domain. Using the image Green function in the whole fluid domain, the integrations on the two lateral walls and bottom are excluded. The simulation results on extreme wave elevations in finite and infinite water-depths are compared with experimental results and second-order analytical solutions respectively. The wave kinematics is also discussed in the present study.     

Numerical simulation of wave-induced laminar boundary layers

The three-dimensional numerical model with σ-coordinate transformation in the vertical direction is applied to the simulation of surface water waves and wave-induced laminar boundary layers. Unlike most of the previous investigations that solved the simplified one-dimensional boundary layer equation of motion and neglected the interaction between boundary layer and outside flow, the present model solves the full Navier–Stokes equations (NSE) in the entire domain from bottom to free surface. A non-uniform mesh system is used in the vertical direction to resolve the thin boundary layer. Linear wave, Stokes wave, cnoidal wave and solitary wave are considered. The numerical results are compared to analytical solutions and available experimental data. The numerical results agree favorably to all of the experimental data. It is found that the analytical solutions are accurate for both linear wave and Stokes wave but inadequate for cnoidal wave or solitary wave. The possible reason is that the existing analytical solutions for cnoidal and solitary waves adopt the first-order approximation for free stream velocity and thus overestimate the near bottom velocity. Besides velocity, the present model also provides accurate results for wave-induced bed shear stress.     

Simulating nonlinear waves on the free surface in surf zones with two-dimensional sloping beach

Unsteady nonlinear wave motions on the free surface in shallow water and over slopes of various geometries are numerically simulated using a finite difference method in rectangular grid system. Two-dimensional Navier–Stokes equations and the continuity equation are used for the computations. Irregular leg lengths and stars are employed near the boundaries of body and free surface to satisfy the boundary conditions. Also, the free surface which consists of markers or segments is determined every time step with the satisfaction of kinematic and dynamic free surface conditions. Moreover, marker-density method is also adopted to allow plunging jets impinging on the free surface. Either linear or Stokes wave theory is employed for the generation of waves on the inflow boundary. For the simulation of wave breaking phenomena, the computations are carried out with various wave periods and sea bottom slopes in surf zone. The results are compared with other existing computational and experimental results. Agreement between the experimental data and the computation results is good.     

A 3-D model based on Navier–Stokes equations for regular and irregular water wave propagation

A spatial fixed σ-coordinate is used to transform the Navier–Stokes equations from the sea bed to the still water level. In the fixed σ-coordinate system only a very small number of vertical grid points are required for the numerical model. The time step for using the spatial fixed σ-coordinate is efficiently larger than that of using a time dependent σ-coordinate, as there is substantial truncation error involved in the time dependent σ-coordinate transformation. There is no need to carry out the σ-coordinate transformation at each time step, which can reduce computational times. It is important that wave breaking can be potentially modeled in the fixed σ-coordinate system, but in a time-dependent σ-coordinate system the wave breaking cannot be modeled. A projection method is used to separate advection and diffusion terms from the pressure terms in Navier–Stokes equations. The pressure variable is further separated into hydrostatic and hydrodynamic pressures so that the computer rounding errors can be largely avoided. In order to reduce computational time of solving the hydrodynamic pressure equation, at every time step the initial pressure is extrapolated in time domain using computed pressures from previous time steps, and then corrected in spatial domain using a multigrid method. For each time step, only a few of iterations (typically six iterations) are required for solving the pressure equation. The model is tested against available experimental data for regular and irregular waves and good agreement between calculation results and the measured data has been achieved.     

波浪与带窄缝多箱体作用共振现象的模拟研究

针对波浪与带有窄缝多箱体结构作用产生的流体共振问题,建立了基于域内源造波技术的二维非线性时域数值波浪水槽模型,其中自由水面满足完全非线性运动学和动力学边界条件,窄缝内流体引入人工阻尼来等效由于涡旋运动和流动分离引起的黏性耗散,计算域边界采用高阶边界元进行离散。通过模拟三箱体间两窄缝内相对波高变化,并与已发表的数值与实验结果对比,验证了本模型的准确性。同时通过大量的数值计算,分析了箱体数量对窄缝内水体共振频率、共振波高以及对结构反射波高和透射波高的影响。     

波浪沿斜面传播过程的数值模拟

精确模拟非线性波沿斜面传播过程非常困难,为此论文从势函数的边界积分方程出发,建立了一种时域内二维波浪模拟的数值模型,主要用来模拟完全非线性波浪的传播变形过程。论文的数值模型使用高阶二维边界元方法,采用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解带自由表面的线性或完全非线性波浪传播问题。在计算区域一端造出线性或非线性的周期性波浪,另一端采用消除反射波的人工粘性吸收边界。通过与现有理论比较证明了论文数值方法所得结果是准确可靠的。     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

A time‐split, forward‐backward numerical model for solving a nonhydrostatic and compressible system of equations

In this paper, we present a numerical procedure for solving a 2‐dimensional, compressible, and nonhydrostatic system of equations. A forward‐backward integration scheme is applied to treat high‐frequency and internal gravity waves explicitly. The numerical procedure is shown to be neutral in time as long as a Courant–Friedrichs–Lewy criterion is met. Compared to the leap‐frog‐scheme most models use, this method involves only two time steps, which requires less memory and is also free from unstable computational modes. Hence, a time‐filter is not needed. Advection and diffusion terms are calculated with a time step longer than sound‐wave related terms, so that extensive computer time can be saved. In addition, a new numerical procedure for the free‐slip bottom boundary condition is developed to avoid using inaccurate one‐sided finite difference of pressure in the surface horizontal momentum equation when the terrain effect is considered. We have demonstrated the accuracy and stability of this new model in both linear and nonlinear situations. In linear mountain wave simulations, the model results match the corresponding analytical solution very closely for all three cases presented in this paper. The analytical streamlines for uniform flow over a narrow mountain range were obtained through numerical integration of Queney's mathematical solution. It was found Queney's original diagram is not very accurate. The diagram had to be redrawn before it was used to verify our model results. For nonlinear tests, we simulated the famous 1972 Boulder windstorm and a bubble convection in an isentropic enviroment. Although there are no analytical solutions for the two nonlinear tests, the model results are shown to be very robust in terms of spatial resolution, lateral boundary conditions, and the use of the time-split scheme.     

斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式

依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍     

A Finite Element Solution of Wave Forces on Submerged Horizontal Circular Cylinders

In this paper, a numerical model is established for estimating the wave forces on a submerged horizontal circular cylinder. For predicting the wave motion, a set of two-dimensional Navier-Stokes equations is solved numerically with a finite element method. In order to track the moving non-linear wave surface boundary, the Navier-Stokes equations are discretized in a moving mesh system. After each computational time step, the mesh is modified according to the changed wave surface boundary. In order to stabilize the numerical procedure, a three-step finite element method is applied in the time integration. The water sloshing in a tank and wave propagation over a submerged bar are simulated for the first time to validate the present model. The computational results agree well with the analytical solution and the experimental data.Finally, the model is applied to the simulation of interaction between waves and a submerged horizontal circular cylinder.The effects of the KC number and the cylinder depth on the wave forces are studied.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depths

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropr...     

北部湾潮汐数值预报及其分析

采用二维非线性运动方程和连续方程组，以多个分潮同步输入，做耦合数值模拟，进行潮汐数值预报。结果表明，K1潮波等振幅线分布与历史潮波图明显不同，此线呈舌状向北伸展，特别在湾的南部海域，其中央区振幅很小，10cm等振幅线可以伸至北纬19°以北。O1分潮的分布与K1分潮类似；北部湾北部以全日潮为主，湾的南部潮汐性质指标数2．0的等值线呈舌状向湾内神，两侧为不正规日潮，中间为不正规半日潮。4．0等值线向北推移，在19°N以南中央海区预报结果与前人的结果不同；南部海区中央一带潮差很小，最大可能潮差不超过2m，这种分布与前人结果有较大差异。