[论文解读] On Overcoming the Transverse Boundary Error of the SU/PG Scheme for Moving Conductor Problems

Sethupathy Subramanian, Udaya Kumar+1 more|arXiv (Cornell University)|Sep 17, 2021

Advanced Numerical Methods in Computational Mathematics参考文献 34被引用 3

一句话总结

本文識別出庫侖規範（∇·A = 0）是移動導體問題中SU/PG格式橫向邊界誤差的根本原因。透過移除規範約束並在弱形式中保留∇·A項，作者推導出一種改良的線性格式，可消除非物理的邊界振盪，同時保持SU/PG方法的計算效率。數值結果顯示，峰值誤差與條件數顯著降低，特別是在高佩克萊特數時效果更為顯著。

ABSTRACT

Conductor moving in magnetic field is quite common in electrical equipment. The numerical simulation of such problem is vital in their design and analysis of electrical equipment. The Galerkin finite element method (GFEM) is a commonly employed simulation tool, nonetheless, due to its inherent numerical instability at higher velocities, the GFEM requires upwinding techniques to handle moving conductor problems. The Streamline Upwinding/Petrov-Galerkin (SU/PG) scheme is a widely acknowledged upwinding technique, despite its error-peaking at the transverse boundary. This error at the transverse-boundary, is found to be leading to non-physical solutions. Several remedies have been suggested in the allied fluid dynamics literature, which employs non-linear, iterative techniques. The present work attempts to address this issue, by retaining the computational efficiency of the GFEM. By suitable analysis, it is shown that the source of the problem can be attributed to the Coulomb's gauge. Therefore, to solve the problem, the Coulomb's gauge is taken out from the formulation and the associated weak form is derived. The effectiveness of this technique is demonstrated with pertinent numerical results.

研究动机与目标

• 識別移動導體問題中SU/PG格式橫向邊界誤差的根本原因。
• 解決SU/PG方法在高運動速度下於介質介面產生的非物理電流與振盪現象。
• 發展一種計算效率高、線性的替代方案，以取代非線性SOLD格式，從而消除邊界誤差而不損失穩定性。
• 證明在弱形式中移除庫侖規範約束（∇·A = 0）可解決誤差問題，同時維持SU/PG方法的效率。

提出的方法

• 作者分析一個二維移動導體在磁場中的模型，專注於導體與空氣之間的介面。
• 他們識別出在SU/PG格式中強制執行∇·A = 0會引入偽造的邊界誤差。
• 透過保留下列∇·A項而非作為約束條件，重新推導安培定律的弱形式。
• 改良的格式保持了原始SU/PG格式的線性結構與計算效率。
• 在二維與三維中進行數值模擬，以驗證所提出方法在不同佩克萊特數下的有效性。
• 使用條件數與峰值誤差指標，將該方法與標準SU/PG、邊元法及極零點抵消方案進行比較。

实验结果

研究问题

• RQ1 為何SU/PG格式在移動導體問題的橫向介面處產生非物理的邊界誤差？
• RQ2 SU/PG格式中橫向邊界誤差的潛在根本原因為何？
• RQ3 是否可避免使用非線性、迭代的SOLD格式，而完全消除邊界誤差？
• RQ4 在弱形式中移除庫侖規範約束（∇·A = 0）是否能帶來穩定、準確且高效的解？
• RQ5 所提出方法在準確性與條件數方面，與SU/PG及邊元法等既有方案相比表現如何？

主要发现

• SU/PG格式中的橫向邊界誤差直接源於庫侖規範（∇·A = 0）的強制執行。
• 在弱形式中保留∇·A項可消除導體-空氣介面處的非物理峰值誤差。
• 在Pe = 1600時，所提出方法將峰值誤差由SU/PG的6.791%降至0.271%，條件數亦顯著降低，顯示穩定性提升。
• 該方法維持線性系統結構與計算效率，與非線性SOLD格式形成對比。
• 在Pe = 1000時，峰值誤差由SU/PG的30.021%降至0.271%，顯示一致的性能提升。
• 所提出方法在準確性與條件數方面均優於SU/PG與極零點抵消方案，特別是在高佩克萊特數時表現更佳。