结垢1.1

import numpy as np
import math
import json

m =(1, 2, 3, 4)
l =(1, 3, 5, 7)
n =(1, 4, 7, 10)
HE_fouling=json.dumps({'HE_inlet_pre':m, 'HE_outlet_pre':l, 'T_volume_flow':n})


def HE_Fouling(HE_fouling1):
    '''
    热交换器结垢
    本函数的作用为求解热交换器壳程阻抗系数和结垢厚度
    Keyword arguments:
        input：
            HE_tube_num 输入换热管数目
            HE_tube_length 输入换热器换热管长度
            HE_tube_diameter 输入换热管内径
            HE_shell_diameter 输入热交换器壳体直径
            Fluid_density 输入介质流体密度
            volume_flow 管程流体体积流量
            inlet_pre 管程进口压力
            outlet_pre 管程出口压力
        mid-parameters
             Fl_res 流动阻抗系数，其值为64/雷诺数，需要查明
             M 定义中间变量
        output：
            HE_shell_impedance 壳程阻抗系数
            HE_T_scaling_thickness1 结垢厚度

    calculate process：
    首先利用管程进口压力、管程出口压力、管程流体体积流量计算壳程阻抗系数。
    然后使用介质流体密度、换热管长度及壳程阻抗系数计算中间变量M。
    最后使用 换热管数目，输入热交换器壳体直径，输入换热管内径，输入介质流体密度，管程流体体积流量计算结垢厚度。
    '''
    HE_tube_num = int(input('换热管数目\t'))
    HE_tube_length = float(input('换热器换热管长度\t'))
    HE_tube_diameter = float(input('换热管内径\t'))
    HE_shell_diameter = float(input('热交换器壳体内径\t'))
    Fluid_density = float(input('热交换器管程流体密度\t'))

    #line 48-52 将HE_fouling1转换为python可读取的字典形式，读取并将数据转换为数组格式。
    Fouling = json.loads(HE_fouling1)
    print(f'管程结垢{Fouling}')
    volume_flow = np.array(Fouling['T_volume_flow'])
    inlet_pre = np.array(Fouling['HE_inlet_pre'])
    outlet_pre = np.array(Fouling['HE_outlet_pre'])

    #line 55-61 计算壳程阻抗系数和结垢厚度
    HE_shell_impedance = ( inlet_pre-outlet_pre) / volume_flow**2
    Fl_res=1 #流动阻抗系数，其值为64/雷诺数，需要查明。管内介质在某温度下流动阻力系数。
    HE_shell_impedance = np.array([1.1,2.2,3.3,4.1])
    M= pow(((8 * Fluid_density * HE_tube_length * Fl_res) / (math.pi * math.pi * HE_shell_impedance)), 1 / 5)
    HE_T_scaling_thickness1 = (pow(
    HE_tube_num **2  * M **2 - 4 * HE_tube_num * M * HE_shell_diameter + 4 * HE_tube_num * HE_shell_diameter **2,
    0.5) - 2 * HE_tube_num * HE_tube_diameter - HE_tube_num * Fluid_density * volume_flow) / (4 * HE_tube_num)

    #将数据转换为元组形式，并存入字典。最后将字典转换为json格式。
    HE_shell_impedance = tuple(HE_shell_impedance)
    HE_T_scaling_thickness1 = tuple(HE_T_scaling_thickness1)
    Fouling_output = {'shell_impedance':HE_shell_impedance, 'caling_thickness':HE_T_scaling_thickness1}
    print(f'结垢厚度：{Fouling_output}')
    Output_fouling = json.dumps(Fouling_output)



HE_Fouling(HE_fouling)
# help(HE_Fouling)
# 上产分析
# 用效率
# 流动阻力损失