基于微机电系统技术的固体微推力器阵列可以为微型卫星提供位置保持、姿态控制、阻力补偿以及轨道调整的冲量和冲量矩,并且有着无可动部件、结构简单、成本低、可靠性高的优点,因此从20世纪90年代末开始受到了广泛的重视。发展固体推进剂微推力器的关键是产生极其精确的推力和元冲量。本文主要研究的是一种固体微推力器的内流场,首先建立了固体微推力器工作过程的数学模型,包括燃面退移模型、壁面滑移模型、离散相颗粒模型以及流固耦合传热模型等。之后建立了二维轴对称几何模型,结合FLUENT软件及用户自定义函数搭建了固体微推力器内流场的数值仿真平台。目前对固体微推力器的数值仿真研究中,较少的是对其瞬态工作过程的研究,故本文采用燃速公式的方法模拟固体推进剂燃面的瞬态退移过程,推进剂采用的是高氯酸铵/聚羟基丁二烯(AP/HTPB)复合推进剂。基于数值仿真平台,对固体微推力器的完整工作过程进行了瞬态仿真,得到了工作过程中燃烧室内温度、压力、速度、颗粒相浓度等参数的分布情况及随时间的变化情况,并计算得到推力器的总冲、比冲等参数,以此来分析固体微推力器的工作性能,研究壁面滑移、壁面热损失、工作环境等因素对其工作性能的影响规律,并初步分析了两相流因素的影响。