多步IMEX局部间断Galerkin方法求解一维非线性四阶演化问题的稳定性与误差分析

《Mathematics》：Stability and Error Analysis of Multistep IMEX Local Discontinuous Galerkin Methods for One-Dimensional Nonlinear Fourth-Order Evolution Problems

　　

研究人员针对一类一维非线性四阶偏微分方程（PDEs），以Kuramoto–Sivashinsky（KS）方程为典型模型，研究了三种多步隐式–显式（IMEX）时间离散格式与局部间断Galerkin（LDG）空间离散的耦合方法。非线性对流项采用显式处理，而线性二阶、三阶和四阶导数项采用隐式推进。在拟一致网格假设以及精确解和数值解的均匀有界性假设下，能量估计给出了三种格式在时间步长限制??≤??₀下的??²稳定性，其中??₀与空间网格尺寸无关。对于所有三种格式，通过耦合LDG椭圆投影和能量估计，在多层次起始值具有相应精度的条件下，获得了最优空间收敛率??(?^??+1)以及时间收敛率??(??)、??(??²)和??(??³)。误差分析采用耦合LDG椭圆投影使离散辅助误差方程保持齐次性；残余的低阶投影项（包括相关项）被显式估计。独立的周期制造解测试验证了时间收敛率分别接近一阶、二阶和三阶。在特定条件下使用三阶格式进行的独立周期空间加密测试恢复了预期的空间收敛行为。对于长时间周期KS测试，在时间步长和网格加密两种条件下比较了动能历史和空间谱。后者被解释为有限观测窗口内的稳健性检验，而非统计平稳混沌状态或逐点长时间收敛的证据。

研究背景方面，一维非线性四阶时变偏微分方程在 interfacial flows、flame fronts 和 thin-film dynamics 等领域具有广泛的应用背景，其中 Kuramoto–Sivashinsky（KS）方程作为典型的非线性耗散–色散模型，以其丰富的斑图形成和确定性混沌行为著称。该类方程中四阶项提供高频耗散，二阶项在系数为负时具有失稳作用，三阶项为色散项，非线性对流项则起到调节作用。然而，对此类高阶非线性方程进行高效的数值求解面临严峻挑战：显式时间推进方法受限于高阶算子的苛刻时间步长约束，而全隐式处理非线性对流项则导致非椭圆型代数系统，迭代求解收敛性差。局部间断Galerkin（LDG）方法虽具有高阶精度、灵活自适应性和局部可解性等优势，但其与高效时间离散格式的耦合在全离散层面的稳定性与误差分析尚不完善，特别是针对同时包含二阶、三阶和四阶线性算子的非线性四阶模型的多步IMEX格式缺乏统一的 fully discrete 分析框架。已有工作分别涉及LDG空间离散或IMEX时间推进，但将多步IMEX与LDG耦合用于该非线性四阶模型并建立网格无关稳定性与全离散误差估计的研究尚属空白。在此背景下，研究人员开展了三种多步IMEX–LDG格式的全离散分析研究。

研究人员的主要工作是：对一维非线性四阶问题（以KS方程为典型模型）建立三种多步IMEX–LDG全离散格式，其中非线性对流项显式处理、线性二阶至四阶导数项隐式处理；在拟一致网格和均匀有界性假设下证明三种格式的能量稳定性，其时间步长上界与空间网格尺寸无关；通过耦合LDG椭圆投影建立最优全离散误差估计，时间精度分别达到一阶、二阶和三阶，空间精度达到??(?^??+1)。研究所获结论表明，三种格式在网格无关时间步长条件下均具有??²稳定性，且误差估计在相应起始值精度条件下达到最优阶。这一工作的重要意义在于首次为同时含二阶、三阶和四阶线性算子的非线性四阶演化问题提供了三种多步IMEX–LDG耦合的统一全离散稳定性与误差分析框架，避免了显式格式的网格相关CFL限制，并为高精度时间推进与LDG空间离散的结合作出了严格的理论保证。

在关键技术方法方面，研究人员采用的主要方法包括：局部间断Galerkin（LDG）空间离散，通过引入辅助变量将四阶导数分解为一阶导数链并采用交替数值迹以保证能量稳定性；三种多步隐式–显式（IMEX）时间离散格式，分别基于前向Euler、Adams–Bashforth二阶（AB2）和Adams–Bashforth三阶（AB3）外推处理非线性对流项并结合相应隐式权重处理线性导数项；耦合LDG椭圆投影技术，用于保持离散辅助误差方程的齐次性；以及能量估计与离散Gronwall不等式用于稳定性与误差分析。数值实验采用周期制造解测试、行波解测试和长时间混沌KS测试进行验证，其中混沌测试的起始历史值由高精度刚性BDF求解器生成以保证多层次起始精度。

