波浪作用下近岸区污染物输移扩散的数学模型及其实验验证

在近岸缓坡浅水海岸,波浪破碎产生沿岸流是近岸海域流场的重要组成部分,它对污染物输移扩散规律的影响重大,在高阶近似抛物化缓坡方程求解大面积波浪场基础上,建立了波浪作用下污染物输移扩散数学模型.计算结果与不同坡度均匀斜坡地形上具有不同波高、周期的规则波及不规则波浪作用下污染物输移扩散实验结果进行了比较,分析了各种因素对波浪作用下沿岸流分布规律影响,所得结论认为地形坡度及入射波高对污染物输移扩散的影响较大,波浪作用将使缓坡海滩上污染物的输移扩散平行岸线方向.     

结合椭圆型缓坡方程模拟近岸波流场

波浪向近岸传播的过程中,由波浪破碎效应所产生的近岸波流场是近岸海域关键的水动力学因素之一.结合近岸波浪场的椭圆型缓坡方程和近岸波流场数学模型对近岸波浪场及由斜向入射波浪破碎后所形成的近岸波流场进行了数值模拟.计算中考虑到波浪向近岸传播中由于波浪的折射、绕射、反射等效应使局部复杂区域波向不易确定,采用结合椭圆型缓坡方程所给出的波浪辐射应力公式来计算波浪产生的辐射应力,在此基础上耦合椭圆型缓坡方程和近岸波流场数学模型对近岸波流场进行数值模拟,从而使模型综合考虑了波浪的折射、绕射、反射等效应且避免了对波向角的直接求解,可以应用于相对较复杂区域的近岸波流场模拟.     

非结构化网格下近岸波生流数值模拟

波浪破碎产生的近岸流是近岸海域关键的水动力因素之一。基于近岸波浪的椭圆型缓坡方程和二维近岸波生流方程,建立了非结构化网格下近岸波浪破碎形成的近岸流数值模型。数值模型中,在空间上采用有限体积法进行数值离散,在时间上采用欧拉向前格式数值离散。数值计算结果表明,该数值模型可以有效地模拟近岸波浪破碎产生的近岸流。     

近岸岛屿周围波生流的数值模拟

基于非结构网格有限体积法的环流模型,建立了相位平均的波流作用模型。模型中考虑了波浪破碎效应的三维辐射应力和波流作用下的扩散系数。通过模拟试验水池中人工岛和天然岛屿周围的波生流,对岛周围波生流的运动特征进行了分析。当波浪垂直入射到均匀斜坡上的人工岛时,岛周围产生回流和对称的环流,流速沿水深分布稍有差异。在天然岛屿中,入射波浪在近岸形成沿岸流。靠近岸线岛屿周围的沿岸流比离岸较远岛屿周围的流速大。这和辐射应力在浅水区域作用有关,岸线区域波生流特征明显。近岸岛屿周围在波浪辐射应力下产生不同的波浪增水,该波浪增水形成压力梯度进而促进岛屿周围的流态变化。     

基于抛物型缓坡方程模拟近岸植被区波浪传播

植被对波浪传播运动有重要影响。考虑近岸波浪在植被区传播中的折射、绕射、破碎及植被引起的波能耗损效应,基于抛物型缓坡方程建立了模拟近岸植被区波浪传播的数学模型,对模型进行了数值模拟验证,采用数值模拟试验分析了植被对波浪传播的影响。数值模拟结果表明,波浪在近岸植被区传播时,随着植被密度和植被高度的增加,波浪传播中的波高衰减增大,波能耗损增加;不同周期波浪在植被区传播中的波高衰减过程也明显不同。     

近岸波浪破碎区不规则波浪的数值模拟

基于近岸不规则波浪传播的抛物型缓坡方程和两类波浪破碎能量损耗因子,对近岸波浪破碎区不规则波浪的波高分布进行了数值模拟,并结合实验结果对数值模拟结果进行了验证分析,结果表明采用两类波浪破碎能量损耗因子所模拟的破碎区波高与实测值均吻合良好,波浪破碎能量损耗因子及波浪破碎指标对破碎区波浪波高分布影响较明显。     

波生沿岸流的数值模拟研究

采用完全非线性Boussinesq方程建立的FUNWAVE模型进行波生沿岸流数值模拟研究,通过对不同斜坡地形和波浪入射条件下波生流的物理模型实验结果进行比较,验证了该数值模型能较准确地计算沿岸流;通过改变波浪的不同入射条件,对不同入射条件的沿岸流数值模拟得出:当其他条件不变时,仅入射波高增大时,沿岸流的流幅和幅值增大,幅值位置向深水移动;仅增大入射周期时,沿岸流的流幅显著增加,幅值的增加较小;斜坡地形坡度的改变能显著影响波生沿岸流的流幅和幅值,但对沿岸流幅值位置的水深影响不大。采用窄缝法处理动边界时,选择合适的窄缝起始水深对沿岸流的准确计算是十分必要的。     

圆形数值波浪水池斜向和多向入射波浪模拟研究

为了研究斜向入射波浪,基于三维不可压缩两相流模型,开发了一套圆形数值波浪水池数值模型。在圆形波浪水池中,通过源项造波法成功生成了任意入射方向的波浪,并且利用人工摩擦项模拟阻尼区以数值耗散反射波浪。模型基于嵌入式多块网格体系,采用FVM法(finite volume method)离散Navier-Stokes方程,VOF法(volume of fluid)追踪自由水面。试验结果表明,斜向入射波浪的模拟结果与理论值基本一致,圆形波浪水池在模拟斜向入射波浪时,有效区域的面积较传统波浪水池显著增大,而且有效区域受波浪入射角度的影响也较小。同时,通过叠加多列斜向入射波浪,模拟出了多向交叉波列,并通过与理论结果对比,发现其具有较高的精度。     

