CN102855396A - Design calculation method of dual-stream spirally-wound tubular heat exchanger - Google Patents

Design calculation method of dual-stream spirally-wound tubular heat exchanger Download PDF

Info

Publication number
CN102855396A
CN102855396A CN201210303321XA CN201210303321A CN102855396A CN 102855396 A CN102855396 A CN 102855396A CN 201210303321X A CN201210303321X A CN 201210303321XA CN 201210303321 A CN201210303321 A CN 201210303321A CN 102855396 A CN102855396 A CN 102855396A
Authority
CN
China
Prior art keywords
flat
spiral pipe
shell
liquid
advances
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.)
Pending
Application number
CN201210303321XA
Other languages
Chinese (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201210303321XA priority Critical patent/CN102855396A/en
Publication of CN102855396A publication Critical patent/CN102855396A/en
Pending legal-status Critical Current

Links

Images

Abstract

The invention relates to a design calculation method for obtaining a dual-stream spirally-wound tubular heat exchanger complex tube bundle model and key parameters through five calculation processes, i.e. dual-stream spirally-wound tubular heat exchanger flow speed calculation, dual-stream spirally-wound tubular heat exchanger Reynolds number calculation, dual-stream spirally-wound tubular heat exchanger Prandtl number calculation, dual-stream spirally-wound tubular heat exchanger total heat transfer calculation and effective heat exchange height calculation, aiming at the internal flow heat transfer technology design calculation process of a dual-stream spirally-wound tubular heat exchanger. In the calculation processes, a logarithmic mean value method is adopted for determining inlet/outlet speed, viscosity, heat conducting coefficient, Reynolds number, Prandtl number, convection heat exchange coefficient and the like in a complex phase-change heat exchange process; a volume fraction method and a mass fraction method are used for solving two-phase average speed and other physical parameters; an effective heat exchange height calculation method is used for providing a simple and convenient calculation method for iterative calculation of the spirally-wound tubular heat exchanger; and a mathematical model for heat exchange technology design calculation of the dual-stream spirally-wound tubular heat exchanger is provided.

Description

Bifilar stream spiral winding tube type heat exchanger design and calculation method
Technical field
The present invention relates to a kind of bifilar stream spiral winding tube type heat exchanger design and calculation method, be mainly used in the low-temperature gas liquefaction separation field, comprise that the liquefaction of-161 ℃ of natural gas in low temperature ,-197 ℃ of air low temperature liquefaction separate ,-197 ℃ of low temperature liquid nitrogens are washed technique ,-70 ℃ of gas low temperatures such as low-temp methanol washing process purify, the low-temperature liquefaction separation technology field.
Background technology
Bifilar stream spiral winding tube type heat exchanger is a kind of Double-bundle helical disk cast heat-exchange apparatus that a kind of heat exchange pipeline forms after twining layer by layer, mainly consisted of by spiral pipe (1), spiral pipe (2), housing (3), core cylinder (4), support critical pieces such as (5), bifilar stream heat interchanger for basic in the wrap-round tubular heat exchanger is mainly used in the fluid heat transfer process that there is the larger temperature difference in tube side.Bifilar stream spiral winding tube type heat exchanger is with its compact conformation, unit volume has larger heat transfer area, but the thermal expansion automatic compensating of heat-transfer pipe, realize easily maximizing, can reduce the advantage such as equipment number of units and become visual plant in low temperature purification, the liquefaction process such as natural gas liquefaction, cryogenic air separation, low-temperature rectisol.Because bifilar stream spiral winding tube type heat exchanger is applied to low temperature environment mostly, internal pipeline twines complicated, there is not general design criterion, there is not unified heat-exchanging process design and calculation method yet, along with technological process or physical parameter characteristics are different and exist than big difference, brought difficulty therefore for bifilar stream spiral winding tube type heat exchanger standardisation process.In addition, because bifilar stream spiral winding tube type heat exchanger spiral pipe (1), spiral pipe (2) winding method are a lot, do not have unified pipeline winding pattern and Theoretical Design computing method to be used for the computer aided calculation process, brought obstacle for the scientific computing process of bifilar stream spiral winding tube type heat exchanger.Be standardization and the science computational problem that solves better bifilar stream spiral winding tube type heat exchanger, the present invention has provided a kind of Simplified Design computing method of bifilar stream spiral winding tube type heat exchanger from bifilar stream spiral winding tube type heat exchanger heat-exchanging process designing and calculating.
Summary of the invention
Bifilar stream spiral winding tube type heat exchanger design and calculation method comprises that bifilar stream spiral winding tube type heat exchanger flow relocity calculation, bifilar stream spiral winding tube type heat exchanger Reynolds number calculate, bifilar stream spiral winding tube type heat exchanger Prandtl number calculates, bifilar stream spiral winding tube type heat exchanger overall heat transfer coefficient is calculated and five main heat-exchanging process computation processes of the effective heat exchange high computational of bifilar stream spiral winding tube type heat exchanger.
Technical solution of the present invention:
1.Bifilar stream spiral winding tube type heat exchanger flow relocity calculation process
Shell-side calculates the average external volume flow
V O is flat=( V O advances- V O goes out)/ln ( V O advances/ V O goes out)
In the formula:
V O is flat--shell-side logarithmic mean flow, m 3/ h;
V O advances--shell-side inlet flow, m 3/ h;
V O goes out--shell-side rate of discharge, m 3/ h;
v O is flat= V O is flat/ (3600 A)
In the formula:
v O is flat--shell fluid mean flow rate, m/s;
A--shell-side passage section area, m 2
The single-phase working fluid volumetric flow rate of shell-side
V Lo is flat=( V Lo advances- V Lo goes out)/ln ( V Lo advances/ V Lo goes out)
In the formula:
V Lo advances--shell-side liquid-inlet flow, m 3/ h;
V Lo goes out--shell-side liquid outlet flow, m 3/ h;
The shell-side flow rate of liquid
v Lo is flat= V Lo is flat/ (3600 A)
In the formula:
v Lo is flat--shell-side liquid mean flow rate, m/s;
Single-phase flowing gas volumetric flow rate
V Vo is flat=( V Vo advances- V Vo goes out)/ln ( V Vo advances/ V Vo goes out)
In the formula:
V Vo advances--shell-side liquid-inlet flow, m 3/ h;
V Vo goes out--shell-side liquid outlet flow, m 3/ h;
The shell-side gas flow rate
v Vo is flat= V Vo is flat/ (3600 A)
In the formula:
v Vo is flat--shell-side gas mean flow rate, m/s;
Mean flow rate during shell-side two-phase mixed flow
v O is flat= R To v Lo is flat+ (1- R To) v Vo is flat
In the formula:
R To--liquid bulk integration rate;
R To= V Lo is flat/ ( V Lo is flat+ V Vo is flat);
V Lo is flat--shell-side liquid average external volume flow, m 3/ h;
V Vo is flat--shell-side gas average external volume flow, m 3/ h;
V Lo is flat=( V Lo advances- V Lo goes out)/ln ( V Lo advances/ V Lo goes out);
V Vo is flat=( V Vo advances- V Vo goes out)/ln ( V Vo advances/ V Vo goes out);
The interior logarithmic mean flow velocity of spiral pipe (1)
V 1i is flat=( V 1i advances- V 1i goes out)/ln ( V 1i advances/ V 1i goes out)
In the formula:
V 1i is flat--spiral pipe (1) side logarithmic mean flow, m 3/ h;
V 1i advances--spiral pipe (1) side-entrance flow, m 3/ h;
V 1i goes out--spiral pipe (1) side outlet flow, m 3/ h;
v 1i= V 1i is flat/ a i n 1
In the formula:
v 1i--spiral pipe (1) inner fluid speed, m/s;
a 1i--spiral pipe (1) interior conduit passage section area, m 2
n 1--spiral pipe (1) number of tubes;
Mean flow rate during the interior two-phase mixed flow of spiral pipe (1)
v 1i is flat= R T1i v L1i is flat+ (1- R T1i) v V1i is flat
In the formula:
R T1i--liquid bulk integration rate;
R T1i= V L1i is flat/ ( V L1i is flat+ V V1i is flat);
V L1i is flat--spiral pipe (1) side liquid average external volume flow, m 3/ h;
V V1i is flat--spiral pipe (1) side gas average external volume flow, m 3/ h;
V L1i is flat=( V L1i advances- V L1i goes out)/ln ( V L1i advances/ V L1i goes out);
V V1i is flat=( V V1i advances- V V1i goes out)/ln ( V V1i advances/ V V1i goes out);
The interior logarithmic mean flow velocity of spiral pipe (2)
V 2i is flat=( V 2i advances- V 2i goes out)/ln ( V 2i advances/ V 2i goes out)
In the formula:
V 2i is flat--spiral pipe (2) logarithmic mean flow, m 3/ h;
V 2i advances--spiral pipe (2) inlet flow rate, m 3/ h;
V 2i goes out--spiral pipe (2) rate of discharge, m 3/ h;
v 2i= v 2i is flat/ a 2i n 2
In the formula:
v 2i--spiral pipe (2) inner fluid speed, m/s;
a 2i--spiral pipe (2) interior conduit passage section area, m 2
n 2--spiral pipe (2) number of tubes;
Mean flow rate during the interior two-phase mixed flow of spiral pipe (2)
v 2i is flat= R T2i v L2i is flat+ (1- R T2i) v V2i is flat
In the formula:
R T2i--liquid bulk integration rate;
R T2i= V L2i is flat/ ( V L2i is flat+ V V2i is flat);
V L2i is flat--spiral pipe (2) side liquid average external volume flow, m 3/ h;
V V2i is flat--spiral pipe (2) side gas average external volume flow, m 3/ h;
V L2i is flat=( V L2i advances- V L2i goes out)/ln ( V L2i advances/ V L2i goes out);
V V2i is flat=( V V2i advances- V V2i goes out)/ln ( V V2i advances/ V V2i goes out);
Adopt the logarithmic mean value method in the computation process by the gas-liquid phase-splitting and mix two-phase flow averaging method calculating flow velocity;
2.Bifilar stream spiral winding tube type heat exchanger Reynolds number computation process
The shell-side Reynolds number
Re o= v O is flat ρ O is flat d o / μ O is flat
In the formula:
Re o---the shell-side Reynolds number;
v O is flat---shell fluid mean flow rate, m/s;
ρ O is flat---shell fluid average density, kg/m 3
d o---outer diameter tube, m;
μ O is flat---shell fluid average viscosity coefficient, Pa.s;
Shell-side gas average density
ρ Vo is flat=( ρ Vo advances- ρ Vo goes out)/ln ( ρ Vo advances/ ρ Vo goes out)
In the formula:
ρ Vo is flat--shell-side gas average density, kg/m 3
ρ Vo advances--shell-side gas feed average density, kg/m 3
ρ Vo goes out--shell-side gas vent average density, kg/m 3
Shell-side liquid average density
ρ Lo is flat=( ρ Lo advances- ρ Lo goes out)/ln ( ρ Lo advances/ ρ Lo goes out)
In the formula:
ρ Lo is flat--shell-side liquid average density, kg/m 3
ρ Lo advances--shell-side liquid-inlet average density, kg/m 3
ρ Lo goes out--shell-side liquid outlet average density, kg/m 3
Average density during the shell-side gas-liquid two-phase
ρ O is flat= R Lo ρ Lo is flat+ (1- R Lo) ρ Vo is flat
In the formula:
R Lo--the liquid phase mass fraction;
R Lo=( V Lo is flat/ ρ Lo is flat)/( V Lo is flat/ ρ Lo is flat+ V Vo is flat/ ρ Vo is flat);
V Lo is flat--shell-side liquid quality flow, m 3/ h;
V Vo is flat--shell-side gas mass flow, m 3/ h;
ρ Lo is flat--shell-side liquid average density, kg/m 3
ρ Vo is flat--shell-side gas average density, kg/m 3
Shell-side gas average viscosity
μ Vo is flat=( μ Vo advances- μ Vo goes out)/ln ( μ Vo advances/ μ Vo goes out)
In the formula:
μ Vo is flat--shell-side gas average viscosity;
μ Vo advances--shell-side gas feed average viscosity, Pa.s;
μ Vo goes out--shell-side gas vent average viscosity, Pa.s;
Shell-side liquid average viscosity
μ Lo is flat=( μ Lo advances- μ Vo goes out)/ln ( μ Lo advances/ μ Lo goes out)
In the formula:
μ Lo is flat--shell-side liquid average viscosity;
μ Lo advances--shell-side liquid-inlet average viscosity, Pa.s;
μ Lo goes out--shell-side liquid outlet average viscosity, Pa.s;
Shell-side gas-liquid two-phase average viscosity
μ O is flat= R Lo μ Lo is flat+ (1- R Lo) μ Vo is flat
In the formula:
R Lo--the liquid phase mass fraction;
Shell-side gas Reynolds number
Re Vo= v Vo ρ Vo is flat d o / μ Vo is flat
In the formula:
Re Vo---manage outer gas Reynolds number;
v Vo---shell-side gas flow rate, m/s;
ρ Vo is flat---shell-side gas average density, kg/m 3
d o---outer diameter tube, m;
μ Vo is flat---manage outer gas average viscosity coefficient, Pa.s;
Shell-side liquid Reynolds number
Re Lo= v Lo ρ Lo is flat d o / μ Lo is flat
In the formula:
Re Lo---manage outer liquid Reynolds number;
v Lo---shell-side flow rate of liquid, m/s;
ρ Lo is flat---shell-side liquid average density, kg/m 3
d o---outer diameter tube, m;
μ Lo is flat---manage outer liquid average viscosity coefficient, Pa.s;
Shell-side gas-liquid two-phase average Reynolds numbdr:
Re O is flat=Re Lo R Lo+ (1- R Lo) Re Vo
In the formula:
Re O is flat---the average Reynolds numbdr during gas-liquid two-phase;
R Lo---the liquid phase massfraction;
The interior gas Reynolds number of spiral pipe (1)
Re V1i= v V1i ρ V1i is flat d 1i / μ V1i is flat
In the formula:
Re V1i---the interior gas Reynolds number of spiral pipe (1);
v V1i---the interior gas flow rate of spiral pipe (1), m/s;
ρ V1i is flat---the interior gas average density of spiral pipe (1), kg/m 3
d i---spiral pipe (1) internal diameter, m;
μ V1i is flat---the interior gas average viscosity of spiral pipe (1) coefficient, Pa.s;
The interior liquid Reynolds number of spiral pipe (1)
Re L1i= v L1 i ρ L1 i is flat d 1 i / μ L1 i is flat
In the formula:
Re L1i---the interior liquid Reynolds number of spiral pipe (1);
v L1i---the interior flow rate of liquid of spiral pipe (1), m/s;
ρ L1i is flat---the interior liquid average density of spiral pipe (1), kg/m 3
d 1i---spiral pipe (1) internal diameter, m;
μ L1i is flat---the interior liquid average viscosity of spiral pipe (1) coefficient, Pa.s;
The interior gas-liquid two-phase average Reynolds numbdr of spiral pipe (1)
Re 1 i is flat=Re L1 i R L1i+ (1- R L1i) Re V1i
In the formula:
Re 1 i is flat---the average Reynolds numbdr during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas average density of spiral pipe (1)
ρ V1i is flat=( ρ V1i advances- ρ V1i goes out)/ln ( ρ V1i advances/ ρ V1i goes out)
In the formula:
ρ V1i is flat--the interior gas average density of spiral pipe (1), kg/m 3
ρ V1i advances--the interior gas feed average density of spiral pipe (1), kg/m 3
ρ V1i goes out--the interior gas vent average density of spiral pipe (1), kg/m 3
Spiral pipe 1 interior liquid average density
ρ L1i is flat=( ρ L1i advances- ρ L1i goes out)/ln ( ρ L1i advances/ ρ L1i goes out)
In the formula:
ρ L1i is flat--the interior liquid average density of spiral pipe (1), kg/m 3
ρ L1i advances--the interior liquid-inlet average density of spiral pipe (1), kg/m 3
ρ L1i goes out--the interior liquid outlet average density of spiral pipe (1), kg/m 3
Average density during gas-liquid two-phase
ρ 1i is flat= R L1i ρ L1i is flat+ (1- R L1i) ρ V1i is flat
In the formula:
R L1i--the liquid phase mass fraction;
The interior gas average viscosity of spiral pipe (1)
μ V1i is flat=( μ V1i advances- μ V1i goes out)/ln ( μ V1i advances/ μ V1i goes out)
In the formula:
μ V1i is flat--the interior liquid average viscosity of spiral pipe (1), Pa.s;
μ V1i advances--the interior liquid-inlet average viscosity of spiral pipe (1), Pa.s;
μ V1i goes out--the interior liquid outlet average viscosity of spiral pipe (1), Pa.s;
The interior liquid average viscosity of spiral pipe (1)
μ L1i is flat=( μ Lo advances- μ Vo goes out)/ln ( μ Lo advances/ μ Lo goes out)
In the formula:
μ L1i is flat--the interior liquid average viscosity of spiral pipe (1), Pa.s;
μ L1i advances--the interior liquid-inlet average viscosity of spiral pipe (1), Pa.s;
μ L1i goes out--the interior liquid outlet average viscosity of spiral pipe (1), Pa.s;
Average viscosity during gas-liquid two-phase
μ 1i is flat= R L1i μ L1i is flat+ (1- R L1i) μ V1i is flat
In the formula:
R L1i--the liquid phase mass fraction;
The interior gas Reynolds number of spiral pipe (2)
Re V2i= v V2i ρ V2i is flat d 2i / μ V2i is flat
In the formula:
Re V2i---spiral pipe (2) gas Reynolds number;
v V2i---spiral pipe (2) gas flow rate, m/s;
ρ V2i is flat---spiral pipe (2) gas average density, kg/m 3
d 2i---spiral pipe (2) internal diameter, m;
μ V2i is flat---spiral pipe (2) gas average viscosity coefficient, Pa.s;
The interior liquid Reynolds number of spiral pipe (2)
Re L2i= v L2i ρ L2i is flat d 2i / μ L2i is flat
In the formula:
Re L2i---the interior liquid Reynolds number of spiral pipe (2);
v L2i---the interior flow rate of liquid of spiral pipe (2), m/s;
ρ L2i is flat---the interior liquid average density of spiral pipe (2), kg/m 3
d 2i---spiral pipe (2) internal diameter, m;
μ L2i is flat---the interior liquid average viscosity of spiral pipe (2) coefficient, Pa.s;
Spiral pipe 2 interior gas-liquid two-phase average Reynolds numbdrs
Re 2 i are flat=Re L2 i R L2i+ (1- R L2i) Re V2i
In the formula:
Re 2i is flat---the average Reynolds numbdr during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas average density of spiral pipe (2)
ρ V2i is flat=( ρ V2i advances- ρ V2i goes out)/ln ( ρ V2i advances/ ρ V2i goes out)
In the formula:
ρ V2i is flat--the interior gas average density of spiral pipe (2), kg/m 3
ρ V2i advances--the interior gas feed average density of spiral pipe (2), kg/m 3
ρ V2i goes out--the interior gas vent average density of spiral pipe (2), kg/m 3
Spiral pipe 2 interior liquid average densities
ρ L2i is flat=( ρ L2i advances- ρ L2i goes out)/ln ( ρ L2i advances/ ρ L2i goes out)
In the formula:
ρ L2i is flat--the interior liquid average density of spiral pipe (2), kg/m 3
ρ L2i advances--the interior liquid-inlet average density of spiral pipe (2), kg/m 3
ρ L2i goes out--the interior liquid outlet average density of spiral pipe (2), kg/m 3
The interior gas-liquid two-phase density of spiral pipe (2)
ρ 2 i are flat= ρ L2 i R L2i+ (1- R L2i) ρ V2i
In the formula:
ρ 2i is flat---the average density during gas-liquid two-phase, kg/m 3
R L2i---the liquid phase massfraction;
The interior gas average viscosity of spiral pipe (2)
μ V2i is flat=( μ V2i advances- μ V2i goes out)/ln ( μ V2i advances/ μ V2i goes out)
In the formula:
μ V2i is flat--the interior liquid average viscosity of spiral pipe (2), Pa.s;
μ V2i advances--the interior liquid-inlet average viscosity of spiral pipe (2), Pa.s;
μ V2i goes out--the interior liquid outlet average viscosity of spiral pipe (2), Pa.s;
The interior liquid average viscosity of spiral pipe (2)
μ L2i is flat=( μ L2i advances- μ V2i goes out)/ln ( μ L2i advances/ μ L2i goes out)
In the formula:
μ L2i is flat--the interior liquid average viscosity of spiral pipe (2), Pa.s;
μ L2i advances--the interior liquid-inlet average viscosity of spiral pipe (2), Pa.s;
μ L2i goes out--the interior liquid outlet average viscosity of spiral pipe (2), Pa.s;
Average viscosity during the interior gas-liquid two-phase of spiral pipe (2)
μ 2i is flat= R L2i μ L2i is flat+ (1- R L2i) μ V2i is flat
In the formula:
R L2i--the liquid phase mass fraction;
Adopt the Reynolds number mean value method in the computation process by the gas-liquid phase-splitting and two-phase flow mass fraction calculating after mixing;
3.Bifilar stream spiral winding tube type heat exchanger Prandtl number computation process
The shell-side Prandtl number
Pr o= C Po μ O is flat/ λ O is flat
In the formula:
Pr o---the shell-side Prandtl number;
C Po---shell-side specific heat at constant pressure, kJ/kg;
μ O is flat---shell-side average viscosity coefficient, Pa.s;
λ O is flat---mean coefficient of heat conductivity, kg/m 3
Shell-side gas average specific heat at constant pressure
C Pvo is flat=( C Pvo advances- C Pvo goes out)/ln ( C Pvo advances/ C Pvo goes out)
In the formula:
C Pvo is flat--shell-side gas average density, kJ/kg;
C Pvo advances--shell-side gas feed average density, kJ/kg;
C Pvo goes out--shell-side gas vent average density, kJ/kg;
Shell-side liquid average specific heat at constant pressure
C PLo is flat=( C PLo advances- C PLo goes out)/ln ( C PLo advances/ C PLo goes out)
In the formula:
C PLo is flat--shell-side liquid average density, kJ/kg;
C PLo advances--shell-side liquid-inlet average density, kJ/kg;
C PLo goes out--shell-side liquid outlet average density, kJ/kg;
Shell-side gas-liquid two-phase average specific heat at constant pressure
C Po is flat= R Lo C PLo is flat+ (1- R Lo) C Pvo is flat
In the formula:
R Lo--liquid quality divides rate;
Shell-side gas average viscosity
μ Vo is flat=( μ Vo advances- ρ Vo goes out)/ln ( μ Vo advances/ μ Vo goes out)
In the formula:
μ Vo is flat--shell-side gas average viscosity;
μ Vo advances--shell-side gas feed average viscosity, Pa.s;
μ Vo goes out--shell-side gas vent average viscosity, Pa.s;
Shell-side liquid average viscosity
μ Lo is flat=( μ Lo advances- ρ Vo goes out)/ln ( μ Lo advances/ μ Lo goes out)
In the formula:
μ Lo is flat--shell-side liquid average viscosity;
μ Lo advances--shell-side liquid-inlet average viscosity, Pa.s;
μ Lo goes out--shell-side liquid outlet average viscosity, Pa.s;
Shell-side gas-liquid two-phase average viscosity
μ O is flat= R Lo μ Lo is flat+ (1- R Lo) μ Vo is flat
In the formula:
R Lo--liquid quality divides rate;
Shell-side gas mean coefficient of heat conductivity
λ Vo is flat=( λ Vo advances- ρ Vo goes out)/ln ( λ Vo advances/ λ Vo goes out)
In the formula:
λ Vo is flat--shell side gas mean coefficient of heat conductivity, W/ (m.K);
λ Vo advances--shell side gas feed mean coefficient of heat conductivity, W/ (m.K);
λ Vo goes out--shell side gas vent mean coefficient of heat conductivity, W/ (m.K);
Shell-side liquid mean coefficient of heat conductivity
λ Lo is flat=( λ Lo advances- λ Vo goes out)/ln ( λ Lo advances/ λ Lo goes out)
In the formula:
λ Lo is flat--shell-side liquid mean coefficient of heat conductivity, W (m.K);
λ Lo advances--shell-side liquid-inlet mean coefficient of heat conductivity, W/ (m.K);
λ Lo goes out--shell-side liquid outlet mean coefficient of heat conductivity, W/ (m.K);
Shell-side gas-liquid two-phase mean coefficient of heat conductivity
λ O is flat= R Lo λ Lo is flat+ (1- R Lo) λ Vo is flat
In the formula:
R Lo--liquid quality divides rate;
Shell-side gas Prandtl number
Pr Vo= C Pvo μ Vo is flat/ λ Vo is flat
In the formula:
Pr Vo---the shell-side Prandtl number;
C Pvo---shell-side specific heat at constant pressure, kJ/kg;
μ Vo is flat---shell-side average viscosity coefficient, Pa.s;
λ Vo is flat---mean coefficient of heat conductivity, kg/m 3
Shell-side liquid Prandtl number
Pr Lo= C PLo μ Lo is flat/ λ Lo is flat
In the formula:
Pr Lo---the shell-side Prandtl number;
C Lpo---shell-side specific heat at constant pressure, kJ/kg;
μ Lo is flat---shell-side average viscosity coefficient, Pa.s;
λ Lo is flat---mean coefficient of heat conductivity, kg/m 3
The average Prandtl number of shell-side gas-liquid two-phase
Pr O is flat=Pr Lo R Lo+ (1- R Lo) Pr Vo
In the formula:
Pr O is flat---the average Prandtl number during gas-liquid two-phase;
R Lo---the liquid phase massfraction;
The interior gas Prandtl number of spiral pipe (1)
Pr V1i= C Pv1i μ V1i is flat/ λ V1i is flat
In the formula:
Pr V1i---the shell-side Prandtl number;
C Pv1i---shell-side specific heat at constant pressure, kJ/kg;
μ V1i is flat---shell-side average viscosity coefficient, Pa.s;
λ V1i is flat---mean coefficient of heat conductivity, kg/m 3
The interior liquid Prandtl number of spiral pipe (1)
Pr L1i= C PL1i μ L1i is flat/ λ L1i is flat
In the formula:
Pr L1i---the shell-side Prandtl number;
C PL1i---shell-side specific heat at constant pressure, kJ/kg;
μ L1i is flat---shell-side average viscosity coefficient, Pa.s;
λ L1i is flat---mean coefficient of heat conductivity, kg/m 3
The average Prandtl number of the interior gas-liquid two-phase of spiral pipe (1)
Pr 1i is flat=Pr L1i R L1i+ (1- R L1i) Pr V1i
In the formula:
Pr 1i is flat---the average Prandtl number during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas average specific heat at constant pressure of spiral pipe (1)
C Pv1i is flat=( C Pv1i advances- C Pv1i goes out)/ln ( C Pv1i advances/ C Pv1i goes out)
In the formula:
C Pv1i is flat--shell-side gas average density, kJ/kg;
C Pv1i advances--shell-side gas feed average density, kJ/kg;
C Pv1i goes out--shell-side gas vent average density, kJ/kg;
The interior liquid average specific heat at constant pressure of spiral pipe (1)
C PL1i is flat=( C PL1i advances- C PL1i goes out)/ln ( C PL1i advances/ C PL1i goes out)
In the formula:
C PL1i is flat--shell-side liquid average density, kJ/kg;
C PL1i advances--shell-side liquid-inlet average density, kJ/kg;
C PL1i goes out--shell-side liquid outlet average density, kJ/kg;
The interior gas-liquid two-phase average specific heat at constant pressure of spiral pipe (1)
C P1i is flat= C PL1i R L1i+ (1- R L1i) C Pv1i
In the formula:
C P1i is flat---the average specific heat at constant pressure during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas average viscosity of spiral pipe (1)
μ V1i is flat=( μ V1i advances- ρ V1i goes out)/ln ( μ V1i advances/ μ V1i goes out)
In the formula:
μ V1i is flat--shell-side gas average viscosity, Pa.s;
μ V1i advances--shell-side gas feed average viscosity, Pa.s;
μ V1i goes out--shell-side gas vent average viscosity, Pa.s;
The interior liquid average viscosity of spiral pipe (1)
μ L1i is flat=( μ L1i advances- ρ L1i goes out)/ln ( μ L1i advances/ μ L1i goes out)
In the formula:
μ L1i is flat--shell-side liquid average viscosity, Pa.s;
μ L1i advances--shell-side liquid-inlet average viscosity, Pa.s;
μ L1i goes out--shell-side liquid outlet average viscosity, Pa.s;
The interior gas-liquid two-phase average viscosity of spiral pipe (1)
μ 1i is flat= μ L1i R L1i+ (1- R L1i) μ V1i
In the formula:
μ 1i is flat---the average viscosity during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas mean coefficient of heat conductivity of spiral pipe (1)
λ V1i is flat=( λ V1i advances- ρ V1i goes out)/ln ( λ V1i advances/ λ V1i goes out)
In the formula:
λ V1i is flat--shell-side gas mean coefficient of heat conductivity, W/ (m.K);
λ V1i advances--shell-side gas feed mean coefficient of heat conductivity, W/ (m.K);
λ V1i goes out--shell-side gas vent mean coefficient of heat conductivity, W/ (m.K);
The interior liquid mean coefficient of heat conductivity of spiral pipe (1)
λ L1i is flat=( λ L1i advances- λ V1i goes out)/ln ( λ L1i advances/ λ L1i goes out)
In the formula:
λ L1i is flat--shell-side liquid mean coefficient of heat conductivity, W/ (m.K);
λ L1i advances--shell-side liquid-inlet mean coefficient of heat conductivity, W/ (m.K);
λ L1i goes out--shell-side liquid outlet mean coefficient of heat conductivity, W/ (m.K);
The interior gas-liquid two-phase mean coefficient of heat conductivity of spiral pipe (1)
λ 1i is flat= λ L1i R L1i+ (1- R L1i) λ V1i
In the formula:
λ 1i is flat---the mean coefficient of heat conductivity during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas Prandtl number of spiral pipe (2)
Pr V2i= C Pv2i μ V2i is flat/ λ V2i is flat
In the formula:
Pr V2i---shell-side Prandtl number, m/s;
C Pv2i---shell-side specific heat at constant pressure, kJ/kg;
μ V2i is flat---shell-side average viscosity coefficient, Pa.s;
λ V2i is flat---mean coefficient of heat conductivity, kg/m 3
The interior liquid Prandtl number of spiral pipe (2)
Pr L2i= C PL2i μ L2i is flat/ λ L2i is flat
In the formula:
Pr L2i---shell-side Prandtl number, m/s;
C PL2i---shell-side specific heat at constant pressure, kJ/kg;
μ L2i is flat---shell-side average viscosity coefficient, Pa.s;
λ L2i is flat---mean coefficient of heat conductivity, kg/m 3
The average Prandtl number of the interior gas-liquid two-phase of spiral pipe (2)
Pr 2i is flat=Pr L2i R L2i+ (1- R L2i) Pr V2i
In the formula:
Pr 2i is flat---the average Prandtl number during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas average specific heat at constant pressure of spiral pipe (2)
C Pv2i is flat=( C Pv2i advances- C Pv2i goes out)/ln ( C Pv2i advances/ C Pv2i goes out)
In the formula:
C Pv2i is flat--the interior gas average density of spiral pipe (2), kJ/kg;
C Pv2i advances--the interior gas feed average density of spiral pipe (2), kJ/kg;
C Pv2i goes out--the interior gas vent average density of spiral pipe (2), kJ/kg;
The interior liquid average specific heat at constant pressure of spiral pipe (2)
C PL2i is flat=( C PL2i advances- C PL2i goes out)/ln ( C PL2i advances/ C PL2i goes out)
In the formula:
C PL2i is flat--the interior liquid average density of spiral pipe (2), kJ/kg;
C PL2i advances--the interior liquid-inlet average density of spiral pipe (2), kJ/kg;
C PL2i goes out--the interior liquid outlet average density of spiral pipe (2), kJ/kg;
The interior gas-liquid two-phase average specific heat at constant pressure of spiral pipe (2)
C P2i is flat= C PL2i R L2i+ (1- R L2i) C Pv2i
In the formula:
C P2i is flat---the average specific heat at constant pressure during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas average viscosity of spiral pipe (2)
μ V2i is flat=( μ V2i advances- ρ V2i goes out)/ln ( μ V2i advances/ μ V2i goes out)
In the formula:
μ V2i is flat--the interior gas average viscosity of spiral pipe (2), Pa.s;
μ V2i advances--the interior gas feed average viscosity of spiral pipe (2), Pa.s;
μ V2i goes out--the interior gas vent average viscosity of spiral pipe (2), Pa.s;
Spiral pipe 2 interior liquid average viscosities
μ L2i is flat=( μ L2i advances- ρ L2i goes out)/ln ( μ L2i advances/ μ L2i goes out)
In the formula:
μ L2i is flat--the interior liquid average viscosity of spiral pipe (2), Pa.s;
μ L2i advances--the interior liquid-inlet average viscosity of spiral pipe (2), Pa.s;
μ L2i goes out--the interior liquid outlet average viscosity of spiral pipe (2), Pa.s;
Spiral pipe 2 interior gas-liquid two-phase average viscosities
μ 2i is flat= μ L2i R L2i+ (1- R L2i) μ V2i
In the formula:
μ 2i is flat---the average viscosity during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas mean coefficient of heat conductivity of spiral pipe (2)
λ V2i is flat=( λ V2i advances- ρ V2i goes out)/ln ( λ V2i advances/ λ V2i goes out)
In the formula:
λ V2i is flat--the interior gas mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ V2i advances--the interior gas feed mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ V2i goes out--the interior gas vent mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
The interior liquid mean coefficient of heat conductivity of spiral pipe (2)
λ L2i is flat=( λ L2i advances- λ V2i goes out)/ln ( λ L2i advances/ λ L2i goes out)
In the formula:
λ L2i is flat--the interior liquid mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ L2i advances--the interior liquid-inlet mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ L2i goes out--the interior liquid outlet mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
The interior gas-liquid two-phase mean coefficient of heat conductivity of spiral pipe (2)
λ 2i is flat= λ L2i R L2i+ (1- R L2i) λ V2i
In the formula:
λ 2i is flat---the mean coefficient of heat conductivity during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
Adopt the Prandtl number mean value method by the gas-liquid phase-splitting and mix the calculating of two-phase flow mass fraction;
4.Bifilar stream spiral winding tube type heat exchanger overall heat transfer coefficient computation process
The shell-side convection transfer rate
h o= Q Always/ ( Q s/ h s+ Q r/ h r)
h o--manage outer convection transfer rate, W/ (m 2.K);
Q Always--overall heat exchange amount, W;
Q s--sensible heat, W;
Q r--latent heat, W;
h s--sensible heat convection transfer rate, W/ (m 2.K);
h s?= R Lo h Lo+(1－ R Lo) ?h vo
In the formula:
R Lo--the liquid phase mass fraction;
h Lo--liquid sensible heat convection transfer rate, W/ (m 2.K);
h Vo--gas sensible heat convection transfer rate, W/ (m 2.K);
h r--latent heat convection transfer rate, W/ (m 2.K);
h r=1000～1100, get h r=1100 W/ (m 2.K);
Convection transfer rate
h Vo=0.296 ( λ Vo is flat/ d o) Re Vo is flat 0.609Pr Vo is flat 0.31
h Lo=0.296 ( λ Lo is flat/ d o) Re Lo is flat 0.609Pr Lo is flat 0.31
Calculate the heat exchange amount
Q Always= C Po is flat mΔ t+ m c r o
In the formula:
m--shell-side oeverall quality flow, kg/s;
m c--gas-liquid phase transition quality, kg/s;
r o--shell-side vaporizing liquid latent heat, kJ/kg;
Q s=( C PLo is flat m Lo is flat+ C Pvo is flat m Vo is flat) Δ t
In the formula:
m Lo is flat--shell-side liquid quality flow, kg/s;
m Lo is flat=( m Lo advances- m Lo goes out)/ln ( m Lo advances/ m Lo goes out)
m Vo is flat--shell-side gas mass flow, kg/s;
m Vo is flat=( m Vo advances- m Vo goes out)/ln ( m Vo advances/ m Vo goes out)
Q r =Q Always -Q s
The interior convection transfer rate of spiral pipe (1)
h 1i =0.037 ( λ Li is flat / d Li) (Re Li is flat 0.75) Pr Li is flat 0.42
Spiral pipe 2 interior convection transfer rates
h 2i =0.037 ( λ 2i is flat / d 2i) (Re 2i is flat 0.75) Pr 2i is flat 0.42
The corresponding overall heat transfer coefficient of spiral pipe (1)
K 1=1/(1/ h o+ R 1o+ δ 1 d 1o/( λ 1 d 1m)+ R 1i d 1o/ d 1i+ d 1o/( h 1i d 1i))
In the formula:
K 1--spiral pipe (1) overall heat transfer coefficient, W/ (m 2.K);
h o--shell-side convection transfer rate, W/ (m 2.K);
h 1i--the inboard convection transfer rate of spiral pipe (1), W/ (m 2.K);
R 1o--the dirty coefficient of shell-side outer wall heat, W/ (m 2.K);
R 1i--spiral pipe (1) inner surface heat dirt coefficient, W/ (m 2.K);
λ 1--spiral pipe (1) tube wall heat conduction coefficient, W/ (m.K);
δ 1--spiral pipe (1) pipeline wall thickness, m;
d 1m--spiral pipe (1) pipeline central diameter, m;
d 1o--spiral pipe (1) outer diameter tube, m;
d 1i--spiral pipe (1) internal diameter of the pipeline, m;
The corresponding overall heat transfer coefficient of spiral pipe (2)
K 2=1/(1/ h o+ R o+ δ 2 d o/( λd m)+ R i d o/ d i+ d o/( h 2i d i))
In the formula:
K 2--spiral pipe (2) overall heat transfer coefficient, W/ (m 2.K);
h o--shell-side convection transfer rate, W/ (m 2.K);
h 2i--the inboard convection transfer rate of spiral pipe (2), W/ (m 2.K);
R 2o--the dirty coefficient of shell-side outer wall heat, W/ (m 2.K);
R 2i--spiral pipe (2) inner surface heat dirt coefficient, W/ (m 2.K);
λ 2--spiral pipe (2) tube wall heat conduction coefficient, W/ (m.K);
δ 2--spiral pipe (2) pipeline wall thickness, m;
d 2m--spiral pipe (2) pipeline central diameter, m;
d 2o--spiral pipe (2) outer diameter tube, m;
d 2i--spiral pipe (2) internal diameter of the pipeline, m;
Adopt the logarithmic mean value method by gas-liquid phase-splitting and sensible heat and latent heat independence computing method calculating convection transfer rate and overall heat transfer coefficient;
5.The effective heat exchange high computational of bifilar stream spiral winding tube type heat exchanger process
Total heat transfer in the spiral winding tube type heat exchanger
Q Always =A Always K AlwaysΔ t
In the formula:
Δ t--log-mean temperature difference, K;
Spiral pipe 1 total heat conduction area
A 1 is total =Q 1 is total/ ( K 1 is totalΔ t 1)
Spiral pipe (1) length
l 1 =A 1 is total/ ( n 1 π d 1m)
In the formula:
l 1--spiral pipe (1) length;
Spiral pipe (2) total heat conduction area
A 2 is total =Q 2 is total/ ( K 2 is totalΔ t 2)
Spiral pipe (2) length
l 2 =A 2 is total/ ( n 2 π d 2m)
In the formula:
l 2--spiral pipe (2) length;
l 1= l 2The time computation process finish;
The effective heat exchange height of heat interchanger
l ?3＝ l 1sin a＝ l 2sin a
In the formula:
a--the spiral pipe ascending angle;
l 1≠ l 2The time readjust the spiral pipe inlet velocity v 1i, v 2iThe value size, and recomputate whole flow process, until l 1= l 2The time finish heat-exchanging process computation process.
Bifilar stream spiral winding tube type heat exchanger housing (3) and other accessory can carry out designing and calculating with reference to design standardss such as GB150, GB151.
The Principle Problems that scheme is related:
Bifilar stream spiral winding tube type heat exchanger is mainly used in low-temperature gas liquefaction to be separated and the gas purification field, such as technical fields such as gas low temperature purification, low-temperature liquefaction separation such as natural gas in low temperature liquefaction, air low temperature liquefaction separation, low-temperature rectisols.Because have fierce phase transformation and polyphasic flow flow process in the low-temperature heat exchange process, the physical parameter amplitude of variation is very large in the heat transfer process, general heat exchange method for designing all is difficult to accurate calculating, does not also have design and calculation method systematic, standard at present; The present invention widely applies logarithmic mean value method, volume fraction method and mass fraction method solution procedure parameter according to thermal conduction study and principles of fluid mechanics, and parameters obtained is applied to overall heat transfer coefficient computation process; Bifilar stream spiral winding tube type heat exchanger heat-exchanging process design and calculation method is by basic engineering parameter makeover process, can be applicable to complicated process of mathematical modeling of twining tube bank, and with the heat-exchanging process design calculation process of building the application of mathematical model in spiral winding tube type heat exchanger, obtain the complicated Pipe bundle structure parameter of bifilar stream spiral winding tube type heat exchanger and relevant heat exchange computation model, design whole bifilar stream spiral winding tube type heat exchanger with this, make bifilar stream spiral winding tube type heat exchanger that a definite design and calculation method be arranged, be conducive to the standardisation process of bifilar stream spiral winding tube type heat exchanger.
Technical characterstic of the present invention:
Be that bifilar stream spiral winding tube type heat exchanger flow relocity calculation, bifilar stream spiral winding tube type heat exchanger Reynolds number calculate, bifilar stream spiral winding tube type heat exchanger Prandtl number calculates, bifilar stream spiral winding tube type heat exchanger overall heat transfer coefficient is calculated and five main heat-exchanging process computation processes of the effective heat exchange high computational of bifilar stream spiral winding tube type heat exchanger by five-step approach, obtain bifilar stream spiral winding tube type heat exchanger complicated tube bank model and key parameter design and calculation method; The velocity computation process of refinement bifilar stream spiral tube side and shell side in the computation process is used the logarithmic mean value method and is determined complicated phase-change heat transfer process import and export speed, viscosity, coefficient of heat conductivity, Reynolds number, Prandtl number, convection transfer rate etc.; Use volume fraction method and mass fraction method and find the solution gas-liquid two-phase mean speed and other physical parameter; Provided the mathematics computing model that calculates bifilar stream spiral winding tube type heat exchanger overall heat transfer coefficient; Using effective heat exchange altitude gauge algorithm provides the simple calculating method of iterative computation spiral winding tube type heat exchanger; The complete bifilar stream wrap-round tubular heat exchanger heat-exchanging process computation model of one cover is proposed, the winding method of bifilar stream spiral winding tube type heat exchanger can be applied to modeling process, and with the building mathematics computing model technology Calculation process that is applied to conduct heat, obtain the bifilar stream spiral winding tube type heat exchanger plumber computation model of planting, design whole bifilar stream spiral winding tube type heat exchanger with this, make bifilar stream spiral winding tube type heat exchanger that a definite design and calculation method be arranged.
Description of drawings
Figure 1 shows that bifilar stream spiral winding tube type heat exchanger critical piece pie graph.
Embodiment
Adopt the logarithmic mean value method by the gas-liquid phase-splitting and mix two-phase flow averaging method calculating flow velocity; Adopt mean value method by the gas-liquid phase-splitting and two-phase flow mass fraction method calculating Reynolds number after mixing; Adopt mean value method by the gas-liquid phase-splitting and mix two-phase flow mass fraction calculating Prandtl number; Adopt mean value method by gas-liquid phase-splitting and sensible heat and latent heat independence computing method calculating convection transfer rate; Divide algorithm to determine total heat transfer, total heat conduction area, spiral pipe (1) length, spiral pipe (2) length according to two plumes; l 1= l 2The time determine the effective heat exchange of spiral winding tube type heat exchanger height; l 1≠ l 2The time readjust the spiral pipe inlet velocity v 1i, v 2iThe value size also recomputates whole flow process, until l 1= l 2The time determine again the effective heat exchange height of spiral winding tube type heat exchanger; With reference to the bifilar stream spiral winding tube type heat exchanger of design standards designing and calculating housing (3) and other accessories such as GB150, GB151 and finish the design calculation process of whole bifilar stream spiral winding tube type heat exchanger.

Claims (6)

1. bifilar stream spiral winding tube type heat exchanger design and calculation method is characterized in that: bifilar stream spiral winding tube type heat exchanger design and calculation method comprises that bifilar stream spiral winding tube type heat exchanger flow relocity calculation, bifilar stream spiral winding tube type heat exchanger Reynolds number calculate, bifilar stream spiral winding tube type heat exchanger Prandtl number calculates, bifilar stream spiral winding tube type heat exchanger overall heat transfer coefficient is calculated and five main heat-exchanging process computation processes of the effective heat exchange high computational of bifilar stream spiral winding tube type heat exchanger.
2. according to claim 1Described bifilar stream spiral winding tube type heat exchanger flow relocity calculation process is characterized in that:
Shell-side calculates the average external volume flow
V O is flat=( V O advances- V O goes out)/ln ( V O advances/ V O goes out)
In the formula:
V O is flat--shell-side logarithmic mean flow, m 3/ h;
V O advances--shell-side inlet flow, m 3/ h;
V O goes out--shell-side rate of discharge, m 3/ h;
v O is flat= V O is flat/ (3600 A)
In the formula:
v O is flat--shell fluid mean flow rate, m/s;
A--shell-side passage section area, m 2
The single-phase working fluid volumetric flow rate of shell-side
V Lo is flat=( V Lo advances- V Lo goes out)/ln ( V Lo advances/ V Lo goes out)
In the formula:
V Lo advances--shell-side liquid-inlet flow, m 3/ h;
V Lo goes out--shell-side liquid outlet flow, m 3/ h;
The shell-side flow rate of liquid
v Lo is flat= V Lo is flat/ (3600 A)
In the formula:
v Lo is flat--shell-side liquid mean flow rate, m/s;
Single-phase flowing gas volumetric flow rate
V Vo is flat=( V Vo advances- V Vo goes out)/ln ( V Vo advances/ V Vo goes out)
In the formula:
V Vo advances--shell-side liquid-inlet flow, m 3/ h;
V Vo goes out--shell-side liquid outlet flow, m 3/ h;
The shell-side gas flow rate
v Vo is flat= V Vo is flat/ (3600 A)
In the formula:
v Vo is flat--shell-side gas mean flow rate, m/s;
Mean flow rate during shell-side two-phase mixed flow
v O is flat= R To v Lo is flat+ (1- R To) v Vo is flat
In the formula:
R To--liquid bulk integration rate;
R To= V Lo is flat/ ( V Lo is flat+ V Vo is flat);
V Lo is flat--shell-side liquid average external volume flow, m 3/ h;
V Vo is flat--shell-side gas average external volume flow, m 3/ h;
V Lo is flat=( V Lo advances- V Lo goes out)/ln ( V Lo advances/ V Lo goes out);
V Vo is flat=( V Vo advances- V Vo goes out)/ln ( V Vo advances/ V Vo goes out);
The interior logarithmic mean flow velocity of spiral pipe (1)
V 1i is flat=( V 1i advances- V 1i goes out)/ln ( V 1i advances/ V 1i goes out)
In the formula:
V 1i is flat--spiral pipe (1) side logarithmic mean flow, m 3/ h;
V 1i advances--spiral pipe (1) side-entrance flow, m 3/ h;
V 1i goes out--spiral pipe (1) side outlet flow, m 3/ h;
v 1i= V 1i is flat/ a i n 1
In the formula:
v 1i--spiral pipe (1) inner fluid speed, m/s;
a 1i--spiral pipe (1) interior conduit passage section area, m 2
n 1--spiral pipe (1) number of tubes;
Mean flow rate during the interior two-phase mixed flow of spiral pipe (1)
v 1i is flat= R T1i v L1i is flat+ (1- R T1i) v V1i is flat
In the formula:
R T1i--liquid bulk integration rate;
R T1i= V L1i is flat/ ( V L1i is flat+ V V1i is flat);
V L1i is flat--spiral pipe (1) side liquid average external volume flow, m 3/ h;
V V1i is flat--spiral pipe (1) side gas average external volume flow, m 3/ h;
V L1i is flat=( V L1i advances- V L1i goes out)/ln ( V L1i advances/ V L1i goes out);
V V1i is flat=( V V1i advances- V V1i goes out)/ln ( V V1i advances/ V V1i goes out);
The interior logarithmic mean flow velocity of spiral pipe (2)
V 2i is flat=( V 2i advances- V 2i goes out)/ln ( V 2i advances/ V 2i goes out)
In the formula:
V 2i is flat--spiral pipe (2) logarithmic mean flow, m 3/ h;
V 2i advances--spiral pipe (2) inlet flow rate, m 3/ h;
V 2i goes out--spiral pipe (2) rate of discharge, m 3/ h;
v 2i= v 2i is flat/ a 2i n 2
In the formula:
v 2i--spiral pipe (2) inner fluid speed, m/s;
a 2i--spiral pipe (2) interior conduit passage section area, m 2
n 2--spiral pipe (2) number of tubes;
Mean flow rate during the interior two-phase mixed flow of spiral pipe (2)
v 2i is flat= R T2i v L2i is flat+ (1- R T2i) v V2i is flat
In the formula:
R T2i--liquid bulk integration rate;
R T2i= V L2i is flat/ ( V L2i is flat+ V V2i is flat);
V L2i is flat--spiral pipe (2) side liquid average external volume flow, m 3/ h;
V V2i is flat--spiral pipe (2) side gas average external volume flow, m 3/ h;
V L2i is flat=( V L2i advances- V L2i goes out)/ln ( V L2i advances/ V L2i goes out);
V V2i is flat=( V V2i advances- V V2i goes out)/ln ( V V2i advances/ V V2i goes out);
Adopt the logarithmic mean value method to calculate by the gas-liquid phase-splitting and mixing two-phase flow averaging method calculating flow velocity in the computation process.
3. according to claim 1Described bifilar stream spiral winding tube type heat exchanger Reynolds number computation process is characterized in that:
The shell-side Reynolds number
Re o= v O is flat ρ O is flat d o / μ O is flat
In the formula:
Re o---the shell-side Reynolds number;
v O is flat---shell fluid mean flow rate, m/s;
ρ O is flat---shell fluid average density, kg/m 3
d o---outer diameter tube, m;
μ O is flat---shell fluid average viscosity coefficient, Pa.s;
Shell-side gas average density
ρ Vo is flat=( ρ Vo advances- ρ Vo goes out)/ln ( ρ Vo advances/ ρ Vo goes out)
In the formula:
ρ Vo is flat--shell-side gas average density, kg/m 3
ρ Vo advances--shell-side gas feed average density, kg/m 3
ρ Vo goes out--shell-side gas vent average density, kg/m 3
Shell-side liquid average density
ρ Lo is flat=( ρ Lo advances- ρ Lo goes out)/ln ( ρ Lo advances/ ρ Lo goes out)
In the formula:
ρ Lo is flat--shell-side liquid average density, kg/m 3
ρ Lo advances--shell-side liquid-inlet average density, kg/m 3
ρ Lo goes out--shell-side liquid outlet average density, kg/m 3
Average density during the shell-side gas-liquid two-phase
ρ O is flat= R Lo ρ Lo is flat+ (1- R Lo) ρ Vo is flat
In the formula:
R Lo--the liquid phase mass fraction;
R Lo=( V Lo is flat/ ρ Lo is flat)/( V Lo is flat/ ρ Lo is flat+ V Vo is flat/ ρ Vo is flat);
V Lo is flat--shell-side liquid quality flow, m 3/ h;
V Vo is flat--shell-side gas mass flow, m 3/ h;
ρ Lo is flat--shell-side liquid average density, kg/m 3
ρ Vo is flat--shell-side gas average density, kg/m 3
Shell-side gas average viscosity
μ Vo is flat=( μ Vo advances- μ Vo goes out)/ln ( μ Vo advances/ μ Vo goes out)
In the formula:
μ Vo is flat--shell-side gas average viscosity;
μ Vo advances--shell-side gas feed average viscosity, Pa.s;
μ Vo goes out--shell-side gas vent average viscosity, Pa.s;
Shell-side liquid average viscosity
μ Lo is flat=( μ Lo advances- μ Vo goes out)/ln ( μ Lo advances/ μ Lo goes out)
In the formula:
μ Lo is flat--shell-side liquid average viscosity;
μ Lo advances--shell-side liquid-inlet average viscosity, Pa.s;
μ Lo goes out--shell-side liquid outlet average viscosity, Pa.s;
Shell-side gas-liquid two-phase average viscosity
μ O is flat= R Lo μ Lo is flat+ (1- R Lo) μ Vo is flat
In the formula:
R Lo--the liquid phase mass fraction;
Shell-side gas Reynolds number
Re Vo= v Vo ρ Vo is flat d o / μ Vo is flat
In the formula:
Re Vo---manage outer gas Reynolds number;
v Vo---shell-side gas flow rate, m/s;
ρ Vo is flat---shell-side gas average density, kg/m 3
d o---outer diameter tube, m;
μ Vo is flat---manage outer gas average viscosity coefficient, Pa.s;
Shell-side liquid Reynolds number
Re Lo= v Lo ρ Lo is flat d o / μ Lo is flat
In the formula:
Re Lo---manage outer liquid Reynolds number;
v Lo---shell-side flow rate of liquid, m/s;
ρ Lo is flat---shell-side liquid average density, kg/m 3
d o---outer diameter tube, m;
μ Lo is flat---manage outer liquid average viscosity coefficient, Pa.s;
Shell-side gas-liquid two-phase average Reynolds numbdr:
Re O is flat=Re Lo R Lo+ (1- R Lo) Re Vo
In the formula:
Re O is flat---the average Reynolds numbdr during gas-liquid two-phase;
R Lo---the liquid phase massfraction;
The interior gas Reynolds number of spiral pipe (1)
Re V1i= v V1i ρ V1i is flat d 1i / μ V1i is flat
In the formula:
Re V1i---the interior gas Reynolds number of spiral pipe (1);
v V1i---the interior gas flow rate of spiral pipe (1), m/s;
ρ V1i is flat---the interior gas average density of spiral pipe (1), kg/m 3
d i---spiral pipe (1) internal diameter, m;
μ V1i is flat---the interior gas average viscosity of spiral pipe (1) coefficient, Pa.s;
The interior liquid Reynolds number of spiral pipe (1)
Re L1i= v L1 i ρ L1 i is flat d 1 i / μ L1 i is flat
In the formula:
Re L1i---the interior liquid Reynolds number of spiral pipe (1);
v L1i---the interior flow rate of liquid of spiral pipe (1), m/s;
ρ L1i is flat---the interior liquid average density of spiral pipe (1), kg/m 3
d 1i---spiral pipe (1) internal diameter, m;
μ L1i is flat---the interior liquid average viscosity of spiral pipe (1) coefficient, Pa.s;
The interior gas-liquid two-phase average Reynolds numbdr of spiral pipe (1)
Re 1 i is flat=Re L1 i R L1i+ (1- R L1i) Re V1i
In the formula:
Re 1 i is flat---the average Reynolds numbdr during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas average density of spiral pipe (1)
ρ V1i is flat=( ρ V1i advances- ρ V1i goes out)/ln ( ρ V1i advances/ ρ V1i goes out)
In the formula:
ρ V1i is flat--the interior gas average density of spiral pipe (1), kg/m 3
ρ V1i advances--the interior gas feed average density of spiral pipe (1), kg/m 3
ρ V1i goes out--the interior gas vent average density of spiral pipe (1), kg/m 3
Spiral pipe 1 interior liquid average density
ρ L1i is flat=( ρ L1i advances- ρ L1i goes out)/ln ( ρ L1i advances/ ρ L1i goes out)
In the formula:
ρ L1i is flat--the interior liquid average density of spiral pipe (1), kg/m 3
ρ L1i advances--the interior liquid-inlet average density of spiral pipe (1), kg/m 3
ρ L1i goes out--the interior liquid outlet average density of spiral pipe (1), kg/m 3
Average density during gas-liquid two-phase
ρ 1i is flat= R L1i ρ L1i is flat+ (1- R L1i) ρ V1i is flat
In the formula:
R L1i--the liquid phase mass fraction;
The interior gas average viscosity of spiral pipe (1)
μ V1i is flat=( μ V1i advances- μ V1i goes out)/ln ( μ V1i advances/ μ V1i goes out)
In the formula:
μ V1i is flat--the interior liquid average viscosity of spiral pipe (1), Pa.s;
μ V1i advances--the interior liquid-inlet average viscosity of spiral pipe (1), Pa.s;
μ V1i goes out--the interior liquid outlet average viscosity of spiral pipe (1), Pa.s;
The interior liquid average viscosity of spiral pipe (1)
μ L1i is flat=( μ Lo advances- μ Vo goes out)/ln ( μ Lo advances/ μ Lo goes out)
In the formula:
μ L1i is flat--the interior liquid average viscosity of spiral pipe (1), Pa.s;
μ L1i advances--the interior liquid-inlet average viscosity of spiral pipe (1), Pa.s;
μ L1i goes out--the interior liquid outlet average viscosity of spiral pipe (1), Pa.s;
Average viscosity during gas-liquid two-phase
μ 1i is flat= R L1i μ L1i is flat+ (1- R L1i) μ V1i is flat
In the formula:
R L1i--the liquid phase mass fraction;
The interior gas Reynolds number of spiral pipe (2)
Re V2i= v V2i ρ V2i is flat d 2i / μ V2i is flat
In the formula:
Re V2i---spiral pipe (2) gas Reynolds number;
v V2i---spiral pipe (2) gas flow rate, m/s;
ρ V2i is flat---spiral pipe (2) gas average density, kg/m 3
d 2i---spiral pipe (2) internal diameter, m;
μ V2i is flat---spiral pipe (2) gas average viscosity coefficient, Pa.s;
The interior liquid Reynolds number of spiral pipe (2)
Re L2i= v L2i ρ L2i is flat d 2i / μ L2i is flat
In the formula:
Re L2i---the interior liquid Reynolds number of spiral pipe (2);
v L2i---the interior flow rate of liquid of spiral pipe (2), m/s;
ρ L2i is flat---the interior liquid average density of spiral pipe (2), kg/m 3
d 2i---spiral pipe (2) internal diameter, m;
μ L2i is flat---the interior liquid average viscosity of spiral pipe (2) coefficient, Pa.s;
Spiral pipe 2 interior gas-liquid two-phase average Reynolds numbdrs
Re 2 i are flat=Re L2 i R L2i+ (1- R L2i) Re V2i
In the formula:
Re 2i is flat---the average Reynolds numbdr during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas average density of spiral pipe (2)
ρ V2i is flat=( ρ V2i advances- ρ V2i goes out)/ln ( ρ V2i advances/ ρ V2i goes out)
In the formula:
ρ V2i is flat--the interior gas average density of spiral pipe (2), kg/m 3
ρ V2i advances--the interior gas feed average density of spiral pipe (2), kg/m 3
ρ V2i goes out--the interior gas vent average density of spiral pipe (2), kg/m 3
Spiral pipe 2 interior liquid average densities
ρ L2i is flat=( ρ L2i advances- ρ L2i goes out)/ln ( ρ L2i advances/ ρ L2i goes out)
In the formula:
ρ L2i is flat--the interior liquid average density of spiral pipe (2), kg/m 3
ρ L2i advances--the interior liquid-inlet average density of spiral pipe (2), kg/m 3
ρ L2i goes out--the interior liquid outlet average density of spiral pipe (2), kg/m 3
The interior gas-liquid two-phase density of spiral pipe (2)
ρ 2 i are flat= ρ L2 i R L2i+ (1- R L2i) ρ V2i
In the formula:
ρ 2i is flat---the average density during gas-liquid two-phase, kg/m 3
R L2i---the liquid phase massfraction;
The interior gas average viscosity of spiral pipe (2)
μ V2i is flat=( μ V2i advances- μ V2i goes out)/ln ( μ V2i advances/ μ V2i goes out)
In the formula:
μ V2i is flat--the interior liquid average viscosity of spiral pipe (2), Pa.s;
μ V2i advances--the interior liquid-inlet average viscosity of spiral pipe (2), Pa.s;
μ V2i goes out--the interior liquid outlet average viscosity of spiral pipe (2), Pa.s;
The interior liquid average viscosity of spiral pipe (2)
μ L2i is flat=( μ L2i advances- μ V2i goes out)/ln ( μ L2i advances/ μ L2i goes out)
In the formula:
μ L2i is flat--the interior liquid average viscosity of spiral pipe (2), Pa.s;
μ L2i advances--the interior liquid-inlet average viscosity of spiral pipe (2), Pa.s;
μ L2i goes out--the interior liquid outlet average viscosity of spiral pipe (2), Pa.s;
Average viscosity during the interior gas-liquid two-phase of spiral pipe (2)
μ 2i is flat= R L2i μ L2i is flat+ (1- R L2i) μ V2i is flat
In the formula:
R L2i--the liquid phase mass fraction;
Adopt the Reynolds number mean value method in the computation process by the gas-liquid phase-splitting and two-phase flow mass fraction calculating after mixing.
4. according to claim 1Described bifilar stream spiral winding tube type heat exchanger Prandtl number computation process is characterized in that:
The shell-side Prandtl number
Pr o= C Po μ O is flat/ λ O is flat
In the formula:
Pr o---the shell-side Prandtl number;
C Po---shell-side specific heat at constant pressure, kJ/kg;
μ O is flat---shell-side average viscosity coefficient, Pa.s;
λ O is flat---mean coefficient of heat conductivity, kg/m 3
Shell-side gas average specific heat at constant pressure
C Pvo is flat=( C Pvo advances- C Pvo goes out)/ln ( C Pvo advances/ C Pvo goes out)
In the formula:
C Pvo is flat--shell-side gas average density, kJ/kg;
C Pvo advances--shell-side gas feed average density, kJ/kg;
C Pvo goes out--shell-side gas vent average density, kJ/kg;
Shell-side liquid average specific heat at constant pressure
C PLo is flat=( C PLo advances- C PLo goes out)/ln ( C PLo advances/ C PLo goes out)
In the formula:
C PLo is flat--shell-side liquid average density, kJ/kg;
C PLo advances--shell-side liquid-inlet average density, kJ/kg;
C PLo goes out--shell-side liquid outlet average density, kJ/kg;
Shell-side gas-liquid two-phase average specific heat at constant pressure
C Po is flat= R Lo C PLo is flat+ (1- R Lo) C Pvo is flat
In the formula:
R Lo--liquid quality divides rate;
Shell-side gas average viscosity
μ Vo is flat=( μ Vo advances- ρ Vo goes out)/ln ( μ Vo advances/ μ Vo goes out)
In the formula:
μ Vo is flat--shell-side gas average viscosity;
μ Vo advances--shell-side gas feed average viscosity, Pa.s;
μ Vo goes out--shell-side gas vent average viscosity, Pa.s;
Shell-side liquid average viscosity
μ Lo is flat=( μ Lo advances- ρ Vo goes out)/ln ( μ Lo advances/ μ Lo goes out)
In the formula:
μ Lo is flat--shell-side liquid average viscosity;
μ Lo advances--shell-side liquid-inlet average viscosity, Pa.s;
μ Lo goes out--shell-side liquid outlet average viscosity, Pa.s;
Shell-side gas-liquid two-phase average viscosity
μ O is flat= R Lo μ Lo is flat+ (1- R Lo) μ Vo is flat
In the formula:
R Lo--liquid quality divides rate;
Shell-side gas mean coefficient of heat conductivity
λ Vo is flat=( λ Vo advances- ρ Vo goes out)/ln ( λ Vo advances/ λ Vo goes out)
In the formula:
λ Vo is flat--shell side gas mean coefficient of heat conductivity, W/ (m.K);
λ Vo advances--shell side gas feed mean coefficient of heat conductivity, W/ (m.K);
λ Vo goes out--shell side gas vent mean coefficient of heat conductivity, W/ (m.K);
Shell-side liquid mean coefficient of heat conductivity
λ Lo is flat=( λ Lo advances- λ Vo goes out)/ln ( λ Lo advances/ λ Lo goes out)
In the formula:
λ Lo is flat--shell-side liquid mean coefficient of heat conductivity, W (m.K);
λ Lo advances--shell-side liquid-inlet mean coefficient of heat conductivity, W/ (m.K);
λ Lo goes out--shell-side liquid outlet mean coefficient of heat conductivity, W/ (m.K);
Shell-side gas-liquid two-phase mean coefficient of heat conductivity
λ O is flat= R Lo λ Lo is flat+ (1- R Lo) λ Vo is flat
In the formula:
R Lo--liquid quality divides rate;
Shell-side gas Prandtl number
Pr Vo= C Pvo μ Vo is flat/ λ Vo is flat
In the formula:
Pr Vo---the shell-side Prandtl number;
C Pvo---shell-side specific heat at constant pressure, kJ/kg;
μ Vo is flat---shell-side average viscosity coefficient, Pa.s;
λ Vo is flat---mean coefficient of heat conductivity, kg/m 3
Shell-side liquid Prandtl number
Pr Lo= C PLo μ Lo is flat/ λ Lo is flat
In the formula:
Pr Lo---the shell-side Prandtl number;
C Lpo---shell-side specific heat at constant pressure, kJ/kg;
μ Lo is flat---shell-side average viscosity coefficient, Pa.s;
λ Lo is flat---mean coefficient of heat conductivity, kg/m 3
The average Prandtl number of shell-side gas-liquid two-phase
Pr O is flat=Pr Lo R Lo+ (1- R Lo) Pr Vo
In the formula:
Pr O is flat---the average Prandtl number during gas-liquid two-phase;
R Lo---the liquid phase massfraction;
The interior gas Prandtl number of spiral pipe (1)
Pr V1i= C Pv1i μ V1i is flat/ λ V1i is flat
In the formula:
Pr V1i---the shell-side Prandtl number;
C Pv1i---shell-side specific heat at constant pressure, kJ/kg;
μ V1i is flat---shell-side average viscosity coefficient, Pa.s;
λ V1i is flat---mean coefficient of heat conductivity, kg/m 3
The interior liquid Prandtl number of spiral pipe (1)
Pr L1i= C PL1i μ L1i is flat/ λ L1i is flat
In the formula:
Pr L1i---the shell-side Prandtl number;
C PL1i---shell-side specific heat at constant pressure, kJ/kg;
μ L1i is flat---shell-side average viscosity coefficient, Pa.s;
λ L1i is flat---mean coefficient of heat conductivity, kg/m 3
The average Prandtl number of the interior gas-liquid two-phase of spiral pipe (1)
Pr 1i is flat=Pr L1i R L1i+ (1- R L1i) Pr V1i
In the formula:
Pr 1i is flat---the average Prandtl number during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas average specific heat at constant pressure of spiral pipe (1)
C Pv1i is flat=( C Pv1i advances- C Pv1i goes out)/ln ( C Pv1i advances/ C Pv1i goes out)
In the formula:
C Pv1i is flat--shell-side gas average density, kJ/kg;
C Pv1i advances--shell-side gas feed average density, kJ/kg;
C Pv1i goes out--shell-side gas vent average density, kJ/kg;
The interior liquid average specific heat at constant pressure of spiral pipe (1)
C PL1i is flat=( C PL1i advances- C PL1i goes out)/ln ( C PL1i advances/ C PL1i goes out)
In the formula:
C PL1i is flat--shell-side liquid average density, kJ/kg;
C PL1i advances--shell-side liquid-inlet average density, kJ/kg;
C PL1i goes out--shell-side liquid outlet average density, kJ/kg;
The interior gas-liquid two-phase average specific heat at constant pressure of spiral pipe (1)
C P1i is flat= C PL1i R L1i+ (1- R L1i) C Pv1i
In the formula:
C P1i is flat---the average specific heat at constant pressure during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas average viscosity of spiral pipe (1)
μ V1i is flat=( μ V1i advances- ρ V1i goes out)/ln ( μ V1i advances/ μ V1i goes out)
In the formula:
μ V1i is flat--shell-side gas average viscosity, Pa.s;
μ V1i advances--shell-side gas feed average viscosity, Pa.s;
μ V1i goes out--shell-side gas vent average viscosity, Pa.s;
The interior liquid average viscosity of spiral pipe (1)
μ L1i is flat=( μ L1i advances- ρ L1i goes out)/ln ( μ L1i advances/ μ L1i goes out)
In the formula:
μ L1i is flat--shell-side liquid average viscosity, Pa.s;
μ L1i advances--shell-side liquid-inlet average viscosity, Pa.s;
μ L1i goes out--shell-side liquid outlet average viscosity, Pa.s;
The interior gas-liquid two-phase average viscosity of spiral pipe (1)
μ 1i is flat= μ L1i R L1i+ (1- R L1i) μ V1i
In the formula:
μ 1i is flat---the average viscosity during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas mean coefficient of heat conductivity of spiral pipe (1)
λ V1i is flat=( λ V1i advances- ρ V1i goes out)/ln ( λ V1i advances/ λ V1i goes out)
In the formula:
λ V1i is flat--shell-side gas mean coefficient of heat conductivity, W/ (m.K);
λ V1i advances--shell-side gas feed mean coefficient of heat conductivity, W/ (m.K);
λ V1i goes out--shell-side gas vent mean coefficient of heat conductivity, W/ (m.K);
The interior liquid mean coefficient of heat conductivity of spiral pipe (1)
λ L1i is flat=( λ L1i advances- λ V1i goes out)/ln ( λ L1i advances/ λ L1i goes out)
In the formula:
λ L1i is flat--shell-side liquid mean coefficient of heat conductivity, W/ (m.K);
λ L1i advances--shell-side liquid-inlet mean coefficient of heat conductivity, W/ (m.K);
λ L1i goes out--shell-side liquid outlet mean coefficient of heat conductivity, W/ (m.K);
The interior gas-liquid two-phase mean coefficient of heat conductivity of spiral pipe (1)
λ 1i is flat= λ L1i R L1i+ (1- R L1i) λ V1i
In the formula:
λ 1i is flat---the mean coefficient of heat conductivity during gas-liquid two-phase;
R L1i---the liquid phase massfraction;
The interior gas Prandtl number of spiral pipe (2)
Pr V2i= C Pv2i μ V2i is flat/ λ V2i is flat
In the formula:
Pr V2i---shell-side Prandtl number, m/s;
C Pv2i---shell-side specific heat at constant pressure, kJ/kg;
μ V2i is flat---shell-side average viscosity coefficient, Pa.s;
λ V2i is flat---mean coefficient of heat conductivity, kg/m 3
The interior liquid Prandtl number of spiral pipe (2)
Pr L2i= C PL2i μ L2i is flat/ λ L2i is flat
In the formula:
Pr L2i---shell-side Prandtl number, m/s;
C PL2i---shell-side specific heat at constant pressure, kJ/kg;
μ L2i is flat---shell-side average viscosity coefficient, Pa.s;
λ L2i is flat---mean coefficient of heat conductivity, kg/m 3
The average Prandtl number of the interior gas-liquid two-phase of spiral pipe (2)
Pr 2i is flat=Pr L2i R L2i+ (1- R L2i) Pr V2i
In the formula:
Pr 2i is flat---the average Prandtl number during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas average specific heat at constant pressure of spiral pipe (2)
C Pv2i is flat=( C Pv2i advances- C Pv2i goes out)/ln ( C Pv2i advances/ C Pv2i goes out)
In the formula:
C Pv2i is flat--the interior gas average density of spiral pipe (2), kJ/kg;
C Pv2i advances--the interior gas feed average density of spiral pipe (2), kJ/kg;
C Pv2i goes out--the interior gas vent average density of spiral pipe (2), kJ/kg;
The interior liquid average specific heat at constant pressure of spiral pipe (2)
C PL2i is flat=( C PL2i advances- C PL2i goes out)/ln ( C PL2i advances/ C PL2i goes out)
In the formula:
C PL2i is flat--the interior liquid average density of spiral pipe (2), kJ/kg;
C PL2i advances--the interior liquid-inlet average density of spiral pipe (2), kJ/kg;
C PL2i goes out--the interior liquid outlet average density of spiral pipe (2), kJ/kg;
The interior gas-liquid two-phase average specific heat at constant pressure of spiral pipe (2)
C P2i is flat= C PL2i R L2i+ (1- R L2i) C Pv2i
In the formula:
C P2i is flat---the average specific heat at constant pressure during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas average viscosity of spiral pipe (2)
μ V2i is flat=( μ V2i advances- ρ V2i goes out)/ln ( μ V2i advances/ μ V2i goes out)
In the formula:
μ V2i is flat--the interior gas average viscosity of spiral pipe (2), Pa.s;
μ V2i advances--the interior gas feed average viscosity of spiral pipe (2), Pa.s;
μ V2i goes out--the interior gas vent average viscosity of spiral pipe (2), Pa.s;
Spiral pipe 2 interior liquid average viscosities
μ L2i is flat=( μ L2i advances- ρ L2i goes out)/ln ( μ L2i advances/ μ L2i goes out)
In the formula:
μ L2i is flat--the interior liquid average viscosity of spiral pipe (2), Pa.s;
μ L2i advances--the interior liquid-inlet average viscosity of spiral pipe (2), Pa.s;
μ L2i goes out--the interior liquid outlet average viscosity of spiral pipe (2), Pa.s;
Spiral pipe 2 interior gas-liquid two-phase average viscosities
μ 2i is flat= μ L2i R L2i+ (1- R L2i) μ V2i
In the formula:
μ 2i is flat---the average viscosity during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
The interior gas mean coefficient of heat conductivity of spiral pipe (2)
λ V2i is flat=( λ V2i advances- ρ V2i goes out)/ln ( λ V2i advances/ λ V2i goes out)
In the formula:
λ V2i is flat--the interior gas mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ V2i advances--the interior gas feed mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ V2i goes out--the interior gas vent mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
The interior liquid mean coefficient of heat conductivity of spiral pipe (2)
λ L2i is flat=( λ L2i advances- λ V2i goes out)/ln ( λ L2i advances/ λ L2i goes out)
In the formula:
λ L2i is flat--the interior liquid mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ L2i advances--the interior liquid-inlet mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
λ L2i goes out--the interior liquid outlet mean coefficient of heat conductivity of spiral pipe (2), W/ (m.K);
The interior gas-liquid two-phase mean coefficient of heat conductivity of spiral pipe (2)
λ 2i is flat= λ L2i R L2i+ (1- R L2i) λ V2i
In the formula:
λ 2i is flat---the mean coefficient of heat conductivity during gas-liquid two-phase;
R L2i---the liquid phase massfraction;
Adopt the Prandtl number mean value method by the gas-liquid phase-splitting and mix the calculating of two-phase flow mass fraction.
5. according to claim 1Described bifilar stream spiral winding tube type heat exchanger overall heat transfer coefficient computation process is characterized in that:
The shell-side convection transfer rate
h o= Q Always/ ( Q s/ h s+ Q r/ h r)
h o--manage outer convection transfer rate, W/ (m 2.K);
Q Always--overall heat exchange amount, W;
Q s--sensible heat, W;
Q r--latent heat, W;
h s--sensible heat convection transfer rate, W/ (m 2.K);
h s?= R Lo h Lo+(1－ R Lo) ?h vo
In the formula:
R Lo--the liquid phase mass fraction;
h Lo--liquid sensible heat convection transfer rate, W/ (m 2.K);
h Vo--gas sensible heat convection transfer rate, W/ (m 2.K);
h r--latent heat convection transfer rate, W/ (m 2.K);
h r=1000～1100, get h r=1100 W/ (m 2.K);
Convection transfer rate
h Vo=0.296 ( λ Vo is flat/ d o) Re Vo is flat 0.609Pr Vo is flat 0.31
h Lo=0.296 ( λ Lo is flat/ d o) Re Lo is flat 0.609Pr Lo is flat 0.31
Calculate the heat exchange amount
Q Always= C Po is flat mΔ t+ m c r o
In the formula:
m--shell-side oeverall quality flow, kg/s;
m c--gas-liquid phase transition quality, kg/s;
r o--shell-side vaporizing liquid latent heat, kJ/kg;
Q s=( C PLo is flat m Lo is flat+ C Pvo is flat m Vo is flat) Δ t
In the formula:
m Lo is flat--shell-side liquid quality flow, kg/s;
m Lo is flat=( m Lo advances- m Lo goes out)/ln ( m Lo advances/ m Lo goes out)
m Vo is flat--shell-side gas mass flow, kg/s;
m Vo is flat=( m Vo advances- m Vo goes out)/ln ( m Vo advances/ m Vo goes out)
Q r =Q Always -Q s
The interior convection transfer rate of spiral pipe (1)
h 1i =0.037 ( λ Li is flat / d Li) (Re Li is flat 0.75) Pr Li is flat 0.42
Spiral pipe 2 interior convection transfer rates
h 2i =0.037 ( λ 2i is flat / d 2i) (Re 2i is flat 0.75) Pr 2i is flat 0.42
The corresponding overall heat transfer coefficient of spiral pipe (1)
K 1=1/(1/ h o+ R 1o+ δ 1 d 1o/( λ 1 d 1m)+ R 1i d 1o/ d 1i+ d 1o/( h 1i d 1i))
In the formula:
K 1--spiral pipe (1) overall heat transfer coefficient, W/ (m 2.K);
h o--shell-side convection transfer rate, W/ (m 2.K);
h 1i--the inboard convection transfer rate of spiral pipe (1), W/ (m 2.K);
R 1o--the dirty coefficient of shell-side outer wall heat, W/ (m 2.K);
R 1i--spiral pipe (1) inner surface heat dirt coefficient, W/ (m 2.K);
λ 1--spiral pipe (1) tube wall heat conduction coefficient, W/ (m.K);
δ 1--spiral pipe (1) pipeline wall thickness, m;
d 1m--spiral pipe (1) pipeline central diameter, m;
d 1o--spiral pipe (1) outer diameter tube, m;
d 1i--spiral pipe (1) internal diameter of the pipeline, m;
The corresponding overall heat transfer coefficient of spiral pipe (2)
K 2=1/(1/ h o+ R o+ δ 2 d o/( λd m)+ R i d o/ d i+ d o/( h 2i d i))
In the formula:
K 2--spiral pipe (2) overall heat transfer coefficient, W/ (m 2.K);
h o--shell-side convection transfer rate, W/ (m 2.K);
h 2i--the inboard convection transfer rate of spiral pipe (2), W/ (m 2.K);
R 2o--the dirty coefficient of shell-side outer wall heat, W/ (m 2.K);
R 2i--spiral pipe (2) inner surface heat dirt coefficient, W/ (m 2.K);
λ 2--spiral pipe (2) tube wall heat conduction coefficient, W/ (m.K);
δ 2--spiral pipe (2) pipeline wall thickness, m;
d 2m--spiral pipe (2) pipeline central diameter, m;
d 2o--spiral pipe (2) outer diameter tube, m;
d 2i--spiral pipe (2) internal diameter of the pipeline, m;
Adopt the convection transfer rate mean value method by gas-liquid phase-splitting and sensible heat and latent heat independence computing method.
6. according to claim 1The effective heat exchange high computational of described bifilar stream spiral winding tube type heat exchanger process is characterized in that:
Total heat transfer in the spiral winding tube type heat exchanger
Q Always =A Always K AlwaysΔ t
In the formula:
Δ t--log-mean temperature difference, K;
Spiral pipe 1 total heat conduction area
A 1 is total =Q 1 is total/ ( K 1 is totalΔ t 1)
Spiral pipe (1) length
l 1 =A 1 is total/ ( n 1 π d 1m)
In the formula:
l 1--spiral pipe (1) length;
Spiral pipe (2) total heat conduction area
A 2 is total =Q 2 is total/ ( K 2 is totalΔ t 2)
Spiral pipe (2) length
l 2 =A 2 is total/ ( n 2 π d 2m)
In the formula:
l 2--spiral pipe (2) length;
l 1= l 2The time computation process finish;
The effective heat exchange height of heat interchanger
l ?3＝ l 1sin a＝ l 2sin a
In the formula:
a--the spiral pipe ascending angle;
l 1≠ l 2The time readjust the spiral pipe inlet velocity v 1i, v 2iThe value size, and recomputate whole flow process, until l 1= l 2The time finish heat-exchanging process computation process.
CN201210303321XA 2012-08-24 2012-08-24 Design calculation method of dual-stream spirally-wound tubular heat exchanger Pending CN102855396A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210303321XA CN102855396A (en) 2012-08-24 2012-08-24 Design calculation method of dual-stream spirally-wound tubular heat exchanger

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210303321XA CN102855396A (en) 2012-08-24 2012-08-24 Design calculation method of dual-stream spirally-wound tubular heat exchanger

Publications (1)

Publication Number Publication Date
CN102855396A true CN102855396A (en) 2013-01-02

Family

ID=47401979

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210303321XA Pending CN102855396A (en) 2012-08-24 2012-08-24 Design calculation method of dual-stream spirally-wound tubular heat exchanger

Country Status (1)

Country Link
CN (1) CN102855396A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104239724A (en) * 2014-09-19 2014-12-24 国家电网公司 Online monitoring and assessment method for heat exchange efficiency of water-cooling system of converter valve
CN106152514A (en) * 2016-09-27 2016-11-23 郑州大学 Novel bifilar stream central flame flame vacuum boiler
CN108595761A (en) * 2018-03-23 2018-09-28 中国科学院理化技术研究所 Formula design of heat exchanger method and device is stuck up based on volume minimization
CN109614712A (en) * 2018-12-12 2019-04-12 东北大学 A kind of spiral winding tube type heat exchanger HEAT EXCHANGE ANALYSIS system
CN110705047A (en) * 2019-09-16 2020-01-17 仲恺农业工程学院 Design calculation method for winding tube type heat exchanger with baffle plate

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101338984A (en) * 2008-08-08 2009-01-07 西安交通大学 Spiral traverse baffle shell type heat exchanger design method
JP2011252626A (en) * 2010-05-31 2011-12-15 Hitachi Cable Ltd Heat transfer tube for double-tube heat exchanger

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101338984A (en) * 2008-08-08 2009-01-07 西安交通大学 Spiral traverse baffle shell type heat exchanger design method
JP2011252626A (en) * 2010-05-31 2011-12-15 Hitachi Cable Ltd Heat transfer tube for double-tube heat exchanger

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
于清野: "缠绕管式换热器计算方法", 《中国优秀硕士学位论文全文数据库（工程科技Ⅰ辑）》 *
彭一川: "气液两相区中平均速度和平均体积分率的全浮力数学模型", 《东北工学院学报》 *
昌锟: "低温翅片管换热器的设计计算研究", 《中国优秀硕士学位论文全文数据库（工程科技Ⅱ辑）》 *
蔡毅: "板式换热器性能的数值模拟和实验研究", 《中国优秀硕士学位论文全文数据库（工程科技Ⅱ辑）》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104239724A (en) * 2014-09-19 2014-12-24 国家电网公司 Online monitoring and assessment method for heat exchange efficiency of water-cooling system of converter valve
CN106152514A (en) * 2016-09-27 2016-11-23 郑州大学 Novel bifilar stream central flame flame vacuum boiler
CN108595761A (en) * 2018-03-23 2018-09-28 中国科学院理化技术研究所 Formula design of heat exchanger method and device is stuck up based on volume minimization
CN109614712A (en) * 2018-12-12 2019-04-12 东北大学 A kind of spiral winding tube type heat exchanger HEAT EXCHANGE ANALYSIS system
CN110705047A (en) * 2019-09-16 2020-01-17 仲恺农业工程学院 Design calculation method for winding tube type heat exchanger with baffle plate
CN110705047B (en) * 2019-09-16 2023-03-07 仲恺农业工程学院 Design calculation method for winding tube type heat exchanger with baffle plate

Similar Documents

Publication Publication Date Title
CN102855396A (en) Design calculation method of dual-stream spirally-wound tubular heat exchanger
CN102261968B (en) Method and device for predicting node temperature of shell and tube heat exchanger
CN105067661B (en) Gas liquid exchanger heat transfer coefficient determining device
CN102810128A (en) Design calculation method for single-strand spiral wound heat exchanger
He et al. Experimental investigation on turbulent heat transfer characteristics of molten salt in a shell-and-tube heat exchanger
CN106932214A (en) A kind of heat exchanger performance and efficiency test platform
CN101414321B (en) Design method for evaporation type cooler/condenser for chemical industry
CN102749104B (en) Method for accurately measuring gas flow and liquid flow in gas and liquid mixed fluid
Zheng et al. Distribution characteristics of gas-liquid mixture in plate-fin heat exchangers under sloshing conditions
Steiner et al. Heat transfer and pressure drop for boiling nitrogen flowing in a horizontal tube: 1. Saturated flow boiling
Zhu et al. Study of convection heat transfer of CO2 at supercritical pressures during cooling in fluted tube-in-tube heat exchangers
Zheng et al. Distribution characteristics of gas-liquid mixture in spiral-wound heat exchanger under sloshing conditions
Ren et al. Effect of frost layer on heat transfer of cryogenic fluid in a finned tube
CN112052550B (en) Design method of shell-side boiling spiral coiled tube heat exchanger
CN106369272B (en) Liquefied natural gas (LNG) gasifier efficient heat-exchanging pipe and open-frame type gasifier
Tong et al. CFD simulation on the operation characteristics of CO2 two-phase thermosyphon loop
CN208350689U (en) A kind of device for evaluating augmentation of heat transfer pipe in-tube condensation characteristic
CN112857108B (en) Design method of wound tube type heat exchanger with shell Cheng Lengning tube pass boiling
Sun et al. Experimental study on improving the performance of dry evaporator with rectifying nozzle type critical distributor
Zhang et al. Experimental and simulation study on the distributor design in plate evaporators
CN108956687A (en) Digitlization heat transfer comprehensive experimental device and method
Zheng et al. The sloshing effects on distribution characteristics of gas–liquid mixture in spiral-wound heat exchanger
Lubis et al. Design of a heat exchanger of three concentric tube layer on contrary flow
CN102999707A (en) Low-temperature methanol-methanol winding tubular heat exchanger design and calculation method
CN110189593B (en) Accurate heat transfer comprehensive experiment device

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20130102