结垢1.1

结垢1.1

import numpy as np
import math
import jsonm =(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-parametersFl_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**2Fl_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)
# 上产分析
# 用效率
# 流动阻力损失