The synthesis of cost optimal heat exchanger networks

Heat RecoverySystems & ClIP Vol. I0, No. 4, pp. 301-328, 1990 Printed in Great Britain

0890.4332/90 $3.00+.00 ~ 1990Pergamon Press plc

REVIEW PAPER* THE SYNTHESIS OF COST OPTIMAL HEAT E X C H A N G E R NETWORKS AN INDUSTRIAL REVIEW OF THE STATE OF THE ART T. Gtn~D~P,,SL~ and L. NA~S Norsk Hydro a.s Research Centre, N-3901 Porsrunn, Norway

(Recked 22 September 1987; receivedfor publication 5 October 1987) Almmzt~Tbe ~ m ~ / h e a t exchanger network synthesis (HENS) probkm is very complex and involves combinatorial in~oblems in the '*matching" between hot and cold streams to enhance heat mmve~,

teml~tme dqm~ant physical and transport properties, the choice of flow conflt~aeon and nme~l,

of ~ o n for the heat esdmnge~ the mmbinalion of"hard" and "soft" problem data (some target tmnptnturm must be met, while others may be varied within limiu if this is of advanta~ for the total prcom ~onomy), various kinds of ~nmalnts (forbidden and oompui~ry marcia) and diffamn typm of ~ (liquid, vaponr and mixed phase). Pressure drop limitatima and the cost of ~ are aim important. The d e a n object/ue inc~,_~_~a quantitative part (cost of heat exchange equipment and aternal utili"t~) and a qualitative part (safety, operability, flexibility and controllability).This makes it d i g i t to eslablish a single objective fon~on to evaluate the design. Due to topological effecU (services are added or removed), the investment cost function exhibits discontinuities since there is a unit cost involved in the equipment. Some of the qualitative aspects mentioned above cannot easily be formulated ahead of time. The cost of flexibility can be calculated, but only for Wen situations (networks). The g/obal optimum is thus hard to guarantee and the engineer hag,to resort to simplifi~tiom of the model and mine rules that will lead towards a near optimal solution. Resear~ is progeeainll along three different lines which are the use of thermodym~¢ mathematk~/method| ami tbe me of/mow/edt,¢ hazed ~jstem| for protein de|ign. In order to mire retl life indutlrial probleme, the engineer should take advantage of all these di~plines. However, the skin and experien~ of the en~l~Der him~lf will m of vital importance. ~ T ~ paper lmmmtg the state of the art of HENS methods together with some applications tO previous literature problems. An evaluation of the various methods is performed from an indmtt~l" point of view. There is also a brief digumion of some of the software tools available to solve such problems. The presentation will emph-,i~ on the design of the heat exchanger network itself althoellh interactions with the rest of the process and the utility system inevitably will be diacumed. The aspects of flexibility and operability will also be briefly mentioned. HENS is the must mature f~eldof proum synthesis when it comes to systematic methods. The increase of eneflly prices drain8 t h e 70s and early 80s has bern the major driving force. As things have bma developefl, however, the emphaah has changed from eneqly optimal (minimal) strmmu~ to coat optimal network. The iatmt developed methods can find the ~ trado.off between invetlm~t oust algl opemtin$ cost for any price 0cenm'io Cmdadiag rqliomd factors such as the coolin$ water unnlm.atore etc.) ahead of design. N ~ all probletm have not yet been solved, and important research is still being conducted addrtming all three areas of HENS, which are targeting, synthesis and optimization. Industry spends a nignificant amount of money to carry out energy analysis of new and existing plants, to support aatdemia in their resusr~ and to [email protected] and acquire accurate and efficient computer took. It is the involvement of our company in these areas that has given us the opportunity, on an industrial barns, to review the field. In the past there have been two schools of HENS. One relies on thermodynamic principles and a few heuristic rules, where the ckaigner manually (or interactively if software is available) s y n ~ the network. The other, more automatic, approach relies on mathematical methods like linear and nonlinear proigranuning. The relative meritt of the~e echoois will he discerned with reference to case studies,There will also be a short presentation of the historical development within these schools. INTIIODUCnON The process design exercise involves several interacting activities, and Li~nnboff (1983) illustrated this by the "onion diagram". The heat exchanger network synthesis (HENS) problem has been viewed as a subproblem of process synthesis, but the i n t e r i ~ o m with the process itself and the utility system must be

taken into consideration. Even in the case of a frozen process and utility system, the industrial HENS problem is very complex. The combinatorial problem when matching streams and sequencin8 heat exchang,~ has been reported for the clauical literature problems with less than ten raeams. A typical (a~tstrtalprocess conUfim

*This review paper has been reproduced with kind permission from Computers & Chemical Enginecrh,tg 12, 503-530 (1988).

30--80 strmms that need heatinl or cooling. TbeB streams may c h a n p phase and the phye/cal and transport properties of relevance to the H E N S prob301

302

T. GUNDEP.S~Nand L. Nx~ss

lem will change with temperature. The heat transfer conditions will be very different in a gas-gas exchanger, in a liquid-liquid exchanger or in a steam boiler. There are many choices with respect to flow configuration and materials of construction for the heat exchangers. Due to the interaction between the heat exchanger network and the rest of the process, it is important to be aware of the fact that some stream parameters are "hard" and some are "soft". Some target temperatures must be met from process requirements, while others may be changed within limits; if this is of advantage for the heat recovery and thus the overall process economy. Finally, there will be qualitative aspects of operability, flexibility, safety and controllability. In many cases this means that nonthermodynamic constraints will he imposed on the model in the form of forbidden, restricted o~ required matches. The work of Westbrook (1961) on the use of dynamic programming and Hwa (1965) on the use of separable programming and a superstructure including promising configurations are among the first reported attempts to systematically solve the HENS problem. The real pioneers in this area were Rudd and coworkers at the University of Wisconsin (Lee et al., 1970; Masso and Rudd, 1969; Rudd, 1968) and Hohmann and Lockhart at the University of Southern California (Hohmann, 1971; Hohmann and Lockhart, 1976) around 1970. Since then, about 200 papers have been published, out of which more than the half are from the last 4 yr. The area has been previously reviewed several times, either stand-alone or as a part of overall process synthesis. Hendry et al. (1973) viewed the HENS problem as a homogeneous subproblem in process synthesis, where the main difficulty involved is the combinatorial matching and sequencing problem. Siirola (1974) classified the methods presented so far and proposed new branching rules to reduce the problem size. He was also among the first to point out that methods should relax the constraint that utilities are handled last and thus placed at extreme temperatures. Rathore and Powers (1975) and Nishida et al. (1977) compared their own methods with previous work in a tabular review, and Hlavacek (1975, 1978) has reviewed the area as part of the overall process design activity including both simulation and synthesis of steady state and dynamic processes. The latest complete review is the excellent work of Nishida et al. (1981) where HENS is one of several process synthesis areas. Since then, Hohmann (1984) presented HENS as a multitier problem and reviewed the latest achievements in using thermodynamic methods. Westerherg and Grossmann (1985) gave a nice and pedagogic introduction to overall process synthesis with emphasis on heat exr,hangar networks and application of mathematical methods. Grossmann (1985a) has reviewed the application of mixed integer linear programming (MILP) in various

process design situations, including HENS, refrigeration systems and total process systems. Norsk Hydro a.s has been active in prototyping software (Naess and Gundersen, 1984) and has access to several computer tools and the latest methods by supporting research at UMIST in Manchester, CMU in Pittsburgh and NTH in Trondheim. There has also been brief contact with the group at CALTECH in Los Angeles. This knowledge has been used in process synthesis studies within the company on existing and new processes, an~l motivated the preparation of this review paper. Initially, reducing the energy consumption was the objective of the research and methods available. This resulted in the discovery of important and fundamental concepts like the h'eat recovery pinch as a bottleneck for reduced energy consumption. The latest methods put emphasis orl' the total cost of energy and equipment and aspects of flexibility and operability are addressed during design rather than as an add on activity. Although retrofit projects are harder to find in periods with low energy prices (such as in 1986), methods are now available for new designs that find the correct trade-off between investment and operating cost for any price scenario.

THE HISTORy OF

When trying to solve the HENS problem systematically, the problem was initially transformed to a mathematical model and solved by numerical methods. The complexity of the industrial HENS problem mentioned above made several simplifications necessary to make these mathematical models manageable. First of all, one could question the industrial usefulness of these simplified models. Secondly, even if accepting these models, the size of the problems that could be solved due to the inherent combinatorial nature was less than 10 streams, which in itself is a severe limitation in industrial applications. A brief review of the first attempts to solve the HENS problem is given below to aid the subsequent presentation of newer methods. A more thorough description of these works are given in earlier review papers mentioned above. To meet the requirements of a review paper while remaining useful for the industry, some methods have been included purely for completeness, other methods for emphasis. Assignment problem o f L P

Kesler and Parker (1969) divided each stream into small heat duty elements ("exchangelets") of equal size and posed the matching between hot and cold elements as an assignment problem. This approach has been improved by Kobayashi et al. (1971) who used the heat content diagram to allow for stream splits and cyclic matches. Nishida et al. (1971) introduced matching rules for minimizing total area

303

The synthesis of cost optimal heat exchanger networks and Cena et al. (1977) allowed for constraints and multiple utilities. Another attempt to solve the matching problem by simultaneous methods is the already mentioned work of Hwa (1965) on separable programming of a superstructure. Methods that decide one match at a time have been named sequential methods, and include tree searching and heuristic methods. These strategies may or may not utilize mathematical methods. .Decomposition or tree search Several early methods were published within the decomposition or tree searching group. Lee et al. (1970) introduced a branch and bound method, but the combinatorial problems were still severe; Siirola (1974) introduced new rules for branching. Pho and Lapidus (1973) used partial enumeration and their "synthesis matrix" was used by Kelahan and Gaddy (1977) who used adaptive random search. The disadvantage of this matrix is the exclusion of cyclic matches and stream splits. Greenkorn et al. (1978) relaxed this constraint and introduced a heat availability function (I-IAF, actually a grand composite curve) to assure good initial solutions. Rathore and Powers (1975) used forward branching to avoid the generation and evaluation of infeasible solutions. Grossmann and Sargent (1978a) combined implicit enumeration with heuristic estimates to solve the configuration problem allowing for constraints on the matches. Finally, Menzies and Johnson (1972) used branch and bound for the synthesis of optimal energy recovery networks including mechanical energy. Heuristic methods The heuristic approach was introduced by Masso and Rudd (1969) who weighted a set of rules according to the adaptive learning during design. Ponton and Donaldson (1974) suggested to match the hot stream with the highest supply temperature with the cold stream with the highest target temperature, an approach that was followed by numerous researchers later. Wells and Hodgkinson (1977) presented an extensive list of heurstic rules for general process synthesis considerations, targets and stream matching.

This fact has slowly been recognized in l~rallef with the failure of the pure mathematical methods. Hohmann's feasibility table was the first rigorous way to establish the minimum utility target ahead of design. The famous ( N - 1) rule givesa close target to the minimum number of exchangers in a n~'twork, and even though Linhoff et al. ('1979) later pointed out that there are cases where this target cannot be reached, Hohmann actually discussed the effect of loops and subgraphs. The minimum heat transfer area target was also addressed in a TQ-diagram by the temperature contention concept, which gave guidance to the splitting of streams to reach the target. If cost data are available, it is possible to determine the optimum utility rates and corresponding minimum temperature of approach ahead of design. A very important part of Hohmann's work is the feasible network solution space demonstrated by the area vs energy diagram in Fig. I. The curve connecting the Am and the E m target for various values of ATm divides the "space" in feasible and infeasible solutions. Hohmann pointed out that this line actually defines an effective maximum number of units necessary to reach the area target. The same diagram was also used to discuss the threshold situation. When decreasing the value of ATm the energy consumption decreases while required heat transfer area increases (A/E trade-off). For some stream systems, there is a limiting value of AT_. where energy remains constant, and the curve turns vertical as indicated in Fig. I. Following the line further increases area and reduces the number of heat exchange units (A/U trade-off). Today we know that this trade-off is more complex, since reduced driving forces have the unfortunate side effect of increasing the number of shells. In his design procedure, Hohmann pointed out that networks with more units than the minimum target include heat load loops which represent degrees of freedom that should be used to minimize area. In this optimization, one should keep in mind the possibility of loop breaking that would reduce the number of units by one and thus harvest a step reduction in the overall cost. Finally, sensitivity fac/

Thermodynamic approach Today, it is amazing how little recognition the work of Hohmann (1971) seemed to get among the researchers in the early 1970s. The explanation often given that the work was little published [only as a Ph.D. extract by Hohnumn and Lockhart (1976)] is only part of the story. Attempts to publish the work were actually turned down twice, and the only reason we can see is a very strong confidence in those days that the HENS problem could be automated and solved by mathematical methods. The work of Hohmann had little contribution to the mathematical approach of that time, but the significance to the solution of industrial HENS problems is paramount.

&Tn| n = 5

,--.,,

~/////J

2T-o, ''o

V/,4_u_rc_~/p'~

, , o, . ". ,o .- -

V//////////J//////////_ l"hre~td

[ too,

Energy

Fig. I. Hohmann's feasible network solution space.

304

T. GUND~t.~Nand L. N~ss

tots were calculated and used to replace rigorous simulation in operability considerations. The additivity of simultaneous changes in inlet temperatures and flowrates was also discussed. The heat recovery pinch

Towards the end of the 1970s, the discovery of the heat recovery pinch as a bottleneck for energy savings resulted in increased efforts in academia and industry in developing and applying systematic methods in retrofitand new designs. Umeda et al. (1978, 1979a,b) draw their composite lines in the available ~,'rgy diagram in such a way that the curves touched at one point (temperature). This point, which forms a bottleneck that prevents further heat integration and thus energy savings, was named the pinch, Huang and Elshout (1976) presented similar ideas using the TQ diagram. The fundamental understanding of the heat recovery pinch, including the decomposition effect was presented by Linnhoff et al. (1979). Linnhoff and coworkers further related excess utility comumption to cross pinch heat flow with contributions from process to proceu heat exchange, heating below and cooling above the pinch. Then came the appropriate placement concept for correct integration of turbines, heat pumps and distillation columns into overall processes, and finally the more general plus/minus principle, that suggests which (and how) process modifications should be made to increase the heat recovery. Since then, Linnhoff and coworkers, first within ICI and later at UMIST have developed the pinch concept into a complete methodology named pinch technology, addressing several aspects of process design that will be discussed in a later section. DESIGN OF HEAT EXC~ANGEI

NETWOItKS In this section, the various contributions will be classified in subsections for: (1) important concepts; (2) targets; (3) synthesis methods; (4) optimization; and briefly (5) flexibility. By concepts we mean both fundamental physical insight and representations that enhance the analysis and understanding. Targets have already been discussed and involve estimates of the theoretically best performance of a heat exchanger network. Synthesis methods include the matching of hot and cold streams and the sequencing of the resulting heat exchangers. Optimization involves both topological and parameter improvmnents that reduce the total annual cost. Finally flexibility is the ability of a heat exchanger network to cope with operating conditions different from design case. Important concepts

The temperatum/eathalpy (TQ) diagram hns boon adopted by several remue.hen. What is new in heat exchanger network duign, is the merging of all hot

streams into a composite heat source curve and all cold streams into a composite heat sink curve. The composite curves were first used by Huang and Elshout 0976). A similar concept is the composite line in the available energy diagram applied by Umeda et al. 0978). The composite curve is an important concept in the work of Linnhoff and coworkers. In the "User Guide" by Linnhoff et aL (1982), the grand composite curve (GCC) is introduced and explained on the basis of the heat cascade. The applications involve multiple utility targeting, flue gas optimization, utility pinches, and to some extent, process modifications. Itoh et ai. (1982) at the same time introduced the heat demand and supply (HDS) diagram, which is equivalent to the Grand Composite. The application was to find load and level for the utilities, to observe options for steam generation and how to integrate different processes on a site (relative position of the individual pinches). Two perhaps less known but early variants of the C-CC are the heat surplus diagram introduced by Flower and Linnhoff (1977) to indicate opportunities for using heat pumps or utility generation (waste heat boilers) and the heat availability function (HAF) introduced by Greenkorn et al. (1978). Its application was to graphically find minimum utility targets and the lowest hot utility temperature level. This work is based on fundamental results obtained by Raghavan (1977). As an aid to make decisions about matching hot and cold streams, Kobayashi et al. (1971) and Nishida et al. (1971) introduced the heat content diagram, which is a heat capacity flowrate vs ternperature diagram. The duty of each stream is represented by boxes in the diagram, and one can allow for temperature d ~ [ ~ d ~ t heat capaciti~. By dividing the box for a stream horizontally or vertically, one can represent multiple matches and stream splits. A similar representation is the heat spectrum diagram used by Nishimura (1980). Linnhoff (1979) introduced two concepts that have proven to be of advantage for HENS, the heat cascade diagram and the stream grid. The heat cascade illustrates the basic idea behind the problem table algorithm 0rrA) and gives rise to the GCC. The stream grid is a clever way to represent hot and cold streams, the process and utility pinches, and most important, the heat exchanger network. The representation makes it easy to redraw other sequences, illustrate head load loops and paths, stream splitting and mixing and has thus been adopted by many other researchers. Another network repres:ntation is the match matrix of Pho and Lapidus (1973), with the inherent disadvantage of not allowing for stream splitting and cyclic matches, a severe limitation to industrial applications. The single most important concept for HENS is the heat recovery pinch that was discovered independently by Umeda et al. (1978) and Linnhoffet al. (1979). The pinch point was also mentioned but not explained by

The synthesis of cost optimal heat exchanger networks Eishout and Hohmann (1979). The same applies to Huang and Elshout (1976) who moved their composite curves towards each other ~ they touch at one point, and the hot composite being above the cold composite at all other places. This position was said to define maximum heat that may be recovered with a corresponding infinite area. The heat recovery pinch forms the basis for a complete methodology named pinch technology developed by Linnhoff and coworkers. The recent progress in the use of automatic and mathematical methods (will be discussed in detail later) is caused by taking the pinch decomposition into consideration. The reason why methods like the one presented by Ponton and Donaldson (1974) may perform poorly on some problems but yield very good solutions in other cases, is the fact that some problems are pinched (even down to AT._ : 0), while o t h e r s a r e threshold problems, which means that for approach temperatures below a certain threshold, only one type of utility (hot or cold) is needed, and that the utility consumption is insensitive to the value of ATm, [see Linnhoff et el. (1979)]. Design should always start where the process is meet coustrained, and for threshold problems with heat surplus, this is at the hot end. Closely linked to the heat recovery pinch is the "plus/m~cus'" principle, which was defined by lJnnhoff ~ l~lr~f (1984) and by l.innhoff and Vredeveld (1984) but had been previously used by Umeda et el. (1979a, b). The idea is to modify the process in such a way, that heat is added to the heat sink above pinch and removed from the heat source below pinch. The appropriate placement conczpt uses the + / - principle to give rules for the integration of heat pumps and turbinm (Townsend and Linnhoff 1983a,b; Linnhoff and Townsend, 1982) and distillation columna (Linnhoff et aL, 1983) in overall processes. A tool for energy management, is the heat path diagram proposed by Westerberg (1983). The utility system and various processes operating on a site are drawn and positioned according to their relative temperatures. Each process is divided by its pinch temperature (if not a threshold problem) into a sink and source part. The relative position of these pinches is important for the possible integration between the processes. The idea to embed several promising networks into one large flowshcet, or superstructure, that was used for the first time by Hwa (1965) has been extended to powerful process synthesis tools during the last years. Papoulius and Grossmann (1983) used MILP to modml and solve superstructures for heat exchanger networks as wull as overall processes. The invention of the atveam ~q~,structure by Floudas et al. (1986) made automatic network generation and optimization possible from the heat load sets given by the MILP solution. Wilcox (1985) also contributed to the automation of the network synthesis by a split-mix-bypau technique.

305

Targets

rUpperand lower bo W were introduced to reduce the combinatorial ~ b l e m limiting some of the first methods. These bounds were later developed further into more or less rigorous targets for energy, heat transfer area, number of beat exchange units and finally total annual cost. I. Energy consumption. Minimum utility targets proposed by Rathore and Powers (1975), Nishida et aL (1977) and Gromanann and Sargent (1978a) are correct for threshold problems, but are too optimistic in pinched situations. These targets were basically used as lower bounds to reduce the search , and as such they were reasonable. Correct energy target has been established for the unconstrained case by several researchers. Hohmann (1971) used the feasibility table, Linnhoff and Flower (1978a) developed the PTA, Umeda et aL (1978) read the utility targets off the TQ diagram and a graphical procedure was also proposed by Greenkorn et al. (1978), who used the HAF (i.e. the grand composite). In constrained situations (forbidden, restricted or required matches) rigorous targets were formulated as an LP transportation problem by Cerda et al. (1983) and as an LP tramshipment problem with reduced size by Papoulias and Grossmann (1983). Recently, Viswanathan and Evans (1987) presented the &,_~_!stream approach to reduce the energy penalty often related to constraints. A hot stream may be heated (then it becomes a cold stream by definition) by the hot part of a forbidden pair, and then substitute the forbidden hot stream to heat up the cold part of the forbidden pair. The problem is formulated as a modified transshipment model and then solved as a network flow problem by the out-of-kilter algorithm. The resulting energy target will then be in the region between the constrained and the unconstrained target, depending on the specific stream data. Both the PTA and LP models may allow for individual stream contributions AT~ to the approach temperature to reflect differences in the heat transfer conditions. Various LP models suggested by Saboo et al. (1986a,b) allow for individual stream contributions, and one may also obtain an energy target for a given structure or a specified total area, a feature that is interesting in retrofit cases. Finally, Doldan et al. 0985), when retrofitting a process with a given utility system for heat and power introduced a secondary bottleneck often encountered in energy studies. This means that there is tittle incentive for reducing the hot utility consumption below the amount available from back pressure turbines, unless the power system is modified. A small nonlinear program (NLP) is used to obtain retrofit energy target. 2. Heat transfer area. When it comes to heat transfer area, one should make the distinction betwcen minimum area procedures which include stra-

306

T.

GUND~

tegies for the network development that minimizes total area, and algorithms for calculating a minimum area target ahead of design. Network generating procedures have been presented by Hohmann (1971) (TQ diagram with temperature and stream contention) and Nishida et al. (1971) who addressed the minimum area task in the heat content diagram matching streams in the order of decreasing inlet temperatures. Nishida et al. (1977) refined this approach by including the utilities in the minimum area algorithm. In this way, utilities are not necessarily placed at extreme temperatures, but located in an "optimal" way in the network. Finally, Umeda et al. (1978) discussed how differences in heat transfer coefficients and unit costs for area could be handled by changing temperatures in the splitting and mixing points in the network. Nishida et al. (1981) in their review paper summarized the work so far on minimum heat transfer area. When all heat transfer coefficients are equal and vertical, counter-current heat transfer is assumed, the following equation for calculating the area target was presented: A~,.

=

U.(Th

-

T,)

(1)

With piecewise linear composite curves as indicated in Fig. 2, this equation can be transformed into a summation: AT~j Townsend and Linnhoff (1984) introduced a new target for area taking into account the individual film transfer coefficients for the streams. The assumption of an average overall U is too crude to be used in cost estimates ahead of design, and equation (3) is thus an important contribution to both area and cost targeting:

zl- q

(3)

The new area target is exact only if all film transfer coefficients are equal and that heat exchange is strictly vertical, but gives good approximations in other situations. The vertical requirement means that all streams, k, in each interval, j in Fig. 2 have to be split and matched according to the total heat capacity flowrate ratio given by the composite curves. The Am network may therefore involve a large number of splits and heat exchangers, and has been named spaghetti design by Litmhoff and o o w o r k e r s . A rigorous target for heat transfer area, even in the case of different film transfer coefficients, can be found by the LP transportation model put forward by Sabot et aL (1986b). Although vertical heat transfer usually leads to minimum total area, there are cases of significant differences in film transfer cotCgg~ents where one could gain from deliberately

and L. NAEss T

Hot streams "qk

I

i+..+t~1~+,.,~

'A

I.]' I

