潮汐流动中污染物排放的一种远区计算模型

本文建立了一个适合于潮汐流动中污染物排放的远区动态平面二维数值计算模型。该模型的计算特点是将流场和浓度场分开计算，即所谓“非耦合解”，运用有限节点法计算了流场，用改进的有限节点法计算了浓度场。该模型应用于长江大厂镇江段，其计算出的流场和该江段典型污染物Ｆ－的浓度场与实测一致，在此基础上，用设计潮计算了一个污水排江方案，确定了其影响范围     

徐六泾控制节点污染物运移轨迹模拟

对三维多功能动力-生态耦合模式(COHERENS)进行二次开发,运用"网格冻结法"实现了漫滩和露滩过程中的干湿交替,突破其为固定边界和限制水深的局限,使COHERENS模式成功应用于长江口浅滩过程的模拟。对潮位和水流流速的模拟结果进行了较好地验证。在斜压流场的基础上对示踪颗粒拉格朗日运移进行追踪,对污染物欧拉输运进行数值模拟。以徐六泾控制节点横断面设置颗粒示踪子和污染物排放点,对排放的悬浮颗粒物质和溶解性污染物的运动轨迹和特征进行了模拟和比较分析,给出了模拟期间的流场特征和污染物输移规律。     

Internal inlet for wave generation and absorption treatment

A new method of implementing, in two-dimensional (2-D) Navier–Stokes equations, a numerical internal wave generation in the finite volume formulation is developed. To our knowledge, the originality of this model is on the specification of an internal inlet velocity defined as a source line for the generation of linear and non-linear waves. The use of a single cell to represent the source line and its transformation to an internal boundary condition proved to be an interesting alternative to the common procedure of adding a mass source term to the continuity equation within a multi-cell rectangular region. Given the reduction of the source domain to a one-dimensional region, this simple new type of source introduced less perturbation than the 2-D source type. This model was successfully implemented in the PHOENICS code (Parabolic Hyperbolic Or Elliptic Numerical Integration Code Series). In addition, the volume of fluid (VOF) fraction was used to describe the free surface displacements. A friction force term was added to the momentum transport equation in the vertical direction, in order to enhance wave damping, within relatively limited number of cells representing the sponge layers at the open boundaries. For monochromatic wave, propagating on constant water depth, numerical and analytical results showed good agreements for free surface profiles and vertical distribution of velocity components. For solitary wave simulation, the wave shape and velocity were preserved; while, small discrepancy in the tailing edge of the free surface profiles was observed. The suitability of this new numerical wave generation model for a two source lines extension was investigated and proven to be innovative. The comparisons between numerical, analytical and experimental results showed that the height of the merging waves was correctly reproduced and that the reflected waves do not interact with the source lines.     

长江口南支潮流的一个数值模拟

为了能够精确地拟合天然水域不规则的水陆边界,并能保持有限差分数值计算在矩形网格中进行,本文利用Thompson的数值网格生成技术,在长江口南支七丫口至横沙河段设计了一个椭圆型边界拟合坐标系,并在这个坐标系下建立了平面二维潮流数学模型。通过数值求解此数学模型,实现了计算区域内现状潮流的一个数值模拟,模拟结果表明,模式的设计是成功的。     

斜压海洋动力学的一种三维数值模式 Ⅰ.动力学方程数值格式

依据自由海面海洋动力学原始方程建立了一种三维有限差分数值模式,可用于潮波、风暴潮和海流的数值模拟和预报。运动方程和连续方程的数值格式采用内、外模态分离的技术。外模态采用交替方向隐格式,用于计算海面高度和垂直平均流速,时间步长不受Ｃｏｕｒａｎｔ－Ｆｒｉｄｅｒｉｃｈｓ－Ｌｅｗｙ条件限制;内模态采用半隐格式,用于计算海流的垂直２颁布,其时间步长可大于外模态时间步长。模式的计算程度比一般显式模式可快１０倍     

2001年春季东中国海三维环流的数值模拟

本文基于2001年4月东中国海区域实测水文资料,应用三维有限元模式FEOM(Finite Element Ocean Model)对东中国海三维环流系统进行了数值计算分析。模式水平网格系统采用单节点线性三角形网格,垂向使用z坐标,观测温、盐度场通过客观分析法插值得到初始条件,分别进行了诊断计算和强诊断计算,计算结果表明:(1)改进逆方法可以很好地反演研究区域流函数和流量分布,为数值模拟提供优质可靠的开边界条件。(2)有限元模式在网格自由度方面和对研究区域的完整覆盖方面优势明显,高分辨率的垂向z坐标也可以较好地拟合海底地形,从而得到较高分辨率的三维数值模拟结果。(3)通过诊断计算,模拟再现了东中国海春季环流的多涡结构,分辨出了台湾暖流、黄海沿岸流、黄海暖流等流系。(4)比较诊断与强诊断两个计算结果,它们在定性上较为一致,在定量上有些差别,总体分布强诊断计算结果更为合理。     

