孤立波在岸礁地形上传播与爬坡数值模拟研究

为准确模拟孤立波在岸礁地形上的传播和爬坡,采用基于完全非线性Boussinesq方程开发的Funwave-TVD模型,探究模型的可行性,并利用验证后的模型进一步研究岸礁各地形因素对孤立波爬高的影响。研究结果表明:模型能准确模拟孤立波在岸礁陡变地形上的传播及变形,摩擦系数对礁前陡坡及礁坪上的波浪传播模拟影响不大,但对爬坡预测的敏感性较强;模型空间步长可适当增大,提高计算效率;随着礁坪宽度的增大以及礁后斜坡的变缓,孤立波爬坡高度下降明显,而礁前陡坡坡度变化对孤立波爬坡高度影响不大。     

Boussinesq modelling of solitary wave and N-wave runup on coast

A total variation diminishing Lax–Wendroff scheme has been applied to numerically solve the Boussinesq-type equations. The runup processes on a vertical wall and on a uniform slope by various waves, including solitary waves, leading-depression N-waves and leading-elevation N-waves, have been investigated using the developed numerical model. The results agree well with the runup laws derived analytically by other researchers for non-breaking waves. The predictions with respect to breaking solitary waves generally follow the empirical runup relationship established from laboratory experiments, although some degree of over-prediction on the runup heights has been manifested. Such an over-prediction can be attributed to the exaggeration of the short waves in the front of the breaking waves. The study revealed that the leading-depression N-wave produced a higher runup than the solitary wave of the same amplitude, whereas the leading-elevation N-wave produced a slightly lower runup than the solitary wave of the same amplitude. For the runup on a vertical wall, this trend becomes prominent when the wave height-to-depth ratio exceeds 0.01. For the runup on a slope, this trend is prominent before the strong wave breaking occurs.     

珊瑚岛礁孤立波爬坡的平面二维数值模拟研究

为了探究岛屿周围珊瑚礁在抵御海啸灾害中的作用,采用激波捕捉类Boussinesq模型FUNWAVE-TVD,对孤立波在理想化三维岛礁地形上的传播及爬坡开展了现场尺度的平面二维数值模拟,分析了入射波高、礁坪水深、礁坪宽度、礁前斜坡坡度、礁后斜坡坡度、珊瑚礁糙率对岛屿四周孤立波爬高分布的影响。结果表明,珊瑚礁的存在总体上可有效降低岛屿四周孤立波的最大爬坡高度;入射波高、礁坪水深、礁坪宽度、珊瑚礁糙率是影响珊瑚岛礁四周孤立波爬坡分布的主要因素,岛礁四周最大爬坡高度会随入射波高和礁坪水深的增大、礁坪宽度和珊瑚礁糙率的减小而不断增大;当礁坪水深增大到一定程度时,珊瑚礁主要会对岛屿背浪面的爬高失去影响,而当礁坪宽度和珊瑚礁糙率减小至一定程度时,会出现岛礁四周最大爬高高于无珊瑚礁时爬高的现象;礁后斜坡的变缓会使岛礁周围的最大爬高有所减小,而礁前斜坡坡度对珊瑚岛礁周围的最大爬高几乎没有影响。     

Numerical modeling of wave runup and forces on an idealized beachfront house

In this paper, a numerical wave model based on the incompressible Reynolds-averaged Navier–Stokes (RANS) and k–ε equations is used to estimate the impact of a solitary wave on an idealized beachfront house located at different elevations on a plane beach. The locations of the free surface are reconstructed by volume of fluid (VOF) method. The model is satisfactorily tested against the experimental data of wave runup, and the analytical solution of wave forces on vertical walls. The time histories of wave profiles, forces, and overturning moments on the idealized house are demonstrated and analyzed. The variations of wave forces and overturning moments with the elevation of the idealized beachfront house are also investigated.     

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.     

Characteristics of the boundary layer at the bottom of a solitary wave

