非线性波传播的新型数值模拟模型及其实验验证

以一种新型的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并通过采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.通过对均匀水深水域内波浪传播的数值模拟说明,模型能较好地模拟大水深水域和强非线性波的传播.通过设置不同的入射波参数来进行潜堤地形上波浪传播的物理模型实验,并将数值解与物理模型实验结果进行了比较.     

强非线性和色散性Boussinesq方程数值模型检验

采用同位网格有限差分法,建立了强非线性和色散性Boussinesq方程数值计算模型。以稳恒波Fourier近似解给定入射波边界条件,对均匀水深深水和浅水域不同非线性的行进波、缓坡地形上深水至浅水域的浅水变形波、以及缓坡和陡坡地形上的波浪水槽实验进行了数值计算,并将计算结果与解析解、解析数值解以及实验值进行了较为详细的比较,从而检验了模型的色散性、非线性以及不同底坡下非线性波的浅水变形性能。     

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

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

Numerical simulation of nonlinear dispersive waves propagating over a submerged bar by IB–VOF model

It is a good test for a numerical model to simulate progressive waves propagating over a submerged bar with a relatively high ratio of slopes. In this paper, the combined IB–VOF model is used to predict nonlinear dispersive waves propagating over a submerged bar with both slopes of 1:2. The predicted free surface elevations are compared with the experimental data and numerical results presented by other researchers. The comparison shows that the IB–VOF model is able to provide satisfactory wave profiles in the shallow water with strong nonlinear effects and in the wave transmitted region with strong wave dispersion in particular. Moreover, the wave evolution behind the submerged bar is described in detail, including the spatial wave profile modulation, spectral analysis of the time-series waves, flow velocity and pressure fields, and kinetic energy distribution. The effect of fluid viscosity on the numerical simulations is also studied, and it is found that the effect on the wave evolution considered in this paper is not significant. Finally, the hydrodynamic force acting on the bar is calculated using the IB–VOF model.     

非线性波传播的新型数值模拟模型及其实验验证引入变换速度变量

以一种新型的含变换速度变量的Boussinesq型方程为控制方程组,采用五阶Runge-Kutta-England格式离散时间积分,采用七点差分格式离散空间导数,并采用恰当的出流边界条件,从而建立了非线性波传播的新型数值模拟模型.对均匀水深水域内波浪传播的数值模拟,说明在引入变换速度后进一步增大了模型的水深适用范围.对潜堤地形上波浪传播的数值模拟说明,在引入变换速度后进一步提高了模型的数值模拟精度.     

A spectral model for unidirectional nonlinear wave propagation over arbitrary depths

A weakly-nonlinear and dispersive wave equation recently developed by the authors is used for formulating a spectral-type unidirectional wave propagation model describing spectral transformations of narrow-band waves travelling over arbitrary depths. The essential characteristics of the model equation are recapitulated first and then the spectral domain representation in terms of spatially varying harmonic amplitudes is presented. The resulting evolution equations are used to simulate the experiments concerning harmonic generation in shallow water and nonlinear random wave transformations over a submerged bar. Furthermore, the spectral model predictions are compared with the field measurements in nearshore with satisfactory results.     

Wave Dispersion Study in the Indian Ocean-Tsunami of December 26, 2004

A numerical study which takes into account wave dispersion effects has been carried out in the Indian Ocean to reproduce the initial stage of wave propagation of the tsunami event that occurred on December 26, 2004. Three different numerical models have been used: the nonlinear shallow water (nondispersive), the nonlinear Boussinesq, and the full Navier-Stokes aided by the volume of fluid method to track the free surface. Numerical model results are compared against each other. General features of the wave propagation agreed very well in all numerical studies. However some important differences are observed in the wave patterns, i.e., the development in time of the wave front is shown to be strongly connected to the dispersion effects. Discussions and conclusions are made about the spatial and temporal distribution of the free surface reaffirming that the dispersion mechanism is important for tsunami hazard mitigation.     

On the modeling of wave propagation on non-uniform currents and depth

