Run-up of solitary waves on slopes with different profiles

The problem of sea-wave run-up on a beach is discussed within the framework of exact solutions of a nonlinear theory of shallow water. Previously, the run-up of solitary waves with different forms (Gaussian and Lorentzian pulses, a soliton, special-form pulses) has already been considered in the literature within the framework of the same theory. Depending on the form of the incident wave, different formulas were obtained for the height of wave run-up on a beach. A new point of this study is the proof of the universality of the formula for the maximum height of run-up of a solitary wave on a beach for the corresponding physical choice of the determining parameters of the incident wave, so that the effect of difference in form is eliminated. As a result, an analytical formula suitable for applications, in particular, in problems related to tsunamis, has been proposed for the height of run-up of a solitary wave on a beach.     

Efficient computation of surf zone waves using the nonlinear shallow water equations with non-hydrostatic pressure

A numerical method for non-hydrostatic, free-surface, irrotational flow governed by the nonlinear shallow water equations including the effects of vertical acceleration is presented at the aim of studying surf zone phenomena. A vertical boundary-fitted grid is used with the water depth divided into a number of layers. A compact finite difference scheme is employed for accurate computation of frequency dispersion requiring a limited vertical resolution and hence, capable of predicting the onset of wave breaking. A novel wet–dry algorithm is applied for a proper handling of moving shoreline. Mass and momentum are strictly conserved at discrete level while the method only dissipates energy in the case of wave breaking. The numerical results are verified with a number of tests and show that the proposed model using two layers without ad-hoc assumptions enables to resolve propagating nonlinear shoaling, breaking waves and wave run-up within the surf and swash zones in an efficient manner.     

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.     

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.     

The influence of seaward boundary conditions on swash zone hydrodynamics

The influence of the seaward boundary condition on the internal swash hydrodynamics is investigated. New numerical solutions of the characteristics form of the nonlinear shallow-water equations are presented and applied to describe the swash hydrodynamics forced by breaking wave run-up on a plane beach. The solutions depend on the specification of characteristic variables on the seaward boundary of the swash zone, equivalent to prescribing the flow depth or the flow velocity. It is shown that the analytical solution of Shen and Meyer [Shen, M.C., Meyer, R.E., 1963. Climb of a bore on a beach. Part 3. Runup. J. Fluid Mech. 16, 113–125] is a special case of the many possible solutions that can describe the swash flow, but one that does not appear appropriate for practical application for real waves. The physical significance of the boundary conditions is shown by writing the volume and momentum fluxes in terms of the characteristic variables. Results are presented that illustrate the dependence of internal flow depth and velocity on the boundary condition. This implies that the internal swash hydrodynamics depend on the shape and wavelength of the incident bore, which differs from the hydrodynamic similarity inherent in the analytical solution. A solution appropriate for long bores is compared to laboratory data to illustrate the difference from the analytical solution. The results are important in terms of determining overwash flows, flow forces and sediment dynamics in the run-up zone.     

Maximum run-up,breaking conditions and dynamical forces in the swash zone: a boundary value approach

On the basis of the approximate analytical solution for the nonlinear shallow water equations of Antuono and Brocchini [M. Antuono & M. Brocchini, The boundary value problem for the nonlinear shallow water equation, Stud. Appl. Maths, 119, 71–91 (2007).], we propose useful regression curves for the prediction of maximum run-up and dynamical forces in the swash zone on a frictionless, uniformly sloping beach. For the first time the dependence of the results on both the wave height and the wave steepness is analyzed in detail providing formulae able to describe a wide class of wave inputs. Finally, the regression formulae are validated through comparison with maximum run-up laws and breaking conditions already available in the literature, the present model results appearing to better account for nonlinear effects.     

基于InterDyMFoam的潜礁孤立波传播数值模拟

建立了基于OpenFOAM动边界的类fixedValuePointPatchVectorField继承的仿物理造波数值波浪水槽。对孤立波在三种不同潜礁地形上的行进、爬坡以及破碎等典型过程进行数值模拟,模拟结果与Boussinesq模型及物理模型试验所得结果进行了对比分析。结果表明,采用动边界进行仿物理造波更适合处理波高水深比较大的孤立波传播问题,可以较好地模拟孤立波在潜礁上传播引起的波浪破碎、水跃等现象。     

Tsunami generation,propagation, and run-up with a high-order Boussinesq model

