搜索

x

留言板

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

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

本振光功率锁定方法应用于激光外差辐射计的研究

李月 李竣 薛正跃 王晶晶 王贵师 高晓明 谈图

罗杨, 陈茂林, 苏冬冬, 许诺, 王忠晶, 韩志聪, 赵豪. 外磁场作用下的磁等离子体动力学过程仿真. 物理学报, 2022, 71(5): 055204. doi: 10.7498/aps.71.20211383
引用本文: 罗杨, 陈茂林, 苏冬冬, 许诺, 王忠晶, 韩志聪, 赵豪. 外磁场作用下的磁等离子体动力学过程仿真. 物理学报, 2022, 71(5): 055204. doi: 10.7498/aps.71.20211383
Luo Yang, Chen Mao-Lin, Su Dong-Dong, Xu Nuo, Wang Zhong-Jing, Han Zhi-Cong, Zhao Hao. Simulation of magnetoplasmadynamic process with applied magnetic field. Acta Phys. Sin., 2022, 71(5): 055204. doi: 10.7498/aps.71.20211383
Citation: Luo Yang, Chen Mao-Lin, Su Dong-Dong, Xu Nuo, Wang Zhong-Jing, Han Zhi-Cong, Zhao Hao. Simulation of magnetoplasmadynamic process with applied magnetic field. Acta Phys. Sin., 2022, 71(5): 055204. doi: 10.7498/aps.71.20211383

本振光功率锁定方法应用于激光外差辐射计的研究

李月, 李竣, 薛正跃, 王晶晶, 王贵师, 高晓明, 谈图

Research on application of local oscillator power locking method to laser heterodyne radiometer

