河口海岸泥沙输移的分组数学模型研究

根据河口海岸水沙输移的特点,建立了一个新的二维分组数学模型,用来预测该区域的水沙输移过程。该模型耦合了水动力模块、泥沙输移模块和床面演变模块。其中水动力模块基于浅水方程组,综合考虑了柯氏力、床面切应力以及表面风应力的影响,引入干湿判断法处理动边界。泥沙输移模块首先将泥沙按照粒径分组,针对不同泥砂性质,对各组泥沙分别进行建模求解。床面演变模块基于质量守恒方程,实时更新床面高程以及床沙级配变化,并传递给水动力模块,更新底部边界。该模型被应用在了英国塞汶(Severn)河口,其预测的泥沙浓度和实测数据以及不分组的模型的预测结果进行了比较,结果显示,文中建立的分组模型预测的结果要明显好于不分组模型。     

非对称歪斜波引起的片流输沙数值研究

本文用了一个可考虑相位差作用和波浪边界层非对称性的瞬态理论模型和一个两相紊流模型共同研究非对称歪斜波引起的片流输沙现象。为了解速度偏度和加速度偏度对输沙通量和输沙率的贡献,两相流模型为理论模型提供了必要的相位超前、瞬时侵蚀深度和边界层的发展过程。理论模型研究显示了由速度偏度和加速度偏度引起的向岸阶段和离岸阶段的泥沙运动非对称性,解释了净输沙的产生原因。在以往的非对称歪斜波片流输沙研究中,净输沙的产生主要被归结于相位差作用。本文的研究则表明了非对称的边界层发展所产生的净流量和动床面效应在净输沙产生过程中的比相位差作用更为重要。     

后台阶下游床面冲刷过程数值模拟研究

以床面瞬时剪应力作为泥沙起动及输运的水动力机制,建立了结构物周围复杂流场下床面局部冲刷的数学模型。并应用大涡模拟方法对后台阶下游三维湍流流动进行数值模拟,得到台阶下游床面瞬时剪应力的分布规律。为了确定床面瞬时剪应力与泥沙上扬通量的关系,先应用数学模型对不同模型参数下,冲刷开始后5分钟时台阶下游床面形状进行试算。通过试算与实验结果的比较,确定床面瞬时剪应力与泥沙上扬通量关系中需要的模型参数。进一步对冲刷开始后30分钟内台阶下游床面演化规律进行模拟,模拟结果与实验结果相符合。     

Simulation of the Fate of Faecal Bacteria in Estuarine and Coastal Waters Based on A Fractionated Sediment Transport Model

A two-dimensional depth-integrated numerical model is refined in this paper to simulate the hydrodynamics, graded sediment transport process and the fate of faecal bacteria in estuarine and coastal waters. The sediment mixture is divided into several fractions according to the grain size. A bed evolution model is adopted to simulate the processes of the bed elevation change and sediment grain size sorting. The faecal bacteria transport equation includes enhanced source and sink terms to represent bacterial kinetic transformation and disappearance or reappearance due to sediment deposition or re-suspension. A novel partition ratio and dynamic decay rates of faecal bacteria are adopted in the numerical model. The model has been applied to the turbid water environment in the Bristol Channel and Severn estuary, UK. The predictions by the present model are compared with field data and those by non-fractionated model.     

水东湾地貌形态演变数值模拟

为了研究泻湖型海湾内经常出现的湾中岛的形成机理,应用水平二维潮流、泥沙输移和地形演变耦合模型,对具有典型的沙坝泻湖地貌形态的水东湾的湾中岛的形成和演化进行了数值模拟,成功地模拟出了涨潮三角洲(大洲岛)的形成过程,计算结果与实际地貌形态总体符合。模型中分别考虑了全沙输沙和推移质输沙两种输沙情况。结果表明,水东湾湾中岛是由涨潮流引起的泻湖内泥沙不断淤积而形成的,是一种涨潮三角洲的地貌形态;全沙输沙模式比推移质输沙模式更适合这一地形演化过程的模拟。模拟结果也再现了湾中岛在25 a期间的地形变化过程,这一研究结果为通过计算数值模拟来研究一般海湾的地貌形态的形成机理和演化过程提供了实际算例。     

High-resolution field measurements and numerical modelling of intra-wave sediment suspension on plane beds under shoaling waves