The results of direct numerical simulations of the boundary layer generated at the bottom of a solitary wave are described. The numerical results, which agree with the laboratory measurements of Sumer et al. (2010) show that the flow regime in the boundary layer can be laminar, laminar with coherent vortices and turbulent. The average velocity and bottom shear stress are computed and the results obtained show that the logarithmic law can approximate the velocity profile only in a restricted range of the parameters and at particular phases of the wave cycle. Moreover, the maximum value of the bottom shear stress is found to depend on the dimensionless wave height only, while the minimum (negative) value depends also on the dimensionless boundary layer thickness. Diagrams and simple formulae are proposed to evaluate the minimum and maximum bottom shear stresses and their phase shift with respect to the wave crest.     

可渗透海床单桩孤立波破碎波浪力数值研究

当波浪传播至近岸浅水区时易发生破碎,波浪破碎后水体直接拍击单桩结构,其波浪作用力显著增大,可能导致结构失稳破坏。首先建立包括斜坡海床的数值波浪水槽,并与已有研究进行对比验证。进而开展考虑斜坡海床可渗透性的孤立波数值模拟,分析孤立波传播与浅水化破碎特征,着重研究竖直单桩上破碎波浪力的特性,及其与单桩位置、海床渗透率的关联性。数值研究发现:对于低渗透率海床,随着单桩位置由深水向岸线位置变动,其所受波浪力先增后减,所受波浪力最大的桩体位置随海床渗透率增加而从波浪破碎点前方移动至破碎点后1D处,而在高渗透率海床上不同位置处桩体所受波浪力均较小;随海床渗透率等梯度增加,海床消波作用逐渐增强,波浪破碎进程延缓,波浪破碎点向岸线方向加速移动,单桩上破碎波浪力呈整体下降趋势,但可能因波浪破碎点的位置变动导致部分位置桩体所受波浪力异常增大。     

Computation of solitary waves during propagation and runup on a slope

A numerical time-simulation algorithm for analysing highly nonlinear solitary waves interacting with plane gentle and steep slopes is described by employing a mixed Eulerian–Lagrangian method. The full nonlinear free surface conditions are considered here in a Lagrangian frame of reference without any analytical approximations, and thus the method is valid for very steep waves including overturning. It is found that the runup height is crucially dependent on the wave steepness and the slope of the plane. Pressures and forces exerted on impermeable walls of different inclinations (slopes) by progressive shallow water solitary waves are studied. Strong nonlinear features in the form of pronounced double peaks are visible in the time history of pressure and force signals with increasing heights of the oncoming solitary waves. The effect of nonlinearity is less pronounced as the inclination of the wall decreases with respect to the bottom surface.     

A modeling study of acoustic propagation through movingshallow-water solitary wave packets

Propagation of 400-Hz sound through continental-shelf internal solitary wave packets is shown by numerical simulation to be strongly influenced by coupling of normal modes. Coupling in a packet is controlled by the mode coefficients at the point where sound enters the packet, the dimensions of the waves and packet, and the ambient depth structures of temperature and salinity. In the case of a moving packet, changes of phases of the incident modes with respect to each other dominate over the other factors, altering the coupling over time and thus inducing signal fluctuations. The phasing within a moving packet varies with time scales of minutes, causing coupling and signal fluctuations with comparable time scales. The directionality of energy flux between high-order acoustic modes and (less attenuated) low-order modes determines a gain factor for long-range propagation. A significant finding is that energy flux toward low-order modes through the effect of a packet near a source favoring high-order modes will give net amplification at distant ranges. Conversely, a packet far from a source sends energy into otherwise quiet higher modes. The intermittency of the coupling and of high-mode attenuation via bottom interaction means that signal energy fluctuations and modal diversity fluctuations at a distant receiver are complementary, with energy fluctuations suggesting a source-region packet and mode fluctuations suggesting a receiver-region packet. Simulations entailing 33-km propagation are used in the analyses, imitating the SWARM experiment geometry, allowing comparison with observations     