By transforming two different time-dependent hyperbolic mild slope equations with dissipation term for wave propagation on non-uniform currents into wave-action conservation equation and eikonal equation, respectively, shown are the different effects of dissipation term on the eikonal equation in the two different mild slope equations. The performances of intrinsic frequency and wave number are also discussed. Thus the suitable mathematical model is chosen in which the wave number vector and intrinsic frequency are expressed both more rigorously and completely. By using the perturbation method, an extended evolution equation, which is of time-dependent parabolic type, is developed from the time-dependent hyperbolic mild slope equation which exists in the suitable mathematical model, and solved by using the alternating direction implicit (ADI) method. Presented is the numerical model for wave propagation and transformation on non-uniform currents in water of slowly varying topography. From the comparisons of the numerical solutions with the theoretical solutions of two examples of wave propagation, respectively, the results show that the numerical solutions are in good agreement with the exact ones. Calculating the interactions between incident wave and current on a sloping beach [Arthur, R.S., 1950. Refraction of shallow water waves. The combined effects of currents and underwater topography. EOS Transactions, August 31, 549–552], the differences of wave number vector between refraction and combined refraction–diffraction of waves are discussed quantitatively, while the effects of different methods of calculating wave number vector on numerical results are shown.     

Wave Dispersion Study in the Indian Ocean-Tsunami of December 26, 2004

A numerical study which takes into account wave dispersion effects has been carried out in the Indian Ocean to reproduce the initial stage of wave propagation of the tsunami event that occurred on December 26, 2004. Three different numerical models have been used: the nonlinear shallow water (nondispersive), the nonlinear Boussinesq, and the full Navier-Stokes aided by the volume of fluid method to track the free surface. Numerical model results are compared against each other. General features of the wave propagation agreed very well in all numerical studies. However some important differences are observed in the wave patterns, i.e., the development in time of the wave front is shown to be strongly connected to the dispersion effects. Discussions and conclusions are made about the spatial and temporal distribution of the free surface reaffirming that the dispersion mechanism is important for tsunami hazard mitigation.     

Three-dimensional nonlinear random wave groups in intermediate water depth

High waves at ocean occur during a complex space–time evolution of wave groups. In this paper the nonlinear structure of three-dimensional sea wave groups at intermediate water depth is investigated. To this purpose, the Boccotti's Quasi-Determinism theory is firstly applied to describe the linear wave groups when a given exceptionally high crest occurs. Then, the second-order correction to the linear solution is derived for the general condition of three-dimensional wave groups, at a finite water depth. Several numerical applications, finally, have been carried out in order to show how both the spectral bandwidth and the directional spreading modify the nonlinear high waves at different water depth.     

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.     

SWAN模型对实验室小尺度浅水波浪模拟研究

基于浅水斜坡地形的物理模型试验数据,考察SWAN模型对实验室小尺度浅水波浪的模拟效果,进而检验其浅水项的模拟精度。模拟中采用直接输入初始测点的实测海浪谱进行造波,重点考察浅水中三波相互作用和变浅破碎两个源项,对不同工况下,SWAN模式在水深条件变化下的有效波高、谱平均周期、海浪谱演化的模拟能力进行研究。研究表明:模拟的有效波高较符合实测波浪的增长和衰减,但谱平均周期计算值明显偏小;海浪谱的能量转移机制同实测有较大区别,频谱模拟结果出现高频高估、低频低估现象。对两个源项进行对比分析得出三波相互作用对海浪谱的能量转换影响远大于变浅破碎耗散。想要提高近岸区谱平均周期和海浪谱的模拟精度则SWAN模型中三波非线性项的计算精确度仍需更多研究和改进。     

Statistical Analyses of Wave Height Distribution for Multidirectional Irregular Waves over A Sloping Bottom

The main objective of this paper is to examine the influences of both the principal wave direction and the directional spreading parameter of the wave energy on the wave height evolution of multidirectional irregular waves over an impermeable sloping bottom and to propose an improved wave height distribution model based on an existing classical formula. The numerical model FUNWAVE 2.0, based on a fully nonlinear Boussinesq equation, is employed to simulate the propagation of multidirectional irregular waves over the sloping bottom. Comparisons of wave heights derived from wave trains with various principal wave directions and different directional spreading parameters are conducted. Results show that both the principal wave direction and the wave directional spread have significant influences on the wave height evolution on a varying coastal topography. The shoaling effect for the wave height is obviously weakened with the increase of the principal wave direction and with the decrease of the directional spreading parameter. With the simulated data, the classical Klopman wave height distribution model is improved by considering the influences of both factors. It is found that the improved model performs better in describing the wave height distribution for the multidirectional irregular waves in shallow water.     

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.     

