Effect of bipolar-plates design on corrosion,mass and heat transfer in proton-exchange membrane fuel cells and water electrolyzers:A review

Attaining a decarbonized and sustainable energy system,which is the core solution to global energy issues,is accessible through the development of hydrogen energy.Proton-exchange membrane water electrolyzers(PEMWEs)are promising devices for hydrogen production,given their high efficiency,rapid responsiveness,and compactness.Bipolar plates account for a relatively high percentage of the total cost and weight compared with other components of PEMWEs.Thus,optimization of their design may accelerate the promotion of PEMWEs.This paper reviews the advances in materials and flow-field design for bipolar plates.First,the working conditions of proton-exchange membrane fuel cells(PEMFCs)and PEMWEs are compared,including reaction direction,operating temperature,pressure,input/output,and potential.Then,the current research status of bipolar-plate substrates and surface coatings is summarized,and some typical channel-rib flow fields and porous flow fields are presented.Furthermore,the effects of materials on mass and heat transfer and the possibility of reducing corrosion by improving the flow field structure are explored.Finally,this review discusses the potential directions of the development of bipolar-plate design,including material fabrication,flow-field geometry optimization using threedimensional printing,and surface-coating composition optimization based on computational materials science.

Numerical Investigations on Fluid Flow and Heat Transfer Characteristics of an Ultra-Thin Heat Pipe with Separated Wick Structures

Thermal and fluid-flow characteristics were numerically analyzed for ultra-thin heat pipes.Many studies have been conducted for ultra-thin heat pipes with a centered wick structure,but this study focused on separated wick structures to increase the evaporation/condensation surface areas within the heat pipe and to reduce the concentration of heat flux within the wick structure.A mathematical heat-pipe model was made in the threedimensional coordinate system,and the model consisted of three regions:a vapor channel,liquid-wick,and container wall regions.The conservation equations for mass,momentum,and energy were solved numerically with boundary conditions by using a code developed by one of the authors.The numerical results with the separated wick structures were compared with those with the centered,which confirmed the effectiveness of the separation of the wick structure.However,the effectiveness of the separation was affected by the position of the separated wick structure.A simple equation was presented to determine the optimum position of the separated wick structures.Numerical analyses were also conducted when the width of the heat pipe was increased with the cooled section,which clarified that the increase in the cooled-section width with the addition of wick structures wasmore effective than the increase in the cooled-section length.A 44%reduction in the total temperature difference of the heat pipe was obtained under the present numerical conditions.Furthermore,a comparison wasmade between experimental results and numerical results.

Analysis of Heat Transfer Performance of Oscillating Heat Pipes Based on a Central Composite Design

Oscillating heat pipes (OHPs) are very promising cooling devices. Their heat transfer performance is af- fected by many factors, and the form of the relationship between the performance and the factors is complex and non-linear. In this paper, the effects of charging ratio, inclination angle, and heat input and their interaction effects on heat transfer performance of a looped copper-water OHP are analyzed. First, suppose that the relationship between the response and the variables approximates a second-order model. And use the central composite design to arrange the ex- periment. Then, the method of least squares is used to estimate the parameters in the second-order model. Finally, multi- variate variance analysis is used to analyze the model. The results show that the assumption is right, that is to say, the re- lationship is well modeled by a second-order function. Among the three main effect variables, the effect of inclination angle is the most significant, but their interaction effects are not significant. In the range of the considered factors, both the optimum charging ratio and the optimum inclination angle increase as the heating water flow rate increases.

Evaporation Heat Transfer Characteristics from a Sintered Powder Wick Structure Sandwiched between Two Solid Walls

An ultra-thin flattened heat pipe has been developed with a centered wick structure.This structure is essential to make the heat pipe thinner.However,the centered wick structure reduces the evaporation and condensation surface areas of the wick structure because it is sandwiched between heat pipe walls.In this study,because detailed discussion has not been made,heat transfer experiments were conducted for the wick structure sandwiched between two solid walls.This study focused on the evaporation heat transfer characteristics from the sandwiched wick structure.The experiments were conducted with three wick structures,that is,strip-shaped sintered copper powders with thicknesses of 0.5,1.0,and 1.5 mm.Water was used as working fluid.The capillary pumping performance,that is,the liquid lifting velocities of the three wick structures were the same.The experimental results of the three wick structures were compared regarding the relation between the evaporation heat transfer rate and the superheat of the working fluid.The heat transfer experiments were also conducted when one of the solid walls was removed from the wick structure.It was confirmed that even if the wick structure was sandwiched between the solid walls,sufficient evaporation of the working fluid occurred from the thin sides of the wick structure.

