A superconvergent hybridizable discontinuous Galerkin method for the convective Cahn--Hilliard equation

數值分析中有一個長期被默認的取捨：迎風穩定化能壓制對流引起的非物理振盪，卻讓分片常數（k=0）近似的 L² 收斂率從最優 O(h²) 退化為 O(h)。四川大學、Buffalo 大學等五位研究者在 2026 年 4 月的 arXiv 新作中，針對對流 Cahn-Hilliard 方程提出超收斂 HDG（可雜化間斷 Galerkin） 格式，以放棄迎風通量、顯式離散對流項為核心策略，同時達成無條件穩定、k≥0 全階最優 L² 收斂以及對稱全局系統三項目標。

從二元合金相分離到多相流模擬

Cahn-Hilliard 方程起源於 1958 年描述二元合金相分離的現象學模型，用序參數 u（局部成分比例）和化學勢 ϕ 兩個場量共同決定界面的形成、粗化與遷移動力學。這個四階偏微分方程今天廣泛用於多相流、聚合物混合與表面潤濕的數值模擬，是計算材料科學與流體力學的核心工具之一。

加入背景流體速度場後，方程納入散度項 ∇·(βu) 形成對流 Cahn-Hilliard 方程，並可進一步與 Navier-Stokes 方程耦合構成 Cahn-Hilliard-Navier-Stokes（CHNS）系統。高 Péclet 數（Pe 大、對流主導）情形下，連續 Galerkin 方法在尖銳梯度附近容易產生非物理振盪；標準 DG 方法雖具局部守恆性與靈活通量設計，全局自由度數量卻遠多於連續方法。這一現實推動了 HDG 方法在 Cahn-Hilliard 類問題中的應用探索。

HDG 方法的架構與迎風通量的 k=0 收斂困境

HDG 方法由 Cockburn、Gopalakrishnan 和 Lazarov 於 2009 年提出，核心思想是在網格面（face）引入 trace 未知量 û 作為全局耦合自由度，而單元內部的通量和標量變量通過局部消去（local static condensation）恢復，將全局系統規模縮減至僅剩面上的未知數。這一設計讓 HDG 同時繼承 DG 的局部守恆性與超收斂後處理能力，並將求解規模控制在接近連續 Galerkin 的水準。

然而，HDG 類方法處理對流項時慣用的迎風穩定化 τ_c(u_h − û_h) 存在一個理論缺陷：它使 k=0 分片常數近似的 L² 收斂率從理論最優 O(h²) 退化為 O(h)。這一問題已在對流-擴散方程、最優控制問題和 Navier-Stokes 方程的 HDG 格式中反覆出現，卻尚無系統性解法。k=0 近似在含尖銳相界面的非線性問題中尤為重要，低次多項式在界面附近能提供更穩健的行為，收斂退化直接影響方法的實用性。

三個設計選擇：顯式對流、凸-凹分裂、對稱 MINRES 求解

論文核心突破在於對流項採用完全顯式離散，且不引入任何迎風穩定化。每個時間步計算中，對流算子 ℬ 直接使用前一步的 (u_h^{n-1}, û_h^{n-1})，從根源上繞開了穩定化引起收斂退化的問題。理論分析表明，這種顯式處理在無任何網格比例約束（不需要 Δt ∝ h² 之類的條件）的前提下，仍能對所有 k≥0 保持最優 L² 收斂率。

時間方向的穩定性由凸-凹分裂（convex-concave splitting）保證：自由能 f(u) = u³ − u 分解為凸部分 u³（隱式離散，保穩定）和凹部分 −u（顯式離散，保精度），使能量在每一步都單調下降，實現無條件能量穩定性。對比之下，Kay et al. 2007 年的 DG 格式對流項隱式離散，僅有條件穩定，要求時間步長 Δt ≤ C·ε²；當界面厚度參數 ε 極小時（如 ε = 0.01），這一限制可讓時間步長縮小至計算不可行的量級。

對流項顯式化還帶來第三項紅利：全局求解系統保持對稱結構。消去單元內部自由度後，全局耦合矩陣仍為對稱形式，可直接用最小殘差法（MINRES）等高效對稱迭代求解器求解，在大規模三維計算中具有顯著效率優勢。

新型耦合 HDG 橢圓投影：k≥0 最優 L² 估計的技術基石

誤差分析方面，論文最重要的技術創新是設計了全新的 HDG 橢圓投影算子。在作者先前無對流的 HDG 格式中，u 和 ϕ 的橢圓投影可分別從兩個解耦 Laplace 方程獲得；引入對流後，ϕ 的投影方程出現依賴 u 的對流-擴散項，必須將兩者聯立為一個耦合系統（一個對流-擴散型加一個 Laplace 型）同時求解。若沿用舊的解耦投影算子，L² 誤差估計將不可避免地退化。

新橢圓投影在 L² 範數下對標量 u、ϕ 均達到 h^{k+2} 逼近階，對通量 q、p 達到 h^{k+1}，對時間導數 ∂_t u 保持相同最優估計。論文結合第三 Strang 引理（量化一致性誤差貢獻）和 Aubin-Nitsche 對偶技術（從 H¹ 提升至 L² 估計），在凸多邊形域的正則性假設下，對所有 k≥0 建立完整最優 L² 收斂理論，誤差常數僅依賴 Pe、ε 和 T，與網格尺寸 h 和時間步長 Δt 無關。

數值驗證：k=0 收斂率對比與 k+2 超收斂確認

論文的數值實驗系統地驗證了三個核心理論預測。迎風通量有無的收斂率對比最為直觀：本文格式（無迎風通量）k=0 的 L² 誤差呈 O(h²) 衰減；切換為帶迎風穩定化的格式，k=0 退化為 O(h)，而 k≥1 時兩種格式均達最優——與理論完全吻合，並由 Table 2 的數值結果清晰呈現。能量曲線實驗顯示，無論時間步長多大，能量在整個演化過程中都單調遞減，驗證無條件穩定性。此外，k+2 超收斂現象也得到數值確認：HDG 格式中 trace 未知量 û 的 L² 誤差比單元內部未知量高一階，達到 O(h^{k+2})，體現了 HDG 方法固有的超收斂機制。

論文最後指出，本工作建立的方法論工具——顯式對流 HDG 離散框架、新型耦合橢圓投影算子與完整的超收斂分析技術——將直接服務於後續構建更複雜 CHNS 方程組的超收斂 HDG 格式，對相場流體模型的高效數值求解具有長期方法論意義。

捨棄迎風穩定化換來三個同時成立的性質：k=0 最優 O(h²) 收斂、無條件能量穩定、對稱 MINRES 求解——對流 Cahn-Hilliard 方程 HDG 離散的一次精確設計示範。