面向空间应用的光伏太阳能电池分析工具开发：基于双二极管模型与光学建模的Python仿真平台

《Materials Today Electronics》：Development of an analytical tool to design photovoltaic solar cells: Analysis in outer space conditions

本文推荐了一款基于Python开发的太阳能电池仿真工具，重点针对空间应用环境（如AM0光谱、电离辐射效应）进行优化。该工具采用面向对象编程，集成双二极管等效电路模型（2M5P）、半导体物理和光学传输矩阵法（TMM），可模拟单/多结电池、抗反射涂层（ARC/DLARC）及保护层（PL）的性能。通过对比商业软件（如COMSOL、PC1D）和实际空间电池（如AZUR Space Si电池）数据，验证了其准确性（效率误差<0.6%）。研究还展示了抗反射涂层设计、电离辐射损伤模拟及多结电池电流匹配优化案例，为空间太阳能电池设计提供了开源、高精度的仿真解决方案。

　　

引言

太阳能电池自1958年首次应用于航天任务（美国Vanguard I卫星）以来，已成为空间能源系统的核心部件。空间环境与地表存在显著差异：太空光谱为未经过大气层吸收的AM0标准（辐照度1366 W/m²），较地表AM1.5光谱强度更高；极端温度波动、等离子体撞击、原子氧侵蚀及高能粒子辐射（如范艾伦带）对电池可靠性提出严苛要求。现有商用仿真软件（如COMSOL、PC1D）多属闭源付费工具，且缺乏针对空间环境的专项优化。本研究开发了一套基于Python的开源分析工具，通过双二极管电学模型、半导体物理参数和光学建模方法，实现空间太阳能电池的精准仿真。

方法论

软件架构与光学模型

程序采用对象导向设计，依赖NumPy数值计算、Matplotlib数据可视化及JSON配置文件管理。光学模型基于垂直入射假设，针对不同层结构采用差异化处理：

• •
  抗反射涂层：应用传输矩阵法计算多层薄膜系统（如TiO₂/Al₂O₃）的干涉效应，其传输矩阵由界面矩阵与传播矩阵乘积构成，可量化反射损失。
• •
  保护层：对玻璃、EVA等厚层材料，假设消光系数趋近零，采用修正TMM处理光强传输。
• •
  结区光传播：通过菲涅尔方程和比尔-朗伯定律循环计算光子通量在结内的透射、反射与吸收过程，累加生成光生电流。

电学模型与参数计算

• •
  光生电流：按发射区、基区和耗尽区分别计算。发射区空穴扩散电流密度公式包含吸收系数α、空穴扩散长度L_p及表面复合速度S_p参数；基区电子扩散电流则依赖电子扩散系数D_n与厚度W_p；耗尽区因强电场作用具有近100%量子效率。
• •
  暗电流：采用双二极管模型，第一二极管模拟扩散电流（理想因子n=1），饱和电流I_diff,0由基区和发射区组分构成；第二二极管表征耗尽区SRH复合电流（n=2），其饱和电流I_dr,0与载流子寿命τ_n、τ_p相关。
• •
  串联电阻：包含基区电阻（与基区厚度W_p和电阻率ρ_b成正比）、发射区电阻（与指栅间距S_f和方块电阻r_□相关）、接触电阻（通过热电子发射/场发射模型计算）及指栅电阻（由银电极几何尺寸决定）。
• •
  电离辐射效应：通过引入陷阱密度N_t和俘获截面σ修正载流子寿命，公式1/τ = 1/τ₀ + σN_tv_th，其中v_th为热速度。
• •
  输出特性求解：采用牛顿-拉夫森法迭代求解包含串联电阻R_s和分流电阻R_sh的隐式电流方程，多结电池需额外进行电流匹配计算。

结果与验证

半导体数据来源

关键材料（Si、Ge、GaAs、Ga_0.51In_0.49P）参数取自《Handbook series on semiconductor parameters》，折射率数据库及Sze等人的电阻率-掺杂浓度曲线（图4）。GaInP的载流子迁移率、寿命等参数通过GaP和InP的线性插值估算。

仿真准确性验证

• •
  与COMSOL对比：Si电池模型在AM1.5光谱下，本工具效率达16.38%，与COMSOL（16.35%）偏差<0.2%；外量子效率曲线在短波区吻合度较高。
• •
  与PC1D对比：GaAs电池（含TiO₂/Al₂O₃ DLARC）效率为22.21%，略低于PC1D（25.37%），差异源于串联电阻模型及光学层简化处理。
• •
  商业电池对标：AZUR Space Si空间电池（S 32型）仿真效率16.33%，与 datasheet标称值（16.9%）误差<3.4%，验证了DLARC模型及接触电阻计算的实用性。

抗反射涂层与保护层研究

以GaAs电池为例，对比不同光学设计：

• •
  裸电池：平均反射率32.46%，效率仅13.78%。
• •
  单层ARC：理想折射率n_ARC = √(n₀n_s)≈1.87（λ=600 nm时）。Al₂O₃层（n=1.65）因低消光系数优于ZnO（n=1.78但k值高），效率提升至19.44%。
• •
  双层ARC：最优结构满足n_top = (n₀²n_s)^1/3≈1.57，n_bot = (n₀n_s²)^1/3≈2.46。SiO₂/TiO₂组合在500–800 nm波段透射率最优，效率达20.36%。
• •
  保护层附加：EVA/玻璃层导致效率降低约1%（至19.32%），体现防护性能与效率的权衡。

电离辐射损伤模拟

基于1 MeV电子辐照实验数据（表12），分析Si电池性能退化：

• •
  单个陷阱效应：空穴陷阱主要降低开路电压（如H1陷阱使V_oc从0.58 V降至0.42 V），电子陷阱同时削弱短路电流（E2陷阱使I_sc减少21%）。
• •
  累积效应：当总注量达10¹⁶ e/cm²时，载流子寿命急剧下降，SRH复合电流激增，导致填充因子（FF）从0.774恶化至0.337，效率仅存0.78%。但模型退化速率较实际电池（文献报道效率保留5.3%）偏快，原因可能未考虑陷阱能级位置对复合效率的影响。

多结电池优化

研究GaInP/GaAs/Ge三结电池的电流匹配问题：

• •
  初始设计：Ge底结电流（35.73 mA）远高于GaAs中结（14.66 mA），导致效率仅22.96%。
• •
  结构优化：将GaInP顶结厚度从0.75 μm减至0.31 μm，掺杂调整为N_d=7×10¹⁸ cm^-3，使GaAs结电流提升至16.97 mA，效率改善至27.73%。外量子效率曲线显示GaAs结在短波区响应增强，而GaInP结吸收减弱，实现更好的光谱分配。

结论

本研究成功开发了一款面向空间应用的太阳能电池仿真平台，通过集成电学、光学与辐射损伤模型，实现了对单结/多结电池的精准分析。工具在抗反射涂层设计、电离辐射效应模拟及多结电池电流匹配等方面展现出实用价值，为航天器能源系统设计提供了开源、低成本的计算基础。未来工作可聚焦于陷阱能级敏感性建模、纹理表面光学优化及实际空间环境验证，进一步提升仿真置信度。