制冷专业毕业设计外文翻译.doc

《制冷专业毕业设计外文翻译.doc》由会员分享，可在线阅读，更多相关《制冷专业毕业设计外文翻译.doc（11页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、中文译文：结霜工况下强制对流换热器的翅片性能摘要 在本文中，经过修改的模拟结霜工况下风机供风翅片式换热器的热交换翅片上霜增长的数值模型已经得到证实。结果发现，制冷换热器上结霜导致翅片热流量、气流速度和翅效率的急剧下降，以及压降的增大。精确地研究显示改变几个设计参数（其中包括风机的型式）的影响。2002年Eldevier科学有限公司和IIR。保留所有权利。1.介绍 对制冷系统、空调设备以及空气对空气型热回收系统来说，翅片式热交换器结霜是一个常见问题。在这些系统中，结霜对制冷换热器来说是最重要的性能影响因素。换热器表面结霜对热交换器的热工性能的影响有好几个方面。 1. 换热器表面结霜增大了翅片与气

2、流间的热阻，从而降低了制冷系统用热交换器的冷却能力。对于气流速度以及气流与翅片表面的温差给定的情况，一定厚度的霜层会降低热流量，这是因为霜层的实际导热系数很小并且霜层间会形成温差。 2. 热交换器表面结霜使得通过换热器的气体量大大减少、气体压降大大上升。根据风机的特性，在连续结霜数小时后，气流的路径会变短甚至会被完全堵塞。主要符号 总表面积， 基管表面积， 最小流通面面积， 翅片面积， 迎风面面积， 压力损失修正系数 摩擦系数 翅的比热容， 霜的比定压热容， 空气的比定压热容， 霜层实际气体扩散速度， 霜表面实际气体速度， 管内径， 穆迪（或达西）摩擦系数 翅高， 对流传热系数， 对流传质系数

3、， 升华潜热， 无霜对流传热系数， 入口流体收缩损失系数 出口气体膨胀损失系数 霜的导热系数， 翅片的导热系数， 翅长， 相变速度， 翅片数 充分发展流的努塞尔数 局部努塞尔数 大气压强， 普朗特数 气体体积流量， 光管热流量， 单翅片热流量， 基面总热流量， 结霜下单翅片热流量, 基面总传热量变化， 以内径为特征尺度的雷诺数 当地舍伍德数 方程15、16的源项， 时间， 翅间距， 温度， 翅间气流平均温度，K 翅表面温度，K 翅间平均气流温度， 基面温度， 入口气流温度， 出口气流温度， 翅间气体平均速度， 迎风面风速， 入口气体湿度， 湿度， 翅片间气体平均湿度， 霜表面湿度， 翅片入口气

4、体含湿量， 坐标， 翅片出口气体含湿量， 坐标， 霜内空气体积分数， 坐标， 霜内冰的体积分数， 霜层厚度， 平均测量霜层厚度， 翅片霜层平均厚度， 无翅区霜层平均厚度， 翅厚， 空气动力粘度， 翅效率 空气密度， 霜密度， 平均测量霜密度， 翅片密度， 冰密度， 蒸汽密度， 翅片压降， 测得翅间平均压降 空气经翅片的温变， 霜层厚度变化， 沿翅片空气的湿度变化， 霜的密度变化， 沿翅片压降的变化， 通过基面的总热流量， 3. 为了维持积霜换热器的良好性能，定期或循环融霜来防止霜的累积是必需的。通常情况下，在除霜周期内，制冷系统将被关闭，同时热量将热交换器加热到冰点以上。Niederer 指出

5、工质携带的热量中只有15至20的用于融霜。其余的热量有的提高了换热器的温度，有的则直接散失到周围的环境中。在融霜周期中蒸发的水通常又会在下一个结霜时间凝结回到热交换器表面。 与无霜工况相比，由于结霜的缘故，相同的冷量下，目前工业用制冷系统的尺寸要大50，同时同一个系统平均要多输入25的能量。 制冷换热器设计人员需要知道霜的累积速度，以及它与气流阻塞、压力损失、热效率的关系。对翅片管换热器结霜过程的研究表明，结霜过程随空间与时间有较大的变化，这些变化由换热器的设计参数与运行工况而决定。在过去超过50年间，发表的关于换热器表面结霜的研究论文近150篇。ONeal、Tree和Padki总结了许多的这

6、类文献。以前的研究大多不涉及表面结霜模型的发展，而是致力于金属平板表面结霜的简单的不变的特性模型，或是监测典型换热器的结霜特性。很少有研究涉及结霜时的压降与气流变化。过去关于热交换器与传热表面结霜的研究为结霜热交换器的设计起到了指导作用，但是不能为设计效果提供详细预测。没有任何一个研究是关于发展翅片表面精确结霜模型的。 数字模型是很重要的，因为它们不仅能模拟一个广泛的运行工况，而且能很好的解释实践中很普遍但实验研究中不切实际的时间因素对运行工况的影响。一经验证，数字模型对于热交换器运行工况设计与优化很有作用。 最近，学者们对结霜工况下热交换器得性能做了更多的研究。一个用于预测热交换器翅片结霜的

7、数学模型已经在Chen等人的研究中被提到。该模型已经经过实验数据证实并且在Thonas与Chen等人发表的文献中也被提到。在这些实验与数值研究中，每次实验时都保持气流速度恒定。这种恒定气流速度的条件在风机驱动的热交换器中是不可能的。 本文中，对于相同布置的热交换器，由Chen发展的数值模型经过修改后比以前的更加符合事实情况。因此，正如Chen的研究中所讨论的那样，这种模拟结果的精确度更高。根据对一种典型风机曲线的研究，恒定的气流经过热交换器后变得各式各样。从而，这种对结霜工况下热交换器特性的模拟结果仍有待保留与商榷。2.数值模拟2.1问题陈述与分析图1显示了系统布置简图。风机布置在热交换器翅侧

8、表面之前，驱动气流穿过结霜的热交换器。图2为低温环境中冷基面（即 40左右）、供冷空气即 20左右以及高相对湿度即 90左右的热交换器翅侧的简明图。对于给定的运行工况与翅的几何尺寸，数值模拟不仅预测了霜的瞬时的累积、翅平面与冷基平面的霜高分布，而且还预测了沿冷基面的热效率及沿翅的气体压力降。表1罗列了一种典型运行工况的参数和翅的几何参数，也就是输入值或模型已知值。 该处呈现的热交换器的布置与翅片管式热交换器的不同，这是因为翅片安装在冷平板基面（如图2所示）。正如Thomas在论文中解释的那样，对于每一个热交换器翅上霜的高度高确定性的的实验测量，这种型式的布置最为方便。这是唯一一种现存的已经证实

9、的结霜数值模型的翅片热交换器布置方式。为了适应与验证复杂几何形状翅片管式热交换器的模型，就需要进行更多的实验与理论模拟研究。同时，验证过的模型能够用于显示结霜工况下典型热交换器的一些重要的特性。 如图2所示，模拟平翅片热交换器的数值模型包括两部分：一个是用于描述翅化表面与非翅化基面上结霜情况的结霜模型；另一个是用于描述每一个铝翅片传热情况的二维热传导模型。读者可以参考附录5、8来了解该模型的所有细节，以及用于该模型的精确实验数据验证、不稳定性数据与特性的分析讨论。该模型的关键部分在于模拟了流动区域。2.2结霜模型 应用于寒冷翅片上每一点的结霜模型都是一维的，并模拟了霜层内传热传质的瞬态情况。结

10、霜模型的核心方程与边界条件是： 能量方程： （1） 水蒸气扩散方程： （2）冰相连续型方程： （3） 霜与气流界面的传质边界条件： （4） 传热边界条件： （5） （6） 式中，是来自于翅的二维热传导模型的翅表面温度。2.3霜与气流间的传热传质 Gnielinski相关用来评估对流传热传质系数，与。为了用Gnielinski相关进一步简化与的计算，假设翅表面霜层增长均匀，即霜层高度相等。同时还假设与只沿气流的y方向变化。 对流传热的Gnielinski相关适用于管内雷诺数在3000到10000之间的管内充分发展紊流。 （7） 雷诺数为： （8）式中，翅间霜自由区域的水力直径。 的计算： （9）

11、平均气流速度的计算： （10） 翅前来流速度。 （11）式中，为翅的来流横截面积。 在Kays与Crawford的论文中，图14-7就是用来确定紊流气流流过翅片区而引起的当地努塞尔数的增强程度。假定随传热传质过程的缓慢变化，传热传质的变化趋势相似 （12）式中，是当地舍伍德数。 因此，与是从与的定义中得到的。 （13） （14）2.4翅的导热 由边界条件可以看出，霜层内传热的模型与翅的导热模型是相似的。翅内传热的核心方程为， （15）式中，是源项，由结霜模型参数决定，例如从霜层到翅表面的热流量。 （16）2.5压降与热流量的计算 由热交换器翅表面结霜而引起的堵塞，会导致沿翅片的气体压降增大。假

12、定气体密度恒定，沿翅的气体压降可由下式计算， （17）式中，是考虑由结霜引起的入口与出口压力损失的实验观测系数。沿每一个结霜翅片的热流密度由下式计算， （18）翅间基平面的热流量由结霜模型得出，由下式计算， （19）式中，是非翅化区基平面总面积。那么，总热流密度为： （20）此处，为翅片总数。2.6气体的平均温度和湿度的计算 方程（4）和（5）中气流温度和湿度定义为：翅化区气体平均温度与湿度。也就是： （21） （22）对翅化区用能量与质量守恒， （23） （24）此处， 空气与翅间的总热流量 霜表面湿度 随着翅片间气流温度的降低，根据方程(21)-(24)的计算湿度的方法，表明气体会达到与气

13、流温度（即伴随着空气中的小水滴或冰晶，相对湿度达到100）相对应的超饱和状态。气流通过翅间的时间非常短（约0.02秒），因此若仅考虑热交换器上的超饱和空气，气流可望保持0.02秒内不在微小尘埃颗粒上凝结与升华。 正如Mao等在论文（附录12）及Chen等在论文（附录13）中所观测的一种特殊情况，即当进入翅间的气源为含霜晶和雾滴的超饱和状态，该状态气体的结霜过程与水蒸气的明显不一样。当气源内含有冰晶时，结霜过程主要受扩散的霜晶以及它们在霜层外边缘积累的影响，这种霜层外边缘除了霜层内水蒸气的扩散外没有其他的扩散。而我们的霜层扩散模型并不适应与这种特殊的物理现象。在这种超饱和情况中，依赖于空气扩散动

14、力的结霜预计进行的将非常的快并且会非常粗糙。因此为了发展超饱和气源结霜的数字模型，研究者们需要进行进一步的实验与理论研究工作。2.7恒定气流速度模型的验证 Chen等通过精确地数据把这种模型的霜层高度分布与一种典型翅片的进行了比较，如恒定流气流速度模型的沿流动方向的霜冻密度的分布，翅面的总热流量及气流的压降等。这些模拟结果与实测数据比较得出的平均差异列于表2。表中的这些数据来源于一系列不同运行工况（如不同的霜厚，霜密度，总热流量，沿翅的压降）的大量的实验数据。 上述数据唯一的不确定性源于实验的平均偏差以及平均模拟数据与测量数据的平均差异。考虑到每一种工况的每一组独立数据之间最大的偏差大约增加一

15、倍，而实验的偏差几乎保持不变。这些验证研究结果与物理、数学模型非常接近，因此现在该模型可以用于模拟与该模型相似物理现象的不同运行工况。2.8风扇供风情况的解决方法结霜工况下风扇供风热交换器的模拟仿真过程与Chen等在在附录5、8中提出的情况很相似；只是，现在的气流速度是随积霜翅间的压降而变化的。 图3示出了用于测量翅间气流速度与压降的轴流风机与离心式风机的共同的典型特性曲线。图3中A点为热交换器表面结霜之前的气流速度状态点。B点为由于热交换器表面结霜气流速度降到最初的50的轴流风机状态点。对于离心风机，B点为风机气流停止不前而气流速度明显下降即将接近到零的状态点。B点和B点为除霜预期的起点。为

16、了简化模型，可以假设翅间空气压降为主导因素。气流循环的其他压降可以忽略。该模型是有气流速度与不同类型风机压降的关系再加上其他的模型来完成的。 图3还显示了该热交换器翅片表面与基表面持续结霜过程的三个不同时间点t1、t2和t3（t3t2t1）的压力流量特性的系统压降曲线。风机特性曲线与系统P-Q特性曲线之间的状态点就是风机特定时间的工作点。这交汇区给出了翅间气流速度与压降。 整个数据模型的求解方法： （1）使用最后时间的翅片表面温度分布，如，对流传热传质系数，如和，及平均气流湿度和温度，如和，作为边界条件解决翅片表面上每个格点的结霜模型。翅片表面霜层高度是随霜的其他特性而得到的。 （2）解决非翅

17、片基表面结霜模型。进而得出霜高，。 （3） 解决翅片热传导模型，其中源项从第一步得出。 （4）运用已知的计算平均霜高。 （5）用方程（11）和（17）求得系统P-Q特性关系。 （6）运用风机曲线与系统P-Q关系计算翅间气流速度的与压降。 （7）运用Gnielinski相关计算传热传质系数。 （8）运用方程（18）、（19）和（20）计算翅片与非翅基面的热流量和总热流量。 （9）运用方程（21）(24)计算平均气流温度与湿度，和。 （10）重复第一至第九步，直至结霜过程结束。3.模拟结果与讨论 模拟结果是建立在表1示出的数据为运行工况的基础之上。图4示出了轴流风机的瞬间压力损失与气流入口速度。图

18、4显示出，在5.5h的时间内，由于热交换器翅表面结霜的缘故，气流速度降低到最初的50，同时压降增加到初始的3.2倍。这些数据意味着随着时间的增长，气流速度与压降间的近似线性的变化关系。附录5、8中Chen等研究了在他们假设的恒定气流速度条件下，压降随时间趋于变快的情况。 图5中部分翅化气流区被霜完全堵塞同时总的翅化气流区无霜的情况被称作霜堵比点。图中显示出，在运行5.5小时后，多余50的翅间气流区被霜层堵塞。Chen等在附录5、8中指出，翅间气流通路的这种堵塞情况主要发生在冷基表面，这种冷基表面可能会被完全的堵塞并且随时间的增长会从基表面向翅尖发展。 图6显示了轴流风机供风的热交换器的每个翅片

19、各个侧面的瞬态热流量。在结霜的最初时刻（t30s），可以观察到热流量的缓慢增长（从3.9-4.1W翅片）。这最初增长的5的热流量（即0.2W翅片）是水蒸气在最初清洁的的翅表面结霜相变传质并直接向翅表面释放潜热的缘故。在这最初的短暂结霜时间之后，每个翅片的热流量都随时间的增长而下降。在结霜5.5h后，气流速度降低到最初的一半，热流量降到最初的11.7（从4.1-2.4W翅片 ）。热流量的下降有两方面的因素：（1）翅表面的霜层像一个保温层一样，因此它增大了翅表面的热阻、降低了热流量；（2）翅表面结霜堵塞了气流通路，导致气流速度下降进而降低了空气和翅表面的对流传热。附录5、8中Chen等指出当在气流

20、恒定和运行工况相似的条件下，翅的热流密度至多降到原来的56。图6的结果表明风机的选择为热流量下降的主要影响因素。 为了表明霜积累是翅片效率随温度的变化关系，引进更改后的新的翅片效率的概念。该翅片效率定义为， （25） 此处 结霜工况下任意时间的翅的热流量，； 翅表面面积，； 无霜下对流传热系数，； 翅间平均气流温度，； 基面温度，。 图7显示了5.5小时以后该翅效率从0.8降到0.5的时间变化。翅效率这种下降的原因与上面提到过的热流量的相似。4.灵敏性研究 该模型的的灵敏性研究是用来探讨各种参数的变化对模型在应用中的影响。 该灵敏性研究中涉及了三个参数。它们是翅间距、翅厚和风机曲线。 表3列举

21、了用于灵敏性研究的参数，而图8-12显示了以轴流风机和表3的数据为基本参数的灵敏性研究结果。图3中当运行工况的气流状态从A变化到B或B时，只有一个参数是随每个结果的基本情况和每个模拟的终止时间而变化的。缩小翅间距使得除霜循环间的热交换器运行时间缩短。同样对于气流速度下降50，基本情况下需要330分钟，而当翅间距降低22时只需要280分钟。缩小翅间距降低了每一个翅的翅效率和热流量。较小的翅间距使得霜堵塞的更快以及空气和翅表面的传热下降的也更加的快。 另一方面，翅间距的缩小将会使得单位气流面积的翅片数增多，以此来提高总的热流密度。因此，这就表明若将此模型建立在各种冷冻食品的融霜循环的现实要求的基础

22、之上，就可以进行翅间距与翅厚的优化。 把翅厚减少50每个翅的热流量就比最初降低25。对于越薄的翅片来说，翅的热流量与效率随时间下降的越慢，这就使得融霜循环间的热交换器运行时间变得更长。如图8-12所示的，在该研究中将轴流风机改为离心风机对翅片的特性并没有显著的影响。而且，图9、10中曲线4的热流量与气流速度下降到零的终端点不变。在附录5中Chen等已经证实，运行工况的进气速度、温度和湿度以及热交换器冷基面的温度都对结霜过程与热交换器的性能有显著的影响。读者可以参考这些结果进行介绍与讨论这些敏感性实验。5.结论 该文已经证明翅片表面的精确结霜模型与选定的风机的特性曲线的耦合对制冷热交换器的优化设计而言必不可少。结霜工况下，热流量、翅片效率以及压降、气流速度都随时间明显降低。风机、翅间距及翅厚的设计与选型将影响热交换器的结霜过程和融霜周期。