CN117371299A - 一种托卡马克新经典环向粘滞力矩的机器学习方法 - Google Patents
一种托卡马克新经典环向粘滞力矩的机器学习方法 Download PDFInfo
- Publication number
- CN117371299A CN117371299A CN202311674834.6A CN202311674834A CN117371299A CN 117371299 A CN117371299 A CN 117371299A CN 202311674834 A CN202311674834 A CN 202311674834A CN 117371299 A CN117371299 A CN 117371299A
- Authority
- CN
- China
- Prior art keywords
- ntv
- tokamak
- input
- viscous torque
- classical
- 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
- 238000010801 machine learning Methods 0.000 title claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 29
- 238000012549 training Methods 0.000 claims abstract description 28
- 230000009466 transformation Effects 0.000 claims abstract description 19
- 230000006870 function Effects 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 14
- 238000013528 artificial neural network Methods 0.000 claims abstract description 12
- 238000003062 neural network model Methods 0.000 claims abstract description 8
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 238000000034 method Methods 0.000 claims description 36
- 238000009826 distribution Methods 0.000 claims description 29
- 238000002474 experimental method Methods 0.000 claims description 13
- 238000005315 distribution function Methods 0.000 claims description 11
- 238000012360 testing method Methods 0.000 claims description 10
- 239000002245 particle Substances 0.000 claims description 9
- 230000003416 augmentation Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 238000011156 evaluation Methods 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 4
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 claims description 3
- 230000008859 change Effects 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 230000010354 integration Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 5
- 238000013461 design Methods 0.000 abstract description 3
- 230000004927 fusion Effects 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
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
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- 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
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/0499—Feedforward networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Artificial Intelligence (AREA)
- Software Systems (AREA)
- Medical Informatics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及可控核聚变技术领域,具体公开了一种托卡马克新经典环向粘滞力矩的机器学习方法,包括以下步骤:S1、获取传统托卡马克新经典环向粘滞力矩计算程序输入输出,构造数据集;S2、根据物理特征设计输入输出变换函数,使之适用于深度神经网络模型训练;S3、设计深度神经网络算法,执行模型训练,得到预测模型;S4、对所述的预测模型性能进行评估及评价。本发明保证了新经典环向粘滞力矩模拟预测的准确性,避免了简化模型可能带来的预测偏差,同时又显著提高了计算效率、降低了计算成本,使得相关模拟分析能够在即时或实时的条件下开展,能够对物理分析、实验方案设计、实验放电控制提供准确且即时的理论参考与建议。
Description
技术领域
本发明涉及可控核聚变技术领域,具体涉及一种托卡马克新经典环向粘滞力矩的机器学习方法。
背景技术
理想的托卡马克磁场位形是环向对称的,然而在真实实验条件下,由于线圈分立性、误差场、外加磁场扰动等因素,托卡马克磁场往往呈现出环向不对称性。这种环向不对称的磁场扰动将通过新经典环向粘滞(NTV)效应产生额外的环向力矩,即NTV力矩,影响等离子体环向旋转,进而影响等离子体不稳定性及装置性能表现。因此,准确且快速的NTV力矩模拟分析有助于增强对实验现象的理解、提升对未来实验的预测和控制能力。
然而,现有的托卡马克新经典环向粘滞力矩模拟方法无法兼顾计算准确性与计算效率:即如果要实现较为准确的NTV力矩模拟,所需要的计算时间无法满足即时甚至实时反馈的要求;采用简化的NTV力矩物理模型能够极大提升计算效率,但简化模型的计算结果可能与真实情况存在较大偏差。计算准确性与计算效率之间的矛盾制约了物理研究的开展与装置性能的优化。
发明内容
本发明的目的在于提供一种托卡马克新经典环向粘滞力矩的机器学习方法,解决上述技术问题。
本发明的目的可以通过以下技术方案实现:
一种托卡马克新经典环向粘滞力矩的机器学习方法,包括以下步骤:
S1、获取传统托卡马克新经典环向粘滞力矩计算程序输入输出,构造数据集;
S2、根据物理特征设计输入输出变换函数,使之适用于深度神经网络模型训练;
S3、设计深度神经网络算法,执行模型训练,得到预测模型;
S4、对所述的预测模型性能进行评估及评价。
作为本发明进一步的方案:在所述的步骤S1中,获取传统托卡马克新经典环向粘滞力矩计算程序输入输出,构造数据集的具体步骤如下所示:
从实验数据库中选取若干具有不同边界安全因子值的若干个放电实验,并获取所述的放电实验中的等离子体平衡位形、安全因子径向剖面分布、电子密度径向剖面分布、电子温度径向剖面分布、离子温度径向剖面分布和等离子体环向旋转频率径向剖面分布;
基于电子密度、电子温度、离子温度、等离子体环向旋转频率4种物理量的径向剖面分布,通过同时改变幅度和剖面梯度的方式进行剖面增广,每种物理量增广后得到5个剖面分布;
对增广后的物理量剖面分布进行组合,共计得到14×54=8750组等离子体参数位形;
针对上述8750组等离子体参数位形的每一种组合,利用传统新经典环向粘滞力矩数值计算方法,执行新经典环向粘滞力矩计算;
传统新经典环向粘滞力矩数值计算方法步骤为,从线性化漂移动力学方程出发:
;
其中f 0为平衡分布函数,为抛射角散射碰撞算符,/>为粒子速度矢量,(ρ,α,θ)为空间磁面坐标系,v d为偏向频率,f为扰动分布函数,即为本方程的解;
进行坐标变换,从坐标系变换到/>坐标系,其中/>为抛射角变量,x为粒子归一化能量,新坐标定义为:
;
;
;
;
其中表示反弹平均,Bθ表示极向磁场,/>,,BM与Bm分别表示磁面上最大和最小磁场强度。
在坐标系中对扰动分布函数和线性化漂移动力学方程进行傅里叶分解,得到:
;
;
其中一些符号定义为:
;
;
;
;
;
;
T为温度,e为电荷,为极向磁通的径向导数,/>为扰动磁场的傅里叶分量,qωE为电漂移频率,ωB为磁漂移频率,DB包含磁漂移频率的极向依赖关系,ωb为粒子反弹频率。
求解得到扰动分布函数后,即可通过如下相空间积分过程得到NTV力矩径向剖面分布,即/>:
;
;
;
;
;
;
;
;
;
;
;
N j 、T j 、e j 分别为第j种粒子的密度、温度、电荷,P为等离子体压强。
经过上述计算步骤,能够得到8750组数据样本。对其按8:2的比例随机划分为训练集和测试集。
作为本发明进一步的方案:在所述的步骤S2中,根据物理特征设计输入输出变换函数,使之适用于深度神经网络模型训练,具体步骤包括:
对每种物理量,使用其典型取值进行无量纲化处理,即:
输入向量中的电子密度=原始电子密度/1019;
输入向量中的电子温度=原始电子温度/103;
输入向量中的离子温度=原始离子温度/103;
输入向量中的等离子体旋转频率=原始等离子体旋转频率/104;
输入向量中的安全因子=原始安全因子/10;
由此能够保证输入物理量取值范围均在1-10量级;
托卡马克典型NTV力矩取值范围约为10-5-1N/m2不等,与输入物理量类似,NTV力矩取值范围同样跨越多个数量级尺度,因此,采用如下变换函数对NTV力矩进行变换:
;
其中T NTV 为变换前的NTV力矩,T' NTV 为变换后的NTV力矩,变换后的NTV力矩取值范围约为1-10的量级。
作为本发明进一步的方案:在所述的步骤S4中,计深度神经网络算法,执行模型训练,得到预测模型,具体包括:
将托卡马克背景磁场B0、变换后的电子密度N e 、电子温度T e 、离子温度T i 、旋转频率、安全因子q组成一个输入向量,按照典型的输入物理量分辨率,输入向量长度为1006;
将变换后的电子NTV力矩、离子NTV力矩组成一个输出向量,按照典型的输出物理量分辨率,输出向量长度为130;
所述深度神经网络包含一个输入层,三个隐藏层,一个输出层,各层之间为全连接,其中输入层包含1006个特征,输出层包含130个特征,三个隐藏层分别包含500、520、260个特征;
设置学习率为0.001,采用Adam优化器,设置epochs为50000;
训练结束得到新经典环向粘滞力矩快速预测模型。
作为本发明进一步的方案:在所述的步骤S4中,模型性能评估与评价,具体包括:
观察损失函数随训练过程的变化;
在测试集上计算R2作为模型性能评估指标,R2定义为:
;
其中y (i)表示第i个样本真实值,y (i)表示第i个样本预测值,表示y (i)的平均值,R2越接近1,模型性能越好;
随机选取测试集中的某些样本,直观比较NTV力矩预测值与真实值的差异;
比较传统NTV力矩模拟方法与本方法计算时间差异。
本发明的有益效果:通过采用上述的技术方案,基于新经典环向粘滞力矩物理模型数据库,利用深度神经网络技术,实现了兼顾计算准确性和计算效率的托卡马克新经典环向粘滞力矩快速模拟预测。通过上述方法,保证了新经典环向粘滞力矩模拟预测的准确性,避免了简化模型可能带来的预测偏差,同时又显著提高了计算效率、降低了计算成本,使得相关模拟分析能够在即时或实时的条件下开展,能够对物理分析、实验方案设计、实验放电控制提供准确且即时的理论参考与建议。
附图说明
下面结合附图对本发明作进一步的说明。
图1是本发明一种实施例中提供的一种托卡马克新经典环向粘滞力矩的机器学习方法的流程示意图;
图2是本发明一种实施例中提供的深度神经网络结构示意图;
图3是本发明一种实施例中提供的模型训练过程示意图;
图4是本发明一种实施例中提供的模型预测表现示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
请参阅图1-4所示,本发明为一种托卡马克新经典环向粘滞力矩的机器学习方法,包括:
S1、获取传统托卡马克新经典环向粘滞力矩计算程序输入输出,构造数据集;具体步骤包括:
S11、从东方超环(EAST)托卡马克实验数据库中选取若干具有不同边界安全因子值的放电实验,下文叙述中以14个放电实验为例;
S12、从实验数据库中获取上述14个放电实验的等离子体平衡位形、安全因子径向剖面分布、电子密度径向剖面分布、电子温度径向剖面分布、离子温度径向剖面分布、等离子体环向旋转频率径向剖面分布;
S13、对于上述14个放电实验,基于电子密度、电子温度、离子温度、等离子体环向旋转频率4种物理量的径向剖面分布,通过同时改变幅度和剖面梯度的方式进行剖面增广,每种物理量增广后得到5个剖面分布;
S14、增广后的物理量剖面分布进行组合,共计得到组等离子体参数位形;
S15、针对上述8750组等离子体参数位形的每一种组合,利用传统新经典环向粘滞力矩数值计算方法,执行新经典环向粘滞力矩计算;
S16、传统新经典环向粘滞力矩数值计算方法步骤为,从线性化漂移动力学方程出发:
;
其中f 0为平衡分布函数,为抛射角散射碰撞算符,/>为粒子速度矢量,(ρ,α,θ)为空间磁面坐标系,v d为偏向频率,f为扰动分布函数,即为本方程的解;
进行坐标变换,从坐标系变换到/>坐标系,坐标变换满足如下关系:
;
;
;
在坐标系中对扰动分布函数和线性化漂移动力学方程进行傅里叶分解,得到:
;
;
其中一些符号定义为:
;
;
;
;
求解得到扰动分布函数后,即可通过如下相空间积分过程得到NTV力矩径向剖面分布,即/>:
;
;
;
经过上述计算步骤,能够得到8750组数据样本。对其按8:2的比例随机划分为训练集和测试集。
S2、根据物理特征设计输入输出变换函数,使之适用于深度神经网络模型训练;具体步骤包括:
S21、托卡马克典型电子密度约为10-19 m -3量级,典型电子温度和离子温度约为10- 3 eV量级,典型等离子体环向旋转频率约为104 rad/s量级,典型安全因子约为1-10量级,输入物理量跨越很大的数量级尺度,同时作为输入会造成输入向量变化剧烈,影响模型训练效果,因此对每种物理量,使用其典型取值进行无量纲化处理,即:
输入向量中的电子密度=原始电子密度/1019;
输入向量中的电子温度=原始电子温度/103;
输入向量中的离子温度=原始离子温度/103;
输入向量中的等离子体旋转频率=原始等离子体旋转频率/104;
输入向量中的安全因子=原始安全因子/10;
由此能够保证输入物理量取值范围均在1-10量级;
S22、托卡马克典型NTV力矩取值范围约为10-5-1N/m2不等,与输入物理量类似,NTV力矩取值范围同样跨越多个数量级尺度,因此,采用如下变换函数对NTV力矩进行变换:
;
其中TNTV为变换前的NTV力矩,T'NTV为变换后的NTV力矩,变换后的NTV力矩取值范围约为1-10的量级。
S3、设计深度神经网络算法,执行模型训练,得到预测模型;具体步骤包括:
S31、将托卡马克背景磁场B0、变换后的电子密度N e 、电子温度T e 、离子温度T i 、旋转频率、安全因子q组成一个输入向量,按照典型的输入物理量分辨率,输入向量长度为1006,如图2中向量/>所示;
S32、将变换后的电子NTV力矩、离子NTV力矩组成一个输出向量,按照典型的输出物理量分辨率,输出向量长度为130,如图2中向量所示;
S33、所述深度神经网络包含一个输入层,三个隐藏层,一个输出层,各层之间为全连接,其中输入层包含1006个特征(即输入向量长度),输出层包含130个特征(即输出向量长度),三个隐藏层(如图2中、/>、/>所示)分别包含500、520、260个特征,网络结构如图2所示;
S34、设置学习率为0.001,采用Adam优化器,设置epochs为50000;
S35、采用一块NVIDIA GeForce RTX 3060 Laptop GPU加速训练;
S36、训练结束得到新经典环向粘滞力矩快速预测模型。
S4、评估模型性能;具体步骤包括:
S41、观察损失函数随训练过程的变化。如图3所示,损失函数随训练过程进行持续减小;
S42、在测试集上计算R2作为模型性能评估指标,R2定义为:
;
其中y (i)表示第i个样本真实值,y (i)表示第i个样本预测值,表示y (i)的平均值,R2越接近1,模型性能越好,对于本实施例,R2≈0.976;
S43、随机选取测试集中的某些样本,直观比较NTV力矩预测值与真实值的差异。图4展示了测试集中编号分别为1、500、1000、1500的样本模型预测结果与真实值的剖面分布比较,可以发现本方法能够较为准确地实现NTV力矩的预测;
S44、比较传统NTV力矩模拟方法与本方法计算时间差异。传统NTV力矩模拟方法单次计算时间约为11.98s,本实施例中的方法单次计算时间约为2.86×10-7,即,计算效率得到了107~108倍的提升。
并且在本方法中,根据一般深度神经网络模型设计惯例,为保证模型预测准确性,将满足R2>0.95的模型作为一级模型,预测准确性较好;将0.95>R2>0.90的模型作为二级模型,预测准确性次之;R2<0.90的模型预测误差较大,不建议使用,需要通过修改模型参数,例如学习率、优化器、epochs等参数进行模型优化与性能提升,直至R2满足相应要求。
本方法所获得的深度神经网络模型在保证模拟准确性(R2)的前提下,预测推理速度极快,克服了传统模拟方法无法兼顾计算准确性与计算效率的难题。若要更进一步提升模型预测准确性,可以从扩大数据集、增加神经网络深度和宽度、增加训练次数、调整学习率等方面开展研究,这可能会增加模型训练时间,但一旦训练结束得到模型,模型预测推理时间几乎不会受到影响。使用者根据自身需求,可以基于上述思路开展更进一步的探索和补充。
以上对本发明的一个实施例进行了详细说明,但所述内容仅为本发明的较佳实施例,不能被认为用于限定本发明的实施范围。凡依本发明申请范围所作的均等变化与改进等,均应仍归属于本发明的专利涵盖范围之内。
Claims (5)
1.一种托卡马克新经典环向粘滞力矩的机器学习方法,其特征在于,包括以下步骤:
S1、获取传统托卡马克新经典环向粘滞力矩计算程序输入输出,构造数据集;
S2、根据物理特征设计输入输出变换函数,使之适用于深度神经网络模型训练;
S3、设计深度神经网络算法,执行模型训练,得到预测模型;
S4、对所述的预测模型性能进行评估及评价。
2.根据权利要求1所述的一种托卡马克新经典环向粘滞力矩的机器学习方法,其特征在于,在所述的步骤S1中,获取传统托卡马克新经典环向粘滞力矩计算程序输入输出,构造数据集的具体步骤如下所示:
从实验数据库中选取若干具有不同边界安全因子值的若干个放电实验,并获取所述的放电实验中的等离子体平衡位形、安全因子径向剖面分布、电子密度径向剖面分布、电子温度径向剖面分布、离子温度径向剖面分布和等离子体环向旋转频率径向剖面分布;
基于电子密度、电子温度、离子温度、等离子体环向旋转频率4种物理量的径向剖面分布,通过同时改变幅度和剖面梯度的方式进行剖面增广,每种物理量增广后得到5个剖面分布;
对增广后的物理量剖面分布进行组合,共计得到14×54=8750组等离子体参数位形;
针对上述8750组等离子体参数位形的每一种组合,利用传统新经典环向粘滞力矩数值计算方法,执行新经典环向粘滞力矩计算;
传统新经典环向粘滞力矩数值计算方法步骤为,从线性化漂移动力学方程出发:
;
其中f 0为平衡分布函数,为抛射角散射碰撞算符,/>为粒子速度矢量,(ρ,α,θ)为空间磁面坐标系,v d为偏向频率,f为扰动分布函数,即为本方程的解;
进行坐标变换,从坐标系变换到/>坐标系,其中/>为抛射角变量,x为粒子归一化能量,新坐标定义为:
;
;
;
;
其中表示反弹平均,Bθ表示极向磁场,/>,/>,BM与Bm分别表示磁面上最大和最小磁场强度;
在坐标系中对扰动分布函数和线性化漂移动力学方程进行傅里叶分解,得到:
;
;
其中一些符号定义为:
;
;
;
;
;
;
T为温度,e为电荷,为极向磁通的径向导数,/>为扰动磁场的傅里叶分量,qωE为电漂移频率,ωB为磁漂移频率,DB包含磁漂移频率的极向依赖关系,ωb为粒子反弹频率;求解得到扰动分布函数/>后,即可通过如下相空间积分过程得到NTV力矩径向剖面分布,即:
;
;
;
;
;
;
;
;
;
;
N j 、T j 、e j 分别为第j种粒子的密度、温度、电荷,P为等离子体压强;经过上述计算步骤,能够得到8750组数据样本,对其按8:2的比例随机划分为训练集和测试集。
3.根据权利要求1所述的一种托卡马克新经典环向粘滞力矩的机器学习方法,在所述的步骤S2中,根据物理特征设计输入输出变换函数,使之适用于深度神经网络模型训练,具体步骤包括:
对每种物理量,使用其典型取值进行无量纲化处理,即:
输入向量中的电子密度=原始电子密度/1019;
输入向量中的电子温度=原始电子温度/103;
输入向量中的离子温度=原始离子温度/103;
输入向量中的等离子体旋转频率=原始等离子体旋转频率/104;
输入向量中的安全因子=原始安全因子/10;
由此能够保证输入物理量取值范围均在1-10量级;
托卡马克典型NTV力矩取值范围约为10-5-1N/m2不等,与输入物理量类似,NTV力矩取值范围同样跨越多个数量级尺度,因此,采用如下变换函数对NTV力矩进行变换:
;
其中T NTV 为变换前的NTV力矩,T' NTV 为变换后的NTV力矩,变换后的NTV力矩取值范围约为1-10的量级。
4.根据权利要求1所述的一种托卡马克新经典环向粘滞力矩的机器学习方法,其特征在于,在所述的步骤S4中,计深度神经网络算法,执行模型训练,得到预测模型,具体包括:
将托卡马克背景磁场B0、变换后的电子密度N e 、电子温度T e 、离子温度T i 、旋转频率、安全因子q组成一个输入向量,按照典型的输入物理量分辨率,输入向量长度为1006;
将变换后的电子NTV力矩、离子NTV力矩组成一个输出向量,按照典型的输出物理量分辨率,输出向量长度为130;
所述深度神经网络包含一个输入层,三个隐藏层,一个输出层,各层之间为全连接,其中输入层包含1006个特征,输出层包含130个特征,三个隐藏层分别包含500、520、260个特征;
设置学习率为0.001,采用Adam优化器,设置epochs为50000;
训练结束得到新经典环向粘滞力矩快速预测模型。
5.根据权利要求1所述的一种托卡马克新经典环向粘滞力矩的机器学习方法,其特征在于,在所述的步骤S4中,模型性能评估与评价,具体包括:
观察损失函数随训练过程的变化;
在测试集上计算R2作为模型性能评估指标,R2定义为:
;
其中y (i)表示第i个样本真实值,y (i)表示第i个样本预测值,表示y (i)的平均值,R2越接近1,模型性能越好;
随机选取测试集中的某些样本,直观比较NTV力矩预测值与真实值的差异;
比较传统NTV力矩模拟方法与本方法计算时间差异。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311674834.6A CN117371299B (zh) | 2023-12-08 | 2023-12-08 | 一种托卡马克新经典环向粘滞力矩的机器学习方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311674834.6A CN117371299B (zh) | 2023-12-08 | 2023-12-08 | 一种托卡马克新经典环向粘滞力矩的机器学习方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117371299A true CN117371299A (zh) | 2024-01-09 |
CN117371299B CN117371299B (zh) | 2024-02-27 |
Family
ID=89396912
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311674834.6A Active CN117371299B (zh) | 2023-12-08 | 2023-12-08 | 一种托卡马克新经典环向粘滞力矩的机器学习方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117371299B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118502964A (zh) * | 2024-07-12 | 2024-08-16 | 安徽大学 | 托卡马克新经典环向粘滞力矩cuda模拟实现方法 |
Citations (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190258904A1 (en) * | 2018-02-18 | 2019-08-22 | Sas Institute Inc. | Analytic system for machine learning prediction model selection |
CN110232205A (zh) * | 2019-04-28 | 2019-09-13 | 大连理工大学 | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 |
US20200074325A1 (en) * | 2018-08-29 | 2020-03-05 | National Chiao Tung University | Systems and Methods for Creating an Optimal Prediction Model and Obtaining Optimal Prediction Results Based on Machine Learning |
US20200184131A1 (en) * | 2018-06-27 | 2020-06-11 | Dalian University Of Technology | A method for prediction of key performance parameter of an aero-engine transition state acceleration process based on space reconstruction |
CN112116153A (zh) * | 2020-09-18 | 2020-12-22 | 上海电力大学 | 一种耦合Copula与堆叠式LSTM网络的园区多元负荷联合预测方法 |
WO2021000362A1 (zh) * | 2019-07-04 | 2021-01-07 | 浙江大学 | 一种基于深度神经网络模型的地址信息特征抽取方法 |
US20210042454A1 (en) * | 2020-04-24 | 2021-02-11 | Harbin Institute Of Technology | Linearization Identification Method for Hysteresis Model of Piezoceramics Based on Koopman Operators |
WO2021082809A1 (zh) * | 2019-10-29 | 2021-05-06 | 山东科技大学 | 一种外汇时间序列预测的训练优化方法 |
CN113011073A (zh) * | 2021-04-28 | 2021-06-22 | 清华大学 | 基于深度学习的一维复杂滞回关系构建与结构模拟方法 |
WO2021123622A1 (fr) * | 2019-12-19 | 2021-06-24 | Extrality | Méthode de simulation numérique par apprentissage automatique |
CN113408852A (zh) * | 2021-05-18 | 2021-09-17 | 江西师范大学 | 基于在线学习行为和深度神经网络的元认知能力评估模型 |
US20210350218A1 (en) * | 2020-05-06 | 2021-11-11 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Bayesian neural network with resistive memory hardware accelerator and method for programming the same |
US20210357555A1 (en) * | 2018-09-14 | 2021-11-18 | Northwestern University | Data-driven representation and clustering discretization method and system for design optimization and/or performance prediction of material systems and applications of same |
CN113780522A (zh) * | 2021-08-27 | 2021-12-10 | 核工业西南物理研究院 | 基于深度神经网络的托卡马克等离子体大破裂预测算法 |
CN113792620A (zh) * | 2021-08-27 | 2021-12-14 | 核工业西南物理研究院 | 基于深度神经网络的托卡马克边缘局域模实时识别算法 |
US20210390230A1 (en) * | 2020-06-16 | 2021-12-16 | Chongqing University | Method for Quickly Optimizing Key Mining Parameters of Outburst Coal Seam |
CN114118361A (zh) * | 2021-10-29 | 2022-03-01 | 北京宇航系统工程研究所 | 一种基于深度学习参数锚定的态势评估方法 |
WO2022151890A1 (zh) * | 2021-01-13 | 2022-07-21 | 中国长江三峡集团有限公司 | 一种风电场发电量评估及微观选址模型建立方法 |
CN114996631A (zh) * | 2022-05-20 | 2022-09-02 | 安徽大学 | 一种托卡马克等离子体平衡位形轻量级重建方法 |
WO2022212313A1 (en) * | 2021-03-30 | 2022-10-06 | Siemens Corporation | Ml-based discrepancy surrogate for digital twin simulation |
CN115204227A (zh) * | 2022-07-12 | 2022-10-18 | 北京航空航天大学 | 基于深度学习的设备故障诊断中不确定性量化校准方法 |
CN115239018A (zh) * | 2022-08-15 | 2022-10-25 | 华南农业大学 | 基于物理信息神经网络的植保无人机下洗风场重建方法 |
CA3159935A1 (en) * | 2021-05-20 | 2022-11-20 | Royal Bank Of Canada | System and method for adversarial vulnerability testing of machine learning models |
WO2022247049A1 (zh) * | 2021-05-24 | 2022-12-01 | 苏州大学 | 基于复数值前向神经网络的风速预测方法 |
CN115563604A (zh) * | 2022-10-19 | 2023-01-03 | 四川大学 | 基于深度神经网络和特征融合的口令强度评估方法及系统 |
CN115615432A (zh) * | 2022-10-11 | 2023-01-17 | 常州大学 | 一种基于深度神经网络的室内行人惯性导航方法 |
DE102022210046A1 (de) * | 2021-09-24 | 2023-03-30 | Robert Bosch Gesellschaft mit beschränkter Haftung | Rechnergestützte neuronales-netz-kraftfeld-trainingsroutinen für molekulardynamik-computersimulationen |
WO2023106990A1 (en) * | 2021-12-10 | 2023-06-15 | Kaaberg Johard Leonard | A modular, variable time-step simulator for use in process simulation, evaluation, adaption and/or control |
CN116415730A (zh) * | 2023-04-12 | 2023-07-11 | 安徽大学 | 一种预测水位的融合自注意力机制时空深度学习模型 |
US20230268084A1 (en) * | 2021-03-15 | 2023-08-24 | Dalian University Of Technology | Simulation method for electron temperature evolution caused by east tokamak radiofrequency wave |
WO2023155231A1 (zh) * | 2022-02-21 | 2023-08-24 | 东南大学 | 一种高度类人的自动驾驶营运车辆安全驾驶决策方法 |
CN116700002A (zh) * | 2023-06-28 | 2023-09-05 | 北京航空航天大学 | 基于深度神经网络的压电执行器Hammerstein建模和预测控制方法 |
CN116720158A (zh) * | 2023-04-26 | 2023-09-08 | 中国科学院合肥物质科学研究院 | 一种带不确定度估计的时间序列回归预测方法及系统 |
-
2023
- 2023-12-08 CN CN202311674834.6A patent/CN117371299B/zh active Active
Patent Citations (33)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190258904A1 (en) * | 2018-02-18 | 2019-08-22 | Sas Institute Inc. | Analytic system for machine learning prediction model selection |
US20200184131A1 (en) * | 2018-06-27 | 2020-06-11 | Dalian University Of Technology | A method for prediction of key performance parameter of an aero-engine transition state acceleration process based on space reconstruction |
US20200074325A1 (en) * | 2018-08-29 | 2020-03-05 | National Chiao Tung University | Systems and Methods for Creating an Optimal Prediction Model and Obtaining Optimal Prediction Results Based on Machine Learning |
US20210357555A1 (en) * | 2018-09-14 | 2021-11-18 | Northwestern University | Data-driven representation and clustering discretization method and system for design optimization and/or performance prediction of material systems and applications of same |
CN110232205A (zh) * | 2019-04-28 | 2019-09-13 | 大连理工大学 | 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法 |
WO2021000362A1 (zh) * | 2019-07-04 | 2021-01-07 | 浙江大学 | 一种基于深度神经网络模型的地址信息特征抽取方法 |
WO2021082809A1 (zh) * | 2019-10-29 | 2021-05-06 | 山东科技大学 | 一种外汇时间序列预测的训练优化方法 |
WO2021123622A1 (fr) * | 2019-12-19 | 2021-06-24 | Extrality | Méthode de simulation numérique par apprentissage automatique |
US20210042454A1 (en) * | 2020-04-24 | 2021-02-11 | Harbin Institute Of Technology | Linearization Identification Method for Hysteresis Model of Piezoceramics Based on Koopman Operators |
US20210350218A1 (en) * | 2020-05-06 | 2021-11-11 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Bayesian neural network with resistive memory hardware accelerator and method for programming the same |
US20210390230A1 (en) * | 2020-06-16 | 2021-12-16 | Chongqing University | Method for Quickly Optimizing Key Mining Parameters of Outburst Coal Seam |
CN112116153A (zh) * | 2020-09-18 | 2020-12-22 | 上海电力大学 | 一种耦合Copula与堆叠式LSTM网络的园区多元负荷联合预测方法 |
WO2022151890A1 (zh) * | 2021-01-13 | 2022-07-21 | 中国长江三峡集团有限公司 | 一种风电场发电量评估及微观选址模型建立方法 |
US20230268084A1 (en) * | 2021-03-15 | 2023-08-24 | Dalian University Of Technology | Simulation method for electron temperature evolution caused by east tokamak radiofrequency wave |
WO2022212313A1 (en) * | 2021-03-30 | 2022-10-06 | Siemens Corporation | Ml-based discrepancy surrogate for digital twin simulation |
CN113011073A (zh) * | 2021-04-28 | 2021-06-22 | 清华大学 | 基于深度学习的一维复杂滞回关系构建与结构模拟方法 |
CN113408852A (zh) * | 2021-05-18 | 2021-09-17 | 江西师范大学 | 基于在线学习行为和深度神经网络的元认知能力评估模型 |
CA3159935A1 (en) * | 2021-05-20 | 2022-11-20 | Royal Bank Of Canada | System and method for adversarial vulnerability testing of machine learning models |
WO2022247049A1 (zh) * | 2021-05-24 | 2022-12-01 | 苏州大学 | 基于复数值前向神经网络的风速预测方法 |
CN113780522A (zh) * | 2021-08-27 | 2021-12-10 | 核工业西南物理研究院 | 基于深度神经网络的托卡马克等离子体大破裂预测算法 |
CN113792620A (zh) * | 2021-08-27 | 2021-12-14 | 核工业西南物理研究院 | 基于深度神经网络的托卡马克边缘局域模实时识别算法 |
DE102022210046A1 (de) * | 2021-09-24 | 2023-03-30 | Robert Bosch Gesellschaft mit beschränkter Haftung | Rechnergestützte neuronales-netz-kraftfeld-trainingsroutinen für molekulardynamik-computersimulationen |
CN114118361A (zh) * | 2021-10-29 | 2022-03-01 | 北京宇航系统工程研究所 | 一种基于深度学习参数锚定的态势评估方法 |
WO2023106990A1 (en) * | 2021-12-10 | 2023-06-15 | Kaaberg Johard Leonard | A modular, variable time-step simulator for use in process simulation, evaluation, adaption and/or control |
WO2023155231A1 (zh) * | 2022-02-21 | 2023-08-24 | 东南大学 | 一种高度类人的自动驾驶营运车辆安全驾驶决策方法 |
CN114996631A (zh) * | 2022-05-20 | 2022-09-02 | 安徽大学 | 一种托卡马克等离子体平衡位形轻量级重建方法 |
CN115204227A (zh) * | 2022-07-12 | 2022-10-18 | 北京航空航天大学 | 基于深度学习的设备故障诊断中不确定性量化校准方法 |
CN115239018A (zh) * | 2022-08-15 | 2022-10-25 | 华南农业大学 | 基于物理信息神经网络的植保无人机下洗风场重建方法 |
CN115615432A (zh) * | 2022-10-11 | 2023-01-17 | 常州大学 | 一种基于深度神经网络的室内行人惯性导航方法 |
CN115563604A (zh) * | 2022-10-19 | 2023-01-03 | 四川大学 | 基于深度神经网络和特征融合的口令强度评估方法及系统 |
CN116415730A (zh) * | 2023-04-12 | 2023-07-11 | 安徽大学 | 一种预测水位的融合自注意力机制时空深度学习模型 |
CN116720158A (zh) * | 2023-04-26 | 2023-09-08 | 中国科学院合肥物质科学研究院 | 一种带不确定度估计的时间序列回归预测方法及系统 |
CN116700002A (zh) * | 2023-06-28 | 2023-09-05 | 北京航空航天大学 | 基于深度神经网络的压电执行器Hammerstein建模和预测控制方法 |
Non-Patent Citations (5)
Title |
---|
XINGTING YAN ET AL.: "Neoclassical toroidal viscosity torque in tokamak edge pedestal induced by external resonant magnetic perturbation", 《PHYSICS OF PLASMAS》 * |
闫星廷, 朱平, 孙有文: "Unified modeling of both resonant and non-resonant neoclassical transport under non-axisymmetric magnetic perturbations in tokamaks", 《PHYSICS OF PLASMAS》, vol. 26, no. 7 * |
闫星廷, 朱平, 孙有文: "共振磁扰动在边界台基区引起的新经典环向粘滞力矩", 《第十八届全国等离子体科学技术会议》 * |
闫星廷: "托卡马克中新经典环向粘滞力矩及其作用下等离子体稳态转动的模拟计算研究.", 《 《中国博士学位论文全文数据库(工程科技Ⅱ辑)》 * |
陈峰;谷俊涛;李玉磊;彭晓溪;韩天甲;: "基于机器视觉和卷积神经网络的东北寒地玉米害虫识别方法", 江苏农业科学, no. 18 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118502964A (zh) * | 2024-07-12 | 2024-08-16 | 安徽大学 | 托卡马克新经典环向粘滞力矩cuda模拟实现方法 |
CN118502964B (zh) * | 2024-07-12 | 2024-09-24 | 安徽大学 | 托卡马克新经典环向粘滞力矩cuda模拟实现方法 |
Also Published As
Publication number | Publication date |
---|---|
CN117371299B (zh) | 2024-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117371299B (zh) | 一种托卡马克新经典环向粘滞力矩的机器学习方法 | |
Alizadeh et al. | Ensemble of surrogates and cross-validation for rapid and accurate predictions using small data sets | |
CN107591815B (zh) | 一种求解含离散控制电力系统无功优化的分解方法 | |
Möller et al. | A method to calculate fission-fragment yields Y (Z, N) versus proton and neutron number in the Brownian shape-motion model: Application to calculations of U and Pu charge yields | |
CN107766656A (zh) | 基于历史数据回归分析的汽轮机阀门流量特性分析方法 | |
CN101968832B (zh) | 基于构造-剪枝混合优化rbf网络的煤灰熔点预测方法 | |
CN112817312A (zh) | 一种基于双重搜索优化算法的路径规划方法 | |
CN107957565A (zh) | 一种核磁共振波谱仪自屏蔽梯度线圈及其设计方法 | |
CN111310348A (zh) | 一种基于pso-lssvm的材料本构模型预测方法 | |
CN101788789A (zh) | 基于混沌混合优化算法的单元机组非线性预测控制方法 | |
CN104750948A (zh) | 一种处理飞行器设计中多极值多约束问题的优化方法 | |
Van Kien et al. | Adaptive inverse multilayer fuzzy control for uncertain nonlinear system optimizing with differential evolution algorithm | |
Ma et al. | Computational intelligence nonmodel-based calibration approach for internal combustion engines | |
CN115169470A (zh) | 一种基于可接受区域的高维小样本数据扩充方法 | |
CN112966399B (zh) | 一种基于机器学习的脉管制冷机工况预测方法及系统 | |
Janiš et al. | Kondo behavior in the asymmetric Anderson model: Analytic approach | |
CN113722853B (zh) | 一种面向智能计算的群智能进化式工程设计约束优化方法 | |
CN111522226B (zh) | 针对伺服转台的多目标优化高型pid最优控制器设计方法 | |
CN106126926B (zh) | 一种能够处理温度分布效应的子群共振计算方法 | |
CN112182850B (zh) | 直线电机模型预测控制权重因子优化方法 | |
Yin et al. | Improved Hybrid Fireworks Algorithm‐Based Parameter Optimization in High‐Order Sliding Mode Control of Hypersonic Vehicles | |
Zhang et al. | The amount prediction and optimization of the returned ore generated from sintering process based on SHAP value and ensemble learning | |
CN109711040B (zh) | 一种基于搜索方向学习的智能工业设计强化学习算法 | |
Kallakuri et al. | Closed loop modeling of the APS-U orbit feedback system | |
Chu et al. | Novel robust approach for constructing Mamdani-type fuzzy system based on PRM and subtractive clustering algorithm |
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 |