应用有限元法与引入迎风思想的有限体积法相结合处理运动电磁问题研究.pdfVIP

应用有限元法与引入迎风思想的有限体积法相结合处理运动电磁问题 应用有限元法与引入迎风思想的有限体积法相结合 处理运动电磁问题 甘艳阮江军张宇 (武汉大学电气工程学院，湖北省武汉市430072) 摘 要 本文提出了应用有限元法与引入迎风思想的有限体积法相结合处理运动电磁问题新方法，讨论了运动电磁闯 题控制方程的特性，详细阐明了一般Galerkin有限元法离散格式的误差估计，从而说明了一般Galerkin有限元 法处理对流占优型对流扩散方程时解的不稳定性。在充分考虑了有限元法与有限体积法各自的特性和优点后，提 出了用有限元法离散控制方程的扩散项而用引入迎风思想的有限体积法离散控制方程的对流项。应用这种方法计 算了一个简单的运动问题，结果表明本文提出的这种方法町以很好的处理运动问题。 关键词：对流扩散方程控制方程Galerkin有限元法有限体积法迎风系数 1引言 目前，工程应用一个比较突出的问题是场域内介质间的相对运动问题，常见于电机、电磁 物体发射及涡流闸等电磁装置。在分析这些问题时，要考虑冈运动引起的切割电势。解决这类问 题的难点主要有二个方面：如何解决数值解的振荡问题以及如何处理介质间的相对运动使得有限 元剖分网格发生变化的问题。目前，讨论介质间相对运动时的网格处理方面的文献较多，关于数 值解振荡的问题也已展开了广泛的研究，并取得了一定的成果。关于解决数值解的振荡问题已达 成的共识是：通过其它的方式将速度因子y融合到刚度矩阵中，避免出现矢量位函数4对空间坐 标的偏导数。文献【Il提出了消除数值解失真振荡的逆风有限元思想，这种处理方法的思想是通过 对权函数的修正、或高斯积分点的偏移使得有限元方程获得稳定的数值解。但是随着导体运动速 度的增加，系数矩阵中速度项所占比例逐渐增人，从而导致对角线元素逐渐减小，甚至会导致系 数矩阵失去“主元占优”的特性。因此这种方法在对对流项(速度项)进行处理的同时还要对扩 散项进行类似的分析、处理。我们知道，当控制方程中没有对流项时，用有限元法求解的稳定性 是可以保证的，导致数值解不稳定的原因是对流项，那么是否可以将对流项与扩散项分开处理， 用常规的有限元法处理扩散项，而用其他的方法来离散对流项，当然，离散对流项的方法应该保 证对对流项处理的稳定性。本文提出了将有限元法与引入迎风思想的有限体积法相结合处理运动 电磁问题的新方法，用常规有限元法处理控制方程中的扩散项而用引入迎风思想的有限体积法处 理对流项。这种方法在不增加计算量的前提下保证了计算精度，消除了数值解的振荡。 2控制方程的特性 描述运动介质间的电磁场的控制方程为： v×(土v×彳)：以+以 (1) ∥ 其中，以为源电流密度，以为涡流密度。 v×(土v×么)一仃(一罢+V×(v×彳))：以 (2) ∥ 研 在二维问题中，引入库仑规范后，方程(1)可变为： !v2A-o．(V．v)A一仃要：一以 (3) 弘 os 式(3)中A是磁场的矢量磁位；矿是介质相对于参考坐标系的运动速度：以是源电流密度； 仃是介质的电导率；∥是介质的磁导率。 式(3)中第一项表示由于介质的特性产生的场，一般称为扩散项，对于运动磁场问题扩散系数 为二。第二项表示由于介质的运动产生的场，一般称J,j对流项，对于运动磁场问题对流系数为导 ∥ 体的电导率与运动速度的乘积。当对流系数远远大于扩散系数时，这类控制方程属于对流占优型 的对流扩散方程。这类问题的控制方程中有一阶导数项，正是这个一阶导数项给运动场的计算增 ／JDT许多复杂因素。在数值计算中所遇到的主要困难就是由这个一阶导数项的离散方式不妥而引 起的121。对流扩散型偏微分方程的共同的特点是：对流的影响远大于扩散的影响，即对流占优性。 具有对流占优性的对流扩