近岸大区域水波数学模型及其数值求解

从水波发展方程出发,导出了大区域缓坡水波数学模型,并采用交错网格系统下的有限差分法对该数学模型进行了数值求解,计算结果表明该数学模型在粗网格下也能得到与实验结果比较一致的结果,从而表明该模型可用于较大区域的水波问题.该模型具有编程简单、求解比较快速、经济的优点.     

Numerical Simulation of Breaking Wave Based on Higher-Order Mild Slope Equation

The "surface roller" to simulate wave energy dissipation of wave breaking is introduced into the random wave model based on approximate parabolic mild slope equation in this paper to simulate the random wave transportation including diffraction, refraction and breaking in nearshore areas. The roller breaking random wave higher-order approximate parabolic equation model has been verified by the existing experimental data for a plane slope beach and a circular shoal, and the numerical results of random wave breaking model agree with the experimental data very well. This model can be applied to calculate random wave propagation from deep to shallow water in large areas near the shore over natu ral topography.     

Numerical Modeling of Wave Diffraction-Refraction in Water of Varying Current and Topography

—A numerical model for wave diffraction-refraction in water of varying current and topogra-phy is proposed,and time-dependent wave mild-slope equation with a dissipation term and correspondingequivalent governing equations are presented.Two different expressions of parabolic approximations forthe case of the absence of current are also given and analyzed.The influence of current on the results ofsimulation of waves is discussed.Some examples show that the present model is better than others in simu-lating wave transformation in large water areas.And they also show that the influence of current shouldbe taken into account,on numerical modeling of wave propagation in water of strong current and coastalareas,otherwise the modeling results will be largely distorted.     

中国海域浮游硅藻钟形中鼓藻Bellerochea orologicalis Stosch初步研究

我国迄今已记录的中鼓藻属中只有一种——锤状中鼓藻Bellerochea malleus(Bright- well)Van Heurck有详细报导,作者在中国海域进行浮游生物取样时采到本属的另一种钟形中鼓藻Bellerochea horologicalis Stosch,1977.对本种与锤状中鼓藻的重要区别,以及本种的细胞形态、结构、生态习性与分布进行了描述。     

综合考虑多种变形因素的近岸规则波传播数学模型Ⅰ. 基本方程的导出

推广了Kirby的有环境水流影响的缓坡方程,得到了综合考虑环境水流(水流因子)、非线性弥散影响(非线性因子)、底摩擦波能损失(底摩擦因子)、非缓坡地形影响(地形因子)、折射、绕射、波浪破碎多种变形因素的波浪传播控制方程,并给出了非线性因子、地形因子、底摩擦因子、水流因子的确定方法。基于导出的方程做进一步推导,得到了波高和波向为变量的综合考虑多种变形因素的波浪传播基本方程,该方程有许多优点:1)其绕开了求解波势函数的困难,将椭圆型方程的边值问题化为初值问题;2)直接求解波高和波向;3)可采用有限差分法离散求解,对空间步长没有限制,适合大面积海区波场计算;4)综合考虑了多种波浪变形因素,方程更为合理,5)容易处理波浪破碎问题。     

结合抛物型缓坡方程计算波浪辐射应力

将波浪辐射应力与抛物型缓坡方程中的待求变量联系起来,提出了一种计算辐射应力的新方法,并用有限差分法对控制方程进行了数值求解。数值结果表明这种方法精度高、编程简单、求解快速,可用于实际大区域波浪辐射应力的计算。     

四叉树网格下的椭圆型缓坡方程数值模型研究

波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。     

Water Wave Refraction and Diffraction over Frictional Topography

- This paper considers the effects of wave energy dissipation induced by sea-bottom friction on the computational results of water wave refraction and diffraction with the parabolic equation method. The presented results show that the friction factor formula adopted in this paper is of higher numerical accuracy than that introduced by Dalrymphe (1984), and it can be used to compute wave propagation over large open areas.     

Large Eddy Simulation for Wave Breaking in the Surf Zone

