（译文）换热器英文参考文献

应用计算数值的方法来研究流体的粘度变化

对板式换热器性能的影响

M.A. Mehrabian and M. Khoramabadi

Department of Mechanical Engineering, Shahid Bahonar University of Kerman,

Kerman, Iran

摘 要

目的--本文的目的是在逆流和稳态条件下，通过数值计算，研究流体粘度的

变化对板式换热器热特性的影响。

设计/工艺/方法--实现这篇文章目的的方法，源于由4部分组成的热量交

换板中间通道中冷热流体的一维能量平衡方程。有限差分法已经用于计算温度分布及换热器的热性能。在侧边通道中，水作为将被冷却的热流体，然而在中央通道中，大量随温度变化同时粘度随之变化剧烈的流体作为将要被加热的冷流体。 发现—这个程序的运行实现了工作流体的结合，例如水与水，水与异辛烷，水与苯，水与甘油和水与汽油等。对于以上所有工作流体的结合，两种流体的温度分布已经沿流动通道划分。总传热系数可以通过冷流体和热流体的温度来绘制。研究发现，若总传热系数呈线性变化，在温度变化范围内既不是冷流体和热流体的温度。当粘度已受温度影响或者冷流体的性质改变时，换热器的影响效果并不是很显著。

创意/价值--对于由2块板为边界的温度控制体来说，本文包含一个可以得

到能量平衡方程数值解的新方法。通过对数值计算结果与实验结果进行比较，验证了这种数值计算方法。

1

关键词:热交换器、热传递、数值分析、有限差分法 研究类型:研究性论文。

2

术 语

A:板传热面积，m2

b:板间距，m C:等式常数

C:热容，W/?C

Cp:定压比热容，J/kg??C

De:当量直径，m

h:对流传热系数，W/m2??Cj:指定轴截面

k:板传导率，W/m??C

L:板长度，m m:粘度修正系数

?m:质量流量，kg/s

n:NuP?1r3与Re之间的斜率

NTU: 传热单元数 Nu: 努塞尔数 Pr: 普朗特数 Q: 传热速率, W Re: 雷诺数

3

r: 方程指数 (8) t: 时间, s T: 温度, ℃ u: 流速, m/s

U：总传热系数，W/m2??C ?U:平均传热系数，W/m2??C

V:通道体积，m3

w: 流动宽度, m x: 横向坐标 y: 轴向坐标

m:流体动粘度系数，kg/m?sr:流体密度，kg/m3

l: 换热器有效性 d: 板厚度, m f: 板投影面积的比值

4

下 标

c : 冷流体 Cv: 控制体 h : 热流体 m : 平均值 min:最小值 w : 板壁

5

数值分析

数值分析法中使换热器分成一些轴向的部分。一个典型的轴向部分都有一个表面积

。对于这个增加的表面积，冷热流体的温度分别是

和

，我们

可以假设总传热系数可以作为这些温度的函数而表示出来。这样：

等式2可以应用在轴截面上，表示为：

等式（3）和（4）也可以运用在换热器相邻通道的两个轴截面上， 可写为：

上述方程的解的获得是当空间导数存在偏差时。以viscosities(Yaws, 2003)为依据的温度数据表被编入计算机程序中，并且这个程序可以表示出每个轴截面上，流体流动时的温度下的黏度。线性插值的操作就开始进行，此时温度数值与表值不一致。像密度、热导率等一些其他的流体性质与温度无关。每种流体的这些特性的数值以平均流体温度来指定，并且作为输入数据。冷热流体的入口温度作为数值分析的边界条件。

板式换热器通道中的流体无量纲传热系数可看成是与热传递相关的一种类型

11

(Rao et al., 2002)：

Shah and Focke (1988)进行了实验研究板式换热器热传递和压降特性。他们注意到，常数C取决于换热板的类型和换热器的几何形状，而常数n取决于流体的流态。

Edwards et al. (1974)研究证明得出，在雷诺数大约小于10时，实验数据是以

标绘的，而不是 APV Junior Paraflow板落在表明典型传热关系

的坡度线1/3处的Re值：

在雷诺数较高时（Re大于10），坡度约为0.7，这样会得出过度条件和湍流条件：

这种关系也可能成为相互距离b的两个平行板之间的湍流类型。可假设为，由于板的褶皱，取决于当量直径

的雷诺数会影响热传递的增加。对于牛顿

流体，Edwards et al得到的结果，表明APV Junior Paraflow板被作为流动通道，并可以推广到任何板的类型，提供常数C来作为修正值。

Mehrabian(1996年)进行了广泛的研究，从实验和理论观点探究流体动力学和板式换热器的热性能。在湍流（Re大于10）条件下，他对APV SR3板提出了以下的关系：

