一种水冷蒸发器（表冷器设计方法）

水冷蒸发器（表冷器）校核计算：蒸发器为铜管铝翅换热器，制冷工质为水，给空气降温。
已知蒸发器进口空气状态参数，假定一个空气出口参数
已知蒸发器水的进口温度，假定水的出口参数，通过迭代比较算出最终的制冷量和出口状态参数。
（1）蒸发器结构参数：采用强制对流的直接换热蒸发器，连续整体铝套片，铜管为 $\phi$ =10×0.5mm正三角形排列，管间距s1=25mm，排间距s2=21.65mm，铝片厚度 $\delta$=0.15mm，片距sf=1.8mm，取空气流动方向管排数为2排，迎面风速为2m/s

（2）空气侧进口参数：空气进口温度为：35℃，相对湿度为：80%（湿工况），空气流量为0.209kg/s
（3）水侧进口参数：水进口温度为15℃，水流量为0.3kg/s
空气和水 的换热关系如下图：

计算流程如下图：

通过python进行计算： 

import math
import matplotlib.pyplot as plt

def Heat_Exchangers():  #定义换热器参数
    s1 = 0.025 #管间距
    s2 = 0.02165 #排间距
    d0 = 0.01 # 管道直径，单位：m
    l = 10 #管长，m
    tt = 0.0005 # 管道壁厚，单位：m
    deltaf = 0.00015  # 翅片厚度，m
    k_cu = 401  # 铜管的导热系数
    sf = 0.0018 #片距
    n = 2 #空气流动方向管排数
    v = 2 #迎面风速，m/s
    af = 2*(s1*s2-3.14*d0*d0/4)/sf #每米管长翅片表面积， m2
    ab = 3.14*d0*(sf-deltaf)/sf  # 每米管长翅片间基管外表面积，m2
    aof=af+ab # 每米管长外总表面积，m2
    a_a = aof*l #风侧换热面积
    ai = 3.14*(d0-tt)*1 #每米管长内表面积，m2
    a_w = ai*l #水侧换热面积
    beta = aof/ai #肋化系数
    alpha = aof/s1 #肋通系数（每排肋管外表面积与迎风面积比）
    epsilon = (s1-d0)*(sf-deltaf)/s1/sf #净面比（最窄流通截面与迎风面积比）

    Heat_Exchangers = [v,s1,sf,deltaf,d0,aof,n,a_a,a_w,s2,l,ai]
    return Heat_Exchangers

def humid_air(t,phi):  #干球温度、相对湿度
    T = t + 273.15 #开尔文温度
    B = 101325 #大气压Pa
    Pqb = 100*math.pow(10,((-7.90298)*(373.16/T-1)+5.02808*math.log10(373.16/(T))+(-1.3816e-7)*(math.pow(10,(11.344*(1-T/373.16)))-1)+0.0081328*(math.pow(10,((-3.49149)*(373.16/T)-1))-1)+math.log10(1013.246)))#饱和水蒸气分压
    Pq = phi*Pqb#水蒸气分压
    rho = 0.003484*B/T-0.00134*Pq/T #空气密度kg/m3
    Cp_a = 1005  # 水比热容，单位：J/kg/K

    Pr_a = 0.7026  # 计算空气侧普朗特数
    vis_a = 1.59e-5  # 空气运动粘度，单位：m2/s
    v = 1/rho #比容

    d = 662*Pq/(B-Pq)
    h1 = 1.005*t+0.001*d*(2501+1.86*t)#(空气焓值：kj/kg干空气)
    td = -35.957-1.8726*math.log(Pq)+1.1689*math.log(Pq)*math.log(Pq) #露点温度

    humid_air = [Pqb,Pq,rho,Cp_a,Pr_a,vis_a,v,d,h1,td]
    return humid_air