I

I

If,\+

~I !

I

l

--

+q"

-, i\Sec,ionj

;rll

....

,..

s++.,,,,,,o,,.,o.o,,

I-I

Oj

Fig. 2. Enthalpy intervals for area calculation. criss crossing as explained by Abroad (1985)and Tjne

and Linnhoff (1986) and illustrated in Fig. 3. The approach taken by Sabot et aL (1986b) is to start with temperature intervals according to the "kinks" on the composite curves and Madually decrease these intervals until either the area t a r p t obtained by solving the LJ) converges, or the size of the LP becomes prohibitively large. 3. Number oft,nits. The ( N - 1) rule where Nis the total number of hot and cold streams including utilities proposed by Hohmann (1971), has been widely used to set target for the minimum number of units in a heat exchanger network. Boland and Linnhoff (1979)extended this rule by Euler's theorem from graph theory to (N - S + L) to account for the number of subsystems (S) and the number of loops (L). This discussion was repeated by Linnhoff eta/. (1979), who also presented a counter example to Hohinann's claim that an (N - 1) network could be found for any stream system. The observation made by Linnhofl" eta/. is correct, but unfortunately the problem they presented (Fig. 4) is not a true counter example. This was shown by Wood et aL (1985), who discussed the application of splitting, mixing and bypassing, including nonisothermai mixing to reduce the number of units. Figure 5 shows how the posed "counter example" actually meets the units target of two units (no utilities needed). In an earlier work, Linnhoff (1979) also discussed unit reduction by a combination of splitting, mixing and bypassing. Vertical heatf

t_



~

_1 f f

*

f ^r~ ""'~"

<^"~ ....

'. . . .

| L

.r~ t

!

i ~

o

I ,

lit=

Fig. 3. Vertical vs "¢riss-croled" heat transfer (Abroad, 1985).

307

The synthesis of cost optimal heat exchanger networks ISO*

~]

_

SSt"

-

Z40" -



.~p:5.O

:

4SO° ~ I~a='f'~ S41.4"

IlO"

IMp=lO

i

104- , ~ , : . o

li] ~ 4 .

] 0 0 " . mep:Ul

E '4so' •

MO°

,~,

324"

\

i

s i o / ° ° ' r ~ ,,.,:~o 112 m61=O~

Fig. 4. Prolxlcl counter example to the ( N - 1) rule (IAnnhoffet al., 1979). An important aspect of the heat recovery pinch is the decomposition of the problem into separate subnetworks, one completely above and one completely below the pinch temperatures. Linnhoff and Turner (1981) introd~__u:e,3_ the max~nmm energy recovery (MER) target for number of units. U,m.un, by applying the ( N - 1) rule above and below pinch. Grimes et al. (1982) introduced a formula that accounts for streams existing in more than one subnetwork (crossing the pinch), which in essence is the same target as the one presented by Linnhoff and Turner. Cerda and Westerberg (1983) formulated an MILP model for the target which was transformed to an LP transportation problem by relaxation in order to avoid solving the original MILP problem. A similar approach was taken by Papoulias and Grossmann (1983) who formulated and solved an MILP model given the constrained utility target from an LP transshipment model. The advantage of using linear programming is that constraints and multiple utilities can be handled, and that necessary loops in the network to achieve MER as well as the existence of subnetworks is detected ahead of design. Figure 6 shows a four stream problem with no need for utilities, where the minimum number of units target of 3 cannot be reached. Cold stream 3 must be matched with both hot streams 1 and 2, to reach the target temperature of 350°C for the hot streams. This gives two exchangers, but since stream 3 is too small in load to completely cool any of the hot streams, they both need another match, and four units is the true minimum. Hohmann's ( N - 1) rule and the extension by Linnhoff and coworkers (that U,m.Mn is achieved by applying the (N - 1) rule for each subnetwork) consider both temperatures and enthalpies, but enthalpy

_ 4oo. 771

I

s~.~ I sto*~i~,:1o --14 !lp=~

4IILO~ .ll~salt I[S4 tlTg

Fig. 6. A true counter example to the (N - 1) rule (Gro~ smann, 19851)).

only on a "composite" basis in the calculation of Q..~. and Qc.-,- The reason why the MILP model gives a rigorous target is that ~u//vi&~7 stream enthalpies and temperatures are considered. Jones and Rippin 0985) used LP to find all heat load distributions (HLD) with the global (disregarding the pinch) minimum number of units, but still achieve MER. This is possible if streams are split in such a way that hot and cold streams being matched in the pinch region have identical heat capacity flowrates. Some exchangers thus operate across the pinch, but without transferring heat across the pinch which would have resulted in a utility penalty. The price for achieving this is a constant approach temperature equal to &Tea, throughout the exchanger, with a corresponding large heat transfcr area. In case~ where networks (or HLDs) with U,m (global) cannot be found, one may relax the number of units, AT.,. or the MER constraint. Colbert (1982) argued that Um is a bad target for networks with shell and tube heat exchangers. Rather than counting the number of units, one should count the minimum number o f shells. Heggs (1985) pointed out that when using the pinch design method, the exchangers operating in the vicinity of the pinch often would have to be strict counter-current units in order to handle the expected duties with the given inlet and outlet temperatures. A graphical stage by stage procedure for estimating the nuraber of shells from the composite curves was proposed by Shiroko and Umeda (1983). Liu et al. (1985) psoposed to increase the number of shells by one until the Fz-factor for the exchanger reaches an acceptable value. Their exl~.,rience, however, was that the question of single or multi-shell is not important for the network topology. As one of the latest developments in pinch technology (will be discussed in a separate section), Ahnmt [ ] 520*C P40* (1985), Ahmad et al. (1988) also include target for the T T t moPu"O minimum number of shells as part of the a , nual.cost ~ 344"C !04" L mop I4JO target. 4. Total annual cost. Having targets available ahead of design obviously was a major step forward. The .'C mop'5. I -- [.'.JH i~l mop .S 0 targets can be used as a motivation or to give the \4.,- Lm'/ ~o / ,oo. designer confidence that a network is close to "optiNO mal". The design procedure for heat exchanger netFig. 5. Network that achievesthe units target (Wood etaL, works after the discovery of the heat recovery pinch 1985). and by using targets is summarized in Fig. 7. The

308

T. GUNDeX~Nand L. Nx~ss

IIItmem 1". rt m.%

_[ -I I

Tatlatlno QhQeUmln

r,,..h

t

IyntheIl$ POM • ~'M|R Opt|mllatllfl E vaU ( vtA

I_

r

Fig. 7. Design procedure based on tarIets.

~jathem method reed ts the ~ dm~a method OPDM) mmated by L m h e ~ tad tendmmh (1983). Other Wnthem methods aim for area, and in that ~ae the area target must he included in the t a r i n g box. The minimum allowed approach temlmmture b obtained from e:perience, but if the elaborate trading off area, enerl~ and units in the optimization box indicates that the chmen AT.,. was inadequate, an outer loop must he performed. A new minimum approach temperature is then applied at the targeting and the synthesis stage. The pinch point may change in the process, resulting in a possibly very different network. Even with expemn~ in designing a curtaintype of pro~ah the ¢hmen value of AT.~. may not he optimal, since the value of the minimum approach temperature that properly trades off the cost factors involved depends on the following three agpects: (1) relative prices for energy and equipment; (2) shape of the composite curves (parallelor wide open); and (3) good or bad heat tramfer condifiom in the pinch

reIion.

The need for a tool to audst the enl;ineer's "experience" was certainly needed, and the answer is cost targttMg, where a total annual COSt is obtained by nddinI the ¢ontdbutiom from the targets for energy, area and units. By calculating the total annual cost for various value, of AT..,., Ahmnd (1985), Abroad

