A Three-dimensional Tidal Model in Boundary-fitted Curvilinear Grids

In hydrodynamic models of marine and estuarine currents the use of boundary-fitted curvilinear grids not only makes the model grids fit to the coastline and bathymetry well, but also makes the kinetic boundary conditions simple and more accurate. Because of these advantages, a three-dimensional tidal model with boundary-fitted curvilinear grids has been developed to simulate both tide and current in estuarine and shelf water. The basic idea is to use a set of coupled σ-stretched and elliptic transformations to map the physical space into a corresponding transformed space such that all boundaries are coincident with co-ordinate lines and the transformed grids are rectangular. The hydrodynamic equations in the transformed space are solved in a rectangular mesh of the transformed grid system. The application of the model in the simulation of the M₂,S₂ , K₁and O₁tidal waves in the Bohai Sea, China, shows that the numerical results are in good agreement with the observations. The method of boundary-fitted curvilinear grids is effective for improving simulation accuracy of current in the estuarine and shallow seas, especially at the coastal regions where the current was usually impractical by a uniform Cartesian grid system with the shoreline and bathymetry represented by numerous stair-steps.     

An Application of Finite Volume WENO Scheme to Numerical Modeling of Tidal Current

1 .IntroductionThe development and application of hydrodynamic models to predict circulation in estuarine ,coastal ,and shelf waters has been extremely active over the past three decades (Spauldinget al .,1992 ,1994) . However , most of these models are b…     

潮流与泥沙数值模拟回顾及进展

回顾和总结了近年来国内外潮流与泥沙数值模拟的进展,评价分析了数学模型的多种数值解法,特别是有关二维、三维潮流数值模拟以及水沙数学模型近年来的进展,并根据目前潮流、泥沙数值模拟技术中存在的不足,讨论了进一步研究的方向和发展趋势。     

A Finite Volume Method with Unstructured Triangular Grids for Numerical Modeling of Tidal Current

The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.     

葫芦岛港柳条沟港区规划潮流泥沙数模研究

建立了基于不规则三角形网格的考虑波浪作用的潮流泥沙数学模型,在对模型进行充分验证的基础上,对葫芦岛海域的潮流、泥沙场的进行了数值模拟,对该海区的潮流扬、泥沙场特征进行了分析;对萌芦岛港柳条沟港区规划方案进行了潮流、泥沙场模拟研究。详细分析了工程前后港池和航道的流扬变化,计算了航道的横流,预报了规划港区港池和航道的泥沙年淤积量,为港区的规划设计提供了科学依据。     

Numerical Simulation of Flow Field in Tidal Estuary Using Body-Fitted Coordinate System

- On the basis of the physical mechanism, a body-fitted coordinate system is developed. By using this system the boundaries in simulation and in real are fitted well, and simulation with great accuracy is achieved. A computation example indicates that compared to traditional two-dimensional computation methods, the body-fitted simulation has an advantange of better coincidence with the real and can be adopted in simulating flow fields in tidal estuaries.     

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

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

Corrected Formula of Bed Resistance Coefficient for Plane Numerical Simulation of Tidal Current

By means of a logarithm law for the velocity profile,a corrected formula of bed resistancecoefficient,which involves many factors such as gradient of still water depth,variation of surfaceelevation,flow direction,and so on,is derived from the 3D governing equations of tidal current by aver-aging over the whole water depth.Theoretical analysis and application have shown that the 2D plane tidalcurrent numerical model would be more reasonable and could be applied to steep bottom topographywhen the corrected bed resistance coefficient is used,therefore the results of reproduction simulation andengineering calculation would be more scientific and reasonable.     

Numerical Study on Tidal Mixing in the Bohai Sea

The spatial and temporal variability of tidal mixing in Bohai Sea is studied using a numerical approach. In calculating tidal mixing, accurate barotropic tidal current is obtained via a harmonic analysis package utilizing the simulated current output from a high-resolution regional ocean model. And a “small-scale” roughness map is adopted to describe the detailed topographic features of Bohai Sea. It is shown that the tidal mixing estimated in Bohai Sea is much higher than the level of global background, and fluctuates considerably at some regions within a single day. In Liaodong Bay, Bohai Bay and Bohai Strait, the mixing varies greatly, with the peak value of O (10^?2) m² s^?1. The order of magnitude of mixing in Laizhou Bay is about O (10^?5～10^?3) m² s^?1. Mixing with background level of O (10^?5) m² s^?1 only appears in central area. Result also shows that rough topography plays relatively a more important role than tidal current in enhancing diapycnal mixing in Bohai Sea. The distributions of tidal mixing in selected sections reveal that the vertical stratification in Bohai Sea is not obvious, generally renders a barotropic structure.     

胜利油田孤东海堤保滩促淤工程二维潮流数值模拟研究

以胜利油田孤东海堤保滩促淤工程为实例，采用二维潮流数值模拟方法，利用实测流速资料对该海区进行潮流计算，对所提出的不同工程方案进行比较，并分析其促淤效果，在此基础上给出工程方案的推荐结论。     

Numerical Modeling of Wave and Current Interaction in the Tidal Inlet Area

This paper investigates the phenomena of wave refraction and diffraction in the slowly varying topography, as well as the current deflection due to wave actions. A numerical model is developed based on depth integrated mean continuity equation and momentum equations, and a 3rd-order wave equation which governs the wave diffraction, refraction and interaction with current. Examples to examine the above model are given comparing with the laboratory data or the numerical results of other researchers. An example simulating the inlet area shows the interesting velocity field which may be used as a pioneer to further study on the nearshore hydrodynamics and sedimentation.     

广西近岸海域潮流数值模拟