Intra-wave sediment suspension is examined using high-resolution field measurements and numerical hydrodynamic and sediment models within 120 mm of a plane seabed under natural asymmetric waves. The detailed measurements of suspended sediment concentration (at 5 mm vertical resolution and at 4 Hz) showed two or three entrainment bursts around peak flow under the wave crest and another at flow reversal during the decelerating phase. At flow reversal, the mixing length was found to be approximately double the value attained at peak flow under the crest. To examine the cause of multiple suspension peaks and increased diffusion at flow reversal, a numerical “side-view” hydrodynamic model was developed to reproduce near-bed wave-induced orbital currents. Predicted currents at the bed and above the wave boundary layer were oppositely directed around flow reversal and this effect became more pronounced with increasing wave asymmetry. When the predicted orbital currents and an enhanced eddy diffusivity during periods of oppositely directed flows were applied in a Lagrangian numerical sediment transport model, unprecedented and extremely close predictions of the measured instantaneous concentrations were obtained. The numerical models were simplified to incorporate only the essential parameters and, by simulating at short time scales, empirical time-averaged parameterisations were not required. Key factors in the sediment model were fall velocities of the full grain size distribution, diffusion, separation of entrainment from settlement, and non-constant, but vertically uniform, eddy diffusivity. Over the plane bed, sediment convection by wave orbital vertical currents was found to have no significant influence on the results.     

适于模拟不规则水域波浪的缓坡方程两种数值模型比较

本文分析比较了适于不规则水域波浪模拟的椭圆型缓坡方程两种数值模型。两种数值模型均采用有限体积法离散,分别基于四叉树网格和非结构化三角形网格建立。首先结合近岸缓坡地形上波浪传播的经典物理模型实验对两种数值模型分别进行了验证,并结合计算结果对比分析了两种模型的计算精度和效率。计算结果表明,两种数值模型均可有效地模拟近岸波浪的传播变形;相对非结构化三角形网格下的模型,基于四叉树网格建立的数值模型在数值离散和求解过程中无需引入形函数、不产生复杂的交叉项,离散简单,易于程序实现,且节约计算存储空间,计算效率高。     

Modelling 3D hydrodynamics governing island-associated sandbanks in a proposed tidal stream energy site

A 3D numerical modelling study to investigate the existing hydrodynamic regime of the Pentland Firth Inner Sound Channel, Scotland, UK is presented. Hydrodynamics that govern some sensitive sedimentary deposits in the Inner Sound Channel are discussed. A 3D hydrodynamic model Delft3D is set up for Pentland Firth, located between Orkney Islands and mainland Scotland and a full sensitivity analysis of the numerical model is carried out. The current model set up sufficiently captures the existing hydrodynamics during a full spring-neap tidal cycle inside Pentland Firth. Using model results, the 3D structure of the dynamics of the tidal flows in the Inner Sound Channel is investigated. The temporal variability of tidal flows, the residual tidal flows in the channel and local flow interactions with the Island of Stroma are described. It is proved that the tidally dominant flows drive the sediment transport gradient model to explain the principle maintenance mechanisms of two island-associated sandbanks present in the Inner Sound. The present study provides detailed information on the physics of the tidal regime in the Inner Sound and explains the presence of sandbanks in an area of high tidal flows. Due to extremely high tidal flows, Inner Sound is considered as one of the most favourable sites for tidal energy extraction in the UK. The findings of this study will be very useful in assessing the significance of impacts of future tidal energy extraction on natural hydrodynamics and sediment dynamics of the area.     

Sand transport under combined current and wave conditions: A semi-unsteady,practical model

For the general purposes of morphodynamic computations in coastal zones, simple formula-based models are usually employed to evaluate sediment transport. Sediment transport rates are computed as a function of the bottom shear stress or the near bed flow velocity and it is generally assumed that the sediment particles react immediately to changes in flow conditions. It has been recognized, through recent laboratory experiments in both rippled and plane bed sheet flow conditions that sediment reacts to the flow in a complex manner, involving non-steady processes resulting from memory and settling/entrainment delay effects. These processes may be important in the cross-shore direction, where sediment transport is mainly caused by the oscillatory motions induced by surface short gravity waves.The aim of the present work is to develop a semi-unsteady, practical model, to predict the total (bed load and suspended load) sediment transport rates in wave or combined wave-current flow conditions that are characteristic of the coastal zone. The unsteady effects are reproduced indirectly by taking into account the delayed settling of sediment particles. The net sediment transport rates are computed from the total bottom shear stress and the model takes into account the velocity and acceleration asymmetries of the waves as they propagate towards the shore.A comparison has been carried out between the computed net sediment transport rates with a large data set of experimental results for different flow conditions (wave-current flows, purely oscillatory flow, skewed waves and steady currents) in different regimes (plane bed and rippled bed) with fine, medium and coarse uniform sand. The numerical results obtained are reasonably accurate within a factor of 2. Based on this analysis, the limits and validity of the present formulation are discussed.     

