A higher-efficient three-dimensional numerical model for small amplitude free surface flows

A higher-efficient three-dimensional non-hydrostatic model is developed to simulate small amplitude free surface flows based on a staggered unstructured grid. In this model, a fractional step algorithm is adopted to solve the Navier-Stokes equations in two major steps. A top-layer pressure method is proposed to minimize the number of vertical layers and subsequently the computational cost. Three classical examples of small amplitude free surface flows are used to demonstrate the capability and efficiency of the model. The satisfactory results demonstrated the capability and efficiency of modelling a range of small amplitude free surface flows with only a small number of vertical layers.     

A Hybrid Finite-Volume and Finite Difference Scheme for Depth-Integrated Non-Hydrostatic Model

A depth-integrated, non-hydrostatic model with hybrid finite difference and finite volume numerical algorithm is proposed in this paper. By utilizing a fraction step method, the governing equations are decomposed into hydrostatic and non-hydrostatic parts. The first part is solved by using the finite volume conservative discretization method, whilst the latter is considered by solving discretized Poisson-type equations with the finite difference method. The second-order accuracy, both in time and space, of the finite volume scheme is achieved by using an explicit predictor-correction step and linear construction of variable state in cells. The fluxes across the cell faces are computed in a Godunov-based manner by using MUSTA scheme. Slope and flux limiting technique is used to equip the algorithm with total variation dimensioning property for shock capturing purpose. Wave breaking is treated as a shock by switching off the non-hydrostatic pressure in the steep wave front locally. The model deals with moving wet/dry front in a simple way. Numerical experiments are conducted to verify the proposed model.     

基于ECOMSED模式的兴化湾三维潮流泥沙数值模拟研究

通过三维ECOMSED模式,对兴化湾内的泥沙变化规律进行了分析论证.模拟结果与实测值的对比分析表明,模拟结果较好,可以较准确地描述兴化湾的潮流泥沙特征.分析发现,兴化湾大部分海域的含沙量很低,兴化湾内普遍呈现冲淤稳定的状态,海区的冲淤变化量维持在1个较小的水平,淤积主要集中在近岸浅滩区域,近岸年淤积变化量普遍呈增大趋势.     

东海舟山群岛海域表层沉积物运移特性

通过潮流和表层沉积物粒度参数等实测资料,对舟山群岛海域表层沉积物的运移趋势、净输运通量以及起动流速进行了初步分析.结果表明,舟山群岛对研究区表层沉积物分布影响明显;舟山群岛海域表层沉积物整体由湾口向外海输运,平均日单宽输沙通量由杭州湾内向外海迅速减小;从杭州湾向外海,泥沙起动越来越难,单纯潮流对表层沉积物再悬浮作用相对较弱.     

日照帆船港工程泥沙运移物模试验研究

通过动床波浪输沙模型试验,研究了日照帆船港3种设计方案的泥沙运移情况.通过分析试验结果,预报了常浪和大浪作用下工程区的冲淤分布.结果表明,日照帆船港工程的建设不会对邻近海岸产生大的影响,航道和港池均不会出现大的淤积.但是,为防止堤脚冲刷,应在挡沙坝堤脚处设置护脚结构.     

A Numerical Model for Natural Backfill of Pipeline Trenches Subjected to Unidirectional／Oscillatory Flows

A numerical model for the self-burial of a pipeline trench is developed. Morphological evolutions of a pipeline trench under steady-current or oscillatory-flow conditions are simulated with/without a pipeline inside the trench. The oscillatory flow in this study represents the action of waves. The two-dimensional Reynolds-averaged continuity and Navier-Stokes equations with the standard k-e turbulence closure, as well as the sediment transport equations, are solved with the finite difference method in a eurvilirrear coordinate system. Both bed and suspended loads of sediment transport are included in the morphological model. Because of the lack of experimental data on the backfilling of pipeline trenches, the numericalmodel is firstly verified against three closely-relevant experiments available in literature. A detailed measurement of the channel migration phenomenon under steady currents is employed for the assessment of the integral performance of themodel. The two experimental results from U-tube tests are used to validate the model‘s ability in predicting oscillatory flows. Different time-marching schemes are employed for the morphological computation under unidirectional and oscillatory conditions. It is found that vortex motions within the trench play an important role in the trench development.     

Nonlinear Model of Mean Free Surface of Waves