def Air_Parameters(t_ain,phi_ain,t_aout,phi_aout):  #定义空气侧物性参数和边界条件
    global Heat_Exchanger,humid_air
    Heat_Exchanger = Heat_Exchangers()
    v,s1,sf,deltaf,d0 =Heat_Exchanger[0],Heat_Exchanger[1],Heat_Exchanger[2],Heat_Exchanger[3],Heat_Exchanger[4]
    t_ain = t_ain  # 空气进口温度，单位：°C
    phi_ain = phi_ain #进口空气想读湿度
    t_aout = t_aout # 空气出口温度，单位：°C
    t_f = (t_ain+t_aout)/2 # 空气平均温度，单位：°C
    Cp_a = humid_air(t_f,phi_ain)[3]  # 水比热容，单位：J/kg/K
    rho_a = humid_air(t_f,phi_ain)[2] # 空气密度
    Pr_a = humid_air(t_f,phi_ain)[4] # 计算空气侧普朗特数
    vis_a = humid_air(t_f,phi_ain)[5]  # 空气运动粘度，单位：m2/s
    v_max = v*s1*sf/(s1-d0)/(sf-deltaf) #最窄空气流速
    Re = v_max*d0/vis_a #雷诺数
    aof,n = Heat_Exchanger[5],Heat_Exchanger[6]
    alpha4 = 0.0014+0.2618 * math.pow(Re,(-0.4)) * math.pow((aof/0.0314),(-0.15))#空气干侧表面传热系数，利用McQuistion通用关联式，4排基管
    alpha2 = alpha4*0.992*math.pow((2.24*math.pow(Re,(-0.092))*math.pow(n/4,(-0.031))),0.607*(n-4))#两排管平均干表面传热系数
    alpha0 = alpha2*rho_a*v_max*Cp_a/(math.pow(Pr_a,0.6667)) #实际干表面传热系数
    Air_Parameters = [t_f,Cp_a,rho_a,Pr_a,vis_a,v_max,Re,alpha0]
    return Air_Parameters

    #确定空气在蒸发器中的状态变化

def Water_Parameters(t_win,t_wout,G_w):  #定义水侧物性参数和边界条件
    global humid_air
    di = Heat_Exchangers()[4]-0.001
    t_win = t_win #进水温度15℃
    t_wout = t_wout  # 出水温度15℃
    t_wf = (t_win + t_wout)/2 # 水侧平均温度
    Cp_w = 4200  # 水比热容，单位：J/kg/K
    G = G_w  # 水入口质量流量，单位：kg/s
    rho_w = 995.7  # 水流体密度
    k_w = 0.618  # 水导热系数，单位：W/m/K
    vis_w = 0.656e-6  # 水运动粘度，单位：m2/s
    u = G / (3.14*di*di/4)  # 计算水入口速度
    a = k_w / rho_w / Cp_w  # 计算换换热器热扩散率
    Re_w = u * di / vis_w  # 计算水侧雷诺数
    Pr_w = vis_w / a  # 计算水侧普朗特数
    h_w = 0.0023* math.pow(Pr_w,0.333)*math.pow(Re_w, 0.8) *0.618/di  # 水侧计算换热系数
    Water_Parameters = [Cp_w,h_w]
    return Water_Parameters

def Calculation(t_ain,phi_ain,t_aout,phi_aout,G_a,t_win,t_wout,G_w):
    global  humid_air,Air_Parameters,Water_Parameters

    #换热器参数
    s1 = Heat_Exchangers()[1]  # 管间距
    s2 = Heat_Exchangers()[9]  # 排间距
    d0 = Heat_Exchangers()[4]  # 管道直径，单位：m
    tt = 0.0005  # 管道壁厚，单位：m
    deltaf = Heat_Exchangers()[3]  # 翅片厚度，m
    k_al = 207 #铝的导热系数
    sf = Heat_Exchangers()[2]  # 片距
    n = Heat_Exchangers()[6]  # 空气流动方向管排数
    v = Heat_Exchangers()[0]  # 迎面风速，m/s
    af = 2 * (s1 * s2 - 3.14 * d0 * d0 / 4) / sf  # 每米管长翅片表面积， m2
    ab = 3.14 * d0 * (sf - deltaf) / sf  # 每米管长翅片间基管外表面积，m2
    aof = Heat_Exchangers()[5]  # 每米管长外总表面积，m2
    A_a = Heat_Exchangers()[7]  # 风侧换热面积
    A_w = Heat_Exchangers()[8]  # 水侧换热面积

    #水侧参数
    t_win = t_win  # 进水温度15℃
    t_wout = t_wout  # 出水温度15℃

    Cp_w = Water_Parameters(t_win,t_wout,G_w)[0]  # 水比热容，单位：J/kg/K
    G = G_w  # 水入口质量流量，单位：kg/s