An analytical solution for the local suspended sediment concentration profile in tidal sea region

The time-averaged and oscillatory solutions of the one-dimensional vertical (1DV) advection–diffusion equation for the suspended sediment have been derived analytically in a tidal sea region of finite water depth. The basic equation assumes constant eddy diffusivity and settling velocity. No net flux condition is set at the sea surface, while a boundary condition with the erosion rate and depositional velocity is prescribed at the sea bottom. The time-averaged solution has been derived in a straightforward manner, while the advection–diffusion equation governing the oscillatory concentration has been first transformed to a simple diffusion equation and then solved using the Galerkin-eigenfunction method. The former is given in a closed form, while the latter is presented in a series solution.A set of calculations has been performed to examine the change in the vertical structure as well as magnitude of the concentration response function. A possible use of the solution to make an estimate of the erosion rate at the sea bottom based on the concentration information at the sea surface is discussed.     

长江口北港口门海域悬沙输运机制分析

根据2012年10月24日至11月1日在长江口北港口门海域获取的沉积动力学数据,采用输运通量分解方法分析水沙输运机制。结果表明,长江口北港口门附近海域涨、落潮期间底部悬沙浓度与近底部流速呈显著线性相关,存在显著的再悬浮作用;潮周期内的悬沙输运具不对称特征,涨潮悬沙浓度大于落潮悬沙浓度,悬浮泥沙有向陆输运的趋势。拉格朗日平流输运是影响悬沙输运的主要贡献项;潮泵效应(尤其是潮汐捕捉效应)以及河口垂向环流作用是两个次要影响因素,在影响程度上前者比后者略大。观测发现,长江口北港口门海域潮流除了具有涨落潮流速、历时等不对称现象外,还具有流速结构不对称的特征,进而导致涨、落潮底部湍流混合程度不对称与输沙不对称,这可能是造成悬沙向河口内输运形成最大浑浊带的重要因素。     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

POM模式在河口湾污染物质输运过程模拟中的应用

具有复杂地形的河口湾,由于冲淡水的影响,湾内的流动具有很强的三维结构.因此要再现河口湾各种过程,必须进行三维数值模拟.在实际应用中,选用水平上采用正交曲线性网格和垂向上采用σ坐标的POM模式作为水动力模式比较合适.本研究在POM模式的基础上添加了示踪粒子三维Lagrange运动轨迹计算及水质模块,用以研究榆林湾水交换能力和污染物质输运过程.     

A Study of Eddy Viscosity Coefficient in Numerical Tidal Simulation

Based on the fluid motion equations, the physical meaning of eddy viscosity coefficient and the rationality of the Boussinesq hypothesis are discussed in this paper. The effect of the coefficient on numerical stability is analyzed briefly. A semi-enclosed rectangular sea area, with an orthogonal spur dike, is applied in a 2-D numerical model to study the effect of horizontal eddy viscosity coefficient (AH). The computed result shows that AH has little influence on the tidal level and averaged flow velocity, but has obvious influence on the intensity and the range of return flow around near the spur dike. Correspondingly, a wind-driven current pool and an annular current are applied in a 3-D numerical model respectively to study the effect of vertical eddy viscosity coefficient (Av). The computed result shows that the absolute value of Av is inversely proportional to that of horizontal velocity, and the vertical gradient value of Av determines the vertical distribution of horizontal velocity. The distrib     

杭州湾的三维水流数值模拟

建立了一个基于σ变换和内外模式分裂技术三维水流数学模型，采用有限节点法（平面）和变步长差分法（σ向）对方程进行了离散；通过风生流和环岛水流两个理论模式对所建模型进行了校验，并将该模型应用到杭州湾的三维水流数值模拟中，效果良好，所建模型垂向分辨率高，简单实用，可应用于河口，海岸，湖泊等大范围水域的水流数值计算。     

象山港水交换数值研究──Ⅰ.对流-扩散型的水交换模式

