近岸波浪和波生流环境中泥沙输运的数值模拟

Owing to lack of observational data and accurate definition,it is difficult to distinguish the Kuroshio intrusion water from the Pacific Ocean into the South China Sea(SCS).By using a passive tracer to identify the Kuroshio water based on an observation-validated three-dimensional numerical model MITgcm,the spatio-temporal variation of the Kuroshio intrusion water into the SCS has been investigated.Our result shows the Kuroshio intrusion is of distinct seasonal variation in both horizontal and vertical directions.In winter,the intruding Kuroshio water reaches the farthest,almost occupying the area from 18°N to 23°N and 114°E to 121°E,with a small branch flowing towards the Taiwan Strait.The intrusion region of the Kuroshio water decreases with depth gradually.However,in summer,the Kuroshio water is confined to the east of 118°E without any branch reaching the Taiwan Strait;meanwhile the intrusion region of the Kuroshio water increases from the surface to the depth about 205 m,then it decreases with depth.The estimated annual mean of Kuroshio Intrusion Transport(KIT) via the Luzon Strait is westward to the SCS in an amount of –3.86×106 m3/s,which is larger than the annual mean of Luzon Strait Transport(LST) of –3.15×106 m3/s.The KIT above 250 m accounts for 60%–80% of the LST throughout the entire water column.By analyzing interannual variation of the Kuroshio intrusion from the year 2003 to 2012,we find that the Kuroshio branch flowing into the Taiwan Strait is the weaker in winter of La Ni?a years than those in El Ni?o and normal years,which may be attributed to the wind stress curl off the southeast China then.Furthermore,the KIT correlates the Ni?o 3.4 index from 2003 to 2012 with a correlation coefficient of 0.41,which is lower than that of the LST with the Ni?o 3.4 index,i.e.,0.78.     

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.     

Numerical Study of Influence of Currents on Nonlinear Focusing Waves

In this paper,a numerical model for nonlinear wave propagation in currents is formulated by a set of enhanced fully nonlinear Boussinesq equations with ambient currents.This model is verified by comparison with the published results.Then the influence of currents on nonlinear focusing waves is studied by use of the numerical model.It is found that the effect of currents on the surface elevations at the focal location is negligible.Following currents can augment the maximum crest of focusing wave while decre...     

Numerical study on water waves and wave-induced longshore currents in Obaky coastal water

In this paper, the water waves and wave-induced longshore currents in Obaky coastal water which is located at the Mediterranean coast of Turkey were numerically studied. The numerical model is based on the parabolic mild-slope equation for coastal water waves and the nonlinear shallow water equation for the wave-induced currents. The wave transformation under the effects of shoaling, refraction, diffraction and breaking is considered, and the wave provides radiation stresses for driving currents in the model. The numerical results for the water wave-induced longshore currents were validated by the measured data to demonstrate the efficiency of the numerical model. Then the water waves and longshore currents induced by the waves from main directions were numerically simulated and analyzed based on the numerical results. The numerical results show that the movement of the longshore currents was different while the wave propagated to a coastal zone from different directions.     

Numerical Simulation of Spatial Lag Between Wave Breaking Point and Location of Maximum Wave-Induced Current

A quasi three-dimensional numerical model of wave-driven coastal currents with the effects of surface rollers is developed for the study of the spatial lag between the location of the maximum wave-induced current and the wave breaking point.The governing equations are derived from Navier-Stokes equations and solved by the hybrid method combining the fractional step finite different method in the horizontal plane with a Galerkin finite element method in the vertical direction.The surface rollers effects are considered through incorporating the creation and evolution of the roller area into the free surface shear stress.An energy equation facilitates the computation process which transfers the wave breaking energy dissipation to the surface roller energy.The wave driver model is a phase-averaged wave model based on the wave action balance equation.Two sets of laboratory experiments producing breaking waves that generated longshore currents on a planar beach are used to evaluate the model's performance.The present wave-driven coastal current model with the roller effect in the surface shear stress term can produce satisfactory results by increasing the wave-induced nearshore current velocity inside the surf zone and shifting the location of the maximum longshore current velocity landward.     

