内孤立波波-波相互作用的数值模拟

基于弱二维的KP方程,并结合南中国海东沙群岛附近内孤立波的观测资料,模拟了内孤立波的波-波相互倌用0数值结果较好的反应了内孤立波的二维特征,同时体现两个内孤立波波-波相互作用的非线性特征,即两波相交处相速随振幅的增大而变大。相比于一维的KdV方程,KP在内孤立波的仿真反演方面具有更大的优势。     

网衣对波浪传播影响的三维数值模拟

为研究网衣对波浪传播的影响,采用多孔介质模型模拟网衣,建立了模拟网衣与波浪相互作用的三维数值波浪水槽模型。基于该数值模型研究了波浪经过网衣作用后波高衰减的变化规律,并与物理模型试验结果进行比较。通过对比数值模拟和物理模型试验的结果,证明了该三维数值模型模拟网衣对波浪传播影响的可行性。     

台湾海峡—种三维风海流数值计算模式

本文根据浅海风海流基础理论，采用数值计算方法，对台湾海峡的风海流流场进行了三维数值模拟，讨论了台湾海峡特征季风情况下风海流流场的水平结构及垂直分布。本模式值得进一步研究和完善。     

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

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

内潮的一种分层三维数值模式

在中国近海的不少海域曾经发现相当强的内潮现象,例如赵俊生(1992)报道了渤海海峡的内潮现象, Yamashiro(1988)在东海陆坡处亦发现内潮的存在。特别是Amoco石油公司在南海陆坡处进行了一年多的海流测量,测得最大的全日潮流振幅可超过60cm/s,这里正压潮流小于10cm/s,故观测到的潮流主要与内潮有关。迄今关于内潮的数值模式大多将表面潮作为已知输入条件、从而使模式的应用受到较多限制(Jiang and Fang，1992)。Heaps(1983)曾发展了一种自由海面多层模式,但他们忽略了非线性平流项,这无疑也会损失许多有用信息。本文将叙述一种自由海面的多层模式,它能够同时模拟表面潮和内潮,并包含了非线性平流项。     

海岸河口三维潮流数学模型

本文采用有限元法建立了适用于海岸河口浅水地区的三维潮流数学模型，垂向采用绝对分层坐标系统，将整个水柱分成若干层，在每层内通过垂向积分平均，将三维问题简化为多个平面二维问题。在求解有限元方程中，引入集中质量矩阵技术，在时间上采用两步ＬａｘＷｅｎｄｒｏｆｆ格式，使有限元方程直接以显式解出，不需联立求解，节省了大量的计算时间和计算内存，通过模拟表面切应力作用在矩形水地上而引起的水流，计算结果与分析解比较一致，并将本模型应用到香港维多利亚水道中，计算结果与实测值亦符合较好，证明本模型是一个实用而有效的三维潮流数学模型。     

河口海岸三维双σ坐标斜压水流数值模型研究

在海底地形陡变、垂向密度分层明显的水域,三维σ坐标模式中会出现一种"伪"水平斜压梯度力,并会引起"伪"密度流,以至于影响模拟的精度。垂向上引入双σ坐标变换,建立河口海岸水域三维斜压水流数值模型。数值试验结果表明,在海底地形陡变水域,双σ坐标模式可以减小水平斜压梯度力处理引起的误差。     

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

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

东海对马暖流水来源的数值模拟

本文通过二维数值模拟对1986年6月～1988年12月东海对马暖流水的来源问题进行了初步探讨,结果得出东海对马暖流水的来源基本上分为三种类型:(1)东海对马暖流水主要为东海黑潮水继续北上部分构成;(2)东海对马暖流水由东海黑潮水、东海陆架水以及东海北部黄海大陆沿岸水几部分混合而成;(3)东海对马暖流水几乎全部由东海北部的黄海大陆沿岸水构成。模拟与实测结果基本一致.     

不完全立波的二维数值波浪水槽模拟

利用二维数值波浪水槽研究波浪在防波堤前的运动状态.模型选用Navier-Stokes方程作为基本方程,以有限元法对基本方程进行离散.以直立堤前完全立波和越浪条件的不完全立波、斜坡堤前的不完全立波等典型情况的模拟为基础,组合这些典型情况模拟削角堤前不完全立波运动,探讨其堤前波面、波节点及波腹点位置、波压力等特性变化.     

A SPH numerical wave basin for modeling wave-structure interactions

Based on a parallel SPH-LES model, a three dimensional numerical wave basin is developed to study wave interaction with coastal structures. The OpenMP programming technology combined with an existing MPI program contained in the parallel version of SPHysics code has been implemented to enable the simulation of hundred millions of particles running on a computer cluster. As part of the numerical basin development work an active absorbing wave maker and a sponge layer are introduced. The dynamic boundary conditions are also corrected to reduce the spurious effects. Wave generation and propagation in the numerical wave basin is first tested and confirmed with analytical results. Then, the model is applied to simulate wave interactions with a vertical breakwater and a vertical cylinder in order to further assess the capability of the numerical wave basin. The predicted free surface elevation near the vertical breakwater is compared with the experimental data while the horizontal forces and overturning moments acting on the vertical cylinder are verified with the analytical results. In all these cases the model results show excellent agreement with the validation data.     