采用有限元三角形网格的分步杂交方法,建立了广西近岸海域的二维潮流数值模型,计算值与实测资料符合较好。采用主要分潮组合输入,模拟了研究海域的平均潮潮流场。模拟结果表明:涨急时,潮流向为东北方向,最大涨潮流速为74cm/s左右;落急时,潮流向为西南方向,最大落潮流速约100cm/s,落潮流速大于涨潮流速。近岸区域潮流为往复流,离岸边越远潮流越接近旋转流。     

非结构网格生成法在二维潮流数学模型中的应用

文章基于Delaunay三角化的思想实现了非结构化数值网格生成技术,并LOP(Local Optimization Procedure)算法对网格进行了优化。用有限体积法在非结构三角形网格上对平面二维潮流方程进行了数值求解。实例应用分析表明,文中的数值网格生成技术和有限体积方法相结合用于求解具有复杂几何边界区域内的潮流问题能取得很好的结果。     

潮流能量的提取及对其流速的影响

本文在经过验证的Quoddy Region海域三维潮流数值模型的基础上,植入了提取潮流能的数值计算模块。确定了Western Passage在3种涡轮机布设方式下能够提取的能量以及提取能量的效率,其中采取3排涡轮机交叉布设的方式时,1个大小潮时间内提取的能量最高,可达1 148.1MWh,提取效率为4.57%,远高于另外2种方式。在单排涡轮机布设的实验中(实验一),通过改变速度衰减系数,得到了可提取能量随速度衰减系数的变化曲线,该曲线与前人在理论分析时得到的曲线类似。同时本文还分析了提取潮流能给周围流场带来的影响。     

胶州湾潮汐潮流高分辨率数值模拟研究

基于采用无结构网格和有限体积方法的FVCOM陆架模式,考虑8个主要的天文分潮,建立胶州湾三维高分辨率数值模型来重现和研究其潮汐潮流变化状况。与实测资料对比验证表明,模拟结果与实测值吻合较好。在此基础上,根据模拟结果计算得到了较以往更为精细的同潮图和潮汐、潮流、余流分布特征。研究结果揭示,最大可能潮流和最大余流都发生在团岛附近,流速分别可达2.14和0.43 m/s;除了湾口附近前人报道过的4个余流系统外,还在中部首次揭示了2个相对较弱的余流系统;潮流能通量在内外湾口呈"左进右出"的结构;胶州湾的平均纳潮量为8.31亿m3;染色试验表明,胶州湾30 d的水交换率为36.8%。     

印度尼西亚近海潮汐潮流的数值模拟

利用FVCOM海洋数值模式计算了印尼近海的M2,S2,K1,O1分潮的分布,计算范围从20°S～20°N,90°～150°E,计算网格分辨率在印尼海域岛屿平均为1/12度,在大陆边界平均为1/5度,在开边界平均为1/2度.计算结果与104个TOPEX/Poseidon卫星高度计交叉点数据和79个验潮站数据进行比较,符合良好;与高度计交叉点比较,M2分潮振幅的均方根差为6 cm,迟角为7°;S2分潮的振幅偏差为3 cm,迟角偏差为8°;K1分潮振幅的偏差为6 cm,迟角偏差为10°;O1分潮振幅偏差为3 cm,迟角偏差为10°.根据计算结果给出了4个分潮的潮汐、潮流、潮余流和潮能通量密度分布图.     

较大尺度下潮流物模试验理论可行性的数值研究

以渤海及其局部海域为例,通过数值模拟量化地转缺失对水动力过程的影响,评价不同空间尺度下实施潮流物理模型试验的理论可行性。结果表明,对整个渤海、辽东湾这么大尺度的海域实施物模试验理论上是不可行的;对90km尺度的渤海湾,地转缺失使模拟结果存在较大偏差;对莱州湾61km尺度的海域可进行物模试验。研究确定60km为实施潮流物模试验一般的理论允许尺度。     

潮流能水平轴水轮机湍流模型研究初探

借助计算流体力学(CFD)方法研究潮流能水平轴水轮机水动力学性能,选择合适的湍流模型至关重要。为此,本文采用计算流体力学软件Fluent,将Spalart-Allmaras湍流模型、标准k-ε湍流模型、RNGk-ε湍流模型、标准k-ω湍流模型和SSTk-ω湍流模型分别应用到水平轴水轮机数值模拟当中。对数值模拟结果进行分析比较,并与模型试验数据进行对比,得出应用SSTk-ω湍流模型可以更好的模拟潮流能水平轴水轮机的水动力学性能。     

水平轴潮流能水轮机数值模拟方法对比分析

叶素动量理论和CFD方法是水平轴潮流能水轮机性能分析中运用较为广泛的数值模拟方法,文中结合小尺寸水轮机模型试验,对比分析了叶素动量理论和CFD方法在水轮机性能分析中的准确性和适用性.验证结果表明:叶素动量理论和CFD方法均能对水轮机的性能进行预报,且CFD精度高于动量叶素理论;大尖速比时,动量叶素理论偏离较高,不再适合性能预报;在小尖速比下,建议采用RNGk-ε模型的CFD方法进行分析计算;动量叶素理论适合设计初期设计方案的对比分析,而CFD方法适合对设计结果的验证校核和详细分析.     

渤海湾三维变动边界潮流数值模拟

基于河口、陆架和海洋模式 (ECOM模型 ) ,引入干湿网格法模拟潮滩涨落的改进 ,建立了渤海湾三维变动边界潮流模型。该模型网格采用为 4 80 m× 6 2 0 m,垂向 4层 ,充分考虑渤海湾潮滩广阔的特点。计算结果与实测值比较符合良好 ,较好地刻画出渤海湾 M2 分潮潮流场的时空分布特点。