Numerical Simulation of the Three-Dimensional Wave-Induced Currents on Unstructured Grid

By coupling the three-dimensional hydrodynamic model with the wave model, numerical simulations of the three-dimensional wave-induced current are carried out in this study. The wave model is based on the numerical solution of the modified wave action equation and eikonal equation, which can describe the wave refraction and diffraction. The hydrodynamic model is driven by the wave-induced radiation stresses and affected by the wave turbulence. The numerical implementation of the module has used the finite-volume schemes on unstructured grid, which provides great flexibility for modeling the waves and currents in the complex actual nearshore, and ensures the conservation of energy propagation. The applicability of the proposed model is evaluated in calculating the cases of wave set-up, longshore currents, undertow on a sloping beach, rip currents and meandering longshore currents on a tri-cuspate beach. The results indicate that it is necessary to introduce the depth-dependent radiation stresses into the numerical simulation of wave-induced currents, and comparisons show that the present model makes better prediction on the wave procedure as well as both horizontal and vertical structures in the wave-induced current field.     

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 morphodynamic modelling of tidal sand waves on the shoreface

An artificial sand wave on the Dutch shoreface of the North Sea has been studied in conditions with relatively strong tidal currents in the range of 0.5 to 1 m/s and sediments in the medium sand size range of 0.2 to 0.5 mm. The sand wave is perpendicular to the tidal current and has a maximum height and length of the order of 5 m and 1 km, respectively. The sand wave is dynamically active and shows migration rates of the order of a few metres per year. A numerical morphodynamic model (DELFT3D model) has been used to simulate the morphological behaviour of the sand wave in the North Sea. This model approach is based on the numerical solution of the three-dimensional shallow water equations in combination with a surface wave propagation model (wind waves) and the advection–diffusion equation for the sediment particles with online bed updating after each time step. The model results show that the sand wave grows in the case of dominant bed-load transport (weak tidal currents; relatively coarse sediment; small roughness height; low waves) and that the sand wave decays in the case of dominant suspended transport (strong currents, relatively fine sediment, large roughness height; storm waves).     

Coupled modeling of currents and wind waves in the Kerch Strait

We present a numerical model of the dynamics of the Kerch Strait allowing one to perform the coordinated analysis of the fields of currents and wind waves. The model includes the spectral wave module and the hydrodynamic block of currents. The influence of waves on the currents is taken into account in the hydrodynamic block both via the surface and bottom tangential stresses and via the radiation stresses. In order to take into account the inverse influence of currents upon the waves, we use the fields of currents and sea level from the hydrodynamic block in the wave module. The specific features of the structure of currents and wind waves in the strait are studied for the typical wave situations. The results of the coupled and separate simulation are compared and the importance of taking into account the mechanisms of interaction between waves and currents in the analysis of the dynamic processes in the strait is demonstrated. __________ Translated from Morskoi Gidrofizicheskii Zhurnal, No. 5, pp. 3–20, September–October, 2007.     

Numerical simulation of wave-induced currents combined with parabolic mild-slope equation in curvilinear coordinates

Researches on breaking-induced currents by waves are summarized firstly in this paper. Then, a combined numerical model in orthogonal curvilinear coordinates is presented to simulate wave-induced current in areas with curved boundary or irregular coastline. The proposed wave-induced current model includes a nearshore current module established through orthogonal curvilinear transformation form of shallow water equations and a wave module based on the curvilinear parabolic approximation wave equation. The wave module actually serves as the driving force to provide the current module with required radiation stresses. The Crank-Nicolson finite difference scheme and the alternating directions implicit method are used to solve the wave and current module, respectively. The established surf zone currents model is validated by two numerical experiments about longshore currents and rip currents in basins with rip channel and breakwater. The numerical results are compared with the measured data and published numerical results.     

Two-dimensional wave-current interaction with a locally enriched quadtree grid system