-A nonlinear model of mean free surface of waves or wave set-up is presented.The model isbased on that of Roelvink(1993),but the numerical techniques used in the solution are based on theWeighted-Average Flux(WAF)method(Watson et al,1992),with Time-Operator-Splitting(TOS)usedfor the treatment of the source terms.This method allows a small number of computational points to beused,and is particularly efficient in modeling wave set-up.The short wave(or primary wave)energy equa-tion is solved by use of a more traditional Lax-Wendroff technique.A nonlinear wave theory(James,1974)is introduced.The model described in this paper is found to be satisfactory in most respects whencompared with the measurements conducted by Stive(1983)except in modeling the mean free surface veryclose to the mean shoreline.     

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.     

贴体正交曲线坐标系下水流泥沙数学模型的构建与应用

为了研究深圳河河口福田保税区围垦填土占用部分河道对深圳河治理工程防洪能力和河床冲淤变化的影响情况，建立了贴体正交曲线坐标系下的平面二维的深圳河河道-河口-海湾整体水域的水流泥沙数学模型，模型中采用双步全隐有限差分法离散求解水流运动方程，悬沙不平衡输运方程采用三阶QUIKEST格式离散后用交错方向法求解。在经实测水文资料验证的基础上，应用数学模型计算了8种洪峰流量和河口高潮位组合4种围垦填土的挖除方案下深圳河沿程水位和河床冲淤变化，讨论了深圳河治理工程防洪能力降低的程度和近河口段洲滩冲刷以及所造成香港一侧红树林的破坏情况，提出了优化挖除方案。     

滩海油田泥沙运移控制初步方案及试验研究

为了降低成本,维持滩海油田的开敞式开发,提出了滩海油田泥沙运移控制初步方案,并进行了物模试验验证。该方案充分利用自然力,通过人工干预促使防护区外泥沙起动,借助海上潮流使高含沙海水流经防护区,实现岸滩蚀退的快速防护,对滩海油田的防护具有指导意义。     

水沙环境中重金属迁移转化模型的两个解析解

选取整体模型,将水体中的溶解态(水相)和吸附态(固相)重金属作为整体,开展了水沙环境中重金属迁移转化数学模型的解析解研究.在一定的简化条件下,得出了泥沙单纯沉积和单纯侵蚀两种特定动床条件下一维重金属迁移转化模型的两个解析解.文中成果可用于数学模型的验证和完善.     

黄河钓口河口行水期泥沙输运过程的三维数值模拟

由于钓口流路行水时期普遍缺乏对河口和近海的系统观测,基于1976年黄河三角洲实测水深数据和汛期水沙数据,采用EFDC三维数值模式对钓口河口泥沙沉积动力过程进行了数值模拟研究。模型结果表明,钓口河口的羽状流和异重流的时空分布具有显著的潮周期变化特征。受与岸线平行的涨落潮流作用,表层羽状流的侧向摆动幅度较大,河口羽状流输沙主要平行于岸线方向,跨等深线的泥沙输运受到限制。汛期黄河入海的高浓度泥沙在钓口河口形成潜没的异重流,在底层沿河口轴线向北输运泥沙至三角洲前缘区域,泥沙输运通量比表层羽状流的输送通量高出一个数量级,且随着离岸距离的增大而呈指数形式快速衰减。随着水深增大,单宽泥沙通量的衰减速率逐渐减小。计算结果显示,由于汛期高浓度泥沙入海,钓口河口的异重流稳定存在,导致约75%的入海泥沙沉积在10 m水深以浅的区域,25%的入海泥沙淤积在三角洲前缘外侧,异重流过程对汛期钓口河口泥沙向深水区输运和沉积具有控制作用。数值模拟结果与前人根据水深变化揭示的三角洲冲淤分布格局一致,与钻孔的沉积记录基本吻合。     

A 2D Mathematical Model for Sediment Transport by Waves and Tidal Currents

In this study, the combined actions of waves and tidal currents in estuarine and coastal areas are considered and a 2D mathematical model for sediment transport by waves and tidal currents has been established in orthogonal curvilinear coordinates. Non-equilibrium transport equations of suspended load and bed load are used in the model. The concept of background concentration is introduced, and the formula of sediment transport capacity of tidal currents for the Oujiang River estuary is obtained. The Dou Guoren formula is employed for the sediment transport capacity of waves. Sediment transport capacity in the form of mud and the intensity of back silting are calculated by use of Luo Zaosen＇ s formula. The calculated tidal stages are in good agreement with the field data, and the calculated velocities and flow directions of 46 vertical lines for 8 cross sections are also in good agreement with the measured data. On such a basis, simulations of back silting after excavation of the waterway with a sand bar under complicated boundary conditions in the navigation channel induced by suspended load, bed load and mud by waves and tidal currents are discussed.     

基于第三代浅水波浪模型的海区泥沙骤淤数值预报及工程应用

