Loading [MathJax]/jax/output/SVG/jax.js

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

高功率高储能脉冲电源中重频电感的设计与分析

王燕 张钦 林福昌 李化

杨宸, 范征锋, 刘杰, 等. 通过理论定标律研究激光间接驱动的参数空间[J]. 强激光与粒子束, 2019, 31: 062001. doi: 10.11884/HPLPB201931.180339
引用本文: 王燕, 张钦, 林福昌, 等. 高功率高储能脉冲电源中重频电感的设计与分析[J]. 强激光与粒子束, 2022, 34: 055002. doi: 10.11884/HPLPB202234.210417
Yang Chen, Fan Zhengfeng, Liu Jie, et al. Parameter spaces of laser indirect-driven target by theoretic scaling laws[J]. High Power Laser and Particle Beams, 2019, 31: 062001. doi: 10.11884/HPLPB201931.180339
Citation: Wang Yan, Zhang Qin, Lin Fuchang, et al. Design and analysis of repetitive frequency inductor in high power and high energy storage pulse power supply[J]. High Power Laser and Particle Beams, 2022, 34: 055002. doi: 10.11884/HPLPB202234.210417

高功率高储能脉冲电源中重频电感的设计与分析

doi: 10.11884/HPLPB202234.210417
详细信息
    作者简介:

    王 燕,yan7037@mail.hust.edu.cn

    通讯作者:

    张 钦,zhangqin@mail.hust.edu.cn

  • 中图分类号: TM83