Wave-induced currents may drive nearshore mixing and transport processes, including coastal pollutant dispersion, littoral drift, and long-term morphological changes through beach erosion and accretion. In this study, a numerical model is newly developed to simulate wave climate and localized currents in complicated coastal environments. The model developed is based on a quadtree grid system. The two-dimensional hydrodynamic governing equations are solved by using an explicit Adams-Bashforth finite difference scheme. Effects of wave breaking, shoaling, refraction, diffraction, wave-current interaction, set-up and set-down, turbulent mixing, bed friction, and shoreline movement are incorporated in the model. Results are presented for set-up at a beach in a flume due to normally incident waves, and longshore currents generated by oblique waves on a plane beach.     

渤海波浪和潮汐风暴潮相互作用对波浪影响的数值研究

基于依赖波浪成长状态波令的表面风应力，提出了一个波浪和风暴潮潮汐运动相互作用的联合数值模式，实现了第三代波浪模式和三维风暴潮潮汐模式联合作用的数值研究，并结合渤海典型天气个例的研究，给出了渤海波浪和风暴潮潮汐相互作用对波浪影响的机制和大小量级的定量估计。研究表明，对不同天气过程，波浪和风暴潮潮汐相互作用对波浪影响的性质和大小不同；对强寒潮过程，对波浪影响主要由风暴潮所支配波高调制可达1m，在黄河口区一般达0．5m；对弱天气过程，对波浪影响主要由潮所控制，波高调制约在0．2m，联合作用模式给出的结果与实测更吻合。     

Tests and Applications of An Approach to Absorbing Reflected Waves Towards Incident Boundary

If the upstream boundary conditions are prescribed based on the incident wave only, the time-dependent numerical models cannot effectively simulate the wave field when the physical or spurious reflected waves become significant. This paper describes carefully an approach to specifying the incident wave boundary conditions combined with a set sponge layer to absorb the reflected waves towards the incident boundary. Incorporated into a time-dependent numerical model, whose governing equations are the Boussinesq-type ones, the effectiveness of the approach is studied in detail. The general boundary conditions, describing the down-wave boundary conditions are also generalized to the case of random waves. The numerical model is in detail examined. The test cases include both the normal one-dimensional incident regular or random waves and the two-dimensional oblique incident regular waves. The calculated results show that the present approach is effective on damping the reflected waves towards the incident wave boundary.     

A coupled numerical model of wave interaction with porous medium

A new coupling model of wave interaction with porous medium is established in which the wave field solver is based on the two dimensional Reynolds Averaged Navier-Stokes (RANS) equations with a kε closure. Incident waves, which could be linear waves, cnoidal waves or solitary waves, are produced by a piston-type wave maker in the computational domain and the free surface is traced through the Piecewise Linear Interface Construction-Volume of Fluid (PLIC-VOF) method. Nonlinear Forchheimer equations are adopted to calculate the flow field within the porous media. By introducing a velocity–pressure correction equation, the wave field and the porous flow field are highly and efficiently coupled. The two fields are solved simultaneously and no boundary condition is needed at the interface of the internal porous flow and the external wave. The newly developed numerical model is used to simulate wave interaction with porous seabed and the numerical results agree well with the experimental data. The additional numerical tests are also conducted to study the effects of seabed thickness, porosity and permeability coefficient on wave damping and the pore water pressure responses.     

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

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

Numerical modeling of nonlinear water waves over heterogeneous porous beds

The transformation of nonlinear water waves over porous beds is studied by applying a numerical model based on Chen's [2006. Fully nonlinear Boussinesq-type equations for waves and currents over porous beds. Journal of Engineering Mechanics, 132:2, 220–230] Boussinesq-type equations for highly nonlinear waves on permeable beds. The numerical model uses a high-order time-marching solution and fourth-order finite-difference schemes for discretization of first-order spatial derivatives to obtain a computational accuracy consistent with the model equations. By forcing the wave celerity and spatial porous-damping rate of the linearized model to match the exact linear theory for horizontal porous bed over a prescribed range of relative depths, the values of the model parameters are optimally determined. Numerical simulations of the damped wave propagation over finite-thickness porous layer demonstrate the accuracy of both the numerical model and governing equations, which have been shown by prior theoretical analyses to be accurate for both nominal and thick porous layers. These simulations also elucidate on the significance of the higher-order porous-damping terms and the influence of the hydraulic parameters. Application of the model to the simulation of the wave field around a laboratory-scale submerged porous mound provides a measure of its capability, as well as useful insight into the scaling of the porous-resistance coefficients. For application to heterogeneous porous beds, the assumption of weak spatial variation of the porous resistance is examined using truncated forms of the governing equations. The results indicate that the complete set of Boussinesq-type equations is applicable to porous beds of nonhomogeneous makeup.     