In this work we extend a high-order Boussinesq-type (finite difference) model, capable of simulating waves out to wavenumber times depth kh < 25, to include a moving sea-bed, for the simulation of earthquake- and landslide-induced tsunamis. The extension is straight forward, requiring only an additional term within the kinematic bottom condition. As first test cases we simulate linear and nonlinear surface waves generated from both positive and negative impulsive bottom movements. The computed results compare well against earlier theoretical, numerical, and experimental values. Additionally, we show that the long-time (fully nonlinear) evolution of waves resulting from an upthrusted bottom can eventually result in true solitary waves, consistent with theoretical predictions. It is stressed, however, that the nonlinearity used far exceeds that typical of geophysical tsunamis in the open ocean. The Boussinesq-type model is then used to simulate numerous tsunami-type events generated from submerged landslides, in both one and two horizontal dimensions. The results again compare well against previous experiments and/or numerical simulations. The new extension compliments recently developed run-up capabilities within this approach, and as demonstrated, the model can therefore treat tsunami events from their initial generation, through their later propagation, and final run-up phases. The developed model is shown to maintain reasonable computational efficiency, and is therefore attractive for the simulation of such events, especially in cases where dispersion is important.     

片状分布植被对孤立波爬高影响的数值模拟研究

近岸植被对波浪爬坡具有一定的衰减作用。在自然界中，由于植被的死亡、再生或人为破坏等原因，近岸植被通常呈片状分布，且其内部分布也是不均匀的。本文以完全非线性Boussinesq方程为基础，引入植被作用项，建立了模拟近岸植被区波浪传播的数值模型，验证了模型可靠性，进而采用该模型模拟分析了片状分布植被对孤立波爬高的影响。数值模拟结果表明，片状分布植被能有效减小孤立波爬高；对于均匀分布的片状植被，高密度片状植被对孤立波爬高的消减效果优于低密度片状植被；对于相同密度、不同分布形式的片状植被，均匀分布的片状植被对孤立波的消减效果优于不均匀分布的片状植被；对于不均匀分布的片状植被，前密后疏的片状植被对孤立波的消减效果优于前疏后密的片状植被。     

Run-up of long solitary waves of different polarities on a plane beach

We study the run-up of long solitary waves of different polarities on a beach in the case of composite bottom topography: a plane sloping beach transforms into a region of constant depth. We confirm that nonlinear wave deformation of positive polarity (wave crest) resulting in an increase in the wave steepness leads to a significant increase in the run-up height. It is shown that nonlinear effects are most strongly pronounced for the run-up of a wave with negative polarity (wave trough). In the latter case, the run-up height of such waves increases with their steepness and can exceed the amplitude of the incident wave.     

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.     

Numerical investigation of tsunami-like wave hydrodynamic characteristics and its comparison with solitary wave

Solitary waves have been commonly used as an initial condition in the experimental and numerical modelling of tsunamis for decades. However, the main component of a tsunami waves acts at completely different spatial and temporal scales than solitary waves. Thus, use of solitary waves as approximation of a tsunami wave may not yield realistic model results, especially in the coastal region where the shoaling effect restrains the development of the tsunami wave. Alternatively, N-shaped waves may be used to give a more realistic approximation of the tsunami wave profile. Based on the superposition of the sech²(*) waves, the observed tsunami wave profile could be approximated with the N-shaped wave method, and this paper presents numerical simulation results based on the tsunami-like wave generated based on the observed tsunami wave profile measured in the Tohoku tsunami. This tsunami-like wave was numerically generated with an internal wave source method based on the two-phase incompressible flow model with a Volume of Fluid (VOF) method to capture the free surface, and a finite volume scheme was used to solve all the governing equations. The model is first validated for the case of a solitary wave propagating within a straight channel, by comparing its analytical solutions to model results. Further, model comparisons between the solitary and tsunami-like wave are then made for (a) the simulation of wave run-up on shore and (b) wave transport over breakwater. Comparisons show that use of these largely different waveform shapes as inputs produces significant differences in overall wave evolution, hydrodynamic load characteristics as well as velocity and vortex fields. Further, it was found that the solitary wave uses underestimated the total energy and hence underestimated the run-up distance.     

Three-dimensional numerical simulation of solitary wave run-up using the IB method

