电推力器等离子体羽流三维非稳态流动模拟.pdfVIP

中国工程热物理学会 传热传质学 学术会议论文 编号：083219 电推力器等离子体羽流电推力器等离子体羽流三维三维非稳态非稳态 电推力器等离子体羽流电推力器等离子体羽流三维三维非稳态非稳态 ** ** 流动模拟流动模拟 流动模拟流动模拟 王平阳 钱中 严立 杜朝辉 (上海交通大学机械与动力工程学院，上海 200030) (Tel: 021Email: wangpy@sjtu.edu.cn) 摘要：以美国格林研究中心设计的脉冲等离子体推力器(PPT)工作时产生的羽流为研究对象，建立三维、 非稳态 PICDSMC 混合方法模型，模拟该羽流的瞬时特性。通过与国外已有实验数据和二维仿真结果 的对比表明，所建模型和所编制的计算机代码基本正确，三维模拟结果在近场明显优于二维模拟，能 更好地展现 PPT 的工作特性以及与实验更广泛地进行比较。本文的研究方法、结果对等离子体羽流的 更深入研究具有参考价值。 关键词：三维 等离子体 羽流 PICDSMC 混合方法 0 前言 脉冲等离子体推力器（Pulsed plasma thruster ，PPT ）具有推力小 （N 级）、控制精 度高、结构简单、易于制造等特点，特别适合微小卫星的姿轨控和深空探测，具有重要 的经济价值。然而，PPT 羽流属于稀薄等离子体超音速射流，其流动情况十分复杂，流 场特性至今未能得到准确测定。美国伍斯特工学院 Gatsonis 等[13]针对格林研究中心的 PPT 样机，运用静电探针法，初步获得了羽流参数分布。同时，在等效面积的基础上， 将三维流场简化为二维，采用粒子单元 （PIC ）法和直接模拟蒙特卡罗 （DSMC ）法对 流场进行了模拟。三维非稳态的等离子体超音速流场的研究至今未见文献报道。 本文将在国外已有研究的基础上，首次建立与推力器结构相符的三维仿真模型，运 用 PICDSMC 法，模拟该 PPT 羽流的流动，并通过与已有实验和二维仿真结果的对比， 验证仿真模型。所得研究结果可为同类问题的研究提供参考。 1 仿真建模 1.11.1 推力器模型 11.1.1 PPT 以脉冲方式工作，采用气体、汞、固态聚四氟 乙烯或聚乙烯作为工质，工质置于两个电极之间，通过 贮能电容器向工质放电产生高温热流使其表面层烧蚀、 离解，形成等离子体。在很薄的电离层内，通过自身磁 压和气体动压的共同作用，加速等离子体并将其喷出而 产生推力。本文研究的PPT 模型由 NASA GRC 设计， 如图 1 所示。它由两个平行的铜质电极 （阴、阳极）、 Teflon 药柱 （由弹簧控制供给）、电容等组成。 图 1 推力器实物模型 *国家自然科学基金资助项目 （No） 1.21.2 电子动量方程及其简化 11.2.2 完整羽流电子动量方程如下 ∂u mn e + mn u ⋅ ∇u = −ne E + u × B − ∇p − mn ν u − u e e e e e