斜坡平台波浪破碎数值模拟

波浪在斜坡地形上破碎,破波后稳定波高多采用物理模型试验方法进行研究,利用近岸波浪传播变形的抛物型缓坡方程和波能流平衡方程,导出了适用于斜坡上波浪破碎的数值模拟方法。首先根据波能流平衡方程和缓坡方程基本型式分析波浪在破波带内的波能变化和衰减率,推导了波浪传播模型中波能衰减因子和破波能量流衰减因子之间的关系;其次,利用陡坡地形上的高阶抛物型缓坡方程建立了波浪传播和波浪破碎数学模型;最后,根据物理模型试验实测数据对数值模拟的效果进行验证。数值计算与试验资料比较表明,该模型可以较好地模拟斜坡地形的波浪传播波高变化。     

波浪增减水的实用数学模型及其数值模拟

提出了一种模拟近岸区波浪增减水的实用数学模型.首先采用考虑波能损失的抛物型缓坡方程数值模拟波浪破碎引起的波浪复振幅变化,接着根据计算得到的波浪复振幅,采用一种新的辐射应力公式计算辐射应力分量,然后采用深度平均方程计算波浪破碎产生的增减水.采用该模型对规则波和不规则波破碎引起的增减水进行了数值模拟,数值模拟结果与实验结果吻合良好,表明该模型可有效模拟近岸区由于波浪破碎引起的增减水.     

Numerical study of pollutant movement in waves and wave-induced long-shore currents in surf zone

Water waves, wave-induced long-shore currents and movement of pollutants in waves and currents have been numerically studied based on the hyperbolic mild-slope equation, the shallow water equation , as well as the pollutant movement equation, and the numerical results have also been validated by experimental data. It is shown that the long-shore current velocity and wave set-up increase with the increasing incident wave amplitude and slope steepness of the shore plane ; the wave set-up increases with the in- creasing incident wave period;and the pollutant morement proceeds more quiekly with the increasing incident wave amplitude and slope steepness of the shore palane. In surf zones, the long-shore currents induced by the inclined incident waves have effectively affected the pollutant movement.     

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.     

Rip currents under obliquely incident wind waves and tidal longshore currents

A field experiment on the nature of rip currents was conducted on the Dutch coast, which differs from previous rip current study sites because it is a wind-sea dominated environment with mostly obliquely incident waves and tidally-driven longshore currents. During the experiment three distinct flow patterns, obtained with GPS tracked drifter instruments, were observed: (1) a locally governed circulation cell, (2) an offshore current that is deflected shore parallel outside the surf zone and (3) a meandering longshore current. The transition from rip currents (flow patterns 1 and 2) to meandering longshore currents (flow pattern 3) occurred gradually within the tidal cycle with longshore currents prevalent at mid to high tide. Rip currents at this site appeared at depressions in the surf zone bar and typically occurred when the water level fell below NAP (equivalent to MSL), even in the presence of obliquely incident waves and tidally driven longshore currents. Hindcast simulations of the drifter experiments were performed with the numerical model XBeach and showed good agreement with field observations. The model was subsequently used to investigate the influence of tidal water level fluctuations, longshore currents and obliquely incident waves on rip currents.Rip currents were initiated when the water level dropped below a specific threshold with the magnitude of the rip current associated with the water level. The strength of the tidal current and its orientation with respect to the incident waves governed the offshore extent and orientation of the rip current. In contrast to other studies that suggest that rip currents solely occur under shore normal (or slightly oblique waves), in this study both observations and numerical model simulations indicate that rip currents can exist under large angles of wave incidence, when the rip channel is sufficiently wide and the wave height is small.     

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.     

Numerical study of wave and longshore current interaction

Wave and longshore current interaction was examined based on the numerical models.In these models,water waves in the presence of longshore currents were modeled by parabolic mild slope equation,and wave breaking induced longshore currents were modeled by shallow water equation.Water wave provided the radiation stress gradients to drive current.Wave and longshore current interactions were considered by cycling the wave and longshore current models to a steady state.The experiments for regular and irregular breaking wave induced longshore currents by Hamilton and Ebersole（2001） and Reniers and Battjes（1997） were simulated.The numerical results indicate that the present models are effective for simulating the interaction of wave and breaking wave induced longshore currents,and the numerically simulated longshore current at wave breaking point considering wave and longshore current interaction show some disagreement with those neglecting the wave-current interaction,and the breaking wave induced longshore current effect on wave transformation is not obvious.     

Longshore currents of regular waves on different beaches

The experiment and numerical computations of longshore currents produced by regularwaves on the two beaches with the slopes of 1:100 and 1:40 are made. The cross-shore distributions of longshore current velocities and wave heights are given and the influences of wave heights, wave periods and beach slopes on the longshore currents are discussed. The discussion is also made for the influences of different eddy viscosity coefficients on the numerical results of longshore current velocities.     

波生流对海岸污染物输移的影响

通过物理模型实验对海岸波浪作用下污染物运动特性进行了分析,重点分析了质量输移流、沿岸流、沿岸流不稳定运动及破波带内旋涡运动等海岸水动力因素的影响.实验中坡度分别取为1:100和1:40,实验中采用CCD摄像机记录墨水的运动轨迹,同步测量流体质点速度以及波面升高.实验表明,在破碎带外污染物主要受波浪非线性引起的质量输移流的影响;在破碎带内主要受沿岸流的影响,同时还受沿岸流不稳定运动及大尺度旋涡运动的影响.     

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.