在简要介绍现有近岸波浪数值模型的基础上,对应用了最新波浪研究成果、适用于海岸、湖泊和河口地区的第三代浅水波浪数值预报模型的研究和应用进行了详细论述。通过该模型数值计算得到黄程山局部区域不同风向的有效波高分布,将模型与海区泥沙骤淤机理相结合,进而得到整个区域平均含沙量的分布,并数值预报近岸区域骤淤量随时间的变化,研究成果可为海区大风浪情况下骤淤预报提供理论指导,本文最后讨论了泥沙骤淤计算模型未来发展的趋势。     

Study of A Geo-Acoustic Model of Gas-Bearing Sediment and Its Application in Sediment with Low Acoustic Veloctiy

A new geo-acoustic model for gas-bearing sediment is proposed based on the work of Dvorkin and Prasad, and Biot theory. Only five geophysical parameters: sediment mineral composition, free gas saturation, tortuosity (also known as the structure factor), permeability, and porosity, are considered in the model. A benefit of this model is that we need only five parameters instead of ten parameters in the Biot's formulas for acoustic velocity and attenuation calculation. Here the model is demonstrated with the in-situ experimental data collected from the Hangzhou Bay, China. The results of this study suggest that free gas content in sediment is the most critical condition resulting in a low acoustic velocity (compressional wave). The respective contributions of the other four parameters in the model are also discussed.     

Study of A Geo-Acoustic Model of Gas-Bearing Sediment and Its Application in Sediment with Low Acoustic Velocity

A new geo-acoustic model for gas-bearing sediment is proposed based on the work of Dvorkin and Prasad, and Biot theory. Only five geophysical parameters: sediment mineral composition, free gas saturation, tortuosity (also known as the structure factor), permeability, and porosity, are considered in the model. A benefit of this model is that we need only five parameters instead of ten parameters in the Biot’s formulas for acoustic velocity and attenuation calculation. Here the model is demonstrated with the in-situ experimental data collected from the Hangzhou Bay, China. The results of this study suggest that free gas content in sediment is the most critical condition resulting in a low acoustic velocity (compressional wave). The respective contributions of the other four parameters in the model are also discussed.     

区域性海洋环流数值模式研究及对南海环流与海峡流量的模拟

基于MOM模式的物理框架,妥善考虑了开边界的物理过程,改造和发展了一个区域海洋数值模式。本模式不仅可以方便地调整开边界条件,使之满足边界的特定物理条件,而且可以方便地做针对性修改,使模式更加可靠。改进后的模式具有MOM模式物理概念明确、公式便于理解、结果便于表达的全部特点,同时克服了MOM模式边界条件不完整、程序不易调整、参数难以改变的缺点。区域性模式比全球模式的计算速度快很多倍,可以成为区域性研究的有效工具。将此模式应用于南海,利用Hellerman＆Rosenstein气候态风应力驱动模式10a,得到与全球模式效果相当的结果。模式模拟结果展现了南海流场的季节特征,在模式分辨率下表现出了多涡结构。根据模拟的流场计算了南海与其它海域的水交换通量。在年平均意义下,外海水通过吕宋海峡进入南海,南海水通过台湾海峡、民都洛海峡和卡里马塔海峡流出南海。各海峡水通量具有明显的季节变化。     

黄岛电厂取水工程潮流泥沙数值模拟

建立了基于不规则三角形网格的考虑波浪及其破碎作用的二维潮流场和泥沙场数学模型,对黄岛电厂取水海域的潮流场、泥沙场进行了数值模拟和分析,对电厂取水口的泥淤积强度进行了计算。计算结果表明,电厂取水口附近海区流弱水清,取水口泥沙淤积轻微,对电厂取水基本没有影响。     

一种新型的非对称台风海面气压场和风场模型

采用台风外围闭合特征等压线等要素,对圆对称的藤田气压模型进行改进,推导出不含有关参数的非对称台风海面气压场和考虑到径向梯度风的风场模型,并利用该模型对给湛江地区造成了严重风暴潮灾的几次台风过程的海面气压场和风场进行了模拟。结果表明,改进后的台风气压场和风场模型更方便准确。     

三维海洋溢油预测模型的建立(英文)

提出一个海洋溢油三维物理归宿和输运的动力学综合模型。该模型包含了一系列的数学公式来描述对流、湍扩散、表面扩展、铅直扩散、乳化和蒸发过程。每一公式的建立是独立的并且与相关过程、环境和其它参数相联系。该模型需要输入流场作为输运的媒介 ,这可从感兴趣区域的三维潮和风驱动的流体动力学模型获得。模型用来预测和后报溢油在海洋环境中的归宿和输移 ,可为溢油应急反应和环境影响评价服务。