等式（9）同时适用于板式换热器冷热流体通道，传热系数分别为截面j的总传热系数是：

12

和h，轴

正如Edwards et al所提出的，在雷诺数低至10的情况，板式换热器流动通道中的流体也可能形成湍流。因此，湍流假设分析是合理的，并且等式（9）适用于水为介质的热流体，同时甘油、苯和辛烷可作为冷流体。这些冷流一定要选择其粘度随温度而变化。

13

结果和讨论

为了得到独立网格数据结果，程序运行时将轴向分为几个不同的部分。将确定的网格点的数值结果与相应的(Haseler et al., 1992)的实验结果进行对比，然后记录两者之间的差值。值得注意的是，增加网格点的数量会减少差值。然而，当轴向部分的数目是17的时候，此时可产生最小差值，当超过17时，差值减小并不明显。通过实验结果表明，轴向部分的数目为17。而3，6,9,15这些数字都是3的倍数。这是在(Haseler et al., 1992).中对比相应点处实验结果的目的。

应该被提及的是，在获得数值解之初，两种流体的出口温度并不清楚，可通过两者的进口平均温度来估算。当得到出口温度时，每种流体的性质可以由其本身的平均温度来估算。(Khoramabadi, 2004).

换热板尺寸和和流体流动条件的数值分析列于表1.这些数值都是(Haseler et al., 1992)中用于试验温度测量的所需尺寸和条件。水的典型温度实验的温度列于表II，以当量直径

的热流体和冷流体的雷诺数分别是879和304。

换热器中间通道的局部温度的数据结果然后与V型区域的实验温度相比较，列于表II并得到数值上的一致。这种误差在

范围内，这样才表明

这种数值法的准确性。由于数值计算程序已经被验证，进口冷热流体之间的变大的温度差会使流体粘度显著影响温度分布、总传热系数和换热器的传热性能。冷流体和热流体的进口温度分别是

这个程序适用于水-水、水-苯、水-

异辛烷四种类型的流体。冷热流体的沿换热器温度分布如图2和3。水-苯与水-苯的温度分布相似，因此，并没有在图2和3中。对于每种类型的工作流体，恒定和变化的流体粘度的温度分布都有所体现。可以注意到，水之间的温度分布的

14

误差大于甘油和异辛烷之间的误差。然而，随着温度改变，水的粘度的变化比甘油和异辛烷较为强烈。这种现象的原因是甘油和异辛烷的对流传热系数小于水之间的对流传热系数，因此，它们控制了总传热系数。

表1 板的尺寸和流动细节

表2 数值比较和实验结果

15

图2 冷流体沿流动通道的温度分布

16

图3 热流体沿流动通道的温度分布

图4 总传热系数与冷流体温度之间关系

17

图5 总传热系数与热流体温度之间的关系

图4和图5分别描绘出了总传热系数与冷热流体之间的变化关系。然而，当假设流体粘度不变时，总传热系数也是不变的。这样使得粘度的变化是控制总传热系数的一个重要的因素。U与T之间的线性关系和液体粘度与温度的关系都在图6中体现。表2中，在在温度范围内粘度显示出了线性变化关系。然而，这种现象并不违背大多数流体以指数函数与温度之间作为的标准。热流体和冷流体传热系数随温度线性变化分别在图7和图8中。

换热器恒定粘度和多变粘度的

关系图在图9中。此图是表示水与异

辛烷之间的关系。然而水-水、水-苯和水-甘油都有相似的图形。换热器的性能并不因温度控制下的粘度影响。

当粘度只受温度限制时，对于不同的工作流体图10描绘了

18

图形。

只要NTU很小（NTU<1），换热器的性能并不取决于流体的性质。但是若NTU增大（NTU>1），换热器的性能就会因不同的流体而改变。与图9和图10有关的性能计算基于NTU和U的定义：

整合起来可以为Tc or Th.

图6 粘度随流体温度的变化

19

图7热端传热系数与热流体温度之间的关系

20

图8 冷端传热系数与冷流体温度的关系

21

图9 可变的U对水-异辛烷中ε-NTU的影响

22

图10 流体种类对可变U的ε-NTU的影响

23

总 结

对于受温度影响的板式换热器流动通道的流体，有限差分数值计算方法已经可以预测轴向温度的变化，像水-水、水-异辛烷、水-苯、水-甘油、水-汽油这样的工作流体的组合已经通过检验。在以上工作流体组合中，换热器侧边通道的水作为将被冷却的热流体。然而，随温度而粘度变化强烈另一种流体，在中间通道中作为将被加热的冷流体。经数值分析的换热器是由四个APV SR3标准板而形成三个通道。这种换热器的几何尺寸和流体流动条件与以前的实验研究结果相符。从本文已提出的数值分析方法获得的预测温度和在相同种类的板和形似条件情况下所作实验的得出的结果相比较，误差值在-2.18和+1.88℃之间。这样，符合相应的条件，证明了数据模型的正确性。经过确认后，程序将以不同的操作条件运行。

