有限体积法及在近岸潮流计算中的应用的研究.pdfVIP

有限体积法及其在近岸潮流计算中的应用研究 捅 要 目前广泛应用于潮流计算的二维潮流方程属浅水方程。齐次浅水方程在数 学上为拟线性双曲型偏微分方程组。浅水方程在计算机技术发展的推动下已发展 了诸多有效的数值计算方法，如有限差分法(FDM)、特征法(MOC)、有限单元法 并取得了巨大成功。有学者很早就发现齐次浅水方程组与正压气流欧拉方程组在 数学上的一致性。随着空气动力学的发展完善，将空气动力学的理论成果和数值 算法应用于计算水动力学已进入实践阶段。本文对有限体积法的研究即是本着这 个思路。有限体积法用于积分形式的守恒律方程。它能结合有限差分法与有限单 元法的优点，如：计算效率高；可应用于无结构网格(有限元网格)：具有激波 捕捉性(shock—capturingproperty)，可以有效地处理间断问题。 本文从探讨一维和二维浅水方程的特征理论和性质出发，详细阐述了应用于 浅水计算的有限体积法一阶及二阶格式(MUSCL途径)的原理及应用方法。应用有 限体积法求解浅水方程的关键问题为求解网元界面处的一维黎曼方程，即确定网 元边界上的数值通量。精确求解黎曼方程非常繁杂且不必要，所以用有足够精度 且简便的近似黎曼解来代替精确的黎曼解。求解近似黎曼解的方法主要有：① Osher格式；②通量向量分裂格式(Fvs)；③通量差分裂格式(FDS)等。 本文还以自动生成(Delaunay原理)及编号的三角形网格为例，详细阐述 了无结构网格的应用方法。本文采用规则矩形网格编制了模拟一维及二维溃坝问 题的测试程序，以检验文中所述方法的正确性和处理间断问题的有效性。最后本 文采用无结构三角形网格并采用有限体积法二阶格式编制了模拟大亚湾潮流运 动的模拟程序。该模型的模拟结果表明了有限体积法用于潮流计算的有效性，其 实现方法也充分体现了有限体积法的上述优势。 关键字：有限体积法，浅水方程，潮流，数值模拟 FiniteVolumeMethodandIts of Application TidaIFlowSimulation Abstract The 2d of flOWsimulationtotheshallow used tidal widely equations belong water turbulentand effectsore shallowwater type．Whca dispersive neglected，the a offirstorder difierential become equations quasi．1inearsystem hyperbolicpartial lotof havebeen the ofthe science．amethods equations．Asdevelopmentcomputer usedto theshallowwater most four solve equationsnumerically．Thesignificant methodsofthemare：FinitcDifference ofCharacteristics Method(FDM)．Method Element FiniteVolume (MOC)。Finite Method(FEM)an