含有紊动黏性项浅水方程的数值求解

紊动黏性项对水流方程数值求解的稳定性至关重要.采用基于无结构网格的有限体积法对方程进行离散,对流项的界面通量采用Roe格式的近似Riemann解计算,紊动黏性项采用单元交界面的平均值估算,底坡源项则采用特征分解法处理.应用该模式模拟了圆柱绕流,模拟结果与实际情形较相符.对港珠澳大桥的海中人工岛水域的潮流运动进行了有效模拟,结果表明该模式能够较好地模拟出水建筑物附近水域复杂的水流结构.

一种新的浅水方程黏性项计算格式的研究与应用

用浅水方程模拟某些紊流特性很强的流动时,紊动黏性项对数值解的稳定性和合理性至关重要。为提高浅水方程中黏性项的计算精度,首先对前人提出的两种计算格式进行精度分析,在此基础上提出了一种具有更高精度的计算格式。这种方法从有限体积法的特点出发,用单元形心的流速直接计算流速梯度,比使用奥高公式具有更高的精度。通过3个典型算例的比较,验证了该方法具有更好的数值稳定性,同时也证明了此格式能比较合理地反映平均水深流动问题中的一些紊流现象,在工程上具有较好的应用价值。