CFD建模研究转速对大型IV型聚合物复合氢气储罐热行为及温度偏移最小化的影响

《Polymers》：CFD Modeling of Rotational Speed Effects on Thermal Behavior and Temperature Excursion Minimization in Large Type IV Polymer Composite Hydrogen Storage Tanks

　　

在快速加注过程中，大型IV型聚合物复合氢气储罐会经历显著的温度梯度，这与气体压缩及焦耳-汤姆逊效应（Joule-Thomson effect）相关，且可能危及容器完整性，严重影响整体安全性。为解决该问题，本研究提出一种新型主动混合方法，即储罐旋转，以增强内部对流传热并最小化温度梯度。研究人员采用Redlich-Kwong真实气体状态方程进行瞬态计算流体力学（CFD）模拟，精确捕捉氢气在高压下的热力学行为。该研究基于对20.56 m³内体积储罐进行1000 s快速加注的工况，分别评估了10、30和50 rad/s的切向旋转速度。结果表明，热性能与转速之间存在强烈的非线性依赖关系。在10 rad/s时，形成了较为均匀的温度分布，且能耗较低。30 rad/s时实现了峰值温度的最显著抑制，即最高效的冷却效果。然而，当转速进一步升高至50 rad/s时，大量的粘性耗散热产生了不利的二次温升，部分抵消了强化混合带来的益处。总体而言，结果表明与30 rad/s更为相关的理想运行窗口能够在大型复合氢气储系统快速加注过程中实现温度均匀性、最高温度限制和能耗之间的最有利折中。

## 一、研究背景、问题与研究意义

氢能作为清洁能源体系的核心组成部分，其储运安全性与效率直接制约着质子交换膜燃料电池（Proton Exchange Membrane Fuel Cells, PEMFCs）等终端应用的规模化发展。在各类高压储氢方案中，IV型储氢罐（type IV hydrogen storage tank）——以内衬高密度聚乙烯（high-density polyethylene, HDPE）、全包裹碳纤维复合材料的结构形式——因其优异的比强度及高压承载能力，已成为大容量氢气储存的最先进技术方案。然而，快速加注过程中氢气经历的近绝热压缩以及显著的焦耳-汤姆逊效应，导致罐内温度急剧攀升，形成极端热-机械应力，严重威胁复合材料层的结构完整性与长期耐久性。

国际技术标准SAE J2601明确规定，复合储罐内部最高温度不得超过85 °C（358 K），以确保材料刚度与系统安全。当前文献中占主导的热管理策略多为被动式方法，包括进站氢气预冷、优化入口扩散器设计、以及罐内置放散热体等。尽管这些方法在一定程度上改善了峰值温度表现，但其效果始终受制于对流换热能力的固有限制，且常伴随额外重量、系统复杂度增加及成本上升等问题，对于质量、体积和成本均受限的大容量氢储存应用尤为不利。在此背景下，无需显著硬件改动、不引入外部 stationary 能耗的主动式罐内对流调控技术成为研究前沿。然而，针对旋转混合这一策略在大规模IV型储氢系统中的系统研究尚属空白。

正是基于上述研究缺口，研究人员开展了这项创新性工作。该研究首次将转动参考系（Moving Reference Frame, MRF）方法系统应用于大容量IV型高压氢气储罐的主动热管理，旨在揭示转速对瞬态温度场的影响规律，识别对流强化与粘性耗散之间的最优平衡窗口，为新一代大容量储氢系统的加注规程设计与主动热管理概念开发提供理论依据。研究成果发表于《Polymers》期刊，对推动高压氢气储运技术的工程实践具有重要参考价值。

## 二、关键技术方法

研究人员基于CFD软件平台开展数值模拟，主要技术方法包括：采用瞬态可压缩湍流模型求解质量、动量和能量守恒方程；运用Redlich-Kwong真实气体状态方程（RK-EOS）精确刻画70 MPa高压范围内氢气的热力学行为，其密度与热力学性质预测精度与NIST REFPROP数据偏差在1–2%以内；选用可实现的k-ε湍流模型（realizable k-ε model）配合增强壁面处理以捕捉旋转诱导的强压力梯度与旋流流动；采用MRF方法在静止罐壁条件下对整个流体域施加恒定角速度（10、30、50 rad/s），引入科里奥利力（Coriolis force）与离心力的源项模拟整体旋转效应；针对20.56 m³大型储罐构建结构化四边形主导网格（约10万单元），并进行50,000至200,000单元的网格独立性验证；边界条件设定为：入口质量流率对应1000 s快速加注工况、入口温度253 K（预冷站端条件）、壁面为绝热无滑移边界、初始状态为2 MPa/298 K的部分耗尽储罐状态；数值格式采用二阶迎风格式与二阶隐式时间推进，收敛残差设定为10^?4（连续性与动量）和10^?6（能量）。模型经Dicken和Mérida的基准实验数据验证，峰值温度偏差小于2.4%。

## 三、研究结果

### 3.1 相对切向速度与流场结构

