R410A和R407C热力性质简化计算

    制冷空调行业的各种探索和研究表明,混合工质在制冷工质替代中具有很大的潜力,其中R410A和R407C作为R22的替代物更是倍受瞩目.为了更好地研究它们对现有制冷系统的影响,计算机仿真是个很好的手段.而热物性程序作为仿真程序的基础部分,对仿真计算的效率和结果有相当的影响.但国内在这方面的研究很少,一般直接采用复杂的状态方程进行迭代计算,这样在相当程度上降低了仿真的速度和稳定性.为了弥补这一不足,本文采用隐式拟合显式计算的方法,参照DuPont公司的数据,对R410A和R407C的热力性质重新拟合,结果可以避免迭代,在显著提高计算速度的同时又能保证所需的精度.

    1　拟合模型

    本文对根据文献[1,2]编制的热力性质程序进行简化,并以该热力性质程序作为拟合的参考数据源和检验简化热力性质的相对精度.饱和热力性质的简化模型采用了文献[3]提出的拟合函数形式:

    式中:u=y-y0;v=x-x0.

    对于过热区的热力性质,文献[4]没能给出形式完全统一的拟合函数.作者统一了过热区简化模型的形式,以便拟合和降低模型的复杂性,具体形式:

    式中:Ts、zs为p所对应的饱和状态参数.

    简化模型是在常用的制冷空调运行工况内进行简化,在饱和区-40～60°C和过热区-40～120°C内保证精度,在温度外推20°C范围内保证变化趋势,以确保仿真计算的正确进行.

    2　隐式拟合的分岔问题和解决方法

    隐式拟合的最大问题就是分岔问题.在隐式拟合方程向显式的计算方程转化时,涉及到根的判别问题.三次方程涉及到3个根,分别代表了不同的根轨迹.但最后需要的可能是其中一条或多条根轨迹组合而成.然而根轨迹的衔接处会出现很小的断裂,这便是分岔现象.断裂处称为分岔点,它造成曲线不连续(分岔点处误差大)和曲线不光滑(分岔点处一阶导数不连续).虽然分岔点的范围很小,在大部分情况对制冷系统的仿真模型不会产生大的影响,但却是个巨大的隐患(可能导致仿真模型计算值的异常).

    通过一系列的尝试,在不改变模型的前提下提出一种解决方法:通过改变拟合数据来调整拟合过程,把分岔点移出拟合范围,同时保证拟合精度.这一方法的数学原理是通过改变拟合数据点可以改变拟合函数的曲率.因此,只要令拟合函数在拟合范围内曲率变化减小,就可使分岔点(即曲率变化最大的点)外移.通过这种方法建立的模型在拟合范围内没有分岔问题,少数模型在外推范围内有分岔,但这对常见制冷空调工况范围内的系统仿真没有影响.具体方法:①改变拟合范围,通常是扩大拟合范围以保证拟合范围内的精度;②改变拟合的点数,大部分情况是减少点数;③用非均匀的数据点拟合(增加某区域内数据个数).相比之下,① 对分岔点的位置影响最大,③则最小.从目前情况看,拟合数据范围、点数和分布的选择对不同的热力参数是不同的,在很大程度上取决于经验.