[論文レビュー] A Stable Weighted Residual Finite Element Formulation for the Simulation of Linear Moving Conductor Problems

Sethupathy Subramanian, Sujata Bhowmick|arXiv (Cornell University)|May 25, 2022

Advanced Numerical Methods in Computational Mathematics参考文献 25被引用数 2

ひとこと要約

本稿では、一次微分項を補助変数を用いて除去することで、上流化を用いず数値安定性を確保する、パラメータフリーで重み付き残差有限要素法を線形移動導体問題に提案する。この手法は、安定かつ高精度な解を得ることができ、1次元、2次元、3次元のシミュレーションにおいて安定化パラメータや反復的微調整を必要とせず、複数の材料配置において最適収束性と頑健性を示す。

ABSTRACT

The finite element method is one of the widely employed numerical techniques in electrical engineering for the study of electric and magnetic fields. When applied to the moving conductor problems, the finite element method is known to have numerical oscillations in the solution. To resolve this, the upwinding techniques, which are developed for the transport equation are borrowed and directly employed for the magnetic induction equation. In this work, an alternative weighted residual formulation is explored for the simulation of the linear moving conductor problems. The formulation is parameter-free and the stability of the formulation is analytically studied for the 1D version of the moving conductor problem. Then the rate of convergence and the accuracy are illustrated with the help of several test cases in 1D as well as 2D. Subsequently, the stability of the formulation is demonstrated with a 3D moving conductor simulation.

研究の動機と目的

• 高周速度下での線形移動導体のシミュレーションにおいて、ガレルキン有限要素法で生じる数値的振動を解消すること。
• 上流化スキームの限界（横方向拡散や境界誤差）を克服し、反復的パラメータ調整を必要としないこと。
• 安定化パラメータ（τ）を含まない、本質的に安定な定式化を構築すること。
• 異なる導電率および透磁率を有する多相材料系において、一貫性と正確性を確保すること。
• 1次元、2次元、3次元の設定、特にTEAM 9aベンチマークのような複雑な幾何形状を含む場合に、手法の安定性と収束性を検証すること。

提案手法

• 磁気ベクトルポテンシャル方程式における一次微分項を、新たな変数 bx = −dAy/dz で置き換えることで、対流項を除去する。
• 重み関数 dN/dz を用いた第二の式を導入し、数値的に不安定な N(dAy/dz) 項を排除する。
• Ay と bx を主変数とする混合ガレルキン弱形式を定式化し、一貫性と安定性を確保する。
• 2階微分項の拡散項に対して部分積分を適用し、弱形式を導出する。
• 1次元の安定性検証にZ変換の極零点解析を用い、不安定モードが存在しないことを確認する。
• 1次元、2次元、3次元のFEMコードに定式化を実装し、3次元ケースでは軸対称および径方向分割を含む。

実験結果

リサーチクエスチョン

• RQ1 上流化を用いず、安定化パラメータを必要としない、パラメータフリーな有限要素法定式化を、線形移動導体問題に開発可能か？
• RQ2 補助変数を用いて一次微分項を除去することで、高周速度シミュレーションにおける数値的振動が解消されるか？
• RQ3 異なる材料特性を有する1次元、2次元、3次元問題において、提案手法の収束速度および正確性はどのように評価されるか？
• RQ4 TEAM 9aベンチマークのような複雑な幾何形状を含む3次元シミュレーションにおいて、定式化は安定性と正確性を維持できるか？
• RQ5 反復的パラメータ調整を必要とせず、上流化スキームに一般的に見られる横方向境界誤差を回避できるか？

主な発見

• Z変換の極零点解析により、1次元においても安定化パラメータを一切使用せずに数値的安定性が確認された。
• 1次元および2次元のテストケースにおいて、∇×A および磁場 b や h の両方で最適収束次数 1 を示した。
• µr = 50 で導電率が変化する2次元シミュレーションにおいて、複数の材料界面を介しても正確性と安定性が維持された。
• TEAM 9a問題の3次元シミュレーションでは、2次元ケースと一貫した結果が得られ、角度方向の離散化を増加させることで正確性が向上した。
• 反復的微調整やパラメータ調整を一切行わず、3次元でも安定かつ正確な結果が得られ、上流化スキームに比べ計算効率に優れた。
• 内部点において、∇×A から得られる反応磁場 br は一貫的かつ正確に計算され、補助方程式に起因する誤差を回避した。