Design and Rating Shell and Tube Heat Exchangers: by John E. Edwards
Design and Rating Shell and Tube Heat Exchangers: by John E. Edwards
Design and Rating Shell and Tube Heat Exchangers: by John E. Edwards
By John E. Edwards
Contents 1.0 2.0 2.1 2.2 3.0 Introduction Fundamentals Basic Theory Heat Transfer Model Selection Design Guidelines
Appendices I II III IV V VI VII VIII IX X XI XII Thermal Design Models Synopsis CC-THERM User Guidelines Thermal Model Selection Shortcut Heat Exchanger Design TEMA Heat Exchanger Layout Designation Typical Overall Heat Transfer Coefficients Typical Resistance Fouling Coefficients LMTD Correction Factor Ft Wolverine Tube General Details Midland Wire Cordage Turbulator Details Tube Dimensional Data Shell Tube Count Data
References 1. 2. 3. 4. 5. 6. 7. 8. [ C20] Hewitt,G.F. et al (1994) Process Heat Transfer, (CRC Press) Perry,R.H. and Green, D. (1984) Perrys Chemical Engineers Handbook, 6th edition (McGraw Hill) Kern,D.Q. (1950) Process Heat Transfer (McGraw Hill) Coulson,J.M. and Richardson,J.F. (1993) Chemical Engineering Vol 1, 4th edition (Pergamon) Skinnet,R.K. (1993) Coulson & Richardsons Chemical Engineering Vol 6, 2nd edition (Pergamon) Chemstations,Inc. CHEMCAD THERM Version 5.1 User Guide Schlunder,E.U. (1993) VDI Heat Atlas (Woodhead Publishing) Seider,D.S., Seader,J.D.Seader and Lewin,R.L. Process Design Principles, (John Wiley & Sons, Inc.) References of this type are to be found in CC-THERM > Help > Appendix
PAGE 2 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
1. 0
Introduction Shell and tube heat exchangers are used extensively throughout the process industry and as such a basic understanding of their design, construction and performance is important to the practising engineer. The objective of this paper is to provide a concise review of the key issues involved in their thermal design without having to refer to the extensive literature available on this topic. The author claims no originality but hopes that the format and contents will provide a comprehensive introduction to the subject and enable the reader to achieve rapid and meaningful results. The optimum thermal design of a shell and tube heat exchanger involves the consideration of many interacting design parameters which can be summarised as follows: Process 1. 2. 3. 4. 5. Mechanical 1. 2. 3. 4. 5. Selection of heat exchanger TEMA layout and number of passes. Specification of tube parameters - size, layout, pitch and material. Setting upper and lower design limits on tube length. Specification of shell side parameters materials, baffle cut, baffle spacing and clearances. Setting upper and lower design limits on shell diameter, baffle cut and baffle spacing. Process fluid assignments to shell side or tube side. Selection of stream temperature specifications. Setting shell side and tube side pressure drop design limits. Setting shell side and tube side velocity limits. Selection of heat transfer models and fouling coefficients for shell side and tube side.
There are several software design and rating packages available, including AspenBJAC, HTFS and CC-THERM, which enable the designer to study the effects of the many interacting design parameters and achieve an optimum thermal design. These packages are supported by extensive component physical property databases and thermodynamic models. It must be stressed that software convergence and optimisation routines will not necessarily achieve a practical and economic design without the designer forcing parameters in an intuitive way. It is recommended that the design be checked by running the model in the rating mode. It is the intention of this paper to provide the basic information and fundamentals in a concise format to achieve this objective. The paper is structured on Chemstations CC-THERM software which enables design and rating to be carried out within a total process model using CHEMCAD steady state modelling software. However the principles involved are applicable to any software design process. In the Attachments a Design Aid is presented which includes key information for data entry and a shortcut calculation method in Excel to allow an independent check to be made on the results from software calculations. Detailed mechanical design and construction involving tube sheet layouts, thicknesses, clearances, tube supports and thermal expansion are not considered but the thermal design must be consistent with the practical requirements. Source references are not indicated in the main text as this paper should be considered as a general guidance note for common applications and is not intended to cover specialist or critical applications. Sources for this paper have been acknowledged where possible. The symbols, where appropriate, are defined in the main text. The equations presented require the use of a consistent set of units unless stated otherwise.
PAGE 3 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
2. 0
Fundamentals The basic layout for a countercurrent shell and tube heat exchanger together with the associated heat curve for a condensing process generated from CHEMCAD are shown below:
PAGE 4 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
2. 1
Basic Theory The fundamental equations for heat transfer across a surface are given by:
Q = U A T lm = w C p ( t ) (t 2 t 1 ) = W C p ( s ) (T1 T 2 )
Where
or
W L
Q heat transferred per unit time (kJ/h, Btu/h) U the overall heat transfer coefficient (kJ/h-m2 oC, Btu/h-ft2-F) 2 2 A heat-transfer area (m , ft ) t lm log mean temperature difference (oC, F) Cp(s) liquid specific heat shell side (kJ/kg-K, Btu/lb-F) Cp(t) liquid specific heat tube side, w tube side flow W shell side flow (kg/h, lb/h)
The log mean temperature difference Tlm (LMTD) for countercurrent flow is given by:
Tlm =
Where T1 T2 t1 t2 inlet shell side fluid temperature outlet shell side fluid temperature inlet tube side temperature outlet tube-side temperature
In design, a correction factor is applied to the LMTD to allow for the departure from true countercurrent flow to determine the true temperature difference. Tm = F t Tm The correction factor is a function of the fluid temperatures and the number of tube and shell passes and is correlated as a function of two dimensionless temperature ratios
(T1 T2 ) R = ( ) t 2 t1
(t 2 t1) S =( ) T1 t1
Kern developed a relationship applicable to any heat exchanger with an even number of passes and generated temperature correction factor plots; plots for other arrangements are available in the TEMA standards. The correction factor Ft for a 1-2 heat exchanger which has 1 shell pass and 2 or more even number of tube passes can be determined from the chart in the Appendix VIII and is given by:
The overall heat transfer coefficient U is the sum of several individual resistances as follows:
1 1 1 1 1 1 = + + + + U h i hfi k x ho hfo
The combined fouling coefficient hf can be defined as follows:
hf =
The individual heat transfer coefficients depend on the nature of the heat transfer process, the stream properties and the heat transfer surface arrangements. The heat exchanger layout depends on the heat transfer area (HTA) so an initial estimate is required based on a trial value of the OHTC. CHEMCAD is used to establish the steady state mass and energy balances across the heat exchanger and typical values of the OHTC are shown in the Attachments. A quick calculation method XLTHERM is also available to assist this procedure. The fouling factors chosen can have a significant effect on the design and again typical values are shown in the Attachments. PAGE 5 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
2.2
Heat Transfer Model Selection The heat transfer model selection is determined by the heat transfer process (sensible, condensing, boiling), the surface geometry (tube-side, shell-side), the flow regime (laminar, turbulent, stratifying, annular), and the surface orientation (vertical, horizontal). A heat transfer model selection flow chart is presented in the Appendix IV to assist in this procedure. This flow chart indicates all the models available in CC-THERM and shows the default selections. A synopsis of the various heat transfer models together with their conditions of application is given in Appendix I. The key features of the models are summarised below: Shellside Film Coefficient Methods for Single Component Condensation in Laminar Flow The Nusselt Method is used for horizontal condensation under stratifying conditions where the liquid film is draining under gravity with minimum influence due to vapour shear. This is the CC-THERM default method. The Eissenberg Method is applicable to condensation over tube banks and considers condensate layer thickening behaviour. This provides the most conservative heat transfer coefficient prediction as compared to the Nusselt and Kern methods for condensation over a single tube. Range of application is for Reynolds Numbers to be in the range 1800 to 2000. The Kern Method Kern adapted the Nusselt equation to allow evaluation of fluid conditions at the film temperature. This method requires the film to be in streamline flow with a Reynolds Numbers range 1800 to 2100. Shellside Film Coefficient Methods for Single Component Condensation in Turbulent Flow The Colburn Method is based on a correlation of industrial data for a wide range of fluids in heat exchangers using standard tube pitch designs. Range of application is for Reynolds Numbers to be in the range 2000 to E06 gives results with a deviation +20% safe. It provides a good method for the verification of computer derived heat transfer coefficients. The McNaught Method takes into account the effects of shear controlled heat transfer and the combination of gravity and shear effects. This is the CC-THERM default method. Tubeside Film Coefficient Methods for Single Component Condensation The Chaddock and Chato adaptation of the Nusselt Method The method is applicable for gravity controlled condensation where the influence of vapour shear is low and we have a liquid film draining under gravity forming a stratified layer at the bottom of the tube The Chemstations Method This is based on Duckler (downflow) and Hewitt (upflow) adaptations to Deissler and von Karman equations. The method is applicable to condensation under shear controlled conditions for vertical and horizontal layouts under laminar or turbulent flow. The influence of gravity is negligible compared to the interfacial shear stress. VDI Film Method The Association of German Engineers (Verein Deutscher Ingenieure, VDI) have developed extensive methods for heat exchanger sizing based on a Heat Atlas method. This method is available as an option in CC-THERM for condensation inside vertical tubes.
PAGE 6 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
Design and Rating of Shell and Tube Heat Exchangers 2.2 Heat Transfer Model Selection
Method for Multi-Component Condensation Silver Bell Ghaly The SBG method is based on the vapor phase condensing / cooling process following the equilibrium integral condensation curve which is met provided the Lewis Number Le, the ratio of Sc to Pr, is close to unity and all the heat released, including that from the liquid phase, passes from the interface to the coolant. Deviations from equilibrium will result in errors in the prediction of vapor temperature. If heat is extracted more rapidly than equilibrium the vapor is super cooled or saturated which can lead to fog formation leading to possible pollution problems. If heat is extracted more slowly than equilibrium the vapor is superheated. Tubeside Film Coefficient Methods for Sensible Heat Transfer in Laminar Flow The Seider Tate Equation is applicable to horizontal and vertical pipes involving organic liquids, aqueous solutions and gases with a maximum deviation 12%. It is not conservative for water. Range of application is for Reynolds Numbers to be in the range 100 to 2100 The VDI-Mean Nusselt Method is applicable to heat transfer behaviour involving tube banks. Correlation constants are available for applications with Reynolds Numbers in the range 10 to 2E06. Tubeside Film Coefficient Methods for Sensible Heat Transfer in Turbulent Flow The Sieder Tate Equation (CC-THERM default) is recommended when heating and cooling liquids involving large temperature differences and when heating gases in horizontal or vertical pipes with a maximum deviation 12%. It is not conservative for water. Application to organic liquids, aqueous solutions and gases with Reynolds Number Re>10000, Prandtl Number 0.7<Pr<700 and L/D>60 (e.g. for L=3 ft, D=0.5in and L>=4ft,D>=0.75), heating or cooling. Colburn Method considers applications with varying heat transfer coefficient (U) by assuming the variation of U to be linear with temperature and by deriving an expression for the true temperature difference accordingly. The Dittus-Boelter Equation is recommended for general use noting the standard deviation 12%. Applicable to both liquids and gases with Reynolds Number Re>10000, Prandtl Number 0.7<Pr<160 and L/D>10 ie suitable for applications with shorter tube lengths. Engineering Sciences Data Unit (ESDU) Method is applicable to both liquids and gases involving Reynolds Number 40000<Re<106 and Prandtl Number 0.3<Pr<300 this method gives more precise calculation. Though not mentioned in the text it is suggested that L/D>60 be used .For Prandtl Numbers <100 the Dittus-Boelter equation is adequate. VDI-Mean Nusselt method determines the average heat transfer coefficient for the whole tube bank, as opposed to a single tube in cross-flow, and has been found to correlate with the maximum velocity between tubes rather than upstream velocity and is of more specific interest to heat exchanger designers.
PAGE 7 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
3. 0
Design Guidelines References: Hewitt et al Process Heat Transfer p267, Kern Process Heat Transfer Chapter 7,p127 and Perry Section 11 p11-0 to p11-19 Definitions Heat exchanger configurations are defined by the numbers and letters established by the Tubular Exchanger Manufacturers Association (TEMA). Refer to Appendix V for full details. For example: A heat exchanger with a single pass shell and multi-pass tube is defined as a 1-2 unit. For a fixed tubesheet exchanger with removable channel and cover, bonnet type rear head, one-pass shell 591mm (231/4in) inside diameter with 4.9m(16ft) tubes is defined SIZE 23-192 TYPE AEL Tube Diameter The most common sizes used are 3/4"od and 1"od Use smallest diameter for greater heat transfer area with a normal minimum of 3/4"od tube due to cleaning considerations and vibration.1/2"od tubes can be used on shorter tube lengths say < 4ft. The wall thickness is defined by the Birmingham wire gage (BWG) details are given in Appendix XI(Kern Table 10) Tube Number and Length Select the number of tubes per tube side pass to give optimum velocity 3-5 ft/s (0.9-1.52 m/s) for liquids and reasonable gas velocities are 50-100 ft/s(15-30 m/s) If the velocity cannot be achieved in a single pass consider increasing the number of passes. Tube length is determined by heat transfer required subject to plant layout and pressure drop constraints. To meet the design pressure drop constraints may require an increase in the number of tubes and/or a reduction in tube length. Long tube lengths with few tubes may give rise to shell side distribution problems. Tube Layout, Pitch and Clearance Definitions and Nomenclature B PT C do D baffle spacing(pitch) tube pitch clearance tube outside diameter shell inside diameter
PT C
Square pitch
PT C
Triangular pitch
PT = do + C
STAGGERED ARRAY INLINE ARRAY FLOW FLOW
Triangular pattern provides a more robust tube sheet construction. Square pattern simplifies cleaning and has a lower shell side pressure drop. Typical dimensional arrangements are shown below, all dimensions in inches. Tube od (in) 5 /8 3 /4 1 1 1 /4 11/2
13 Note 1 For shell 12 pitch(square) /16 15 Note 2 For shell 12 pitch(square) /16
Table above uses minimum pitch 1.25 times tube diameter ie clearance of 0.25 times tube diameter. Smallest pitch in triangular 30 layout for turbulent / laminar flow in clean service. For 90 or 45 layout allow 6.4mm clearance for 3/4 tube for ease of cleaning. PAGE 8 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
Design and Rating of Shell and Tube Heat Exchangers 3. 0 Design Guidelines
Shell Diameter The design process is to fit the number of tubes into a suitable shell to achieve the desired shell side velocity 4ft/s(1.219m/s) subject to pressure drop constraints. Most efficient conditions for heat transfer is to have the maximum number of tubes possible in the shell to maximise turbulence. Preferred tube length to shell diameter ratio is in the range 5 to 10. Tube count data are given in Perry Table 11-3 where the following criteria have been used 1) Tubes have been eliminated to provide entrance area for a nozzle equal to 0.2 times shell diameter 2) Tube layouts are symmetrical about both the horizontal and vertical axes 5 3 3) Distance from tube od to centreline of pass partition 7.9mm( /16 ) for shell id <559mm (22in) and 9.5mm ( /8) for larger shells. Heat Transfer Area Using the maximum number of tubes, subject to adequate provision for inlet nozzle, for a given shell size will ensure optimum shell side heat transfer in minimizing tube bundle bypassing. The heat transfer area required design margin is then achieved by adjusting the tube length subject to economic considerations. On low cost tube materials it may be more economical to use standard lengths and accept the increased design margin. It is a common practice to reduce the number of tubes to below the maximum allowed particularly with expensive tube material. In these situations the mechanical design must ensure suitable provision of rods, bar baffles, spacers, baffles to minimize bypassing and to ensure mechanical strength. Baffle Design - Definitions Shellside cross flow area aS is given by Where D B C PT shell i.d. baffle spacing clearance between tubes tube pitch
as =
DC B PT
Minimum spacing (pitch) Segmental baffles normally should not be closer than 1/5th of shell diameter (ID) or 50.8mm(2in) whichever is greater. Maximum spacing (pitch) Spacing does not normally exceed the shell diameter. Tube support plate spacing determined by mechanical considerations e.g. strength and vibration. Maximum spacing is given by
0.75 B = 74 do
Most failures occur when unsupported tube length greater than 80% TEMA maximum due to designer trying to limit shell side pressure drop. Refer to attachments. Baffle cut Baffle cuts can vary between 15% and 45% and are expressed as ratio of segment opening height to shell inside diameter. The upper limit ensures every pair of baffles will support each tube. Kern shell side pressure drop correlations are based on 25% cut which is standard for liquid on shell side When steam or vapour is on the shell side 33% cut is used Baffle pitch and not the baffle cut determines the effective velocity of the shell side fluid and hence has the greatest influence on shell side pressure drop. Horizontal shell side condensation require segmental baffles with cut to create side to side flow To achieve good vapour distribution the vapour velocity should be as high as possible consistent with satisfying pressure drop constraints and to space the baffles accordingly. Baffle clearances The edge distance between the outer tube limit (OTL) and the baffle diameter has to be sufficient to prevent tube breakthrough due to vibration. For example fixed tube-sheet clearances are shown below. Refer to Perry p11-11 for floating head clearances. Shell inside diameter mm (in) Clearance shell id and OTL mm(in) 7 254(10) to610(24) 11( /16) 1 635(25) 13( /2) PAGE 9 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
Design and Rating of Shell and Tube Heat Exchangers 3. 0 Design Guidelines
Db = do (N t K 1)
1n
Where
do Nt
tube o.d. number of tubes 1 0.319 2.142 2 0.249 2.207 Triangular Pitch pt = 1.25 do 4 6 0.175 0.0743 2.285 2.499 Square Pitch pt = 1.25 do 4 6 0.158 0.0402 2.263 2.617 8 0.0365 2.675
Number Passes K1 n
1 0.215 2.207
2 0.156 2.291
8 0.0331 2.643
Typical fouling coefficients are shown in Appendix VII. It can be shown that the design margin achieved by applying the combined fouling film coefficient is given by:
Af = + UC 1 AC hf
where AC is the clean HTA , Af is the dirty or design HTA and UC is the clean OHTC. Fouling Resistances Inside Outside 0.002 0.001 0.002 0.001 0.002 0.002 0.001 0.001 Corrosion Fouling Heavy corrosion can dramatically reduce the thermal performance of the heat exchanger. Corrosion fouling is dependent on the material of construction selection and it should be possible to eliminate altogether with the right choice. However if economics determine that some corrosion is acceptable and no data is available from past experience an allowance of 1 /16 in (1.59 mm) is commonly applied. Design Margin The design margin to be applied to the design is based on the confidence level the designer has regarding the specific application and the future requirements for multipurpose applications. Design of condensers for multipurpose use, where a wide possible variation in flow conditions can exist, provides a particular problem in this regard. It is standard practice to apply a design margin of 15% to the design (dirty) heat transfer area with the result that this is applied to the design margin resulting from the application of the fouling film coefficients discussed previously giving an added safety factor. Pressure Drop
(8)
Results for Typical Fouling Coefficients (British Units) Fouling Coefficients Clean OHTC Inside Outside Combined 500 1000 333 50 500 1000 333 100 500 500 250 50 1000 1000 500 50
For process design using a simulation the following preliminary conservative estimates are given for pressure drops due to friction. Note an additional pressure change occurs if the exchanger is placed vertically. Initial Process Design Pressure Drop Estimates Process Description Pressure Drop (psi) Liquid streams with no phase change 10 Vapor streams with no phase change 2 Condensing streams 2 Boiling streams 1 PAGE 10 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK Pressure (kPa) 70 14 14 7
www.pidesign.co.uk
APPENDIX I HEAT TRANSFER MODEL SYNOPSIS Shellside Film Coefficient Methods for Single Component Condensation in Laminar Flow Horizontal condenser subcoolers are less adaptable to rigorous calculation but give considerably higher overall clean coefficients than vertical condenser subcoolers which have the advantage of well defined zones. The Nusselt Method (Hewitt et al p590)[C20] The mean heat transfer coefficient for horizontal condensation outside a single tube is given by the relationship developed by Nusselt. This correlation takes no account of the influence of vapour flow which, in addition to the effect of vapour shear, acts to redistribute the condensate liquid within a tube bundle.
3 k L L (L G )g ho = 0.725 L do (Tsat Tw ) 0.25
The Kern Method(Kern p263)[S2] Kern adapted the Nusselt equation to allow evaluation of fluid conditions at the film temperature
3 2 k f f g = 0.943 ho f do t f
0.25
For horizontal tube surfaces from 0 to 180 the above equation can be further developed to give
3 2 k f f g = 0.725 ho f do t f
0.25
McAdam extended the above equation to allow for condensate film and splashing affects where the loading per tube is taken to be inversely proportional to the number tubes to the power of 0.667.
3 2 k f f g = 1.51 ho 2 f
0.33
4W 0.667 L N t f
0.33
This equation requires the film to be in streamline flow corresponding to Reynolds Numbers in range 1800 to 2100 The Eissenberg Method (Hewitt et al p660)[C20] Horizontal shell side condensation involving multiple tubes in the presence of vapour is much more complex than the Nusselt single tube correlation, as the flow of condensate from one tube to another results in the condensate layer thickening on the lower tubes decreasing the heat transfer coefficient. For a bank of n tubes the heat transfer coefficient determined by the Nusselt Method above is modified by the Eissenberg expression given below
0.25 hn = ho 0.6 + 0.42 n
hn = ho n
0.167
The Eissenberg correction is more conservative than that due to Kern with Nusselt method being the most conservative ie the highest film coefficient. Shellside Film Coefficient Methods for Single Component Condensation in Turbulent Flow McNaught Method (Hewitt et al p661)[C21] This method is probably the best available at the moment as it takes into account the effects of shear controlled heat transfer and the combination of gravity and shear effects.
PAGE 11 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
APPENDIX I HEAT TRANSFER MODEL SYNOPSIS Tubeside Film Coefficient Methods for Single Component Condensation Kern Modification of Nusselts equation (Perry 10-21, equation 10-105) Laminar Flow This stratified flow model represents the limiting condition at low condensate and vapor rates Horizontal condensation inside tubes based on do
3 k L L (L G )g ho = 0.815 L do (Tsat TW ) 0.25
3 L k L L (L G )g ho = 0.761 WT L
0.25
Where WT is total vapor condensed in one tube A simplification can be made by setting G = 0 in the above correlations. The Nusselt Method (Hewitt et al p594) Chaddock and Chato adaptation for gravity stratifying flow For horizontal condensation inside tubes there are two extreme cases 1) Gravity controlled where the influence of vapour shear is low and we have a liquid film draining under gravity forming a stratified layer at the bottom of the tube 2)Shear controlled where a uniform annular film is formed. The influence of gravity is negligible compared to the interfacial shear stress. For horizontal condensation under stratifying conditions (case 1) the mean coefficient for the whole circumference is given by
ho =
0 0.72 G.75
3 k L L (L G )g hLg L do (Tsat TW )
0.25
The Chemstations Method (Hewitt et al p580-p589 and Perry 10-21)[C23] Duckler (downflow) and Hewitt (upflow) adaptations to Deissler and von Karman equations For condensation under shear controlled conditions for vertical and horizontal conditions the methods for laminar and turbulent flow uses the following procedure for determining the heat transfer coefficient can be summarised : a) The interfacial shear stress is calculated. b) The condensate flow per unit periphery and the Reynolds Number for the liquid film Ref is calculated. c) Estimate which is a function of Ref and
+
e) Calculate the local liquid film heat transfer coefficient from the following relationship
hi =
CpL (L o )
+
0.5
An alternative and more simple method due to Boyko and Kruzhilin is available but not used in CC-THERM Boyko and Kruzhilin adaptation of the Mikheev correlation Vertical condensation inside tubes Mikheev correlation
hLO = 0.021
L hi = hLO 1 + x 1 G
where x
PAGE 12 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
APPENDIX I HEAT TRANSFER MODEL SYNOPSIS Tubeside Film Coefficient Methods for Single Component Condensation VDI Film Method ( VDI Heat Atlas 1992 pJa6- pJa8) [C24] The Association of German Engineers (Verein Deutscher Ingenieure, VDI ) have developed extensive methods for heat exchanger sizing based on a Heat Atlas method. This method is available as an option in CC-THERM for condensation inside vertical tubes. Method for Multi-Component Condensation Silver Bell Ghaly (SBG) (Hewitt et al p635-p636 ) [C1] [C2] The SBG method is based on the following assumptions Vapor phase condensing / cooling follows the equilibrium integral condensation curve (i.e.,Tv = TE ) This condition is met provided the Lewis Number Le is close to unity, where
Le = Sc Pr
All the heat released, including that from the liquid phase passes from the interface to the coolant The heat transfer dQ in an increment of exchanger area comprises heat extracted due to latent heat dQl and sensible heat in the gas dQG and liquid dQL phases giving
dQG = hG (TE Ti ) dA
We define a parameter Z where
Z =
dQG dA dQG = dQ dA dQ
QT 0
A =
(1 + Z U
i
hG dQ
U (TE TC )
Deviations from equilibrium will result in errors in the prediction of vapor temperature. If heat is extracted more rapidly than equilibrium leads to the vapor temperature being less than TE the vapor is super cooled or saturated which can lead to fog formation leading to possible pollution problems. If heat is extracted more slowly than equilibrium giving a vapor temperature greater than TE the vapor is superheated.
PAGE 13 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
APPENDIX I HEAT TRANSFER MODEL SYNOPSIS Tubeside Film Coefficient Methods for Sensible Heat Transfer in Laminar Flow Seider-Tate Equation (Kern p103) Application 100<Re<2100 in heating or cooling applications and in horizontal / vertical pipes involving organic liquids, aqueous solutions and gases with maximum deviation 12%.It is not conservative for water.
0.33
B W
0.14
Tubeside Film Coefficient Methods for Sensible Heat Transfer in Turbulent Flow Seider-Tate Equation (Perry 10-16)[S1] Sieder-Tate applies a viscosity correction factor when heating/cooling liquids with large temperature differences or when heating gases as heat transfer is reduced (ie TB/TW<1). Correction is not required when cooling gases even with large temperature differences. Application to organic liquids, aqueous solutions and gases with Reynolds Number Re>10000, Prandtl Number 0.7<Pr<700 and L/D>60 (e.g. for L=3 ft, D=0.5in and L>=4ft,D>=0.75), heating or cooling and horizontal / vertical pipes with maximum deviation 12%. It is not conservative for water.
Nu = 0.023 Re Pr
0.8
0.33
B W
0.14
Colburn Method (Hewitt et al p105) [S2] Applying the analogy between heat transfer and friction to the friction factor for turbulent flow gives
f 0 = 0.046 Re
0.2
The Colburn equation for turbulent heat transfer in smooth pipes is derived
0. 8 0.33 0.33 Nu = (f 0 2 )Re Pr = 0.023 Re Pr
Colburn also developed a method (Kern p 94 and Fig17) for applications with varying heat transfer coefficient (U) by assuming the variation of U to be linear with temperature and by deriving an expression for the true temperature difference accordingly. Dittus-Boelter Equation (Hewitt et al p105)[S2] Application to both liquids and gases with Reynolds Number Re>10000, Prandtl Number 0.7<Pr<160 and L/D>10 (ie less stringent than Sieder-Tate above) This is recommended for general use bearing in mind standard deviation error of 13%
0.8 n Nu = 0.023 Re Pr
Engineering Sciences Data Unit (ESDU) Method (Hewitt et al p 105) Application to both liquids and gases with Reynolds Number 40000<Re<10 and Prandtl Number 0.3<Pr<300 this method gives more precise calculation. Though not mentioned in the text it is suggested that L/D>60 be used .For Prandtl Numbers <100 the Dittus-Boelter equation is adequate.
0.795 0.495 exp 0.0225 (ln Pr Nu = 0.0225 Re Pr
6
)2]
PAGE 14 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
APPENDIX I HEAT TRANSFER MODEL SYNOPSIS Tubeside Film Coefficient Methods for Sensible Heat Transfer in Turbulent Flow VDI-Mean Nusselt (Hewitt et al p 73-79)[S19] This method determines the average heat transfer coefficient for the whole tube bank, as opposed to a single tube in cross-flow, and has been found to correlate with the maximum velocity between tubes rather than upstream velocity and is of more specific interest to heat exchanger designers. Most cross-flow tube banks are arranged either in in-line arrays or staggered arrays as shown below
STAGGERED ARRAY INLINE ARRAY FLOW FLOW
is the mean Nusselt Number is the Reynolds Number is based on the maximum flow velocity Reynolds Number is calculated using Vmax formulae given in Hewitt Table 2.4 p76 correlation constants correction factors for surface to bulk physical property variations and for the effect of the number of tube rows in the array respectively where F1 is given by
F1 =
This relationship is valid for Pr < 600 and Re > 10
Pr B Pr W
0.26
Where the number of cross-flow tube rows nr >10 F2 1 and for nr = 4 F2 0.9 Values of a and m correlation constants for
p1 D = 1.2 to 4
and
P 2 D 1.15
Table 2.4 p76 for further details re tube bank layouts. In-Line Banks a 0.742 0.211 0.116 m 0.431 0.651 0.700 a 1.309 0.273 0.124 Staggered Banks m 0.360 0.635 0.700
Reynolds Number 10 to 300 300 to 2.0 E05 2.0 E05 to 2.0 E06.
PAGE 15 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
APPENDIX I HEAT TRANSFER MODEL SYNOPSIS Shellside Film Coefficient Methods for Sensible Heat Transfer in Turbulent Flow Stream Analysis (CC-THERM default) This method balances the pressure drop across the baffles for each of the possible flow paths. These include the spaces between the tube od and the baffle hole, between the shell id and the OTL, shell id and baffle od, pass clearance lanes and across the tube bundle. Bell-Delaware Method (Hewitt et al p 275 to p 277) This method incorporates correction factors for the following elements 1.Leakage through the gaps between the tubes and the baffles and the baffles and the shell. 2.Bypassing of the flow around the gap between the tube bundle and the shell 3.Effect of the baffle configuration recognising that only a fraction of the tubes are in pure cross flow. 4.Effect of adverse temperature gradient on heat transfer in laminar flow (Re < 100) but is considered of doubtful validity. The first step is to calculate the ideal cross flow heat transfer coefficient using the VDI-Mean Nusselt The maximum velocity is calculated using flow area calculations depending on tube layout and pitch, baffle spacing, shell diameter and tube bundle diameter. Correction factors are applied to the calculated heat transfer coefficient for baffle configuration, for leakage related to shell to baffle and tube to baffle, and for bypass in the bundle to shell gap. Kern Method due to Colburn (Kern p137) Based on a correlation of industrial data for hydrocarbons, organic compounds, water and aqueous solutions and gases when the bundle employs baffles having acceptable clearances between baffles/tubes and baffles/shell and tube pitches (in) shown below. Range of application is for Reynolds Number 2000<Re<106 gives results with deviation +20% ie safe Tube od (in) 1 1 1 Square Pitch (in) 1 1 9 1 /16 1 7/8 Triangular Pitch (in) 15 /16 or 1 1 9 1 /16 1 7/8
0.14
B W
PAGE 16 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
APPENDIX I HEAT TRANSFER MODEL SYNOPSIS Shellside Film Coefficient Methods for Single Component Condensation in Laminar Flow Horizontal condenser subcoolers are less adaptable to rigorous calculation but give considerably higher overall clean coefficients than vertical condenser subcoolers which have the advantage of well defined zones. The Nusselt Method (Hewitt et al p590)[C20] The mean heat transfer coefficient for horizontal condensation outside a single tube is given by the relationship developed by Nusselt. This correlation takes no account of the influence of vapor flow which, in addition to the effect of vapor shear, acts to redistribute the condensate liquid within a tube bundle.
3 k L L (L G )g = 0.725 ho L do (Tsat Tw )
0.25
The Kern Method(Kern p263)[S2] Kern adapted the Nusselt equation to allow evaluation of fluid conditions at the film temperature
3 2 k f f g ho = 0.943 f do t f
0.25
For horizontal tube surfaces from 0 to 180 the above equation can be further developed to give
3 2 k f f g ho = 0.725 f do t f
0.25
McAdam extended the above equation to allow for condensate film and splashing affects where the loading per tube is taken to be inversely proportional to the number tubes to the power of 0.667.
3 2 k f f g ho = 1.51 2 f
0.33
4W 0.667 L N t f
0.33
This equation requires the film to be in streamline flow corresponding to Reynolds Numbers in range 1800 to 2100 The Eissenberg Method (Hewitt et al p660)[C20] Horizontal shell side condensation involving multiple tubes in the presence of vapor is much more complex than the Nusselt single tube correlation, as the flow of condensate from one tube to another results in the condensate layer thickening on the lower tubes decreasing the heat transfer coefficient. For a bank of n tubes the heat transfer coefficient determined by the Nusselt Method above is modified by the Eissenberg expression given below
0.25 hn = ho 0.6 + 0.42 n
hn = ho n
0.167
The Eissenberg correction is more conservative than that due to Kern with Nusselt method being the most conservative ie the highest film coefficient. Shellside Film Coefficient Methods for Single Component Condensation in Turbulent Flow McNaught Method (Hewitt et al p661)[C21] This method is probably the best available at the moment as it takes into account the effects of shear controlled heat transfer and the combination of gravity and shear effects.
PAGE 17 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
APPENDIX II CC-THERM USER GUIDELINES Design Optimisation CC-THERM always searches from a small size to a large size which ensures the minimum possible excess area consistent with satisfying the user specified shell side and tube side pressure drop and velocity design constraints. If design is pressure drop or velocity limited leading to an oversized area the user can relax the pressure drop and/or the velocity design constraint and possibly adjust tube pitch or diameter to make the design a heat transfer area limited design. CC-THERM issues a message at the end of its search advising if the design is pressure drop, velocity or area limited to assist in the optimization process. The heat exchanger design can be forced by setting design limits to constrain certain parameters. For example restricting tube length to meet an installation constraint will result in an increase in the number of tubes and hence shell diameter. Standard shell sizes are used so an increase in diameter from 8 to 10 could lead to an oversize of 56% derived from the increase in shell area ratio. To achieve final design optimisation the user should switch to the rating mode and adjust tube length until the desired area safety margin has been achieved. Tube Counts For a selected shell diameter, tube design parameters (diameter, pitch, layout) and clearances there is a limit to the number of tubes that can fit determined by the outer tube limits (OTL). Standard tube count tables are provided in Perry Table 11-3 and CC-THERM will always use these values if standard tube sizes are specified in Imperial units. If the design is based on Metric units the user should ensure a practical design has been achieved in regards to tube counts. The table value can be achieved by entering the Imperial size exactly in Metric e.g. entered as 19.05mm not 19mm. LMTD When running UnitOp HEATEX in CHEMCAD the LMTD is based on the inlet and outlet temperatures. CC-THERM LMTD is based on a zone by zone computation resulting in an overall LMTD being a weighted mean average by zone heat load hence the two values will be different. Heat Exchanger Layout When specifying multiple pass configurations in CHEMCAD UnitOp HEATEX this information is not passed on to CCTHERM; the user needs to re-enter this information. User Specified Components For a new component the designer is normally provided with physical properties at the inlet and outlet conditions only. Pure regression can be carried out using two data points only for viscosity, specific heat and thermal conductivity. Density regressions will sometimes require forcing (set weighting at high value e.g. 106 for a given data point) or to change the library equation in the density parameter to a simpler form e.g. linear between close limits and set the data limit values.
PAGE 18 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK
www.pidesign.co.uk
Nomenclature D shell diameter do tube outside diameter B baffle spacing PT tube pitch C clearance Where PT = do + C
Using CC-THERM specify Type, tube size, tube layout, material Assign fluids to shell or tube side
Shortcut design method Estimate heat transfer coefficients Calculate area required
Specify thermal design criteria Shell and tube side heat transfer model Shell and tube side fouling coefficients
Set tube size, length, layout Calculate number of tubes Estimate shell diameter
Specify key design criteria Shell and tube side dp allowable Shell and tube side velocity allowable
Set baffle cut and spacing Check shell side velocity Check tube side velocity
Decide upper and lower design limits Tube length and shell diameter Baffle cut and spacing
L/D 5 to 10 PT 1. 25 do
Using CC-THERM calculate Number of tubes and shell diameter Confirm optimization process valid
Check CC-THERM results for validity Heat transfer coefficients, design margin Tube, shell and baffle details
Shell dp
Tube dp
Optimize for shell dp specification 1 Shell diameter limit not constraining 2 Adjust baffle cut, spacing, tube pitch
Optimize for design margin required 1 Set shell diameter adjust tube length 2 Set tube length adjust pitch or shell
Optimize for tube dp specification 1 Adjust tube length, diameter, passes 2 Shell diameter limit not constraining
Using CC-THERM optimize design Validate with shortcut techniques Document and plot results
PAGE 19 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
CONDENSATION
TUBESIDE
SHELLSIDE GEOMETRY
MULTI
SINGLE COMPONENT
STREAM ANALYSIS
SHELLSIDE GEOMETRY
TUBESIDE
HORIZONTAL
VERTICAL
ORIENTATION
SIEDER TATE
KERN METHOD
NUSSELT
MCNAUGHT METHOD
COLBURN METHOD
THE EISSENBERG
NUSSELT
DITTUS BOELTER
KERNS METHOD
TABOREK METHOD
ESDU METHOD
GRAVITY STRATIFYING
FLOW
SHEAR ANNULAR
PAGE 21 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
G Split Flow
J Divided Flow
PAGE 22 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
X Cross Flow
U U-Tube Bundle
TEMA CLASS R C B
APPLICATION Severe requirements of petroleum and related process applications Moderate requirements of commercial and general process applications Chemical process service
PAGE 23 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
2000
Outside (Btu/ft hdegF) Cold fluid Water Aqueous solutions(1) Organic solvents Light oils Medium organics (1) Light organics(1) Heavy organics(1) Heavy organics(1) Gases Water (1) Water (1) Water Water(1) Water Water(1) Water Water Brine Brine Brine Water
2000
Aqueous solutions <2.0 cp (1) Aqueous solutions >2.0 cp (1) Organic solvents Light organics/oils Medium organics (1) Heavy organics/oils Gases Heavy oils Gases Steam Hydrocarbon vapors Water Water Water Water Aqueous solutions Light organics Heavy organics
PAGE 24 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
2.00E-03 1.00E-03 2.00E-03 1.00E-03 5.00E-04 1.00E-03 3.00E-03 5.00E-04 1.00E-03 3.00E-03 3.00E-03 2.00E-03 3.00E-03 5.00E-04 1.00E-03 3.00E-03
2.00E-03 5.00E-04 1.00E-03 1.00E-03 5.00E-04 1.00E-03 3.00E-03 5.00E-04 1.00E-03 3.00E-03 2.00E-03 1.00E-03 2.00E-03 5.00E-04 1.00E-03 3.00E-03
2.00E-03 1.00E-03 3.00E-03 2.00E-03 5.00E-04 2.00E-03 5.00E-03 5.00E-04 1.00E-03 5.00E-03 4.00E-03 3.00E-03 4.00E-03 1.00E-03 2.00E-03 5.00E-03
2.00E-03 1.00E-03 2.00E-03 2.00E-03 5.00E-04 2.00E-03 4.00E-03 5.00E-04 1.00E-03 5.00E-03 3.00E-03 2.00E-03 3.00E-03 1.00E-03 2.00E-03 4.00E-03
1.00E-03 2.00E-03 1.00E-03 1.00E-02 1.00E-03 1.00E-02 1.00E-02 5.00E-03 2.00E-03 1.80E-03 5.00E-04 1.00E-03 3.00E-03 3.00E-03 2.00E-03 1.00E-03 2.00E-03 2.00E-03 1.00E-03 2.00E-03 2.00E-03 5.00E-04 1.00E-03 3.00E-03 1.00E-03 2.00E-03 5.00E-03 4.00E-03 1.00E-03
PAGE 25 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
PAGE 26 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
PAGE 27 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
APPENDIX X MIDLAND WIRE CORDAGE TURBULATOR DETAILS The information presented here has been downloaded from www.midlandwirecordage.co.uk/htdivision
To calculate hi for fluid flowing through tubes using HTD removable turbulators proceed as follows: 1. 2. Obtain Reynolds Number in plain tube: Rei = (Gi x Di)/ Select insert to be considered in the design and use the appropriate performance curve to determine the values for heat transfer factor (JH) and friction factor (f) corresponding to the Rei value calculated in step 3 3. 4. Calculate hi = JH x (k/Di) x Pr
1/3
x (/w)
0.14
Calculate pressure drop through tubes from: P = (Z x f x L x Np x Gi2)/(g x x Di x (/w)0.14 Where Z = 9.807 x 10 for SI units, giving P in bar Z = 5.36 x 10-10 for English units, giving P in Ib/in2
-5
5.
Iterate design to optimise the relationship between through and over tubes performance, taking into account any variations of external surface which may be appropriate.
CC-THERM provides the facility for taking into account enhanced performance due to turbulators. Enter data under Tube specification or force the inside film coefficient by entering hi determined in Step 3 above.
PAGE 28 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
PAGE 29 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk
APPENDIX XII SHELL TUBE COUNT DATA th These tables are presented for thermal design guidance only. Perry 7 Edition and onwards have removed this table and show methods by calculation. In any event final layout is subject to detailed mechanical design.
PAGE 30 OF 30 MNL 032A Issued 29 August 08, Prepared by J.E.Edwards of P & I Design Ltd, Teesside, UK www.pidesign.co.uk