Numerical simulation of oscillatory flows over a rippled bed by immersed boundary method

In this paper, a well-developed numerical model based on the immersed boundary (IB) method is used to study oscillatory flows over a bed with large-amplitude ripples in a systematic manner. The work shows that the complex flow over the rippled bed can be numerically dealt with in Cartesian coordinate by the IB method and that the IB method is able to provide main features of the flows near the ripples. An accurate simulation of vortices generation as a result of flow separation at the rippled bed is obtained. It is found that the oscillatory flows start to separate during the flow deceleration when the Keulegan–Carpenter (KC) number is small. The steady streaming for various ripple steepness is simulated and the criterion for separating the single and double structure streaming is also discussed. Moreover, a new type of steady streaming which consists of a pair of embedded recirculations in the vicinity of the ripple trough is obtained for relatively steep ripples in this work. The numerical results, including the steady streaming in particular, may be helpful to improve the understanding of the sediment transport and the seabed evolution with natural ripples under sea waves.     

Numerical study on wave-induced filtration flow across the beach face and its effects on swash zone sediment transport

A numerical model, coupling an analysis of beach groundwater flow with an analysis of swash wave motion over a uniform slope, is presented. Model calculations are performed to investigate the variations of swash-induced filtration flows across the beach face for different input parameters. Swash zone sediment transport under the influence of such filtration flow across the beach face is investigated through modification of effective weight of sediment particle and modification of swash boundary layer thickness. These effects are quantified based on a bed load transport model with a modified Shields parameter.     

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.     

人工岬角和海滩养护下泥沙输运与海床演变

实施人工岬角和海滩养护的软硬措施相结合的方法是探索解决海滩侵蚀问题的新举措之一。基于验证良好的水动力模型建立了泥沙输运和海床演变模型,对秦皇岛市新开河口至南山岸线海滩经整治修复工程实施后的泥沙输运和海床演变规律进行模拟分析。结果表明:大潮与常浪耦合作用下,研究区域周围海域流速整体上落急时刻大于涨急时刻,沙坝掩护水域的流速较弱,波浪在沙坝处增强并在沙坝后破碎并减弱;在沙坝处出现高悬沙浓度区域,悬沙浓度分布主要由底床切应力控制,部分受水动力影响。大潮与强浪耦合作用下,在沙坝处出现远大于常浪时的增强水流(波生流),方向与强浪向一致,波浪在沙坝处已破碎且沙坝对强浪的削减效率大于常浪;与常浪作用下不同,悬沙浓度场在强浪作用下沙坝处没有出现明显的高浓度区域。波流耦合作用下,人工岬角单独存在时护岸效果有限,工程实施后在人工岬角和人工沙坝的配合下,沙坝处侵蚀而沙坝后侧掩护区域淤积,沙坝起到了有效养护海滩的作用。     

四叉树网格下的椭圆型缓坡方程数值模型研究

波浪是近岸海域关键的水动力因素之一。考虑到近岸地形复杂、波浪演化显著的特点,建立了四叉树网格体系下的椭圆型缓坡方程数值模型,采用有限体积法对模型进行数值离散,应用GPBiCG(m, n)算法求解离散后的控制方程。模型中根据波浪波长布局计算网格,生成多层次四叉树网格,对复杂计算域有较好的适应性,并且在离散和方程求解中无需引入形函数、不产生复杂的交叉项,节约了存储空间和计算时间。将模型成功应用于物理模型实验及Acapulco海湾的波浪场数值模拟,结果表明该模型能够准确、高效地模拟近岸波浪场,可为近岸波浪场的模拟提供一定的理论和技术支持。     

Numerical modelling of suspended sediment fluxes in estuarine waters

