求解具有吸引型内部转向点和间断源项的一维奇异摄动抛物型系统的分裂方法

《Journal of Computational and Applied Mathematics》：Splitting Methods for Solving one Dimensional Parabolic Singularly Perturbed Systems with an Attractive Interior Turning Point and a Discontinuous Source Term

　　

本文研究并求解了含有两个对流-扩散型方程的一维奇异摄动系统（singularly perturbed system, SPS），其中对流项导致吸引型内部转向点（attractive interior turning point）；此外，研究人员假设源项在该点存在跳跃间断；在此条件下，系统的解通常包含难以解析的内部边界层（internal boundary layer）。研究人员首先关注所有方程中扩散参数（diffusion parameter）相同的简单情形。所提出的数值方法将构造于特殊非均匀Shishkin网格（Shishkin mesh）上的迎风格式（upwind scheme）与标准一阶分裂方法（splitting method）相结合。研究人员证明，该组合得到的数值算法关于扩散参数一致收敛（uniformly convergent），时间方向达到一阶精度，空间方向达到几乎一阶精度。该方案的一个显著优点是易于推广至含更多分量的系统；此外，与经典方法相比，该算法显著降低了计算成本。针对不同测试问题的数值结果验证了该方法的实际一致收敛性和效率，与理论结果吻合。同时，研究人员仅通过数值实验表明，该技术可推广至系统中各方程扩散参数量级不同的情形。

本文解读基于发表在《Journal of Computational and Applied Mathematics》上的论文《Splitting Methods for Solving one Dimensional Parabolic Singularly Perturbed Systems with an Attractive Interior Turning Point and a Discontinuous Source Term》。研究背景涉及一类奇异摄动系统（singularly perturbed system, SPS），其由两个对流-扩散型方程耦合而成，在空间域内部存在一个吸引型转向点（attractive interior turning point），且源项在该点处具有跳跃间断。此类问题在磁流体动力学、湍流边界层分离、非均匀介质中光传播等物理现象中具有重要应用，但现有文献对含转向点及间断源项的耦合系统研究不足。目前存在的主要问题包括：解在转向点附近出现内部边界层（internal boundary layer），难以用传统数值方法高效解析；系统耦合导致经典隐式方法需要求解大型带状线性系统，计算成本高昂。为此，研究人员开展了本研究，旨在设计一种均匀收敛、计算高效的数值格式。

研究人员提出了一种结合空间离散与时间离散的数值算法。主要关键技术方法包括：第一，采用迎风格式（upwind scheme）在特殊的非均匀Shishkin网格（Shishkin mesh）上进行空间离散，该网格在转向点附近局部加密以捕捉强内部层；第二，采用分数步隐式欧拉方法（fractional implicit Euler method）结合分量分裂（splitting by components）进行时间离散，在每个时间层仅需求解三对角线性系统，从而大幅降低计算复杂度。该方法无需额外试剂或样本队列，仅基于数学推导与数值实验验证。

研究结果按论文各小节呈现如下：

**2. 连续问题精确解的渐近行为**：研究人员通过将解分解为常规分量（regular component）和内部层分量（interior layer component），并利用扩展区域技术，推导出各分量导数的精确估计。结果表明，当扩散参数ε充分小时，解在x=0处出现陡峭的内部边界层，且层宽度由ε^1/2控制。这一分析为后续网格设计提供了理论依据。

**3. 基于Shishkin网格的空间离散**：在非均匀Shishkin网格上应用迎风格式，研究人员构造了半离散问题。通过逆正性引理和障碍函数技术，证明了空间离散误差满足一致界‖[u(x,t)]_N?U_N(t)‖_ΩN≤C N^-1ln N，即空间离散具有几乎一阶的一致收敛性。

**4. 全离散格式与一致收敛性**：采用分数步隐式欧拉方法结合分量分裂对时间进行离散，得到全离散格式。研究人员证明了该格式的无条件稳定性（收缩性）和一致一阶时间收敛性，并最终通过三角不等式导出全局误差估计：max_0≤m≤M‖U_N^m?[u(x,t_m)]_Ω￣N‖_Ω￣N≤C(N^-1ln N+M^-1)，表明格式在空间和时间上均为ε-均匀收敛。

**5. 数值结果**：针对含两个和三个方程的系统进行数值实验，利用双网格原理估计误差。结果表明，对于不同ε取值，最大误差随N和M增加而减小，数值收敛阶接近一阶（空间接近一阶，时间一阶）。计算CPU时间对比显示，分裂算法相比经典隐式方法加速约1.2至3.3倍，且随方程数增加优势更明显。此外，对扩散参数量级不同的情况，研究人员通过构造更复杂的Shishkin网格（包含多个过渡参数）进行了数值测试，同样观察到几乎一阶的均匀收敛。

**6. 结论**：研究人员设计并分析了求解一类具有吸引型内部转向点和间断源项的一维奇异摄动抛物型系统的均匀收敛数值方法。该方法将基于Shishkin网格的迎风格式与分数步隐式欧拉分裂方法相结合，在时间上达到一阶精度，空间上几乎一阶精度，且关于扩散参数一致收敛。算法效率优于经典隐式方法，尤其当系统方程数增加时优势更为显著。数值算例验证了方法的均匀收敛性和低计算成本。此外，研究还展示了该离散技术可推广至多方程和多扩散参数的系统。