Comprehensive Molten Salt Storage Shell and Supporting Structure Design and Analysis-Part I:A Conductive and Convective Theoretical Heat Transfer Analysis for Molten Salt Cylindrical Shells at 565℃ and 700℃

In this paper a full theoretical thermal analysis of a large molten salt container,80-foot in diameter and 46-foot high,including a four-foot elliptic shell roof,is presented for two temperatures,the standard 565℃ and a futuristic 700℃,which substantially improves the efficiency of the molten salt containers through the use of a highly stable chloride salt called SS700(SaltStream 700).The theoretical analysis includes conductive and convective heat transfer analysis in the steel container,elliptic roof shell,the fiberglass insulation,and firebrick insulation,and includes thermal insulation designs to safeguard against energy losses at high temperatures.The underlying soil and the high temperature concrete foundation were analyzed theoretically using conductive heat transfer,however the area surrounding the soil surface around the bottom of the molten salt storage tank had convective heat transfer analysis included.The final designs presented in this paper seek to limit heat losses to a maximum of 250 W/m^(2) while being able to operate at a minimum external ambient temperature of-10℃,which determines the thicknesses of the fiberglass and firebrick insulation.

Conceptual design and heat transfer performance of a flat-tile water-cooled divertor target

The divertor target components for the Chinese fusion engineering test reactor(CFETR)and the future experimental advanced superconducting tokamak(EAST)need to remove a heat flux of up to20 MW m-2.In view of such a high heat flux removal requirement,this study proposes a conceptual design for a flat-tile divertor target based on explosive welding and brazing technology.Rectangular water-cooled channels with a special thermal transfer structure(TTS)are designed in the heat sink to improve the flat-tile divertor target’s heat transfer performance(HTP).The parametric design and optimization methods are applied to study the influence of the TTS variation parameters,including height(H),width(W*),thickness(T),and spacing(L),on the HTP.The research results show that the flat-tile divertor target’s HTP is sensitive to the TTS parameter changes,and the sensitivity is T>L>W*>H.The HTP first increases and then decreases with the increase of T,L,and W*and gradually increases with the increase of H.The optimal design parameters are as follows:H=5.5 mm,W*=25.8 mm,T=2.2 mm,and L=9.7 mm.The HTP of the optimized flat-tile divertor target at different flow speeds and tungsten tile thicknesses is studied using the numerical simulation method.A flat-tile divertor mock-up is developed according to the optimized parameters.In addition,high heat flux(HHF)tests are performed on an electron beam facility to further investigate the mock-up HTP.The numerical simulation calculation results show that the optimized flat-tile divertor target has great potential for handling the steady-state heat load of 20 MW m-2under the tungsten tile thickness<5 mm and the flow speed7 m s^(-1).The heat transfer efficiency of the flat-tile divertor target with rectangular cooling channels improves by13%and30%compared to that of the flat-tile divertor target with circular cooling channels and the ITER-like monoblock,respectively.The HHF tests indicate that the flat-tile divertor mock-up can successfully withstand 1000 cycles of20 MW m-2of heat load without visible deformation,damage,and HTP degradation.The surface temperature of the flat-tile divertor mock-up at the 1000th cycle is only930℃.The flat-tile divertor target’s HTP is greatly improved by the parametric design and optimization method,and is better than the ITER-like monoblock and the flat-tile mock-up for the WEST divertor.This conceptual design is currently being applied to the engineering design of the CFETR and EAST flat-tile divertors.

Optimal Experiment Design for the Identification of the Interfacial Heat Transfer Coefficient in Sand Casting

