[論文レビュー] 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

ひとこと要約

本稿は、移動導体問題におけるSU/PG法の横方向境界誤差の根本原因として、クーロンゲージ（∇·A = 0）を特定した。ゲージ制約を排除し、弱形式に∇·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法の横方向境界誤差の根本原因を特定すること。
• 高周速度下で生じる物質界面における非物理的電流および振動を解消すること。
• 非線形SOLDスキームに代わる、計算効率が高く線形な代替手法を開発し、境界誤差を解消しながら安定性を損なわないこと。
• 弱形式においてクーロンゲージ制約（∇·A = 0）を排除することで、誤差が解消されつつSU/PG法の効率が維持されることを示すこと。

提案手法

• 磁場中を移動する導体の2次元モデルを分析し、導体と空気の界面に注目する。
• SU/PG定式化において∇·A = 0 を強制することで、誤った境界誤差が生じることを特定する。
• アマルの法則の弱形式を再導出する際、制約としてではなく∇·A 項を保持する。
• 修正された定式化は、元のSU/PG法と同様に線形性と計算効率を維持する。
• さまざまなペクレ数における2次元および3次元の数値シミュレーションを実施し、提案手法の妥当性を検証する。
• 条件数とピーク誤差の指標を用いて、標準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法およびポールゼロキャンセレーション手法を上回り、精度と条件数の両面で優れた性能を示した。