乐于分享
好东西不私藏

航天器等离子体相互作用软件参数设置教程(ii)

航天器等离子体相互作用软件参数设置教程(ii)

点击上方蓝字

关注

Opac Channel

SPIS(航天器等离子体相互作用软件)是一款功能强大且先进的三维建模工具,用于模拟航天器与等离子体的相互作用,能够模拟与航天器充电相关的多种现象。SPIS 包含一个基于非结构化网格的先进三维静电等离子体求解器,能够动态模拟复杂且真实的几何形状周围的精细鞘层结构。SPIS 集成了多种模型(混合模型/全 PIC 模型/回溯模型等),能够处理大规模的等离子体条件和物理过程。

多年来,SPIS 已成为欧洲航天器充电分析的事实标准参考工具。

如今它已演变成几种不同的版本:

  • SPIS-SC 用于表面电荷分析。它包含了多年来不断改进的多项功能:

  • GEO/LEO,增加 GEO/LEO 任务的表面充电分析和快速工程分析;

  • SCI,为表面等离子体相互作用领域的高级研究和科学应用提供支持;

  • DUST,增加了尘埃等离子体及其与着陆器相互作用的建模,并考虑了月球或小行星的地面地形。

  • EP,增加了对电推进器对空间系统影响的建模。

  • SPIS-IC 专注于内部充电分析

  • SPIS-ESD 专门用于对太阳能电池阵列上的静电放电进行建模

本推送将着重讲述在SPIS 6.2版本中进行Global parameters的设定,

自定义global parameters参数各选项的核心思路:

文件名
内涵
适合做什么
对于“空间天气—航天器电势”课题的价值
globalParameters-SPISGEO-ECSSGeoWorstCase-EclipseExit.xml
GEO 轨道、ECSS 标准体系下的地球同步轨道最坏充电工况,且包含 Eclipse Exit 食出状态
模拟航天器刚出地影、光照恢复、光电子发射突然变化时的充电过程
很重要。适合做“食出后推进器点火,羽流/CEX 离子如何重新分布,敏感元器件表面电位是否突变”
globalParameters-SPISGEO-ECSSGeoWorstCase.xml
GEO 轨道、ECSS 标准下的 GEO 最坏充电工况,不强调食出瞬态
作为欧洲标准体系下的 GEO charging worst case 基准
最适合作为第一套空间天气高风险基准工况
globalParameters-SPISGEO-NASAWorstCase.xml
GEO 轨道、NASA 标准体系下的最坏充电工况
与 ECSS worst case 交叉比较
很适合做稳健性验证:如果 ECSS 和 NASA worst case 下,某个元器件都高风险,结论更有说服力
globalParameters_standardVacuumChamber.xml
标准真空舱环境参数
地面实验、真空舱测试、羽流实验室复现实验
不适合直接代表空间天气,但适合模拟地面 plume test,或对比“地面真空舱结果”和“在轨环境结果”
SPIS-DUST_MS_parameters.xml
SPIS-DUST 模块参数。从命名看,MS 多半对应 magnetospheric / magnetosphere scenario,具体含义要打开 XML 确认
带电尘埃、尘埃粒子在磁层等离子体中的充电和运动
研究溅射产物、尘埃污染、二次颗粒,否则先不用
SPIS-DUST_PS_parameters.xml
SPIS-DUST 模块参数。PS 的确切含义需要看 XML 内容,可能对应 planetary surface / plasma sheath 一类 dust 场景
尘埃环境、表面附近带电颗粒问题
目前不建议用于电推进羽流侵蚀主线
SPIS-DUST_SW_parameters.xml
SPIS-DUST 模块参数。SW 通常可理解为 solar wind 场景
太阳风中的尘埃充电与运动
如果研究深空/太阳风环境下的尘埃,可用;若研究 GEO 附近 Hall 推进器羽流,暂时不用
globalParameters-5.1.0_modified.xml
SPIS 5.1.0 版本某个默认 global parameter 的修改版
个人或软件包自带的改动模板
它可能改过时间步、粒子数、输出监控、环境参数。没有对比 XML 前,不建议作为论文级仿真基准

以GlobalParameters-5.1.0为例,我们介绍几项参数的含义和常用值:

1. densitiesMapsMonitorStep = -5 或 -10

位置:

Global parameters → Outputs

作用:控制 Volume Densities 里密度场的输出次数。

按 SPIS 这种写法:

0     = 不输出密度图
-5    = 整个仿真过程中输出 5 次密度图
-10   = 整个仿真过程中输出 10 次密度图

要看羽流形态,就必须保留这个参数。
如果设成 0,跑完以后可能没有这些:

log10 of improved Thruster0_3_Xe+ charge density at t = ...
log10 of improved Thruster0_4_Xe+ charge density at t = ...

建议:

densitiesMapsMonitorStep = -10

如果想减少输出、加快一点,可以用:

densitiesMapsMonitorStep = -5

2. densityLogPlotFlag = 2

位置:

Global parameters → Outputs

作用:控制是否输出 log10 密度图

建议保留:

densityLogPlotFlag = 2

3. energyMapMonitorStep = 0

位置:

Global parameters → Outputs

作用:控制是否输出体域能量图,比如粒子能量分布之类。

设:

energyMapMonitorStep = 0

这不会关闭粒子运动,也不会关闭物理计算,只是不输出能量图,减少文件和后处理负担。


4. numericsMapsMonitorStep = 0

位置:

Global parameters → Outputs

作用:控制数值诊断类图,比如宏粒子数、数值状态图等。

它对应 Data Mining 里看到的:

Volumic Number Of Macro Particles

未完待续......

 Produced by Opac Channel 

欢迎扫描上方二维码,或搜索

微信号|Oppaic

加入Opac Channel!

点击赞和在看,您的分享对我们的帮助真的很大!