The interfacial heat transfer coefficient(IHTC)is one of the main input parameters required by casting simulation software.It plays an important role in the accurate modeling of the solidification process.However,its value is not easily identifiable by means of experimental methods requiring temperature measurements during the solidification process itself.For these reasons,an optimal experiment design was performed in this study to determine the optimal position for the temperature measurement and the optimal thickness of the rectangular cast iron part.This parameter was identified using an inverse technique.In particular,two different algorithms were used:Levenberg Marquard(LM)and Monte Carlo(MC).A numerical model of the solidification process was associated with the optimization algorithm.The temperature was measured at different positions from the mould/metal interface at d=0 mm(mould/metal interface),30 mm,60 mm and 90 mm.the thicknesses of the cast part were:L1=40 mm,60 mm and 80 mm.A comparative study on the IHTC identification was then carried out by varying the initial value of the IHTC between 500 Wm^(-2)K^(-1) and 1050 Wm^(-2)K^(-1).Results showed that the MC algorithm used for estimating the IHTC gives the best results,and the optimal position was at d=30 mm,the position closest to the mould/metal interface,for the lowest thickness L1=40 mm.

A Comprehensive Review on Heat Transfer Design Parameters in MEMS through Experimentations and Numerical Modeling

Micro electro mechanical systems are studied in this paper along to comparing the design, fabrication, and operation of such devices especially in the electronic industry. The advancements of such systems are important as scale reduction is becoming more economical and also in heat dissipation efficiencies. Furthermore, the fabrication techniques are becoming more of a challenge in today’s industrial market. In this paper we will introduce a comparative study of experimental research in the micro scale heat exchanger with varieties of industrial application primarily cooling integrated circuits. Significant effect in the process of cooling electronic devices has showed experimentally that scale variations impact into this improvement. Numerical studies have been a leading technique in the development of thermal systems especially at the micro or small scales which helps the industry to cut cost and provide innovative products at high efficiencies. There are so many approaches to numerical analysis and research, where the availability of ready to use software and data acquisition modules is variant. Material selection is also a variable in the efficiency and operations of these micro channel thermal system and heat exchangers.

Comprehensive Molten Salt Storage Shell and Supporting Structure Design and Analysis-Part II:A Conductive and Convective Numerical Finite Element Heat Transfer Analysis for Molten Salt Cylindrical Shells at 700℃,and Comparison with Theoretical Analysis

In this paper a finite element thermal analysis model-using COMSOL-of a large molten salt container,80-foot in diameter and 46-foot high that includes a four-foot elliptic shell roof,is presented for a futuristic 700℃ design,which uses a highly stable chloride salt called SS700(SaltStream 700)that improves the efficiency of the tank when compared to the traditional 565℃.The FEA(finite element analysis)includes conductive and convective heat transfer analysis in the steel container,elliptic roof shell,the fiberglass insulation,and firebrick insulation,and includes thermal insulation designs to safeguard against energy losses at high temperatures.The underlying soil and the high temperature concrete foundation were analyzed by finite element using conductive heat transfer,however the area surrounding the soil surface around the bottom of the MS storage tank had convective heat transfer analysis included.The finite elements analyses presented are to verify the final fiberglass and firebrick insulation designs,which seeks to limit heat losses to a maximum of 250 W/m^(2) while being able to operate at a minimum external ambient temperature of-10℃.These results are also compared to previously calculated theoretical results.

Optimal design of the separate type heat pipe heat exchanger

Separate type heat pipe heat exchangers are often used for large-scale heat exchanging. The arrangement of such a heat exchanger conveniently allows heat input to and output from the heat exchanger at remote locations. The traditional method of designing an ordinary HPHE (heat pipe heat exchanger) is commonly applied in the separate type exchanger design, but the calculations have to be carried out separately, which makes it very complicated. In this work, the ε-NTU (effectiveness-Number of Transfer Units) method was applied for optimization analysis of single- or multi-level separate type heat pipe heat exchangers. An optimizing formula for single-level separate type heat pipe heat exchangers was obtained. The optimizing principles of effec- tiveness-NTU and heat transfer rate by the equal distribution method for multi-level separate type heat pipe heat exchanger are presented. The design of separate type heat pipe heat exchangers by the optimizing method is more convenient and faster than by the traditional method.

基于场协同理论的热空气吹掠型太阳能蒸馏器对流传热性能优化分析

以热空气吹掠型太阳能蒸馏器为研究对象,利用数值模拟结果和均匀设计法分析蒸馏器内流动与传热场协同角的变化规律和影响因素。结果表明,影响场协同角的主要因素是热空气进口速度,热空气进口温度和水温的影响不显著。在蒸馏水温30~70℃、进口空气温度40~90℃,进口空气速度为8 m/s时平均传热场协同角最小,为73.14°,传热系数最大,可达到同等条件下自然对流传热系数的4倍。

