一种有限体积法三维水流悬沙数学模型的研究.pdf

摘要 河口、海岸潮流泥沙运动一般具有三维特征，故需要三维数学模型才能合理 反应其运动过程。本文在对河口、海岸潮流泥沙数学模型的研究现状进行综合分 析的基础上，结合已有工作，探讨了一种基于无结构三角网格有限体积法的三维 浅水流动及悬沙运动数学模型。水流模型的主要特征有：(1)采用模式分裂技术， 内、外模分别计算，以提高求解效率；(2)采用k一占二方程紊动模型计算垂向紊 动粘滞系数，物理意义更加明确；(3)在计算方法上，变量垂向分布用隐格式求 解，以增强数值稳定性；(4)平面上采用无结构三角形网格，以适应天然水流不 规则水陆边界；(5)垂向上采用盯坐标系，以简化表面和底部边界条件；(6)采用 有限体积法离散方程，并采用性能优越的Ostler格式计算单元交界面通量，使水 流模型具有很好的守恒性、逆风性、单调性和高分辨率地捕提间断的能力，并且 在边界的处理上不需引入冗余的数值边界条件；(7)用干湿网格法决定移动边界 variation 位置，使模型适应缓坡地形情况。悬沙输运方程采用TVD(total diminishing)格式进行离散，并引入高阶数值黏性项，使模型既可高精度捕捉间断， 又有良好的稳定性。 将水流模型应用到黄河贵德水文站下游某河段，模拟的流速、水位与实测资 料吻合良好；流量沿程分布结果表明模型具有良好的物理量守恒性；给出了流速 的垂向分布、紊动能量横断面分布，力河道三维水流结构分析提供了依据。 在一个理想矩形水槽中，运用水流模型进行了潮流计算，结果表明模型具有 合理的动边界模拟效果。计算潮位及涨落潮流场亦是合理的，为以后模拟真实潮 波运动奠定了基础。 在一个理想矩形水槽中对三维悬沙悬扬、扩散进行了模拟，结果令人满意， 表明该模型提出的三维水流悬沙扩散数值计算方法是可靠的。 关键词： 有限体积法Osher格式七一占紊动模型三维浅水模型泥沙运动 ABSTRACT in Themotionoftidalcurrentandsedimentestuarineandcoastalareasis be characterized3D shouldsimulated3Dmathematicmodel． by variation,which by availablemathematical Basedonreviewofthe modelsoftidalflowandsediment the andcoastal 3Dnumerical water for estuarine region，a modelforshallow transport isestablishedonthefinite flowand sediment based volume suspended transport methodandunstructured is characterized triangulargad，which by：(1)the solution is inthe theinnermode applied mode—splittingtechnique model，namely mode the