Wind wave breaking under the conditions of inhomogeneous currents and of the atmospheric boundary layer

As known fromin situ observations, inhomogeneities of flows and of the atmospheric boundary layer produce variations of the intensity of wind wave breaking. A relevant phenomenological model is suggested here, describingin situ data on the breaking of waves in the presence of internal waves. The response of the wave breaking to the flow's inhomogeneity enhances with the growth of its spatial or temporal scale. For the mesoscale (10–100 km) inhomogeneities, the model is essentially simplified—wave breakings depict the local energy inputs to wind waves. The model allows us to compute currents of various type in the wave breaking intensity field. The results may have practical implications, in terms of remote sensing of the ocean. Translated by Vladimir A. Puchkin.     

Boundary integral equation solutions for solitary wave generation, propagation and run-up

The boundary integral equation method (BIEM) is developed as a tool for studying two-dimensional, nonlinear water wave problems, including the phenomena of wave generation, propagation and run-up. The wave motions are described by a potential flow theory. Nonlinear free-surface boundary conditions are incorporated in the numerical formulation. Examples are given for either a solitary wave or two successive solitary waves. Special treatment is developed to trace the run-up and run-down along a shoreline. The accuracy of the present scheme is verified by comparing numerical results with experimental data of maximum run-up.     

海堤上波浪爬高与越浪计算分析

基于线性波理论,采用VOF法、RNG k-ε的湍流模型,通过FLUENT软件平台的二次开发功能实现了质量源造波。利用多孔介质模型,在动量方程里添加多孔介质源项,对斜坡堤上铺设异型块体时的爬高与越浪进行了数值模拟,并通过模拟爬高值率定出FLUENT中多孔阻力系数与多种异形块体糙渗系数之间的关系,进而建立了可供实际应用的海堤上波浪爬高与越浪数值模型。     

Effect of surface wave breaking on the surface boundary layer of temperature in the Yellow Sea in summer

Near-surface enhancement of turbulent mixing and vertical mixing coefficient for temperature owing to the effect of surface wave breaking is investigated using a two-dimensional (2-D) ocean circulation model with a tidal boundary condition in an idealized shelf sea. On the basis of the 2-D simulation, the effect of surface wave breaking on surface boundary layer deepening in the Yellow Sea in summer is studied utilizing a 3-D ocean circulation model. A well-mixed temperature surface layer in the Yellow Sea can be successfully reconstructed when the effect of surface wave breaking is considered. The diagnostic analysis of the turbulent kinetic energy equation shows that turbulent mixing is enhanced greatly in the Yellow Sea in summer by surface wave breaking. In addition, the diagnostic analysis of momentum budget and temperature budget also show that surface wave breaking has an evident contribution to the turbulent mixing in the surface boundary layer. We therefore conclude that surface wave breaking is an important factor in determining the depth of the surface boundary layer of temperature in the Yellow Sea in summer.     

Numerical modeling on the interaction of internal solitary wave with slope-shelf and modal analysis

On the basis of a nonhydrostatic numerical model,the interaction of internal solitary wave with slope-shelf was studied.The breaking and polarity transformation were analyzed.A "kink" structure,due to shoaling topography and higher nonlinear effect,was found to be generated by the leading wave before breaking.Coherent vortex shedding behind the leading wave was presented.The evolution characteristics of the modal structure were analyzed based on the empirical orthogonal function method.The modal structure was complicated due to the effect of the variable topography,especially when breaking occurred.In the performed experiments,the contributions to the total variance from higher mode jumped from no more than 20% to over 40%.     

Reflection of large-scale inhomogeneities in the air-sea boundary layer on the intensity of surface wave breaking in the ocean

