基于聚集式太阳能发电的流化床颗粒接收器降阶建模研究学术报告
一、 主要作者、机构与发表信息
本研究由来自美国科罗拉多矿业大学机械工程系的Keaton J. Brewster，美国国家可再生能源实验室的Janna Martinek，以及Brayton Energy LLC等多家机构的联合研究团队共同完成。通信作者为科罗拉多矿业大学的Gregory S. Jackson教授。该研究以论文“Reduced order modeling of a fluidized bed particle receiver for concentrating solar power with thermal energy storage”的形式，于2025年2月12日在线发表在学术期刊《Solar Energy》第289卷上。
二、 学术研究背景
本研究属于清洁能源技术领域，具体聚焦于下一代聚集式太阳能发电技术中的高温颗粒接收器与热储能系统的设计与优化。
研究背景与动因：传统的聚集式太阳能电站通常使用熔融硝酸盐作为传热流体和热储能介质，但其工作温度被限制在570°C以下，且存在腐蚀性挑战，制约了电站发电效率的进一步提升。为了提高效率，下一代超临界二氧化碳布雷顿循环的涡轮机进口温度需要超过700°C。因此，开发能够在750°C以上温度运行的高温传热与储能介质至关重要。氧化物颗粒因其成本低廉、耐高温（可达1200°C）、适合长期储热等优点，被视为理想的候选介质。然而，颗粒介质的低体积热导率使得高效地将聚集的太阳能热量传递至颗粒流成为一项关键技术挑战。
核心问题与目标：为了解决这一问题，研究团队提出了一种间接式流化床颗粒接收器的概念。与直接将颗粒暴露于太阳辐射的“直接式”接收器不同，间接式接收器使用金属壁面（如高温合金）来容纳颗粒流。这种设计可以允许使用成本更低、太阳能吸收率更低但尺寸更小的颗粒（如硅砂），并能防止颗粒逸散。然而，间接式接收器的核心挑战在于，必须确保颗粒流与容纳壁面之间具有极高的换热系数，才能在极高的太阳辐射通量下，将壁面温度控制在结构材料的极限（如低于950°C）以内，同时将颗粒加热至750°C以上的目标温度。流态化技术是实现高壁面-床层换热系数的关键途径。
研究目的：本研究旨在开发一个降阶模型，用于模拟和分析一种逆流、窄通道、鼓泡流化床颗粒接收器的热工水力性能。该模型集成了经过实验校准的壁面换热和颗粒轴向分散关联式，以探索在不同操作参数（如颗粒质量通量、气体质量通量、太阳辐射通量、颗粒分散特性）下，流化床如何实现高效的间接式颗粒接收，并为兆瓦级接收器的设计和放大提供理论基础。
三、 详细研究流程与方法
本研究采用基于物理定律和实验关联式的数值建模方法，构建了一个准二维降阶模型来模拟接收器中心通道的性能。整个研究流程可概括为以下步骤：
1. 模型系统定义与几何构建 * 研究对象：研究的核心是一个原型间接式腔体接收器的单个中心通道。该通道被设计为具有浅角度的腔体壁，以降低壁面上的有效法向太阳通量（将980 kW/m²的孔径通量降低至壁面处低于250 kW/m²）。通道内，硅砂颗粒（Sauter平均直径dp = 225 μm）在重力作用下向下流动，而流体化空气从底部向上逆流流动，形成鼓泡流化状态。 * 几何分区与离散化：模型将接收器通道在深度方向（x方向）划分为三个区域（Zone 1, 2, 3），以对应不同的太阳通量分布条件。每个区域在垂直高度方向（y方向）又进一步被离散为多个计算节点（通常使用51个节点），从而形成一个准二维的计算网格。每个节点求解气固两相及壁面的状态变量。
2. 控制方程与物理模型建立 * 控制方程体系：模型建立了稳态下的质量、动量和能量守恒方程组，分别针对气相、固相（颗粒）以及接收器壁面的内部和外部。这些方程以微分形式列出，涵盖了气体和颗粒的流动、压力分布、体积分数、温度分布以及壁面导热和辐射换热。 * 关键输运系数关联式（本研究核心创新与校准部分）： * 床层-壁面换热系数：采用了团队先前通过实验标定的半经验关联式。该关联式基于Molerus的理论，通过无量纲阿基米德数和过量流化速度的函数来预测对流换热努塞尔数，并叠加了辐射换热贡献。对于本研究使用的225微米硅砂颗粒，在温和鼓泡流化条件下，预测的ht,w可超过1000 W m⁻² K⁻¹。 * 轴向分散系数：引入了由佩克莱特数关联的轴向热分散模型。该模型将颗粒轴向热分散系数与过量流化速度线性关联，其佩克莱特数值（Pe=13.02）通过对实验室窄通道流化床热交换器的实验数据进行拟合获得。此分散效应模拟了上升气泡导致的颗粒上下混合，从而平滑了床层内的垂直温度梯度。 * 横向分散系数：由于缺乏直接的实验测量，模型中引入了一个可变的横向分散系数Dxx,s，并将其与轴向分散系数Dyy,s的比值作为关键参数进行研究。文献表明，横向分散通常小于轴向分散，其比例在1%到50%之间变化。本研究中系统性地探讨了该比值（从0.01到1.0）对接收器性能的影响。 * 边界条件与外部辐射： * 太阳辐射通量输入：外部壁面吸收的太阳辐射通量不是简单假设，而是通过专业的光学模拟软件SolTrace计算得到。该软件对一个极坐标定日镜场聚焦到接收器腔体孔径上的过程进行了蒙特卡洛射线追迹模拟，从而获得了腔体各壁面上真实、非均匀的局部吸收通量分布。此通量分布随后被空间平均并映射到模型的三个区域及各垂直节点上。 * 壁面辐射交换：模型计算了外部壁面之间以及壁面向环境的红外辐射热交换。通过蒙特卡洛方法计算了各区域壁面之间的视角因子，并将其用于辐射换热计算中。
3. 数值求解与模型验证 * 数值方法：将连续的控制方程在离散节点上采用有限体积法进行离散化。对流项采用迎风格式，扩散项采用中心差分。最终形成的大型非线性代数方程组，使用MATLAB中的‘fsolve’求解器（采用信赖域狗腿算法）进行求解。 * 模型验证：为确保模型的可靠性，研究团队将模型预测结果与其实验室规模的窄通道流化床实验数据进行了对比。在相同的颗粒质量通量、气体流速和外部加热条件下，模型成功预测了实验测量的颗粒垂直温度分布和床层压力梯度分布，验证了模型核心的传热和流动物理模型的准确性。
4. 参数化研究与性能分析 * 参数设置：在验证模型的基础上，对原型接收器进行了广泛的参数化研究。研究的变量包括：颗粒质量通量（10-50 kg m⁻² s⁻¹）、气体质量通量（0.01-0.25 kg m⁻² s⁻¹）、峰值吸收太阳通量（100-400 kW m⁻²）、接收器高度（1.1-4.6 m）以及横向与轴向分散系数之比（0.01-1.0）。 * 性能指标：主要评估指标包括：接收器热效率（定义为颗粒吸收的热量与壁面吸收的太阳能之比）、颗粒出口温度（目标>750°C）、以及最关键的外壁面最高温度（需低于材料极限，如950°C）。此外，还分析了不同区域之间温度和热通量的均匀性。
四、 主要研究结果
1. 基准工况性能分析：在基准工况下（孔径峰值通量980 kW/m²，对应壁面吸收峰值188 kW/m²，颗粒通量25 kg m⁻² s⁻¹，气体通量0.1 kg m⁻² s⁻¹，Dxx,s/Dyy,s = 0.4），模型预测表明： * 温度分布不均匀性：接收器前沿区域（Zone 1）由于承受最高的太阳通量，其颗粒温度和壁面温度显著高于后方区域（Zone 2, 3）。Zone 1的最大外壁面温度接近970°C，而Zone 2和3则低得多。这揭示了间接式接收器面临的核心挑战：局部热点。 * 流化床的作用：即使气体质量通量仅为颗粒的0.5%，鼓泡流化产生的高床层-壁面换热系数（>1000 W m⁻² K⁻¹）使得在如此高的壁面热通量下，平均颗粒出口温度能够达到约740°C，满足了高温储能的要求。轴向分散效应有效平滑了垂直方向的温度梯度。
2. 气体质量通量的影响：增加流体化气体质量通量对性能有双重积极影响。 * 提升换热系数：在一定范围内（约 ?̂? =40-60），增加气体通量可提高ht,w，从而降低在相同热通量下的壁面与颗粒的温差。 * 增强混合，降低不均匀性：气体通量增加同时增强了轴向和横向（由于模型假设Dxx,s与Dyy,s相关）的颗粒热分散。结果表明，当气体通量从接近最小流化状态增加到0.1 kg m⁻² s⁻¹时，三个区域间的颗粒出口温度差异从超过80°C缩小到30°C以内，且最大壁面温度可被抑制在950°C以下。热效率也随之提高。
3. 横向分散的关键作用：参数研究表明，横向颗粒热分散是决定接收器壁面温度均匀性和峰值的关键因素。 * 当横向分散系数与轴向分散系数之比从0.01增加到1.0时，前沿区域（Zone 1）的最高壁面温度可降低约60°C，同时三个区域的颗粒出口温度变得更加均匀。 * 这一发现具有重要的工程意义：它表明，通过改进接收器设计（如采用多向气体喷射等方式）来主动增强横向混合，可以有效缓解局部过热问题，降低对前沿区域壁面材料的苛刻要求，并减少因温差过大引起的热应力。
4. 太阳辐射通量与接收器尺寸的权衡：研究探索了在维持颗粒出口温度~750°C的前提下，改变太阳通量浓度和接收器高度的权衡关系。 * 提高峰值太阳通量可以减少所需接收器尺寸，从而提高热效率（因为相对热损失面积减小）。然而，这会导致前沿区域的壁面温度急剧升高。例如，将吸收峰值通量提高到250 kW/m²以上，会使Zone 1壁温超过1000°C，超出了当前高温合金的承受范围。 * 相反，降低太阳通量（增大接收器尺寸）虽然能降低壁温，但会导致热效率下降（低于65%），并增加成本。 * 核心结论：对于当前基于平滑壁面和高强度合金的设计，存在一个可行的操作窗口。基准工况（孔径通量~1000 suns，壁面吸收峰值~190 kW/m²）接近该窗口的上限，壁温已接近材料极限。
5. 性能提升途径探索：模型分析指出了未来改进的方向。 * 强化换热：模拟显示，若能将有效床层-壁面换热系数从1000提升至1500 W m⁻² K⁻¹，则可在保持相同壁温（如950°C）的前提下，允许壁面热通量增加约50%；或在相同热通量下，显著降低壁温。这为在壁面增加扩展表面等强化换热技术提供了理论依据。 * 降低热损失：当前设计的接收器热效率（70-78%）低于85%的常规目标。主要热损失来自前沿区域向环境的大幅辐射。优化腔体几何形状以增强“光陷阱”效应，或采用腔中腔设计，是未来提高效率的重要途径。
五、 研究结论与价值
结论：本研究成功开发并应用了一个经过实验校准的降阶模型，系统评估了逆流流化床间接式腔体颗粒接收器的可行性。模型证实，采用硅砂颗粒和温和鼓泡流化技术，可以实现足够高的床层-壁面换热系数，从而在接近1000 suns的高孔径太阳通量下，将颗粒加热至750°C以上，同时将合金壁面最高温度控制在950°C左右。研究明确了操作参数（尤其是气体流率和横向混合强度）对接收器性能（效率、温度均匀性、峰值壁温）的显著影响。
科学价值： 1. 提供了关键的设计工具：所建立的降阶模型计算高效，能够进行大范围的参数化研究，弥补了全尺度计算流体力学模拟耗时过长的不足，为接收器的初步设计和优化提供了强有力的分析工具。 2. 揭示了核心物理机制：定量分析了流化床中高换热系数与颗粒热分散（尤其是横向分散）对缓解间接式接收器局部热点、提高运行均匀性的核心作用，深化了对该复杂多相流传热系统的理解。 3. 明确了性能边界与权衡关系：清晰勾勒了在给定材料和设计下，太阳通量浓度、接收器尺寸、热效率与壁面温度之间的相互制约关系，为工程决策提供了科学依据。
应用价值： 1. 指导兆瓦级示范装置设计：研究结果直接为后续建造更大规模的流化床颗粒接收器原型机指明了关键设计准则和操作策略，例如需要重点优化气体分布以增强横向混合，以及探索强化换热表面。 2. 推动下一代CSP技术发展：该研究工作为将低成本颗粒热储能与高效超临界二氧化碳发电循环相结合的技术路线扫清了一个关键部件——高温颗粒接收器——的理论障碍，有助于加速这一具有前景的、可调度可再生能源技术的商业化进程。
六、 研究亮点
七、 其他有价值内容
本研究的数据与软件将根据请求向对应作者提供，体现了研究的可重复性与开放性。此外，论文中详细列出了所使用硅砂颗粒和壁面材料的关键物性参数，为其他研究者进行类似模拟提供了便利。研究也得到了美国国家科学基金会和美国能源部太阳能技术办公室的资助，显示了该课题在清洁能源领域的重要战略地位。