稀薄等离子体稳态流动的数值模拟论文.pdf

中固丁程热物理学会 传热传质学 第1．-一届学术会议 编譬：063265 稀薄等离子体稳态流动的数值模拟宰 王平阳 钱中 杨帆杜朝辉 (．【=海交通人学机械oj动力T程学院，j：海200030) 联系电话：021 E．mail：wangpy@sjtu．edu．an 器(SPT-70)为背景，进行丫数值模拟。输入条件米自神：真空舱内对SPT一70的实验测量，数值模拟结 果．0实验结果进行丫比较，表明本文建讧的数学模型和编制的软件足町信的，可以为电推力器性能评 价及潜41的荆流污染订估提供何效途径。 关键词：Hall推力器等离了体羽流PIC／DSMC方法 1引言 高效率、高比冲、轻重量和长寿命推力器是深空探测和微小卫星姿轨控所急需的先 进推进系统，并已成为航天领域的重要研究方向和热点之一。目前已经发展了电热、电 磁和静电等二人类，十多种不同型式电推力器。其中稳态离子推力器(SPT)‘是目前研 究得最成功、应用最多的电磁式推力器。它具有以下几个特点【o】：(1)结果比离子发动 机简单，可靠性高；(2)无空间电荷效应，推力密度比离子发动机高，体积小；(3)比 冲(约1600s)正好处予目前深空探测和微小卫星所需的比冲范围之内。SPT已经在卫 星轨道维持和姿态控制等方面得到一定程度的应用。 然而，目前SPT推力器的稀薄等离子体羽流对飞行器敏感表面的潜在污染比中性气 体羽流具有更严重的污染问题IlJ，冈为除了中性气体羽流的热污染、光污染、’机械污染 和化学污染外，还有电磁污染。为了有效抑制等离子体羽流的各种污染、进一步改善推 力器性能，有必要对羽流参数进行详细研究，研究途径通常有实验和数值模拟两种，前 者由丁．费Hj人，并且无法获得等离子体羽流的全部信息，冈此，采用实验获得部分信息， 然后作为输入条件，进行数值模拟获得其它信息是一条可行的途径。但由于稀薄等离子 体的复杂特性，采用传统动力论方法、单一粒子模拟方法都不能有效模拟其流动。冈此， 由丁电子质量远小于离子，直接模拟其与各种离子和原子的碰撞显然非常没有效率。 本文利H{实验获得的输入条1：，I：，并假没电子为没有质域但具有动量的流体，再结合 PIC／DSMC方法对稀薄等离子体流动进行数值模拟。模拟结果与实验结果进行了比较。 2仿真模型的建立 2．1计算区域及边界条件 由丁．涉及的问题是二维轴对称的，故只需取一、卜区域作为计算区域，其边界条件如 1．2m， 高O．6m。 Xe+，Xe2+，电子e；由丁． 不考虑溅射的影响，即羽流中包含：中性粒子Xe，离f 高价离子数日比例通常不足1％，故忽略：．二价以上离子。 ·网家臼然科学桀会资助项日(No 1199 (2)同壁面及对对称面的处理 对于同壁面。计算中假定其温度为500K，仿真粒子与壁面接触采用漫反射计算模 型。带电离子，与壁面接触后被中和，为保证等离子体的电中性，假定同时一个电子到 底壁面并且消失。壁面的电势都假定为0。 对称面采用镜反射模型，且假定其径向电场E=0。 ChamberWall＆Er=0 图l计算区域及边界条件 (2)羽流出口(即计算区域入口) SPT羽流计算的出口条件如表所示 表1等离子体出口参数 ．粒子种类Xe Xe+Xe计 出口速度／m·S-1 300 15000 15000