Although the finite difference method is computationally efficient, it is acknowledged to be inferior when dealing with flow-over on structures with a complex geometry because of its rectilinear grid system. Therefore, we developed a numerical procedure that can cope with flow over structures with complex shapes while, at the same time, retaining the simplicity and efficiency of a rectilinear grid system. We used the immersed boundary method, which involves application of immersed boundary forces at solid boundaries rather than conventional boundary conditions, to investigate wave interactions with coastal structures in a three-dimensional numerical wave tank by solving the Navier–Stokes equations for two-phase flows. We simulated the run-up of a solitary wave around a circular island. Maximum run-up heights were computed around the island and compared with available laboratory measurements and previous numerical results. The three-dimensional features of the run-up process were analyzed in detail and compared with those of depth-integrated equations models.     

孤立波与带窄缝双箱相互作用模拟研究

针对孤立波与带窄缝双箱的作用问题,应用时域高阶边界元方法建立了二维数值水槽。其中,自由水面满足完全非线性运动学和动力学边界条件,对瞬时自由表面流体质点采用混合欧拉-拉格朗日法追踪,采用四阶龙格库塔法对下一时刻的自由水面的速度势和波面升高进行更新。采用加速度势法求解物体湿表面的瞬时波浪力。采用推板方法生成孤立波。通过模拟孤立波在直墙上的爬高以及施加在直墙上的波浪力,并与已发表的实验和数值结果对比,验证本数值模型的准确性。通过数值模拟计算研究了窄缝宽度、方箱尺寸对波浪在箱体迎浪侧爬高,窄缝内波面升高,箱体背浪侧透射波高及箱体受波浪荷载的影响。同时研究了有一定时间间隔的双孤立波与带窄缝双箱系统作用问题。     

Run-up of tsunamis and long waves in terms of surf-similarity

In this paper we review and re-examine the classical analytical solutions for run-up of periodic long waves on an infinitely long slope as well as on a finite slope attached to a flat bottom. Both cases provide simple expressions for the maximum run-up and the associated flow velocity in terms of the surf-similarity parameter and the amplitude to depth ratio determined at some offshore location. We use the analytical expressions to analyze the impact of tsunamis on beaches and relate the discussion to the recent Indian Ocean tsunami from December 26, 2004. An important conclusion is that extreme run-up combined with extreme flow velocities occurs for surf-similarity parameters of the order 3–6, and for typical tsunami wave periods this requires relatively mild beach slopes. Next, we compare the theoretical solutions to measured run-up of breaking and non-breaking irregular waves on steep impermeable slopes. For the non-breaking waves, the theoretical curves turn out to be superior to state-of-the-art empirical estimates. Finally, we compare the theoretical solutions with numerical results obtained with a high-order Boussinesq-type method, and generally obtain an excellent agreement.     

A Lagrangian model for boundary layer growth and bed shear stress in the swash zone

A new model for the boundary layer development and associated skin friction coefficients and shear stress within the swash zone is presented. The model is developed within a Lagrangian reference frame, following fluid trajectories, and can be applied to both laminar flow and smooth turbulent flow. The model is based on the momentum integral approach for steady, flat-plate boundary layers, with appropriate modifications to account for the unsteady flow regime and flow history. The model results are consistent with previous measurements of bed shear stress and skin friction coefficients within the swash zone. These indicate strong temporal and spatial variation throughout the swash cycle, and a clear distinction between the uprush and backwash phase. This variation has been previously attributed the unsteady flow regime and flow history effects, both of which are accounted for in the new model. Fluid particle trajectories and velocity are computed using the non-linear shallow water wave equations and the boundary layer growth across the entire swash zone is estimated. Predictions of the bed shear stress and skin friction coefficients agree reasonably well with direct bed shear stress measurements reported by Barnes et al. (Barnes, M.P., O’Donaghue, T., Alsina, J.M., Baldock, T.E., 2009. Direct bed shear stress measurements in bore-driven swash. Coastal Engineering 56 (8), 853–867) and, for a given flow velocity, give stresses which are consistent with the bias toward uprush sediment transport which has consistently been observed in measurements. The data and modelling suggest that the backwash boundary layer is initially laminar, which results in the late development of significant bed shear during the backwash, with a transition to a turbulent boundary layer later in the backwash. A new conceptual model for the boundary layer structure at the leading edge of the swash is proposed, which accounts for both the no-slip condition at the bed and the moving wet–dry interface. However, further development of the Lagrangian Boundary Layer Model is required in order to include bore-generated turbulence and to account for variable roughness and mobile beds.     

