数理模型在长江口水动力过程的研究中的应用.pdfVIP

数理模型在长江口水动力过程研究中的应用 袁政涛 李祥 邢少佳 上海市青浦区水务局水文勘测队 上海201700 我国采用数理模型研究长江口及河口工程的可行性已愈加普遍。从总体上而言，进入80年代后，除继续应 用推广一维潮流模型外，大多已采用二维潮流数学模型，并按需要配以泥沙、温度、盐度和营养盐等物质输移模型， 三维模型也已逐步展开。本文从长江口潮流场的数值模拟方法进行了较为系统的归纳总结，较全面地论述了我国20 多年来在长江口数值模拟研究地进展，认为从80年代以来，我国将数学模型应用于河口研究已取得显著进展，它已 成为研究河口的有力工具。 数 学模 型 ;长 江 口;水 动力 垂向仃坐标；动量和标量的平流项可采用TvD格式； 长江口南槽各参量横向差异较小．潮流主要呈 25阶湍流闭合模型求解湍流系数；水平扩散系数基 往复流形态。在该区域，横向积分的垂向二堆模式显 然是合理的。基于潮流有限差分方程，黄世昌等【4l 于Smagorinsky参数化方法；分裂算子技术——内模( 斜压)外模(正压)模态分开计算；考虑热力学因素。 2．2ECOMSED模式 黏滞系数。同时考虑由盐度差引起的密度异重流，采 用Munk．At)【te瑚n的衰减函数来模拟长江口南槽的 ECOMSED模式是目前国际上较为先进的水流、 泥沙、水质模型之一，具有计算三维水流、水质动床泥 垂向二维结构。模拟结果揭示了南槽既是涨潮优势 沙等诸多功能。该模式主要由水动力模块、泥沙输运模 流的转换地带，又是最大流速区。时钟等l呓003年应 块、波浪模块、热通量模块和水质跟踪模块组成，各模 用变网格有限元方法，采用任意四边形等参单元，离 块之间相互耦合。其中水动力模块是一个基于静压和 向二维数学模型．对长江口北槽水域两个观测站的 Boussi啪8q近似的三维斜压原始方程模式；自由海表 面。水平网格采用的是曲线正交坐标系统；垂直方向采 潮流水位、流速垂线分布和整个北槽潮流流速纵向 用R坐标．变量空间配置Arakawa_C格式。2．5阶湍流分布进行了模拟。潮流水位、流速垂线分布的模拟值 闭合模型求解垂向混合系数和扩散系数。水平湍流粘 分别与观测站的实测值可以较为吻合。 3．3三维数学模型 滞和扩散系数基于Smgorinsky参数化方法计算。 2．3 FVCOM模式 易家豪、叶雪祥1983年。赵士清1985年，宋元 为研究近海和河口环流，美国麻州大学提出了一 平等1990年，徐贵泉、褚君达1992年。刘子龙等 个三角形网格、有限体积、三维(3D)原始方程组海洋 1996年对长江口潮流三维数学模型进行了研究。根 VolumeCo鹅t锄dOce锄Mod．据潮汐非恒定流的三维运动方程和连续方程。运用 模式(FVCOM，Finite e1)。该模式最大的特色是结合了有限元法易拟合边 有限差分和有限元相结合的方法，易家豪、叶雪祥 界、局部加密的优点和有限差分法便于离散计算海洋 1983年闱求解了长江河口潮汐河道中的三维流速分 原始方程组的优点。有限元法采用三角形网格，给出线 布和水位随时间的变化过程。求解时。平面上运用了 性无关的基函数，求其待定系数，特点是三角形网格易 Galerkin方法．采用具有三角形单元的有限元公式。 拟合边界、局部加密；而有限差分法直接离散差分海洋 垂直方向和时间上均采用有限差分公式。在垂直方 海洋原始方程组．特点是动力学基础明确、差分直观、 向引进了无尺度坐标。从而解决了自由水面的计算。 计算高效。FVCOM结合两者的优点，数值计算采用方此外，考虑了柯氏力的作用，并将河床阻力作为边界 程积分形式和更好的计算格式。使动量、能量和质量具 条件