螺旋半圆管夹套式余热排出装置传热分析及设计程序开发

提出一种适用于小型反应堆的螺旋半圆管夹套式非能动余热排出装置,利用螺旋半圆管夹套具有的传热效率高、结构紧凑等优点,与反应堆主容器形成一体化的紧凑设计。基于此,对螺旋半圆管夹套的传热过程进行了分析,提出传热计算模型并基于此开发了热工设计程序,通过与公开文献中试验数据的对比分析,验证了程序的正确性,可以用于半圆管夹套式余热排出装置的工程设计。

面向金属泡沫散热器设计的自然对流拓扑优化

针对传统方块型金属泡沫散热器散热效率不高,难以满足高功率电子器件散热需求的问题,发展了基于格子Boltzmann和水平集方法的三维拓扑优化方法。首先,采用格子Boltzmann方法求解多孔介质内的流动和传热过程;然后,采用伴随格子Boltzmann方法计算梯度;最后,根据反应-扩散方程更新水平集函数。通过对底部放置发热元件方腔的自然对流过程开展拓扑优化,得到特征格拉晓夫数为2.4×10^(2)~1.2×10^(5)的4种优化散热器结构,定量评估了散热器的强化换热性能并分析相关机理。研究结果表明:随着格拉晓夫数的增大,散热主导机制由热传导转变为热对流,优化结构由伸展的树枝状向收缩的花朵状结构转变,以留出更多空间使得流动涡旋充分发展;与传统的金属泡沫方块结构和开槽的金属泡沫块结构相比,所提优化结构表现出优异的散热性能,其分支结构和中空结构能够优化方腔内的流动,从而将散热效率提升了29.7%以上,表明了所提拓扑优化方法的有效性。研究工作可为新型金属泡沫散热器的设计提供理论指导。

汽轮机蒸汽管道与阀门及汽缸外缸外表面传热系数的计算模型和计算方法

构建了汽轮机蒸汽管道、主汽调节阀与汽缸外缸的外表面传热系数计算模型与计算方法。给出了汽轮机蒸汽管道、主汽调节阀、高压外缸、中压外缸与低压外缸的内表面传热系数以及保温结构外表面与低压外缸外表面的复合传热系数计算公式。建立了汽轮机蒸汽管道、高压外缸与中压外缸的金属内层壁与保温结构外层壁的双层圆筒壁、主汽调节阀金属内层壁与保温结构外层壁的双层球壁以及低压外缸单层圆筒壁的传热过程计算模型。考虑到金属内层壁内表面温度与外表面温度和保温结构外层壁外表面温度均为待定温度,采用迭代法确定这些汽轮机部件的表面温度、传热过程的传热系数与热流密度。完成了汽轮机蒸汽管道、主汽调节阀、高压外缸、中压外缸与低压外缸的外表面传热系数计算分析的应用实例,得出了这些汽轮机部件传热过程金属内层壁外表面的传热系数、热流密度与等效表面传热系数的计算结果,为这些汽轮机部件温度场与应力场的有限元数值计算提供了传热边界条件。

增程式电动汽车动力系统一体化冷却传热分析

以一款增程式轻型卡车为研究对象,对增程式动力系统的发动机和电机采用一体化冷却系统设计。在预设的运行工况,通过试验测试,计算出工况点增程器和电驱系统的散热需求,对冷却液和冷却水泵进行选型匹配。结合试验测试建立了增程器和电驱系统冷却回路一维仿真模型,对比分析了电驱系统冷却管道串联与并联的冷却效果与能耗情况,验证了驱动电机冷却液帮助发动机冷起动的可行性。研究结果表明:通过试验测量散热量选择到的电子水泵,能更加高效准确地控制冷却液流量,冷却效果更好。冷却系统一维仿真分析中能从温升曲线直观看到并联冷却管路冷却效果更好,并联管路中的电磁阀根据工况不同控制开口大小,比串联管路统一控制冷却液流量冷却更高效。并联冷却管路中驱动电机冷却时快速产生热量,在与发动机冷起动时进行热交互可以帮助发动机快速暖机,改善发动机冷起动排放差、动力弱的情况。

变流量超临界CO_(2)传热特性实验设计