A two-dimensional finite difference numerical model, capable of predicting depth-averaged tidal flow fields in coastal and estuarine waters, has been extended to include tide-induced non-cohesive sediment transport processes. The partial differential equations governing the conservation of mass, momentum and suspended sediment in an incompressible turbulent flow are included in a depth-integrated form in the model. For the representation of the processes of erosion and deposition of sediment from the bed an empirically based source-sink term was refined, based on the results of three mobile bed flume studies. The model has been tested by simulating tidal flows and suspended sediment fluxes in two estuaries, with particular application to the Humber estuary in the U.K. The model was calibrated and found to produce an encouraging degree of agreement between the numerical predictions and corresponding field measurements for this estuary. Furthermore, the predicted gross deposition and erosion features of the estuary were found to be in close agreement with interpretations from Eulerian tidal residual predictions.     

Numerical modelling of sand beach evolution around coastal structures

I~crIOWIn the coastal area, especially at the sandy seashore, wave and nearshore current are the major factors which affect sediment transPOrt and the motyhChdynamics.The numerical models of predicting the beach evolution can be classified intO the medi~term and long-term models according to their space and time scales (De Briend et al., 1993;Watanabe, 1990; Watanabe et al., 1986; Tao, 1996). In the medium-term model the effects ofwave, nearshore current and sediment transport are conside…     

近岸波、流作用下结构物附近海岸演变的数值模拟

针对与砂质海岸在波浪作用下的演变有关的波浪、近岸流及输沙问题进行了系统的研究,并对结构物附近海岸演变进行了数值模拟。考虑了波浪折射-绕射及波浪破碎的综合作用,在近岸流场的模拟中用沿水深积分形成的K方程模型确定涡粘系数。计算岸滩地形变化时,综合波浪、近岸流作用的底沙和悬沙输沙率,并考虑波浪对泥沙作用的影响。模型对防波堤和近岸沉船附近地形变化进行了模拟,效果良好。     

The Sediment Source and Transport Trends around the Abandoned Yellow River Delta,China

Sediment source and transport trends are influenced by various hydrodynamic factors, and thus play important roles in sedimentary evolution and coastal stability. To examine sediment transport trends around the abandoned Yellow River delta promontory and its erosion mechanism, we employ empirical orthogonal function (EOF) analysis to study sedimentary characteristics and transport trends of the abandoned Yellow River delta in northern Jiangsu Province, China. The results showed that: (1) the main sediment source in the abandoned Yellow River delta is the submarine coastal slope and both sides of the abandoned Yellow River Delta; (2) the main hydrodynamics controlling sediment transport is the current that runs along the shore, coupled with waves, especially southward currents; (3) the sediment of the study area was redistributed under hydrodynamics; coarse sediments were eroded and broadly transported to the south. Therefore, it is concluded that the sediment sources and transport have important influence on coastal evolution: the sediment source area shows mass loss of deposits and erosion; deposits in the submarine coastal slope provide the source and were continuously eroded to provide materials to other places as a sediment source.     

Morphodynamics of ebb-tidal deltas: a model approach

The results of 2DH numerical models of the Frisian Inlet (located in the Dutch Wadden Sea) are discussed to gain further knowledge about the physical mechanisms causing the presence of both ebb-tidal deltas and of channels and shoals in tide-dominated inlet systems. A hydrodynamic model, extended with sediment transport formulations, was used to verify earlier conceptual models that deal with ebb-tidal delta characteristics. The model does not confirm their hypothesis concerning the observed spatial asymmetry of ebb-tidal deltas and suggests that long-term morphological simulations are needed to understand this aspect. Furthermore, the model indicates that the initial formation of the ebb-tidal delta is mainly due to convergence of the tidally averaged sediment flux related to residual currents, whilst the net sediment transport in the basin is mainly caused by tidal asymmetry. A second model (accounting for feedbacks between tidal motion and the erodible bottom) was used to simulate the long-term bathymetric evolution of the Frisian Inlet under fair weather conditions. This model reproduces the gross characteristics of the observed morphology: the presence of a double-inlet system with two distinct ebb-tidal deltas having different sizes and the presence of channels and shoals. The role of the ‘Engelsmanplaat’, a consolidated shoal in the middle of the Frisian Inlet, was not found to be crucial for the morphodynamic stability of this inlet system.