这份题为《面向太阳能塔式电站系统级分析的落粒式粒子接收器热模型》的研究论文，发表于Elsevier旗下期刊《Solar Energy》的第268卷（2024年）。论文第一作者及通讯作者单位均为意大利米兰理工大学能源系。研究旨在开发一种计算高效、适用于系统级优化的落粒式粒子接收器热模型，并对不同接收器配置进行性能分析与优化。
学术背景与研究目标
太阳能热发电（Concentrating Solar Power, CSP）是应对气候变化、实现能源可持续发展的重要技术路径。在“第三代CSP技术路线图”中，美国能源部将基于粒子的接收器技术列为实现更高运行温度和更低成本目标的最有前景方案。落粒式粒子接收器作为其中一种关键技术，具备直接加热、可承受超过1000°C高温、且储热介质成本低廉等显著优势。
然而，要充分发挥该技术的潜力并优化整个电站系统的设计与运行，需要一个能够准确描述接收器在不同工况下行为、且计算速度足够快的热模型。尽管文献中存在一些复杂的计算流体动力学（CFD）模型或简化模型，但前者计算成本过高，不适用于系统级的重复优化计算；而后者往往忽略了许多关键物理效应，精度有限。因此，本研究的目标是：在已有简化模型的基础上，集成一系列关键物理特征，开发一个更为精细、准确且计算高效的落粒式粒子接收器热模型，并应用该模型进行接收器尺寸优化、非设计工况及全年性能分析，比较不同配置（如自由落体与多级落体）的性能差异。
详细研究流程与方法
本研究的工作流程是一个从模型开发、验证到系统应用的综合过程，主要包括以下几个步骤：
1. 热模型开发与增强 本研究以González-Portillo等人（2021年）提出的“基础模型”为起点，该模型是一个沿粒子下落方向的一维稳态模型，并已通过与CFD结果的对比得到初步验证。为了提升模型的精度和适用性，本研究在基础模型上集成了多项关键改进： * 二维/三维离散化：将粒子帘的离散化从一维扩展到二维（宽度和下落方向），从而可以处理非均匀的聚光热流密度图。将后墙的离散化从一维扩展到三维（宽度、高度和厚度方向），能够模拟后墙内的三维导热，并允许设置多层墙体结构（如内部耐高温层、中间隔热层和外部保护层）。 * 阻力效应：在粒子帘的动量平衡方程中引入了基于Hruby等人工作的阻力模型，以考虑颗粒与腔内空气之间的拖曳力对颗粒速度的影响。这是对基础模型的重要补充。 * 可变光学与热物性参数： * 采用González-Portillo等人开发的光学模型，根据颗粒体积分数、帘厚度等参数动态计算粒子帘的吸光率、透光率和反射率。 * 粒子的比热容设定为温度的函数，采用由美国桑迪亚国家实验室拟合的经验公式。 * 后墙反射率根据辐射来源（太阳光谱或热辐射光谱）采用不同的值，以更真实地模拟实际接收器中为减少吸热和热损失而设计的墙面特性。 * 改进的对流换热与视因子计算：采用更符合物理的对流换热系数关联式来评估后墙向环境的自然对流与强制对流损失。此外，引入了三维计算的等效视因子，以更精确地描述粒子帘与接收器空腔开口之间的辐射换热交互作用。
模型的物理域包括粒子帘和接收器后墙。控制方程主要包括粒子帘的质量、动量与能量守恒方程，以及后墙的三维导热方程。这些方程在离散化的网格上联立求解，最终得到粒子温度、速度、体积分数以及后墙温度场等关键参数的分布。
2. 模型验证 为了确保所开发模型的可靠性，研究团队进行了两个层次的验证： * 阻力模型验证：将集成的阻力模型计算结果与Hruby等人文献中的模型结果及实验测量数据进行比较。对比显示，在高温条件下（773K），模型预测的颗粒速度与实验数据吻合良好。尽管在某些低温、低质量流率工况下存在高估，但考虑到商业应用的目标质量流率（10-20 kg/s/m）远高于验证所用值，且高温工况更具代表性，该模型被认为是可接受的。 * 整体热模型验证：将开发的热模型与Mills等人文献中的CFD模拟结果进行对比。对比了三种模型变体：不含任何改进的“基础模型”、包含所有改进但不含阻力效应的“无阻力”模型、以及包含所有改进的“完整”模型。结果显示，“完整”模型与CFD结果符合得最好，其总热效率的绝对误差为0.8%，优于“基础模型”的1.0%。这表明所引入的各项改进有效提升了模型的预测精度，同时保持了较低的计算成本（单工况模拟约4分钟），验证了其适用于系统级分析的可行性。
3. 案例研究与系统级分析 在验证之后，研究团队将模型应用于一个具体的100 MWe超临界二氧化碳（sCO2）发电塔式电站案例中，地点选在美国加利福尼亚州达格特。分析流程如下： * 光学分析：使用SolarPilot软件生成特定太阳能场布局下的聚光热流密度图，并计算不同开口尺寸下的光学效率。 * 设计点优化：针对不同的接收器配置（自由落体、2级、3级、4级、5级落体以及带再循环的自由落体），使用热模型计算其热效率。通过变化接收器开口面积，寻找光学效率与热效率的乘积——“光-热转换效率”的最大值，从而确定每种配置的最优尺寸。瞄准策略以保证后墙内层最大热流密度不超过2 MW/m²为约束。 * 非设计工况性能评估：针对优化后的配置，通过改变直接法向辐照度（DNI）来模拟不同太阳辐照条件下的接收器性能，绘制热效率随入射热功率变化的曲线，并评估最低运行功率限制。 * 全年性能分析：将接收器的非设计性能曲线与太阳能场的小时光学性能、天气数据、储热系统（TES）和sCO2动力循环（PB）模型集成，进行全年的逐时动态模拟。运行策略为优先保证PB在额定热功率下运行，不足时由TES补充，过剩时存入TES。同时，评估了再循环策略在部分负荷下切换使用的经济性（基于净热效率比较）。
主要研究结果
1. 设计点优化结果 * 随着级数增加，最优接收器开口尺寸增大（从自由落体的28m×28m增至五级落体的29.5m×29.5m），光-热转换效率也相应提高。 * 自由落体接收器的设计点热效率为81.1%，光-热转换效率为49.5%。五级落体接收器的设计点热效率最高，达到83.1%，光-热转换效率为51.5%。 * 研究发现，在颗粒离开收集槽时假设其为理想混合（温度均匀）与无混合两种极端情况下，热效率差异极小（83.1% vs. 83.0%），因此后续分析采用理想混合假设是合理的。
2. 模型内部物理场分析 对最优尺寸下的自由落体和五级落体接收器内部物理场分析显示： * 温度分布：自由落体接收器由于聚光热流不均匀，粒子帘宽度方向温度差异超过200°C。五级落体接收器因每一级后的混合作用，温度分布更为均匀。 * 速度与体积分数：颗粒在下落过程中因温度升高导致空气粘度增加、密度减小，拖曳力增大，会出现减速现象，而非达到恒定终端速度。颗粒体积分数在起始处较高，但随着帘厚增加迅速下降。多级配置通过在每一级起始处“重置”帘厚，使体积分数周期性出现峰值，提高了帘的平均不透明度。 * 光学性质：粒子帘的反射率主要取决于体积分数，体积分数高时反射率也高。
3. 非设计工况与全年性能结果 * 非设计性能曲线：多级配置在所有入射功率水平下都展现出比自由落体更高的热效率，且在低负荷时优势更明显，因为帘的不透明度更高。五级落体配置具有最宽的操作范围（最低入射功率要求153 MW）。 * 再循环策略效果有限：虽然再循环（将粒子帘在宽度方向分为两半并循环下落）能略微扩展自由落体接收器的操作范围并提升低负荷效率，但其带来的全年净收益非常有限。 * 全年性能汇总： * 自由落体接收器的全年平均热效率为72.0%，年净发电量567 GWh，容量因子64.7%，系统整体（太阳能到电能）年效率为17.6%。 * 五级落体接收器性能最优，全年平均热效率为75.3%，年净发电量597 GWh，容量因子68.2%，系统整体年效率为18.5%。 * 再循环配置相比自由落体，年净发电量仅增加3 GWh（570 vs. 567 GWh），考虑到额外的设备成本和控制复杂性，其应用价值可能不高。
研究结论与价值
本研究成功开发并验证了一个适用于太阳能塔式电站系统级分析的落粒式粒子接收器热模型。该模型在计算效率与物理细节之间取得了良好平衡，能够模拟二维非均匀热流、三维墙体导热、阻力效应、可变物性等关键现象。
通过案例研究，论文得出以下重要结论： 1. 多级落体接收器性能更优：与自由落体接收器相比，采用多级落体结构能有效提升接收器在设计点和非设计点的热效率，且级数越多，性能越好。五级落体配置在全年性能上表现最佳。 2. 明确了性能提升潜力与局限：即使采用最优的五级落体设计，其年热效率（75.3%）仍未达到美国“SunShot”计划设定的90%目标，表明该技术仍有进一步研发和改进的空间。论文提出的模型为此类优化研究提供了有力工具。 3. 评估了运行策略：研究指出，再循环作为一种部分负荷运行策略，对全年性能的提升效果有限，可能不足以证明其增加的资本和运营成本的合理性。
研究亮点
其他有价值的内容
论文在附录中进行了详尽的敏感性分析，探讨了颗粒与空气速度比、对流换热系数等关键参数的不确定性对模型结果的影响，增强了研究结论的稳健性。此外，附录还讨论了阻力模型修正因子、后墙加权反射率计算方法等细节，体现了研究的严谨性。论文最后还展望了未来可能的模型改进方向，如开发基于局部速度的对流换热关联式、考虑风的影响、以及对多级接收器中收集槽进行更详细的建模等。