以溶解态的保守性物质作为湾内水的示踪剂,建立了对流－扩散型的海湾水交换数值模型。数值模型使用参数化的方法把重力环流和潮振荡的垂向剪切作用的水平混合效应包纳在水平二维的示踪剂对流－扩散方程中。在空间网距较小时,模型的稳定性和守恒性均可满足海湾水交换研究的需要。     

胶州湾及其邻近海域潮流和污染物扩散的数值模拟

采用曲线网格技术，将ECOMSED模式应用于胶州湾的潮流和污染物COD的稀释扩散研究，建主了一个三维保守污染物COD输运的对流扩散数值模型，对胶州湾M2分潮和COD污染状况进行了模拟。结果表明，胶州湾及其附近海域存在许多大小、强弱和旋转方向不同的余流涡；湾内COD浓度东部比西部高，湾口外的COD浓度等值线成舌状分布，向东伸展；余流和COD浓度垂向分布差异并不显。模拟结果与胶州湾的实测资料符合较好。     

单点系泊潜标系统动力学建模与姿态分析

通过对单点系泊潜标进行动力学分析,利用集中质量法建立系统的动力学模型,编制Matlab仿真程序,通过数值模拟获取不同海流条件下,分析潜标系统的姿态。文中将一个典型潜标系统的浮球、锚链、重力锚等部件,利用集中质量法简化为节点,通过模拟真实环境中的海流流速在垂直方向的分布,设计1～5 m/s的流速剖面分布。将节点的物理参数和海流参数代入仿真程序中,得到各节点的姿态、节点间张力以及重力锚与海底摩擦力等数据。分析仿真结果,海流流速与节点的偏降、节点间的张力以及重力锚与海底的摩擦力成正比,仿真结果与实际情况相符。文中设计的建模方法和仿真程序可以为单点系泊潜标的总体设计、姿态估算、锚系配置提供工程参考。     

Modeling of coastal effluent transport: an approach to linking far-field and near-field models

One of the challenges in effluent transport modeling in coastal tidal environments is the proper specification of initial dilution in connection with the far-field transport phenomena. An approach of external linkage of far-field and near-field effluent transport models is presented, and applied to simulating the effluent transport in the Port Angeles Harbor, Washington in the Strait of Juan de Fuea. A near-field plume model was used to calculate the effluent initial dilution and a three-dimensional (3-D) hydrodynamic model was developed to simulate the tidal circulation and far-field effluent transport in the Port Angeles Harbor. The hydrodynamic model was driven by tides and surface winds. Observed water surface elevation and velocity data were used to calibrate the model over a period covering the neap-spring tidal cycle. The model was also validated with observed surface drogue trajectory data. The model successfully reproduced the tidal dynamics in the study area and good agreements between model results and observed data were obtained. It is demonstrated that the linkage between the near-field and far-field models in effluent transport modeling can be achieved through iteratively adjusting the model grid sizes such that the dilution ratio and effluent concentration in the circulation model grid cell match the concentration calculated by the near-field plume model.     

迎风有限元法在三维潮流数值模拟中的应用

为研究和开发迎风有限元法在河口海岸三维波数值模拟中的应用，并为建立三维物质扩散模型，最终为解决实际工程问题打下基础，以非线性浅水波方程为基础，采用一咱平面迎风有奶元与垂向工有限差分相结合地数值计算方法建立三维潮流数学模型，基地实际的物理过程，在计算中将三维流动分成外重力波和内重力波，耦联求解潮位和流速的空间分布。应用本模式计算了天津新港附近渤海海域的三维潮流运动，以较高的分辩率揭示了潮流空间结构特     

一种污水排海的近远区计算与控制方法──以大亚湾水容量计算及污水排海为例

针对点容量计算中混合区限制的要求，采用了动量积分法、经验公式法和有限差分法，分别对污水排放后可能引起的近区、过渡区和远区的污水浓度变化进行了模拟预测，并就一定的限制性条件对容量点近区、高浓度混合区及容量点远区影响范围进行了分析。在大亚湾水容量计算及污水排海规划中的应用表明，该法具有计算机时较省、精度相对较高等优点。     

洋浦经济开发区污水排海总量控制

研究建立洋浦近岸海域岸界自适应网格变动边界的潮流模型及物质输运模型。对开发区提出的多种污水排海方案进行污水质点拉格朗日追踪、污水排海的水质预测及工程投资分析 ,推荐了最佳排污口位置 ,最后计算了污水最大排放量及 8种污染物控制指标排放总量 ,为洋浦开发区的环境管理与保护提供决策依据。