Li Yue, Li Jun, Xue Zheng-Yue, Wang Jing-Jing, Wang Gui-Shi, Gao Xiao-Ming, Tan Tu
Article Text (iFLYTEK Translation)
PDF
HTML
导出引用
  • 激光外差辐射计具备成本低、体积小、光谱分辨率高等优势, 可扩展现有的地面碳测量网络, 验证卫星观测结果, 并能在卫星观测区域外提供数据覆盖. 本文在现有的激光外差辐射计的基础上, 报道了基于掺铒光纤放大器的可实现本振光功率锁定的近红外激光外差辐射计原型机. 该激光外差辐射计利用一个中心波长为1.603 μm的分布反馈式半导体激光器作为本振光源, 采用掺铒光纤放大器放大本振光功率, 并利用自动功率控制电路实现掺铒光纤放大器输出端光功率的锁定, 消除了由本振光功率变化引起的基线斜率, 从而实现免基线拟合的整层大气透过率谱的测量. 详细评估了基于掺铒光纤放大器的高度集成化的激光外差辐射计的仪器性能, 并在合肥市科学岛(31.9°N, 117.2°E)地区进行了整层大气CO2透过率谱的测量. 在一天的测量时间内得到6组大气CO2透过率谱, 与大气辐射模型模拟结果进行比对, 测量结果一致. 实验结果表明, 掺铒光纤放大器的应用可以提高激光外差辐射计的性能, 优化其结构, 进而为实现无人值守的长期大气CO2浓度观测和构建全面的碳观测网络提供仪器设备的补充.
    Laser heterodyne radiometer has the advantages of low cost, small size, and high spectral resolution. It can expand the existing ground carbon measurement network, verify satellite observation results, and provide data coverage outside the satellite observation area. Using the existing laser heterodyne radiometer, is presented a prototype of near-infrared laser heterodyne radiometer based on the erbium-doped fiber amplifier that can realize local oscillator power locking. In the laser heterodyne radiometer a distributed feedback semiconductor laser with a center wavelength of 1.603 μm is used as a local oscillator light source. The erbium-doped fiber amplifier is used to enhance the local oscillator power, and the automatic power control circuit is adopted to lock the output optical power of the erbium-doped fiber amplifier. The baseline slope caused by the change of the local oscillator power is eliminated, and the whole layer atmospheric transmittance spectrum without baseline fitting is measured. The instrument performance of a highly integrated laser heterodyne radiometer based on an erbium-doped fiber amplifier is evaluated in detail, and the atmospheric CO2 transmittance spectrum is measured in the Science Island (31.9°N, 117.2°E) of Hefei. Six groups of atmospheric CO2 transmittance spectra are obtained during the measurement period of one day, which are compared with the simulation results from an atmospheric radiation model, showing that they are consistent with each other. The experimental results show that the application of erbium-doped fiber amplifier can improve the performance of laser heterodyne radiometer, optimize its structure, and provide equipment supplement for realizing unattended long-term atmospheric CO2 concentration observation and building a comprehensive carbon observation network.
      通信作者: 谈图, tantu@aiofm.ac.cn
    • 基金项目: 国家自然科学基金(批准号: 42075128)和国家重点研发计划(批准号: 2022YFF1300102-1)资助的课题.
      Corresponding author: Tan Tu, tantu@aiofm.ac.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 42075128) and the National Key Research and Development Program of China (Grant No. 2022YFF1300102-1).

    磁等离子体动力学推力器(magnetoplasmadynamic thruster, MPDT)是通过洛伦兹力产生推力的航天推进装置, 也称洛伦兹力加速器, 可实现兆瓦级功率和百牛级推力性能, 是高功率、大推力电推力器的重要发展方向. MPDT在美国、前苏联/俄罗斯、日本、德国、意大利、中国均得到了广泛关注和重点研究[1,2], 实验室内实现了数千瓦至数兆瓦功率的MPDT测试, 最大测试推力达100 N以上.

    研究人员对MPDT进行了大量的理论分析、仿真模拟和实验测试等工作. 在仿真模拟研究方面早期主要采用流体模型: Brushlinskii和Morozov[3,4]最早于20世纪60年代采用流体模型对大功率MPDT开展数值模拟, 获得了MPDT轴向密度分布与电流的关系; Kimura等[5]提出了单温度的二维计算模型, 并使用该模型对MPDT的工作过程进行模拟; Niewood[6]和Sheppard等[7]提出了双温度轴对称数值模型来研究计算推力器内部的流场分布, 考虑了非平衡电离过程和中性气体的影响; Caldo等[8]采用双温度模型来研究了反常运输对MPDT流动的影响. Rudolph[9]在MPDT流体模拟过程中采用多重网格和多重时间步长法, 解决了时间依赖性问题. Heiermann[10]、赵博强等[11]、Sary等[12]、Tahsini[13]采用流体模型对MPDT内部的电离过程和多种输运过程进行了详细建模仿真, 实现了对MPDT推力器的整体性能计算.

    流体模型适用于MPDT的高密度等离子体特性, 具有计算适用性强、计算约束条件少、计算速度快的特点. 但是流体作为基于宏观的平均化方法, 流体模型无法捕捉到等离子体的一些瞬时现象. 粒子云(particle in cell, PIC)[14]等基于动力学过程的等离子体粒子模拟方法是捕捉MPDT工作过程中等离子体演化过程和瞬时现象的优选方法. 由于MPDT高等离子体密度特性对PIC仿真时间步长和空间步长的限制, 使得其对计算资源和计算机时需求巨大, 故早期粒子模拟方法在MPDT仿真中应用较少. Tang等[15-17]开展了PIC方法在MPDT模拟中的应用探索, 采用改变重粒子质量的加速计算方法减少了计算机时, 实现了低功率下MPDT工作过程模拟; Li等[18]采用PIC方法开展了无外加磁场作用下的MPDT工作过程模拟. 本文采用PIC方法结合基于自相似准则的缩比模型进行外加磁场作用下的MPDT工作过程模拟, 解决MPDT极高等离子体密度对模拟时间步长、空间步长的限制, 在减小计算量的同时, 进行真实工况下MPDT工作过程模拟和性能影响因素分析.

    根据磁场产生方式的不同, MPDT分为两大类: 自身场磁等离子体动力学推力器(self-field MPDT, SF-MPDT)和附加场磁等离子体动力学推力器(applied-field MPDT, AF-MPDT). 本文主要开展AF-MPDT工作过程的仿真研究, 其主要结构包括阴极、阳极和电磁线圈, 如图1所示.

    图 1 AF-MPDT结构和工作原理示意图\r\nFig. 1. Schematic of AF-MPDT structure and operating principle.
    图 1  AF-MPDT结构和工作原理示意图
    Fig. 1.  Schematic of AF-MPDT structure and operating principle.

    AF-MPDT 主要有3种加速机制产生推力: 气动加速、自身场加速和附加场加速, 附加场加速又包括霍尔加速和涡旋加速.

    1)气动加速. 气动加速分量是通过电弧欧姆加热将电能转换为热能来提高推进剂的总焓, 进一步喷管将热能转换为动能而产生推力. 当 MPDT 在高功率工况工作时, 电磁加速占主导, 气动加速分量可忽略.

    2)自身场加速. 图2(a)为自身场加速原理. MPDT工作时, 电流轴向分量会产生一个周向自感磁场Bθ, 该磁场又与电流相互作用产生洛伦兹力: 其中径向电流jr与周向自感磁场Bθ, 产生一个轴向的力(jr×Bθ), 该作用力直接对推力产生贡献; 而轴向的电流jz与自感磁场Bθ相互作用产生径向的洛伦兹力(jz×Bθ), 该作用力会导致轴线附近的等离子体密度增大, 并产生一个不平衡的压力, 间接对推力产生贡献. 由于感应磁场Bθ与电流线性相关, 自身场加速产生的推力与电流表现为平方关系. 自身场加速机制在SF-MPDT中贡献较大, 而在AF-MPDT中贡献较小.

    图 2 MPDT加速机理 (a)自身场加速; (b)霍尔加速; (c)涡旋加速\r\nFig. 2. MPDT acceleration mechanism: (a) Self-field acceleration; (b) hall acceleration; (c) vortex acceleration.
    图 2  MPDT加速机理 (a)自身场加速; (b)霍尔加速; (c)涡旋加速
    Fig. 2.  MPDT acceleration mechanism: (a) Self-field acceleration; (b) hall acceleration; (c) vortex acceleration.

    3)霍尔加速. 图2(b)为霍尔加速原理. 在AF-MPDT中, 外部电磁线圈会在阴阳极间产生轴向磁场Bz和径向磁场Br. 根据欧姆定律, 轴向磁场Bz和径向电场作用下, 等离子体会感生出周向电流jθ, 与自身场机制相似, jθ会与附加磁场产生箍缩力(jθ×Bz)和轴向力(jθ×Br).

    4)涡旋加速. 图2(c)为涡旋加速原理. 在AF-MPDT中, 阴阳极间的轴向和径向电流分量, 也会与外部电磁线圈产生的径向磁场Br和轴向磁场Bz相互作用, 产生周向电磁力(jr×Bz)(jz×Br), 使等离子体产生角向涡旋, 宏观上表现为自旋的等离子体径向膨胀, 在物理/磁喷管的作用下再转化为轴向作用力, 产生推力. 涡旋加速产生的推力与电流和外加磁场的乘积IB成正比, Fradkin[19]的实验也验证了该推力关系.

    为了便于仿真结果与实验结果对比分析, 采用Myers团队设计的AF-MPDT推力器[20]作为计算模型, 其几何结构如图3所示, 图中标注单位为mm, 虚线框内为仿真模拟区域.

    图 3 MPDT几何结构与计算区域\r\nFig. 3. MPDT geometry and simulation area.
    图 3  MPDT几何结构与计算区域
    Fig. 3.  MPDT geometry and simulation area.
    2.3.1   PIC方法

    PIC方法是等离子体粒子模拟的常用方法, 其核心过程包括电磁场求解、带电粒子在电磁场作用下运动、以及粒子碰撞三部分[21].

    MPDT工作过程中的电磁场分布及演化可以通过麦克斯韦方程描述, 由于MPDT工作过程中电磁场的时变量较小, 可以忽略, 故可将麦克斯韦方程简化为

    {×B=μ0J,×E=ρε0,
    (1)

    式中, B=×A为磁感应强度, E=ϕ为电场强度, μ0为真空磁导率, J为电流, ρ为电荷密度, ε0为真空介电常数, A为矢量磁位, ϕ为电势. 进一步将(1)式改写为泊松方程形式, 便于数值求解:

    {2A=μ0J,2ϕ=ρε0.
    (2)

    带电粒子的运动通过牛顿-洛伦兹方程进行描述, 如下:

    {mdvdt=e(E+v×B),dsdt=v,
    (3)

    式中, m为粒子质量, v为粒子速度, e为单位电荷, s为粒子位置, t为时间.

    PIC方法处理粒子间的碰撞一般使用蒙特卡罗碰撞(Monte Carlo collision, MCC)方法[22]. 在MPDT模拟过程中, 考虑如下碰撞类型: 电子与中性分子间的弹性碰撞、激发碰撞、电离碰撞, 离子与中性分子间的弹性碰撞和电荷交换(charge exchange, CEX)碰撞, 电子-电子、电子-离子、离子-离子间的库仑碰撞, 以及电子的Bohm碰撞.

    2.3.2   边界条件设置

    仿真中, 主要包含阴极、阳极、绝缘壁面、对称轴和自由边界5种边界类型. 阴极边界设置为Neumann边界, 固定法向场强E=σ/ε0, 切向场强为E//=0, 其中σ为阴极表面电荷密度; 固定阴极的轴向电流密度为jz,cathode=I/(π R2cathode), 径向电流密度和周向电流密度均为0, 其中, I为阴极电流, Rcathode为阴极半径; 阴极作为粒子吸收边界, 粒子碰撞阴极后注销. 考虑离子碰撞阴极壁面的二次电子发射过程, 对于钨阴极, 当氩离子的入射能量介于10—100 eV之间时二次电子发射系数为0.095左右.

    阳极设置为零电势的Drichlet边界和粒子吸收边界, 粒子碰撞阳极后注销. 阴极阳极间的绝缘壁面为Neumann边界, 可认为壁面表面的场强是由累积电荷得到, 且垂直于壁面, 即法向场强E=σ/ε0, 切向场强为E//=0. 对于离子, 绝缘壁面为吸收边界, 离子碰撞绝缘壁面后注销; 对于电子, 分别考虑在壁面沉积、反射、激发一个二次电子、以及激发两个二次电子的情况, 4种情况的概率分别为

    Pdep(ξ)=0.5eξ243.45922,Pref(ξ)=0.5eξ2302,P2se(ξ)=0.5eξ2127.89582,P1se(ξ)=1Pdep(ξ)Pref(ξ)P2se(ξ),
    (4)

    式中, ξ为电子能量, Pdep(ξ)为电子沉积壁面的概率, Pref(ξ)为电子发射反射的概率, P1se(ξ)为电子撞击出一个二次电子的概率, P2se(ξ)为电子撞击出两个二次电子的概率.

    自由边界主要是粒子的出口, 粒子运动越过自由边界后注销; 计算模型是二维轴对称模型, 中心轴线为对称边界, 粒子运动出对称边界时进行镜面反射处理; 电磁场方程均使用Neumann边界为自由边界, 可认为电场、磁场在自由边界处的法相分量为零.

    PIC仿真过程中, 空间步长需要小于德拜长度, 时间步长需要小于等离子体频率, 并同时满足CFL条件. 对于MPDT, 电子数密度一般大于1020 m–3, 对应德拜长度一般小于10–6 m, 单个自由度的仿真网格数一般需要数万至数十万, 等离子体频率约1012 Hz, 对应时间步长约10–12 s. 要同时满足空间步长和时间步长需求, 也会给计算量带来极大挑战, 因此需要一些加速方法在保证仿真准确性的同时, 减少计算量.

    目前对全粒子PIC模拟有两大类加速方法, 一种是直接改变介电常数以及重粒子质量来加速计算的方法[23], 另一种是Khayms[24]和Taccogna等[25]在霍尔推力器仿真中使用基于自相似准则的缩比模型方法, 其基本特征是保证推力器性能参数成比例缩放.

    鉴于MPDT和霍尔推力器工作过程和物理特征相似, 本文采用基于自相似准则的缩比模型方法, 结合PIC粒子模拟方法, 进行MPDT仿真模型构建和加速计算. 为了获得准确的推力器性能参数, 缩比模型中, 相似准则选取如下: 1)推力器的比冲Isp不变; 2)磁场对电子的约束效果不变; 3)电子平均自由程与推力器特征尺寸比值不变.

    缩比模型中, 将推力器特征尺寸缩小至1/τ, 则物理量电流I、质量流量˙m、数密度ne, ni、电场E、磁场B、推力F、比冲Isp和效率η的缩比系数分别为1/τ, 1/τ, 1/τ, τ, τ, τ, 1, 1和1/τ. 其中自变量为计算输入时需要调整的量, 因变量表示根据调整后的输入量而发生改变的物理量. 缩比模型将推力器尺寸缩小为原来的1/τ后, 电子数密度同样也扩大为原来的τ倍, 使得德拜长度缩小1/τ1/2, 可以认为缩比系数τ对单方向的网格数量的缩小量为1/τ1/2, 对于轴对称模型, 网格数量缩小1/τ. 本文仿真中τ取为8 × 106.

    中性气体流场分布一般采用直接蒙特卡罗模拟(direct simulation Monte Carlo, DSMC)方法求解, 将等离子体放电和输运过程模拟的PIC方法与中性气体流场模拟的DSMC方法耦合, 会大幅增大计算需求. 本文将中性气体流场与等离子体流场解耦处理, 即先单独采用DSMC方法模拟MPDT中的中性气体流场, 然后采用电推进仿真中常用的PIC/MCC方法处理带电粒子与中性气体原子间的相互作用, 即将具备一定分布的气体原子作为背景, 参与带电粒子的碰撞处理, 但其分布不随模拟过程变化, 这样可在保证较好的计算精度的同时, 大幅缩减计算机时.

    针对Ar工质, 流量为˙m = 0.1 g/s的情况下的气体初始流场, Ar原子数密度分布和流线如图4所示, 流场内中性粒子数密度量级约在1020— 1021 m–3量级.

    图 4 Ar原子数密度分布和流线\r\nFig. 4. The number density and streamlines of Ar.
    图 4  Ar原子数密度分布和流线
    Fig. 4.  The number density and streamlines of Ar.

    通过推力器内部等离子体特性参数的演化, 可以分析推力器内部放电和带电粒子的输运过程. 针对外磁场B = 0.03 T, 阴极电流I = 1000 A工作参数下的MPDT的工作过程进行仿真, 在缩比模型中, 计算时间步长为1 × 10–18 s, 对应实际情况的时间步长8 × 10–12 s.

    3.2.1   电子数密度

    MPDT的放电由阴极发射的电子电流主导, 故首先分析放电过程中的电子数密度分布随时间的演化过程, 不同时间步的电子数密度分布仿真结果如图5所示.

    图 5 电子数密度分布\r\nFig. 5. Distribution of electron number density.
    图 5  电子数密度分布
    Fig. 5.  Distribution of electron number density.

    MPDT点火启动过程中, 电子数密度分布随时间的变化可初步分为3个阶段: 放电初期, 电子从阴极发射后, 受外加磁场约束, 沿磁力线螺旋运动, 呈扩张型分布; 放电中期, 电子在阴极尖端下游堆积, 逐渐形成一个高密度分布区域, 并向阳极方向延伸, 构成一条连通阴阳极的高密度电子分布带, 形成放电电弧; 放电末期, 放电电弧逐渐稳定, 等离子体区域进一步扩张, 并在阴阳极下游形成完整的高密度的稳态分布, 推力器进入稳定工作状态. 推力器在达稳态工作状态后, 在推力器外部, 电子束流分为两个部分, 一个集中在阴极尖端的轴线处向下游扩散, 另一部分电子随磁力线向下游扩散, 这和实验观测到的现象一致[1]. 稳态时, MPDT阳极尖端的电子数密度达1021 m–3, 阴阳极下游的电子数密度达1020 m–3.

    3.2.2   离子数密度

    阴极发射电子在电磁场共同作用下运动的同时, 会碰撞中性气体分子发生电离, 产生新的电子和离子. 由于离子完全由电离过程产生, 离子数密度分布随时间的演化过程更真实地反映MPDT的电离过程, 不同时间步的离子数密度分布仿真结果如图6所示.

    图 6 离子数密度分布\r\nFig. 6. Distribution of ion number density
    图 6  离子数密度分布
    Fig. 6.  Distribution of ion number density

    图6可知, 在放电初期, 离子数密度同样由小增大, 并随着放电进行, 其分布也逐渐稳定. 在放电初期, 离子数密度的分布与电子分布相似, 其原因是离子由电子电离碰撞产生的, 在高电子数密度处离子数密度高; 随着放电进行, 离子数密度的分布逐渐偏离电子分布, 呈现出高密度区域逐渐由阴极尖端附近向整个下游区域扩张并向轴线附近聚焦的趋势, 离子分布偏离电子分布的原因是离子质量大, 受磁场的约束远小于电子, 故运动轨迹和最终分布也与电子不同. 稳态工作时, 离子数密度最高达到1021 m–3, 主要仿真区域离子数密度约1020 m–3, 与电子数密度相当.

    外加磁场是影响MPDT推力器性能的重要参数, 在阴极电流I = 1000 A, 流量为˙m = 0.1 g/s工况下, 开展外磁场B = 0.008, 0.017, 0.025, 0.034, 0.052和0.069 T情况下的MPDT工作过程仿真, 分析不同磁场强度下的电子数密度分布和推力器性能参数, 并进行仿真结果与实验结果对比, 验证仿真模型和方法的有效性.

    3.3.1   不同外磁场的电子数密度分布

    不同外磁场作用下, MPDT稳态工作时的电子数密度分布如图7所示.

    图 7 不同磁场强度下的电子数密度分布\r\nFig. 7. Electron number density distribution with different applied magnetic field.
    图 7  不同磁场强度下的电子数密度分布
    Fig. 7.  Electron number density distribution with different applied magnetic field.

    图7可知, 在B = 0.008 T情况下, 电子仅在阴极下游的轴线附近形成高密度区域, 未能形成高效的阴阳极电弧放电; 随着磁场强度的增大, 中轴线高密度区域逐渐向径向移动, 在外磁场强度B = 0.017—0.340 T之间时, 电子数密度出现两个高密度区域, 分别位于轴线附近和推力器出口下游区域, 构成稳定的放电电弧; 外磁场强度继续增大, 轴线附近的高密度区域逐渐缩小, 当磁场B > 0.052 T时, 中心轴线的高密度区域基本不存在. 分析认为: 在外磁场强度较小时, 推力器的自感磁场比重较大, 自感磁场主要为周向磁场, 约束电子在轴线附近聚集; 随着外磁场的增加, 自感磁场的比重降低, 外磁场对电子的约束成为主要部分, 外磁场将使得电子沿着磁力线向下游运动.

    3.3.2   外磁场对推力器性能的影响

    根据不同外磁场作用下的仿真结果, 计算了MPDT稳态工作时的放电电压U、推力F和效率η, 如图8所示.

    图 8 MPDT性能随外磁场的变化\r\nFig. 8. MPDT performance changes with applied magnetic field.
    图 8  MPDT性能随外磁场的变化
    Fig. 8.  MPDT performance changes with applied magnetic field.

    图8可知, 当B < 0.053 T时, 放电电压随外磁场增大而降低, B ≥ 0.052 T时放电电压逐渐变大. 其原因是阴极表面电荷密度决定了阴极电势即放电电压的大小, 而阴极表面电子的运动受磁场影响严重: 当外磁场较小时, 自感磁场对电子的约束占据主导, 电子向阴极聚集, 增大了阴极表面电荷密度, 放电电压升高; 随着外加磁场增大, 电子逐渐跳出自感磁场约束, 远离阴极, 使得放电电压降低; 当外磁场进一步增大时, 外加磁场对电子约束增强, 电子难以跨越阴极磁力线沿径向扩散, 使得阴极表面电荷密度增大, 放电电压再次升高.

    不同磁感应强度条件下的推力结果表现为随磁感应强度B近似线性增大, 与Myers[20]的实验结论一致. 在B = 0.069 T情况下的推力值为0.98 N, 略小于Myers[20]的实验结果1.16 N, 误差产生的原因可能是本文仿真忽略了气动推力的贡献; 推力器的效率同样表现为随磁感应强度B增大而增大的趋势.

    3.3.3   与实验测试数据对比分析

    Myers[20]实验测量了不同外加磁场情况下的MPDT推力器的比冲、效率等性能参数, 为了进一步验证本文仿真程序的准确性, 将仿真结果与Myers[20]的实验结果进行对比, 如图9所示.

    图 9 仿真与实验结果对比\r\nFig. 9. Comparison of simulation and experimental results.
    图 9  仿真与实验结果对比
    Fig. 9.  Comparison of simulation and experimental results.

    图9可知, 仿真获得的比冲、效率等性能参数随外加磁场强度的变化趋势与实验结果一致. 计算比冲值低于实验结果约20%, 其原因是计算比冲由Isp=F/˙m计算得来, 在本文仿真中未考虑气动力对总推力的贡献, 使得计算推力性能低于真实情况, 并导致计算比冲低于实验结果. 计算效率值相比实验结果有约50%的误差, 其原因可能是计算效率与推力和比冲的乘积成正比, 计算推力和计算比冲均低于实验结果, 导致计算效率误差进一步增大.

    虽然推力、比冲和效率的仿真结果因忽略气动加速贡献略小于真实值, 但误差不是很大, 并且能准确反映工况对性能的影响, 故认为, 本文建立的模型和方法在高磁场强度工况下有较好的准确性和可信度.

    阴极电流是影响MPDT推力器性能的另一个重要参数, 在外磁场B = 0.025 T, 流量˙m = 0.1 g/s工况下, 开展阴极电流I = 750, 1000, 1250, 1500, 1750和2000 A情况下的MPDT工作过程仿真, 分析阴极电流对电子数密度分布和推力器性能的影响.

    3.4.1   电子数密度分布

    不同阴极电流情况下, MPDT稳态工作时的电子数密度分布如图10所示.

    图 10 不同阴极电流下的电子数密度分布\r\nFig. 10. Distribution of electron number density under different cathode current
    图 10  不同阴极电流下的电子数密度分布
    Fig. 10.  Distribution of electron number density under different cathode current

    图10可知, 推力器内部的电子数密度随着电流的增大而增大; 对于B = 0.025 T, 流量˙m = 0.1 g/s工况, 750—2000 A的阴极电流情况, 都在阴极尖端下游位置都出现了电子高密度分布区域. 阴极电流增大, 高密度区域也逐渐扩大, 具体表现为随着电流增大, 外磁场对放电等离子体束的约束效果逐渐减弱, 低电流工况的等离子体束流发散角明显小于高电流工况, 其原因是电流增大, 电子数密度增大, 电子间的库仑碰撞增多, 会使得更多的电子沿跨越磁力线沿径向扩散.

    3.4.2   阴极电流对推力器性能的影响

    根据不同阴极电流条件下的仿真结果, 计算了MPDT稳态工作时的放电电压U、推力F和效率η, 如图11所示.

    图 11 MPDT性能随阴极电流的变化\r\nFig. 11. MPDT performance changes with cathode current.
    图 11  MPDT性能随阴极电流的变化
    Fig. 11.  MPDT performance changes with cathode current.

    图11可知, 放电电压随阴极电流的增大而增大, 原因是阴极电流越大, 阴极表面的电子数密度越大, 阴极表面电荷增加, 放电电压增大; 另外, 自感磁场随阴极电流增大而增大, 并进一步约束电子向阴极聚集, 使得阴极表面电荷增加, 放电电压增大.

    在750—2000 A阴极电流范围内, 推力随电流增大而近似线性增大. 在I = 750—2000 A的电流工况下, 计算区域内的自感磁场强度为0.001—0.040 T, 量级上与外磁场相当, 但自感磁场主要集中于阴极附近, 在其余大部分区域外磁场明显强于自感磁场, 故此时MPDT的推力由涡旋加速机制产生, 正比于电流大小.

    推力器的效率同样表现为随阴极电流I增加而增大的趋势, 这与MPDT低功率工况效率低, 高功率工况效率高的实验结果一致.

    本文建立了基于缩比模型的MPDT工作过程的粒子仿真模型, 研究了MPDT点火启动时的放电电弧建立过程, 以及外磁场强度和阴极电流对MPDT性能参数的影响, 并通过仿真结果与已有实验数据对比, 验证了仿真模型的准确性. 点火启动过程研究表明, 放电电弧建立后, 在阴极下游和阳极附近形成两处高电子数密度区并连通以维持稳定的放电电弧, 离子主要集中在阴极下游并沿轴向喷射形成等离子体喷流. 不同外磁场强度下的结果表明, 提高外磁场强度有利于构建阴阳极放电电弧, 随着外磁场强度由0.008 T提升至0.069 T, MPDT的推力产生机制逐渐由自身场加速转变为涡旋加速, 推力和比冲性能随外磁场强度增加近似线性增大. 不同阴极电流的结果表明, 在B = 0.025 T, I = 750—2000 A条件下, 在推力器内部绝大部分区域, 外磁场强度大于自感磁场, 推力产生机制为涡旋加速, 推力随电流增加而近似线性增加. MPDT效率总体表现为随外磁场强度和阴极电流增大而增大的趋势, 符合其适用于高功率工况的特性.

    [1]

    查玲玲, 王薇, 谢宇, 单昌功, 曾祥昱, 孙友文, 殷昊, 胡启后 2022 光谱学与光谱分析 42 1036Google Scholar

    Cha L L, Wang W, Xie Y, Shan C G, Zeng X Y, Sun Y W, Yin H, Hu Q H 2022 Spectrosc. Spect. Anal. 42 1036Google Scholar

    [2]

    Wang J, Sun C, Wang G, Zou M, Tan T, Liu K, Chen W, Gao X 2020 Opt. Lasers. Eng. 129 106083Google Scholar

    [3]

    Wang J, Wang G, Tan T, Zhu G, Sun C, Cao Z, Chen W, Gao X 2019 Opt. Express. 27 9610Google Scholar

    [4]

    李竣, 薛正跃, 刘笑海, 王晶晶, 王贵师, 刘锟, 高晓明, 谈图 2022 物理学报 71 074204Google Scholar

    Li J, Xue Z Y, Liu X H, Wang J J, Wang G S, Liu K, Gao X M, Tan T 2022 Acta Phys. Sin. 71 074204Google Scholar

    [5]

    Lu X, Huang Y, Wu P, Liu D, Ma H, Wang G, Cao Z 2022 Remote Sens. 14 1489Google Scholar

    [6]

    薛正跃, 李竣, 刘笑海, 王晶晶, 高晓明, 谈图 2021 物理学报 70 217801Google Scholar

    Xue Z Y, Li J, Liu X H, Wang J J, Gao X M, Tan T 2021 Acta Phys. Sin. 70 217801Google Scholar

    [7]

    Deng H, Li R, Liu H, He Y, Yang C, Li X, Xu Z, Kan R 2022 Opt. Lett. 47 4335Google Scholar

    [8]

    Sappey A, Masterson B, Howell J 2021 Appl. Opt. 61 2697Google Scholar

    [9]

    卢兴吉, 曹振松, 谈图, 黄印博, 高晓明, 饶瑞中 2019 物理学报 68 064208Google Scholar

    Lu X J, Cao Z S, Tan T, Huang Y B, Gao X M, Rao R Z 2019 Acta Phys. Sin. 68 064208Google Scholar

    [10]

    Xue Z, Shen F, Li J, Liu X, Wang J, Wang G, Chen W, Gao X, Tan T 2022 Opt. Express 30 31828Google Scholar

    [11]

    孙春艳, 王贵师, 朱公栋, 谈图, 刘锟, 高晓明 2020 物理学报 69 144201Google Scholar

    Sun C Y, Wang G S, Zhu G D, Tan T, Liu K, Gao X M 2020 Acta Phys. Sin. 69 144201Google Scholar

    [12]

    Clarke G B, Wilson E L, Miller J H, Melroy H R 2014 Meas. Sci. Technol. 25 055204Google Scholar

    [13]

    Wilson E L, DiGregorio A J, Riot V J, Ammons M S, Bruner W W, Carter D, Mao J, Ramanathan A, Strahan S E, Oman L D, Hoffman C, Garner R M 2017 Meas. Sci. Technol. 28 035902Google Scholar

    [14]

    Wilson E L, DiGregorio A J, Villanueva G, Grunberg C E, Souders Z, Miletti K M, Menendez A, Grunberg M H, Floyd M A M, Bleacher J E, Euskirchen E S, Edgar C, Caldwell B J, Shiro B, Binsted K 2019 Appl. Phys. B 125 211Google Scholar

    [15]

    Zenevich S, Gazizov I, Churbanov D, Plyashkov Y, Spiridonov M, Talipov R, Rodin A 2021 Remote Sens. 13 2235Google Scholar

    [16]

    王晶晶 2021 博士学位论文 (合肥: 中国科学技术大学)

    Wang J J 2021 Ph. D. Dissertation (Hefei: University of Science and Technology of China) (in Chinese)

    [17]

    Deng H, Yang C, Xu Z, Li M, Huang A, Yao L, Hu M, Chen B, He Y, Kan R, Liu J 2021 Opt. Express 29 2003Google Scholar

    [18]

    Huang J, Huang Y, Lu X, Liu D, Yuan Z, Qi G, Cao Z 2022 Front. Phys. 10 835189Google Scholar

    [19]

    Hoffmann A, Huebner M, Macleod N, Weidmann D 2018 Opt. Lett. 43 3810Google Scholar

    [20]

    Parvitte B, Zéninari V, Thiébeaux C, Delahaigue A, Courtois D 2004 Spectrochim. Acta, Part A 60 1193Google Scholar

    [21]

    Nakazawa M 2014 Opt. Rev. 21 862Google Scholar

    [22]

    Shen F, Wang G, Wang J, Tan T, Wang G, Jeseck P, Te Y, Gao X, Chen W 2021 Opt. Lett. 46 3171Google Scholar

    期刊类型引用(3)

    1. 罗杨,陈茂林,许诺,高辉,叶剑民,韩志聪. 场致发射电推力器羽流分布特性仿真. 推进技术. 2023(06): 201-209 . 百度学术
    2. 吴建军,欧阳,张宇,赵元政,李宇奇. 激光-电磁复合推力器研究现状与展望. 推进技术. 2023(06): 58-77 . 百度学术
    3. 李平,彭超,陶华杨,张博文,陈思乐,陈兆权. 磁场增强DBD等离子体改善煤炭燃烧的特性. 高电压技术. 2023(10): 4436-4444 . 百度学术

    其他类型引用(4)

  • 图 1  激光外差辐射计结构示意图. DL, 半导体激光器; EDFA, 掺铒光纤放大器; DFB, 分布反馈式半导体激光器; PD, 光电探测器; Amp, 放大器; Schottky diode, 肖特基二极管

    Fig. 1.  Structure diagram of laser heterodyne radiometer. DL, diode lasers; EDFA, erbium-doped fiber amplifier; DFB, distributed feedback; PD, photodetector; Amp, Amplifier.

    图 2  不同输出电压下激光波数

    Fig. 2.  Laser wavenumber under different output voltages.

    图 3  EDFA原理图、实物图及激光器输出功率变化图 (a) EDFA原理图(ISO, 光隔离器; EDF, 掺铒光纤; WDM, 波分复用器); (b) EDFA实物图; (c) 激光功率(红色实线), 经自动功率控制锁定后的EDFA输出功率(黑色实线)

    Fig. 3.  Schematic diagram and physical diagram of EDFA, and power variation diagram of laser: (a) The schematic diagram of EDFA (ISO, isolator; EDF, erbium-doped fiber; WDM, wavelength division multiplexer); (b) the physical diagram of EDFA; (c) laser power (red solid line), output power of EDFA locked by automatic power control (black solid line).

    图 4  激光外差辐射计实物图

    Fig. 4.  Physical picture of laser heterodyne radiometer.

    图 5  信号功率谱

    Fig. 5.  Signal power spectrum.

    图 6  激光波数及稳定性分析 (a) 激光波数实际值(黑色曲线)与设定值(红色曲线); (b) 激光波数稳定性

    Fig. 6.  Laser wavenumber and stability analysis: (a) The actual value of laser wavenumber (black curve) and the set value (red curve); (b) laser wavenumber stability.

    图 7  外差信号原始数据及艾伦方差曲线函数 (a) 外差信号原始数据; (b) 艾伦方差曲线函数

    Fig. 7.  Original data of heterodyne signal and Allen variance curve function: (a) The original data of heterodyne signal; (b) Allen variance curve function

    图 8  仪器装置的信噪比测量

    Fig. 8.  SNR measurement of instrument device.

    图 9  使用EDFA前(a)后(b)外差信号及基线变化

    Fig. 9.  Heterodyne signal and baseline changes before (a) and after (b) using EDFA.

    图 10  测量期间信号光功率

    Fig. 10.  Power of signal during measurement.

    图 11  实测大气CO2透过率谱与大气辐射模型模拟结果比较 (a) 大气辐射模型模拟结果; (b)实测大气CO2透过率谱

    Fig. 11.  Comparison between the measured atmospheric transmittance spectrum of CO2 and the simulation results of atmospheric radiation model: (a) The simulation results of atmospheric radiation model; (b) the measured atmospheric transmittance spectrum of CO2.

  • [1]

    查玲玲, 王薇, 谢宇, 单昌功, 曾祥昱, 孙友文, 殷昊, 胡启后 2022 光谱学与光谱分析 42 1036Google Scholar

    Cha L L, Wang W, Xie Y, Shan C G, Zeng X Y, Sun Y W, Yin H, Hu Q H 2022 Spectrosc. Spect. Anal. 42 1036Google Scholar

    [2]

    Wang J, Sun C, Wang G, Zou M, Tan T, Liu K, Chen W, Gao X 2020 Opt. Lasers. Eng. 129 106083Google Scholar

    [3]

    Wang J, Wang G, Tan T, Zhu G, Sun C, Cao Z, Chen W, Gao X 2019 Opt. Express. 27 9610Google Scholar

    [4]

    李竣, 薛正跃, 刘笑海, 王晶晶, 王贵师, 刘锟, 高晓明, 谈图 2022 物理学报 71 074204Google Scholar

    Li J, Xue Z Y, Liu X H, Wang J J, Wang G S, Liu K, Gao X M, Tan T 2022 Acta Phys. Sin. 71 074204Google Scholar

    [5]

    Lu X, Huang Y, Wu P, Liu D, Ma H, Wang G, Cao Z 2022 Remote Sens. 14 1489Google Scholar

    [6]

    薛正跃, 李竣, 刘笑海, 王晶晶, 高晓明, 谈图 2021 物理学报 70 217801Google Scholar

    Xue Z Y, Li J, Liu X H, Wang J J, Gao X M, Tan T 2021 Acta Phys. Sin. 70 217801Google Scholar

    [7]

    Deng H, Li R, Liu H, He Y, Yang C, Li X, Xu Z, Kan R 2022 Opt. Lett. 47 4335Google Scholar

    [8]

    Sappey A, Masterson B, Howell J 2021 Appl. Opt. 61 2697Google Scholar

    [9]

    卢兴吉, 曹振松, 谈图, 黄印博, 高晓明, 饶瑞中 2019 物理学报 68 064208Google Scholar

    Lu X J, Cao Z S, Tan T, Huang Y B, Gao X M, Rao R Z 2019 Acta Phys. Sin. 68 064208Google Scholar

    [10]

    Xue Z, Shen F, Li J, Liu X, Wang J, Wang G, Chen W, Gao X, Tan T 2022 Opt. Express 30 31828Google Scholar

    [11]

    孙春艳, 王贵师, 朱公栋, 谈图, 刘锟, 高晓明 2020 物理学报 69 144201Google Scholar

    Sun C Y, Wang G S, Zhu G D, Tan T, Liu K, Gao X M 2020 Acta Phys. Sin. 69 144201Google Scholar

    [12]

    Clarke G B, Wilson E L, Miller J H, Melroy H R 2014 Meas. Sci. Technol. 25 055204Google Scholar

    [13]

    Wilson E L, DiGregorio A J, Riot V J, Ammons M S, Bruner W W, Carter D, Mao J, Ramanathan A, Strahan S E, Oman L D, Hoffman C, Garner R M 2017 Meas. Sci. Technol. 28 035902Google Scholar

    [14]

    Wilson E L, DiGregorio A J, Villanueva G, Grunberg C E, Souders Z, Miletti K M, Menendez A, Grunberg M H, Floyd M A M, Bleacher J E, Euskirchen E S, Edgar C, Caldwell B J, Shiro B, Binsted K 2019 Appl. Phys. B 125 211Google Scholar

    [15]

    Zenevich S, Gazizov I, Churbanov D, Plyashkov Y, Spiridonov M, Talipov R, Rodin A 2021 Remote Sens. 13 2235Google Scholar

    [16]

    王晶晶 2021 博士学位论文 (合肥: 中国科学技术大学)

    Wang J J 2021 Ph. D. Dissertation (Hefei: University of Science and Technology of China) (in Chinese)

    [17]

    Deng H, Yang C, Xu Z, Li M, Huang A, Yao L, Hu M, Chen B, He Y, Kan R, Liu J 2021 Opt. Express 29 2003Google Scholar

    [18]

    Huang J, Huang Y, Lu X, Liu D, Yuan Z, Qi G, Cao Z 2022 Front. Phys. 10 835189Google Scholar

    [19]

    Hoffmann A, Huebner M, Macleod N, Weidmann D 2018 Opt. Lett. 43 3810Google Scholar

    [20]

    Parvitte B, Zéninari V, Thiébeaux C, Delahaigue A, Courtois D 2004 Spectrochim. Acta, Part A 60 1193Google Scholar

    [21]

    Nakazawa M 2014 Opt. Rev. 21 862Google Scholar

    [22]

    Shen F, Wang G, Wang J, Tan T, Wang G, Jeseck P, Te Y, Gao X, Chen W 2021 Opt. Lett. 46 3171Google Scholar

  • [1] 沈晓阳, 成一灏, 夏林. 紧凑型冷原子高分辨成像系统光学设计. 物理学报, 2024, 73(6): 066701. doi: 10.7498/aps.73.20231689
    [2] 黄一帆, 邢阳光, 沈文杰, 彭吉龙, 代树武, 王颖, 段紫雯, 闫雷, 刘越, 李林. 亚角秒空间分辨的太阳极紫外宽波段成像光谱仪光学设计. 物理学报, 2024, 73(3): 039501. doi: 10.7498/aps.73.20231481
    [3] 吴长茂, 唐熊忻, 夏媛媛, 杨瀚翔, 徐帆江. 用于空间相机设计的高精度光线追迹方法. 物理学报, 2023, 72(8): 084201. doi: 10.7498/aps.72.20222463
    [4] 侯晨阳, 孟凡超, 赵一鸣, 丁进敏, 赵小艇, 刘鸿维, 王鑫, 娄淑琴, 盛新志, 梁生. “机器微纳光学科学家”: 人工智能在微纳光学设计的应用与发展. 物理学报, 2023, 72(11): 114204. doi: 10.7498/aps.72.20230208
    [5] 李竣, 薛正跃, 刘笑海, 王晶晶, 王贵师, 刘锟, 高晓明, 谈图. 激光外差光谱仪模拟风场探测. 物理学报, 2022, 71(7): 074204. doi: 10.7498/aps.71.20211252
    [6] 邱乙耕, 范元媛, 颜博霞, 王延伟, 吴一航, 韩哲, 亓岩, 鲁平. 光声光谱仪用三维扩展光源光场整形系统设计与实验. 物理学报, 2021, 70(20): 204201. doi: 10.7498/aps.70.20210691
    [7] 薛正跃, 李竣, 刘笑海, 王晶晶, 高晓明, 谈图. 基于激光外差探测的大气N2O吸收光谱测量与廓线反演. 物理学报, 2021, 70(21): 217801. doi: 10.7498/aps.70.20210710
    [8] 许祥馨, 常军, 武楚晗, 宋大林. 基于双随机相位编码的局部混合光学加密系统. 物理学报, 2020, 69(20): 204201. doi: 10.7498/aps.69.20200478
    [9] 冯帅, 常军, 胡瑶瑶, 吴昊, 刘鑫. 偏振成像激光雷达与短波红外复合光学接收系统设计与分析. 物理学报, 2020, 69(24): 244202. doi: 10.7498/aps.69.20200920
    [10] 孙春艳, 王贵师, 朱公栋, 谈图, 刘锟, 高晓明. 基于高分辨率激光外差光谱反演大气CO2柱浓度及系统测量误差评估方法. 物理学报, 2020, 69(14): 144201. doi: 10.7498/aps.69.20200125
    [11] 徐平, 杨伟, 张旭琳, 罗统政, 黄燕燕. 集成化导光板下表面微棱镜二维分布设计. 物理学报, 2019, 68(3): 038502. doi: 10.7498/aps.68.20181684
    [12] 操超, 廖志远, 白瑜, 范真节, 廖胜. 基于矢量像差理论的离轴反射光学系统初始结构设计. 物理学报, 2019, 68(13): 134201. doi: 10.7498/aps.68.20190299
    [13] 刘飞, 魏雅喆, 韩平丽, 刘佳维, 邵晓鹏. 基于共心球透镜的多尺度广域高分辨率计算成像系统设计. 物理学报, 2019, 68(8): 084201. doi: 10.7498/aps.68.20182229
    [14] 卢兴吉, 曹振松, 谈图, 黄印博, 高晓明, 饶瑞中. 激光外差光谱仪的仪器线型函数研究. 物理学报, 2019, 68(6): 064208. doi: 10.7498/aps.68.20181620
    [15] 沈本兰, 常军, 王希, 牛亚军, 冯树龙. 三反射主动变焦系统设计. 物理学报, 2014, 63(14): 144201. doi: 10.7498/aps.63.144201
    [16] 裴琳琳, 吕群波, 王建威, 刘扬阳. 编码孔径成像光谱仪光学系统设计. 物理学报, 2014, 63(21): 210702. doi: 10.7498/aps.63.210702
    [17] 任洪亮. 有限远共轭显微镜光镊设计和误差分析. 物理学报, 2013, 62(10): 100701. doi: 10.7498/aps.62.100701
    [18] 邱巍, 马英驰, 吕品, 刘典, 徐晓娟, 张程华. 室温掺铒光纤放大器中实现参量控制无损耗光速减慢传输. 物理学报, 2012, 61(9): 094204. doi: 10.7498/aps.61.094204
    [19] 王 方, 朱启华, 蒋东镔, 张清泉, 邓 武, 景 峰. 多程放大系统主放大级光学优化设计. 物理学报, 2006, 55(10): 5277-5282. doi: 10.7498/aps.55.5277
    [20] 孙 强, 于 斌, 王肇圻, 母国光, 卢振武. 谐衍射双波段红外超光谱探测系统研究. 物理学报, 2004, 53(3): 756-761. doi: 10.7498/aps.53.756
  • 期刊类型引用(3)

    1. 罗杨,陈茂林,许诺,高辉,叶剑民,韩志聪. 场致发射电推力器羽流分布特性仿真. 推进技术. 2023(06): 201-209 . 百度学术
    2. 吴建军,欧阳,张宇,赵元政,李宇奇. 激光-电磁复合推力器研究现状与展望. 推进技术. 2023(06): 58-77 . 百度学术
    3. 李平,彭超,陶华杨,张博文,陈思乐,陈兆权. 磁场增强DBD等离子体改善煤炭燃烧的特性. 高电压技术. 2023(10): 4436-4444 . 百度学术

    其他类型引用(4)

计量
  • 文章访问数:  4625
  • PDF下载量:  62
  • 被引次数: 7
出版历程
  • 收稿日期:  2023-02-22
  • 修回日期:  2023-03-13
  • 上网日期:  2023-03-21
  • 刊出日期:  2023-05-05

/

返回文章
返回