基于格子玻尔兹曼高阶展开的有限差分格式：提升对流-扩散方程计算精度的新方法

《Computers & Mathematics with Applications》：Lattice-Boltzmann-inspired finite-difference schemes for the convection-diffusion equation

本研究针对标准格子玻尔兹曼方法（LBM）在求解对流-扩散方程时精度受限的问题，提出两种新策略。其一，通过六阶展开推导最优松弛时间以消除高阶误差；其二，基于此展开构建了三级和四级有限差分格式（TLFD/FLFD）。理论分析与高斯-丘陵基准测试表明，在最优松弛时间下，新方法在粗网格上可达四阶精度，其相对误差可比标准LBM低两个数量级，显著提升了数值求解的准确性与效率。

　　

• •
  基于对格子玻尔兹曼方程进行六阶展开，推导出了可消除特定高阶误差项的最优松弛时间。
• •
  受此展开启发，开发了两种全新的、仅用平衡分布函数表达的多级有限差分格式。
• •
  理论分析证实，在粗网格上，这两种格式在各自的最优松弛时间下均可实现四阶精度。
• •
  通过高斯-丘陵基准测试验证，数值计算的最优松弛时间与收敛速率与理论预测完美吻合。
• •
  结合狄利克雷和诺伊曼边界条件的高阶外推公式，四级有限差分格式实现了最低的相对误差，比标准格子玻尔兹曼模型低两个数量级。

<2. 标准格子玻尔兹曼方法中的高阶截断误差>

本研究聚焦于以下对流-扩散方程：

?_tφ + ?·(φu) = ?·(D?φ)，

其中φ代表一个标量参数，u是对流速度，D是一个各向异性扩散张量。为了在格子玻尔兹曼方法的框架内求解此方程，方程(1)可重写为：

?_tφ + ?·(φu) = ?·(α?·K)，

此处K = Dφ/α，而α是一个扩散参数。在不失一般性的前提下，我们采用D2Q9格子模型。后续分析旨在通过查普曼-恩斯科格展开识别标准LBM中的高阶截断误差，并推导出能够抵消主导误差项的最优松弛时间表达式，从而为构造高精度有限差分格式奠定基础。

<3. 受格子玻尔兹曼启发的有限差分格式>

<3.1. 三级有限差分格式>

受格子玻尔兹曼启发的有限差分格式的核心思想，是用有限差分公式来近似方程(18)等号右侧的主导阶项。我们假设方程(18)可以重写为如下形式：

g_i(x, t) ≈ K₀?_i^eq(x, t) + K₁?_i^eq(x ? e_iδ_t, t ? δ_t)

+K₂?_i^eq(x ? 2e_iδ_t, t ? 2δ_t) + K₃?_i^eq(x ? 3e_iδ_t, t ? 3δ_t)，

其中系数K_i可以通过系数匹配法确定。通过巧妙地选择系数，可以构建出仅依赖三个或四个时间层（即“三级”或“四级”）的有限差分格式。这些格式完全基于平衡分布函数，绕过了传统LBM中分布函数的演化步骤，在计算上更为直接。理论分析表明，在精心选择的最优松弛时间下，三级有限差分格式（TLFD）和四级有限差分格式（FLFD）均能以四阶精度恢复目标对流-扩散方程。

<3.2. 四级有限差分格式与边界处理>

四级有限差分格式（FLFD）在三级格式的基础上，通过引入更多的时间层信息，进一步提升了格式的精度和稳定性。对于复杂的物理模拟，边界条件的正确处理至关重要。为了充分发挥高阶格式的潜力，我们为FLFD格式开发了相应的高阶外推公式，分别用于处理狄利克雷边界条件（指定边界上的值）和诺伊曼边界条件（指定边界上的法向导数）。这些高阶边界处理技术确保了即使在计算域的边缘，数值解的精度也能与内部区域保持一致，从而显著降低了整体误差。

<4. 数值测试>

<4.1. 高斯-丘陵基准测试>

我们采用高斯-丘陵问题作为基准测试，以验证所提出的模型及预测的最优松弛时间。计算域[-1,1]×[-1,1]被离散为N×N的网格，所有边界均施加周期性边界条件。对于一个给定的扩散张量D，其解析解为：

φ^ana(x, t) = [φ₀/(2π√(‖σ_t‖))] exp(-(1/2)σ_t^-1:(x-ut)(x-ut))，

其中σ_t= σ₀²I + 2tD。我们设置φ₀= 2πσ₀²，使得初始分布由下式给出：

φ(x, 0) = exp(-|x|²/(2σ₀²))。

通过对比不同模型（标准BGK-LBM、TLFD、FLFD）在不同网格密度和松弛时间下的数值解与解析解，我们评估了它们的收敛速度和精度。结果表明，在理论推导出的最优松弛时间下，TLFD和FLFD格式确实达到了四阶收敛精度，且FLFD格式结合高阶边界处理后，其相对误差比标准LBM低了整整两个数量级，优势非常显著。

<4.2. 复杂边界条件下的误差评估>

除了周期边界，我们还测试了所提方法在具有狄利克雷和诺伊曼边界条件的更实际场景下的表现。通过将FLFD格式与专门设计的高阶外推边界格式相结合，即使在计算域边界附近，数值解的精度也得以有效保持。误差分析表明，这种组合策略能够显著抑制由边界处理引入的额外误差，从而在各种复杂边界条件下都能获得高精度的全局数值解，验证了该方法在实际应用中的鲁棒性。

<5. 总结与结论>

本文受格子玻尔兹曼方程高阶展开的启发，发展出适用于对流-扩散方程的多级有限差分格式。所提出的三级（TLFD）和四级（FLFD）格式能够以二阶精度恢复目标方程。为了进一步提升数值精度，我们通过高阶截断误差分析推导了BGK、TLFD和FLFD模型的最优松弛时间，并随后利用高斯-丘陵基准测试进行了验证。数值实验证实，在最优松弛时间下，TLFD和FLFD格式在粗网格上均能达到四阶精度。更重要的是，结合高阶边界外推公式的FLFD格式，其相对误差可比标准格子玻尔兹曼模型低两个数量级。这项工作为高效、高精度地求解对流-扩散方程提供了一类新的、有潜力的数值工具。