CN112214893B - 一种换热器用单侧传热系数剥离方法 - Google Patents

一种换热器用单侧传热系数剥离方法 Download PDF

Info

Publication number
CN112214893B
CN112214893B CN202011089182.6A CN202011089182A CN112214893B CN 112214893 B CN112214893 B CN 112214893B CN 202011089182 A CN202011089182 A CN 202011089182A CN 112214893 B CN112214893 B CN 112214893B
Authority
CN
China
Prior art keywords
coefficient
heat transfer
cold
hot
heat exchanger
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011089182.6A
Other languages
English (en)
Other versions
CN112214893A (zh
Inventor
马石磊
杨志刚
牛忠华
杨洁
董昊然
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Electric Power Generation Equipment Co Ltd
Original Assignee
Shanghai Electric Power Generation Equipment Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Electric Power Generation Equipment Co Ltd filed Critical Shanghai Electric Power Generation Equipment Co Ltd
Priority to CN202011089182.6A priority Critical patent/CN112214893B/zh
Publication of CN112214893A publication Critical patent/CN112214893A/zh
Application granted granted Critical
Publication of CN112214893B publication Critical patent/CN112214893B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Heat-Exchange Devices With Radiators And Conduit Assemblies (AREA)
  • Heat-Pump Type And Storage Water Heaters (AREA)

Abstract

本发明提供了一种换热器用单侧传热系数剥离方法,所述换热器的冷热侧流道结构相同,冷热两侧的传热准则关联式中的各个系数的相应关系均未知,换热过程中,所述换热器的冷热两侧普朗特数Pr变化范围不大。本发明可以应用于冷热侧流道结构相同,冷热两侧的传热准则关联式中的各个系数的相应关系均未知,换热过程中冷热两侧Pr变化范围不大,且不方便单独研究单侧的传热系数的换热器的换热性能研究,拓宽换热器传热性能研究的实施途径。

Description

一种换热器用单侧传热系数剥离方法
技术领域
本发明涉及一种换热器用单侧传热系数剥离方法,用于剥离冷热侧介质流动和传热形式相似的单侧传热系数,属于换热器传热研究领域。
背景技术
目前针对换热器单侧传热系数的剥离方法主要有两种,分别是威尔逊图解法和修正威尔逊图解法。
参考史美中、王中铮编著的《热交换器原理与设计》,应用威尔逊图解法研究管壳式换热器的传热关系式时,当管内介质处于旺盛湍流时,管内传热系数与管内流速的0.8次方成正比,即:
Figure BDA0002721411780000011
式(1)中,αi为管内传热系数,ci为传热公式中待定系数,wi为管内介质流速。
总传热系数可以表示成:
Figure BDA0002721411780000012
式(2)中,KO为换热器总传热系数,αO为管外传热系数,rw为换热管金属热阻,rs为污垢热阻,FO为以换热管外表面为基准的换热面积,Fi为以换热管内表面为基准的换热面积。
试验过程需保证式(2)右边前三项不变,在不同管内流速wi下分别测出相应的KO,则式(2)变为:
Figure BDA0002721411780000013
式(3)就相当于一个直线方程y=ax+b,表示在y-x的直角坐标系中,为一条截距为a、斜率为b的直线,则有:
Figure BDA0002721411780000021
截距a代表定数
Figure BDA0002721411780000022
系数
Figure BDA0002721411780000023
进而得到换热器的管内传热计算公式。
威尔逊图解法适用条件如下:(1)所需研究测定的一侧传热系数与试验变量如流速的方次关系已知;(2)同一组试验中必须保证另一侧流体传热情况基本不变;(3)在同一组试验中污垢热阻基本不变。如果两侧传热基本规律均不确定,且试验过程中不能保证一侧流体传热情况恒定,则威尔逊图解法将不再适用。
因此出现了修正威尔逊图解法,具体使用过程如下。
两侧流体的传热准则关联式分别为:
Figure BDA0002721411780000024
Figure BDA0002721411780000025
式(5)、(6)中,Nu1、Nu2表示冷热两侧的努塞尔数,
Figure BDA0002721411780000026
c1、c2表示试验待定系数,Re1、Re2表示冷热侧的雷诺数,
Figure BDA0002721411780000027
Figure BDA0002721411780000028
ρ1和ρ2为冷热侧介质的密度,Pr1、Pr2表示冷热侧的普朗特数,
Figure BDA0002721411780000029
μ1、μ2表示冷热侧介质在定性温度下的动力粘度,cp1、cp2表示冷热侧介质在定性温度下的比定压热容,λ1、λ2表示冷热侧介质在定性温度下的导热系数,,μw1、μw2表示表示冷热侧介质在壁温下的动力粘度。
将式(5)、(6)改写为:
Figure BDA00027214117800000210
Figure BDA0002721411780000031
式(7)、(8)中,
Figure BDA0002721411780000032
采用平均面积计算总传热系数K如下:
Figure BDA0002721411780000033
以下标“i”表示试验点的序号,可得:
Figure BDA0002721411780000034
可以改写为:
Figure BDA0002721411780000035
式(11)中,Ki表示第i次试验的总传热系数,Re1,i、Re2,i表示第i次试验中冷热侧的雷诺数,B1,i、B2,i表示第i次试验的B1和B2值,
Figure BDA0002721411780000036
表示第i次试验冷热侧介质在壁温下的动力粘度。
式(11)也相当于一个直线方程y=a+bx,截距
Figure BDA0002721411780000037
及斜率
Figure BDA0002721411780000038
可通过线性回归求得。式(11)中的每一个试验点的值相应为:
Figure BDA0002721411780000039
Figure BDA00027214117800000310
运用该方法时,由于包含了三个未知数c1、c2和m2,所以必须选择其中某一个数,假设它的初值,通过试算来求解。其步骤为:
(1)假设c1的初始值为c10
(2)确定壁温tw1,i及tw2,i,壁温可以通过牛顿迭代法来确定。
(3)由式(9)求出α2,再利用式(8),通过线性回归求取m2
对式(8)两边取对数,得:
Figure BDA0002721411780000041
式(14)中,α2,i表示第i次试验中一侧的传热系数。
式(14)中下标i相当于某一个试验点,式(14)可以改写为:
Figure BDA0002721411780000042
此即相当于一直线方程y=a'+b'x,式中a'=lgc2,为了与下面由式(11)所得的c2相比较,令在此所得c2为c20,即c20=c2=10a',m2=b'。
(4)求c1
因m2已经得出,由式(11)线性回归得:
Figure BDA0002721411780000043
(5)比较c1与c10,c2与c20是否满足
|c1-c10|<ξ1,|c2-c20|<ξ2 (16)
式(16)中,ξ1、ξ2表示预先规定的所允许的循环收敛数值。
可以看出,修正威尔逊图解法工作量较大,一般需要编写相应的计算机程序来实现。另外还需注意如果假定的一侧系数与指数均未知,则修正威尔逊图解法也难以求解。
发明内容
本发明要解决的技术问题是:威尔逊图解法和修正威尔逊图解法存在适用范围限制。
为了解决上述技术问题,本发明的技术方案是提供了一种换热器用单侧传热系数剥离方法,所述换热器的冷热侧流道结构相同,冷热两侧的传热准则关联式中的各个系数的相应关系均未知,换热过程中,所述换热器的冷热两侧普朗特数Pr变化范围不大,其特征在于,包括以下步骤:
步骤1、设换热器的冷热侧传热公式相同,均如下式(17)所示:
Nu=c0 Rec1Pr1/3 (17)
式(17)中,Nu表示传热系数,c0、c1表示待定系数,Re表示雷诺数;
步骤2、假定待定系数c1为常数c11;
步骤3、在分别计算总换热系数和冷热侧与换热计算相关的参数后,得:
Figure BDA0002721411780000051
式(18)中,K表示总传热系数,rw为换热管金属热阻,Reh、Rec表示热侧和冷侧的雷诺数,Prh、Prc表示热侧和冷侧的普朗特数,λh、λc表示热侧和冷侧介质在定性温度下的导热系数,Ah、Ac表示以热侧和冷侧换热表面为基准的换热面积,d表示热侧和冷侧的水力直径;
式(18)变形为下式(19):
Figure BDA0002721411780000052
假定
Figure BDA0002721411780000053
基于一些列试验数据,通过线性回归可得到一个拟合程度最高的待定系数c0的值;
步骤4、将待定系数c1假定为另外一个值,返回步骤3,重新得到一个拟合程度最高的待定系数c0的值;
步骤5、通过步骤3及步骤4得到一系列待定系数c0的值,在这一系列的待定系数c0的值中,取其中线性回归可决系数R2最逼近1时的c0的值和c1的值为最终的系数,数据拟合完毕。
优选地,所述常数c11的取值范围为0.4~1。
本发明可以应用于冷热侧流道结构相同,冷热两侧的传热准则关联式中的各个系数的相应关系均未知,换热过程中冷热两侧Pr变化范围不大,且不方便单独研究单侧的传热系数的换热器的换热性能研究,拓宽换热器传热性能研究的实施途径。
具体实施方式
下面结合具体实施例,进一步阐述本发明。应理解,这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解,在阅读了本发明讲授的内容之后,本领域技术人员可以对本发明作各种改动或修改,这些等价形式同样落于本申请所附权利要求书所限定的范围。
本发明提出了一种适用于换热器冷热两侧的传热准则关联式中的各个系数均未知,且两侧介质流动和传热形式相似时的单侧传热系数剥离方法,拓宽换热器传热性能研究的实施途径。
本发明所涉及的换热器冷热侧流道结构(包括流道横截面形状和流道走向形式等)相同,冷热两侧的传热准则关联式中的各个系数的相应关系均未知,换热过程中冷热两侧普朗特数Pr变化范围不大,且不方便单独研究单侧的传热系数。假定待研究换热器的冷热侧传热公式相同,均为:
Nu=c0 Rec1Pr1/3 (20)
式(20)中,Nu表示传热系数,c0、c1表示待定系数,Re表示雷诺数。
换热管/壁壁温计算方法与修正威尔逊图解法相同,不作赘述。
本发明提供的一种换热器用单侧传热系数剥离方法具体包括以下步骤:
步骤1、假定待定系数c1为常数c11,常数c11的取值范围大致为0.4~1;
步骤2、在分别计算总换热系数和冷热侧与换热计算相关的参数后,得:
Figure BDA0002721411780000061
式(21)中,K表示总传热系数,rw为换热管金属热阻,Reh、Rec表示热侧和冷侧的雷诺数,Prh、Prc表示热侧和冷侧的普朗特数,λh、λc表示热侧和冷侧介质在定性温度下的导热系数,Ah、Ac表示以热侧和冷侧换热表面为基准的换热面积,d表示热侧和冷侧的水力直径;
式(21)变形为下式(22):
Figure BDA0002721411780000071
假定
Figure BDA0002721411780000072
基于一些列试验数据,通过线性回归可得到一个拟合程度最高的待定系数c0的值;
步骤3、在0.4~1的范围内,常数c11取为另外一个值,将待定系数c1假定为新的常数c11,返回步骤2,重新得到一个拟合程度最高的待定系数c0的值;
步骤4、通过步骤2及步骤3得到一系列待定系数c0的值(可以设置一个迭代阈值,当超过该迭代阈值时,停止执行步骤2及步骤3),在这一系列的待定系数c0的值中,取其中线性回归可决系数R2最逼近1时的c0的值和c1的值为最终的系数,数据拟合完毕。

Claims (2)

1.一种换热器用单侧传热系数剥离方法,所述换热器的冷热侧流道结构相同,冷热两侧的传热准则关联式中的各个系数的相应关系均未知,换热过程中,所述换热器的冷热两侧普朗特数Pr变化范围不大,其特征在于,包括以下步骤:
步骤1、设换热器的冷热侧传热公式相同,均如下式(17)所示:
Nu=c0Rec1Pr1/3 (17)
式(17)中,Nu表示传热系数,c0、c1表示待定系数,Re表示雷诺数;
步骤2、假定待定系数c1为常数c11;
步骤3、在分别计算总换热系数和冷热侧与换热计算相关的参数后,得:
Figure FDA0002721411770000011
式(18)中,K表示总传热系数,rw为换热管金属热阻,Reh、Rec表示热侧和冷侧的雷诺数,Prh、Prc表示热侧和冷侧的普朗特数,λh、λc表示热侧和冷侧介质在定性温度下的导热系数,Ah、Ac表示以热侧和冷侧换热表面为基准的换热面积,d表示热侧和冷侧的水力直径;
式(18)变形为下式(19):
Figure FDA0002721411770000012
假定
Figure FDA0002721411770000013
基于一系 列试验数据,通过线性回归可得到一个拟合程度最高的待定系数c0的值;
步骤4、将待定系数c1假定为另外一个值,返回步骤3,重新得到一个拟合程度最高的待定系数c0的值;
步骤5、通过步骤3及步骤4得到一系列待定系数c0的值,在这一系列的待定系数c0的值中,取其中线性回归可决系数R2最逼近1时的c0的值和c1的值为最终的系数,数据拟合完毕。
2.如权利要求1所述的一种换热器用单侧传热系数剥离方法,其特征在于,所述常数c11的取值范围为0.4~1。
CN202011089182.6A 2020-10-13 2020-10-13 一种换热器用单侧传热系数剥离方法 Active CN112214893B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011089182.6A CN112214893B (zh) 2020-10-13 2020-10-13 一种换热器用单侧传热系数剥离方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011089182.6A CN112214893B (zh) 2020-10-13 2020-10-13 一种换热器用单侧传热系数剥离方法

Publications (2)

Publication Number Publication Date
CN112214893A CN112214893A (zh) 2021-01-12
CN112214893B true CN112214893B (zh) 2022-09-06

Family

ID=74053288

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011089182.6A Active CN112214893B (zh) 2020-10-13 2020-10-13 一种换热器用单侧传热系数剥离方法

Country Status (1)

Country Link
CN (1) CN112214893B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113791115B (zh) * 2021-09-15 2023-10-20 湘潭大学 一种板式换热器传热性能测试方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111274672A (zh) * 2020-01-07 2020-06-12 浙江工业大学 2-4型管壳式换热器的显式传热计算方法
CN111428184A (zh) * 2020-04-30 2020-07-17 庆安集团有限公司 一种板翅式换热器芯体尺寸计算方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102135511B (zh) * 2010-12-29 2013-08-28 浙江银轮机械股份有限公司 一种板翅式换热器翅片表面传热性能的测试方法及装置
CN106484968B (zh) * 2016-09-23 2019-04-09 西安交通大学 一种基于响应面的换热压降关联式的快速预测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111274672A (zh) * 2020-01-07 2020-06-12 浙江工业大学 2-4型管壳式换热器的显式传热计算方法
CN111428184A (zh) * 2020-04-30 2020-07-17 庆安集团有限公司 一种板翅式换热器芯体尺寸计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
波纹板强化换热装置中带水膜饱和烟气对流传热的实验研究;郭勇超等;《热能动力工程》;20171120(第11期);全文 *

Also Published As

Publication number Publication date
CN112214893A (zh) 2021-01-12

Similar Documents

Publication Publication Date Title
Verma et al. ANN: Prediction of an experimental heat transfer analysis of concentric tube heat exchanger with corrugated inner tubes
Shah Assessment of modified Wilson plot techniques for obtaining heat exchanger design data
Zimparov Enhancement of heat transfer by a combination of three-start spirally corrugated tubes with a twisted tape
Rad et al. Heat transfer enhancement in shell-and-tube heat exchangers using porous media
Ameel et al. Optimization of X-shaped louvered fin and tube heat exchangers while maintaining the physical meaning of the performance evaluation criterion
CN112214893B (zh) 一种换热器用单侧传热系数剥离方法
Roy et al. Advances in heat intensification techniques in shell and tube heat exchanger
Youn et al. An experimental investigation on the airside performance of fin-and-tube heat exchangers having sinusoidal wave fins
Liu et al. Experimental and numerical investigation on the flow and heat transfer characteristics of the tube with an integrated internal longitudinal fin
Sobota Experimental and numerical analysis of heat transfer in the helically coiled heat exchanger
Kim Comparison of optimal thermal performances of finned tube annuli with various fin shapes
Conté et al. Numerical and experimental investigations of heat transfer performance of rectangular coil heat exchangers
CN116150993A (zh) 一种基于频域特性的换热器设计方法及换热器
Iqbal et al. Analysis of thermally developing laminar convection in the finned double-pipe heat exchanger
Abdelaziz et al. Modeling heat exchangers under consideration of manufacturing tolerances and uncertain flow distribution
Liu et al. Influence of fin arrangement on fluid flow and heat transfer in the inlet of a plate-fin heat exchanger
Xie et al. Numerical investigation of heat transfer and fluid flow characteristics inside a wavy channel
Zhao et al. Influence of tube arrangement on the thermal hydraulic performance of a membrane helical-coil heat exchanger
Zhang et al. Heat transfer and performance enhancement investigation of biomimetic honeycomb gas coolers in transcritical CO2 heat pumps
Al-Zahrani Heat transfer characteristics of innovative configurations of double pipe heat exchanger
CN114925598A (zh) 一种基于克里金法的管状换热器换热误差可靠性分析方法
Al-Lami et al. New pressure drop and heat transfer correlations for turbulent forced convection in internally channeled tube heat exchanger ducts
Mo et al. Study on heat transfer and fluidity of three-dimensional invisible finned-tube
Zheng et al. Heat exchanger specification: Coupling design and surface performance evaluation
CN112035957A (zh) 一种空空中冷器性能预测方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant