基于lsm的钻孔热阻并验证模型.docxVIP

基于lsm的钻孔热阻并验证模型 摘要：基于Comsol有限元软件建立垂直地埋管二维导热数值模型，在综合考虑支管半径、钻孔半径、两支管间距、土壤导热系数和回填材料导热系数的基础上进行仿真实验，得出不同热物性参数值下的钻孔热阻并验证模型的可靠性。基于田口方法，建立L 0 引言地源热泵系统因具备使用年限长、能源利用率高、环境友好等特性而发展前景广阔，但其初投资较高，发展速度因而受到限制本文在上述研究基础上，利用COMSOL Multiphysics 5.3a有限元软件建立钻孔内二维导热问题的数值模型，求得不同热物性参数值下的钻孔热阻值，并基于田口方法得出最优参数组合和参数贡献率，进而得出最优钻孔热阻值。1 钻孔二维导热数值模型1.1 二维传热模型本文建立钻孔内二维稳态传热模型，忽略钻孔的轴向导热支管内流体的对流换热以及支管的自身热阻，考虑在垂直于钻孔轴线的横截面上由于支管位置、大小、形状等因素的影响而引起的二维温度场。模型如图1所示。控制方程：钻孔内回填材料的导热微分方程式为：钻孔外土壤的导热微分方程式为：钻孔热阻的计算公式初始条件：土壤及回填材料初始温度：T边界条件：当土壤计算区域无限大时，边界区域受埋管的影响可忽略，在工程中，埋管间距为3～6 m，因此取距钻孔中心长度为6 m的圆边界作为数值模拟的远边界，且其温度值同土壤初始温度相等：进出水支管的外壁平均温度为内边界条件：T1.2 热阻仿真实验对传热模型进行网格划分，网格划分示意图如图2所示。文献将T1.3 网格无关性检验为了消除因网格和探针数量过少而使仿真结果出现偏差的情况，本文基于不同数量的网格和探针，进行仿真实验，热物性参数的设定值如表1所示。测试结果如图3所示。定义变动率为：当ROC0.05%时，可认为钻孔壁温的变化可忽略不计。由图3及式（6）可知，当设定的域点探针数量由70增至100时，钻孔壁温度趋于平缓，变动率ROC=0.027%；另一方面，当网格数量由654增至938时，钻孔壁温度无变动，变动率ROC=0，而运算时间会有所增加。因此，在仿真实验中，拟将整个区域划分为大小不等的654个网格，在钻孔壁上设定70个域点探针和边界探针。1.4 模型可靠性分析为了验证模型的可靠性，本文进行5次仿真实验，并将实验结果与学者Paul、Hellstrom、Bennet和Sharqawy得到的结果进行对比，每次仿真实验的热物性参数设定值如表2所示，各种方法得到的钻孔热阻值如图4所示。由图4可知，由本模型仿真得到的热阻值与由Paul、Hellstrom、Bennet和Sharqawy求得的计算值有相近的变化趋势。本文利用仿真值与计算值间的平均绝对误差MAE和均方根误差RMSE来确定仿真值与计算值之间的差异，进而验证数值模型的准确性：由式（8）和式（9）计算得出的RMSE值和MAE值如表3所示。由表3可知，仿真值和计算值之间的误差小、可靠性高。仿真值同Hellstrom的计算值偏差程度最小，RMSE和MAE值分别为0.13%和0.10%，其次是Sharqawy、Bennet和Paul。值得注意的是，仿真值同Paul的计算值偏差程度较大，RMSE和MAE值分别为1.13%和0.86%，这主要是因为Paul的钻孔热阻公式是基于实验求得的，热阻值的准确与否与实验次数有很大关联。2 田口方法2.1 求解信噪比田口方法能通过正交实验以及各种评价指标有效地评定各参数对目标指标的影响，并能给出最优参数组合影响钻孔热阻的5种热物性参数：（A）支管半径r通常，田口方法的最小实验次数N本文所选参数数目N表5列出了各参数在正交表中的排列位置以及每次数值模型得到的热阻值和相应的信噪比。2.2 确定最优参数组合对每次数值模型求得的信噪比进行求解，可得如表6所示的信噪比响应表。表6包含5种参数3类水平下的平均信噪比以及每个参数的最大平均响应值和最小平均响应值之差Delta，并根据Delta的大小进行相应排秩。通常来说，信噪比越大品质越好，而且参数的最大平均响应值和最小平均响应值的差值越大，说明该参数对结果的影响作用越大。由表6可知，钻孔半径、回填导热系数和支管半径的Delta值较大，说明这些参数对结果的影响作用较大。对比各参数各水平下的平均信噪比，应选择A2.3 参数的方差分析为了量化各参数对钻孔热阻值的影响程度，本文基于方差分析法求得各种参数的平方和、方差、P值和F值，进而求得各参数对钻孔热阻值的贡献率，用以确定参数的相对重要程度。信噪比方差分析数据如表7所示，F值是组间平均方差和组内平均方差的比，用以量化参数对对钻孔热阻值的的贡献程度。P值可确定钻孔热阻值与每个参