基于贝叶斯优化的高能质子空间传输精确控制.pdf
《基于贝叶斯优化的高能质子空间传输精确控制.pdf》由会员分享,可在线阅读,更多相关《基于贝叶斯优化的高能质子空间传输精确控制.pdf(7页珍藏版)》请在咨信网上搜索。
1、粒子束及加速器技术基于贝叶斯优化的高能质子空间传输精确控制*申诗雨1,杨晓虎1,张国博1,赵子琦1,马燕云2(1.国防科技大学理学院核科学与技术系,长沙410073;2.国防科技大学前沿交叉学科学院,长沙410073)摘要:综合考虑地磁场、高能质子的相对论效应以及同步辐射的影响,建立了质子在空间传输的单粒子运动模型。基于该模型,提出利用贝叶斯优化方法,实现给定质子能量下,质子从空间初始点传输到目标点的精确控制方法,获得了出射方向随出射位置的变化规律,当位置径向角取 0和 180时,位置轴向角的取值不会改变粒子的最优出射方向。研究结果为质子束在空间环境中的长程传输提供理论支撑。关 键 词:高能质
2、子;传输控制;空间传输;贝叶斯优化中图分类号:O53文献标志码:Adoi:10.11884/HPLPB202335.230231Precisecontrolofhigh-energyprotonstransportinspaceenvironmentbyusingbayesianoptimizationShenShiyu1,YangXiaohu1,ZhangGuobo1,ZhaoZiqi1,MaYanyun2(1.Department of Nuclear Science and Technology,National University of Defense Technology,Chan
3、gsha 410073,China;2.College of Advanced Interdisciplinary Studies,National University of Defense Technology,Changsha 410073,China)Abstract:Consideringthegeomagneticfield,therelativisticeffectandbremsstrahlungradiationofhigh-energyprotons,asingleparticlemotionmodelofprotontransportinthespaceenvironme
4、ntisestablished.Basedonthismodel,theBayesianoptimizationmethodisproposedtorealizetheprecisecontrolofprotonstransportfromtheinitialpositiontothetargetunderagivenprotonenergy.Thedependenceoftheprotonlaunchangleonthelaunchheightisobtained,thatis,whenthecoordinateradialangleis0and180,thevalueofthecoordi
5、nateaxialanglewillnotchangetheoptimalemissiondirectionoftheparticles.Theresultscanprovidetheoreticalreferencesforlong-distancetransportofprotonbeamsinthespaceenvironment.Keywords:highenergyproton,transportcontrol,spacetransmission,Bayesianoptimization在目前的空间主动实验中,已经完成的形式包括发射带电粒子束、释放等离子体、释放化学物质和核爆炸等。其
6、中发射带电粒子对于辐射防护以及空间碎片清除等方面具有工程价值。由于探寻带电粒子在空间传输中的精确控制实验成本昂贵,因此有必要利用理论方法、数值模拟等简易、经济的研究手段,对质子束在空间环境中的长程传输控制开展预研和探索。由于高能电子束加速相比起高能质子束加速要相对容易,在国内外开展的研究中主要集中在电子束传输研究,对高能质子的研究相对较少。Neubert 等1-2使用蒙特卡罗模拟方法,分析了同向地磁场对束包络的影响;Mironychev 等3则结合蒙特卡罗方法计算了电子束水平射入大气环境的情况,指出垂直方向的地磁场会改变电子束的传输方向;Porazik 等4-5研究了电子束损失锥的计算以及电子
7、束注入位置的选择,之后 Powis6利用单粒子运动模型模拟了电子束从轨道卫星发射后,沿磁力线传输并最终抵达地球附近的全过程,发现磁层环境的变化会显著影响电子束在大气层内束斑的落点。Hao 等7-8利用二维 PIC 方法,建立了适用于真空环境内相对论电子束长程传输的计算模型,总结了地磁场及自生电磁场对电子束长程传输的影响。姚海波等人的研究表明通过调制相对论电子束的束/背景等离子体密度比值、相对论因子、束包络和密度分布等,可以抑制电子束传输中的束流不稳定*收稿日期:2023-07-26;修订日期:2023-09-19基金项目:国家自然科学基金项目(12175309、11975308、1200529
8、7);中国科学院战略先导专项(XDA25050200);湖南省自然科学基金项目(CX20190001);国防科技大学青年创新奖资助课题(20180104)联系方式:申诗雨,。通信作者:杨晓虎,。第35卷第10期强激光与粒子束Vol.35,No.102023年10月HIGHPOWERLASERANDPARTICLEBEAMSOct.,2023104005-1性,如成丝不稳定性和水龙管不稳定性,保证电子束稳定的长距离传输9。然而,人们对高能质子的研究较少,而高能质子的传输对空间武器和空间垃圾清理等的研究有重要作用,故有必要对高能质子的空间传输精确控制进行研究。高能质子在空间中的传输与地磁场的强度及
9、方向、本身能量等因素密切相关,已有的研究并未充分考虑地磁场对于高能质子的偏转,磁场中运动粒子的同步辐射,高能质子的相对论效应对传输过程的影响,以及高能质子传输的准确控制等。本文细致研究了不同能量的高能质子在空间传输的精确控制问题,综合考虑地磁场、相对论以及同步辐射等影响,建立了高能质子的单粒子运动模型。给定初始能量,利用贝叶斯优化方法,实现了高能质子在空间中从起始点到目标点的精确传输。1高能质子在地磁场中的传输首先,我们建立地磁场的笛卡尔坐标系,认为地磁北极与地理北极重合,以地球球心为圆心,以从球心向地磁北极的连线为 z 轴方向,向赤道与零度经线交点的连线为 x 轴方向建立右手笛卡尔坐标系。高
10、能质子在空间中传输时受到地磁场的作用,由于受到洛伦兹力而产生速度的偏转,在某种情况下地球可认为是一个巨大的磁偶极子。所以,可以用偶极场来近似地磁场,其模型为10Bdip(r)=04r33(M r)r M(1)r=x x+y y+z z r=r/rMM=M zx=ReRey=z=0B0=3.07105T0M/4=B0R3e式中:为矢径,为单位矢量,为地球的磁矩。对于地球,反平行于 z 轴。在地球赤道处(为地球半径),地磁场强度约为,得。于是在笛卡尔坐标中,磁场的表达式为Bdip=B0R3er53xz x+3yz y+(2z2x2y2)z(2)R2=0.973 9选取零度经线与赤道交点处高度为 0
11、km 至 1000km每隔 10km 取一个值的数据与 IGRF 模型数据比较,如图 1所示,可以看出理论模型计算所得数据与 IGRF 数据库11数据拟合效果较好,拟合优度计算得,验证较好,可以代替实际地磁场使用。m0近地空间中,带电量为 q、质量为、速度为 v 的质子,在电场强度为 E、磁感应强度为 B 的场中运动时,其所受的磁场作用力由洛伦兹力公式给出,即dm0vdt=qE(r)+qvB(r)(3)1/10当粒子能量大于本身静止能量的时,应当考虑相对论对于粒子速度、质量等的影响,相对论对于粒子的影响根据洛伦兹变换公式给出,即dmvdt=qE(r)+qvB(r)(4)=(1v2/c2)m=m
12、0/(1v2/c2)式中:为相对论因子,。给定质子能量为 100MeV,出射位置为高度 500km、经度 0、北纬 30处,出射角度为(60,0),即出射径向角为 60,出射轴向角为 0,运动时间为 0.01s。是否考虑相对论的高能质子在地磁场中运动轨迹图如图 2 所示。可以发现,质子在地磁场中呈回旋运动,总体轨迹呈沿磁场方向螺旋前进的螺旋线,而由于相对论的作用,质子质量发生变化,变化约为 11%,最终作用于洛伦兹力偏转半径时相差量级约为 100m,虽然运动时间很短,但由于质子速度接近光速,导致实际运动距离较长,最终导致差距较大。如图 2 所示,在运动时间为 0.01s 时考虑相对论效应和未考
13、虑相对论效应对于质子运动有较大影响,经过计算得同一时刻下两种情况间最大距离可达 20km 以上。2同步辐射对高能质子传输的影响磁场中的高能粒子运动时受到洛伦兹力的作用产生垂直于运动方向的加速度,从而发生同步辐射,粒子在发02004006008001 000h/km2.02.53.03.5B/(105 T)theoretical model dataIGRF database dataFig.1MagneticinductionintensitydiagramoftheoreticalmodelandIGRFdatabase图1理论模型与 IGRF 数据库磁感应强度图强激光与粒子束104005-
14、2射同步辐射后损失能量,影响后期洛伦兹偏转半径和运动轨迹,故需予以考虑。粒子在磁场中发射的同步辐射功率由下式确定12P(t)=e2 v260c31(1v2/c2)2(5)给定质子能量为 100MeV,出射位置为高度 500km、经度 0,北纬 30处,出射方向为(60,0),即出射径向角为 60,出射轴向角为 0,运动时间为 0.01s。是否考虑同步辐射的高能质子在地磁场中运动轨迹图如图 3 所示。可以看出由于地磁场的磁镜效应,质子出现弹跳运动;同步辐射对质子运动轨迹的影响与相对论的影响相似,是否考虑同步辐射对于质子的总体运动轨迹影响不大,其对于质子的影响在运动初期并不明显,但由于质子的速度接
15、近光速,仅 0.01s 的时间两者之间的距离已经相差了 80km 以上,相差了将近一个螺旋,可以说明同步辐射的影响会随着运动距离的增加而逐渐增加,予以考虑。3空间中的单粒子运动模型建立及应用3.1空间中的单粒子运动模型的建立由于粒子在运动过程中同时受到地磁场偏转、相对论效应以及同步辐射的影响,因此在综合计算中将时间离散化后需要在一个时间步长里同时进行三个效果的计算。由于时间步长取得很小,在某一个时间步长内的计算可认为粒子的状态不发生变化。首先使用上一时刻的最终速度计算由相对论效应引发的质量变化;之后由于地磁场的作用发生洛伦兹偏转及同步辐射:受到洛伦兹偏转会影响粒子的加速度,受到同步辐射会影响粒
16、子能量进而影响速度,且偏转和辐射应是同步发生,故选择上一时刻的速度及本时刻计算所得的质量同时计算偏转和辐射,在计算的过程中需将粒子的运动在三维进行展开来分别计算每个方向的速度及加速度。再使用通过偏转计算得到的加速度以及通过辐射计算得到的速度来计算并保存本时刻的位移、计算并保存本时刻的最终速度并作为下一时刻的初始速度。3.2空间中的单粒子运动模型的应用X0Xtdt=108s建立笛卡尔坐标系,并通过坐标变换将给定坐标中的经纬度换算为笛卡尔坐标系下的坐标,给定初始点的坐标、目标点的坐标、质子的初始能量及出射角度,将质子的运动在时间上进行离散化,以为一个时间步长,共走 N 步,应用上述空间中的单粒子运
17、动模型,计算每个时刻粒子的坐标及速度并保存,最后得到粒子运动坐标(Xp)及运动速度(vp)的 3*N 维矩阵。XpXtN通过计算中每个坐标点与目标点的距离,找到其中的最小值,得到此出射速度方向下与目标点的最小距离;通过使用每时刻粒子的速度即可计算得到每时刻粒子的能量,与初始能量相比即可得到每时刻的能量损耗;再找到距离最小值的坐标后,其对应的步数即对应粒子运动的时间。4基于贝叶斯优化方法对高能质子传输的优化控制高能质子在空间传输中由于受到地磁场影响,其运动轨迹比较复杂,难以实现准确控制。因此,我们提出利用500y/km3 4203 4403 460503 480 x/km3 5005 4003
18、520z/km5 5003 5405 6003 5605 7003 5805 8003 6003 6205 9006 000without relativistic effectwith relativistic effectstarting pointfinal point without relativistic effectfinal point with relativistic effectFig.2Protontrajectorywithandwithoutrelativisticeffect图2考虑相对论效应和未考虑相对论效应的运动轨迹图0200y/km5 0006 0007 0
19、002 0008 000 x/km4009 00010 0000z/km2 0004 000without synchrotron radiationwith synchrotron radiationstarting pointFig.3Protontrajectorywithandwithoutsynchrotronradiation图3考虑同步辐射和未考虑同步辐射的质子运动轨迹图申诗雨,等:基于贝叶斯优化的高能质子空间传输精确控制104005-3P(A|B)=P(B|A)*P(A)/P(B)P(A|B)P(B|A)P(A)P(B)贝叶斯优化方法对高能质子空间传输进行精确控制。贝叶斯优化是
20、一种非常有效的全局优化算法,可以快速高效的找到相对最优的输入参数。目标函数通常是一个黑盒函数,它的精确表达式或导数是未知。高斯过程是贝叶斯优化的基础13。一个高斯过程由其均值函数和协方差函数唯一定义14。高斯过程假设给定参数为正态分布,并将具有不确定性的未知函数建模为替代函数。通过应用贝叶斯公式,可以从先验分布得到后验分布。表示为后验分布,代表需要优化的分布。同时,代表目标函数,是未知函数的先验分布,代表证据或边际可能数据15。经过多次迭代,基于高斯过程的贝叶斯优化器会以一定的置信水平逼近函数,并逐渐增强置信水平,收敛到最优解。y=f(X)以到目标点距离和质子出射角度的优化为例。将目标点坐标,
21、出射角度及其对应的路径中每个点的坐标之间的关系建模为,其中 X 代表相关参数X=X1,Xn(6)y=(XtXp)2设为目标点与路径点之间的距离,可以将问题转化为寻找出射角度对应到目标点最小距离的优化问题 A,即X=argminXAf(X)(7)X0Xf(X)贝叶斯优化需要初始采样来获得数据集。给定均值函数和核函数,高斯过程回归然后可以通过使用均值向量和协方差矩阵来提供数据集的简洁表示。X0(X)使用均值和核函数获得初始函数 f1,但由于有限的数据和高度不确定性,它可能不是最佳的。为了改进估计,除了先前的数据点,还引入了额外的观察数据点。这降低了不确定性,并提供了更好的设计变量。可以使用各种方法
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 贝叶斯 优化 高能 质子 空间 传输 精确 控制
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。