Design and analysis of repetitive frequency inductor in high power and high energy storage pulse power supply

  • 摘要: 针对高功率脉冲电源集成系统连续放电的需求,研制了一种用于该系统的重频水冷电感。水冷电感在系统中既调节电源的电流波形,又能起到隔离作用。该重频水冷电感耐压高、通流大、充电间隔周期6 s,可连续工作10 次放电。针对连续放电的需求,通过去离子循环水对重频电感进行降温。现以单模块能量334 kJ、电感量30 µH、通流 100 kA为例进行设计分析,建立重频水冷电感温度场仿真模型,通过ANSYS仿真软件对该电感内部温度瞬态特性进行分析。结果表明:水冷电感通过去离子水冷却效果好,电感每次工作后温度最高达到47 ℃,在下一个工作点到来之前电感温度又恢复到41 ℃;同时该电感在没有加循环水的情况下通过了112 kA的电动力考核。试验结果与理论分析吻合较好,重频水冷电感运行稳定,从而验证了理论分析和设计的正确性。
  • 惯性约束聚变(ICF)[1-4]的目标是约束氘氚聚变燃料并点火,实现可控核聚变反应并产生可观的能量。激光间接驱动[5-6]是目前最流行的方式之一。将激光注入装有氘氚燃料靶丸的金或铀黑腔中转换为X射线辐射,辐射均匀地辐照到靶丸上,烧蚀驱动靶丸内爆实现氘氚燃料的高密度压缩,并在靶丸芯部产生高温热斑,达到热核点火条件,并且在燃料飞散前充分燃烧、产生足够的聚变能量。2018年美国劳伦斯利弗莫尔国家实验室(LLNL)在国家点火装置(NIF)[7]上实现了聚变能输出大于内爆动能的实验[8]。我国2016年实现了间接驱动内爆集成实验[9],2018年在黑腔中进行了内爆实验[10]。激光聚变间接驱动目前的研究难点是控制黑腔内的激光等离子体不稳定性(LPI)、控制内爆中的流体力学不稳定性以及解决内爆混合[11-14]的问题。黑腔中靶丸加持膜[15]和充气管引起的不稳定性[16]也限制了内爆压缩达到高增益点火条件。

    直接数值模拟是点火靶物理研究与设计的重要手段,为此美国发展了LASNEX[5], HYDRA[17]等辐射流体力学模拟程序,我国也发展了系列辐射流体力学程序[18]。这些程序物理建模复杂、代码量庞大、网络分辨率高,具有较高的模拟置信度。然而,对数值模拟的过多依赖导致ICF内爆过程中的基本物理机制无法简明扼要地获得,不利于系统地分析内爆过程中影响靶丸能量增益的基本物理机制和规律。

    本文的工作是用定标方法建立简单、易懂的内爆动力学模型,分析影响内爆和能量增益的基本物理机制,在物理规律上研究设计参数对靶丸性能的影响。内爆加速阶段使用火箭模型[1, 2, 19]描述,减速阶段依靠劳森判据和定标律[20]描述,热斑形成与点火采用一维流体力学开源程序MULTI-IFE[21]辅助定标。由此获得靶丸设计所需的参数空间,研究参数空间在不同初始条件下的变化规律。在一定的靶丸性能要求下,给出了靶丸基本的设计参数范围。当改变靶丸的驱动辐射温度、飞行熵增因子、整形速度和烧蚀体材料的条件下,靶丸的装量、半径以及壳层厚度应当作相应地调整。

    图 1(a)所示,实际的间接驱动点火靶丸由同心的球型壳层组成。壳层分为外层烧蚀体和内层主燃料两部分。外层烧蚀体由聚苯乙烯塑料(CH)或金属铍(Be)等低原子序数物质制成,为了防止预热、控制流体力学不稳定性等,烧蚀体一般由不同的掺杂层构成,掺杂的元素有硅、锗、钨等。内层主燃料为冷冻氘氚(DT-ice)固体。壳层内部的空腔中充满稀薄的氘氚蒸汽(DT-gas)。壳层厚度比靶丸半径小很多,经过整形驱动压缩后,壳层变得更薄,因此靶丸壳层可以假设成一个具有初始装量M0且位置在R0处的薄壳,如图 1(b)所示。在薄壳假设的基础上,1.1至1.3部分分阶段对靶丸内爆过程进行物理建模。

    图  1  实际靶丸构型与简化模型
    Figure  1.  The real target configuration and the simplified model

    当靶丸受到辐射烧蚀时,壳层的一部分烧蚀体向外喷射,反作用力推动壳层聚心内爆。辐射对靶丸壳层产生的质量烧蚀速率和烧蚀压为$ \dot m$=aTr3pa=bTr3.5,其中$ \dot m$的单位是g·cm-2·μs,pa的单位是Mbar,辐射温度Tr的单位是heV。烧蚀定标律系数a, b由烧蚀体的材料决定:(1)文献[5]式(6-2)(6-3)给出了CH的烧蚀定标律系数;(2)文献[22]给出了Ge-doped CH、HDC和Cu-doped Be的质量烧蚀速率与辐射温度之间的实验数据图并拟合其数据给出这三种材料的系数a;文献[22]也给出了HDC和Cu-doped Be的烧蚀压与辐射温度之间的数据图,拟合可得bHDC=5;bBe=6.31,但是未给出Ge-doped CH烧蚀压数据。为了方便计算和讨论,我们假设有b/a=10。由此获得表 1中的数据。

    表  1  不同材料的烧蚀定标律系数
    Table  1.  The parameters for different ablator materials
    materials a/(g·cm-2·μs·heV-3) b/(Mbar·heV-3.5)
    CH 0.30 3.0
    Ge-doped CH 0.35 3.5
    HDC 0.50 5.0
    Cu-doped Be 0.75 7.5
    下载: 导出CSV 
    | 显示表格

    驱动壳层内爆的辐射由整形脉冲和主脉冲两部分构成。本文从主脉冲驱动的加速阶段开始研究,假设整形脉冲结束后,壳层获得初始内爆速度为v0。壳层的初始半径位置设为R0,壳层的初始装量(包括烧蚀层和DT主燃料)设为M0。随后,辐射主脉冲烧蚀靶丸壳层并驱动其加速内爆,根据质量守恒、牛顿第二定律可得[19]

    {dM(t)dt=4πR2(t)˙mM(t)dv(t)dt=4πR2(t)padR(t)dt=v(t) (1)

    由式(1)有dv(t)=-(pa/$ \dot m$)(dM(t)/M(t)),对其积分可得加速到t时刻的内爆速度v=v0+(pa/$ \dot m$)ln(1/δt)。物理上,随着壳层向内聚心,靶丸体积逐渐变小,内部气体压力逐步上升,当内部气体压力等于靶丸驱动压力时,靶丸壳层获得最大内爆速度,之后壳层开始减速。最大内爆速度时刻,记壳层剩余质量份额为δ=ms/M0,其中ms为剩余壳层质量,则最大的内爆速度为vimp=v0+(pa/$ \dot m$)ln(1/δ)。

    仍由式(1)得

    C0(1δ)+δlnδ=˙m2paM04π3R30=D0;C0=v0˙mpa+1 (2)

    用初始参量表达最大内爆速度有vimp=v0-10lnδTr0.5=v0δ-1+[10Tr0.5(δ-1-1)-($ \dot m$/M0)(4π/3)R03δ-1]。进一步简化计算,将此加速过程近似为均加速的,使用v≈(v0+vimp)/2可得壳层的加速时间τ=R0/v=R0/(v0-5lnδTr0.5)和平均加速度g=2(R0-v0τ)/τ2=10(5lnδTr0.5-v0)lnδTr0.5/R0

    假设加速过程中密度为常数,壳层压力等于烧蚀压,加速过程飞行熵增因子为αif,可得壳层密度为ρ=[pa/(αifAdeg)]3/5=Adeg-0.6b0.6αif-0.6Tr2.1,其中Adeg=2.17×1012 (erg/cm3)/(g/cm3)5/3是氘氚固体的特征常数[2]。壳层的初始厚度为ΔR0=M0/(ρR02)=(4π)-1Adeg0.6b-0.6αif0.6Tr-2.1M0R0-2,靶丸的初始形状因子为Ar0=R0R0=(4π)1Adeg-0.6b0.6αif-0.6Tr2.1M0-1R03

    在靶丸加速过程中,烧蚀面上压力梯度和密度梯度反向,引发烧蚀瑞利—泰勒不稳定性[23](A-RTI)增长,此时界面附近的物理量可表示为$ f(x, y ; t)=f_0(x)+\sum\limits_k f_k(x) \mathrm{e}^{\mathrm{i} k y} \mathrm{e}^{\gamma_t t}$,其中

    γ1=0.9k1g1+k1L1.4k1va (3)

    是A-RTI线性增长率[24-25]。不稳定性增长率由RT增长项$ 0.9 \sqrt{k_1 g /\left(1+k_1 L\right)}$和烧蚀致稳项-1.4klva组成,其中波数为kl=2π/λ=l/R=2lR0-1; R=R0/2。式(3)中L表示密度梯度标长[24],本文用辐射自由程[2]来近似,即L=ρ/(dρ/dr)=1/(κpρ),其中κp=0.43(Z3/A2)ρTr-7/2=0.43Z3A-2Adeg-0.6b0.6αif-0.6Tr-1.4表示辐射的平均光厚[2]。烧蚀致稳项中的烧蚀速度可定义为

    va=M0(1δ)4πρ(34R0)2τ=49πA0.6deg(v05lnδT0.5r)(1δ)b0.6α0.6ifT2.1rM0R40 (4)

    这里近似加速阶段壳层的表面积从4πR02收缩到0,故va中有R=3R0/4。A-RTI发展的时间就是壳层加速内爆的时间τ。则A-RTI线性增长的指数因子为

    γlτ=0.94l1+4.65lZ3A2A1.2degb1.2α1.2ifT0.7rR105lnδT0.5r(5lnδT0.5rv0)0.396lA0.6deg(1δ)b0.6α0.6ifT2.1rM0R40 (5)

    l模扰动线性增长的放大倍数为Gfl=eγlτ,记Gfl=1~100阶模放大倍数的最大值。设计上要求A-RTI的幅值不超过壳层厚度,则所容忍的最大初始扰动幅值可以近似为ηR0/Gf

    减速过程[19, 26]中,壳层对气体做功,内部气体的压力和能量升高,壳层的动能减小。当壳层动能消耗殆尽时,达到阻滞时刻,之后靶丸开始向外飞散。阻滞时刻,壳层的熵增因子记为αs(称为阻滞熵增因子),其与飞行熵增因子的关系为$ \alpha_{\mathrm{s}}=\alpha_{\mathrm{if}} \sqrt{M}=\alpha_{\mathrm{if}}^{0.85} v_{\mathrm{imp}}^{0.5} p_{\mathrm{a}}^{-0.1} \gamma^{-0.25} A_{\mathrm{deg}}^{-0.15}$,其中$ M=v_{\mathrm{imp}} / c_{\mathrm{if}}=v_{\mathrm{imp}} / \sqrt{\gamma p_{\mathrm{a}} / \rho}$是最大内爆速度时刻的马赫数[26, 27]cif是最大内爆速度时刻壳层的声速。又记阻滞时刻靶丸的半径为rs、壳层的密度和压力分别为ρsp。并且我们假设阻滞时刻:(1)壳层和热斑等压,热斑的质量可忽略;(2)壳层压力可用阻滞熵增因子与壳层密度表达;(3)用理想气体的内能描述壳层的内能。我们还假设最大内爆速度时刻壳层的动能可以全部转化为阻滞时刻壳层的内能。由这些假设可得

    {ms=4π3r3sρsp=αsAdegρ5/3s4π3r3spγ1=12msv2imp (6)

    化简后有

    {p=(γ12)5/2v5impα3/2sA3/2degρs=(γ12)3/2v3impα3/2sA3/2degrs=(4π3)1/3(γ12)1/2m1/3sv1impα1/2sA1/2deg (7)

    再将阻滞熵增αs换为飞行熵增αif,可得

    {p=(γ12)2.5γ0.375A1.275degv4.25impα1.275ifp0.15aρs=(γ12)1.5γ0.375A1.275degv2.25impα1.275ifp0.15ars=(4π3)1/3(γ12)0.5γ0.125A0.425degm1/3sv0.75impα0.425ifp0.05a (8)

    所以,我们获得了壳层内爆到阻滞时刻的物理量p, ρs, rs, αs对于最大内爆速度时刻的物理量ms, vimp, αif, pa和常数Adeg, γ的定标律。

    在理想情况下,阻滞或阻滞时刻之前,壳层内部(即靶丸中心)已形成点火热斑,引发聚变点火与燃烧。这样的燃烧从中心热斑开始,逐渐向外传播并引燃DT燃料壳层,进而引发所有的燃料燃烧。对于燃烧中的DT燃料等离子体,热斑和壳层等压且有均匀的压强p,此时约束热斑和壳层等离子体的是物质惯性。随后,这团等离子体将以稀疏波的形式向外膨胀,膨胀速度为阻滞时刻的声速cs,则惯性约束时间[4]可表示为

    τc=rscs=rsγp/ρs (9)

    用最大内爆速度时刻的物理量表示它,得到

    τc=(4π3)1/3(γ12)1γ0.625A0.425degm1/3sv1.75impα0.425ifp0.05a (10)

    我们需要进一步描述靶丸燃烧后的聚变反应产能和能量增益,这需要知道中心热斑和燃料壳层的总面密度。对于一个聚变点火并燃烧的靶丸,要求中心点火热斑的面密度(ρr)h>0.3 g/cm2、离子温度Th>4 keV且燃料壳层的面密度(ρr)s>1 g/cm2[28]。为了得到总面密度(包括热斑面密度和壳层面密度)的定标律[3],我们用MULTI-IFE程序模拟了从最大内爆速度时刻开始到靶丸最后飞散的过程中,热斑和壳层各个物理量的演化。我们在如下范围内,随机生成500组初始参数:壳层的初始质量ms=200~400 μg、密度ρs=10~30 g/cm3、温度Ts=1~31 eV、速度vimp=200~400 km/s,热斑的半径rh=50~300 μm、密度ρh=0.01~0.26 g/cm3、温度Th=0.1~1.6 keV。统计MULTI-IFE的输出结果后,得到了中间变量和总面密度ρr之间的数据,中间变量有壳层的质量ms,动能ks=msvimp2/2,烧蚀压pa,飞行熵增因子αif,热斑的体积Vh=4πrh3/3,质量mh=ρhVh,压力ph=ΓBρhTh,其中ΓB=7.66×1014 erg/(g·keV)是氘氚气体常数。拟合数据得到

    ρr=(4.5927×106)m0.476sk0.619sp0.0325aα0.369ifV0.603hm0.0975hp0.274h (11)

    图 2所示是拟合结果的准确性效果,横坐标是MULTI-IFE计算得到的总面密度,纵坐标是用拟合后的定标律给出的数值,可见点火面密度的定标律式(11)拟合的很好。再忽略式(11)中次方小的项,可以得到ρr=(8.2968×10-5)ms0.476ks0.624αif-0.342Vh-0.694pa-0.327,由于加速阶段的火箭模型无法确定最大内爆速度时刻热斑的参数,所以取Vh为500组随机Vh的平均值,最后得总面密度的定标律为

    ρr=(0.087325)m0.476sk0.624sα0.342ifp0.327a=(0.056663)(M0δ)1.100v1.248impα0.342ifp0.327a (12)
    图  2  500组随机初始参数拟合的定标面密度效果
    Figure  2.  The scaling MULTI-IFE data vs of areal density

    文献[3]从广义劳森判据推得了总面密度定标律,文献[25]从流体动力学相似律等价理论出发推导得到了壳层面密度和热斑面密度的定标律。

    由总面密度,我们可知氘氚燃料被点燃的比例即燃烧损耗[2, 25]

    Φ=ρrρr+HB (13)

    其中HB=6 g/cm2是燃烧参数,氘氚燃料反应后的能量产额

    Q=QDTms5mpΦ=(3.37184×1018erg/g)msΦ (14)

    其中QDT=17.6 MeV是一对氘氚原子发生核反应放出的能量,mp是单位原子质量。最后,靶丸内爆的能量增益为

    G=Q2.746F (14)

    其中,已将通过靶丸初始表面的辐射能F=4σTr4πR02τ=4πσTr4R03/(v0-5lnδTr0.5)耦合为间接驱动给与的激光能量Elaser=2.746F,即激光转换效率是36.42%,这个效率参考了文献[5, 30]的理论与数据。

    我们还需要判断靶丸内爆点火的裕量,即是否能顺利点火。由劳森判据[2-4],阻滞时刻壳层的压强p和约束时间τc乘积必须满足

    pτc=(4π3)1/3(γ12)1.5γ0.25A0.85degm1/3sv2.5impα0.85ifp0.1a>C (16)

    所以对于能顺利点火的靶丸,其壳层加速完成时的质量

    ms>mmin=C3(4π3)[(γ12)4.5γ0.75A2.55deg]v7.5impα2.55ifp0.3a (17)

    故最小点火能量

    Eigt=C32(4π3)[(γ12)4.5γ0.75A2.55deg]v5.5impα2.55ifp0.3a (18)

    定义出一维点火阈值因子

    ITF1D=12msv2impEigt=[C3(4π3)1(γ12)4.5γ0.75A2.55deg]msα2.55ifv7.5impp0.3a (19)

    ITF1D>1的靶丸是可以点火的靶丸。再使用ms=M0δ,并且在一个典型的点火点(ms=0.3 mg; vimp=360 km/s; αif=1.4; pa=131.6 Mbar; ITF1D=3.0)上确定常系数后有

    ITF1D=(2.9145×1057)M0δv7.5impα2.55ifp0.3a (20)

    ITF1D描述了靶丸能否顺利点火。文献[3]从广义劳森判据推得了点火阈值因子,文献[29]讨论了考虑热斑扰动和热斑纯度后的点火阈值因子。

    通过上述内爆动力学模型,我们获得了靶丸设计参量与性能参量之间的定标关系。关键的靶丸设计参量有:靶丸烧蚀体的材料系数a, b,靶丸的装量M0与初始半径R0及其关联的剩余质量份额δ与形状因子Ar0,靶丸所处环境参量即黑腔的辐射温度Tr,靶丸的内爆过程参量飞行熵增因子αif与初始内爆速度v0。关键的靶丸性能参量有:点火阈值因子ITF1D,核反应能量增益G,通过靶丸初始表面的辐射能F,靶丸的不稳定性增长因子Gf和允许的最大初始扰动η。通过分析这些设计参量与性能参量之间的变化关系,我们研究了靶丸设计参数空间的变化规律。

    靶丸的性能参量随靶丸装量——半径的变化规律是基本的关注点。在确定的靶丸烧蚀材料CH,飞行熵增因子αif=1.4,辐射温度Tr=3 heV和初始内爆速度v0=0 km/s的条件下,我们获得图 3。由图 3(a),对于确定的靶丸装量,靶丸的初始半径越大,被烧蚀的外表面越大,壳层加速完成时剩余质量份额越小。为了保留足够的燃料质量,需要δ>0.1。由图 3(b),点火阈值因子ITF1D>1才能顺利点火,考虑一定裕量,取ITF1D>2。δ>0.1且ITF1D>2的区域即为靶丸点火所考虑的范围,如图 3(c)以黄色所示,即点火岛区域。

    图  3  剩余质量比δ和点火阈值因子ITF1D在靶丸装量——半径参数空间中的变化及点火岛区域
    Figure  3.  The contour of residual mass fraction δ and ignition threshold factor ITF1D in the capsule mass-radius parameter space and the volcanic islands

    考虑点火岛中的性能参量。如图 4所示,初始装量为3 mg,初始半径为1 mm是主要参考点。该点靶丸需要消耗1.49 MJ激光能量转换为0.542 MJ通过靶丸初始表面的X射线能量,预计会产生8.46 MJ的聚变能得到5.69倍的能量增益。并且该靶丸在获得386.5 km/s的最大内爆速度时,壳层剩余质量份额为10.7%。靶丸的点火阈值因子ITF1D=5.59,形状因子是41.76。加速阶段产生的烧蚀瑞利—泰勒不稳定性增长因子为346.3,故最大初始扰动不能超过0.07 μm。

    图  4  点火岛中的参量
    Figure  4.  The parameters in the volcanic islands

    辐射温度Tr由3 heV增大到3.3 heV,质量烧蚀速率和烧蚀压增大,对装量相同的靶丸,为了保证烧蚀效果不变,必须减小靶丸受到烧蚀的面积,使得壳层所有的质量不被辐射烧蚀掉。如图 5所示,点火岛位置将向靶丸初始半径减小的方向移动,点火岛区域大小变化不显著。

    图  5  辐射温度对点火岛的影响
    Figure  5.  The effect of radiation temperature on the volcanic islands

    图 6所示,增大辐射温度后,虽然质量烧蚀速率和烧蚀压增大,但是在相同剩余质量份额的条件下,加速时间减小,使得通过靶丸初始表面的辐射能基本不变,需要的激光能量不变。靶丸产生的聚变能量与辐射温度无关,致使提高辐射温度后能量增益不变。增大辐射温度后,不稳定性的烧蚀致稳作用增强,不稳定性增长率和发展时间都变小,靶丸的稳定性变好,可以容忍更大的初始扰动。设计上,在靶丸装量不变的条件下,需要减小靶丸的半径,同时保证形状因子不变,所以靶丸的厚度也相应需要减小。

    图  6  不同辐射温度下性能参量的比较
    Figure  6.  The comparisons of parameters in different radiation temperatures

    飞行熵增因子αif由1.4增大到2.0,虽然不影响壳层的最大内爆动能,但是阻滞时刻可以获得的压力变小密度变小,最小点火能量的增大使得点火阈值因子变小,最终导致点火岛变窄。如图 7所示,点火岛的位置不变,但是范围变小。

    图  7  飞行熵增对点火区域的影响
    Figure  7.  The effect of in-flight adiabat on the volcanic islands

    图 8所示,增大飞行熵增因子后,虽然通过靶丸初始表面的辐射能不变,需要的激光能量不变,但是点火面密度略微减小,产能因此减小,所以能量增益略微减小。增大飞行熵增因子后,加速阶段烧蚀体被压缩的密度减小,其平均光厚减小,密度梯度标长增大,烧蚀速度增加,故此烧蚀致稳作用增强。不稳定性增长因子变小,允许的最大初始扰动变大,内爆稳定性变好。另外,由于减小了壳层的压缩密度,设计上需要增加靶丸的初始厚度,形状因子因此变小。

    图  8  不同飞行熵增下性能参量的比较
    Figure  8.  The comparisons of parameters in different in-flight adiabat

    整形速度v0由0 km/s增大到40 km/s,根据式(2)在靶丸装量不变的条件下,相同剩余质量份额对应的靶丸初始半径增大,又由于整形速度对最大内爆速度的贡献,点火阈值因子提高。如图 9所示,点火岛的位置向靶丸初始半径增大的方向移动,点火岛的区域范围变大。

    图  9  整形速度对点火区域的影响
    Figure  9.  The effect of shaping velocity on the volcanic islands

    图 10所示,增大初始整形速度后,虽然最大内爆速度增大,导致最大内爆动能增大。但是装量相同的靶丸初始半径增大,导致通过靶丸初始表面的辐射能增大,需要的激光能量增大。虽然内爆速度的增大对点火面密度有贡献,致使产能提高,但是权衡也增大的激光能量后,能量增益基本不变。增大初始整形速度后,内爆不稳定性增长因子的增加不超过$ \sqrt{v_{{\rm{i m p}}} /\left(v_{{\rm{i m p}}}-v_0\right)}$倍,对vimp=300~400 km/s, v0=40 km/s,不稳定性增长因子增加不会超过5.4%~7.4%,即内爆稳定性变化不显著。另外,由于相同装量下靶丸初始半径增大,设计上壳层初始厚度变小、形状因子变大。壳层厚度的变小使得允许的最大初始扰动变小。

    图  10  不同整形速度下性能参量的比较
    Figure  10.  The comparisons of parameters in different shaping velocity

    当改变靶丸的烧蚀材料,提高质量烧蚀速率与烧蚀压时,根据式(2)对于相同的δ参数D0是不变的,质量烧蚀速率$ \dot m$的提高使得R0减小,有R0~a-1/3b~R0-3。点火阈值因子不受影响,只是跟着δ移动位置。如图 11所示,点火岛向靶丸初始半径减小的方向移动,点火岛的区域范围大小不变。

    图  11  靶丸烧蚀材料对点火区域的影响
    Figure  11.  The effect of ablator material on the volcanic islands

    图 12所示,当改变靶丸的烧蚀材料,提高质量烧蚀速率与烧蚀压时,由于初始半径R0明显减小, 通过靶丸初始表面的辐射能明显减小,需要的激光能量明显减小。点火面密度不受影响产能不变,由于节约了驱动能量,所以能量增益变大。另外,虽然烧蚀材料烧蚀性能的改变,增大了烧蚀体的密度ρ~b0.6,使得烧蚀体的平均光厚增大κp~b0.6,密度梯度标长减小L~b-1.2,但是初始半径减小的效应更强。最终效果是增大了烧蚀速度va~b-0.6R0-4~b0.73,增加了烧蚀致稳~b1.06的效果,不稳定性增长因子明显减小,可以容忍的初始扰动增大,稳定性明显变好。设计上,对于装量相同的靶丸由于初始半径的减小,初始厚度需要增加,靶丸形状因子减小。

    图  12  不同靶丸烧蚀材料下性能参量的比较
    Figure  12.  The comparisons of parameters in different ablator material

    本文通过理论方法和一维流体力学模拟,给出了激光间接驱动点火靶丸内爆过程的定标律公式。通过这些定标律公式获得了在给定黑腔辐射温度、飞行熵增因子、整形速度和烧蚀材料的条件下,靶丸装量——半径参数空间的点火岛区域。研究了靶丸性能参数随辐射温度、飞行熵增因子等的变化规律:(1)当靶丸所处黑腔辐射温度升高时,内爆的稳定性将变好;设计上在靶丸装量不变的条件下,靶丸半径需要减小。(2)当靶丸的飞行熵增因子增大时,内爆增益略微减小,内爆稳定性变好;但是点火阈值因子减小导致点火岛的区域变窄。(3)当靶丸的整形速度增大时,点火岛的区域略微变大,内爆稳定性变化不显著;设计上在靶丸装量不变的条件下,需要增大靶丸半径,这会导致靶丸壳层形状因子变大。(4)当改变靶丸烧蚀材料,提高质量烧蚀速率与烧蚀压时,能量增益变大且稳定性增强;设计上在靶丸装量不变的条件下,需要减小靶丸半径。

    通过这些分析,我们获得了改进靶丸性能的几种方案:增大辐射温度、增大飞行熵增因子可增强A-RT不稳定性的烧蚀致稳作用,减弱内爆不稳定性;增大整形速度,可提高产能,但是并不能有效提高增益;改变靶丸的烧蚀材料提高质量烧蚀速率与烧蚀压,可以节约激光能量提高能量增益并控制不稳定性。在后续工作中,我们将进一步精细化分析影响靶丸增益和内爆不稳定性的物理机制,寻找更佳的优化设计方案。

  • 图  1  单电源模块结构图

    Figure  1.  Single module structure diagram

    图  2  重频水冷电感外形图

    Figure  2.  Contour diagram of the repeatitive frequency water-cooled inductor

    图  3  重频水冷电感三维图

    Figure  3.  3D diagram of the repeatitive frequency water- cooled inductor

    图  4  重频水冷电感线圈

    Figure  4.  Diagram of repeatitive frequency water-cooled inductor coil

    图  5  1/4 匝导线模型

    Figure  5.  1/4 turn traverse model

    图  6  重频水冷电感温度分布云图

    Figure  6.  Cloud map of temperature distribution of the repeatitive frequency water-cooled inductor

    图  7  电感线圈内部导体上的应力

    Figure  7.  Stress on conductor inside inductor coil

    图  8  电感线圈导体表面电磁密度

    Figure  8.  Surface electromagnetic density of inductor coil conductor

    图  9  重频水冷电感通流试验原理图

    Figure  9.  Schematic diagram of repeatitive frequency water-cooled inductor test circuit

    图  10  重频水冷电感通流试验接线图

    Figure  10.  Wiring diagram of repeatitive frequency water cooled inductor flow test

    图  11  重频水冷电感通过了112 kA电流峰值的波形图

    Figure  11.  Waveform of repeatitive frequency water-cooled inductor passing a peak current of 112 kA

    表  1  导线线径与温升的关系

    Table  1.   Relationship between wire diameter and temperature rise

    numberФo/mmФi/mmenergy/kJtemperature/℃
    114626752.96
    下载: 导出CSV

    表  2  趋肤深度与铜管水温、频率的关系

    Table  2.   Relationship between skin depth and temperature and frequency

    numberfrequency/Hztemperature/℃temperature of coefficient/Kskin depth/mm
    120047/40
    47/40
    1.1154/1.0853
    1.1154/1.0853
    3.707
    2166.74.117
    下载: 导出CSV

    表  3  电感流量特性实验数据表

    Table  3.   Experimental data of inductance flow characteristics

    flow rate/(L/min) pressure loss/MPa
    5 0.257
    5.5 0.309
    6 0.364
    6.5 0.413
    7 0.512
    下载: 导出CSV
  • [1] Akiyama H, Sakugawa T, Namihira T, et al. Industrial applications of pulsed power technology[J]. IEEE Transactions on Dielectrics and Electrical Insulation, 2007, 14(5): 1051-1064. doi: 10.1109/TDEI.2007.4339465
    [2] 戴宇峰, 鲁军勇, 张晓, 等. 脉冲功率电源连续发射水冷模拟负载[J]. 国防科技大学学报, 2016, 38(6):6-11. (Dai Yufeng, Lu Junyong, Zhang Xiao, et al. Water-cooling simulated resistance for continuously launching pulsed power supply[J]. Journal of National University of Defense Technology, 2016, 38(6): 6-11 doi: 10.11887/j.cn.201606002
    [3] Hundertmark S, Liebfried O. Options for an electric launcher system[J]. IEEE Transactions on Plasma Science, 2019, 47(10): 4433-4438. doi: 10.1109/TPS.2018.2888899
    [4] Li Zhenxiao, Zhang Yazhou, Wu Jinguo, et al. Design of a 30MJ capacitor-based pulsed power supply for EML[C]//Proceedings of the 21st International Conference on Pulsed Power. Brighton, UK: IEEE, 2017: 1-5.
    [5] Lehmann P, Peter H, Wey J. First experimental results with the ISL 10 MJ DES railgun PEGASUS[J]. IEEE Transactions on Magnetics, 2001, 37(1): 435-439. doi: 10.1109/20.911871
    [6] Hundertmark S, Schneider M, Simicic D, et al. Experiments to increase the used energy with the PEGASUS railgun[J]. IEEE Transactions on Plasma Science, 2014, 42(10): 3180-3185. doi: 10.1109/TPS.2014.2314131
    [7] 张龙霞, 李碧清, 霍敏. 国外电磁炮发展概述[J]. 飞航导弹, 2011, 40(10):23-27. (Zhang Longxia, Li Biqing, Huo Min. Overview of foreign electromagnetic gun development[J]. Maneuverable Missile, 2011, 40(10): 23-27
    [8] 范晶, 宋朝文. 舰载电磁轨道炮用高功率脉冲电源研究进展[J]. 电气技术, 2010(s1):70-72. (Fan Jing, Song Chaowen. Research progress of the high-powered pulsed power supply of naval electromagnetic rail gun[J]. Electrical Engineering, 2010(s1): 70-72
    [9] 何巧, 马游春, 马子光, 等. 高功率脉冲电源电压测试系统设计[J]. 仪表技术与传感器, 2020(2):46-49,113. (He Qiao, Ma Youchun, Ma Ziguang, et al. Design of high power pulse power supply voltage test system[J]. Instrument Technique and Sensor, 2020(2): 46-49,113 doi: 10.3969/j.issn.1002-1841.2020.02.011
    [10] Conway J T. Inductance calculations for circular coils of rectangular cross section and parallel axes using Bessel and Struve functions[J]. IEEE Transactions on Magnetics, 2010, 46(1): 75-81. doi: 10.1109/TMAG.2009.2026574
    [11] Conway J T. Inductance calculations for noncoaxial coils using Bessel functions[J]. IEEE Transactions on Magnetics, 2007, 43(3): 1023-1034. doi: 10.1109/TMAG.2006.888565
    [12] Conway J T. Analytical solutions for the self- and mutual inductances of concentric coplanar disk coils[J]. IEEE Transactions on Magnetics, 2013, 49(3): 1135-1142. doi: 10.1109/TMAG.2012.2229287
    [13] Luo Yao, Chen Baichao. Improvement of self-inductance calculations for circular coils of rectangular cross section[J]. IEEE Transactions on Magnetics, 2013, 49(3): 1249-1255. doi: 10.1109/TMAG.2012.2228499
    [14] Babic S, Salon S, Akyel C. The mutual inductance of two thin coaxial disk coils in air[J]. IEEE Transactions on Magnetics, 2004, 40(2): 822-825. doi: 10.1109/TMAG.2004.824810
    [15] 卡兰塔罗夫, 采依特林. 电感计算手册[M]. 陈汤铭, 刘保安, 罗应立, 等译. 北京: 机械工业出版社, 1992

    КалантаровП. Л. , ЦейтлинЛ. А. Расчет индуктивностей справочаи книга[M]. Chen Tangming, Liu Baoan, Luo Yingli, et al, trans. Beijing: China Machine Press, 1992
    [16] 罗垚, 陈柏超, 袁佳歆, 等. 倾斜轴空心矩形截面圆柱线圈互感计算[J]. 电工技术学报, 2012, 27(5):132-136. (Luo Yao, Chen Baichao, Yuan Jiaxin, et al. Mutual inductance calculations of inclined axial air-core circular coils with rectangular cross-sections[J]. Transactions of China Electrotechnical Society, 2012, 27(5): 132-136
    [17] 罗垚, 陈柏超. 空心矩形截面圆柱线圈自感计算的新方法[J]. 电工技术学报, 2012, 27(6):1-5. (Luo Yao, Chen Baichao. New method for self-inductance calculations of air-core circular coils with rectangular cross-sections[J]. Transactions of China Electrotechnical Society, 2012, 27(6): 1-5
    [18] 罗垚, 陈柏超, 周洪. 有限长圆柱磁屏同轴线圈电感计算方法[J]. 电工技术学报, 2016, 31(14):122-129. (Luo Yao, Chen Baichao, Zhou Hong. Approach for inductance calculations of coaxial circular coils shielded by cylindrical magnetic screen of finite length[J]. Transactions of China Electrotechnical Society, 2016, 31(14): 122-129 doi: 10.3969/j.issn.1000-6753.2016.14.014
  • 期刊类型引用(1)

    1. 陈毅,黄凯,闫杰,李菊香,刘承东. 脉冲电源输出电流精度影响因素分析. 火炮发射与控制学报. 2023(03): 93-100 . 百度学术

    其他类型引用(2)

  • 加载中
图(12) / 表(3)
计量
  • 文章访问数:  963
  • HTML全文浏览量:  398
  • PDF下载量:  112
  • 被引次数: 3
出版历程
  • 收稿日期:  2021-09-23
  • 修回日期:  2021-12-10
  • 录用日期:  2022-01-04
  • 网络出版日期:  2022-01-22
  • 刊出日期:  2022-05-15

目录

/

返回文章
返回