一种中性气体释放诱发航天器放电的数值模拟方法，在初始或特定条件下给出场致电流、热电流等种子电流，并考虑周围等离子体以及材料放气的初始状态，然后分别通过粒子模拟方法和蒙特卡洛法求解空间中的粒子轨迹和粒子与气体碰撞反应过程，研究这一过程中气压和放电激烈度的变化，从而得到不同初始条件下，中性气体释放诱发航天器放电的影响程度和规律。本发明能够辅助分析中性气体释放在航天器表面充放电过程中的影响机制，通过计算分析不同气压不同结构的航天器带电部件在充放电过程中可能出现的非自持到自持，再到雪崩放电的转换阈值，并且可以作为基础方法研究多类复杂环境下共同诱发航天器充放电的过程，以便科学有效的进行航天器设计及加固。