室温磁致冷的仿真与解决方案.doc

室温磁致冷的仿真与分析 ? ? 摘　要:室温磁致冷技术尚处于样机研制阶段,仿真分 析可以简化该过程。借助有限元分析工具计算磁场分布,并将工质细化然后结合磁工质材料磁热特性计算工质微元完整循环中的热量变化,最后求和近似得到磁致冷 装置的整体运行特性。针对一种简化室温磁致冷原理样机的Ericsson循环制冷过程进行详细的计算。仿真结果表明提高循环的制冷温跨可增大总回热量,但 是制冷效率将会降低。该方法有助于综合评定室温磁致冷装置的设计效果、指导室温磁致冷机的研制过程。 ? ??室温磁致冷由于其高效和环保,具有广阔的发展前景,目前这项研究工作主要集中于新的磁工质材料和磁致冷装置研制。由于可能应用于室温磁致冷的材料的磁 热效应有限,必须借助较高强度的外磁场变化才可能获得较好的制冷效果,而强磁场源的设计实现非常困难,因此距实用化尚有一定距离。 ? ? 国内外对各种磁热力循环、特别是Ericsson循环性能的研究做了大量的工作[1,2],本文在总结这些研究成果的基础上,提出了综合磁场、磁热效应、 磁致冷循环和传热等多个环节对室温磁致冷过程进行全面仿真的分析方法,然后对一个改进设计的室温磁致冷装置进行了计算分析,验证了该方法的可行性和有效 性。 ? ? 1　室温磁致冷技术仿真分析方法原理 ? ? 与传统气体压缩制冷过程相比较,磁致冷的最大特点是其所采用的制冷工质为固体磁性材料,这决定了磁致冷装置在结构上与气体压缩制冷机有着根本区别。通常, 室温磁致冷装置主要由磁工质、磁场发生源、传热流体循环单元、蓄冷器以及冷热源热交换器等部件组成,对磁致冷装置的仿真分析必须从这些主要模块入手。磁工 质是磁致冷装置中完成吸热和散热作用的核心单元,对磁工质在循环过程中的各个状态变化进行计算是室温磁致冷仿真分析的核心问题。由于工质磁热效应的发挥受 到磁场源产生的磁场分布影响,计算中须将磁工质单元细化成微元,然后根据材料的磁热特性和场强分布情况计算工质微元在循环过程中的热量变化,在此基础上求 和即得整个磁工质单元的制冷效果和性能。 ? ? 1.1　磁场区域场强分布计算 ? ? 实际过程中无论磁致冷装置的磁场源设计如何优化,所产生的磁场分布都无法达到理想的均匀分布;处于工作磁场区域内的磁工质在循环过程的有限时间内,由外磁 场改变而引起状态和热量的变化也不能达到完全均一,必须按磁场分布做相应计算才能获得对整个磁工质工作状况的估计。磁场分布的计算可以借助有限元计算工具 如ANSYS来完成。 对于含有恒磁场过程的循环方式来说,其磁化或去磁过程可利用绝热或等温的条件进行分析,又由于通常情况下循环的较低恒磁场过程即取零磁场状态近似,所以室 温磁致冷的分析中真正需要做磁场分布分析的仅是较高恒磁场的过程。计算后得到工质微元所处磁场的分布情况及网格剖分结果,按磁工质微元的顺序表示即为: ? ? 式中:H2表示网格剖分后所有微元的平均场强,H2(k)表示每个磁工质微元上的平均场强,ρ表示所选用的磁工质材料的密度,v(k)为网格剖分后得到的各个子单元的体积参数,m表示网格剖分后所有微元的质量。 ? ? 1.2　磁工质材料磁热特性分析 磁性材料磁热特性的研究不仅是选择磁工质的重要依据,还是室温磁致冷技术仿真分析的首要前提。恒场热容CH,等温磁熵变ΔSM和绝热温度变化ΔTad是材料磁热特性的表征方法[3],均为磁场和温度的函数。通过测量及计算可得到材料磁特性的离散值: ? ? 式中Hi和Tj分别为测量恒场热容时的外磁场强和温度。对CH、ΔSM和ΔTad的离散值采用曲线拟合和插值的方法可近似得到在一定温度和磁场强度变化范 围内随温度连续变化的磁热效应曲线,以便于后面的分析计算。恒磁场下,材料的熵还包含有晶格熵SL和电子熵SE[4],室温时材料晶格熵较大,必须加以考 虑,利用由Pierre-Louis Dulong定律得到晶格比热容[5],即 cL= 24.94×103/M, (3) ? ? 式中M为材料的原子量或分子量。相比较而言,材料的电子熵较低[4],本文暂不计其作用效果。 ? ? 1.3　磁致冷循环过程的计算方法推导 ? ? 由于室温时材料晶格熵大到不能忽略,因而只有具有回热特性的磁Brayton、磁Ericsson和磁Stirling循环才适用于室温磁致冷[3]。分 析磁致冷循环须首先分析其过程的特点如恒场、绝热、恒温等以确定其循环方式,然后采用相应不同的计算方法。3种可行的室温磁致冷循环中, Ericsson循环虽不具有理想回热,但由于由两个等温和两个恒磁场过程组成,可获得相对较大的温跨,因而是一种重要的循环方式[1,2]。下面主要对 Ericsson循环进行工质热量变化计算方法的讨论。 ? ? 磁工质的最低和最高工作温度分别是Tc和Th。等温磁化过程ΔT=0,即晶格熵引