CN110232205A - 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 - Google Patents
用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 Download PDFInfo
- Publication number
- CN110232205A CN110232205A CN201910348249.4A CN201910348249A CN110232205A CN 110232205 A CN110232205 A CN 110232205A CN 201910348249 A CN201910348249 A CN 201910348249A CN 110232205 A CN110232205 A CN 110232205A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- flux
- tokamak
- magnetic flux
- formula
- 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.)
- Granted
Links
Classifications
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Plasma Technology (AREA)
Abstract
本发明属于磁约束受控核聚变领域托卡马克装置放电的数值模拟,涉及一种用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法。根据托卡马克上线圈组中的电流数据,利用毕奥‑萨伐尔定律求解出感应磁场;根据大环径比近似下,积分得到感应磁场对应的极向磁通的表达式;利用傅里叶变换伪谱方法,把磁通变换成不同模数的分量,得到其在托卡马克中位形下的三维分布,并耦合到磁流体方程中,利用Crank‑Nicolson方法求解存在共振磁扰动下的磁场随时间的演化,以此来研究其对新经典撕裂模的影响。本发明能模拟出新经典撕裂模锁模、抑制和场穿透等物理过程,计算效率高、数值稳定性较好,是一种高效且能准确模拟实验的数值方法。
Description
技术领域
本发明属于磁约束受控核聚变领域托卡马克装置放电的数值模拟,特别涉及一种用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法。
背景技术
近几十年来由于人类大量的使用化石燃料造成了世界范围内严重的环境问题,比如全球变暖,雾霾和光污染等。特别是在以中国为代表的发展中国家,这些问题更加严重。经济的快速发展是牺牲环境所带来的。因此,发展清洁的可再生能源成为了当前中国的迫切需求。科学家们相信核聚变能是最有希望从根本上解决能源问题的途径之一。其中,磁约束核聚变被认为是目前最可行的方法。磁约束聚变研究目前遇到的最大困难就是要解决聚变装置托卡马克装置运行是存在的各种不稳定性。新经典撕裂模就是其中一种非常危险的宏观不稳定性,外加螺旋共振磁扰动是一种能有效控制新经典撕裂模的办法。实验中主要通过在托卡马克外部设置鞍形线圈组,通上交流电,产生感应磁场来和托卡马克芯部的等离子体相互作用。由于托卡马克装置巨大,聚变反应物理过程十分复杂,实验一次的费用极其昂贵且周期很长,实验操作不当还会损坏装置,造成巨大的经济损失。因此,能够对其中关键物理问题进行数值模拟研究的软件和方法就显得尤为重要了。
发明内容
为了解决上述问题,本发明提供一种用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法。本发明能模拟出新经典撕裂模锁模的抑制和场穿透等物理过程,采取的数值算法稳定性好、并行效率高,是一种高效且能准确模拟实验的数值方法。
一种用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法,具体步骤如下:
步骤1:根据所要模拟的托卡马克装置的大环的半径R0和小圆截面的半径a参数,建立三维坐标系(r,θ,z),构建三维的计算网格,其中,r、θ和z分别为径向、极向和轴向坐标;
步骤2:在托卡马克装置的外部设置多组磁扰动线圈,每个磁扰动线圈可实现正向通电、反向通电、不通电的不同工作方式,通过各工作方式之间的不同组合,利用毕奥-萨伐尔定律计算线圈在等离子体区域所产生的磁场B(r),计算公式如下:
其中,μ0为真空磁导率,C为积分路径,I为线圈电流强度,为源电流的微小线元,为电流元指向待求场点的单位向量,r为待求场点的坐标;
步骤3:将步骤2得到的磁场B(r),在托卡马克小圆截面上通过坐标变换转换为径向磁场Br和极向磁场Bθ的形式:
在大环径比近似下,托卡马克中的磁场写成如下形式:
其中,B为总磁场强度,B0为轴向磁场强度,为轴向单位矢量,ψ为极向磁通,用于表征极向磁场的物理量;
通过公式(2)得到径向磁场Br和极向磁场Bθ的表达式如下:
利用公式(3)和(4),通过对Br和Bθ进行积分计算求解出所等价的极向磁通ψ;
步骤4:把步骤3得到的等价的极向磁通ψ进行傅里叶展开,得到不同模式的磁通形式,即任意模式的共振磁扰动的磁通在三维空间上的分布,具体表示如下:
ψ(r,θ,z)=ψ0(r)+∑ψm,n(r)ei(mθ+nφ) (5)
其中,r、θ和z分别为径向、极向和轴向坐标,ωm,n为不同模式的磁通的本征函数,m和n分别为极向和环向模数;ψ0(r)表示模数为0的磁通,称为平衡磁通;i为虚数单位;
步骤5:将步骤4得到的共振磁扰动的磁通叠加到计算区域的边界,即计算网格的最外层磁面上,耦合到磁流体方程中,如下所示:
其中,u为涡量,ψ为磁通,p为等离子体压强,为等离子体电流,jb为自举电流,φ为流函数,R为雷诺数,SA为磁雷诺数,χ||和χ⊥分别为平行和垂直输运系数,Ez0和S0分别为电磁和压强源项。
采用Crank-Nicolson方法对涡量u、磁通ψ和等离子体压强p进行时间推进计算,得到在存在共振磁扰动下,涡量u、磁通ψ和等离子体压强p在计算网格上的空间的分布;
Crank-Nicolson方法的计算公式如下:
其中,U和F为任意函数,x为空间坐标,t为时间坐标,Δt为时间间隔。
为提高计算精度,对Crank-Nicolson方法进行改进,具体过程如下:
以磁通的时间推进为例:
假设已知t时刻磁通值ψt,需要求得t+Δt时刻的磁通ψt+Δt,步骤如下:
第一步:预计
根据公式(7)和ψt,求得下标rhs代表公式(7)中的右端项;通过公式(10),取时间间隔为计算得到时刻的磁通利用公式(7)和求得时刻的右端项
第二步:校正
回滚到t时刻,取时间间隔为Δt,利用公式(11)和第一步预计得到的右手项求得t+Δt时刻的磁通ψt+Δt。
步骤6:保存步骤5中得到的涡量u、磁通ψ和等离子体压强p的三维空间分布信息;
步骤7:重复进行步骤2-6直至时间t等于设定的最大计算时间Tmax,即得到任意时刻的涡量u、磁通ψ和等离子体压强p等物理量的变化。
本发明的有益效果:本发明提供了一种能够模拟外加共振磁扰动对新经典撕裂模影响的方法,能模拟出新经典撕裂模锁模,抑制和场穿透等物理过程。采用了傅里叶谱方法,与传统的直接在边界上增加磁扰动的方法相比,能够避免边界处理不当带来的数值不稳定性,实现长时间的非线性演化。同时,在托卡马克实验上通常也是采用谱分析的方法,这样的数值算法能更好的结合实验,计算效率高,数值稳定性较好,是一种的高效且有准确的数值模拟方法。
附图说明
图1是本发明所适用的托卡马克实验装置的示意图。
图2是本发明用于模拟计算采用的网格。
图3是三维磁岛结构示意图。
图4是本发明模拟的2/1共振磁扰动的三维分布图。
图5(a)、图5(b)和图5(c)分别是本发明关于新经典撕裂模抑制、小磁岛锁模和大磁岛锁模的模拟结果图。
图6是本发明的方法流程图。
具体实施方式
下面结合附图和技术方案,进一步说明本发明的具体实施方式。
托卡马克装置是一种通过强磁场把热等离子体约束在一个环形容器内发生反应的装置,其结构如图1所示,其中外部结构就是外加的共振磁扰动线圈。首先,根据真实的装置位形建立三维坐标系(r,θ,z),r、θ和z分别为径向、极向和轴向坐标,构建三维的计算网格,如图2所示。托卡马克中初始的平衡磁面是层层嵌套的圆或有一定形变的圆,新经典撕裂模会撕裂磁面,形成磁岛结构,如图3所示。然后,根据实验中磁扰动线圈中所通交流电的数据,求解出其所产生的感应磁场。接下来,积分得到感应磁场对应的极向磁通的表达式。利用傅里叶变换伪谱方法,把磁通变换成不同模数的分量,通过一定的选择和滤波,得到其在托卡马克中位形下的三维分布,图4即为求得的一例三维共振磁扰动磁通的示意图。再把扰动磁通添加到计算区域的边界即最外层磁面上,耦合到磁流体方程中,计算出存在扰动磁场时,下一时刻各个物理量的值。重复之前的步骤多次,以此来得到任意时刻各物理量随时间的非线性演化,从而分析得到其对新经典撕裂模的影响和其中的物理过程。部分模拟结果如图5(a)、图5(b)和图5(c)所示,分别为新经典撕裂模抑制、小磁岛锁模和大磁岛锁模。
具体实施步骤如下:
步骤1:根据所要模拟的托卡马克装置的大环的半径R0,小圆截面的半径a等参数,建立三维坐标系(r,θ,z),r、θ和z分别为径向、极向和轴向坐标,构建三维的计算网格,如图2所示;
步骤2:在托卡马克装置的外部设置多组磁扰动线圈,每个磁扰动线圈可实现正向通电、反向通电、不通电的不同工作方式,通过各工作方式之间的不同组合,利用毕奥-萨伐尔定律计算出线圈在等离子体区域所产生的磁场,计算公式为:
其中,μ0为真空磁导率,C为积分路径,I为线圈电流强度,为源电流的微小线元,为电流元指向待求场点的单位向量,r为待求场点的坐标;
步骤3:把步骤2得到的磁场在托卡马克小圆截面上通过坐标变换转换成Br和Bθ的形式:
在大环径比近似下,托卡马克中的磁场可以写成如下形式:
其中,B为总磁场强度,B0为轴向磁场强度,为轴向单位矢量,ψ为极向磁通,用于表征极向磁场的物理量;
通过公式(2)可以得到径向磁场Br和极向磁场Bθ的表达式如下:
利用公式(3)和(4),通过对Br和Bθ进行积分计算求解出所等价的极向磁通ψ;
步骤4:把步骤3得到的等价的极向磁通进行傅里叶展开,得到不同模式的磁通形式,即任意模式的共振磁扰动的磁通在三维空间上的分布,具体表示如下:
ψ(r,θ,z)=ψ0(r)+∑ψm,n(r)ei(mθ+nφ) (5)
其中,r、θ和z分别为径向、极向和轴向坐标,ψm,n为不同模式的磁通的本征函数,m和n分别为极向和环向模数;
步骤5:将步骤4得到的共振磁扰动的磁通叠加到计算区域的边界,即计算网格的最外层磁面上,耦合到磁流体方程,如下:
其中,u为涡量,ψ为磁通,p为等离子体压强,为等离子体电流,jb为自举电流,φ为流函数,R为雷诺数,SA为磁雷诺数,χ||和χ⊥分别为平行和垂直输运系数,Ez0和S0分别为电磁和压强源项。
采用Crank-Nicolson方法对涡量u、磁通ψ和等离子体压强p进行时间推进计算,得到在存在共振磁扰动下,涡量u、磁通ψ和等离子体压强p在计算网格上的空间的分布;
Crank-Nicolson方法的计算公式如下:
其中,U,F为任意函数,x为空间坐标,t为时间坐标,Δt为时间间隔。
为提高计算精度,对Crank-Nicolson方法进行改进,具体过程如下:
以磁通的时间推进为例:
假设现在已知t时刻磁通值ψt,想要求得t+Δt时刻的磁通ψt+Δt,第一步:预计
根据公式(7)和ψt,求得下标rhs代表公式(7)中的右端项。通过公式(10),取时间间隔为计算得到时刻的磁通利用公式(7)和求得时刻的右端项
第二步:校正
回滚到t时刻,取时间间隔为Δt,利用公式(11)和第一步预计得到的右手项便可求得t+Δt时刻的磁通ψt+Δt。
步骤6:保存步骤5中得到的涡量u、磁通ψ和等离子体压强p的三维空间分布信息;
步骤7:重复进行步骤2-6直至时间t等于设定的最大计算时间Tmax,即得到任意时刻的涡量u、磁通ψ和等离子体压强p等物理量的变化。
Claims (1)
1.一种用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法,其特征在于,具体步骤如下:
步骤1:根据所要模拟的托卡马克装置的大环的半径R0和小圆截面的半径a参数,建立三维坐标系(r,θ,z),构建三维的计算网格,其中,r、θ和z分别为径向、极向和轴向坐标;
步骤2:在托卡马克装置的外部设置多组磁扰动线圈,每个磁扰动线圈可实现正向通电、反向通电、不通电的不同工作方式,通过各工作方式之间的不同组合,利用毕奥-萨伐尔定律计算线圈在等离子体区域所产生的磁场B(r),计算公式如下:
其中,μ0为真空磁导率,C为积分路径,I为线圈电流强度,为源电流的微小线元,为电流元指向待求场点的单位向量,r为待求场点的坐标;
步骤3:将步骤2得到的磁场B(r),在托卡马克小圆截面上通过坐标变换转换为径向磁场Br和极向磁场Bθ的形式:
在大环径比近似下,托卡马克中的磁场写成如下形式:
其中,B为总磁场强度,B0为轴向磁场强度,为轴向单位矢量,ψ为极向磁通,用于表征极向磁场的物理量;
通过公式(2)得到径向磁场Br和极向磁场Bθ的表达式如下:
利用公式(3)和(4),通过对Br和Bθ进行积分计算求解出所等价的极向磁通ψ;
步骤4:把步骤3得到的等价的极向磁通ψ进行傅里叶展开,得到不同模式的磁通形式,即任意模式的共振磁扰动的磁通在三维空间上的分布,具体表示如下:
ψ(r,θ,z)=ψ0(r)+∑ψm,n(r)ei(mθ+nφ) (5)
其中,r、θ和z分别为径向、极向和轴向坐标,ψm,n为不同模式的磁通的本征函数,m和n分别为极向和环向模数;ψ0(r)表示模数为0的磁通,称为平衡磁通;i为虚数单位;
步骤5:将步骤4得到的共振磁扰动的磁通叠加到计算区域的边界,即计算网格的最外层磁面上,耦合到磁流体方程中,如下所示:
其中,u为涡量,ψ为磁通,p为等离子体压强,为等离子体电流,jb为自举电流,φ为流函数,R为雷诺数,SA为磁雷诺数,χ||和χ⊥分别为平行和垂直输运系数,Ez0和S0分别为电磁和压强源项;
采用Crank-Nicolson方法对涡量u、磁通ψ和等离子体压强p进行时间推进计算,得到在存在共振磁扰动下,涡量u、磁通ψ和等离子体压强p在计算网格上的空间的分布;为提高计算精度,对Crank-Nicolson方法进行改进,具体过程如下:
以磁通的时间推进为例:
假设已知t时刻磁通值ψt,需要求得t+Δt时刻的磁通ψt+Δt,步骤如下:
第一步:预计
根据公式(7)和ψt,求得下标rhs代表公式(7)中的右端项;通过公式(10),取时间间隔为计算得到时刻的磁通利用公式(7)和求得时刻的右端项
第二步:校正
回滚到t时刻,取时间间隔为Δt,利用公式(11)和第一步预计得到的右手项求得t+Δt时刻的磁通ψt+Δt;
步骤6:保存步骤5中得到的涡量u、磁通ψ和等离子体压强p的三维空间分布信息;
步骤7:重复进行步骤2-6直至时间t等于设定的最大计算时间Tmax,即得到任意时刻的涡量u、磁通ψ和等离子体压强p物理量的变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910348249.4A CN110232205B (zh) | 2019-04-28 | 2019-04-28 | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910348249.4A CN110232205B (zh) | 2019-04-28 | 2019-04-28 | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110232205A true CN110232205A (zh) | 2019-09-13 |
CN110232205B CN110232205B (zh) | 2020-08-25 |
Family
ID=67860276
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910348249.4A Active CN110232205B (zh) | 2019-04-28 | 2019-04-28 | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110232205B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259599A (zh) * | 2020-01-07 | 2020-06-09 | 大连理工大学 | 一种用于分析托卡马克中高能量粒子约束性能的测试粒子模拟方法 |
CN112309588A (zh) * | 2020-10-12 | 2021-02-02 | 中国科学院合肥物质科学研究院 | 一种适用于磁约束聚变装置的共振磁扰动线圈及实现方法 |
CN112632892A (zh) * | 2020-11-27 | 2021-04-09 | 浙江大学 | 一种等离子体中具有刮削层特性的偏滤器平衡位形的模拟方法 |
CN112733363A (zh) * | 2021-01-11 | 2021-04-30 | 大连理工大学 | 一种用于east托卡马克中三维全空间等离子体响应的模拟方法 |
CN112766386A (zh) * | 2021-01-25 | 2021-05-07 | 大连理工大学 | 一种基于多输入多输出融合网络的广义零样本学习方法 |
CN112786273A (zh) * | 2020-12-31 | 2021-05-11 | 中国科学院合肥物质科学研究院 | 一种基于立方体永磁体块的仿星器磁体及其优化布置方法 |
CN112992385A (zh) * | 2021-05-06 | 2021-06-18 | 西南交通大学 | 一种准环对称仿星器磁场位形的设计方法 |
CN113283135A (zh) * | 2021-05-17 | 2021-08-20 | 核工业西南物理研究院 | 一种托卡马克真空室电磁仿真等离子体电流数据提取方法 |
CN113536211A (zh) * | 2021-07-24 | 2021-10-22 | 大连理工大学 | 一种用于east托卡马克磁信号快速提取的数值方法 |
WO2022193041A1 (zh) * | 2021-03-15 | 2022-09-22 | 大连理工大学 | 一种用于east托卡马克射频波导致电子温度演化的模拟方法 |
CN115479900A (zh) * | 2022-03-24 | 2022-12-16 | 核工业西南物理研究院 | 一种聚变装置玻璃观察窗污染防护装置 |
CN116167247A (zh) * | 2023-04-21 | 2023-05-26 | 中国科学院合肥物质科学研究院 | 一种基于冯·哈格诺夫方法的gs方程数值计算方法 |
CN117371299A (zh) * | 2023-12-08 | 2024-01-09 | 安徽大学 | 一种托卡马克新经典环向粘滞力矩的机器学习方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5120487A (en) * | 1990-03-01 | 1992-06-09 | The United States Of America As Represented By The United States Department Of Energy | Method of measuring the dc electric field and other tokamak parameters |
CN104332184A (zh) * | 2014-08-30 | 2015-02-04 | 中国科学院等离子体物理研究所 | 适用于east托卡马克装置全钨偏滤器的靶板探针系统 |
CN107068204A (zh) * | 2011-11-14 | 2017-08-18 | 加州大学评议会 | 用于形成并维持高性能frc的系统和方法 |
CN108288509A (zh) * | 2018-01-10 | 2018-07-17 | 成都大学 | 基于气体注入增强磁扰动抑制逃逸电流的方法及装置 |
CN108733946A (zh) * | 2018-05-24 | 2018-11-02 | 大连理工大学 | 一种用于托卡马克中电子回旋电流剖面演化的模拟方法 |
-
2019
- 2019-04-28 CN CN201910348249.4A patent/CN110232205B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5120487A (en) * | 1990-03-01 | 1992-06-09 | The United States Of America As Represented By The United States Department Of Energy | Method of measuring the dc electric field and other tokamak parameters |
CN107068204A (zh) * | 2011-11-14 | 2017-08-18 | 加州大学评议会 | 用于形成并维持高性能frc的系统和方法 |
CN104332184A (zh) * | 2014-08-30 | 2015-02-04 | 中国科学院等离子体物理研究所 | 适用于east托卡马克装置全钨偏滤器的靶板探针系统 |
CN108288509A (zh) * | 2018-01-10 | 2018-07-17 | 成都大学 | 基于气体注入增强磁扰动抑制逃逸电流的方法及装置 |
CN108733946A (zh) * | 2018-05-24 | 2018-11-02 | 大连理工大学 | 一种用于托卡马克中电子回旋电流剖面演化的模拟方法 |
Non-Patent Citations (1)
Title |
---|
WEIKANG TANG, ET AL.: "Effects of resonant magnetic perturbation on locked mode of neoclassical tearing modes", 《PLASMA SCIENCE AND TECHNOLOGY》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259599B (zh) * | 2020-01-07 | 2022-09-20 | 大连理工大学 | 一种用于分析托卡马克中高能量粒子约束性能的测试粒子模拟方法 |
CN111259599A (zh) * | 2020-01-07 | 2020-06-09 | 大连理工大学 | 一种用于分析托卡马克中高能量粒子约束性能的测试粒子模拟方法 |
CN112309588A (zh) * | 2020-10-12 | 2021-02-02 | 中国科学院合肥物质科学研究院 | 一种适用于磁约束聚变装置的共振磁扰动线圈及实现方法 |
CN112309588B (zh) * | 2020-10-12 | 2022-10-21 | 中国科学院合肥物质科学研究院 | 一种适用于磁约束聚变装置的共振磁扰动线圈及实现方法 |
CN112632892B (zh) * | 2020-11-27 | 2022-03-25 | 浙江大学 | 一种等离子体中具有刮削层特性的偏滤器平衡位形的模拟方法 |
CN112632892A (zh) * | 2020-11-27 | 2021-04-09 | 浙江大学 | 一种等离子体中具有刮削层特性的偏滤器平衡位形的模拟方法 |
CN112786273A (zh) * | 2020-12-31 | 2021-05-11 | 中国科学院合肥物质科学研究院 | 一种基于立方体永磁体块的仿星器磁体及其优化布置方法 |
CN112733363A (zh) * | 2021-01-11 | 2021-04-30 | 大连理工大学 | 一种用于east托卡马克中三维全空间等离子体响应的模拟方法 |
CN112733363B (zh) * | 2021-01-11 | 2022-09-20 | 大连理工大学 | 一种用于east托卡马克中三维全空间等离子体响应的模拟方法 |
WO2022148333A1 (zh) * | 2021-01-11 | 2022-07-14 | 大连理工大学 | 一种用于east托卡马克中三维全空间等离子体响应的模拟方法 |
CN112766386A (zh) * | 2021-01-25 | 2021-05-07 | 大连理工大学 | 一种基于多输入多输出融合网络的广义零样本学习方法 |
CN112766386B (zh) * | 2021-01-25 | 2022-09-20 | 大连理工大学 | 一种基于多输入多输出融合网络的广义零样本学习方法 |
WO2022193041A1 (zh) * | 2021-03-15 | 2022-09-22 | 大连理工大学 | 一种用于east托卡马克射频波导致电子温度演化的模拟方法 |
CN112992385A (zh) * | 2021-05-06 | 2021-06-18 | 西南交通大学 | 一种准环对称仿星器磁场位形的设计方法 |
CN112992385B (zh) * | 2021-05-06 | 2021-08-03 | 西南交通大学 | 一种准环对称仿星器磁场位形的设计方法 |
CN113283135A (zh) * | 2021-05-17 | 2021-08-20 | 核工业西南物理研究院 | 一种托卡马克真空室电磁仿真等离子体电流数据提取方法 |
CN113283135B (zh) * | 2021-05-17 | 2023-02-21 | 核工业西南物理研究院 | 一种托卡马克真空室电磁仿真等离子体电流数据提取方法 |
CN113536211A (zh) * | 2021-07-24 | 2021-10-22 | 大连理工大学 | 一种用于east托卡马克磁信号快速提取的数值方法 |
CN113536211B (zh) * | 2021-07-24 | 2024-02-06 | 大连理工大学 | 一种用于east托卡马克磁信号快速提取的数值方法 |
CN115479900A (zh) * | 2022-03-24 | 2022-12-16 | 核工业西南物理研究院 | 一种聚变装置玻璃观察窗污染防护装置 |
CN116167247A (zh) * | 2023-04-21 | 2023-05-26 | 中国科学院合肥物质科学研究院 | 一种基于冯·哈格诺夫方法的gs方程数值计算方法 |
CN117371299A (zh) * | 2023-12-08 | 2024-01-09 | 安徽大学 | 一种托卡马克新经典环向粘滞力矩的机器学习方法 |
CN117371299B (zh) * | 2023-12-08 | 2024-02-27 | 安徽大学 | 一种托卡马克新经典环向粘滞力矩的机器学习方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110232205B (zh) | 2020-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110232205A (zh) | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 | |
Powell et al. | A solution-adaptive upwind scheme for ideal magnetohydrodynamics | |
Srinivasan et al. | Magnetic field generation in Rayleigh-Taylor unstable inertial confinement fusion plasmas | |
CN104331544B (zh) | 一种基于eic原理的三相三柱变压器建模方法 | |
CN102346442B (zh) | 适于节点分析框架的微网暂态仿真并行计算方法 | |
Boyer et al. | Central safety factor and βN control on NSTX-U via beam power and plasma boundary shape modification, using TRANSP for closed loop simulations | |
Lore et al. | Implementation of the 3D edge plasma code EMC3-EIRENE on NSTX | |
CN105160115B (zh) | 基于近似与灵敏度分析的反射面天线机电集成优化设计方法 | |
Stanier et al. | Solar particle acceleration at reconnecting 3D null points | |
McGreivy et al. | Optimized finite-build stellarator coils using automatic differentiation | |
Gourgouliatos et al. | Magnetar activity via the density–shear instability in Hall-MHD | |
Gaye et al. | H∞ stabilization of the current profile in tokamak plasmas via an LMI approach | |
CN103150581B (zh) | 基于最小二乘支持向量机燃烧模型的锅炉优化方法和装置 | |
Romero et al. | Plasma internal inductance dynamics in a tokamak | |
Sharma et al. | Modeling and control of TCV | |
Ambrosino et al. | Optimization of experimental snowflake configurations on TCV | |
Gu et al. | AC losses in HTS tapes and devices with transport current solved through the resistivity-adaption algorithm | |
Beach | Master equation approach to computing RVB bond amplitudes | |
Suzuki et al. | Monte Carlo simulations of the two-dimensional quantal and classical spin systems—A new type of phase transition with vortices | |
Coleman et al. | The design and optimisation of tokamak poloidal field systems in the BLUEPRINT framework | |
Gerard et al. | Optimizing the HSX stellarator for microinstability by coil-current adjustments | |
CN106777472A (zh) | 基于拉盖尔多项式的减少分裂误差的完全匹配层实现方法 | |
Gonzalez-Martin et al. | Modelling the Alfvén eigenmode induced fast-ion flow measured by an imaging neutral particle analyzer | |
Hyatt et al. | Designing, constructing, and using plasma control system algorithms on DIII-D | |
Zhang et al. | Analysis of the impacts of the J-TEXT TBM mock-up on the equilibrium magnetic field |
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 |