太阳风磁帆推进性能数值研究_蔡静媛.pdf
《太阳风磁帆推进性能数值研究_蔡静媛.pdf》由会员分享,可在线阅读,更多相关《太阳风磁帆推进性能数值研究_蔡静媛.pdf(10页珍藏版)》请在咨信网上搜索。
1、利用考虑行星际磁场作用的磁流体动力学模型,建立了磁帆三维数值模拟方法,对计算方法的可靠性进行了验证,发现了线圈尾部的磁重联现象,研究了太阳风来流速度、等离子体离子数密度以及攻角对磁帆推进性能的影响。得出以下结论:不同速度、不同离子数密度的太阳风主要通过改变 z 方向电流的大小改变洛伦兹力,进而影响磁帆的推进性能:太阳风离子数密度恒定时,随着来流速度由 30km/s 逐渐增大至75km/s,z 方向电流最大值由 4205A/m2增至 14709A/m2,磁帆所受推力由 3.39N 增至 13.40N;太阳风来流速度恒定时,随着离子数密度由 1.81019m3增大至 4.51019m3,z 方向电
2、流最大值由 6039A/m2增至 10585A/m2,磁帆所受推力由 6.62N 增至 12.27N。磁帆攻角变化,主要通过磁场构型的变化影响磁帆推进性能:攻角为 0和 90时的磁层半径分别为 0.14m 和 0.18m,磁帆所受推力分别为 6.62N 和 11.03N,由此推测实际应用中保持线圈轴线与太阳风来流方向平行,可获得更大推力。系统研究了相关因素对磁帆推进性能的影响,可为磁帆的推力调节研究提供参考和支持,对未来磁帆的深入研究具有重要的参考价值。关键词:磁帆;磁流体动力学(MHD)推进;太阳风速度;离子数密度;攻角中图分类号:V439.4文献标志码:ANumericalstudyonp
3、ropulsionperformanceofsolarwindmagneticsailCAIJingyuan,LILai,ZHUGuiping(CollegeofAstronautics,NanjingUniversityofAeronauticsandAstronautics,Nanjing210016,China)Abstract:The three-dimensional numerical simulation of the solar wind magnetic sail wasestablishedwiththemagnetohydrodynamic(MHD)modelconsid
4、eringtheinterplanetarymagneticfield.Theverificationofthecalculationmethodwasaccomplishedbycomparisonwithexperimentaldata.Inaddition,theobservationandconfirmationofthemagneticreconnectionwereachievedatthetailofthecoil.Thepropulsionperformanceofthemagneticsailwasstudiedintermsoftheincomingvelocity,the
5、plasma ion number density and the attack angle of of the solar wind.The solar wind with differentvelocitiesanddifferentionnumberdensitiesmainlyinfluencedtheLorentzforcebyaffectingthecurrentinthezdirection,whichfurtheraffectedthepropulsionperformanceofthemagneticsail.Asincomingvelocityincreasedfrom30
6、km/sto75km/sofsolarwindwithfixedionnumberdensity,themaximumcurrentinzdirectionincreasedfrom4205A/m2to14709A/m2,andthethrustofmagneticsailincreasedfrom3.39Nto13.40N.Astheionnumberdensityincreasedfrom1.81019m3to4.51019m3ofsolarwindwithfixedincomingvelocity,themaximumcurrentinzdirectionincreasedfrom603
7、9A/m2to10585A/m2,andthethrustincreasedfrom6.62Nto12.27N.Thevariationoftheattackangleaffectedthepropulsionperformanceofmagneticsailbyinfluencingtheconfigurationofmagneticfield.Withtheattack angle of 0 and 90,the radius of the magnetic cavity was 0.14 m and 0.18 m,respectively.收稿日期:2021-09-21作者简介:蔡静媛(
8、1998),女,硕士生,主要从事新型推进技术研究。通信作者:朱桂平(1983),女,副教授,博士,主要从事新型推进技术研究。E-mail:引用格式:蔡静媛,李来,朱桂平.太阳风磁帆推进性能数值研究J.航空动力学报,2023,38(4):850-859.CAIJingyuan,LILai,ZHUGuiping.NumericalstudyonpropulsionperformanceofsolarwindmagneticsailJ.JournalofAerospacePower,2023,38(4):850-859.第38卷第4期航空动力学报Vol.38No.42023年4月JournalofA
9、erospacePowerApr.2023Correspondingly,thethrustofthemagneticsailwas6.62Nand11.03N,respectively.Itwasinferredthatlargerthrustcanbeobtainedbykeepingtheaxisofthecoilparalleltothedirectionofthesolarwindinpracticalapplication.Theinfluenceofrelevantfactorsonthepropulsionperformanceofmagneticsailwasstudieds
10、ystematically,whichcanprovideareferenceandsupportfortheresearchofthrustregulationofthesail,andhasimportantreferencevalueforthefurtherstudyofmagneticsail.Keywords:magneticsails;magnetohydrodynamic(MHD)propulsion;solarwindvelocity;ionnumberdensity;attackangle磁帆技术最早由 Zubrin 和 Andrews 提出1-2,它利用太阳风中粒子与磁线
11、圈生成的磁场相互作用产生推力。太阳风是从恒星上层大气射出的超声速等离子体带电粒子流,以 200800km/s 的速度向外辐射。航天器携带低密度超导材料制作的线圈,发射入轨之后展开线圈并通电,线圈中的电流将在航天器周围形成一个磁场区域称作人工磁腔,通过该磁场与空间中的太阳风等离子场相互作用产生动量交换从而为航天器提供动力。磁帆仅需电源系统即可提供充足能源,是一种无工质损耗、高质量利用率的深空探测动力方式,开展太阳风磁帆推进系统研究对未来空间飞行具有深远意义。自磁帆的概念提出之后,各国学者便逐步开展该项技术的相关理论及实验研究。Nishida 等3运用磁流体动力学(MHD)模型模拟了太阳风和磁帆偶
12、极子磁场之间的相互作用,验证了太阳风和航天器之间的动量传递过程,阐明了磁帆获得推力的原理及方式。Akita 和 Suzuki4采用网格质点法(particle-in-cell,PIC)进行数值计算,得到了流场和磁帆周围电磁场的基本特征,以及磁帆所受作用力。Kajimura 等5-6采用混合网格质点法(hybridparticle-in-cell,HPIC)进行三维仿真,研究了磁帆周围的等离子体流特征,探讨了磁帆的电磁推力特性,如阻力、升力、横向力的值等,研究了攻角和推力大小、磁矩和磁帆姿态之间的关系。Ashida 等7提出了一个包括离子的有限拉莫尔半径效应的数值模型,采用 flux-tube
13、模型求解离子的粒子运动轨迹,并验证了该计算方法的可靠性;运用 PIC 模拟程序对小型磁层开展了三维数值模拟8,研究了小型磁帆的推力特性。Ueno等9对推力大小为 1N 左右的磁帆比例模型进行了实验研究,获得了离子尺度下磁帆的推力特性;并采用平行四边形摆法10对空间中磁帆的推力特性进行了实验测试。上述仅由一个线圈组成的磁帆通常称为纯磁帆。根据 Zubrin 的理论,若要太阳风粒子与磁场相互作用产生推力进行加速,磁线圈的直径必须达到 100km 级以上,但这在实际发射和太空展开中都很难实现。为了克服这个问题,Winglee 等11-13在原有研究的基础上提出了将磁帆与等离子体射流结合起来的微磁层等
14、离子推进技术(minimag-netosphericplasmapropulsion,M2P2)。众多学者对这种自带等离子体的磁帆模型进行了数值与实验研究14-17,但其推进原理与原始磁帆并无二致,且仍需消耗大量的等离子体才能使磁层膨胀到所需大小。虽然经过 30 多年的发展,磁帆推进技术的数值模拟方法已趋于成熟,线圈参数对磁帆推进性能的影响也已经研究得比较透彻,但太阳风变化对其性能的影响机制还尚未明确。太阳风是速度和导电特性多变的等离子体,也是磁帆推进系统主要的动力来源,研究其参数对磁帆推进性能的影响,有助于进一步理解太阳风与磁场相互作用原理、完善磁帆推进性能影响机制。因此本文在数值计算方法得
15、到验证的基础上,系统地研究了太阳风来流速度、等离子体离子数密度以及攻角对磁帆推进性能的影响,可为在不同性质太阳风来流下磁帆的推力调节研究提供参考和支持。1计算方法本研究采用理想磁流体动力学(MHD)控制方程组,规范化后的方程组如下18:t+(V)=0(1)(V)t+(VV+pI)=J B(2)Bt(VB)=0(3)Et+(E+p+B22)VB(BV)=0(4)J其中 为电流密度,由诱导磁场法可得第4期蔡静媛等:太阳风磁帆推进性能数值研究851J=10B(5)E 为能量密度,其表达式为E=12V2+p1+B22(6)pVB5/3式中 是等离子体密度,为压强,为太阳风速度矢量,其在 x、y 和 z
16、 方向上的分量分别记作 u、v和 w。为磁场强度。由于太阳风主要组成为质子和电子,所以这里 设置为单原子理想气体的比热值。当太阳风等离子体流经线圈产生的磁场时,带电粒子受到洛伦兹力作用反射回去,磁帆受太阳风反作用力推动前行。线圈磁场产生的磁压力与外部太阳风来流的动压在磁场某处达到一种平衡即(2Bm)220=nimiV2sw(7)2Bm0nimiVsw式中为平衡处的磁场强度,为真空磁导率,为离子数密度,为质子质量,为太阳风来流速度。L平衡处边界至线圈中心的距离 即磁腔半径,其表达式为L=(0M2d82nimiV2sw)1/6(8)MdS=L2L式中为线圈产生的磁矩大小。平衡处所围成的磁场区域称为
17、磁腔,磁腔的大小。距离线圈中心距离为 处的磁场强度为Bm=0Md4L3(9)riL磁腔的边界称为磁层顶,磁层顶处离子的拉莫半径为riL=miVswe2Bm(10)e式中 为元电荷。由式(8)式(10)可以进一步推导拉莫尔半径的另一表达式riL=mie(20nimi)1/2(11)riLLriL/L推力系数 Cd的值与拉莫尔半径和 的比值有关,其表达式为Cd=|3.4riL/Lexp0.22(riL/L)2riL/L 13.6exp0.28(riL/L)2riL/L 1(12)F磁帆所受推力 的公式为F=Cd12V2swS(13)V2sw/2S式中为太阳风的动压,为上述磁腔面积。2计算模型为了考
18、察算法和程序的准确性,本文将数值模拟结果同文献 9 中的实验结果进行了对比。实验采用内径为 0.05m、外径为 0.054m、宽度为0.01m 的线圈(图 1)。计算模型如图 2 所示,内径为 0.05m,外径为 1.0m,宽度为 0.01m,前表面与后表面采用周期边界条件。网格总量在 300000、360000 和 420000 时,流速相差在 0.15%以内,因此为了同时兼顾计算的准确性和时效性,选用网格数 360000 进行计算,其中 z=0m 平面的网格如图 3 所示,图中 x 和 y 表示坐标。0.01 m0.05 m0.054 m图1实验用线圈示意图9Fig.1Schematicd
19、iagramoftestcoil9Computational zonexzy1.0 m0.05 m0.01 mCoilVswBack boundaryFront boundary图2计算模型示意图Fig.2Schematicdiagramofcalculationmodel当计算区域半径远大于线圈半径时,环形导电线圈产生的磁场可近似为偶极子磁场19。图 4给出了 z=0m 平面磁场模型。Vswni采用文献 9 中的边界条件对文中方法进行验证。太阳风沿 y 轴负方向流入,线圈攻角为 0(即线圈轴线与来流方向垂直),如图 4 所示。太阳风来流速度=45km/s、离子数密度=1.8852航空动力学报
20、第38卷miBc1019m3、质子质量=1.6721027kg、线圈中心磁场强度为 1.5T。V图 5 给出了太阳风等离子体与线圈磁场作用后的速度()和磁场流线分布与文献实验结果对比,从图中可观察到,由于受到洛伦兹力阻碍作用而形成的激波层与实验中的激波层位置一致。L数值计算结果显示磁腔半径 约为 0.14m,由式(11)式(13)计算得到的结果与文献 9 中实验测得的结果见表 1。对比发现,理论计算得到的推力值大于实验测得的推力值,这主要是因为理论计算的推力系数是在无碰撞条件下得到的,而实验中所采用的等离子体是部分电离,且实验设备尺寸相对较小导致无法忽略带电粒子之间和带电粒子与中性粒子之间的碰
21、撞,因而造成了数值模拟结果和实验结果之间在推力和推力系数上的偏差20。但数值模拟得到的磁腔半径、理论计算得到的拉莫尔半径与磁腔半径之比与实验结果具有较好的一致性,验证了计算方法的可靠性。3计算结果及分析本文采用理想氢气等离子体作为工质。由于采用 MHD 模型计算,因此所讨论的流场主要表现在宏观特征,即流体受到洛伦兹力的作用。由动量方程(式(2)可知磁帆所受洛伦兹力在 y 方向上的分量为Fy=Jz(B0 x+bx)Jx(B0z+bz)(14)JxJzB0 xB0zbxbz(B0 x+bx)其中和分别为 x、z 方向上的电流密度分量,和分别为外加磁场在 x、z 方向上的分量,和分别为诱导磁场在 x
22、、z 方向上的分量。图 6(a)给出了 x=0m 中线上各方向磁场分量值对比,相较于 x 方向总的磁场分量,z 方向上的磁场分量几乎为零,因此洛伦兹力 y 方向分0.500.500.0300.250.250.025000.0200.250.250.0150.500.0300.0200.010 x/mx/my/my/m图3z=0m 平面区域计算网格示意图Fig.3Schematicdiagramofcalculationgridforz=0mplane0.500.500.2500.2500.250.250.50 x/my/m图4 z=0m 平面外加磁场流线分布Fig.4Streamlinedis
23、tributionofexternalmagneticfieldonz=0mplane09000180002700036000450000.30.40.20.100.10.2y/m4050022500090004500Bow shock180003150027000450001350036000V/(m/s)图5数值模拟结果与文献 9 中实验结果的对比图Fig.5ComparisonofnumericalresultswithexperimentalresultsinRef.9表1计算结果与文献 9 中实验结果Table1Calculationresultsandexperimentalres
24、ultsinRef.9参数计算值实验值L/m0.140.15riL/m0.038riL/L0.270.24Cd3.527F/N6.621.64第4期蔡静媛等:太阳风磁帆推进性能数值研究853Jz(B0 x+bx)量的大小主要由式(14)中的第一项即决定,图 6(b)给出了 x 方向上磁场分布。当线圈磁场不变时,影响洛伦兹力 y 方向分量的主要因素是 z 方向上的电流大小。0.060.500.500.050.250.250.04000.030.250.250.020.500.500.2500.250.500.010B/Ty/my/mx/mB0 x+bxB0y+byB0z+bz(b)x方向磁场分量
- 配套讲稿:
如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。