在先前(Mehrabian》,2003)的第一作者对板式换热器总传热系数的变化的影响已经研究过。在这项调查研究中，总传热系数四种不同的变化模式被提出。其中一种是U和T的直线关系。在采用研究的温度范围内，目前获得的数据结果看出，U和Tc、U和Th也都呈直线关系。通过两种不同方法所获得的结果一致性对于设计工程师和换热器生产商是非常关注的。这就意味着，在 U和Tc或Th之间随粘度呈指数变化而与温度呈线性关系的流体类别存在比较缓和的操作温度下。

分析中所获得的结果表明，粘度的变化对换热器的换热效率只是有轻微的影响。这可以得出这样的结论，若

或

修正系数被纳入等式（8）

中，可以影响粘度的变化。基于现在的计算机程序数值预算可以运用估算指数m。

24

结 论

板式换热器热量分析中的平均粘度与局部粘度的影响，通过计算数值的方法进行研究。对于粘度随温度变化的流体，恒定粘度和可变粘度的流体温度分布之间的差别要比粘度受温度影响较弱的流体要小。，对于粘度可变化的流体，当温度变化较小且缓和时，总传热系数随冷流体或热流体温度呈线性变化。若应用粘度受温度控制的流体，则对换热器的热性能影响并不显著。当NTU较小时，流体的性质并不改变换热器的性能。

25

参 考

Colburn, A.P. (1933), “Mean temperature difference and heat transfer coefficient in liquid heat

exchangers”, Industrial and Engineering Chemistry, Vol. 29 No. 8, pp. 873-6.

Dunkley, W.L. et al. (1961), “Variation in heat treatment of individual product particles in plate

heat exchangers”, Food Technology, Vol. 15, pp. 400-3.

Edwards, M.F. et al. (1974), “Heat transfer and pressure drop characteristics of a plate heat

exchanger using Newtonian and non-Newtonian fluids”, The Chemical Engineer, May, pp. 286-8.

Foote, M.R. (1967), “Effective mean temperature differences in multi-pass heat exchangers”,

Report No. 303, National Engineering Laboratory, pp. 1-22.

Haseler, L.E. et al. (1992), “Flow distribution in a plate and frame heat exchanger”, Paper 42,

paper presented at 3rd UK National Heat Transfer Conference, Birmingham,Vol. 1, pp. 361-7.

Kern, D.Q. (1950), Process Heat Transfer, McGraw-Hill, New York, NY.

Khoramabadi, M. (2004), “Influence of variable viscosity and axial heat conduction on the

thermal performance of plate heat exchangers”, MSc thesis, Department of Mechanical Engineering , Shahid Bahonar University of Kerman, Kerman.

Mazza , V. (1984), “Plate heat exchangers for HVAC applications”, Heating, Piping and Air

Conditioning, Vol. 56, pp. 59-63.

Mehrabian , M.A. (1996), “Experimental, theoretical and computational modeling of flow in

corrugated channels to investigate thermal and hydrodynamic characteristics of plate heat

26

exchangers”, PhD thesis, Department of Mechanical Engineering, University of Bristol, Bristol.

Mehrabian, M.A. (2003), “Longitudinal temperature changes in plate heat exchangers:

analytical-numerical approach”, Iranian Journal of Science and Technology, Vol. 27, B3, pp. 583-98.

Rao, B.P. et al. (2002), “Effect of flow distribution to the channels on the thermal performance of a

plate heat exchanger”, Chemical Engineering and Processing, Vol. 41, pp. 49-58. Reppich, M. (1999), “Use of high performance plate heat exchangers in chemical and process

industries”, International Journal of Thermal Sciences, Vol. 38, pp. 999-1008.

Schlunder, E.U. (1989), Heat Exchanger Design Handbook, Hemisphere Publishing Corporation,

New York, NY.

Shah, R.K. and Focke, W.W. (1988), “Plate heat exchangers and their design theory”,

in Shah, R.K., Subbarao, E.C. and Mashelkar, R.A. (Eds), Heat Transfer Equipment Design, Hemisphere Publishing, Washington, DC, pp. 227-54.

Yaws, C.L. (2003), Yaws’ Handbook of Thermodynamics and Physical Properties of Chemical

Compounds, Knovel, Norwich, NY.

27

通 讯

M.A. Mehrabian can be contacted at: MA_Mehrabian@yahoo.com

28

本文来源：https://www.bwwdw.com/article/wqrx.html