三维浅海潮流模型及其混合有限分析解法研究.pdfVIP

三维浅海潮流模型及其混合 有限分析解法 姜国强1 陶建华1 槐文信2 (1．天津大学机械工程学院天津300072I 2．武汉大学水利水电学院，湖北武汉430072) 摘要本文以具有自由表面的三维非线性Na啊er-stokes方程组为基本控制方程，垂向 采用。坐标变换后得到三维浅海潮流数学模型，并采用干湿判别动边界技术处理潮滩问题。 数值求解方法采用具有自动迎风、精度高及稳定性好的混合有限分析法，最后，采用三个算例 验证了本文模型及其混合有限分析解法模拟潮流、风生流及处理湖滩动边界问题的有效性和 正确性。 关键词 三维数学模型，混合有限分析解法，干湿判别法 1 引言 随着计算机技术和数值模拟方法的发展，以及在研究泥沙底沙起动，疏浚抛泥、河口盐水 楔入侵等这些工程问题时的实际需要，自20世纪70年代中期浅海三维水动力学数值模拟的 研究开始逐步发展起来。 对于各种三维模式，按照不同垂向分层的方法有：由Leendertse[1]等人提出的固定分层模 式，它是在垂向上以分层平均的方法筒略地将水体分为若干层，在每层中将流动基本控制方程 沿垂向积分，从而将每层水体转化成为平面二维问题求解；B1un南erg[2]等人通过将垂向坐标 无量纲化得出相对分层模式，该模式使得整个计算区域具有相同的垂向计算网格数，对浅海计 算区域具有更高的分辨率，并且层层之间存在质量交换，对流动的模拟更符合实际情况； Heaps[3]建立的浅海三维水动力学的谱方法模式也是颇具特色的一种方法，该方法将水平流 速的垂向变化以本征方法做成谱展式，把一个三维问题变成一系列的二维问题。 对流扩散方程的混合有限分析求解方法具有自动迎风、精度高及稳定性好的特点H]。为 此，本文通过混合有限分析法求解垂向采用d坐标变换后得到三维浅海潮流数学模型，潮滩动 边界采用干湿判别方法进行。并通过经典算例验证了本文模型及其混合有限分析解法模拟风 生流、潮流及处理潮滩动边界问题的有效性和正确性。 2 口坐标系下的三维浅海潮流数学模型 2．1基本方程及定解条件 ’ 于三维潮波流动水平尺度一般远大于垂向尺度，丽虽水质点的垂向加速度远小于重力加速度， 依据静压假定垂直方向的动量方程可退化为静压方程。 对于河口或海湾的水流运动，在其物理区域内，存在着自由水面随时间的不断变化，同时 实际情况中，水底地形也沿水平方向呈不规则变化，使得水深在空间和时间上不断变化。为 此，引入新的坐标系，通过对基本方程进行。坐标变换，将自由水面和水底变换成为口坐标平 面，使整个计算水域具有相同的网格单元数目。它与直角坐标系的变换关系为 ￡2 2 z；z ￡=z7；z=z7；y=y7；z=P(z7，了7，口，￡7)=口(五+咖+叩z；y=y。；z=P(z’，了‘，口，￡+)=口(h十叩)十叩 那么，在a坐标系下，各项微分有如下表达式： a aP， a —a a午l a一． 一==——’…’ a￡ a￡’ ad az az’ P。’ a口9。’ 一：=’_—，一● a a aP’， a a 1 a了 a9 a6午o’ az ao甲， 一==———_’一·—一==。一—一 式中算子 红=雾；％=磬；纷一多；吼=筹 由此得到d坐标系下的基本方程 连续方程： 。 ‘ 。 挈+宅丝+宅丝+鸯孕：o