[发明专利]一种数值模拟磁感应现象的积分方程法

说明书

本发明公开了一种采用积分方程法数值模拟磁感应现象的方法，该方法首先把偏微分方程形式的控制方程变为积分方程形式，积分方程中所有积分的积分区域均局限于可导体占据的有限区域，因此避免了对无限大区域的离散，从而减少了计算工作量和对内存的需求。由于积分方程形式的磁感应方程中不包含任何关于未知函数的导数，因此非常适用于求解在局部区域具有不连续或不光滑的磁场的磁感应现象。由于积分方程对于被积函数的要求低，利用积分方程形式的磁感应方程便于求解磁导率或电导率不连续的问题。本发明可用于任何形状的可导体的磁感应现象的数值模拟。特别是可用于研究快速增殖核反应堆的安全问题。

技术领域

本发明涉及一种采用积分方程法进行数值模拟磁感应现象的方法。

背景技术

在快速增殖反应堆和正在研究中的磁约束受控核聚变反应堆中，液态金属均用于冷却剂。在这些工况下，液态金属的流动通常具有大的磁雷诺数，且有可能出现磁流体自激放大发电现象，影响反应堆的安全。同时在反应堆中通常使用具有高磁导率的磁性材料，其有可能引起大的磁感应磁场的产生，影响反应堆的安全。本发明拟提出一种高效和应用范围广的数值模拟这些现象的积分方程法。

通常磁感应现象可由如下pre-Maxwell方程组

和如下欧姆定律

j(r,t)＝σ(r)[E(r,t)+F(r,t)] (4)

描述。其中，r是位置矢量，t是时间，E是电场，F是电动势，如果可导流体的流动是层流，则有

F(r,t)＝u(r)×B(r,t) (5)

其中，u是流体流动的速度，如果可导流体的流动是湍流，则有：

其中，包含α的项表示了α效应，包含β的项表示了β效应。设可导流体约束在具有有限体积D的容器中，容器外的无限大区域记为D′。由上述方程可得以磁场为未知函数的磁感应方程如下：

在有限体积D内,磁场满足:

在容器外无限大区域D′内,磁场还满足:

在所有区域上,磁场均满足:

在具有不同的电导率或磁导率的区域的交界面上磁场满足如下连续性条件:

r＝r→∞当时

B＝O(r^-3) (10)

其中,B是磁场，u是速度，μ是磁导率，σ是电导率，n是交界面的单位法线矢量。

数值模拟上述方程存在两大困难，一是如何处理可导流体外的无限大区域D′，二是往往会出现两个区域具有差别非常大的磁导率，磁导率在其交界面上的巨大跳跃严重影响现有许多方法的稳定性。目前对第一个问题，一种策略是采用“Vertical fieldcondition”近似的边界条件，从而避免对无限大区域的划分。另一种方法是围绕可导流体做一足够大的球，设在可导流体和球面之间充满低电导率的材料。以上两种处理方法和真实的边界条件相差甚远，必然导致大的计算误差。另外一类方法是对可导性流体的流动采用有限体积或有限元法，而对外面的无限大区域采用边界元，但这一类方法存在如何衔接两种不同的方法的问题，因此至今没有有关和现有实验比较的结果，有待实验验证。