[论文解读] 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

一句话总结

该论文提出了一种无参数的、加权残差的有限元公式，用于线性移动导体问题，通过引入辅助变量消除一阶导数项，从而在无需迎风法的情况下确保数值稳定性。该方法在1D、2D和3D模拟中均实现了稳定且精确的解，无需稳定参数或迭代优化，展示了在多种材料配置下的最优收敛性和鲁棒性。

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.

研究动机与目标

• 解决在高阶速度下模拟线性移动导体时，Galerkin有限元方法出现的数值振荡问题。
• 克服迎风法方案的局限性，如横向扩散和横向边界误差，这些误差通常需要迭代调整参数。
• 开发一种本质上稳定且无需稳定参数（τ）的公式，避免计算开销。
• 在电导率和磁导率各异的多材料系统中确保一致性和准确性。
• 验证该方法在1D、2D和3D配置中的稳定性与收敛性，包括TEAM 9a基准的复杂几何结构。

提出的方法

• 将磁矢势方程中的一阶导数项替换为新变量 bx = −dAy/dz，以消除对流项。
• 引入使用权函数 dN/dz 的第二个方程，以消除数值不稳定的 N(dAy/dz) 项。
• 采用 Ay 和 bx 作为主要变量，构建混合Galerkin弱形式，确保一致性和稳定性。
• 对二阶导数扩散项应用分部积分，推导弱形式。
• 使用Z变换极点-零点分析进行1D稳定性验证，确认无不稳定模态。
• 在1D、2D和3D FEM代码中实现该公式，包括3D情况下的轴向和径向离散化。

实验结果

研究问题

• RQ1 能否开发一种无参数的有限元公式，在无需迎风法的情况下稳定Galerkin方法用于线性移动导体问题？
• RQ2 通过引入辅助变量消除一阶导数项，是否能消除高速模拟中的数值振荡？
• RQ3 该公式的收敛速率和准确性在1D、2D和3D问题中，针对不同材料属性的表现如何？
• RQ4 该公式在涉及复杂几何结构（如TEAM 9a基准）的3D模拟中是否能保持稳定性和准确性？
• RQ5 该方法是否避免了迎风法中常见的横向边界误差，且无需迭代调整参数？

主要发现

• 所提出的公式在1D中无需任何稳定参数即可实现数值稳定性，经Z变换极点-零点分析验证。
• 在1D和2D测试案例中，该方法对∇×A以及磁场b或h均表现出最优的收敛阶数1。
• 在µr = 50且电导率变化的2D模拟中，该公式在多个材料界面间保持了准确性和稳定性。
• TEAM 9a问题的3D模拟结果与2D情况一致，且随着角向离散化程度提高，精度进一步提升。
• 该公式在3D中无需迭代优化或参数调整即可产生稳定且精确的结果，在计算效率上优于迎风法。
• 在内部点处，反应磁场 br 均能从 ∇×A 一致且准确地计算得出，避免了由辅助方程引入的误差。