《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)可以重寫為如下形式：
gi(x, t) ≈ K0?ieq(x, t) + K1?ieq(x ? eiδt, t ? δt)
+K2?ieq(x ? 2eiδt, t ? 2δt) + K3?ieq(x ? 3eiδt, t ? 3δt)，
其中系數Ki可以通過系數匹配法確定。通過巧妙地選擇系數，可以構建出僅依賴三個或四個時間層（即“三級”或“四級”）的有限差分格式。這些格式完全基于平衡分布函數，繞過了傳統LBM中分布函數的演化步驟，在計算上更為直接。理論分析表明，在精心選擇的最優松弛時間下，三級有限差分格式（TLFD）和四級有限差分格式（FLFD）均能以四階精度恢復目標對流-擴散方程。
<3.2. 四級有限差分格式與邊界處理>
四級有限差分格式（FLFD）在三級格式的基礎上，通過引入更多的時間層信息，進一步提升了格式的精度和穩定性。對于復雜的物理模擬，邊界條件的正確處理至關重要。為了充分發揮高階格式的潛力，我們為FLFD格式開發了相應的高階外推公式，分別用于處理狄利克雷邊界條件（指定邊界上的值）和諾伊曼邊界條件（指定邊界上的法向導數）。這些高階邊界處理技術確保了即使在計算域的邊緣，數值解的精度也能與內部區域保持一致，從而顯著降低了整體誤差。
<4. 數值測試>
<4.1. 高斯-丘陵基準測試>
我們采用高斯-丘陵問題作為基準測試，以驗證所提出的模型及預測的最優松弛時間。計算域[-1,1]×[-1,1]被離散為N×N的網格，所有邊界均施加周期性邊界條件。對于一個給定的擴散張量D，其解析解為：
φana(x, t) = [φ0/(2π√(‖σt‖))] exp(-(1/2)σt-1:(x-ut)(x-ut))，
其中σt= σ02I + 2tD。我們設置φ0= 2πσ02，使得初始分布由下式給出：
φ(x, 0) = exp(-|x|2/(2σ02))。
通過對比不同模型（標準BGK-LBM、TLFD、FLFD）在不同網格密度和松弛時間下的數值解與解析解，我們評估了它們的收斂速度和精度。結果表明，在理論推導出的最優松弛時間下，TLFD和FLFD格式確實達到了四階收斂精度，且FLFD格式結合高階邊界處理后，其相對誤差比標準LBM低了整整兩個數量級，優勢非常顯著。
<4.2. 復雜邊界條件下的誤差評估>
除了周期邊界，我們還測試了所提方法在具有狄利克雷和諾伊曼邊界條件的更實際場景下的表現。通過將FLFD格式與專門設計的高階外推邊界格式相結合，即使在計算域邊界附近，數值解的精度也得以有效保持。誤差分析表明，這種組合策略能夠顯著抑制由邊界處理引入的額外誤差，從而在各種復雜邊界條件下都能獲得高精度的全局數值解，驗證了該方法在實際應用中的魯棒性。
<5. 總結與結論>
本文受格子玻爾茲曼方程高階展開的啟發，發展出適用于對流-擴散方程的多級有限差分格式。所提出的三級（TLFD）和四級（FLFD）格式能夠以二階精度恢復目標方程。為了進一步提升數值精度，我們通過高階截斷誤差分析推導了BGK、TLFD和FLFD模型的最優松弛時間，并隨后利用高斯-丘陵基準測試進行了驗證。數值實驗證實，在最優松弛時間下，TLFD和FLFD格式在粗網格上均能達到四階精度。更重要的是，結合高階邊界外推公式的FLFD格式，其相對誤差可比標準格子玻爾茲曼模型低兩個數量級。這項工作為高效、高精度地求解對流-擴散方程提供了一類新的、有潛力的數值工具。