研究结果方面，首先，在稳定性分析中，通过能量估计与离散Gronwall不等式，研究人员证明了三种多步IMEX–LDG格式在时间步长??≤??₀（??₀与网格尺寸h无关）条件下的??²稳定性，其中非线性对流项通过单调数值通量的单元熵不等式控制，低阶隐式项通过辅助变量能量恒等式控制。其次，在误差估计中，通过构造耦合LDG椭圆投影使离散辅助误差方程保持齐次性，并对残余低阶投影项进行显式估计，获得了三种格式的最优全离散误差界，即空间??(?^??+1)分别配合时间??(??)、??(??²)和??(??³)，多步格式的结果以相应精度的起始值为条件。第三，在独立时间收敛性验证中，采用周期制造解测试固定空间离散仅加密时间步长，验证了三种格式的时间收敛率分别接近一阶、二阶和三阶，其中三阶格式在 finest refinement 上达到约三阶的观测收敛率。第四，在独立空间收敛性验证中，采用三阶格式在周期制造解上固定小时间步长进行空间加密，观测到空间收敛率基本达到三阶，与??=2时的??(?³)理论预期一致。第五，在行波解测试中，对KS方程及含三阶导数项的KS方程在非均匀网格上进行了空间–时间联合加密测试，??²和??^∞误差随加密一致衰减，激波剖面传播保持尖锐且无可见伪振荡，但研究人员明确指出这些表格应解释为联合空间–时间收敛数据而非纯时间阶验证。第六，在长时间周期KS混沌测试中，通过动能历史和直接DG谱诊断评估了时间步长和网格加密的稳健性：时间步长从??=0.01减至??=0.005时，时间平均谱的相对??²差异仅为0.41%；空间加密从??=128至??=256时，谱的相对??²差异为1.06%，耗散尾部衰减至约10^-5，这些结果被解释为有限观测窗口内的分辨率检验而非统计平稳或逐点长时间收敛的证据。

讨论部分，研究人员总结了本研究的主要贡献与局限性。研究结论部分翻译如下：研究人员研究了非线性四阶时变问题的三种多步IMEX–LDG离散格式，以Kuramoto–Sivashinsky方程为典型模型。非线性对流项显式处理，而线性二阶、三阶和四阶导数项隐式处理。在拟一致网格假设和均匀有界范围假设下，稳定性分析给出了时间步长条件??≤??₀下的??²界，其上界与网格尺寸无关，因此估计不要求网格相关的CFL限制。利用耦合LDG椭圆投影，研究人员获得了分别对应一阶、二阶和三阶格式的误差界，多步结果以相应精度的起始值为条件。数值实验现在区分了时间和空间验证。周期制造解时间测试恢复了分别约为一阶、二阶和三阶的 finest-refinement 收敛率，而三阶格式在特定条件下的空间测试给出了基本等于三的观测收敛率。行波表格仍为有用的联合空间–时间数据，但其三阶列不用于验证三阶误差定理。对于周期KS诊断，以某基线配置为基准。将时间步长从某值减半至另一值时，时间平均直接DG谱在观测窗口内的相对??²范数变化仅为0.41%。在另一波数处，空间加密使相对谱差异在相同频带上为1.06%，两个谱在特定高波数处均衰减至约10^-5。这些结果被解释为有限窗口的时间和空间分辨率检验，而非逐点长时间收敛或统计平稳不变测度的证据。本分析有意局限于一维情形。在二维张量积网格上，LDG导数链和交替迹原则上可以逐方向构造，一维投影估计的张量积类似物也可获得。主要的额外困难在于多维离散能量中方向界面项与非线性扰动的同步控制。在非结构三角网格上，推广更为显著：失去简单的一维耦合椭圆投影/交替迹结构，面方向必须一致地纳入交替迹，投影/误差估计必须在单元面上重建。对于二维KS或Cahn–Hilliard等模型，这些问题因多维非线性算子和（后者情形下的）化学势公式而更加复杂。因此，严格的二维推广留待未来工作，而非从现有一维理论中推断。