计算两股不同流量的气体，通过换热器换热后，高温气体的出口温度

在这里插入图片描述

# -*- coding: utf-8 -*-
"""
Created on Thu Nov 30 11:23:12 2023
计算两股不同流量的气体，通过换热器换热后，高温气体的出口温度
(煤烟二级，计算煤烟二级热侧出口温度)
------------------------------------------------
已知:冷测烟气量、出入口温度,热侧烟气量、热侧入口温度
求热侧出口温度
保证热出比冷进高15℃
------------------------------------------------
@author: xujifa
"""
import sympyclass Com_two_hot_out_temperature():# 已知物理参数、边界条件等#冷测入口烟气量 Nm3/h，#冷测入口温度 ℃，#冷测出口温度 ℃，#热测入口烟气量 Nm3/h，#热测入口温度 ℃def __init__(self, cold_in_fl, cold_in_tem,cold_out_tem,hot_in_fl,hot_in_tem,name):self.cold_in_fl = cold_in_fl  #冷测入口烟气量 Nm3/hself.cold_in_tem = cold_in_tem  #冷测入口温度 ℃self.cold_out_fl = cold_in_fl  #冷测出口烟气量 Nm3/hself.cold_out_tem = cold_out_tem  #冷测出口温度 ℃self.hot_in_fl = hot_in_fl  #热测入口烟气量 Nm3/hself.hot_in_tem = hot_in_tem  #热测入口温度 ℃self.hot_out_fl = cold_in_fl  #热测出口烟气量 Nm3/hself.name = str(name)self.cha = 15 #保证热出比冷进大于15℃print(f"👇👇👇以下计算为 {self.name} 的参数👇👇👇")#定义热损失，换热器效率，根据热测入口温度 ℃来定义def ex_eff(self):if self.hot_in_tem > 400:return 0.94else:return 0.965#fl 风量Nm3/h、 tem温度℃，求携带热量def com_heat(self,fl,tem):fl = fl/0.0224 #气体摩尔量 mol/htem = tem + 273.15 #温度转化为k，因为比热容的方程是根据温度k来计算的cp = 7.22e-10*pow(tem,3) - 5.937e-06*pow(tem,2) + 0.01707 *tem + 26.52 #比热容 J/mol-Kheat =  cp*fl*tem #计算携带的热量  J/h# print(f"---携带热量模块---{heat} j/h")#f{}格式化字符串return heat#换热计算,高温在前，低温在后,风量Nm3/h、 tem温度℃def ex_heat(self):cold_in_heat = self.com_heat(self.cold_in_fl,self.cold_in_tem)cold_out_heat = self.com_heat(self.cold_out_fl,self.cold_out_tem)heat = cold_out_heat - cold_in_heat #计算进出口的热量差值  J/h,1kj/h = 1/3600kwprint(f"***{self.name} @ 冷测热量：{heat} J/h",)print(f"***{self.name} @ 冷测热量：{heat/3600/1000} KW",)return heat#计算高温出口热侧所需热量，设置换热效率def com_hot_out_heat(self):hot_in_heat = self.com_heat(self.hot_in_fl,self.hot_in_tem)  #热侧进口热量print(f"+++{self.name} @ 热侧入口携带热量：{hot_in_heat} J/h")cold_ex_heat = self.ex_heat()/self.ex_eff() #冷测的换热量,换热效率hot_out_heat = (hot_in_heat - cold_ex_heat)/(self.hot_in_fl/0.0224) #热侧出口所需热量print(f"+++{self.name} @ 热侧出口所需热量：{hot_out_heat} J/h")return hot_out_heat#求解一元四次函数，tem为变量def sol_tem(self,h):tem = sympy.Symbol("tem")s = sympy.solve(7.22e-10*pow(tem,4) - 5.937e-06*pow(tem,3) + 0.01707 *tem*tem + 26.52*tem-h)print(f'↓↓↓***{self.name} @ 计算结果温度 K-待筛选***↓↓↓')print(s)print(f'↑↑↑***{self.name} @ 计算结果温度 K-待筛选***↑↑↑')return s#对sol求解的函数进行后处理def post_sol(self,data):for num in data:if "I" not in str(num):if num > 0:print(f"👉👉👉{self.name}@热侧出口温度 {round(num-273.15,2)} ℃")print(f"👆👆👆以上计算为 {self.name} 的参数👆👆👆")return round(num-273.15,2)else:print("solve-复数")return -911111111  #说明是复数def getdata(self):hot_out_tem = self.post_sol(self.sol_tem(self.com_hot_out_heat()))#直接调用函数print(f"热侧出口温度为{hot_out_tem}-经过方程求解")# if hot_out_tem - self.cold_in_tem < self.cha and hot_out_tem>0:#     hot_out_tem = self.cold_in_tem + self.cha# else:#     hot_out_tem = -99999999print("---------------------------------")print(f"{self.name}边界参数")print(f"{self.name}-冷测：\n1、流量{self.cold_in_fl}Nm3/h\n2、入口温度{self.cold_in_tem}℃\n3、出口温度{self.cold_out_tem}℃\n4、温差{self.cold_out_tem-self.cold_in_tem}℃")print(f"{self.name}-热侧：\n1、流量{self.hot_in_fl}Nm3/h\n2、入口温度{self.hot_in_tem}℃\n3、出口温度{hot_out_tem}℃\n4、温差{self.hot_in_tem-hot_out_tem}℃")print("---------------------------------")return self.cold_out_fl,self.cold_out_tem,self.hot_out_fl,hot_out_temif __name__ == "__main__":Com_two_hot_out_temperature(122000, 240, 280, 15000,200,"二级换热器").getdata()