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

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

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.     

A Study on the Spectral Form of Nearshore Water Waves

- The spectral form of wind waves is investigated based on the ocean wave data observed at three nearshore stations of Taiwan. In this study, the generalized forms of Pierson-Moskowitz spectrum and JONSWAP spectrum are used to describe the local wave spectrum by selecting suitable spectral form parameters. It is shown that, at a specific site, the similarity of wave spectral form exists. Thus it is possible to use a representative spectral form for a given nearshore region to describe the wave spectrum at this nearshore. On the other hand, the effects of relative water depth on spectral form are examined. The feasibility of two spectral models in finite water depth is evaluated by using the same field wave data.     

波浪作用下孔隙海床-管线动力相互作用分析

波浪作用下海床中的孔隙水压力与有效应力是影响海底管线稳定性的主要因素。然而,在目前的海床响应分析中一般将管线假定为刚性,并不能合理地考虑海床与管线之间的相互作用效应,同时也没有考虑土体和管线加速度对海床动力响应的惯性影响,从而无法确定由此所引起的管线内应力。为此考虑管线的柔性,分别采用饱和孔隙介质的Biot动力固结理论和弹性动力学理论列出了海床与管线的控制方程,进而采用摩擦接触理论考虑海床与管线之间的相互作用效应,基于有限元方法建立了海床-管线相互作用的计算模型及其数值算法。通过变动参数对比计算讨论了管线几何尺寸、海床土性参数对波浪所引起的管线周围海床孔隙水压力和管线内应力的影响。     

复合加载条件下吸力式沉箱基础承载特性数值分析

吸力式沉箱基础的承载特性是海洋工程结构设施建造与设计中的一个关键问题。这种新型的深水海洋基础型式,通常承受竖向上拔荷载与水平荷载的共同作用,其工作性能与设计理论远远不能满足工程实践的需要。本文采用有限元分析方法对吸力式沉箱基础的极限承载特性进行数值计算。以大型通用有限元分析软件ABAQUS为平台,通过二次开发,数值实现了Swipe试验加载方法和固定位移比分析方法,针对不同的沉箱长径比、土的强度折减系数,探讨了沉箱基础在垂直上拔荷载和水平荷载单调联合作用下的极限承载力,通过对不同荷载组合的数值计算构造了复合加载条件下沉箱基础破坏包络面。     

钛型吸附剂的表面特性

用自装的BET容量法装置测定了几种粒状钛型吸附剂的比表面及其孔径分布。测定是在常温减压下进行的。结果表明:钛型吸附剂在加铀海水中的吸铀量与样品的比表面及孔径分布有一定的内在联系。这种表面性质的研究对吸附剂的筛选有重要的指导意义。     

有限尺度李雅普诺夫指数与海表温度梯度相关性的初步分析

为了探究海表温度和海面高度之间的瞬时相关性,介绍了一种卫星测高数据的拉格朗日分析指数--有限尺度李雅普诺夫指数(Finite Size Lyapunov Exponent,FSLE),以黑潮延伸体区域的涡旋和南大西洋的亚南极锋为例,通过对观测、模式结果和融合产品结果的分析,探讨了该指数与海表温度梯度(Sea Surface Temperature Gradient,SSTG)之间的相关性。比较FSLE图像和SSTG图像发现,FSLE与SSTG均呈丝状结构,对海洋表层水体结构描述具有一致性,尤其在温度梯度大和地转流强的区域更为一致。二者的一致性要远好于其他常用方式,比如全流速、OW参数涡旋识别方法和Winding-Angle涡旋识别方法。不同区域FSLE与SSTG之间的相关性表现不同,黑潮延伸体区域相关系数存在显著的季节变化,而南大西洋亚南极锋区域季节内变化突出。     

非常规管节点疲劳寿命分析与计算

对非常规管节点的疲劳寿命进行了分析研究，用SESAM计算疲劳载荷，应用精细有限元分析计算热点应力，用规范的S—N曲线计算管节点的疲劳寿命。通过实例计算表明，这些处理对非常规管节点疲劳寿命分析与计算是很有效的。     

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

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

Computation of Viscous Uniform and Shear Flow over A Circular Cylinder by A Finite Element Method

The incompressible viscous uniform and shear flow past a circular cylinder is studied. The two-dimensional Navier-Stokes equations are solved by a finite element method. The governing equations are discretized by a weighted residual method in space. The stable three-step scheme is applied to the momentum equations in the time integration. The numerical model is firstly applied to the computation of the lid-driven cavity flow for its validation. The computed results agree well with the measured data and other numerical results. Then, it is used to simulate the viscous uniform and shear flow over a circular cylinder for Reynolds numbers from lO0 to lO00. The transient time interval before the vortex shedding occurs is shortened considerably by introduction of artificial perturbation. The computed Strouhal number, drag and lift coefficients agree well with the experimental data. The computation shows that the finite element model can be successfully applied to the viscous flow problem.