Modelling Solitary Waves and Its Impact on Coastal Houses with SPH Method

The interaction between solid structures and free-surface flows is investigated in this study. A Smoothed Particle Hydrodynamics (SPH) model is used in the investigation and is verified against analytical solutions and experimental observations. The main aim is to examine the effectiveness of a tsunami-resistant house design by predicting the wave loads on it. To achieve this, the solitary wave generation and run-up are studied first. The solitary wave is generated by allowing a heavily weighted block to penetrate into a tank of water at one end, and the near-shore seabed is modelled by an inclined section with a constant slope. Then, the SPH model is applied to simulate the three-dimensional flows around different types of houses under the action of a solitary wave. It has been found that the tsunami-resistant house design reduces the impact force by a factor of three.     

Moving shoreline boundary condition for nearshore models

The paper develops and analyzes two fully nonlinear boundary conditions that incorporate the motion of the shoreline in nonlinear time domain nearshore models. A moving shoreline essentially means the computational domain is changing with the solution of the flow. The problem is solved in two steps. The first is to establish an equation that determines the motion of the shoreline based on the local momentum balance. The second is to develop and implement into a shoreline model the capability of accommodating a changing computational domain. The two models represent two different ways of addressing this step: one is to track the position of the shoreline in a fixed grid by establishing a special shoreline point which generally is not a fixed grid point. The second is by a coordinate transformation that maps the changing domain onto a fixed domain and solves the basic equations in the mapped domain. The two shoreline conditions are tested against three known solution for nonlinear shoreline motion. Two are the 1-D solutions to the nonlinear shallow water (NSW) equations by Carrier and Greenspan [J. Fluid Mech. 4 (1958) 97], one representing the response to a transient change in the offshore water level, the other the motion due to a periodic standing wave, both on slopes steep enough to allow full reflection. The third is the 2-D horizontal (2DH) computational solution by Zelt [Coast. Eng. 15 (1991) 205] for the run-up of a solitary wave on a cusped beach. In all cases, both models are shown to behave well and give high accuracy results for suitably chosen grid and time spacings.     

基于Boussinesq水波模型的聚焦波模拟

基于最高导数为3阶的单层Boussinesq方程,建立了聚焦波的时域波浪计算模型。数值模型求解采用了预报?校正的有限差分法。对于时间差分格式,预报和校正分别采用3阶Adams-Bashforth格式和4阶Adams-Moulton格式。首先,针对不同水深条件下水槽中传播的强非线性波进行模拟,并将数值结果与流函数的数值解析解进行了比较,结果表明无论是波面位移、波面处的水平速度和垂向速度均与解析解符合较好,最大波峰面的速度分布伴随水深的增加与解析解吻合程度变差,非线性速度分布的适用范围与线性解析解适应范围kh<3.5基本一致。其次,对深水聚焦波演化进行了模拟研究,研究中聚焦波的生成采用在边界点累加不同频率线性规则波的方法。应用聚焦波物理模型实验结果验证模型,计算聚焦位置处的波面位移和沿水深的速度分布与实验结果的对比表明,波面位移吻合程度较好,垂向的水平速度分布基本吻合。最后,保持中心频率(周期)不变,数值模拟了周期范围变化下最大聚焦波峰面以及波峰面水平速度的变化趋势,结果表明波峰面值和波峰面水平速度随着周期范围缩小而增大。     

On the run-up and back-wash processes of single and double solitary waves — An experimental study

The run-up and back-wash processes of single and double solitary waves on a slope were studied experimentally. Experiments were conducted in three different wave flumes with four different slopes. For single solitary wave, new experimental data were acquired and, based on the theoretical breaking criterion, a new surf parameter specifically for breaking solitary waves was proposed. An equation to estimate maximum fractional run-up height on a given slope was also proposed. For double solitary waves, new experiments were performed by using two successive solitary waves with equal wave heights; these waves were separated by various durations. The run-up heights of the second wave were found to vary with respect to the separation time. Particle image velocimetry measurements revealed that the intensity of the back-wash flow generated by the first wave strongly affected the run-up height of the second wave. Showing trends similar to that of the second wave run-up heights, both the back-wash breaking process of the first wave and the reflected waves were strongly affected by the wave–wave interaction. Empirical run-up formula for the second solitary wave was also introduced.