#空气侧参数
    t_ain = t_ain  # 空气进口温度，单位：°C
    phi_ain = phi_ain  # 进口空气想读湿度
    t_aout = t_aout  # 空气出口温度，单位：°C
    phi_aout = phi_aout  # 进口空气想读湿度
    G_a = G_a #空气质量流量，kg/s


    cur_eps = 1
    x = []
    y = []
    i = 1

    if t_win >= humid_air(t_ain,phi_ain)[9]: #干工况
        while (cur_eps > 0.01 and i < 100):
            xi=1 #析湿系数
            t_f = (t_ain + t_aout) / 2  # 空气平均温度，单位：°C
            Cp_a = humid_air(t_f, phi_ain)[3]  # 空气比热容，单位：J/kg/K
            alpha0 = Air_Parameters(t_ain, phi_ain, t_aout, phi_aout)[7]
            m = math.sqrt(2 *alpha0*xi*1/k_al/deltaf)
            rho1 = 1.27*s1/d0*math.sqrt(0.7)
            hh = d0/2*(rho1-1)*(1+0.35*math.log(rho1)) #肋片折合高度
            etaf = math.tanh(m*hh)/(m*hh) #翅片效率
            alphaj = xi*alpha0*((etaf*af+ab)/(af+ab)) #空气侧当量表面传热系数
            h_w =Water_Parameters(t_win,t_wout,G_w)[1] # 水侧计算换热系数
            K = 1 / ((A_a / A_w) * 1 / h_w + 1 / alphaj+tt/k_al)  # 计算当量传热系数,237为铝的导热系数
            thta = abs(((t_ain - t_win) - (t_aout - t_wout)) / (math.log((t_ain - t_win) / (t_aout - t_wout))))  # 对数温差
            Heat_Flow = K * A_a * thta
            t_wnext = t_win + Heat_Flow / (Cp_w * G)
            t_anext = t_ain - Heat_Flow / (Cp_a * G_a)
            t_wout = (t_wout + t_wnext) / 2
            t_aout = (t_wout + t_anext) / 2
            Heat_Flow = K * A_a * thta
            Heat_Flow1 = Cp_w * G * math.fabs(t_wout - t_win)
           # Heat_Flow2 = Cp_a * G_a * math.fabs(t_aout - t_ain)
            cur_eps = math.fabs(Heat_Flow1 - Heat_Flow)
            i = i + 1
            x.append(i)
            y.append(cur_eps)
            plt.plot(x, y, color='b')  # 线条颜色 bluelinewidth=1,marker='o',  # 折线图中添加圆点markersize=3,  # 点的大小
            plt.xlabel('迭代次数')
            plt.ylabel('误差值')
            # 显示图形
            plt.pause(0.1)
            plt.show()
        print('干工况','\n''出水温度：',t_wout,'\n''出风温度：',t_aout,'\n''换热量：',Heat_Flow,'\n''迭代次数：',i)


    else:   #湿工况
        while (cur_eps > 0.1 and i < 60):
            t_f = (t_ain + t_aout) / 2  # 空气平均温度，单位：°C
            Cp_a = humid_air(t_f, phi_ain)[3]  # 水比热容，单位：J/kg/K
            rho_a = humid_air(t_f, phi_ain)[2]  # 空气密度
            Pr_a = humid_air(t_f, phi_ain)[4]  # 计算空气侧普朗特数
            vis_a = humid_air(t_f, phi_ain)[5]  # 空气运动粘度，单位：m2/s
            v_max = v * s1 * sf / (s1 - d0) / (sf - deltaf)  # 最窄空气流速
            Re = v_max * d0 / vis_a  # 雷诺数
            #析湿系数
            d1 = humid_air(t_ain, phi_ain)[7]
            d2 = humid_air(t_aout, phi_aout)[7]
            dm = (d1+d2)/2
            # 对流换热系数计算
            # xi =1+ 2.501*(d1-dm)/(1.005*(t_ain-t_aout))
            # alpha4 = 0.0014 + 0.2618 * math.pow(Re, (-0.4)) * math.pow((aof / 0.0314),(-0.15))  # 空气干侧表面传热系数，利用McQuistion通用关联式，4排基管
            # alpha2 = alpha4 * 0.992 * math.pow((2.24 * math.pow(Re, (-0.092)) * math.pow(n / 4, (-0.031))),0.607 * (n - 4))  # 两排管平均干表面传热系数
            # alpha0 = alpha2 * rho_a * v_max * Cp_a / (math.pow(Pr_a, 0.6667))  # 实际干表面传热系数
            # m = math.sqrt(2 * alpha0 * xi * 1 / k_al / deltaf)
            # rho1 = 1.27 * s1 / d0 * math.sqrt(0.7)
            # hh = d0 / 2 * (rho1 - 1) * (1 + 0.35 * math.log(rho1))  # 肋片折合高度
            # etaf = math.tanh(m * hh) / (m * hh)  # 翅片效率
            # alphaj = xi * alpha0 * ((etaf * af + ab) / (af + ab))  # 空气侧当量表面传热系数
            #
            #
            # h_w = Water_Parameters(t_win,t_wout,G_w)[1]  # 水侧计算换热系数
            #
            # K = 1 / ((A_a / A_w) * 1 / h_w + 1 / alphaj + tt / k_al)  # 计算当量传热系数,237为铝的导热系数
            # thta = abs(((t_ain - t_win) - (t_aout - t_wout)) / (math.log((t_ain - t_win) / (t_aout - t_wout))))  # 对数温差
            # Heat_Flow = K * A_a * thta
            Heat_Flow2 = (Cp_a * G_a * (t_ain - t_aout)) + (d1 - d2) * 250.1 * G_a
            Heat_Flow1 = Cp_w * G * math.fabs(t_wout - t_win)


            t_wnext = t_win + Heat_Flow2 / (Cp_w * G)
            t_anext = t_ain - ((Heat_Flow1-250.1*(d1-d2)*G_a) / (Cp_a * G_a))

            #d2_next = (d1 - (Heat_Flow1 - (Cp_a * G_a * (t_ain - t_anext))) / 250.1 / G_a)

            d2_next =(humid_air(t_anext, phi_aout)[8]-1.005*t_anext)/(2501+1.86*t_anext)/0.001
            print(250.1*(d1-dm)*G_a)
            t_wout = (t_wout + t_wnext) / 2
            t_aout = (t_wout + t_anext) / 2
            d2 = (d2_next + d2) / 2
            phi_aout = d2*101325/(662+d2)/(humid_air(t_anext, phi_aout)[0])
            #Heat_Flow = K * A_a * thta
            Heat_Flow2 = (Cp_a * G_a * (t_ain - t_aout)) + (d1 - d2) * 250.1 * G_a
            Heat_Flow1 = Cp_w * G * math.fabs(t_wout - t_win)
            w = (d1-d2)* G_a #除湿量,单位g

            cur_eps = math.fabs(Heat_Flow1 - Heat_Flow2)

            i = i + 1
            x.append(i)
            y.append(cur_eps)
            plt.plot(x, y, color='b')  # 线条颜色 bluelinewidth=1,marker='o',  # 折线图中添加圆点markersize=3,  # 点的大小
            plt.xlabel('迭代次数')
            plt.ylabel('误差值')
            # 显示图形
            plt.pause(0.1)
            plt.show()
        print('湿工况', '\n''出水温度：', t_wout, '\n''出风温度：', t_aout, '\n''换热量：', Heat_Flow1,'W','\n''除湿量：',w,'g' ,'\n''迭代次数：', i)

#输入边界条件：
#空气侧参数
t_ain = 35  # 空气进口温度，单位：°C
phi_ain = 0.8  # 进口空气想读湿度
t_aout = 25  # 空气出口温度，单位：°C
phi_aout = 0.8  # 进口空气想读湿度
G_a = 0.209 #空气质量流量，kg/s
# 水侧参数
t_win = 15  # 进水温度15℃
t_wout = 16  # 出水温度15℃
G_w = 0.3  # 水入口质量流量，单位：kg/s

Calculation(t_ain,phi_ain,t_aout,phi_aout,G_a,t_win,t_wout,G_w)

参考资料：《小型制冷装置设计指导 》《制冷空调系统仿真原理与技术》《制冷技术与装置设计》