In this paper, (he large eddy simulation method is used combined with the marker and cell method to study the wave propagation or shoaling and breaking process. As wave propagates into shallow water, the shoaling leads lo the increase of wave height, and then at a certain position, the wave will be breaking. The breaking wave is a powerful agent for generating turbulence, which plays an important role in most of the fluid dynamic processes throughout the surf zone, such as transformation of wave energy, generation of near-shore current and diffusion of materials. So a proper numerical model for describing the turbulence effect is needed. In this paper, a revised Smagorinsky subgrid-scale mode! is used to describe the turbulence effect. The present study reveals that the coefficient of the Smagorinsky model for wave propagation or breaking simulation may be taken as a varying function of the water depth and distance away from the wave breaking point. The large eddy simulation model presented in this pape     

Numerical study of vegetation damping effects on solitary wave run-up using the nonlinear shallow water equations

Vegetation damping effects on propagating water waves have been investigated by many researchers. This paper investigates the effects of damping due to vegetation on solitary water wave run-up via numerical simulation. The numerical model is based on an implementation of Morison's formulation for vegetation induced inertia and drag stresses in the nonlinear shallow water equations. The numerical model is solved via a finite volume method on a Cartesian cut cell mesh. The accuracy of the numerical scheme and the effects of the vegetation terms in the present model are validated by comparison with experiment results. The model is then applied to simulate a solitary wave propagating on a plane slope with vegetation. The sensitivity of solitary wave run-up to plant height, diameter and stem density is investigated by comparison of the numerical results for different patterns of vegetation. The numerical results show that vegetation can effectively reduce solitary wave propagation velocity and that solitary wave run-up is decreased with increase of plant height in water and also diameter and stem density.     

Wave-Current Propagation over a Frictional Topography

—In this paper the parabolic approximation model based on mild-slope equation is used tostudy wave propagation over a slowly varying and frictional topography under wave-current interaction.A governing equation considering the friction effects is derived by the authors for the first time.A simpli-fied form for the rate of wave energy dissipation is presented on the basis of the wave-current action conser-vation equation and the bottom friction model given by Yoo and O＇connor(1987).Examples reveal thatthe present computational method can be used for the calculation of wave elements for actual engineeringprojects with large water areas.     

Calculation of Wave Radiation Stress in Combination with Parabolic Mild Slope Equation

A new method for the calculation of wave radiation stress is proposed by linking the expres-sions for wave radiation stress with the variables in the parabolic mild slope equation.The governing equa-tions are solved numerically by the finite difference method.Numerical results show that the new methodis accurate enough,can be efficiently solved with little programming effort,and can be applied to the calcu-lation of wave radiation stress for large coastal areas.     

木本植被覆盖岸坡上波浪爬升过程的数值模拟研究

近岸木本植物构成的生态缓冲带作为新型的海岸软防护结构，兼具功能性和生态友好性，在沿海工程建设中愈发受到关注，如何深入开展其防护效果的机理研究是目前亟待解决的问题。本文采用数值模拟方法，在N-S方程中分别考虑树枝和树干的拖曳力影响，提出了木本植被作用下波浪沿斜坡爬升的表面波衰减的连续介质等效模型，并采用MAC法来跟踪自由曲面上的水颗粒轨迹。本文以波浪沿1/30的斜坡爬升为算例，对比讨论了有无植被作用下波浪的传播过程，并将算例结果与以往试验结果规律进行对照，验证了数值模型的有效性。最后，分别讨论了植物枝干的高度、密度、树枝倾斜角度等植被特性和波浪因素对植被消浪效果的影响，得到植被消浪的基本规律。文中的计算结果也可为实际的护岸工程和生态景观设计提供参考。     

A new form of generalized Boussinesq equations for varying water depth

A new set of equations of motion for wave propagation in water with varying depth is derived in this study. The equations expressed by the velocity potentials and the wave surface elevations include first-order non-linearity of waves and have the same dispersion characteristic to the extended Boussinesq equations. Compared to the extended Boussinesq equations, the equations have only two unknown scalars and do not contain spatial derivatives with an order higher than 2. The wave equations are solved by a finite element method. Fourth-order predictor–corrector method is applied in the time integration and a damping layer is applied at the open boundary for absorbing the outgoing waves. The model is applied to several examples of wave propagation in variable water depth. The computational results are compared with experimental data and other numerical results available in literature. The comparison demonstrates that the new form of the equations is capable of calculating wave transformation from relative deep water to shallow water.     

A numerical model for wave propagation in curvilinear coordinates

Using the perturbation method, a time dependent parabolic equation is developed based on the elliptic mild slope equation with dissipation term. With the time dependent parabolic equation employed as the governing equation, a numerical model for wave propagation including dissipation term in water of slowly varying topography is presented in curvilinear coordinates. In the model, the self-adaptive grid generation method is employed to generate a boundary-fitted and varying spacing mesh. The numerical tests show that the effects of dissipation term should be taken into account if the distance of wave propagation is large, and that the outgoing boundary conditions can be treated more effectively by introduction of the dissipation term into the numerical model. The numerical model is able to give good results of simulating wave propagation for waters of complicatedly boundaries and effectively predict physical processes of wave propagation. Moreover, the errors of the analytical solution deduced by Kirby et al. (1994) [Kirby, J.T., Dalrymple, R.A., Kabu, H., 1994. Parabolic approximation for water waves in conformal coordinate systems. Coastal Engineering 23, 185–213.] from the small-angle parabolic approximation of the mild-slope equation for the case of waves between diverging breakwaters in a polar coordinate system are corrected.     

Wash Waves Generated by High Speed Displacement Ships

It is difficult to compute far-field waves in a relative large area by using one wave generation model when a large calculation domain is needed because of large dimensions of the waterway and long distance of the required computing points. Variation of waterway bathymetry and nonlinearity in the far field cannot be included in a ship fixed process either. A coupled method combining a wave generation model and wave propagation model is then used in this paper to simulate the wash waves generated by the passing ship. A NURBS-based higher order panel method is adopted as the stationary wave generation model; a wave spectrum method and Boussinesq-type equation wave model are used as the wave propagation model for the constant water depth condition and variable water depth condition, respectively. The waves calculated by the NURBS-based higher order panel method in the near field are used as the input for the wave spectrum method and the Boussinesq-type equation wave model to obtain the far-field waves. With this approach it is possible to simulate the ship wash waves including the effects of water depth and waterway bathymetry. Parts of the calculated results are validated experimentally, and the agreement is demonstrated. The effects of ship wash waves on the moored ship are discussed by using a diffraction theory method. The results indicate that the prediction of the ship induced waves by coupling models is feasible.     

Modeling of propagation and transformation of transient nonlinear waves on a current

A novel theoretical approach is applied to predict the propagation and transformation of transient nonlinear waves on a current. The problem was solved by applying an eigenfunction expansion method and the derived semi-analytical solution was employed to study the transformation of wave profile and the evolution of wave spectrum arising from the nonlinear interactions of wave components in a wave train which may lead to the formation of very large waves. The results show that the propagation of wave trains is significantly affected by a current. A relatively small current may substantially affect wave train components and the wave train shape. This is observed for both opposing and following current. The results demonstrate that the application of the nonlinear model has a substantial effect on the shape of a wave spectrum. A train of originally linear and very narrow-banded waves changes its one-peak spectrum to a multi-peak one in a fairly short distance from an initial position. The discrepancies between the wave trains predicted by applying the linear and nonlinear models increase with the increasing wavelength and become significant in shallow water even for waves with low steepness. Laboratory experiments were conducted in a wave flume to verify theoretical results. The free-surface elevations recorded by a system of wave gauges are compared with the results provided by the nonlinear model. Additional verification was achieved by applying a Fourier analysis and comparing wave amplitude spectra obtained from theoretical results with experimental data. A reasonable agreement between theoretical results and experimental data is observed for both amplitudes and phases. The model predicts fairly well multi-peak spectra, including wave spectra with significant nonlinear wave components.     

A numerical model for combined refraction-diffraction of short waves on an island

A combined wave refraction-diffraction numerical model was developed to predict wave conditions around an arbitrary island. The methodology was based on the mild-slope equation, solved using a finite difference scheme with a marching procedure. The new model reduced the computer's memory demand considerably in comparison with finite-element, parabolic, error vector propagation and other finite difference approaches, and could therefore predict wave conditions for a large coastal area under given offshore boundary-wave conditions. Laboratory data on wave conditions under submerged circular and elliptical shoal conditions were selected to validate the numerical results. Good agreement was observed in all cases. Wave characteristics around an island were predicted using this model with the given deep-water wave condition. The model can predict wave conditions for any island with a mild-slope coastline.     

Parabolic Approximation of the Weakly Nonlinear Mild Slope Equation with Bottom Friction

This paper presents a refined parabolic approximation model of the mild slope equation to simu-late the combination of water wave refraction and diffraction in the large coastal region.The bottom frictionand weakly nonlinear term are included in the model.The difference equation is established with the Crank-Nicolson scheme.The numerical test shows that some numerical prediction results will be inaccurate in com-plicated topography without considering weak nonlinearity;the bottom friction will make wave height damp-ing and it can not be neglected for calculation of wave field in large areas.     

数值波浪水槽中完全非线性浅水波仿真特性研究

应用基于势流理论的时域高阶边界元方法,建立一个完全非线性的三维数值波浪水槽,通过实时模拟推板造波运动的方式产生波浪。通过混合欧拉-拉格朗日方法和四阶Runge-Kutta方法更新自由水面和造波板的瞬时位置。利用所建模型分别模拟了有限水深波和浅水波,与试验结果、相关文献结果和浅水理论结果吻合较好,且波浪能够稳定传播。系统地讨论造波板的运动圆频率、振幅和水深等对波浪传播和波浪特性的影响,并对波浪的非线性特性进行分析,研究发现造波板运动频率、运动振幅以及水深均将对波浪形态和波浪非线性产生显著影响。结果为真实水槽造波机的运动控制以及波浪生成试验提供了依据,便于实验室设置更合理的参数来准确模拟不同条件下的波浪。