An Explicit High Resolution Scheme for Nonlinear Shallow Water Equations

The present study develops a numerical model of the two-dimensional fully nonlinear shallow water equations （NSWE） for the wave run-up on a beach. The finite volume method （FVM） is used to solve the equations, and a second-order explicit scheme is developed to improve the computation efficiency. The numerical fluxes are obtained by the two dimensional Roe＇ s flux function to overcome the errors caused by the use of one dimensional fluxes in dimension splitting methods. The high-resolution Godunov-type TVD upwind scheme is employed and a second-order accuracy is achieved based on monotonic upstream schemes for conservation laws （MUSCL） variable extrapolation; a nonlinear limiter is applied to prevent unwanted spurious oscillation. A simple but efficient technique is adopted to deal with the moving shoreline boundary. The verification of the solution technique is carried out by comparing the model output with documented results and it shows that the solution technique is robust.     

A numerical tidal and storm surge model of the North Sea

Abstract

A two‐dimensional, vertically integrated, numerical model has been used to compute the M₂ and S₂ tidal components for the North Sea. The model has also been applied for simulation of two storm surge situations. The results of the computations are compared with sea level observations from coastal stations, and for one of the simulations a comparison with current measurements is also performed.     

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

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

Numerical modeling of tides in the Gulf of Suez,Egypt

A vertically integrated two‐dimensional (2‐D) and a five‐layer three‐dimensional (3‐D) numerical models were developed to compute the tides in the Gulf of Suez, Egypt. The computational grid used to schematize the Gulf has a horizontal resolution of 3 × 3 km and consists of a lattice of 23 × 100 points in the 2‐D model and five such lattices in the 3‐D model. Both 2‐D and 3‐D simulations clearly revealed the Kelvin wave nature of the tide with partial reflection. The M₂ simulations also showed a strong tidal signature in the southern sector as compared to the northern part. For the 3‐D simulations, the horizontal and vertical eddy viscosity coefficients and the coefficient of bottom friction were respectively set to 3 × 10⁶ cm²/s, 200 cm²/s, and 0.001. The tidal range decreases from the entrance of the Gulf of Suez toward the Bank of Tor where it reaches a small value and then increases again to about 1.5 m at Suez. A difference of six hours exists between the times of high water at the southern and northern ends of the Gulf. Quantitatively the agreement between the observed and computed tide is better in the 2‐D simulation than in the 3‐D simulation. However, the counterclockwise gyres in the top three layers of the 3‐D model (upper 30 m) during slack water following the flood tide do not appear in the 2‐D simulation.     

徒骇河感潮河段冲淤变化的数值模拟

研究天然河流,特别是感潮河口段河流的泥沙冲淤问题时,由于其复杂性和多变性,使得数学模型成为一种重要研究手段。本文采用天然冲积河道中缓变非恒定一维水流泥沙方程组建立了研究徒骇河感潮河段泥沙冲淤变化的数学模型。数值模拟时对河道进行了梯形断面概化。通过过河道子测站的实测资料对数值计算结果进行验证,二者吻合较好,在此基础上对徒骇河上游坝上处建闸后该河段在各级径流量影响下的冲淤变化趋势作出了模拟。     

NUMERICAL SIMULATION FOR THE SEASONAL MOVEMENT OF SUBTROPICAL HIGH RIDGE LINE Ⅰ. HEATING EFFECT OF SOLAR RADIATION

Jn this paper, the numerical simulation of the latitude position variation of subtropical high ridge line is made by using the two dimensional zonal average models. When the solar radiation heating is considered in the model, the jumping process of the subtropical high ridge takes place during its seasonal northward movement. There would be no such a jumping process if the radiation heating is absent. So, it may be concluded that the basic character of seasonal movement of subtropical high depends mainly on the solar radiation heating.     

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

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

渤海潮波系统的数值模拟

利用二维非线性潮波方程组,讨论了渤黄海主要分潮(全日潮、半日潮及浅水分潮) 数值模拟中的有关问题。数值模拟中同时考虑了4个主要分潮(M2,S2,K1,O1)和两个浅水分潮(M4,MS4)。分析表明,在渤黄海潮波系统数值模拟中,稳定后选取14 d的数值模拟结果进行调和分析能够取得最佳(最合理)的调和分析结果。计算出调和常数的模拟值与实测值之差的绝对平均值:M2分潮的振幅差为4cm,迟角差为3.3°,S2分潮的振幅差为2cm,迟角差为4.2°,K1 分潮的振幅差为1cm,迟角差为3.7°,O1分潮的振幅差为2 cm,迟角差为5.5°。实验结果较好地体现了渤黄海潮波系统的特征。     

A dynamic response analysis of tension leg platforms including hydrodynamic interaction in regular waves

A numerical procedure is described for predicting the motion and structural responses of tension leg platforms (TLPs) in waves. The developed numerical approach, in a TLP is assumed to be flexible instead of rigid, is based on a combination of the three dimensional source distribution method and the finite-element method. The hydrodynamic interactions among TLP members, such as columns and pontoons, are included in the motion and structural response analysis. Numerical results are compared with the experimental and numerical ones. The results of comparison confirmed the validity of the proposed approach.