Evolution of waves and currents over a submerged laboratory shoal

The vertically-integrated effect of interaction between waves and wave-induced currents on wave transformation over a submerged elliptic shoal was investigated based on numerical simulations of the Vincent and Briggs experiment [Vincent, C.L., Briggs, M.J., 1989. Refraction- diffraction of irregular waves over a mound. Journal of Waterway, Port, Coastal and Ocean Engineering, 115(2), pp. 269–284.]. The numerical simulations were performed using two numerical wave-current model systems: one, a combination of the wave model SWAN and the current model SHORECIRC, and the other, a combination of the wave model REF/DIF and the same current model. A time-dependent, phase-resolving wave and current model, FUNWAVE, was also utilized to simulate the experiment. In the simulations, the developed wave-induced currents defocused waves behind the shoal and brought on a wave shadow zone that showed relatively low wave height distributions. For the breaking case of monochromatic waves, the wave heights computed using FUNWAVE showed good agreement with the measurements and the resulting wave-induced currents showed a jet-like velocity distribution in transverse direction. And the computed results of the two model combinations agreed better with the measurements than the computed results obtained by neglecting wave-current interaction. However, it was found that for the case in which transverse interference pattern caused by refracted waves was strong, REF/DIF-SHORECIRC did not correctly evaluate radiation stresses, the gradients of which generate wave-induced currents. SWAN-SHORECIRC, which cannot deal with the interference patterns, predicted a jet-like wave-induced current. For breaking random wave cases, the computed results of the two model combinations and FUNWAVE agreed well with the measurements. The agreements indicate that it is necessary to take into account the effect of wave-induced current on wave refraction when wave breaking occurs over a submerged shoal.     

New model to determine forces at on-bottom slender pipelines

The present paper proposes a numerical model to determine horizontal and vertical components of the hydrodynamic forces on a slender submarine pipeline lying at the sea bed and exposed to non-linear waves plus a current. The new model is an extension of the Wake II type model, originally proposed for sinusoidal waves (Soedigdo et al., 1999) and for combined sinusoidal waves and currents (Sabag et al., 2000), to the case of periodic or random waves, even with a superimposed current. The Wake II type model takes into account the wake effects on the kinematic field and the time variation of drag and lift hydrodynamic coefficients. The proposed extension is based on an evolutional analysis carried out for each half period of the free stream horizontal velocity at the pipeline. An analytical expression of the wake velocity is developed starting from the Navier–Stokes and the boundary layer equations. The time variation of the drag and lift hydrodynamic coefficients is obtained using a Gaussian integration of the start-up function. A reduced scale laboratory investigation in a large wave flume has been conducted in order to calibrate the empirical parameters involved in the proposed model. Different wave and current conditions have been considered and measurements of free stream horizontal velocities and dynamic pressures on a bottom-mounted pipeline have been conducted. The comparison between experimental and numerical hydrodynamic forces shows the accuracy of the new model in evaluating the time variation of peaks and phase shifts of the horizontal and vertical wave and current induced forces.     

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.     

Wave Numerical Model for Shallow Water

The history of forecasting wind waves by wave energy conservation equation is briefly des-cribed.Several currently used wave numerical models for shallow water based on different wave theoriesare discussed.Wave energy conservation models for the simulation of shallow water waves are introduced,with emphasis placed on the SWAN model,which takes use of the most advanced wave research achieve-ments and has been applied to several theoretical and field conditions.The characteristics and applicabilityof the model,the finite difference numerical scheme of the action balance equation and its source termscomputing methods are described in detail.The model has been verified with the propagation refractionnumerical experiments for waves propagating in following and opposing currents;finally.the model is ap-plied to the Haian Gulf area to simulate the wave height and wave period field there,and the results arecompared with observed data.