随着我国“双碳”目标的确立,二氧化碳作为具有较优环保性能的自然工质,是制冷剂替代工作中被广泛关注的研究对象。搭建了变流量超临界二氧化碳(CO_(2))传热特性实验台,对不同入口压力、质量通量、管径下水平微通道管内换热特性进行了探究。结果表明,超临界CO_(2)换热系数随压力的升高而减小,随着质量通量的增大而增大,随着管径的增大而减小。分析可知,影响换热的主要因素是物性的变化,换热系数由于在临界点附近物性的突变存在峰值,实验中获得了不同参数对超临界CO_(2)换热的影响规律。本实验装置和实验内容能够较好地用于相关专业的实验教学,也为超临界CO_(2)换热设备和冷热系统的设计提供指导。

夏热冬冷地区动态传热保温设计室外典型气象参数确定方法

GB50176-2016《民用建筑热工设计规范》中围护结构外墙保温设计采用稳态方法,但夏热冬冷地区由于特有的建筑运行模式,使外墙传热方向发生逆转,故稳态方法难以满足冬季室内热环境水平提升的实际需求。提出一种基于动态传热的冬季外墙保温设计室外计算温度确定方法,探究室外干球温度和太阳辐射逐时波动的最不利组合,将保温设计的目标由“日”提升到“时”。基于典型城市2007-2017年逐时气象数据,以外墙内表面温度的累计不保证天数作为热工设计指标,将累计不保证天数恰好不为0(指各朝向外墙内表面温度低于室内基本热舒适条件下所能允许温差的天数)时作为典型日挑选的依据,并将典型日出现频率最高的日期作为室外计算温度选取的依据,获得冬季保温设计典型日,并获得了典型日计算时长的影响。可为夏热冬冷地区围护结构外墙保温设计室外计算参数的挑选提供理论支撑。

高温环路热管设计方法研究与传热极限分析

高超声速飞行器在以高超声速飞行时尾喷管会面临严酷的燃烧热考验,结合高温热管导热性强和环路热管传热距离长等优势,高温环路热管可以作为尾喷管热防护的新方法.根据液态钠高导热系数及高表面张力的物性设计了一种新型蒸发器,提出了一套高温环路热管的设计流程.使用这套设计流程设计了一根高温环路热管,结合热管传热极限的计算方法,对钠热管进行仿真计算,分析了温度和结构对热管传热极限的影响.研究发现在高温环路热管的运行温区内,热管的传热性能随着温度升高一般由声速极限、毛细极限和沸腾极限依次决定.结合研究结论提出了热管结构优化方法,为热管设计提供了理论依据.

基于Fluent流体传热仿真的毛细管风道设计优化

针对毛细管冷却风道内分离段毛细管柱轴向散热不均匀问题,设计一种变截面送风道。利用Fluent仿真软件对风道设计过程进行建模和数值模拟,优化改进风道冲击角、出入口形状与风道长度尺寸等参数,并通过对风道出风口风速与风道内柱面温度的实测验证设计的可靠性。结果表明:毛细管冷却风道经优化后,出风口的风速更加均匀,柱面轴向温差减小,温度均匀性得到改善。仿真结果较好地拟合了实验结果,4个出风口的风速在10.38~10.67m·s^(-1),柱面温差缩小至±1℃以内。优化后的毛细管冷却风道设计对提升分离段冷却效率具有指导意义。

基于Kriging模型的螺旋波纹管流动换热特性及结构优化

在螺旋管式换热器的基础上结合波纹结构,设计出一种新型的螺旋波纹管换热器.通过数值计算研究了不同波纹深度和波纹节距对其流动换热的影响.以数值计算结果训练Kriging模型,训练后的Kriging模型可预测大量不同螺旋波纹结构特征样本的流动换热特性.通过筛选摩擦系数较低且换热效率较高的结构参数以实现螺旋波纹管换热特性优化.结果表明,波纹结构产生的流速突变和螺旋结构产生的二次流提高了换热效率和流动阻力.通过与数值计算结果的对比,Kriging模型预测换热特性具有较高的效率和精度.在所选定的参数范围内,当波纹深度、波纹节距、螺旋中径、螺旋节距分别为1.32,27.99,86.49,57.85 mm时,螺旋波纹管具有较低的摩擦系数和较高的换热效率.