CN107391434B - 一种数值模拟磁感应现象的积分方程法 - Google Patents
一种数值模拟磁感应现象的积分方程法 Download PDFInfo
- Publication number
- CN107391434B CN107391434B CN201710590450.4A CN201710590450A CN107391434B CN 107391434 B CN107391434 B CN 107391434B CN 201710590450 A CN201710590450 A CN 201710590450A CN 107391434 B CN107391434 B CN 107391434B
- Authority
- CN
- China
- Prior art keywords
- equation
- integral
- magnetic induction
- magnetic
- magnetic field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种采用积分方程法数值模拟磁感应现象的方法,该方法首先把偏微分方程形式的控制方程变为积分方程形式,积分方程中所有积分的积分区域均局限于可导体占据的有限区域,因此避免了对无限大区域的离散,从而减少了计算工作量和对内存的需求。由于积分方程形式的磁感应方程中不包含任何关于未知函数的导数,因此非常适用于求解在局部区域具有不连续或不光滑的磁场的磁感应现象。由于积分方程对于被积函数的要求低,利用积分方程形式的磁感应方程便于求解磁导率或电导率不连续的问题。本发明可用于任何形状的可导体的磁感应现象的数值模拟。特别是可用于研究快速增殖核反应堆的安全问题。
Description
技术领域
本发明涉及一种采用积分方程法进行数值模拟磁感应现象的方法。
背景技术
在快速增殖反应堆和正在研究中的磁约束受控核聚变反应堆中,液态金属均用于冷却剂。在这些工况下,液态金属的流动通常具有大的磁雷诺数,且有可能出现磁流体自激放大发电现象,影响反应堆的安全。同时在反应堆中通常使用具有高磁导率的磁性材料,其有可能引起大的磁感应磁场的产生,影响反应堆的安全。本发明拟提出一种高效和应用范围广的数值模拟这些现象的积分方程法。
通常磁感应现象可由如下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”近似的边界条件,从而避免对无限大区域的划分。另一种方法是围绕可导流体做一足够大的球,设在可导流体和球面之间充满低电导率的材料。以上两种处理方法和真实的边界条件相差甚远,必然导致大的计算误差。另外一类方法是对可导性流体的流动采用有限体积或有限元法,而对外面的无限大区域采用边界元,但这一类方法存在如何衔接两种不同的方法的问题,因此至今没有有关和现有实验比较的结果,有待实验验证。
发明内容
本发明要解决的技术问题是克服现有技术的缺陷,提供一种采用积分方程法进行数值模拟磁感应现象的方法。
为了解决上述技术问题,本发明提供了如下的技术方案:
本发明提供一种采用积分方程法进行数值模拟磁感应现象的方法,其包括以下步骤:
考虑可导区域D,其由m个子区域Di(i=1,2,3,…,m)构成,子区域的表面记为Si。设Di子区域的导电率和导磁率分别是σi和μi。Sij表示Di和Dj的交界面。nij表示垂直于Sij表面的单位法矢量,且从Di指向Dj。设流场是稳态的,则磁场对时间的依赖关系可表示为:~eλt。应用积分理论和边界元法,磁感应方程可表示为如下积分方程组:
把上述积分方程中的积分利用复合梯形公式离散,方程(11)对应的离散方程写成分量形式为:
Br=LrkBk+λPrjAj+Nrlφl (14)
方程(12)对应的离散方程为:
0.5φl+Elmφm=HlkBk+λDljAj (15)或者为
Glmφm=HlkBk+λDljAj (16)
其中,Glm=0.5δlm+Elm。方程(13)对应的离散方程为:
Aj=QjkBk (17)
由方程(14)-(17)可得以磁场的分量在节点的值为未知数的代数方程组为:
Br=LrkBk+λPrjQjkBk+Nrm(G-1)mlHlkBk+λNrl(G-1)lmDmjQjkBk (18)
求解代数方程组(18)可得磁感应产生的磁场。
本发明所达到的有益效果是:
积分方程形式的磁感应方程中所有积分的积分区域均局限于可导体占据的有限区域,因此避免了对无限大区域的离散,从而减少了计算工作量和对内存的需求。
由于积分方程形式的磁感应方程中不包含任何关于未知函数的导数,因此非常适用于求解在局部区域具有不连续或不光滑的磁场的磁感应现象。
由于积分方程对于被积函数的要求很低,利用积分方程形式的磁感应方程便于求解磁导率或电导率不连续的问题。
本发明可用于任何形状的可导体的磁感应现象的数值模拟。特别是可用于研究快速增殖核反应堆的安全问题。
Giesecke等提出的数值模拟快速增殖反应堆的液态钠在磁性材料的影响下的流动可能引起的磁场自激放大发电现象时,只考虑了局部磁雷诺数小于等于20,相对磁导率小于20的情况。本软件将采用的积分方程法,可用于磁雷诺数和相对磁导率均大于20的情况。
具体实施方式
以下对本发明的优选实施例进行说明,应当理解,此处所描述的优选实施例仅用于说明和解释本发明,并不用于限定本发明。
一种采用积分方程法进行数值模拟磁感应现象的方法,其包括以下步骤:
考虑可导区域D,其由m个子区域Di(i=1,2,3,…,m)构成,子区域的表面记为Si。设Di子区域的导电率和导磁率分别是σi和μi。Sij表示Di和Dj的交界面。nij表示垂直于Sij表面的单位法矢量,且从Di指向Dj。设流场是稳态的,则磁场对时间的依赖关系可表示为:~eλt。应用积分方程理论和边界元法,磁感应方程可表示为如下积分方程组:
把上述积分方程中的积分利用复合梯形法离散,方程(19)对应的离散方程写成分量形式为:
Br=LrkBk+λPrjAj+Nrlφl (22)
方程(20)对应的离散方程为:
0.5φl+Elmφm=HlkBk+λDljAj (23)
或者为
Glmφm=HlkBk+λDljAj (24)
其中,Glm=0.5δlm+Elm。方程(21)对应的离散方程为:
Aj=QjkBk (25)
由方程(22)-(25)可得以磁场的分量在节点的值为未知数的代数方程组为:
Br=LrkBk+λPrjQjkBk+Nrm(G-1)mlHlkBk+λNrl(G-1)lmDmjQjkBk (26)
求解代数方程组(26)可得磁感应产生的磁场。
积分方程形式的磁感应方程中所有积分的积分区域均局限于可导体占据的有限区域,因此避免了对无限大区域的离散,从而减少了计算工作量和对内存的需求。
由于积分方程形式的磁感应方程中不包含任何关于未知函数的导数,因此非常适用于求解在局部区域不连续或不光滑的磁场。
由于积分方程对于被积函数的要求很低,利用积分方程形式的磁感应方程便于求解磁导率或电导率不连续的问题的求解。
本发明可用于任何形状的可导体的磁感应现象的数值模拟。特别是可用于研究快速增殖核反应堆的安全问题。
Giesecke等提出的数值模拟快速增殖反应堆的液态钠在磁性材料的影响下的流动可能引起的磁场自激放大发电现象时,只考虑了局部磁雷诺数小于等于20,相对磁导率小于20的情况。本软件将采用的积分方程法,可用于磁雷诺数和相对磁导率均大于20的情况。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种采用积分方程法进行数值模拟磁感应现象的方法,其特征在于,包括以下步骤:
考虑可导区域D,其由m个子区域Di,其中i=1,2,3,…,m构成,子区域的表面记为Si;设Di子区域的导电率和导磁率分别是σi和μi;Sij表示Di和Dj的交界面;nij表示垂直于Sij表面的单位法矢量,且从Di指向Dj;设流场是稳态的,则磁场对时间的依赖关系表示为:~eλt,即磁场可表示为B(r,t)=B(r)eλt,λ的实部代表增长率,虚部代表频率;应用积分理论和边界元法,磁感应方程表示为如下积分方程组:
Br=LrkBk+λPrjAj+Nrlφl (4)
其中,Br代表磁场分量在r单元的离散值,方程(4)的右边的第一至第三项分别代表方程(1)右边的第一至第三个积分的离散形式,利用了爱因斯坦求和约定,重复的下标表示求和;同样方程(2)对应的离散方程为:
0.5φl+Elmφm=HlkBk+λDljAj (5)
或者为
Glmφm=HlkBk+λDljAj (6)
其中,Glm=0.5δlm+Elm;方程(3)对应的离散方程为:
Aj=QjkBk (7)
由方程(4)-(7)可得以磁场的分量在节点的值为未知数的代数方程组为:
Br=LrkBk+λPrjQjkBk+Nrm(G-1)mlHlkBk+λNrl(G-1)lmDmjQjkBk (8)
求解代数方程组(8)可得磁感应产生的磁场。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710590450.4A CN107391434B (zh) | 2017-07-18 | 2017-07-18 | 一种数值模拟磁感应现象的积分方程法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710590450.4A CN107391434B (zh) | 2017-07-18 | 2017-07-18 | 一种数值模拟磁感应现象的积分方程法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107391434A CN107391434A (zh) | 2017-11-24 |
CN107391434B true CN107391434B (zh) | 2020-10-30 |
Family
ID=60335921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710590450.4A Expired - Fee Related CN107391434B (zh) | 2017-07-18 | 2017-07-18 | 一种数值模拟磁感应现象的积分方程法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107391434B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105487025A (zh) * | 2015-11-27 | 2016-04-13 | 宁波南车时代传感技术有限公司 | 基于磁场积分方程法的磁通门瞬态分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8797702B2 (en) * | 2011-06-29 | 2014-08-05 | Toyota Motor Engineering & Manufacturing North America, Inc. | Focusing device for low frequency operation |
-
2017
- 2017-07-18 CN CN201710590450.4A patent/CN107391434B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105487025A (zh) * | 2015-11-27 | 2016-04-13 | 宁波南车时代传感技术有限公司 | 基于磁场积分方程法的磁通门瞬态分析方法 |
Non-Patent Citations (2)
Title |
---|
A Model for the Forward Problem in Magnetic Induction Tomography Using Boundary Integral Equations;M. H. Pham et al.;《IEEE Transactions on Magnetics》;20081031;第44卷(第10期);第2262-2267页 * |
大地电磁三维矢量有限元正演;严家斌 等;《吉林大学学报(地球科学版)》;20160930;第46卷(第5期);第1538-1549页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107391434A (zh) | 2017-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sun et al. | Diffusion-based coarse graining in hybrid continuum–discrete solvers: Applications in CFD–DEM | |
Kar et al. | The potential energy of interaction between dissimilar electrical double layers | |
Ohshima | Modified Henry function for the electrophoretic mobility of a charged spherical colloidal particle covered with an ion-penetrable uncharged polymer layer | |
Ohshima | Sedimentation potential in a concentrated suspension of spherical colloidal particles | |
Chandran et al. | Unsteady hydromagnetic free convection flow with heat flux and accelerated boundary motion | |
Bochkov et al. | Solving elliptic interface problems with jump conditions on Cartesian grids | |
Mazaheri et al. | Coupling behavior of the pH/temperature sensitive hydrogels for the inhomogeneous and homogeneous swelling | |
Kim et al. | Systematic analysis for electrical conductivity of network of conducting rods by Kirchhoff's laws and block matrices | |
Amani et al. | Mixed discrete least squares meshless method for planar elasticity problems using regular and irregular nodal distributions | |
Wang et al. | Spontaneous detachment of colloids from primary energy minima by Brownian diffusion | |
Saleh et al. | Plasma-waves evolution and propagation modeled by sixth order Ramani and coupled Ramani equations using symmetry methods | |
Guittet et al. | A Voronoi interface approach to cell aggregate electropermeabilization | |
CN107391434B (zh) | 一种数值模拟磁感应现象的积分方程法 | |
Korpusov et al. | Blow-up of solutions of a full non-linear equation of ion-sound waves in a plasma with non-coercive non-linearities | |
Gedik et al. | Unsteady flow of two-phase fluid in circular pipes under applied external magnetic and electrical fields | |
Haverkort et al. | Magnetohydrodynamic effects on insulating bubbles and inclusions in the continuous casting of steel | |
Johnson et al. | An analysis of electrophoresis of concentrated suspensions of colloidal particles | |
Talaat et al. | Simulation of the electric field and the GPR resulting from vertical-driven rods earthing system in a multi-layers earth structure | |
Abade et al. | Dynamics of permeable particles in concentrated suspensions | |
Oulaid et al. | Accurate boundary treatments for lattice Boltzmann simulations of electric fields and electro-kinetic applications | |
CN114417748B (zh) | 一种水上飞机的水气交混状态湍流场的模拟方法及装置 | |
Moleon Baca et al. | Calculation of ζ potential by non-equilibrium double layer theory in positive polystyrene model colloids | |
Roa et al. | Ion size effects on the electrokinetics of salt-free concentrated suspensions in ac fields | |
Cheng et al. | Electrophoresis of a soft particle normal to a plane | |
Ding et al. | Sedimentation velocity and potential in a suspension of charge-regulating colloidal spheres |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201030 Termination date: 20210718 |