and Linnhoff (1984, 1986), Unnhoff and Abroad (19g(ahb) can predict an optimal valne of the apimaw,h umlmatmu which will lead to mar-m/nimum coot networks. The aml piece of this "Immle" was the area target from Townsend and Linahoff (1984). As dismIed in a wevious met/on this t a q ~ is not r ~ but ~ mine ~ must he made about the area dimibulion in the ~ thatisnot ~ ~ RI ~ q u ~ wh~ I ~y should go for the risorom target obtained by the LP model put forward by Saboo et o2. (IW~b). 1 ~ does not p~'vent that such a target may find its applica-

tions, but when it comes to cost targeting and design initialization, AT.,, is the crucial parameter. If errors in area target and distribution counteract (experience shows that this in fact is the case) to make a close cost target, we are willing to sacrifice science for industrial applicability. By introducing cost data earlier in the design procedure and applying the total annual cost target to find the optimum ATIa ~ of design (referred to as "supertarfeting"; Linnhoffand Abroad 1986a,b; Naess and Gundersen, 1984), Fig. 8 is obtained which indicates a significant reduction in the design work, especially during optimization. The way targets for energy, area and units play together to obtain the total annual cost target is shown in Fig. 9, where the discontinuities are caused by the kinks on the composite curves. When targeting for total annual cost ahead of design, the investment or capital cost is much harder to predict than the operating cost, since investment is more dependent on the network structure. Energy can he targeted for in a rigorous manner, and experience shows that designs meeting this target are nearoptimal also in total cost. The idea of preoptimization or trading off energy and capital ahead of design has recently been discussed in a few other works, apparently without the necessary rigor in their capital cost targets. Ray and Fonyo (1983) give sparse information about how the capital cost i~ obtained. Li and Motard (1986) proposed an interesting analytical approach to the preoptimization, but ,their method relies on the limiting assumption of an averal~, heat transfer coetf~ent U for all exchangen. A small error is also introduced, since the utility exchangers are not included in the total area calculation. A more promising approach is takeri by Akeelvoli and Lzken (1987) who allows', for temperaturedependent individual film transfer coefficients' for the streams. Vertical heat transfer is assumed, but a heat transfer utilization factor may be specified by the engineer to cope with the ideal counter-cummt

auunlption.

I Itraam data T s T t m.Op

_ _ vJl~ !



"gYps" mrg~I~ =>&TI~ I

Synthesis POMf~MIrR

I

Optimisu zatl~Eb, I Fig.8. New designprocedurewith costmrsets.

The synthesis of cost optimal heat exchanger networks

yf

/

ilm. L ATI

°'"

-AT,

i '''n ATs

ATm~

Fig. 9. The four targets as a function of Tu.

Several contributions have been made to add rigor to the capital cost target. In the previous section on target for minimum number of heat exchange units, the shell vs unit aspect was discussed. Other aspects that might be considered are pumping and piping costs. Finally, some works on exchanger type and materials of construction will be mentioned. Hall (1985) and Ahmad (1985) both discuss how to incorporate exchanger type and materials of construction into the capital cost target. Their suggestion is to adjust the film heat transfer coefficients to account for exotic materials of construction and exchanger configurations different from I-2 shell and tube. Detailed algorithms are given in these works. The most extensive work on capital cost targeting and also on how to design for this capital cost has been given by Ahmad (1985). One of the important conclusions in this work is that if capital cost is wrong, it does not matter as long as the slope (sensitivity) is nearly correct. In that case the predemgn trade-off between capital and energy is correct resulting in a correct & T i and thus a proper initialization of the design. After all, this is more important than the target for capital cost itself. This also means that it is only in ex~e,me cases that exchanger type and materials of construction play a significant role. A rule of thumb is that the manipulated film coefficients will have to differ by an order of magnitude from the original ones. Synthesis In this section we will briefly discuss methods proposed for establi~ing the heat exchanger network. The presentation will attempt to give a complete picture, while the most interesting methods from an industrial point of view will be presented in a later section on Advances in HENS Methods. Some of the early methods were discussed in a previous section and will be omitted here, like the usigmnent methods (Cena et aL, 1977; Kesler and Parker, 1969; Kobayash/ eta[., 19/I; Nishida et ¢d., 1971), the tree mmrching methods (Orem&om et al., 1978; Orommmnn and Sarll~t, 19/8a; Kelahan and Cruddy, 19/7; Lee et aL, 19/0; Mmudes and Johnson, 19/2; Pho and Lapidus, 1973; l~thore and Poweil, 19/$; Sfirola, 1974) and the heuristic methods (Masso and

309

Rudd, 1969; Ponton and Donaldson, 1974; Wells and Hodgkinson, 1977). Most of the early methods used minimum area as the design criterion and the guide for selecting matches. Due to increased energy prices, development of rigorous utility targets and the discovery of the heat recovery pinch, the methods developed towards the end of the 70s and beginning of the 80s were geared towards maximum energy recovery (MER). The heuristic rule accepted and used by most methods today is that networks achieving MER in the fewest number of units are near-optimal, and that area minimization within such a structure is applied as the second step. While total cost considerations earlier w~e handled in an optimization step following the ,ynthmis, the ability to target for capital cost and thus trade-off energy and capital, triggered the development of methods to design for capital as well as energy. Finally, some methods are aimed at reducing the number of units, some methods' are exerjy, based and some methods are in effect multi-objective. 1. Design for area. Minimum area s~jnihesis methods were presented by Nishida et al. (1971), Hohmann (1971), Nishida et,al. (1977) and Umeda et aL (1978). The fast heuristic method suggested by Ponton and Donaldson (1974), also makes sure that driving forces are used properly resulting in low area requirements. In addition, the method will lead to low utility requirements at least for threshold problems. Improvements of Ponton and DonaldsOn's rules have been presented in Zachoval and Konecny (1982), Zachoval eta/. (1984) to handle cross-flow exchangers. Rev and Fonyo (1983) combined the pinch design method with a modification of Ponton and Donaldson's fast method to automate the design. Hama and Matsumura (1982) presented a five-stage procedure for designing heat exchanger networks, where the engineer is allowed to interact between the stages. The synthesis task was performed by dynamic programming to reduce the area to a minimum. As already discussed, area should not be the first objective in the network generation. Later Hama (1984) extended this method to automatic or interactive synthesis, forbidden or imposed matches and minimum allowed exchanger size. A limitation with this method is that utilities are placed last, thus at extreme temperatures. A very interesting approach that simultaneously addresses the task of reducing area and number of units is the "compound" structure (a combination of series and parallel matching) put forward by Rev and Fonyo (1982). Again, since the method does not take the pinch decomposition into account, it will work best in threshold cases, especially those with one large hot (cold) stream and several cold (hot) streams. In Fig. 10, one cold stream is heated by several hot streams. The approach is to split the cold stream in two branches, find one long range stream (LRS),

310

T. GUNDER.~EN and L. NA~SS

[3S] -

HTG

I

. / , o o ; , , , . 7 ,oo 1Io.

(zo]/s~e~s"

I,m~,zo'=

S~ [sT] / o . % 1 T°='" /,/.,.o ~'zso. [ 3 : Heot copocity ftowrOte

Fig. 10. The compound structure proposed by Rev and Fonyo (1982).

Fig. 11. R ~ and Fonyo's solution with a total cost of 19.847 $ yr-i

meaning a stream with a large AT between the supply and target temperature, and classify the remaining hot streams in a low temperature group (LTG) and a high temperature group (HTG). This procedure may be applied recursively to LTG and HTG. 2. Desig~ for units. One of the problems often encountesed when designing for area only is that the number of units tend to be large. This is due to the inherent feature of the HENS problem that ininimnm area requires an effective maximum number of units ("spa~tti" design). Figure 11 shows a five-unit network resulting from Rev and Fonyo's method which is cheaper in total annual cost than a ten-unit solution from Ponton and Donaldson (1974) and a six-unit network from Nishida et al. (1977). The minimum number of units for the problem is four, since it is of the threshold type. When solving this problem with the combined efforts of Papoulias and Grossmann (1983) (MILP model) and Floudas et al. (1986) (stream superstructure and minimum investment cost), a four-unit network (Fig. 12) may be found. With the cost data used by Rev and Fonyo, this network is only marginally cheaper (2.4%) than gev and Fonyo's five-unit network. The advantage of saving one unit has been redL_,O~__by an increase in total heat transfer area. Jones and Pippin (1985) found all heat load distributions satisfying the global minimum number of units without energy relaxation by an LP model. Since Urn.usa normally is greater than Um in pinched cases, heat load loops exist. The traditional approach is to break these loops and pay a penalty in energy (to reestablish A T 1 ) or area (keeping the reduced ATm). An alternative is to split streams passing through the pinch in such a way that hot and cold stream branches matched in an exchanger have identical heat capacity flowrates. Two exchangers on each side of the pinch may be merged into one, thus operating across pinch but without transferring heat across the pinch. This approach may find its applications, but should be used with care. The merged pinch heat exchangers will operate at a constant approach temperature equal to A T u and thus need large heat transfer area. With almost parallel composite curves,

the driving forces are limited anyway, and the approach may find cost effective networks that would be very hard to discover by hand. If, however, the composite curves are opening up away from the pinch temperature, the resulting use of driving forces lead to a total area in excess of what is necessary, and trmy completely reduce the saving of one or more units. Figure 13 shows an MER design developed by the

pinch d r a m method with six units, which is oa¢ t m than the U,~Mam target, since there is an ideal match above pinch between stream 1 and 4. This example was discussed by Jones and Pippin (1985) and the stream data were taken from the Uxr Ouide (l.innhoff et aL, 1982). The dashed line i n d k a ~ the heat recovery pinch. Jones and Pippin demonstrated that a Um solution of five units can be found by a spec~l split stream arrangement as shown in Fig, 14. It is also possible to start with the design in Fig. 13 and use the concspt of loops and paths (Litw,hoffet

aL, 1982) to reduce the number of units without increasing enm'gy commmption, but at the expense of e reduced value for &Tin. The result is the uetwot'k shown in Fig. 15, and it is easy to add up the ne~ssary U. A for the exchangers and compare the two alternative ways to reach U,mo. The three networks have the following total U.A, when utility exchangers are ignored (they carry the same load and operate at identical temperatures for the three cases): PDM solution 6 units (U. A) = 24.75 i Jones and Pippin 5 units (U -.4 ) -- 29.342 Evolved PDM 5 units (U.A)--25.660

[20~ [10.23

28* (?o)

rjLos)

500" 4SSA"

.e/40u-

[S~[]400" 3'3" 48~040'.4" (25.?,] 1~100."ll*,

t

410.__

f,.L,S

[ 3 :Hoot copocity ftowrote Fig. 12. Network from automatic synthesis (MILl)) with a total cost of 19.366 $ yr -t.

The synthesis of cost optimal heat exchanger networks 1"to* [~}

1

m5

135"

~ i

0

~s



~

, .

. l

LO

i

0 90

140

SO* mc e :

~ i

311

30" ml+p: 'l.S

,so.

20 e mep:2.0

,ss.

I

90 30 .L Ilk'lithe ~ . - - a , ~ !,~0 mct:.,.u

--

,~

,o

,io/;

]

.

~

~.

. [] Io. ,,,,,to

.o.-

1;'~'le°

Fig. 13. PDM solution for a 4-stream example from the "User Guide".

Fig. 15. Evolved network from PDM solution by loop breaking.

It should be obvious that the network in Fig. 15 is the best solution, since the area is significantly less than for the network in Fig. 14, and actually only 3.6% more than the six-unit network in Fig. 13. Usually a larger increase in area will be experienced, but in this case the match between stream 1 and 3 is favoured by the merging of the two exchangers between stream 2 and 3. Finally, it should be noted that the composite curves in this case are moderately parallel, which should favour Jones and Rippin's approach. It should be noted that allowing a smaller value of ATm in Fig. 15 adds degrees of freedom to the problem as far as network topology is concerned. It has been shown, by using the MAGNETS computer program (MILP and NLP models), that when the approach temperature is allowed to drop to AT.,, = 5°C, while keeping the energy consumption unchanged from Fig. 13, one of the configurations that can be found is the network of Fig. 15. The same thing would be true for Jones and Rippin's approach. The concept of loops and paths gives, however, the engineer easy guidance in the optimization from the MER design into the non-MER design. A variant of Ponton and Donaldson's rule aimed at units rather than area was presented by Rathore (1982), who ranked the hot and cold streams according to their heat capacity flowrate before matching, in order to reduce the number of heat exchange units. Mocsny and Govind (1984) used a decomposition strategy, looking for approximate subsets in the problem, which will reduce the number of units towards the global minimum. A similar approach was presented by Li and Motard (1986), where AT~r.~ was adjusted until subsystems were obtained resulting in reduced number of units. 3. Design for shells. The difference between exchanger units and actual number of shells has been

reported by Mitson (1984) in some cases to have a significant effect on capital cost. Kardos and Streiow (1983) pointed out that a network solution of the synthesis problem based on purely countercurrent heat exchangers only will remain optimal in practice if each unit can be realized by one exchanger. This is the exception rather than the rule. A two-unit example was presented, where the optimal network used sequential matching. Given the fact that each of these units actually required four shells, a different network structure would have been optimal, indicating that one should consider shells rather than units at the synthesis stage. The topic of shells vs units was treated under targeting and will be further discussed in the section on designing for capital cost. 4. Designfor energy. The first systematic approach to obtaining MER designs is the temperature interval (TI) method of Linnhoff and Flower (1978a). Matching hot and cold streams within temperature intervals, the pinch decomposition is obeyed, since the pinch will coincide with one of the interval temperatures. Liunhoff and Flower (19781)) also proposed an evolutionary development (ED) method as a second step to reduce the number of units. Since stream splitting during the TI method will distort the ED stage, it is recommended not to split streams unless necessary to avoid very complex networks. Even though cyclic matches and stream splitting may substitute each other to obey AT.,, and reduce area requirement, there are cases where streams have to he split to obtain MER. Thus, the nonsplit TI method cannot guarantee MER. A method with some similarities to the TI method was presented by Naka and Takamatsu (1982). Rather than dividing into temperature intervals, they divided the problem at the kinks of the composite curves into enthalpy intervals. An explicit pinch decomposition was made, grouping the intervals above and below pinch. Using enthalpy rather than temperature intervals has the advantage of vertical heat transfer and no heaters and coolers inside the network that must be shifted towards the end or combined to heat exchangers as the case is with the TI/ED method. With the discovery and full understanding of the heat recovery pinch, it became clear that the design exercise should start where the problem is most constrained. For pinched problems, this is at some

~0•

170e U

15o. D

-

--

~" [oJs]

~ = so

GO* mcp:3.0

so..,~,:l.s

|o.s] 140"

~

I"

OOe mgt=4. 0

240

Fig. 14. Jones and Rippin's minimum number of units solution. HR$ IO/+--B

312

T. GUNDV.R~Nand L. NAFAS

intermediate pinch temperature, and for threshold problems it is at the highest (heat surplus) or lowest (heat deficit) temperature. An outline of the PDM is given in Linnhoff and Turner (1981) and in the User Guide (Linnhoff et al., (1982). The comprehensive description is by Linnhoff and Hindmarsh (1983). The PDM gives guidelines for matching streams as well as when and how to split streams in order to achieve MER. In order to reduce the number of units to the fewest possible, each exchanger is made as large as possible,hopefully to satisfyeitherthe hot or the cold stream with respect to enthalpy change. This heuristic rule is referred to as tick-off. The stream termination concept was also used in previous methods, but since these methods did not obey the pinch decomposition, the rule sometimes worked across the pinch, thus resultingin utilitypenalties. Grimes et al. (1982) also decomposed the problem at the pinch in a method intended for application by hand. The method follows some of the ideas presented by Greenkorn et al. (1978), like the "tickoff" heuristic. Reducing units while having a high degree of heat recovery is the main objective of their method, which also can handle multiple utilitiesand constraints on the matches. W h e n designing for energy, the P D M makes sure that no heat is transferred across the pinch and that proper matches are chosen in the thermodynamically most constrained part of the process. However, a more general statement can be made about exchanging heat in various temperature regions of the process. Ahmad (1985) pointed out that the energy target will be violated if matches occur that transfer more heat across an interval temperature than cascaded by the PTA. The heat recovery pinch is the most constrained interval temperature, where no heat is allowed to be transferred. A similar statement was made by Rev and Fonyo 0986), who encountered hidden or pseudo-pinches when designing away from the pinch. The tick-off rule was constrained by the total allowable heat load assigned to a match based on the available heat flow in the cascade. Having established that MER designs are nearoptimal in total cost, especially if the heat recovery level is achieved in the fewest number of units, it is also important to be aware of the fact that a large number of MER designs exist for each stream data set. These designs have different structures and may thus have very different operability and flexibility characteristics. They also have different scope for reducing the number of units by energy (nonMER) or area relaxation (AT < &T,m~). Flower and Linnhoff(1980) approached the HENS problem in a combinatorial manner, but using thermodynamic criteria and topological arguments to reduce the problem size. For classical literature examples like 6SPI, the number of possible networks is approximately 3.6 × I 0~ [(nh' no)!]. There will be 924

U u networks and no more than 6 feasible designs which achieve MER with the minimum number of units. The reduction in number of possible solutions is enormous, and makes it possible to list, evaluate and rank all networks. The TC (thermodynamic combinatorial) method does, however, not allow for stream splits. A similar approach is presented by Govind et al. (1986). An automated approach is also described by Akselvoll and Lzken (1987) which is based on the PDM and presents a predefined number of promising networks. The level of energy recovery is found at the targeting stage, trading-off energy and capital. With the corresponding value of A Tin, networks are created and ranked according to their investment cost. The method allows for temperature dependent physical properties, constraints and different cost equations depending on materials of construction. A partial enumeration of the tree of all possible stream matchings is performed, bounded by the investment cost multiplied with a factor, in order to find a set of near-optimum solutions. The authors also claim that the strategy for stream splitting guarantees minimum number of splits, however, without documentation or presentation of their algorithm. The problem of finding the minimum number of stream splits to obtain MER designs have been solved as an MILP problem by Hanson and Cornish (1984), who included constraints on the matches and the number of split branches. Thermodynamic considerations about the heat recovery pinch is built into the model. Saboo et al. (1986a) allowed the engineer to put constraints on the number and position of stream splits, as well as the number of heat exchangers on a a specific stream. 5. Design for capital cost. A new and important guiding device for the PDM in order to consider area and thus capital is the driving force plot (DFP) presented by Linnhoff and Vredeveld (1984), where each exchanger can be evaluated according to its use of temperature driving forces. In order to reduce total area, it is important that each exchanger uses exactly the amount of driving forces which are available in that temperature region of the process. Compromises must often be made between the DFP and the tick-off rule, since the minimum number of units is not compatible with the minimum total area. The driving force plot may also recommend the splitting of streams in cases to save area, even though not strictly required to achieve MER. The most comprehensive work on designing for capital (as well as energy) was presented by Ahmad 0985). The objective of that research was to add capital cost guidance to the PDM which steers the designer towards M ER in the fewest number of units. With the ability to target for capital cost, two new methods were introduced in order to try to meet this target. In addition to the above mentioned DFP, Ahmad 0985) introduced the remaining problem

The synthesis of cost optimal heat exchanger networks ana/ysb (RPA) for energy and area. Ahmad also applied the energy/area target plot, which is similar to Hohnumn's feasible solution space (Fig. 1), as a map for designs. A network on the target line is optimal with respect to area and energy, but may have a nonoptimal trade-off (wrong AT.,.). The pinch concept requires no heat transfer across the pinch (vertical heat transfer) and limited heat transfer across other interval temperatures (some crisscrossing is allowed). The total heat transfer area will he most sensitive to the heat transfer taking place in the most constrained part of the process, the pinch region. In most cases, minimum area is achieved by vertical heat transfer, which means that the pinch concept and the crin-cross idea actually merge in the pinch region. Away from pinch some amount of nonvertical heat transfer is allowed. The RPA is used to check whether a proposed match will be in conflict with the targeted energy or area. Another important contribution from Ahmad's work is the treatment of near-pinches as regular pinches when the composite curves are narrow. An outline of the des/gu procedure ,is given in Fig. 16. 6. Design f o r exergy. Reducing the exergy loss in a process has always been a sound strategy for design. Liunhoff (1979)concluded, however, that second law analysis in the context of chemical process design is both diffgult to produce and difficult to interpret. The literature on heat exchanger networks show that simplified thermodynam/c methods have replaced entropy or exergy functions. A few approaches will be mentioned that use second law analysis to study some aspects of designing heat exchanger networks. Schmidt (1984) used exergy considerations to deterrains the optimal approach temperature for a single heat exchanger, trading off the cost of exergy loss and the cost of equ/pment. Exergy was evaluated by using the electricity price. The conclusion drawn that was illustrated by a case study is that the optimal approach temperature depends on the actual tern. perature level. The two established techniques of either using a single value of AT.,. for all exchangers or using a single value for the exergetic efficiency were both shown to be wrong. Irazoqui (1986) used entropy to simultaneously study the ability of the process to generate heat and shaft power, giving the heat recovery from hot to cold streams highest priority. Problem independent charts of optimal operating lines were used to derive the network. Liu (1982) and Pchler and Liu (1983) approached the HENS problem as a multiobjective problem (exergy and eapital) and used a multistep evolutionary method. Their thermoeconomic concept made it possible during synthesis to evaluate the designs without cost calculations. The important observation that the loss of available energy actually is independent of the heat exchanger network with fixed utility consumption, was made by Umeda et el. (1977). Exergy loss in heat

313

Optimum eneaDy target I Optimum units target l Optimum area target Optimum tOtaL COSttarget

rl

I.__[

Cemposlte curves

,

,,r, ving

I

~

[

I-"

tl . .

.to,

II

I * rmolmO probtem Ii I analysis II

* driving fame I~Ot III==

=

"ll

!j

DesglnmeM~l

Iv. uttte I> Nonitemflve : n~tch ptecament

Fig. 16. Design for energy and capital (Ahmad, 1985).

exchangers is 5nked to heat transfer with driving forces larger than zero. The trade-off between exergy loss and area was the bas/s for a method of Nishitani et el, (1982). Vinngrad et el. (1983) improved Umeda et al.'s method to handle multiple utilities in an exergy driven method. The matching is performed by graphical tools and dividing the problem into zones by the various utility temperatures. Some of the problems with traditional 2nd law analysis was discussed by Linnhoff and Carpenter (1981), the most important one being to distinguish between inevitable and avoidable losses. In a nitric acid plant, for instance, the converter accoonts for more than half of the exergy loss in the process, but most of it is inevitable. Other problems exist when designing heat exchanger networks. Linnhoff (1986b) pointed out that 2nd law analysis is unit operation oriented, which means that network interactions are lost. Linnboff actually showed that pinch technology (PT) is firmly based on the 2nd law of thermodynamics. The advanmile with PT is that it performs true t y n ~ no base c u e daign is needed as in 2nd law mudym. Linnhoff (1986b) also showed that, the inevitable exergy loss can be found (approx.) from the composite curves with A T . , . - 0 . Alternatively, a practical and economic exergy lois can be a t a ~ from the composite curves and the optimum value of AT., (supertargeting). Finally, Jakob et el. (1986) dltcmsed the practical aspects of 2nd law analysis, and presented a com-

314

T. G t n ~ v ~ and L. N ~

purer program FLEXER to find the location and reason for exergy losses. Automatic synthesis

Since the first formulation of the HENS problem, researchers have tried to automate the design stage. While the majority of these methods used some kind of mathematical programming, a few methods will be mentioned first that try to automate methods based on thermodynamics and heuristic rules on a computer. Rev and Fonyo (1983, 1986) combined the strict rules of the PDM in the pinch region with modifications of the fast method of Ponton and Donaldson (1974) to automate the synthesis of the network away from the pinch. Anothex automated approach is described by AkselvoU and Laken (1987) and by Govind et a/. (1986). When it comes to mathematical methods to automate the synthesis stage, one of the most popular approaches has been to formulate the matching problem as an assignment task of LP. Some of these methods were discussed under the historical section. More recently, the assignment algorithm approach has been used and improved by Jezowski (1982b) who increased the engineer's involvement and applied heuristic rules based on thermodynamics. Jezowski (1982a) also approached the inherent problem with the assignment algorithm that a large number of small exchangers arc created. By introducing preferred matches as those who create neighbours, the evolution phase was simplified by the merging of such units. Later Jezowski and Halat (1986) incorporated the PTA into the optimal ~_~__~_'_'gnmentproblem which was solved repeatedly to avoid excess utility consumption. The solution of the HENS problem by the assignment algorithm has also been reported by Ivakhnenko et al. (1982), Kafarov et al. (1982a,b,c) Kafarov and Schmidt (1982) and Ostrovsky et al. (1985). A number of contributions have been given to the solution of the automatic synthesis problem with more potential for industrial use than the above mentioned mathematical methods. Fundamental to the success of this approach is the use of MILP and taking basic thermodynamic coneepts like the heat recovery pinch into consideration. Papouliu and Grossmann (1983) first solved the minimum utility cost problem (optimizing tbe selection among utilities) by an LP trannhipment model, and then the minimum number of units problem as an MILP model. The resulting beat load sets (or distributions) can be used to derive the network by hand, but for problems of industrial size, this has been reported by Hartmann (1984) to be very difficult even with good knowledge of the actual proccu. The stream superstructure introduced by Fioudas et al. (1986) automates the network generation given the heat load sets, and at the same time minimizes the

investment cost for the given utility consumption and number of units. Cerda and Westerberg (1983) also contributed to the solution of the automated network generation. The minimum utility problem was solved by the LP transportation model presented by Cerda et al. (1983) and the minimum number of units by an MILP model that was transformed and solved repeatedly as an LP problem. The basic difference from the approach taken by Papoulias and Grossmann is that transportation models are used rather than transshipment, the latter being significantly smaller in size. Optimization

The example discussed through Figs 13--15, leads naturally into the optimization section. Optimizing the network involves both topology and parameter changes of the initially synthesized dt~ign in order to minimize the total annual cost. Several methods presented in the 705 have been named evolutionary and may qualify as optimization methods (Linnhoff and Flower, 1978b; McGalliard and Westerberg, 1982; Nishida et aL, 1977; Shah and Westerber8, 1975; Umeda ei a/., 1978). The PDM as presented by Linnhoff and Hind. marsh (1983) describes how MER networks can be evo/ved towards non-MER designs with i ~ energy consumption, but reduced number of exchangers. This is achieved by sequentially breaking loops and reestablishing the driving forces by incmunng utility cousumption through a *'path" in the network. The new concept of supertargeting actually means optimizing the heat recovery level ahead of design by trading off capital cost and operating cost using predesign targets. Naka and Takamatsu (1982) evolved their network established by the EI method by three basic operations (rearrange, merge and shift). Muraki and Hayakawa (1982) introduced evolutionary rules for reducing the number of units from a good initial network obtained by pinch decomposition, the engineer's judgement to split streams, and a new matching rule. Another interesting feature is that available utilities are treated as streams which means that utility placement may be optimized. Shiroko and Umeda (1983) used evolution to simplify complex minimum area networks and finally l~hler and Liu (1984) applied thermoeconomic evolution. By the term thermoeconomic the authors mean that tbermodynamic considerations like available energy is used to replace traditional economic calculations, during network evolution. The target for their evolutionary approach is to simultaneo~ly reduce the number of units and loss of available energy. Evolution was also used by Ghmes et ai. 0982) who first generate a network through a procedure with several similarities to the PDM (targ~ing, pinch decomposition and tick-off ~le to reduce the number of units). Stream splits are introduced only if both Qs.mi~ and U~,.ME~ are "threatened". An interesting

The synthesis of cost optimal heat exchangcr networks feature of their evolutionary approach is that loops are not only broken to reduce units, but sometimes created by adding infinitely small exchangers and then shifting load to march through different topologies. Su and Motard (1984) presented a loop breaking procedure for evolving networks from MER designs established by the TI method towards optimum or near-optimum networks with close to the minimum number of units (without increasing utiSty consumption). A primary level did not allow for stream splits, which was introduced at the second level if U.,~Mn could not be reached in a split-free design. Basically, this is an automated routine which goal is the same as the manual PDM. The networks established with the TI method are complex, and there will be many possible ways to break loops to reach minimum number of units. The approach gives no guidance as to which loop to break apart from the requirement of thermodynamic feasibility. Finally, Saboo et al. (1986a,b) introduced an MILP with constraints on the number and location of stream splits and exchangers. In retrofit cases, the engineer can incrementally add these constraints to the MILP to semi-automatically evolve a network structure more compatible with the existing network. The previous discussion has mainly been concerned with topological optimization, although the concept of loops and paths in the PDM simultaneously addresses network contignration and parameters like duties, areas and temperatures. The optimization of the parameters in a heat exchanger network with a given structure will not be discussed any further, but many of the computer programs mentioned in a later section have facilities for such optimization. Flexibility

Various terms have been used in the literature to describe other features of a heat exchanger network than the steady state operation with given flowrates, supply and target temperatures. Operability normally includes flexibility, controllability, reliability and safety, while flexibility is limited to the steady-state feasible (thermodynamic) operation of the HEN for different operating modes. Another frequently used terra is resilience, with a similar meaning as operability, and which has been divided into dynamic and steady state resilience. In this paper, emphasis will be on flexibility aspects. It is of course desirable that the entire process is as operable or resilient as possible, but Colberg and Morari (1988) stated that in a tightly energyintegrated plant it is especially important that the heat exchanger network is resilient. Nash eta/. (1978) pointed out the feature of HENS that several nearminimum cost topologies exist. When choosing among these, flexibility and operability should be emphasized. Early works addressing flexibility of heat exchanger networks introduced the concept of semi-

315

twity analysis (Hohmann, 1971; and Takamatsu, et a/., 1970). In the latter work, sensitivity tables were used to determine the overdesign factors for each exchanger rather than adding an equal contingency to all units. Grossmarm and Sargent (1978b) introduced probability distribution functions to replace empirical overdesign factors. Finally, Bingzben et al. (1982), applied sensitivity analysis to check the flexibility of the network in retrofit situations. After reporting failures of existing synthesis techniques with respect to handling process uncertainties and variations, Mar~lle et a[. (1982) defined the res0~nt HENS problem. In addition to traditional economic requirements, like maximum energy recovery and minimum number of units, qualitative aspects of flexibility, operability, controllability and safety (named resilience) should be taken into account. Properties of resilience and some advances towards systematic methods were made. Four "worst cases" were designed separately and combined into a design with the flexibility of handling all these situations. I.Analys~. The comer-point theorem was introd~__eed by Saboo and Morari (1984). It claims that a network which is feasible for the corner points (boundary values of the disturbance range for the parameters) will operate for any intermediate situation. Unfortunately, this is only true for a limited class of problems, and Saboo and Morari discussed what has been referred to as the problem of nonconvexity. The reasons for nonconvexity have also been discussed by Colberg and Morari (1988) and by Caland~nil and Stephanopoulos (1986), and involve change of pinchpoint in the process, flowrate variations, stream spfits, temperature dependent heat capacities and phase changes. Saboo et al. (1985) introduced the resilience index (RI) as a measure for the flexibility of heat exchanger networks. This index is obtained through simple physical considerations or a rigorous numerical algorithm, and its application is to detect bottlenecks in existing networks or to guide the selection between alternatives in new design. The key issue is to determine which exchanger that limits, the RI by its load or driving force. In their approach, simplifying assumptions were made that reduce the industrial applicability (no boiling or condeming strea_m~ constant heat capacities, no mA~ flowrate variations and no stream splits). These assumptiops are addressed in current research at CALTECH where some of, the results are in the process of publication (Baboo et al., 1987a,b). Swaney and Grossmann (1985a,b) introduced a different resilience measure called the flexibility index (FI). While the RI is the largest total disturbance load a network is able to handle, the FI determines the largest variation of one parameter independent of other disturbances. 2. Design. The importance of the structure as opposed to the unit operations themselves has been

316

T. G

~

increasingly appreciated during the last few years when it comes to economy, exergy, controllability and also flexibility. Townsend and Morari (1984), among others, have stated that oversizing is an unreliable approach to resilience. The need to consider the aspect of network flexibility at the synthesis stage is thus obvious. Townsend and Morari (1984) introduced a target for the RI to guide the synthesis. They also found that increased network resilience is compatible with minimum heat transfer area for a specific value of ATm~a. Since minimum area requires more than the minimum number of units (actually an effective maximum), there is a trade-off between resilience and number of units. While moderate nonvertical heat transfer in most cases has marginal effect on the total area, it may have a significant effect on the resilience of the network. Colberg et aL (1986) later introduced a similar FI target for heat exchanger networks. Parkinson et aL 0982) applied results from Raghavan (1977), using the heat availability function. Their synthesis procedure for resilient HENs was based on the "worst average" HAF and flexibility tests for energy and temperature. A more mathematical approach was taken by Ikautyman and Cornish (1984) who designed flexible heat exchanger networks based on computer simulation and optimization. Floudas and Groumann (1986) managed to extend the MILP models of Papoulias and Grossmann (1983) to handle the multiperiod situation, which is another aspect of the network flexibility problem. In this case, a predefined number of periods with given duration and process conditions are known, and the "optimal" network is found by taking into consideration equipment cost and utility costs for all periods. The approach is able to handle pinch changes from one period to another, and the automatic network generation is obtained by extending the techniques (NLP formulation based on the stream superstructure) from single-period problems (Floudas and Grossmann, 1987a). Later Fioudas and Grossmann (1987b) showed how the MILP and NLP models from previous works could be effectively coupled with flexibility analysis. Specified uncertainties in flowrates and inlet ternperatures are handled in a two-stage design procedure. A quite differetapproach is taken by Linnhoff and Kotjabasakis 0984) who discussed operability considerations by using the stream grid diagram. Using concepts likedownstream paths (the propagation of a "disturbance" through the heat exchanger network to a "controlled" target temperature) it is possible at the design stage to take operability considerations into account. When breaking a loop in order to reduce the number of units, an additional effect can be that a disturbance path is broken. Another aspect is the placement of utility heaters and coolers and bypasses on heat exchangers. Linnhoff and Smith (1985) stated

and L. NAe&s that the flexibility problem is economically undefined, and that rigorous techniques like the RI concept is of little help in the design process. An extended path analysis is proposed to handle disturbance propagations not only within the heat exchanger network, but through the overall process. Linnhoff and Kotjabasakis (1986) refined these ideas (Kotjabasakis and Linnhoff, 1986; Linnhoffand Kotjabasakis, 1986) into a methodology that makes it possible to simultaneously study the three-way trade-off between energy, capital cost and flexibility. The use of sensitivity tables replaces rigorous simulation to study the passive response of internal network temperatures and target temperatures to changes in supply temperatures, heat capacity flowrates and effective ( U . A ) values for the heat exchangers. A three-stage procedure is suggested in their works. By the concept of downstream paths one can analyze whether a modified flowrate or inlet temperature will affect a specific target temperature. If it will, the stream grid is applied to find possible topological modifications to either break the downstream path or to introduce means to reestablish the target temperature. Finally, the question of how much flexibility and which action to choose can be considerect. Kotjabamkis and Linnhoff (1986) have pointed out that the relative cost-effectiveness of alternative design options depends on the degree of flexibility required. This is illustrated in their work by a graph which indicates the investment of each contingency action as a function of the new target temperature. Backing off on a target temperature to save capital may in some cases also involve choosing another network manipulation. In more complex networks, combined actions can be considered using the sensitivity tables. An interactive approach using object-oriented concepts derived within AI was tak'~'n by Calandranis and Stephanopoulos (1985, 1986). They used a s~nilar approach as the research at UMIST referred to above, with the stream grid as an excellent tool to study structural aspects of operability, and flexibility. Two important concepts in their work are the disturbance propagation network and the freedom of each exchanger, defined by the ability to tolerate a load shift. This concept was used for the first time in HENS by Linnhoff and Flower (1978b). Calandranis and Stephanopoulos (1986) also contributed to the physical understanding of the operability problem, and the reasons for nonconvexities were grouped into intrinsic and pinch associated mechanisms. ADVANCES IN HENS METHODS

In this section, we will present the advances that have been made along three linesof research that we consider to be of importance for the solution of industrial HENS problems today and in the future.

The synthesisof cost optimal heat exchanser networks These lines will be referred to as pi.nch technology (based on thermodynamic c o n c e p ~ ) ~ t h ~ f i c l t l programming (trying to automate subtasks of the design) and knowledge based systems using AI techniques. Pinch technology

Linnhoff and Vredeveld (1984) stated that pinch technology "has come of age", referring to a complete methodology for the design of heat exchanger networks, utility systems, heat and power systems and integration of distillation columns with the background process. They also pointed out the importance of considering process modifications to improve heat recovery before the heat exchanger network problem is addressed. The best introduction to pinch technology is given by Linnhoffet al. (1979), I.,innhoff and Turner (1981), Linnhoff et ai. (1982) and Linnhoff and Hindmarsh (1983). A new view on energy savings in distillation is given by Linnhoff et al. (1983) and the integration and appropriate placement of heat pumps and turbines is given by Townsend and Linnhoff (1983a,b). 1. Applications. There is little doubt that so far, pinch technology has had a lead among the three mentioned research lines when it comes to industrial applications. Linnhoff and Turner (1980) reported applications within ICI and concluded that three beliefs often encountered in industry when it comes to heat integration are quite unfounded: (1) to design better processes, one needs more design time: (2) process integration for energy conservation leads to control problems; and (3) energy conservation always costs money (area vs energy trade-off). Most processes are highly integrated, and pinch technology is only a tool to make sure that this integration is done properly. Dyson and Kenny (1982) reported how the PDM has been used within Davy McKee, with special emphasis on crude preheat trains. Some interesting practical aspects of using such systematic methods are given. Industrial applications are also reported by Boland (1983), Boland and Hindmarsh (1984) and Linnhoff and Vredeveld (1984) together with some considerations on technology transfer within large companies like ICI and Union Carbide. Steinmetz and Chancy (1985) discussed total site energy integration by using pinch technology and Westerberg's heat path diagram. The importance of the relative positions of the single process pinch points for the total plant integration was discussed. Hindmarsh et al. (1985) described how the "Tensa" Group in ICI addresses the integration of diesel engines, turbines and refrigeration systems into the overall process. Heat load and level information from the GCC is used to exploit heat interactions between a process and a power system to maximize cogenerated power. Finally, O'Reilly (1986) reported experiences with the pinch technology by warning against the indis-

317

criminate use of systematic methods. He stated that design conStrmnts vail be encountered and must be properly dealt with. Since 1984, important new developments have taken place within pinch technology. The specific areas are capital cost targeting and design, retrofit targeting and design, flexibility considerations at the design stage and finally total site beat and power integration. 2. Capital. The ability to target and design for capital cost as well as for energy has added a new dimension to pinch technology. Even though MER designs in the fewest number of units by experience are near-optimal in total cost, there are obviously cases where total area can be improved to reduce the investment. In cases where no significant improvement in total annual cost is observed, targeting and design for capital and energy will at least reduce the time consuming optimization f A / E ind U / E trade-offs) and also give the engineer confidence that the design is near-optimal. A spin-off from this rescarcn has been increased understanding of the mechanisms in HEN design. This includes concepts like topology traps and crisscrossing, which will be briefly discussed,in the following. The area target based on individual film beat transfer coefficients from Townsend amt Linnhoff (1984) is fundamental to this progress. Ahmad and Linnhoff(1984) reported how targets can be obtained for total annual cost. Tiffs made it possible.to obtain a near optimal value of &Tin ahead of design, which reduced downstream optimization. Later Ahmad (1985) introduced a new method for targeting the number of shells as opposed to units. The last contribution along this line is the concept of topology traps by Linnhoff and Ahmad (1986a,b). When choosing a certain AT, t, one also indirectly makes topological decisions. In a design optimization there may be idiosyncrasies in the topology that make evolution from an initial network to the global optimum impossible. With price differences or local specialities like cooling water temperature, the optimum AT,,~ will be different, possibly leading to different topologies rather than one "world-wide" optimum design. With the driving force plot and the remaining problem analysis, Ahmad 0985) showed how the PDM can be extended to simultaneously design for energy and capital. While the cross pinch concept relates to energy, criu~rossing is linked to heat transfer area. These two mechanisms converge in the most constrained part of the process. 3. Retrofit. Tjoe and IAnnhoff (1984) defined the retrofit HENS problem. In this case, the optimal minimum approach temperature is a function of the energy/capital trade-off and in addition the existing structure of the network and the expected payback time. They outlined a retrofit procedure where heat

318

T. GUNDF.S~Nand L. N ~ s

exchangers transferring heat across the pinch are removed and that compatibility with the existing network is searched by manipulating heat load loops and paths. Linnhoffand Tjoe (1985) presented a procedure for retrofit targeting based on the idea that new exchangers or shells installed will be utilized with the same area efficiency. They presented a graph relating investment cost and utility savings. Depending on expected payback times or an upper limit on the investment, one may obtain an estimate for AT,,, to be used in the network manipulation together with target for energy reductions. Tjoe and Linnhoff (1986) presented a complete methodology for retrofit targeting and design. An existing network is usually nonoptimai in both structure and area/energy tradeoff. They pointed out that the smallest approach temperature only indicates wrong use of driving forces, and that this value should not be used in the retrofit design procedure. An in depth treatment of this topic, including industrial case studies, is given by Tjoe (1986). Linnhoff and Witherell (1986) reported application of these latest developments when retrofitting parts of an ethylene plant. The importance of correct data extraction was emphasized. 4. Flexibility. With the work of Kotjabasakis and Linnhoff (1986), Linnhoff and Kotjabasakis 0984, 1986), pinch technology is able to address flexibility considerations at the design stage. The concepts of downstream paths and sensitivity tables applied in the stream grid give the engineer a "bird's eye view" when designing cost optimal networks that also are flexible. 5. Total site. Morton and Linnhoff (! 984) discussed total site integration by using the GCC and the relative location of the individual process pinches. Linnhoff (1986a) introduced the balanced composite curves to handle the situation where the use of intermediate level utilities may introduce new pinch points called utility pinches. The paper discussed the process/utility interface, and the relative merits of the CC and the GCC are explained. These principles were further discussed by Linnhoff and Ahmad (1986b). Current research in pinch technology is progressing along several lines and Linnhoff (1985) introduced a variant of the "Rubic Cube" to describe three basic dimensions of the development. From heat exchanger networks, the technology has been extended to handle heat and power networks, total processes and overall plants. After starting with energy, capital has been taken into consideration and by the work of Kotjabasakis and Linnhoff (1986), Linnhoff and Kotjabasakis (1984, 1986), operability is handled at the design stage. The next step will he raw materials utilization. Finally, the last dimension of the picture is the development from new designs to also handling retrofit cases. The techniques have also been briefly applied to batch processes. Linnhoff (1985) gave this overview as a reply to O'Reilly (1985) who suggested

that the pinch concept had been developed beyond its capabilities. Although this section has been devoted to the research at UMIST in Manchester (Linnhoff and coworkers), the contribution of Hohmann and Lockart and the Japanese researchers, first of all Umeda and coworkers, have been of paramount importance in the development of the thermodynamic approach as a whole. Mathematical programming

In this section, we will describe the .latest developments of the various approaches that apply mathematical methods in order to make parts of, or the whole design exercise automatic:, The considerable efforts in this area during the 70s did not result in many appfications in the industry. There has, however, been a considerable progress in the last few years, when some basic thermodyuan~c concepts like the heat recovery pinch have been included in the models. The activity on using various LP models (including MILP) is fundamental in this picture, and important contributions have been given by Westerherg, Grossmann and coworkers at CMU in Pittsburgh and Morari and coworkers at CALTECH Pasadena. In addition, the works of Su and Motard (1984) who automated the PDM, and the minimum number of units approach of Jones and Rippin 0985) should be mentioned. Rigorous targets can be found for the minimum utility cost; allowing for multiple utilities, all kinds of constrained matching and stream individual contributions to the minimum approach temperature. A rigorous alternative to the (N - 1) rule as target for the minimum number of units, again allowing for constraints, can be found by MILP tranuhipment models. A rigorous area target has also been proposed, although one may question its importance. Automatic and semi-automatic network generation have been reported based on the heat load distributions that obey the targets. This allows for area or cost optimization, by finding optimal values of stream splits and bypass streams. The networks in Figs 6 and 12 illustrate the power of MILP models and automatic synthesis. First, the true minimum number of units can be obtained ahead of design, but more important is the automatic discovery and optimization (investment cost) of networks achieving this target. The MILP approach alto has the advantage that the model for the heat exchanger network easily can be incorporated in an MILP model for the overall process. A mathematical approach to the simultaneously optimization and heat integration of chemical processes was presented by Duran and Grossmann (1986), who reported considerable savings when compared with the sequential approach. A new pinch location method was developed to cope

The synthesis of cost optimal heat exchanger networks

319

with the situation of not having p ~ t a b l i s h e d tem- algorithm to solve the HENS problem Grimes et aL, perature intervals. This is the case because the tern. 1982, 9anares-Alcantara et al., 1985). Calandranis and Stephanopoulns (1985) used a peratures of the hot and cold streams are among the parameters that are allowed to vary in the opti- Symbolies 3640 computer with the Flavors environment for the operability amz/y~s of MER heat mization. Finally, the latest developments in using various exchangernetworks. After the decision on the range LP models to handle interesting aspects of the retrofit of operating parameter variations for a given netsituation should he mentioned. Saboo et al. (1986a,b) work is made, the program reduces the problem space described procedures for finding the most economic by employing the underlying stream structure to way to add area to installed exchangers. A new assess the number of varying parameters. The postrade-off introduced is the number of exchangers to sible nonconvexity of the network is investigated he modified vs the total additional area. The existing using several theorems that are actually production utility consumption can be evaluated, not only rules of a KBS. A set of eight base cases is used for against the traditional minimum target for a given the analysis of the network. The base cases depend ATm and all structures, but also against a new target upon pinch point changes ("swaps"), the number of for the given structure and AT,,~ =0. The latter disturbances to the network and flowrate variations. situation is interesting, since adding area to an ex- The program does not only tell whether the network changer is fairly simple and may be done during is operable, but also which heat exchanger that ~ should be further investigated together with proper normal plant down time. guidance for retrofit action. The program is based on Knowledge based systems object oriented concepts that enable a transparent, Recently there has been a very rapid growth in the graphics-based, interaction with the designer. By using the proposed theorems, the problem application of knowledge based systems (KBS), or expert systems (ES), in many areas. This has not yet, under consideration can he broken down into simpler however, had any profound effect on the activity problems and solved. Strictly speaking, the theorems related to HENS. When Niida et ai. (1986) described presented are not rigorous insights, but used in experiences in using expert systems in process systems connection with the chosen inference mechanism, a engineering within a major Japanese engineering/ practical result is obtained within the given domain. construction company, there is typically very few This approach shows that it is possible to reduce the applications reported in heat exchanger network de- problem space and that big combinatorial problems sign, and more use in for instance separation systems. may he solved to a "sufficient" degree. Note that Jezowski and Kueiel (1979) approached the HENS optimality still is the designer's responsibility. Finally, in order to combine several techniques problem by generating a search tree of various structures. Jeeowski (1981) then applied AI techniques like for the synthesis of heat exchanger networks, ordered search to find the best network, and this Viswanathan and Evans (1986, 1987) used an expert system to generate and control the various modules work is continued by Jezowski et al. (1983) and involved in HENS. They used a modified transJezowski (1983). Hartmann et al. (1986) pointed out that a large shipment model for the minimum utility connumber of feasible structures exist for the HENS sumption, a modified transportation model for the problem, and with multiple criteria. Since some of minimum number of units and a new transshipment these evaluation criteria to some extent are "fuzzy", model using the flexible "out-of-kilter" algorithm for they devised a simple interactive expert system com- the forbidden match problem. The expert system bing heuristic rules and f u z z y theory. This approach controls the interaction between the various problem is also described in two other papers from the same aspects and generates the models to be solved. The research group (Zeising et al., 1985; Zheljev et al., apparent advantage of using knowledge based tech1985). Unfortunately, both the works of Jezowski niques in this manner is the flexibility in applying, or and coworkers and Hartmann and coworkers have to combine, several methods for subproblems within the disadvantage that the pinch decomposition is not HENS. It is easy to modify or improve the capabilities of the system by expanding the contents of the included, i.e. more "knowledge" is needed. Grimes et al. (1982) presented the HEATEX pro- knowledge base. gram which synthesized heat exchanger networks with minimum utility requirements. A heuristic apCOMPUTER TOOLS FOR HENS proach was used, and the rules were represented as a Norsk Hydro a.s has access to the following comset of 115 production rules. By using OPS3 RX, an early version of the today well known OPS5 knowl- puter programs for designing heat-exchanger netedge engineering tool, the performance of HEATEX works: ADVENT, HEXTRAN, INTERHEAT, was too slow to he really useful (Chowdhury, 1985). MAGNETS, RESHEX and SUPERTARGET, the Interestingly, by formulating the knowledge base for last one only as an early prototype version. We have HEATEX, a more detailed insight in the HENS also tried a copy of HENSYN from Rev and Fonyo problem was obtained and formed the basis for a new (1984), available in universities from the EURECHA

320

T. GUNDF.IL~Nand L. NAF.SS

committee. These programs will be briefly discussed together with other commercial or research type programs referenced in the literature. Two programs included in this review are based on the use of LP techniques to handle several aspects of HENS. Both RESHEX and MAGNETS are university developed, prototype software aimed at research studies, but with significant potential for future industrial applications. HEXTRAN obviously has some nice features for rigorous simulation, optimization and monitoring of heat exchanger networks, and is used worldwide in a large number of industrial companies. When it comes to synthesis, however, our company has decided to use ADVENT with its user-friendly, colour graphics, workstation-based approach as the primary tool. The core of this program is based upon the latest results in pinch technology. HEXTRAN

Being introduced in 1980, and described by Challand and O'Reilly (1981) HEXTRAN from Simulation Sciences Inc. was for a long time the only commercially available computer program addressing the design of heat exchanger networks on a general basis. This batch-oriented proram handles targeting, synthesis, optimization and rigorous simulation and rating of the network. The targeting and synthesis part is based on the work of Linnlioff and Flower (1978a), with the loop breaking algorithm of Su and Motard (1984) to reduce the number of heat exchangers while achieving MER. Challand et al. (1981) and Colbert (1982) introduced the idea of a double temperature approach (DTA), where the heat recovery approach temperature (HRAT) determines the relative position of the composite curves and thus the level of heat integration, while the exchanger minimum approach temperature (EMAT) is allowed to have values below HRAT, however, without changing the utility consumption. O'Reiily (1985) suggested that this was contradictory to the pinch principle which says that whenever heat is transferred across the pinch, there will he a double penalty on the utility side. What really happens in HEXTRAN is that an equal amount of heat is transferred from above to below pinch (with AT > AT=in) and from below to above (AT < AT,a). In pinch technology, this has been named criss-crossing (Ahmad, 1985), and may be used to optimize the total heat transfer area without changing the utility consumption. Another approach to this situation which actually describes the industrial fact that streams may have quite different heat transfer conditions, is the concept of individual stream contributions to AT,~. Applications of HEXTRAN were reported by Kleinschrodt and Hammer (1983) who used the DTA approach to handle multisheli situations in crude preheat trains. Recently, Jones et al. (1986) described how the program can handle retrofit situations in-

eluding targeting for payback time and synthesis through constraints to evolve the network for compatibility. A major advantage with HEXTRAN is a rigorous handling of streams (compositional) and their physical properties (data base) as well as the heat-exchanger calculations. ADVENT

An interactive, pinch technology based, tool was introduced commercially by Union Carbide Corp. in 1985. ADVENT operates in a workstation environment using colour graphics, multiple windows, full-screen panels and mouse to enhance the man-machine interface. Gautam and Linnhoff (1985) emphasized that the very nature of process synthesis requires interactive and flexible tools to give the engineer full control of the network development. Gautam and Smith (1985) presented the modules of ADVENT, which include HENS, process modification, optimization, heat and power, separation synthesis, operability/flexibility and furnaces. The program also has certain simulation capabifities, first of all meant to supply the synthesis activities, but also for detailed heat exchanger calculations. Running on a Unix-based process engineering workstation (Dickert et aL, 1985; Smith, 1986), Gautam and Smith claimed that ADVENT combines state-ofthe-art in computing technology and process synthesis. The program certainly has an edge when it comes to user interface and the application of pinch technology. First, the user can address the aspects of process modification, utility allocation and finding the optimal level of heat recovery by graphical aids like the composite and grand composite curves and an implementation of the "supertargeting" features described in the section on pinch technology. Then the network can be interactively generated in the stream grid environment, with additional tools like the driving force plot to evaluate each heat exchanger thermodynamically. An automatic synthesis algorithm based on the TI method with a subsequent unit reduction is also available, or the program may "fill in the rest" after the engineer has decided the crucial near-pinch matches. Finally network evolution includes graphical and caiculational tools for loop and path manipulations. INTERHEAT

As the name suggests, this program is basically interactive, but with some automated features for network generation and parameter optimization. The program was initially developed by professor Lzken at the Norwegian Institute of Technology in Trondheim, but an extension to the program, marketed under the name of HEATNET, has been developed in cooperation with National Engineering Laboratory (NEL) in the U.K. This also includes links to the PPDS physical property data base and

The synthesis of cost optimal heat exchanger networks to detailed heat exchanger calculations by HTFS software. The features of the program including temperature dependent physical properties (enthalpy and heat transfer coefficients), parameter variations, multivariable optimization (up to 9 parameters) and heat pump and turbine models are described in Lzken, (1984, 1985, 1986). Akselvoll and Lzken (1987) describe the most recent enhancements of the program on automatic network generation, where a number of near-optimai solutions are listed by treating the problem in a combinatorial way, obeying the pinch decomposition and introducing the minimum number of split streams that is necessary to obtain MER.

321

the traditional energy target obtained by PTA or LP models, which gives the minimum utility consumption for a specified value of ATm and all possible structures. RESHEX calculates the minimum energy requirements given the existing structure and ATm ffi 0. This may help in deciding between repiping and area addition to the existing structure. In the latter case, the program finds the minimum additional area necessary in order to achieve a certain reduction of the energy bill. The new trade-off between the number of exchangers to be modified and the corresponding total additional area can be addressed by penalizing area addition on selected heat exchangers. MAGNETS

The activity at Carnegie-Mellon University in Professor Grossmann's group on heat exchanger netLinnhoff March, process integration consultants, work synthesis and the application of MILP models have marketed a pure targeting program named has been implemented in an automatic synthesis TARGET II, based on university software developed program named MAGNETS (Grossmann, 1985b). in professor Linnhoff's group at UMIST, ManThe LP and MILP transshipment models of chester (see Linnhoff and Senior, 1983). The program Papoulias and Grossmann (1983) give rigorous concan handle constrained targeting to find the energy penalty of a forbidden match, and results are given strained targets for energy and minimum number of graphically (composite and GCCs, energy target plot units. The constraints take the form of forbidden, and stream grid with stream population above and required or restricted matches. The automatic network generation was made posbelow pinch). sible by the clever stream supersUucture described The idea of establishing the optimum value of AT..~ ahead of design presented by Ahmad and by Floudas et al. (1986). Based on MER heat load Linnhoff, (1984, 1986), Ahmad (1985) and Linnhoff distributions for the minimum number of units, the and Ahmad (1986a,b) forms the basis of a new optimal match sequence and stream split ratios are found which minimizes the investment cost through program named SUPERTARGET. In addition to the solution of a N I P problem. One advantage of the these extended targets for heat exchanger networks, program is that heat transfer coefficients and heatthe program gives the user some guidance in the match selection when establishing the network. Im- exchanger cost data are specified at a stage of development where the retches are known, but the netportant tools are the CP-table and the DFP. The work topology and stream structure still free to be philosophy taken is that the initialization through optimized. Good or bad heat transfer conditions as annual cost targeting to find a close to optimal A T ~ well as the possible need for exclusive materials of is the crucial task. The network development may construction can thus be accounted for. then be done by hand. When the minimum utility cost problem has been solved by the LP transshipment model, the user has RESHEX the option to obey or disregard one or more of the RESHEX was initially developed to aid the re- possible pinch points in the process, which means search on resilience of heat exchanger networks in that designs with units in the range between Um.u~a Professor Morari's group at the University of Wis- and Um can be studied. An outer (at present) manual comin and later CALTECH (see Saboo and Morari, iteration loop is necessary to address the problem of 1984). The latest developments of the program have cost-optimal A T e . A fimitation with the current been described by Saboo et al. (1986a,b). RESHEX approach is that nonlinear heat capacities cannot be can give rigorous comtmiued energy and area targets handled, but this is subject to further research at allowing for individual stream contributions to CMU. In pinch technology one obvious problem is that AT.~,. The resilience of a network structure can be tested for a specified disturbance range, or the re= there are many MER networks that can be developed silienee index can be computed. Automatic network by PDM. When optimizing these initial networks by generation is based on a modification of the M I I P loop-breaking and path manipulation, each network model of Papoulias and Grossmann (1983), where the will follow different patterns, and an exhaustive MER requirement can be relaxed in order to reduce search is n __ee~___ry to guarantee the global "optimum". A ~milar situation arises when using MILP the number of units. Some fentures are included in the program that are models. There will be many possible HLDs resulting of special interest in retrofit situations. In addition to in the minimum number of units, and each of these TARGET ll/SUPERTARGET

322

T. Gu~3~,~,s~Nand L. N~ss

will lead to different network topologies and stream splits. Other programs

A number of computer tools have been briefly described in the literature, but without enough information to make a proper review possible. The majority of these programs have been developed to support research in the HENS field. EAPEX, which was presented by Ono et al. (1982), is a large computer package involving both process simulation and heat integration, where the modules of the program communicate through files. The heat integration part is based on the approach of Naka and Takamatsu (1982). HENSYN is one half of the program SYNSET developed to aid teaching of process synthesis, which is available for univer~ties through the EURECHA committee. HENSYN is developed by Rev and Fonyo 0984) and combines the PDM and a modification of the fast matching algorithm of Ponton and Donaldson (1974) to develop nearoptimal initial MER networks with the fewest possible number of units. The modification from Ponton and Donaldson is simply that the minimum utility target for the remaining problem is calculated after each match to check whether MER is prevented by the last match. HENS from the Research Institue of Chemical Equipment in Brno, Czechoslovakia, is also based on a combination of methods from Linnhoff and coworkers and Ponton and Donaidson's matching rule. Since this procedure is very fast, one can afford to obtain networks for several values of AT,u, and thus find the optimum trade-off between capital and energy. The features of the program are described by Klemes and Ptacnik (1984, 1985) and include a link to the simulation program SIPRO. Later, Ptacnik and Klemes (1987) introduced HENS-II to handle multipass exchangers. The pinch decomposition is included and an inverse Ponton and Donaldson rule is applied below pinch, meaning that the design starts from the cold end. A microcomputer implementation of a synthesis algorithm based on the thermodynamic combinatorial method of Flower and Linnhoff (1980) has been presented by Govind et al. (1986). The simplifications involved in this approach reduce industrial applicability, but this does not prevent its usefulness as a teaching program. Jabali (1984) implemented and evaluated several thermodynamicbatedmethods in prototype software including the aspect of resilien~ and flexibility. Jezowski et al. (1983) developed computer programs for both tree.~,arch and assignment task algorithms. Depth-first search was applied for mudl-scale and ordered search for large-scale problems~ The assignment algorithm was implemented such that some of the inherent drawbacks were removed.

Finally, a tailor made program for crude oil preheat trains, named EXTRAIN, has been marketed by Pace Consultants in Houston, Texas (member of the Jacobs Engineering Group). The program handles heat exchanger rating and preheat train optimization. Simulation programs

Some of the programs mentioned above have simulation capabilities ranging from rigorous calculations to shortcut methods aimed at supporting the synthesis activity. In addition, there are several general-purpose simulation and design programs available on the market that will not be mentioned here. Rather, we will refer to a few contributions that attempt to exploit the special characteristics of heat exchanger networks in the simulation procedure. While, for example, EAPEX is a sequential modular program, Shindo et al. (1982) argued that the integrated nature of heat exchanger networks will lead to a high number of information recycle loops. The special purpose simulation, rating and design program for heat exchanger networks, HENS, that was presented in their work thus uses the equationbased solution strateg¢. Another equation-based program is the EROS flowsheeting package for evaluation and optimization of heat exchanger networks including mixers and splitters, presented by Shah and Westerberg (1980). Even though Simulation Sciences Inc. uses sequential modular techniques in their general-purpose flowsheeting program, PROCESS, they have chosen to use matrix methods in HEXTRAN. Klemes and Vasek (1983) suggested that in order to handle industrial heat exchanger systems as opposed to ideal counter current formulations, integrated software for the synthesis and simulation is required. Their HENS synthesis program is thus integrated with the SIPRO simulation program, which uses ~In equation-based solution approach. Finally, Asbj~rnsen and Haug (1985) discussed the choice of variables and approximations in the optimization of heat exchanger networks. RECOMMENDATIONS & FUTURE TRENI~

Systematic methods and computer tools are available today which can be used to address and solve the i ~ t r i o J heat e x c ~ network synthesis problem. The proper solution of the various trade-offs involved in heat integration has significant impact on the process economy. Large energy savings with short payback times that meet todays short-term based economic reality have been reported. In some cases, proper intewation can reduce both energy consumption and capital investment and even produce more operable processes. To remain competitive, industry will have to take advantage of the blend of methods presented in this review, i So far, pinch technology has produced the most extensive list of good projects in industry for ira-

The synthesis of cost optimal heat exchanger networks proved heat integration. H o w e v e r fundamental t o pinch technology, mathematical methods and knowledge based systems is the skill and experience of the engineer himself. Industrial companies should therefore put emphasis on technology transfer to increase the understanding of the technical and economical mechanisms involved. Universities should reduce their emphasis on unit operations when teaching process design and put more emphasis on the network structure and related 'important trade-offs in economy, operability, controllability, safety and flexibility. Process synthesis with its blend o f heuristic and systematic methods is fundamental in this respect. New "revolutions" in the design o f heat exchanger networks, resulting from discoveries like predesigu targets for best performance and the heat recovery pinch, are not expected. The use of thermodynamic principles in other areas of process synthesis, however, is still to he explored. Linnhoff uses the picture of the Onion Diagram and the Rubic Cube to describe the extensions and new areas for research within pinch technology. There is still significant scope for improvements in the mathematical methods that have been employed to automate the design of heat-exchanger networks, both in solution techniques, in added realism of the models to satisfy industrial needs, and in flexibility to make sure that the process synthesis activity remains designer driven. The use o f sophisticated hardware and software will be of importance for future development of HENS tools. The requirement for some high level software to integrate the various approaches to HENS is one incentive, another is the need to accommodate design engineers with tools that support the design activity, are user-friendly and have fast response. The use o f knowledge based techniques in some high-level environment seems inevitable in this respect, and will be of advantage in several areas. Pinch technology could gain from "intelligent" programs which behave differently according to the nature of the stream data and the environment of the heat exchanger network. Mathematical techniques that use some kind of superstructure will benefit from sophisticated software which ensure that the optimum solution is contained within the superstructure, while excluding unnecessary possibilities, thus reducing the problem space. Some of the qualitative considerations included in the optimization (safety, environmental aspects, company policy etc.) can be programmed using knowledge based techniques. Finally, we will bring across some of the ideas presented by Stephanopoulos (1986) on expert systems and computing environments for process design activities. Today's programming tools like F O R T R A N do not have the flexibility to meet tomorrow's needs. However, the advanced and flexible tools for program development now available within AI, are too slow for engineering applications. It is thus expected that design tools are developed on one

323

machine, but used in applications in a faster environ. ment on another maehme. The strong argument is made that the culture of personal computers has misdirected the capabilities in computer aided engineering. Stephanopoulos (1986) further suggested that the expert systems built for computer aided process design should take the form of an expert assistant or consultant rather than the expert solver.

Acknowledgements--The authors are grateful to Norsk Hydro a.s, for the opportunity to explore this industrially important area, and for supporting this publication. We also wish to express our thanks and admiration to the _rese__~,r.h groups that we have been in contact with during the last few years. NOMENCLATURE A ffi Heat tranffer at'ca Am •ffi Minimmn heat tranlfer area AI - Artif~ial intelligence C = Cooler (Figs 13-15) DFP = Driving force plot DTA ffi Double temperature approach E ffi Energy ED ffi Evolutionary development method EMAT = Exchanger rain approach temperature Em ffi Minimum energy consumption ES ffi Expert system FI ffi Flexibility index F r ffi Heat exchanger configuration factor GC'C ffi Grand composite curve h ffi Film heat tranffer coeffaciem H ffi Enthalpy or heater (Figs 13-15) HAF = Heat availability function HDS ffi Heat demand and supply diagram HENS ffi Heat exchanger network synthesis HI.X) ==Heat load distribution HRAT ffi Heat recovery approach temperature HTG ffi High temperature group (Fig. 10) I-ITS ffi High temperature streams (Fig. 10) KBS ffi Knowledge based system L ffi Number of heat load loops LP Linear pmgramming~ LRS ffi Long range stream (Fig. 10) LTG ffi Low temperature group (Fig. 10) LTS ffi Low temperature streams (Fig. 10) m,cp : Heat capacity flowrate : MER ffi Maximum energy recovery MILP = Mixed integer linear programming nh -- Number of hot streams nc ffi Number of cold streams N ffi Number of streams and utilities NLP = Nonlinear programming, PDM ffi Pinch design method PT = Pinch technology PTA ffi Problem table algorithm Q ffi Enthalpy AQ ffi Enthalpy change Qh ==Hot utility consumpiton Qx,m ffi Minimum hot utility Qc ffi Cold utility consumpiton Qc.m ffi Minimum cold utility RI ffi Resilience index RPA ffi Remaining problem analysis S = Number of subsystems ST ffi Steam (Fig. I 1) T = Temperature T~ffi Temperature of a cold stream Th ffi Temperature of a hot stream

324

T. GU~D~N

7", = Supply (start) temperature Tt = Target (end) temperature TC = Thermodynamic combinatorial method TI ==Temperature interval method TO =, Temperature enthalpy AT = Temperature difference ATI = Stream contribution to &T..~n &TLM= Log mean temperature difference ATI = Minimum approach temperature ATo~ = Optimum temperature difference (Fig. 16) A T ~ = Approach temperature at the pinch U = Heat transfer coefficient U = Number of units U m - Minimum number of units (global) UmiLMER ms Minimum number of units CMER)

itEFgEENCES Abroad S. Heat excl~amgm" networks: cost trade-offs in energy and capital, Ph.D. Thesis, University of Manchester, Inet of Sci. and Tee,hnol. (1985). Abroad S. and B. Linnhoff. Overall rest ~ for heat exdumger networks. IChemE Annl Res. Mtf, Bath, (t984). Abroad S. and B. IAnnhoff. SUPERTARGL~I': optimisation of a chemical solvents plant=-differentpmcem structures for differenteconomics. A S M £ Mtg, Anaheim, California (1986). Aimmd S., B. Limmhoffand R. Smith. D a p of multipau heat exchangers: an altenuttive atlqpmaeh.Journal of Heat Transfer(ASME) In press(1988). Akselvoll K. and P. A. L~um. Automatic heat eva.hanger network synthem. CEF#7--7be U~ of Computers Chem. J ~ l ' q . Taormina, Italy (1987). A_sbjmmen O. A. and M. Haug. Choice of variables and approximations in the optimization of heat exchanger networks. Proc. Int. Conf. Process Systems F2~£ (PSE8S), pp. 197-208, Cambridge, U.K. (1985). Banarm-Alcantara R. et o.1. KnowladF treed expert systems for CAD. Chem. F.a~£ ProS,. 81, 25-50 (1985). Beautyman A. C. and A. R. H. Cornish. The design of flexible heat exchanger networks. Pro¢. 1st U,K. Nat. Conf. Heat Transfer, Vol. 1, pp. 547-565, Leeds, U.K. (1984). Bingzhen C. ¢t al., ~ synthesis and semitivity analysis of heat e x c ~ networks. Proc. CHEMCOMP-82 Antwerp, Belgium (1982). Boland D. Energy ~ t : empham in the 80's. The Chem. Fmfr March 24--28 (1983). Boland D. and E. Hindmarsh. Heat e x c ~ network

improvements. Chem. &g,~ ProS. St, 47-84 (1984), Bohmd D. and B. IAnnhoff. The ix~iminary design of networks for heat m : h e n ~ by symmmtic methods. The C&nn. F.a,gr. April 9-15 (1979). Calamtranis I. and O. Stt,phanopuulm. Structural operability analysis of heat exchanger networks. AIChE Mt£. Houston (1985). Calandranis J. and G. Stephanopuulos. Structural operability analysis of beat exchangm- networks. IChemE, Chem. F,ngng Res. Des. 64, 347-364. Cena V., C. Mustacchi and F. Natali. Synthesis of heat exchange networks: a non-iterative approach. Chem. F.~,ng. Sci. 32, 1227-1231 (1977). Cerde J. and A. W. Westerberg. Synthesizing heat exchanger networks having restrkted stream/stream

mm~ ruing t m a , p u m e ~ prot/ma form~eom. Ch,m. ,.~i. ~ 1723-1740 (1983). C..erda J., A. W. Westerberg, D. Mason and B. Linnhoff. Minimum utility usage in heat exohanger network synth¢sis--a transportation problem. Chem. F.ng~ Sci. 38, 373-387 (1983).

and L. N~a.ss Challand T. B. and M. G. O'Reilly. New engineering software for energy conservation projects. Proc. 2nd Int. Conf. Engng Software, London, pp. 306-316 (1981). Challand T. B., R. W. Colbert and C. K. Venkatesh. Computerized heat exchanger networks. Chem. Engng Prog. 77, 65-71 (1981). Chowdhury J. Expert systems gear up for process synthesis jobs. Chem. Engng 92, 17-23 (1985). Colberg R. D. and M. Morari. Analysis and synthesis of resilient heat exchanger networks. Adv. Chem. Engmg. In preparation (1988). Colberg R. D., M. Morari and D. W. Townsend. Design of resilient processing plants--a flexibility index target for heat exchanger networks. In preparation (1986). Colbert R. W. Industrial heat exchange networks. Chem. Engng Prog. 78, 47-54 (1982). Dickert B. F., J. S. Newell and Y. P. Tang. An integrated process engineering software system in a workstation environment. AIChE Mtg. Chicago (1985). Doidan O. B., M. J. Baiptjewicz and J. C.erda. Designing heat exchanger networks for ex'mtin8 chemical plants. Comput. chem. Engng 9, 483-498 (1985). Duran M. A. and I. E. Grossmann. Simultaneous optimization and heat integration of chemical processes. A1ChE JI 32, 123-138 (1986). Dyson A. E. S. and P. J. Kenny. The optimization of oil refinery heat exchanger networks. Proc. lOth Aust. Chem. F.ngng Conf. Sidney, pp. 248-252 (1982). Elshout R. V. and E. C. Hohmann. The heat excha_~,u~_r network simulator. Chem. Engng Prog. ?S, 72-77 (1979). Floudas C. A. and I. E. Grossmann. Synthesis of flexible heat exchanger networks for mnitipefiod operation. Comput. chem. Engnf 1O, pp. 153-168 (1986). Floudes C. A. and I. E. Gr0~smann, Automatic gmefation of multiperiod heat exchanger network coe£w~ttione. Comput. chem. Engng II, 123-I'42 (1987a). Floudas C. A. and I. E. Groesmann~ Synthesis of flexible heat exchanger networks with uncertain flowrates and temperatures. AICAE Mtg. Houston, Texas (198To). Floudas C. A., A. R. Ciric and I. E. Groumann. AuWmatic synthesis of optimum heat e ~ network configurations. AIChE JI 32, 276-290 (1986). Flower J. R. and B. Linnhoff. Computer aided synthesis of energy recovery schemes. IChemE 'Conf. Comput. Chem. ~gng, Edinburgh (1977). Flower J. R. and B. Linnhoff. A thermodynamic combinatorial approach to the design of optimum heat exchanger networks. AIChE/I 26, 1-9 (1980). Gautam R. and B. Linnhoff. Computing tw.hnoiogy and engineering dialogue--interfaces as if people mattered. AIChE Mtg, Houston (1985). Gautam R. and J. A. Smith. A computer-aided system for process synthesis and optimization. AIChE Mtg, Houston (1985). Govind R., D. Mocsny, P. Cosson and J. Klei. E.xchah4ler network synthesis on a microcompuL-r. Hydrocarbon Process ~ , 53-57 (1986). Greenkorn R. A., L. B. Koppel and S. Raghavan. Heat exchanger network synthesis---a thermodynamic approach. AIChE Mtf, Miami (1978). Grimes L. E., M. D. Rychener and A. W. Westerberg. The synthesis and evolution of networks of heat exchange that feature the minimum number of units. Chem. Commun. 14, 339-360 (1982). Groumann I. E. Mixed-integer p m g r a ~ approach for the synthesis of integrated proems flowsheets. Comput. chem. Fa~gng9, 463-482 (1985a). Grossmann I. E. Personal Commtm/eation, CarneEieMellon University, Pittsburgh (1985b). Grossmann I. E. and R. W. H. Sargent. Optimum design of heat exchanger networks. Comlmt. e&nm. Eng~ 2, 1-7 (1978a). Grossmann I. E. and R. W. H. Sargent. Optimum d e ~ of

The synthesis of cost optimal heat exchanger networks chemical plants with uncertain parameters. AICA£ JI 24, 1021-1028 (1978b). Hall S. G. Capital cost targets for heat exchanger networks: differing materlah of comtrnetion and differing heat exchanger types. M.Sc. Thesis, University of Manchester, Inst. of Sci. and Technol. (1985). Hama A. Computer-aided synthesis of heat exchanger networks; an iterative approach to global optimum networks. Proc. First U.K. Nat. Conf. Heat Transfer, Vol. I, pp. 567-582, Leeds (1984). Hama A. and M. Matsumura. Stage-by-stage synthesis of ~ t op timal heat exchanger networks. Pro¢. (Technical Sessions) Int. Conf. Proc. Syst. Engng (PSE-82), pp. 106-113, Kyoto (1982).

325

orcim~ se~tch, inz. Chem. i Proc. (Poland) 2, 45-58 (1981). Jezowski J. Selected probim of heat exchani~ network synth~. Hung. J./m/. Chem. 10, 345-356 (1982a). Jexowski J. Heat exchanger network ffnthm/s by combined intepr prolirsmmi~--thc~modynall~ approach. Proc. CHEMCOMP-82, Antwerp, Belgium (1982b).

Jezowski J. Heat exchange network synthesis---method of generating solutions. Proc. 4th Conf. Appl. Chem. Unit Oper. Proc., Vol. 1, pp. 59-64, Verzprem, Hungary, (1983). Jexowski J. and A. Halat. A ~mpie assignment task based method of heat exchanger network synthesis. Proc. MATCHEM, Vol. II/1, pp. 241-250, Balatonfured, Hanson K. W. and A. R. H. Cornish. An algorithm to Hungary (1986). develop heat exchanger networks with the minimum Jezowski J. and E. Kuciel. Method of synthesis of the optimal network of heat exchangers. In=. Chem. (Poland) number of stream splits and maximum energy recovery. 9, 621--630 (1979). Proc. First U.K. Nat. Conf. Heat Transfer, Vol. I, Jezowski J., J. Stanislawski and A. Hakai. A comparison of pp. 583-598, Leeds (1984). some methods of heat exchanger ttetwork synthesis. Proc. Hartmsnn C. Process synthesis through appfication of mixed integer linear programming. M.Sc. Thesis (in Nor3rd Int. Conf. Inf. $ci. Che~n. Engng Vol. 2, Paris (1983). wqian). The Norwegian Inst. of Technol., Trondhelm Jones S. A. and D. W. T. Pippin: The genergtion of heat (1984). load distributions in heat exchanger network synthesis. Hartmann K., W. Kauschus and M. Wagenknecht. Optimal Proc. Int. Conf. Proceu $ysten~ Enh,ng (PSE-&,¢), design of chemical process systems by fuzzy methods. pp. 157-177, Cambridge (1985). Proc. MATCH£M, Vol. !, pp. 156-165, Balatonfured, Jones D. A., A. N. Yilmaz and B. E. Tilton. Synthesis Hungary (1986). techniques for retrofitting heat recovery systems. Chem. Hellgs P. J. Heat exchangers for heat recovery. 17th Int. Engng Prog. 82, 28-33 (1986). Syrup. High Temp.. Heat Exch. Dubrovnik, Yugoslavia Kafarov V. V. and L. A. Shmidt. Synthesis of an optimal heat-exchange system with consideration of the errors in 0985). Hendry J. E., D. F. Rndd and J. D. Sender. Synthesisin the the mathematical models of its elements. Doklady Chem. design of chemical processes. AIChE JI 19, 1-15 (1973). TeclmoL, Proc. Aead. ScL USSR 262-264,' 10-14 1982. Hindmarth E., D. Boland and D. W. Townsend. Maxi- Kafarov V. V. et al. Breakdown method for the automated mizing energy savings for heat engines in process plants. complex synthesis of heat-exchange systems. Doklady Chem. F~gng 92, 38--47 (1985). Chem. Technol. Proc. Acad. Sci. USSI~ 26~--264, 64-67 Hlavaeek V. Analysis and synthesis of complex plants. (1982a). S!~dy state and tramient behaviour. Proc. Comput. Des. Kafarov V. V. et aL Breakdown-variation method for the £rr. Chem. Plants, Vol. 3, pp. 903-972, Kariovy Vary automted complex synthesis of heat-exelmnge systems. 0975). Doklady Chem. Technol. Pro¢. Aead. Sci. USSR 265-267, Hlavacek V. Synthem in the des/gn of chemical proceau. 99-102 (1982b). Comput. chem. ~ 2, 67-75 (1978). Kafarov V. V. et al., An investigation of the properties of Hohnumn E. C. Optimum networks for heat exchange. the effgiency criterion used in the complex problem of Ph.D. Thetis, Univ. of Southern California (1971). designing heat-exchange systems with the objective of Hohmann E. C. Heat-Exchange Technology, Network Syndeveloping search algnrithms for the global extremum. thesis. Kirk-Otlener F.nc~lopnedia of Chemical TechDoMady Chem. Technol. Proc. Aead. $cL USSR 265.-267, nology, 3rd FAn, Suppl. pp. 521-545. Wiley New York l l l - l l 5 (1982c). (1984). Kardos J. and O. Strelow. Structural synthesis of heat Hohmann E. C. and F. J. Lockhart. Optimum heat exexchanger networks with standardized exchangers. Proc. changer network synthesis. AIChE Mtg. Atlantic City 4th Conf Appl. Chem., Unit Operation Proc., Vol. I, (1976). pp. 65-70, V ~ p m n , Hungary (1983). Huang F. and g. V. Elshout. Optimizing the heat recovery Kelahan R. C. and J. L. Gaddy. Synthesis of heat exchange of crude units. Chem. £ngng/hog. 72, 68-74 (1976). networks by mixed integer optimization. AIChE JI 23, Hwa C. S. Mathematical formulation and optimization of 816-822 (1977). heat exchanger networks using separable programming. Kesler M. G. and R. O. Parker. Optimal networks of heat AIChE-I Chem E Syrup. Ser. 4, 101-106 (1965). exchange. Chem. Engng Prog. Syrup. Ser. 6& i!1-120 Irazoqni H. A. Optimal thermodynamic synthesis of ther(1969). mal energy recovery systems. Chem. Engng ScL 41, Kleinschrodt F. and G, A. Hammer. Exchanger networks 1243-1255 (1986). for crude units. Chem. FJqp,g Prof. 79, 33-38 (1983). Itoh J.,K. Shiroko and T. Umeda. Extensive application of Kiemes J. and R. Ptacnik. Heat-exchange network synthesis the T--Q dingzam to heat integrated system synthesis. of production plants. CHISA.84~ Prague (1984). P~c. (T~mk=l S m ~ m ) Int. Coq[. Proc. Syst. Klemes J. and R. Ptacmk. Computer-aided synthesis of heat (PSE-82~ pp. 92-09. Kyoto (1982). exchange network. J. Heat Recmwry Systems $, 425--435 Ivakhaenko V. I., 0. M. Ostrovskfi ~nd T. A. Ber,~intkli. (1985). A method for optimal synthem of heat-exchun~ systems. Klemes J. and V. Vmek. Application ofa flowsheet program Tkee,'. Foum/at. Chem. ~ lf~ 250-255 (1982). in an interactive synthesis of heat exchanger networks. Jahali S. Contribution a la syntham optlmale des reseaux, Proc. 3rd Int. Co~. lnfo. $cL Chem. Engng Vol. I, d'edhangeun de ehalem'. Ph.D. Thm/s 0n French), Ecole pp. 82/!-6, Paris (1983). Ceatrnie des Arts et Manufactures, Soutenne, France Kobayashi S., T. Umeda and A. Ichikawa. Synthesis of (1984). optimal heat exchange systema--an approach by the Jakob K. et aL Second law analym of industrial piants-optimal auignment problem in linear programming. aspects, Pro¢. MATCH£M, Vol. II/I, pp. 220Chem. F,ngng SoL 26, 1367-1380 (1971). 230, Balatonfured, Hungnry (1986). Kotjabasakis E. and B. Linnhoff. Sensitivity tables in the Jezowski J. Heat exchanger network synthesis algorithms of design of flexible processes: !. How much contingency in

326

T. G U N D ~ N

heat exchanger networks is cost-effective? IChemE Chem. Engng Res. Des. 24, 197-211 (1986). Lee K. F., A. H. Masso and D. F. Rudd, ~amch and bound synthesis of integrated process designs. Ind. En8~ Chem. Fundam. 9, 48-58 (1970). Li Y. and R. L. Motard. Optimal pinch approach temperature in heat exchanger networks. Ind. Engng Chem. Fundam. 25, 577-581 (1986). Linnhoff B. Thermodynamic analysis in the design of process networks. Ph.D. Thesis, University of Leeds (1979), Linnhoff B. New concepts in thermodynamics for better chemical process dmgn. Trans IChemE, Chem. Eng. Res. Des. 61, 207-223 (1983). Linnhoff B. Who feels the pinch. Personal view. The Chem. F,,,~gr. 62-63, February (1985). Linnhoff B. The process/utility interface. ~ Int. Mtg National Use of Energy, Liege, Belgium (1986a). Linnhoff B. Pinch technology for the synthesis of optimal heat and power systems. ASME Mtg Anah~ra, California (1986b). Linnhoff B. and S. Ahmad. SUPERTARGETING, or the optimisation of heat exchanger networks prior to design. World Cong. III,Chem. Enos, Tokyo (1986a). Linnhoff B. and S. Ahmad. SUPERTARGETING: optimum synthesis of energy management systems. ASME Mtg, Anahaim, California (1986b). Linnhoff B. and K. J. Carpenter. Energy conservation by exergy analysis---the quick and simple way. Proc. 2rid World Cong. Chem. Engng Montreal, Canada (1981). Linnhoff B. and J. R. Flower. Synthesis of heat exchanger networks: I. Systematic generation of m~,rgy optimal networks. AIChE JI 24, 633--642 (1978a). Linnhoff B. and J. R. Flower. Synthesis of heat exchanger networks: II. Evolutionary generation of networks with various criteria of optimality. AIChE JI 24, 642-654 (1978b). Linnhoff B. and E. Hindmarsh. The pinch design method for heat exchanger networks. Chem. F,ngng Sei. 38, 745-763 (1983). Linnhoff B. and E. Kotjabasakis. Design of operable heat exchanger networks. Proc. U.K. Nat. Conf. Heat TraMfer, VoL 1, pp. 599--618, Leeds (1984). Linnhoff B. and E. Kotjahasakis. Downstream paths for operable process design. Chem. Engnf Prof. 82, 23-28 (1986). Linnhoff B. and S. Parker. Heat exchanger networks with process modification, lChemE Annl Res. Mtg Bath (1984). Linnhoff B. and P. Senior. Energy ~ clarify scope for better heat integration, Proc. F~gnf March, 29-33 (1983). IAnnhoff B. and R. Smith. Design of flexible plant--how much capital for how much flexibility. AIChE Mtg, Houston (1985). Linnhoff B. and T. N. Tjoe. Pinch technology retrofit: setting targets for existing plants. AIChE Mtg, Houston (1985). Linnhoff B. and D. W. Townsend. Designing total energy systems. Chem. Enfnf Prof. 78, 72-80 (1982). Linnhoff B. and J. A. Turner. Simple concepts in process synthem give energy savings and elegant detigns. The Chem. F,~fr. December, 742-746 (1980). Linnhoff B. and J. A. Turner. Heat-recovery networks: new insights yield big savings. Chem. Engnf 8& 56-70 (1981). Linnhoff B. and D. R. Vredeveld. pinch technology has come of age. Chem. Engng Prof. gO, 33--40 (1984). Linnhoff B. and W. D. Withereil. pinch technology guides retrofit. O//& Gas J. 114, 54-65 (1986). Linnhoff B. et al. User Guide on Process Integration for the Efficient Use of Energy. Inst. Chem. Engrs, U.K. 0982). Linnhoff B., H. Dunford and R. Smith. Heat integration of distillation columns into overall processes. Chem. Engng Sci. 38, 1175-1188 (1983).

and L. NA~S Linnhoff B., D. R. Mason and I. Wardle. Understanding heat exchanger networks. Comput. chem. Engng 3, 295-302 (1979). Liu Y. A. A practical approach to the multiobjective synthesis and optimizing control of resilient heat exchanger networks. Proc. Am. Control Conf, Vol. 3, pp. I 115-1126, Arlington,. U.S.A. (1982). Liu Y. A., F. A. Pehler and IX R. Cahela. Studies in chemical process design and synthesis. Vll: systematic synthesis of multipass heat exchanger networks. AIChE JI 31, 487--491 (1985). Lzken P. A. Process integration of heat pumps. 2rid Int. Syrup. Large Scale Applic. Heat Pumps, York (1984). Lzken P. A. Interactive computer program used on the retrofit of a dewatering process. AIChE Mtg, Houston (1985). .. Leken P. A. Energy conservation in an industrial .process. World Cong. 111 Chem. En£ng Tokyo (1986). MarseUe D. F., M. Morari and D. F. Rudd. Design of resilient processing plants: II. Design and control of energy management systems. Chem. Engng $ci. 37, 259-270 (1982). Masso A. H. and D. F. Rudd. The synthesis of system designs. II. Heuristic structuring. AIChE JI 15, 10-17 (1969). McGallard R. L. and A. W. Westerberg. Structural sensitivity analysis in design synthesis. Chem. Engng J. 4, 127-138 (1972). Menzies M. A. and A. 1. Johnson. Synthesis of optimal energy recovery networks using discrete methods. Can. J. Chem. Engng 50, 290-296 (1972). Mitson R. J. Number of shells vs number of units in heat. exchanger network design. M.Sc. Thesis, Univ. of Manchester, Inst. of Sci. and Technol. (1984). Mocsny D. and R. Govind. Decompos/tion strategy for the synthesis of minimum unit heat exchanger networks. AIChE Jl 30, 853-856 (1984). Morton R. J. and B. Linnhoff. Individual process improvements in the context of site-wide interactions. IChemE Aani ICes. Mt£, Bath (1984). Muraki M. and T. Hayakawa. Practical syntham method for heat exchanger network. J. Chem. Engng Japan IS, 136-141 (1982). Naess L. and T. Gundersen. Supertarget~prototype software for heat exchanger network design. Confidential Norsk Hydro Report, Porsgrunn (1984). Naka Y. and T. Takamatsu. Design of heat exchanger networks in the heat availability diagram. Proc. CHEMCOMP-82. Antwerp, Belgium (1982). Nash D. B., E. C. Hohmann, J. Beckman and L. Dobrzanski. A simplified approach to heat exchanger network analysis. AICh£ Mtg, Philadelphia (1978). Nida K. et al. Some expert system experiments in process engineering. IChemE, Chem. Engng Res. Des. 64, 372-380 (1986). Nishida N., S. Kobayashi and A. lchikawa. Optimal synthesis of heat exchange systems. Necessary conditions for minimum heat transfer area and their application to systems synthesis. Chem. F,nl~g Sci. 7.6, 1841-1856 (1971). Nishida N., Y. A. Liu and L. Lapidus. Studies in chemical process design and synthesis: I!I. A simple and practical approach to the optimal synthesis of heat exchanger networks. AIChE J! 23, 77-93 (1977). Nishida N., G. Stephanopoulos and A. W. WestergerS. A review of process synthesis. AIChEJ127, 321-351 (1981). Nishimura H. A theory for the optimal synthesis of heat exchanger systems. J. Opt. Theory Appfic. 3a, 423-450 (1980). Nishitani H., E. Kunugita and L. T. Fan. On the vector optimization of heat exchange. J. Chem. Engng Japan IS, 475-480 (1982). Ono K. et al. A computer program for energy analysis and

The synthesis of cost optimal heat exchanger networks

327

Saboo A. K., M. Morari and D. C. Woodcock. Da/gn of resilient processing plants: VIII. A resilience index for heat exchanger networks. Chem. F.nffof Sci., 40, O'ReiUy A., Experiences in protein integration. The Chem. Emgr. May, 56-59 (1986). 1553-1565 (1985). O'Reilly M. Personal view. The Chem. Engr. January, Schmidt A. Exergy cotmideratiom in the dmdp of heat exchangers. Proc. 4tk ltai.-Yus-Amtr. Chem. Eng. Conf. 46-47 (1985). Grado, Italy, Vol. 1, pp. 7-13 (1984). Ortrovsky G. M. el al. Synthesis of heat exchanger Shah J. V. and A. W. Westerber8. Evolutionary synthesis networks. Hung. J. Ind. Chem. 13, 107-119 (1985). of heat-cachanger networks. AIChE Mtg, Los Aagglm Paponlias S. A. and I. E. Groumann. A structural optimization approach in process synthesis---II. Heat recov(1975). ery networks. Comput. chem. Engof 7, 707-721 (1983). Shah J. V. and A. W. Westerberg. EROS: a program for quick evaluation of en~gy recovery systems. Compul. ParkJrmon A. R. et al. The optimal design of resifient heat exchanger networks. AIChE Syrup. Set. 78, 85-98 (1982). chem. Eofng 4, 21-32 (1980). Pehler F. A. and Y. A. Liu. Thermodynamic availability Shindo A., C. Nagamwa and T. Maejima. Computer aided analysis in the synthesis of optimum energy and miniheat exchange system design. Proc. CHEMCOMP-82, mum cost heat exchanger networks• ACS Syrup. Ser. 235, Antwerp, Belgium (I 982). 161-178 (1983). Shiroko K. and T. Umeda. A practical approach to the Pehler F. A. and Y. A. Liu. Studies in chemical process optimal design of heat exchange systems. Proc. Econ. Int.' design and synthesis: VI. A thermodynamic approach to 3, 4449 (1983). the evolutionary synthesis of heat exchanger networks. S'drola J. J. Status on heat exchanger network ryathem. Chem. Engng Commun. 2S, 295-310 (1984). AICkE Mrs. Tulm (1974). Pho T. K. and L. Lapidus. Topics in computer-aided Smith J. A. User and program interface for an integrated UNIX based prcgeu engineering workstation. AIChE design: II. Synthesis of optimal heat exchanger networks Mtg, New Orleans (1986). by tree searching algorithms. AIChE JI 19, 1182-1189 Steinmetz F. J. and M. O. Cinmey. Total plant process (1973). energy integration. Chem. Engng Prof. 81, 27-32(1985). Ponton J. W. and R. A. B. Donaldson. A fast method for the synthesis of optimal heat exchanger networks• Chem. Stephanopoulos G. Expert systems and computing environments for process systems engineering. Corsmmm. AICh£ Engng Sci. 29, 2375-2377 (1974). 9, 8-17 (1986). Ptacnik R. and J. Klemes. An improved approach to synthesis of multipass heat exchanger networks. Proc. Su J.-L. and R. L. Motard. Evolutionary synthesis of heat' exchanger networks. Comput. chem. Engng 8, 67-80 CEF-87. Use Comput. Chem. Engng. Taormina, Italy (1984). (1987). Raghavan S. Heat exchanger network synthesis. Ph.D. Swaney R. E. and I. E. Grossrnann. An index for operational flexibility in chemical process design. I. FormuThesis, Purdue University, U.S.A. (1977). lation and theory. AIChE JI 31, 621-630 (1985a). Rathore g. N. S. Prcom reNqnencing for energy conserSwaney R. E. and I. E. Grossmann. An index for opervation. Chem. Engng Prof. 78, 75-82 (1982). ational flexibility in chemical process design• II. ComRathore R. N. S. and G. J. Powers. A forward branching putational algorithms. AICh£ J! 31, 631--641 (1985b). scheme for the synthesis of energy recovery systems. Ind. Takamatsu T., I. Hushimoto and H. Ohno. Optimal design Engng Chem. Process Des. Dev. 14, 175-181 (1975). of a large complex system from the viewpoint of sensiRev E. and Z. Fonyo. Synthesis of heat exchanger nettivity analysis. Ind. Emgng Chem. Process Des. Dev. 9, works. Chem. Engng Commun. 18, 97-106. (1982). 368-379 (1970). Rev E. and Z. Fonyo. Process synthesis of heat exchanger networks. Proc. 4th Conf. Appi. Chem., Unit Operations Tjoe T. N. Retrofit of heat exchanger networks. Ph.D. Proc., VOl. 1, pp. 71-76, Vet-zpt~n, Hungary (1983). Thesis, Univ. of Manchester, Inst. of Sci and Tedmol. Rev E. and Z. Fonyo. SYNSET~Process Synthesis Teach(1986). ing Programs, User Manual. Technical Univ. of Tjoe T. N. and B. Linnhoff. Heat exchanger network retrofits, lChemE Anni ICes. Mtg, Bath. (1984). Budapest, Hungary 0984). Rev E. and Z. Fonyo. Hidden and pseudo pinch phenom- Tjoe T. N. and B. Linnhoff. Using pinch technology for process retrofit. Chem. Engng 93, 47--60 (1986). ena and relaxation in the synthesis of heat-exchange networks. Comput. chem. Engng 10, 601-607 (1986). Townsend D. W. and B. IAnnhoff. Heat and power Rudd D. F. The synthesis of system designs: I. Elementary networks in process design. I: criteria for placement of decomposition theory. AIChE JI 14, 343-349 (1968). heat engines and heat pumps in process networks. AIChE JI 29, 742-748 (1983a). Saboo A. K. and M. Morari. Design of resilient processing plants: IV. Some new results on heat exchanger network Townsend D. W. and B. Linnhoff. Heat and power networks in process design. I1: design procedures for synthesis. Chem. Engng Sci. 39, 579-592 (1984). equipment selection and process matehing. AIChE J129, Saboo A. K., M. Morari and R. D. Colberg. RESHEX-an interactive software package for the synthesis and 748-771 (1983b). analysis of resilient heat exchanger networks---l. Program Townsend D. W. and B. Linnhoff. Surface area targets for heat exchanger networks. IChemE Ansi Res. Mtg, Bath description and application. Comput. chem. Engng I0, 577-589 (1986a). (1984). Saboo A. K., M. Morari and R. D. Colberg. RESHEX--an Townsend D. W. and M. Morari. ~ of heat interactive software package for the synthesis and anal- • e~hanl~" networks: an objective compatible with minimum cost. AIChE Mtg, San Francisco (1984). ysis of resilient heat exchanger networks---ll. Discussion of area targeting and network synthesis algorithms. Umeda T., J. ltoh and K. Shiroko. Heat exchange system synthesis by thermodynamic approach. Proc. Paeif. Comput. chem. Engng 10, 591-599 (1986b). Saboo A. K., M. Morari and R. D. Colberg. Resilience Chem. Engng Cong. 2, 1216-1223 (1977). analysis of heat exdumger networks--l. Temperature Umeda T., J. ltoh and K. Shiroko. Heat exchange system synthesis. Chem. F.mgng Prof. 74, 70--76 (1978). dependent heat capacities. Comput. chem. Engng I1, Umeda T., T. Harada and K. Shiroko. A thermodynamic 399--408 (1987a). approach to the synthesis of heat integration systems in Saboo A. K., M. Morari and R. D. Colberg. Resilience chemical processes. Comput. chem. F~gng 3, 273-282 analysis of heat exchanger networks--II. Stream splits and flowrate variations. Comput. chem. Engof I1, (1979a). Umeda T., K. Niida and K. Shiroko. A thermodynamic 457--468 (1987b).

synthesisof chemical ~ with the available energy. Proc. CHEMCOMP-82, Antwerp, Belgium (1982).

HRS I0t4--C

328

T. GUNDV.tEN and L. N~ess

approach to heat integration in distillation systems. AIChE J! 2S, 423-429 (1979b). Vinograd D. L. et al. Synthesis of heat-exchange networks by a thermodynamic method in the presence of several heating agents and refrigerants. DoMady Chem. Technol. Proc. Acad. Sci. U$$R 2t[&-270, 19--23 (1983). Viswanathan M. and L. B. Evans. Development of an expert system to synthesize heat integration systems for process plants. Energy Engng Sci. Syrup. Chicago (1988a). Viswanathan M. and L. B. Evans. Studies in the heat integration of chemical process plants. AIChE JI 1781-1790 (1987). Wells O. and M. Hodllkinmn. The heat content diageam way to heat eae~0mqer networks. Proc. Eqvq[ Aulgtmt, 59--63. (1977). Westbrook G. T., Use this method to size each ~ for best operation. Hy~ocarb. Process Petrol. Refm. 40, 201-206 (1961). Weeterberg A. W. The heat path diagram for energy mmummamt. AICItE ~ Jubilee Mtg, Waslfinllton, D.C. (1983). Westerberli A. W. and I. E. Grmmuan. Process synthesis techniques in the chmaieal industry and their impact on eneqly u,e. Report. Chem. Engng Dept, CMU. Pittsburgh (1985).

Wilcox R. Synthem of heat and power systems in clmmcal plants. Ph.D. Thins, Massachusetts Inst. of Technol. (1985). Wood R. M., R. J. Wilcox and I. E. Groesmann. A note on the minimum number of units for heat exchanger network synthesis. Chem. Engng Commun. 39, 371-380 0955). Zachoval J. and Z. Konecny. An improved method for the generation of heat exchanger networks, espec~lly suitable for the networks formed of ctomnow exclumgers. Proc. CHEMCOMP-82, Antwerp, Belllium (1982). Zar.hoval J., Z. Konecay and O. Navmtil. New ways in heat-exchange network design. CHISA.M, Prague (I9S4). Zehting O., M. Walpmknw,ht and K; Hartma-n Fin Unzeharfer Aisorithmtm zur Syuthme you Warmeubemqpmlpyetmnen. Wlum,.sc& ~ t dev Tecbn. HocKuqnde "Carl Sc/uw/enmt,r" VoL 27, pp. 62--69. Leuna-Men~ D.D.R. 0985). Zheljev T., M. Wallmzkne~t, K. Hartmann and W. Kauschus. Vmlleichende Analym yon Synthemmthnden zur Strukturienms von Wsunmmlmmqpm|-sysmn~ (WUS). Wtumuc*. Ze/ucb#~ ,br T~Jm. Hod,Jdude "Carl $chor/emmer", Vol. 27, pp. 70-89. Leuna. Merseburg, D.D.R. (1985).