河口三维流动数学模型研究评述.pdfVIP

河口三维流动数学模型研究评述术 l-l-嗣●，，1，{，j 刘桦 何友声 (上海交通大学，上海200030) 摘要 本文回顾了近年来国内外河口与海岸三维盛垫盎兰数学模型研究的新进展t着 重讨论r三维水动力学数学模型研究中的若干基率问题．本文认为，无论是正交苗卡尔 坐标系还是一般曲线坐标系，河口三维流动问题的基本方程己经解决；尽管已有一些成 功的有限差分模型，但是寻求更确效的高精度数值格式仍值得我们去努力．河口大觏摸 动边界和漫滩水流的模{c【有待作迸～步研究。 关链词河口，三维数学模型，{平述 1 前言 河口与海岸环境中的流动具有明显的空间三维特性，从理论上讲，可以从PANS方 N．S ： 程(ReynoldAveragedequations)出发建立三维水动力学数学模型。尽管有各种湍流 模式或大涡模拟(LES)技术可用，但是由于受计算机容量和速度的限制，对这类大水体 的运动进行精细的三维数值模拟在目前看来仍是困难的．考虑到近海潮流的水平空间尺 度远大于垂直空间尺度，因此可以看作一种准平行流动，这时水质点运动的垂直加速度 比重力加速度太得多，垂直方向的流体运动方程退化为静压方程，即所谓的静压假定。 J 这样，我们可以用水位变量f(x，y，f)代替压力变量，使得需求解的数学问题得到简化。 Leendertse f1973)的工作具有开创性，他在垂直方向采用固定分层法，即将计算水域划分 为固定的多层．在每层中沿水深积分使之成为二维问题，并用ADI格式进行数值离散． Kim等(1994)应厢固定分层方法建立了海湾三维潮流、盐度模型。为了更好地模拟河 床地形的变化．人们将Ph[1ips({957)提出的a坐标变换应用到河口与海岸三维模型中。 为了较严格地确定涡粘性系数和扩散系数，湍流模式理论也在潮流模型中得到了应用。 OceanModel)。该 于三维数值模型的开发与应用研究，其代表性软件为POM(Princeton 模型在垂直方向采用盯坐标系，在水平方向采用正交曲线坐标系及ArakawaC差分网格， 但数值求解过程仍在物理平面上进行。为了提高计算效率，POM采用丁模态分离法，将 带自由面的三维流动问题分成表面波的传播问题(外模态)和内波的传播问题(内模态)， and 对外模态仍用显式差分．有关POM的详细介绍可在BlumbergMellor(1983)， Me!lor(1996)@找到。Sheng(1987J建立了一般曲线坐标下的三维水动力学模型(CH3D)。 该模型也采用盯坐标系，水平方向的运动方程采用水平流速矢量的逆变分量来表示。Wang 三维潮流数学模型的研究进展作了回顾。 为了满足我国大型河口研究和重大整治工程建设的需要，国内学者近年来对河El三 ·84· 海域的三维潮流进行了数值模拟。韩国其等(1989)采用O-坐标系和算子分裂技术建立 盯坐标系和模态分裂法对渤海湾的三维潮流作了数值模拟。宋志尧等(1998)基于模态 分裂法和ADI格式建立了三维潮流的计算模式，并应用于海岸辐射沙洲的潮流场分析。 建立了计及河口密度分层效应的三维潮流、盐度数学模型，垂直方向为d坐标系，水平 方向为直角笛卡儿坐标系，经模态分离后，外模态用改进的ADI法求解。为了合理确定 垂直涡粘性和涡扩散系数，应用了与POM类似的非均质流体的两方程湍流模式。该模型 应用大涡模拟技术建立了近海三维流动数学模型，较好地反映了台风过境时中国南海的 风吹流流场。沈焕庭(1997)曾对我国的河口数学模型作了回顾，对三维模型也有