浅水波方程的高阶保正Well-Balanced ADER间断Galerkin格式.pdf

Advances in Applied Mathematics 应用数学进展应用数学进展,2023,12(8),3728-3743 Published Online August 2023 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2023.128367 文章引用文章引用:周翔宇,张志壮,钱守国,李刚.浅水波方程的高阶保正 Well-Balanced ADER 间断 Galerkin 格式J.应用数学进展,2023,12(8):3728-3743.DOI:10.12677/aam.2023.128367 浅水波方程的高阶保正浅水波方程的高阶保正Well-Balanced ADER间间断断Galerkin格式格式 周翔宇，张志壮，钱守国，李周翔宇，张志壮，钱守国，李 刚刚*青岛大学数学与统计学院，山东 青岛 收稿日期：2023年7月26日；录用日期：2023年8月18日；发布日期：2023年8月24日 摘摘 要要 本文针对具有不规则几何形状和非平坦底地形的浅水本文针对具有不规则几何形状和非平坦底地形的浅水波波方程，引入了保正高阶方程，引入了保正高阶ADER间断间断Galerkin方法，方法，该方法能准确地保持静水的稳态。为了该方法能准确地保持静水的稳态。为了满足满足well-balanced的性质，我们提出了的性质，我们提出了well-balanced的数值通的数值通量，并基于分解算法将数值解分解为两部分，构造了一种新的源项近似，并相应地将源项近似分解为两量，并基于分解算法将数值解分解为两部分，构造了一种新的源项近似，并相应地将源项近似分解为两部分。此外，部分。此外，还还引入了一个简单的保正限制器引入了一个简单的保正限制器，从而在干湿锋面附近提供高效和鲁棒从而在干湿锋面附近提供高效和鲁棒性性的模拟。大量的的模拟。大量的数值实验数值实验也表明，所得到的也表明，所得到的格式格式s能够准确地捕捉静止稳定状态下湖泊的小扰动，保持水面高度的非负性，能够准确地捕捉静止稳定状态下湖泊的小扰动，保持水面高度的非负性，同时保持光滑解的真正同时保持光滑解的真正高阶精度。高阶精度。关键词关键词 浅水波方程，浅水波方程，ADER方法，间断方法，间断Galerkin格式，保证格式，微分变换过程，全离散格式，保证格式，微分变换过程，全离散 High Order Well-Balanced ADER Discontinuous Galerkin Scheme for Shallow Water Wave Equations Xiangyu Zhou,Zhizhuang Zhang,Shouguo Qian,Gang Li*School of Mathematics and Statistics,Qingdao University,Qingdao Shandong Received:Jul.26th,2023;accepted:Aug.18th,2023;published:Aug.24th,2023 Abstract In this paper,we introduce positivity-preserving high order ADER discontinuous Galerkin Me-*通讯作者。周翔宇 等 DOI:10.12677/aam.2023.128367 3729 应用数学进展 thods for the shallow water equations with irregular geometry and a non-flat bottom topogra-phy,which maintain the still water steady state exactly.To achieve the well-balanced property,we propose the well-balanced numerical fluxes,and construct a novel source term approxima-tion by decomposing the numerical solutions into two parts based on decomposition algorithm,and resolving the approximation to the source term into two parts accordingly.Moreover,a simple positivity-preserving limiter is introduced in one dimension and then extended to two dimensions to provide efficient and robust simulations near the wetting and drying fronts.Ex-tensive one and two dimensional simulations also indicate that the resulting schemes enjoy the ability to accurately capture small perturbations to the lake at rest steady state,maintain the non-negativity of the water height,and keep the genuine high order accuracy for smooth solu-tions at the same time.Keywords Shallow Water Wave Equations,ADER Approach,Discontinuous Galerkin Scheme,Positivity-Preserving Method,Differential Transformation Procedure,Fully-Discrete Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 在本文中，我们主要的工作是构造非平坦底部地形上的浅水波方程(SWEs)的保正 well-balanced 性质的间断 Galerkin 格式，该方程具体形式如下1 2 3 4：()()220,1.2txxtxhhuhuhughghb+=+=(1)在上式中，h 为水深，u 为水深平均速度，b 为水底地形。g 是重力加速度常数。浅水波方程(1)属于双曲守恒律的一种，并且其静水稳态解由下式给出：0,const,uhb=+=(2)其中源项与通量梯度完全平衡。这种平衡律的数值模拟，包括非平坦底部地形上的非线性浅水波方程，有一个主要的问题：标准数值方法可能无法在离散水平上保持稳态下通量梯度和源项之间的平衡，并且可能在稳态附近产生伪振荡。Well-balanced 方法旨在克服这一问题，并且在稳态（或接近稳态）解决格式下表现良好。对于非线性 SWEs，已经有多种方法被提出并应用，其中包括5-12及其参考文献。在相关问题的数值模拟中，另一个经常遇到的问题是对少水或无水地区的润湿和干燥情况的处理。典型的问题包括溃坝问题、洪波和爬高现象。数值上，在计算过程中可能会产生负值，这可能会带来额外的困难。有许多现有的保正方法可以克服这一困难，关于这些方法我们参考了文献13-18。上述对浅水波方程的所有工作都是在有限差分法或有限体积法的基础上进行的。在过去的几十年中，高阶有限元 Galerkin(DG)方法在求解包含双曲型守恒律的偏微分方程方面得到了广泛的关注。近十多年来，高阶时间和空间任意导数方法在求解对流占优问题中得到了广泛的关注。Toro 和合作Open AccessOpen Access周翔宇 等 DOI:10.12677/aam.2023.128367 3730 应用数学进展 者最初在双曲守恒律的有限体积格式背景下提出了 ADER 格式19 20 21。格式19 20 21是单步、单级和全离散的。后来，Dumbser 和 Munz 22 23 24提出了一种采用单阶段 ADER 方法19 20 21的 DG 方法，用于时间离散化。随后，Dumbser 等人对该方法进行了一系列研究25 26 27 28。DG 方法25 26 27 28的关键是构造高阶数值通量。网格间状态以时间泰勒级数展开的形式表示，其中时间导数由使用控制 PDE 的空间导数表示。该程序被称为 Cauchy-Kowalewski 程序(在文献中也称为Lax-Wendroff 程序29)。此外，还需要求解网格间的广义黎曼问题，以获得状态本身及其空间导数。然而，Cauchy-Kowalewski 方法非常繁琐，尤其是在求解高维问题的时候。因此，CauchyKowalewski 程序的替代或简化是非常有意义的。值得注意的是，Dumbser 和 Munz 30提出了一种高效的 Cauchy-Ko-walewski 过程算法。最近，Dumbser 等人31，Duan 和 Tang 32采用局部连续时空 Galerkin 预测方法33来取代 DG 方法框架下的 Cauchy-Kovalewski 过程。Galerkin(DG)方法是一类以高阶间断分段多项式空间作为测试函数空间的有限元方法34 35 36。半离散DG方法(也称为RKDG方法)使用间断的分段多项式空间作为数值解和测试函数空间，并通过TVD Runge-Kutta 格式进行稳定的高阶时间积分37-42。Qian 43设计了几何形状不规则的明渠浅水波方程的高阶保正平衡方法。最近，ADER 方法由 Toro 和 Titarev 引入到有限体积格式中，该方法通过完全离散的一步公式获得空间和时间上任意高阶精度。原始版本的 ADER 方法在空间和时间上使用截断的 Taylor 展开，并通过Cauchy-Kovalevskaya(C-K)过程44将时间导数替换为空间导数，但是该过程不能处理刚性源项。为了克服这些问题，ADER-DG 方法用局部时空 DG 预测器取代 C-K 过程45。因此，单步公式只需要逐点计算通量和源项，这样就极大地减少了空间元素之间的模板依赖，并允许几乎完美的并行化。在本文中，我们构造了湖处于静止稳态的 SWEs 的高阶良好平衡的 ADER-DG 方法，并用微分变换代替 C-K 过程。为了达到良好的平衡性能，我们建议结合基于流体静力重建思想的数值通量，以及一种新的源项近似和分解算法。所得到的 DG 方法在保持湖泊静止稳态的 well-balanced 特性的同时，又保持了光滑解的高阶精度。本文结构如下：在第 2 节中，我们提出了保正平衡的 ADER-DG 方法。第 3 节给出了一些数值算例来验证所提出的 ADER-DG 方法的性能。最后，第 4 节给出结论。2.ADER-DG 方法方法 在本节中，我们构造了浅水波方程的高阶良好平衡的 ADER-DG 方法(1)。为了准确地保持静水稳态(2)，我们基于流体静力重建技术构建了 well-balanced 的数值通量，并提出了一种新的基于分解算法的源项近似。2.1.符号记法符号记法 首先将空间区间,a b离散成 N 个单元1122,1,2,jjjIxxjN+=。每个单元格的中心是()112212jjjxxx+=+，单元格jI的大小记为1122jjjxx+=，最大网格大小记为1maxjj N=。符号1,nnjjItt+=表示时空控制体积，1,nntt+为时间间隔。在这里，我们取近似解空间如下：()()(),:,jkkjVx tx tP=(3)其中()kjP表示j上 k 次以下的时空多项式集合。注意，kV中的函数允许在元件界面上具有间断。此周翔宇 等 DOI:10.12677/aam.2023.128367 3731 应用数学进展 外，我们使用()111222,2jjjvvv+=+表示函数v在元素界面12jx+处的算术平均值。2.2.微分变换过程微分变换过程 在本节中，我们将介绍一种基于将微分方程转化为其时空序列展开解的系数递推关系的方法，称为微分变换法(DTM)46。值得注意的是，该方法在评估 C-K 过程时更加简洁高效47。为方便起见，我们首先考虑一维平流方程()()00,0.txqqq xqx+=(4)将(),0q x的初始条件展开为泰勒多项式()()()0,0,0.xxkxjkqxq kxx=(5)接下来，我们假设解可以表示为每个时空域内的以下多项式()()()()()00,.txxtkknxtjjkkqx tq k kxxttx t=(6)微分变换q 的形式是这样的()(),1,0,1,2,.!xtxtnjkkxtxtkkxtx xt tq x tq kkkkkkxt+=(7)将(6)代入平流式(4)，可得递归关系如下：()()1,11,.1xxtxttkq k kq kkk+=+(8)因此，从式(5)中的初始微分变换(),0 xq k开始，反复重复式(8)，得到任意时刻的微分变换。因此，我们可以通过泰勒级数展开(6)得到平流方程(4)的解，即任意高阶解。随后，我们考虑了浅水波方程的微分变换方法(13)。我们将(),0U x的初始条件投影到分段多项式空间kV中，得到(),0Ux，其形式为()()()0,0,0.xxkkxjkUxU kxx=(9)与 Runge-Kutta DG 方法不同，我们通过 k 次的时空多项式在每个时空域j内寻找浅水波方程(13)的数值解，如下所示()()()()00,.txxtkkkknxtjkkUx tU k kxxtt=(10)与 2.2 节类似，我们将(10)代入(45)，得到浅水波方程(12)的递推关系如下：()()()()()()()()1231,11,1111,11,1,1,.111xxtxttxxxxtxtxtxttttkh k khukkkkkkhuk kG kkGkkGkkkkk+=+=+(11)其中 周翔宇 等 DOI:10.12677/aam.2023.128367 3732 应用数学进展 ()()()()()()()()()()()()()()()()()()11,1,0020033,01,002,1,0,2,1,0,0 xtxtxxtkkxtabxtrskkxtxtrskxttaxrkkaxtxtrskbxtbxtsGkkGr s Gkr ksgGkkh r s h kr ksGkkh r kGkrGkkhur shukr ksGkkh r s Gkr ksh=()()()03,11.xtkraxxxGkg kb k=+因此，将初始微分变换(),0 xU k和(),0,xxb kkk=代入(11)，得到微分变换(),xtU kk，并利用(10)得到在时空域j上的近似(),U x t。2.3.ADER-DG 方法的构建方法的构建 在本节中，我们将利用 2.2 节中的微分变换过程构建 SWEs(1)的 ADER-DG 方法。首先，我们将原控制方程(1)改写为()()()()()2220,11.22txxtxxhhuhuhughg bg hb xh+=+=+(12)我们将源项xghb按照48 49中的方法改写为动量方程中的等效形式()()212xxg bg hb b+。为了便于表示，我们引入了可以将(12)改写为紧致向量形式的符号()(),tUF US U+=(13)其中()()TT221,2Uh huF Uhu hugh=+和()()T210,2xxS Ug hb bgb=+分别表示保守变量、物理通量和几何源项。式(13)与任意空间测试函数()x的乘法并使用j上的部分积分产生以下弱形式()()()()()()()()()()()()()11122111221,d,d,d,d,d dd dfor1,2,.nnjjnnjjtnnjjIIttjjtU x txxU x txxF xttxF xttxF x txx tS Utxx tjN+=这里，为了便于表示，我们将物理通量写作一个二元函数()()(),F x tF U x t=。求解(13)的全离散 ADER-DG 方法定义如下：对于所有测试函数()()xkjxxx=且0,1,xkk=，从kV出发，解(),kUx tV满足 周翔宇 等 DOI:10.12677/aam.2023.128367 3733 应用数学进展 ()()()()()()()()()()()111122221,d,d,d dd d1,2,.jjjjnnjjjjIIUx txxUx txxFxFxFx txx tS Utxx tjN+=(14)这里用数值通量12jF+来近似网格之间物理通量的时间积分，即()112,dnntjtF xtt+。其中()()()()()111222,jjjF Utf UtUt+=数值通量与简单的 Lax-Friedrichs 通量一致()()()()()()()()()()111111222222,1,.2jjjjjjf UtUtF UtF UtUtUt+=+(15)其中()()11max,dnntntx IU x ttt+=.(16)其中()U为雅可比矩阵()F U的特征值，其具体形式为()()()11max,dnntntx Iux tghx ttt+=+.(17)2.4.平衡良好的数值通量平衡良好的数值通量 请注意，在前一小节中提出的 ADER-DG 方法(45)不具有良好的平衡特性。在本节中，我们提出了对数值通量的修正，目的是准确地保持稳态(2)。为了保持正性(如第 2.14 节所述)，我们还在时间步长nt上引入了新的网格边界值：()1111122222*,max 0,max,njjjjjhhbbb+=+(18)其中 1122*,*,jjhh+=在(2)中满足稳定性。现在，让我们讨论平衡的数值通量。在(15)定义的 Lax-Friedrichs 数值通量12jF+中，附加项()1122,nnjjUU+对数值黏度起到了作用，这是非线性守恒律所必需的。然而，它们可能破坏稳态的良好平衡特性。因此，我们将该通量(15)修改为()()()()()()()()()()()()()()()()()()()1111222211112222112211221212*,*,2222,*,1,11222jjjjjjjjjjjjjjhuthutf UtUthuthutg htg htututhhhut+=+()()1212*,.jjhut+(19)周翔宇 等 DOI:10.12677/aam.2023.128367 3734 应用数学进展 这里我们将通量项的第一个分量从()()12,jhut+修改为()()1122*,*,jjht ut+以便在第 2.14 节中解释保正性。2.5.源项近似源项近似 接下来，我们给出源项积分的高阶平衡近似，其形式为()()()221d dd dd d,2jjjxxSt v x tghbt b v x tgb v x t=+(20)其中2S表示源项的第二分量。我们可以进一步把这个积分分解为()()()()()()()()()()()()()()()()()()()()()111221112222d dd d1d dd d2d ddd1dd,2jjjjjnnjnnjxjxxjxjtxjjjtItxjjtISt v x tghbthbtb v x tghbt b v x tgb v x tghbthbtb v x tghbttbvxbvxbvxgtbvxbvxbvx+=+=+(21)其中项()()()()2,jb bhbthbt+分别用它们的数值近似值()()2,jb bhbt+代替。此外，将2,b b的边界值替换为它们在网格面处的平均值，记为 2,bb，与稳态时的数值通量(51)一致。2.6.斜率限制器斜率限制器 对于双曲守恒律系统，当解包含间断点时，ADER-DG 方法通常与斜率限制程序相结合，以压缩间断点附近可能出现的振荡。文献中有许多不同的斜率限制器选择，本文我们考虑经典的全变分有界(TVB)限制器50。标准限制器可能与良好平衡性质相冲突，我们根据51中提出的思想提出了一个良好平衡限制器程序。未知U上的标准 TVB 限制器包括两个步骤。第一步是根据网格平均值,1,jjUU和网格边界值,1 2,1 2,jjUU+来检查单元jI中是否需要任何限制。第二步是对该单元jI中的变量U应用 TVB 限制器。在提出的平衡斜率限制程序中，我们将首先根据单元平均值()(),1,jjhbhb+和()(),1 2,1 2,jjhbhb+来检查是否需要限制。如果单元jI被标记为需要限制，则实际的 TVB 限制器像往常一样应用于U。注意()hb+在静水稳态(2)时变为常数，因此在达到稳态时不施加限位器，保持了良好的平衡性质。2.7.均衡格式概要均衡格式概要 我们提出的浅水波方程(1)的良好平衡 ADER-DG 方法由式(45)给出，其中的数值通量定义在式(51)中，源项近似在上文中已经给出。综合前面各小节的结果，可以很直观地证明如下结果：命题命题 1 上述浅水波方程(1)的 ADER-DG 方法可以保持稳态解(2)的 well-balanced 性。证明 在稳态(2)时，我们有(),0.hbC u+=我们只需要证明包含源项的动量方程的平衡性。将方程中的通量项用()2F表示。对于式(51)中给出的数值通量，式(45)左侧的数值近似为 周翔宇 等 DOI:10.12677/aam.2023.128367 3735 应用数学进展 ()()()()()()()()()()()()()()11111122221111112222222222,ddd d111ddd d.222nnnnjnnnnjttxjjjjttttxjjjjttFUtt vFUtt vFU tvx tg htt v xg htt v xght vx t+=(22)由于()()()(),jhbthbtC+=+，源项近似变成()()1111122221111122222222d ddd1dd.2nnjjnnjtxjjjjtItxjjtIjjSt v x tgCtbvbvb vxgtbvbvb vx+=+注意这个等式()()()()()()()22221111.2222xxxxghtg hbtg hbt bgbgC bg b=+=+并由此我们得出结论，通量和源项近似在静水稳态时相互平衡，从而得到期望的良好平衡特性。2.8.Positivity-Preserving 方法方法 在52中，Zhang 和 Shu 提出了设计双曲守恒律的高阶最大原理保持方法的框架。从那时起，该方法得到了许多关注，并被广泛应用，包括53中的浅水波方程。结果表明，在合适的 CFL 条件下，能保持水高的非负性而不影响质量守恒，并能保持通解的高阶精度。在这里，我们将探讨这种方法在第 2 节中为浅水波方程提出的平衡良好的 ADER-DG 方法中的应用。如54所述，实现这一目标的关键组成部分是以下两项：该方法一阶版本的正性，以及与高阶方法耦合的简单保正限制器。在建立之后，我们只考虑简单的欧拉时间离散化，同样的结果可以推广到多步和TVD 高阶龙格库塔方法。考虑平衡良好的 ADER-DG 方法(45)与(51)中定义的数值通量，取测试函数1v=导致湿截面的单元平均更新如下()()11221111,nnjjjjhhFF+=(23)其中()()()()1111111122222221*,*,*,*,1d.2nntjjjjjjjtFhuthuthht+=+(24)并且有11x=。为方便起见，设1nt=。为了不丧失一般性，我们在本节中忽略下标。该方法的一阶形式为()()11221111.nnjjjjhhFF+=(25)其数值通量定义为()()()()111111122222221*,*,*,*,1,2njjjjjjjFthuthuthh+=+(26)其中()max ugh=+。关于这种一阶方法的保正性，我们有如下引理：引理引理 2 在 CFL 条件1下，考虑具有数值通量(57)的一阶格式(56)。如果1,nnjjhh是非负的，那么1njh+周翔宇 等 DOI:10.12677/aam.2023.128367 3736 应用数学进展 也是非负的。证明 格式(56)可写成 11Nrrw=()()()()()()()112211221122*,*,1*,*,11*,*,11111112222111.22jjnnnnnnjjjjjjnnjjjjnnnjjjjjhhhuuhuhuhhhhuhuh +=+因此，1njh+是1njh，njh和1njh+的线性组合，并且所有系数都是非负的，因为12*,0jh+，12*,njjhh+。因此，10njh+。接下来，我们讨论高阶格式。我们参考54了解细节，这里只介绍主要思想。我们在区间jI上引入 N点(23Nk)Legendre gaas-lobatto 正交规则，用1122121,NNjjjjjjjSxxxxxx+=，表示这些正交点，在区间1 2,1 2上对应的正交权值rw满足。由于求积对于 k 次多项式是精确的，我们有()()11221121d.jNNnnnrnrjjrjjrjjrNjjIrrhhxw hxw hxw hw hx+=+(27)按照55的方法，我们得到的结果是：命题命题 3 对于浅水波方程(1)的格式(45)，我们设()njhx为单元格jI中的 DG 多项式。如果1122,jjhh+和()()2,1nrjjhxrN=都是非负的，则在 CFL 条件下1njh+也是非负的：1.w (28)证明 我们把(58)代入(54)，重写(54)通过加减()()()()1111122221*,*,dnntjjjjtfhuthutt+：()111 12,NnnrjrjjNNrhw hxw hw h+=+(29)以及()()()()()()()()111111111122222222211*,*,*,*,11,d,nntjjjjjjjjjtrhhfhuthutfhuthuttw+=(30)()()()()()()()()111111111122222222211*,*,*,*,1,d.nntNjjjjjjjjjtNhhfhuthutfhuthuttw+=(31)注意到()()()()()()()11111111222222211111122222212112212*,*,*,*,111*,*,*,*,*,112d11d2nnnntjjjjjjjtjjjjjjtjnjjtjhhhuthuthhwhuthuthhthhuttwwht+=+=+()()()11111222111111122222*,*,1*,*,1111d2111dd.22nnnnnntnjjjtttnnjjjjjtthutt hwthutt hutt hwwwtt+(32)周翔宇 等 DOI:10.12677/aam.2023.128367 3737 应用数学进展 可以清楚地看到1h是1122*,jjhh+和12*,jh的线性组合，并且所有系数都是非负的，因为1122*,0jjhh+(这也证明了我们在11 w的合适 CFL 条件下分别定义 和为(49)和(50)的原因。因此，我们得到10h，同样有0Nh。由于1Nww=，这两个 CFL 条件相同，也就是(59)。因此，10njh+，因为它是1,Nh h和()()2,1nrjjhxrN=的凸组合。我们在 DG 多项式()()()()()T,nnnjjjUxhxhux=上采用了以下保正极限器56 57，()()(),min 1,njnnnnjjjjnjjhUxUxUUhm=+=(33)并有()()1,minmin.jnnrjjjjx SrNmhxhx=(34)很容易观察到()()01,nrjjhxrN=，该极限保持变量()njUx的局部守恒。我们计算修正的多项式()njUx，并在均衡方法中使用()njUx代替()njUx(45)。根据56 57中的证明，我们可以验证在 CFL 条件下，与该保正限制器耦合的平衡方法是高阶精度、保正和质量守恒的 1.w (35)为了提高效率，只有当对下一个时间步长的初步计算产生负湿截面时，我们才实现时间步长限制(66)。3.数值结果数值结果 为了证明 well-balanced 的 ADER-DG 格式的良好性能，我们给出了几个标准的浅水波方程的数值算例的结果。在所有的算例中，CFL 数和引力常数 g 分别取 0.16 和 9.812。3.1.平坦底部的黎曼问题平坦底部的黎曼问题 在本节中，我们考虑包含平坦底部的干面积(即()0b x)的黎曼问题。这些例子在58中使用过，这 Figure 1.The numerical and exact solutions of the first Riemann problem at different time with 200 cells in Section 3.1:left:the water height h;right:the discharge hu 图图 1.第 3.1 节中 200 个网格的第一黎曼问题在不同时间的数值解和精确解水高 h(左)和水体动量 hu(右)周翔宇 等 DOI:10.12677/aam.2023.128367 3738 应用数学进展 里选择它们是为了证明我们的方法的保正能力。该例的计算域设置为300,300，初始条件为()()10 if0,00 and,00otherwise.xhu xu x=左侧为静水高度设为 10，右侧为干区。该问题的解析解可在59中找到。我们用我们的平衡保正方法计算了这个问题，该方法具有简单的传输边界条件和 200 个均匀单元。4,8,12t=时刻的解如图1 所示。我们还在这些图中绘制了精确的解，以便进行比较。我们可以观察到，数值结果很好地反映了精确解。3.2.抛物线碗抛物线碗 Sampson 等人60推导出了具有抛物型底地形的一维浅水波方程的解析解。这为我们的数值方法提供了一个很好的测试用例。这个例子已在61中用于具有摩擦源项的浅水波方程。我们取抛物线底部：Figure 2.The water surface level in the parabolic bowl problem at different time in Section 3.2:top left:t=1000;top right:t=2000;bottom left:t=3000;bottom right:t=4000 图图 2.第 3.2 节中抛物面碗问题不同时刻的水面水平：左上：t=1000；右上：t=2000；左下：t=3000；右下：t=4000 周翔宇 等 DOI:10.12677/aam.2023.128367 3739 应用数学进展 ()()20,b xhx a=常数0h和 a 下文定义。计算域设置为5000,5000。对于不含摩擦源项的浅水波方程，解析水面由式给出()()()()22008,cos 2cos,442gBBBxh x tb xhttggag+=其中02,constantgha B=。干湿锋面的确切位置如图 2 所示：()200cos.2B axtagh=3.3.干燥河床的生成问题干燥河床的生成问题 我们首先考虑62中提出的平底()0b x=的数值实验。初始条件为()()0.3 if5,00.1,00,10.0.3 otherwise,xh xhu xx=在这个例子中，干燥的河床会在两个方向相反的稀有波的中间产生。干层的生成使得这个例子很难用数值计算。将 t=1 时刻得到的结果进行比较，正如图 3 可以观察到的那样，目前的格式在两个稀疏波之间留下了一个小的湿区。通过图中绘制出的精确解我们可以清晰的观察到，数值结果具有良好的保正性，并且很好地反映了精确解。Figure 3.Water depth h(left)and water discharge hu(right)at t=1 图图 3.t=1 时水深 h(左)和水体动量 hu(右)4.结论结论 本文设计了一种有效的保正高阶 ADER-DG 方法，用于求解具有非平坦底地形的一维浅水波方程。所提出的方法具有两个良好的特点：对静水稳态解具有良好的平衡性和在干湿锋附近保持正性。通过一种新颖的源项分裂和每个分裂源项的适当的良好平衡近似来实现良好平衡特性。采用63中的一个简单的保正限幅器，以确保所得方法保持截面湿面积的非负性。我们进行了相应的数值模拟实验，结果表明本文所提出的方法具有良好的平衡性，在稳态解附近的小扰动测试中有效，在干湿前沿附近保持正性以及周翔宇 等 DOI:10.12677/aam.2023.128367 3740 应用数学进展 高阶精度，并且在连续解和间断解上都表现良好。基金项目基金项目 本研究得到了山东省自然科学基金面上项目(No.ZR2021MA072,ZR2023MA012)的资助。参考文献参考文献 1 Vzquez-Cendn,M.E.(1999)Improved Treatment of Source Terms in Upwind Schemes for the Shallow Water Equ-ations in Channels with Irregular Geometry.Journal of Computational Physics,148,497-526.https:/doi.org/10.1006/jcph.1998.6127 2 Garcia-Navarro,P.and Vzquez-Cendn,M.E.(2000)On Numerical Treatment of the Source Terms in the Shallow Water Equations.Computers&Fluids,29,951-979.https:/doi.org/10.1016/S0045-7930(99)00038-9 3 Balbas,J.and Karni,S.(2009)A Central Scheme for Shallow Water Flows along Channels with Irregular Geometry.ESAIM:Mathematical Modelling and Numerical Analysis,43,333-351.https:/doi.org/10.1051/m2an:2008050 4 Hernndez-Dueas,G.and Karni,S.(2011)Shallow Water Flows in Channels.Journal of Scientific Computing,48,190-208.https:/doi.org/10.1007/s10915-010-9430-x 5 Audusse,E.,Bouchut,F.,Bristeau,M.-O.,Klein,R.and Perthame,B.(2004)A Fast and Stable Well-Balanced Scheme with Hydrostatic Reconstruction for Shallow Water Flows.SIAM Journal on Scientific Computing,25,2050-2065.https:/doi.org/10.1137/S1064827503431090 6 Xing,Y.and Shu,C.-W.(2005)High Order Finite Dif