紊流数值模拟中的图形动态显示

目前 ,紊流数值模拟程序几乎都是ＦＯＲＴＲＡＮ语言编写的 ,由于用此语言编写图形处理程序的能力较弱 ,所以现今绝大部分数值计算均是先计算出结果 ,然后将其输入图形处理软件 ,再进行显示和输出 ,这样一来 ,周期长且效率往往不高。这是因为数值计算的成果很多 ,仅凭某种法则难以判断计算是否成功 ,有时看似成功的计算 ,当进行图形处理后 ,发现计算成果与实际不符。其原因可能是程序中的某个地方有误 ,或者法则本身有误。这种有误的程序有时也能得出收敛的结果 ,但其结果显然与实际相差较远 ,以致编程者不得不重新检查程序和重新计算。对于…     

台湾东侧黑潮区域海水的T-σt关系及其在动力计算中的应用

迄今为止的调查研究表明,海洋中的海流,特别是那些流速较大的强流:大都以地转流为其主要成份。因此,可以根据地转关系,应用观测所得的海水质量分布(即温、盐度或密度分布)计算出与实际较为相符的流速流量。这种过程就是海流的动力计算。 但是,由于温度资料比盐度资栂精确且易得,因此,人们在根据海水质量分布计算海流时,曾进行了种种简化,试图只根据温度分布及变化来研究海流。这大致可分为下列两方面： 其一,认为海水密度分布主要决定于温度,如果仅欲求得海流分布变化的粗略结果,可以不考虑盐度,直接根据单一的温度分布来计算流速流量。近年来,国内外所开展的一些工作,已为这种简化提供了初步的但较有力的证据。例如, Bryden(1974)通过对中大洋动力实验(MODE)区域两个锚锭站的温度和海流实测资料的分析,得出了水平流速的垂直切変和温度水平梯度的垂直积分之间颇为一致的关系。管秉贤在研究通过苏澳一与那国岛断面的黑潮的变异时,找到了100-200米水层处的温度水平梯度最大值与流量之间较好的关系。 其二,为了求得较为精确的结果,引入了平均的T-S关系曲线来考虑盐度分布对密度的影响。例如, Emery和Wert(1976)曾利用太平洋上每10°方格海域中的平均T-S曲线,对大量的温度断面进行了动力高度计算。所得结果,除了个別区域外,与用温盐资料算得的动力高度相比,误差都比较小。另外,在中大洋动力实验中,对许多只有温度而没有盐度资料的观测站,也是应用该海区的平均T-S曲线进行动力计算的。 本文应用国际黑潮合作调查(CSK)资料,对台湾东侧黑潮区域的水温(T)和条件密度(σt)进行比较,得出了这一区域的T-σt的相关关系式。运用这一关系式对该海区的温度断面进行动力计算,所得流速与用经典的动力计算方法得到的地转流速相比,其方向和量级都是一致的。另外,对东海的黑潮区域的海水,也用本关系式进行了讨论。     

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

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

渤海、黄海强温跃层数值预报模式

文中对前人关于水温垂直结构计算、模拟和预报研究成果进行了评价，并在此基础上建立了考虑海面吸收辐射和透射辐射、地形、风、海流、界面摩擦及其引起的混合作用的强温跃层三维数值预报模式。该模式运行的驱动量仅为风和气温场资料，因而具有较好的实用性；此外，从试报的结果看，效果是令人满意的。     

贡湖及其相关水域风生流模拟研究

在太湖风生流的数值模拟研究的基础上，设计模拟贡湖及其相关水域风生流的套网格数值模式。应用该模式研究不同风情下，梅梁湖，太湖主体湖水对贡湖的影响。     

南黄海潮流和潮余流的数值计算

一、引言南黄海是一个平均深度不超过45米的陆架浅海。其独有的地理条件,使该海域潮运动独具特征。早从本世纪三十年代开始,日本学者就已对黄海的潮汐进行了研究。近年来,我国、苏联、日本和南朝鲜等国学者又运用数值计算等方法进一步探讨了该海域的潮汐、潮流运动。然而,这些研究多偏重于揭示海域潮波运动的规律。对潮流的计算基本上未经实测资料校核。尤其是,对于由潮流的非线性效应所产生的潮汐余流研究甚少。为了更     

用流函数理论计算缓坡上波浪的破碎位置

本文通过理论计算及实验验证发现,在计算缓坡上波浪的变形、特别是确定波浪破碎点时,流函数波浪理论明显优越于线性波液理论。     

复杂网络河的数值模拟

依据复杂网络河的特点及其内在联系，本文建立了其特征矩阵及由Newton-Raphson格式表述的隐式圣－维南差分方程组的系数矩阵，经输入边界条件及解方程组，可得到不同各支段上，下游端点上的水位和流量值。     

Modified PIC Method for Sea Ice Dynamics

The sea ice cover displays various dynamical characteristics such as breakup, rafting, and ridging under external forces. To model the ice dynamic process accurately, the effective numerical modeling method should be established. In this paper, a modified particle-in-cell （PIC） method for sea ice dynamics is developed coupling the finite difference （FD） method and smoothed particle hydrodynamics （SPH）. In this method, the ice cover is first discretized into a series of lagrangian ice particles which have their own sizes, thicknesses, concentrations and velocities. The ice thickness and concentration at Eulerian grid positions are obtained by interpolation with the Gaussian function from their surrounding ice particles. The momentum of ice cover is solved with FD approach to obtain the Eulerian cell velocity, which is used to estimate the ice particle velocity with the Gaussian function also. The thickness and concentration of ice particles are adjnsted with particle mass density and smooth length, which are adjusted with the redistribution of ice particles. With the above modified PIC method, numerical simulations for ice motion in an idealized rectangular basin and the ice dynamics in the Bohai Sea are carried out. These simulations show that this modified PIC method is applicable to sea ice dynamics simulation.