在10 rad/s转速下，储罐内形成高度相干的旋涡流动模式，呈现强环状涡旋结构，最大相对切向速度约6.36 m/s，近壁区域显著降低。流场呈现对称的双涡旋模式，中部存在滞止区，无速度间断与回流，流动组织良好，有利于抑制剪切生热并保证加注过程中温度均匀性。体积平均相对切向速度约2.5–3.0 m/s，足以提供对流输运以破坏热分层，同时避免过度粘性耗散。

### 3.2 湍流动能分布

30 rad/s时湍流动能（Turbulence Kinetic Energy, TKE）呈现强轴对称分布，峰值达38.3 m²/s²，出现在双环形高剪切区，中心区域相对平静。体积平均TKE较非旋转工况增加18–22 m²/s²，表明旋转对湍流生成的贡献显著增强。该转速下湍流输运足以破坏热分层并加速热移除，同时峰值TKE尚未引发主导的粘性耗散与反作用生热。

### 3.3 流函数比较

10 rad/s时流函数呈对称、有组织的环状涡旋结构，值约631 kg/s，径向混合充分；50 rad/s时流函数值升至1040 kg/s，上部形成高浓度流线聚集区，下部相对平静（200–300 kg/s）。体积平均流函数从10 rad/s的约250–300 kg/s增至50 rad/s的450–550 kg/s，表明旋转速度对整体循环强度的显著影响。中间转速在强化混合与控制剪切生热之间取得更好平衡。

### 3.4 静态温度分布

1000 s加注结束后，10 rad/s时中心存在明显高温区（峰值约271 K），温度梯度向壁面平滑过渡，但未完全消除热分层；30 rad/s时温度 profile 显著改善，峰值降至约256 K，中心热点大幅减弱，温度过渡更为平滑，实现三者中最均匀的温度分布与最低最高温度；50 rad/s时温度分布恶化，峰值回升至253 K以上，上部出现显著热点，对称性破坏，表明粘性耗散开始抵消混合收益。

### 3.5 压力分布

10 rad/s时压力场相对均匀（0.08–0.09 MPa），旋转对压力均匀化影响有限；30 rad/s时峰值压力 arousal 至约2.88 MPa，分布改善且更为对称，高压中心区域扩展；50 rad/s时局部压力峰值达1.144 MPa（上部集中），对称性受损，动量分布不均。平均压力场呈非单调发展，30 rad/s为压力均衡化的最优工况。

### 3.6 速度幅值分布

10 rad/s时速度温和有序，峰值约54 m/s，呈中心环状结构；30 rad/s时速度幅值显著增长至约186 m/s，高速区扩大且近壁速度提升至50–80 m/s，对流混合增强；50 rad/s时速度急剧变化，峰值接近297 m/s，近壁形成高强度剪切层，局部剪切率超过200 m/s，粘性耗散成为显著因素。

### 3.7 转子焓分布

10 rad/s时转子焓（rothalpy）呈对称穹顶状，中心高（约3.33×10⁵ J/kg）、壁面低；30 rad/s时径向更均匀，峰值约4.62×10⁵ J/kg，核心区扩大且过渡平缓，能量周向传输更有效；50 rad/s时峰值回升至约4.59×10⁵ J/kg且局部化于上部，对称性丧失，局部能量累积抵消低温优势。

## 四、讨论与结论总结

研究人员在讨论中深入分析了旋转热管理的物理机制与工程可行性。数值模拟明确揭示热性能与转速之间存在非单调关系：适度旋转通过增强湍流动能与对流传输有效降低峰值温度（30 rad/s时较非旋转工况降低约10–15 K），但过度旋转导致近壁高剪切区的粘性耗散热成为主导因素，引发局部二次温升。这一竞争机制——湍流驱动热移除与剪切诱导生热之间的权衡——与旋转换热器、搅拌化学反应器等领域的经典物理一致。

研究进一步探讨了实际工程应用路径，强调MFR模拟作为数值代理手段，并不意味着储罐整体必须物理高速旋转。可行的替代方案包括：（i）由入流动能或小型辅助电机驱动的内部旋转扩散器或推进器；（ii）固定安装中形成大规模涡旋结构的涡旋生成入口设计；（iii）外部低频振荡或摇动机理。对于车辆应用，整体机械旋转面临轴承设计、复合层疲劳管理、动态平衡、系统复杂度与重量增加、以及离心力下内衬完整性维持等显著工程挑战，更适用于固定储运单元。

在安全合规方面，全部转速工况（10–50 rad/s）均满足SAE J2601的85 °C温度上限要求。研究确定10–50 rad/s为关键操作窗口：低于此范围则对流混合不足，高于50 rad/s则机械应力与能耗增加而无对应冷却效率增益。

研究结论指出，旋转加注相较于传统静止框架加注具有本质优势，可有效避免入口射流冲击与局部过热问题，但其成功局限于特定速度范围。针对所考察的大规模储罐设置与加注条件，30 rad/s被确认为"最佳平衡点"（sweet spot），在温度均匀性、峰值温度限制与粘性耗散控制之间实现最优折中。该研究充分验证了旋转混合作为一种无需硬件改动的主动热管理方案在高压氢储存中的前景，同时为未来涡旋生成入口、内部低速转子或振荡机制的具体硬件设计提供了目标性能参数（速度、TKE、耗散率）。