Temporal and azimuthal dependence of sound propagation in shallowwater with internal waves

The short time scale (minutes) and azimuthal dependence of sound wave propagation in shallow water regions due to internal waves is examined. Results from the shallow water acoustics in random media (SWARM-95) experiment are presented that reflect these dependencies. Time-dependent internal waves are modeled using the dnoidal solution to the nonlinear internal wave equations, so that the effects of both temporal and spatial variability can be assessed. A full wave parabolic equation model is used to simulate broadband acoustic propagation. It is shown that the short term temporal variability and the azimuthal dependence of the sound field are strongly correlated to the internal wave field     

Nonlinear saturation-based whitecapping dissipation in SWAN for deep and shallow water

This study investigates the effectiveness of a revised whitecapping source term in the spectral wind wave model SWAN (Simulating WAves Nearshore) that is local in frequency space, nonlinear with respect to the variance density and weakly dependent on the wave age. It is investigated whether this alternative whitecapping expression is able to correct the tendency towards underprediction of period measures that has been identified in the default SWAN model. This whitecapping expression is combined with an alternative wind input source term that is more accurate for young waves than the default expression. The shallow water source terms of bottom friction, depth-induced breaking and triad interaction are left unaltered. It is demonstrated that this alternative source term combination yields improved agreement with fetch- and depth-limited growth curves. Moreover, it is shown, by means of a field case over a shelf sea, that the investigated model corrects the erroneous overprediction of wind-sea energy displayed by the default model under combined swell-sea conditions. For a selection of field cases recorded at two shallow lakes, the investigated model generally improves the agreement with observed spectra and integral parameters. The improvement is most notable in the prediction of period measures.     

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.     

The applicability of the shallow water equations for modelling violent wave overtopping

The shallow water equations (SWE) have been used to model a series of experiments examining violent wave overtopping of a near-vertical sloping structure with impacting wave conditions. A finite volume scheme was used to solve the shallow water equations. A monotonic reconstruction method was applied to eliminate spurious oscillations and ensure proper treatment of bed slope terms. Both the numerical results and physical observations of the water surface closely followed the relevant Rayleigh probability distributions. However, the numerical model overestimated the wave heights and suffered from the lack of dispersion within the shallow water equations. Comparisons made on dimensionless parameters for the overtopping discharge and percentage of waves overtopping between the numerical model and the experimental observations indicated that for the lesser impacting waves, the shallow water equations perform satisfactorily and provide a good alternative to computationally more expensive methods.     

Numerical simulation of wave overtopping at coastal dikes and low-crested structures by means of a shock-capturing Boussinesq model

In this paper, a shock-capturing numerical model, based on the combined solution of Boussinesq and nonlinear shallow water equations is applied to the simulation of wave runup, overtopping and wave train propagation over impermeable, emerged and low-crested, structures. In order to improve the performances of the scheme at wet–dry interfaces, a numerical treatment is introduced to provide well-balancing of advective fluxes and source terms. One- and two-dimensional test cases are presented to validate the performances of the model under both regular and irregular wave conditions.     

Modeling of nonlinear wave propagation over fringing reefs

The applicability of existing nonlinear (triad) spectral models for steep slopes (0.1–0.2) characteristic of reef environments was investigated, using both deterministic (phase-resolving) and stochastic (phased-averaged) formulations. Model performance was tested using laboratory observations of unidirectional wave transformation over steep and smooth bathymetry profiles. The models, developed for mild slopes, were implemented with minimal modifications (the inclusion of breaking parametrizations and linear steep-slope corrections) required by laboratory data. The deterministic model produced typically more accurate predictions than the stochastic one, but the phase averaged formulation proved fast enough to allow for an inverse modeling search for the optimal breaking parametrization. The effects of the additional assumptions of the stochastic approach resulted in a slower than observed evolution of the infragravity band. Despite the challenge posed by the fast wave evolution and energetic breaking characteristic to the steep reef slopes, both formulations performed overall well, and should be considered as good provisional candidates for use in numerical investigation of wave–current interaction processes on steep reefs.