An overview of experiments is given on the observation of the dependence of the intensity of wind wave breaking on large-scale currents' inhomogeneities and atmospheric boundary layer stratification. The data were interpreted using a model in which the variance of the wind wave breaking intensity reflects fluctuations in the influx of energy to wind waves due to various factors. Translated by V. Puchkin.     

Validation of a new generation system for bottom boundary layer beneath solitary wave

Understanding of sea bottom boundary layer characteristics, especially bottom shear stress acting on the sea bed, is an important step needed in sediment transport modeling for practical application purposes. In the present study, a new generation system for bottom boundary layer under solitary wave is proposed. Applicability of this system is examined by comparing measured and numerical solution velocities. Moreover, transitional behavior from laminar to turbulence was investigated. It is concluded that the critical Reynolds number in the experiments shows good agreement with DNS result of Vittori and Blondeaux (2008) and laboratory data of Sumer et al. (2010), indicating validity of the generation system. Since the present generation system enables continuous measurement to obtain ensemble averaged quantities, it can be effectively utilized for future experimental studies on solitary wave boundary layers, including sediment transport experiments with movable bed.     

Numerical simulation of wave transformation and runup incorporating porous media wave absorber and turbulence models

A numerical wave tank is first established using the Navier–Stokes equations and the VOF method assuming laminar flow. The standard k–ε, realizable k–ε and RNG k–ε turbulent models are then incorporated to the numerical tank. An effective numerical method for wave absorption utilizing the energy-dissipating property of porous media is also included. To validate the accuracy of the proposed models, the propagation of a solitary wave, where analytical solution is available for comparison, is first simulated. This is followed by the simulation of irregular wave runup on a composite seawall, wave propagation over submerged bars and wave refraction and diffraction over an elliptic shoal, where experimental data are available for comparison. All computed results agree well with either the analytical solution or the experimental data.     

The role of bottom friction in models of nonbreaking tsunami wave runup on the shore

The role of bottom friction in the runup of nonbreaking long waves on the shore is analyzed. The case of the normal incidence of monochromatic waves is considered. The relief of the model region consists of an even horizontal bottom area conjugated with a flat slope. The energy dissipation is estimated as the work of bottom friction forces over the wave field obtained using the known analytical solution based on the Carrier-Greenspan transforms. Estimates of energy losses for waves whose periods are typical for tsunami waves have been obtained. The energy dissipation is shown to be not concentrated in the shore line area as a rule. The question about the practicability of using partially reflecting boundary conditions on the coast to take into account the bottom friction in large-scale models of tsunami propagation is considered.     

Drift velocity at the wave breaking point

Wave-tank studies were conducted on the measurement of the drift velocity at the breaking point under different types of breaking waves on a rigid, plane beach. The drift velocity has onshore direction near the surface and close to the bottom; in the main flow column, the drift velocity is always offshore. The offshore drift velocity shows a more uniform vertical distribution than that in the offshore region. The experimental data are compared with theoretical values of three different second-order constant-depth wave theories. Comparisons with data from other sources are also made.     

Intermittency of wind wave group breaking

Knowledge on intermittency of wave breaking is so far limited to a few summary statistics, while the probability distribution of time interval between breaking events can provide a full view of intermittency. Based on a series of experiments on wind wave breaking, such probability distributions are investigated. Breaking waves within a wave group were taken as a single breaking event according to recent studies. Interval between successive wave groups with breaker is the focus of this paper. For intervals in our experiments with different fetch and wind conditions, their distributions are all skewed and weighted on small intervals. Results of Kolmogorov-Smirnov tests on time series of these intervals indicate that they all follow gamma distribution, and some are even exponential type. Average breaking-group-interval decreases with friction velocity and significant steepness until the wind is strong enough;most of them are more than 10 times the dominant wave period. Group breaking probability proposed by Babanin recently and the average number of breaking waves in wave groups are also discussed, and they are seemingly more reasonable and sensitive than traditional breaking probability defined in terms of single wave.