22

Multiple Effect Evaporator System

Keywords: Multiple Effect Evaporator, Feed Arrangements, Optimum No. of Effects
The chief factor influencing the economy of an evaporator system is the number of effects. By
increasing the number of effects we can increase the economy of an evaporator system. The first effect of
a multiple effect evaporator is the effect to which the raw steam is fed, vapors obtained from first effect
act as a heating medium for another effect.

Different types of feed arrangement of multiple effect evaporators -:

1)FORWARD FEED ARRANGEMENT

In this arrangement the feed and steam introduced in

the first effect and pressure in the first effect is highest and pressure in last effect is minimum, so
transfer of feed from one effect to another can be done without pump.

Vapor

Vapor

Vapor

P1,T1
F
V
P0,T0

Liquid
Condensate

Liquid
Condensate

Liq. Concentrate
Condensate

Fig. 22.1: Forward feed arrangement in Multi effect Evaporator

2) BACKWARD FEED ARRANGEMENT: In this arrangement feed is introduced in last effect and
steam is introduced in first effect. For transfer of feed, it requires pump since the flow is from low
pressure to higher pressure. Concentrated liquid is obtained in first effect.

Vapor

Vapor

Vapor

P1,T1

F
V

P0,T0

Condensate

Liquid
Condensate

Liquid
Condensate

Fig. 22.2: Backward feed arrangement in Multi effect Evaporator

Condensate

3) MIXED FEED ARRANGEMENT: In this arrangement feed is introduced in intermediate effect,

flows in forward feed to the end of the series and is then pumped back to the first effect for final
concentration. This permits the final evaporation to be done at the highest temperature.

Vapor

Vapor

Vapor

P1,T1

P0,T0

Liquid
Liquid
Condensate

Condensate

Liq. Concentrate

Condensate

Fig. 22.3: Mixed feed arrangement in Multi effect Evaporator

ANALYSIS OF MULTIPLE EFFECT EVAPORATOR SYSTEM :

With the increasing trend in the cost of coal, fuel oil it becomes to use the vapors of previous effect in
the steam chest of following effect. This requires the multiple effect evaporator system .As the number of
effect increases the steam economy increases on the other side capital cost will be more.
There is economic balance between the fixed cost and the operating cost so that one can select the
optimum number of effects.
Fixed cost :
According to Coston and Lindey,the annual fixed cost of a multiple effect evaporator is approximately
proportional to the 0.75 power of the number of effects.

V1 = C1/A* N0.75
The estimated cost of a single effect (C1) can be obtained from a number of sources once the heat transfer
surfaces requirements are known.
Operating cost :
Operating costs can be divided in to steam cost and all operating costs (labour, cooling water,power and
maintenance) such that
V1 = h*W*C2/S + V0
H = operating time (hr/yr)
W = evaporation rate (Kg/hr)
C2 = cost of steam (Rs/Kg)
S = steam economy
V0 = all operating cost other than the cost of steam (Rs/Yr)
Steam economy can be expressed as:
S = S1+S 1S2+S 1S22+..+S 1S2N-1
S = S1 (1-S2N)/ (1-S2)
Thus V2 becomes
V2 = (1-S2)h*W*C2/S1 (1-S2N) + V0
Total cost
VT = V1 + V2
VT = = C1/A* N0.75 + (1-S2)h*W*C2/S1(1-S2N) + V0
For minimum cost VT/N = 0
0 = C1/N{(N+1)0.75 (N0.75)} + (1-S2)h*W*C2/S1{[1/(1-S2N+1)]-[1/(1-S2N)]}
Ahwc2/C1 = {[(N+1)0.75(N0.75)] * S1(1-S2N)(1-S2N+1)}/(1-S2)2 *S2N

Modeling of multiple effect evaporator system

Multiple effect evaporators involve a large number of state and design variables. A change in any
variable can upset the operation of the evaporator. To achieve the goal of energy conservation in multiple
effect evaporators, it is necessary to know how does steam economy alter with changes in operating

variables for a given end product concentration . To quantise the changes in steam economy , a functional
relationship correlating it with variables should be developed. For this, it is necessary to identify all the
variables which affect the steam economy of a multiple effect evaporator.
Variables of a multiple effect evaporator :
In a evaporator , the variables can be classified as geometrical-operating, and self balancing variables. As
regards the geometrical variable , it is the area of heat transfer surface in each effect of an evaporator .
Hence, N-effect evaporator will have N number of geometrical variables.
From industrial practices .we know that there are some operating variables which plant engineer can
change them independently to annual any imbalance in the operation of an evaporator . They include:
feed temperature, feed concentration ,feed flow rate,and steam temperature (pressure), saturation
temperature (pressure) in the last effect. Feed arrangement (forward/backward/mixed) is also one of the
operating variables. Thus , total number of operating variables is six.
As regards the vapour and liquid streams from each effect of a multiple effect evaporator, they cannot be
changed independently by a plant engineer. Therefore, they are self balancing streams. The variables
associated with these streams are: flow rate, temperature and concentration of liquid streams ; and
saturation temperature(pressure) of each effect. However, temperature of vapour stream equals to the
temperature of liquid stream. In this way , for N-effect effect evaporator the number of self balancing
variables becomes 5N. It is important to point out here that the saturation temperature(pressure) of the
last effect, has already been taken in to account as an operating variable . Therefore, it cannot be
considered as a self- balancing variable. Flow rate of steam to the first effect is the another self-balancing
variable whose value is usually not altered. Thus the total number of net self balancing variables for Neffect evaporator, becomes 5N [=5N-1+1].
The summation of geometrical, operating and self -balancing variables gives the total number of variables
in an evaporator. They are equal to 6N+6[=N+6+5N].
Mathematical model:
A mathematical model of a multiple effect evaporator is a relationship amongst the geometrical, operating
and self -balancing variables. This can be obtained from the equations of material balance, energy
balance, heat transfer rate, and boiling point rise.
For the simplicity of the mathematical model, following assumptions have been made in this analysis:
1.
2.
3.
4.
5.
6.

The vapours entering in to steam chest of respective effects are at their saturation temperature.
There is no sub cooling of the condensate from different steam chests.
Condensation of vapour in steam chest occurs at constant pressure.
There is no carry- over of liquid droplets with vapors leaving the respective effects.
There is no heat dissipation to surroundings.
Heat transfer surface does not undergo fouling.

Design of multiple effect evaporator without boiling point elevation for forward feed :
Equations are developed for the case where boiling point elevations are negligible ,also the effect of
composition on liquid enthalpy is neglected. The equations so obtained are generalized for the case where
boiling point elevations cannot be neglected. For definiteness , forward feeds are employed.
Specifications: F , XF,TF,To,Po,P3(ort3),X3(orl3),U1,U2,U3,equal areas, forward feed, negligible boiling
point elevations.
To find : V0, T1,L1,T2,L2 and A
Actually ,four additional dependent variables exist; namely: V1,V2,X1,X2

V1=F-L1(vapour)

V2=L1-L2 (vapour)

V3=L2-L3(vapour)

P31,T11

P3,T2

P3,T3

P0,T0
P0,T0

P1,T1

P2,T2

TF
V0

L1,T1

L2,T2

V0 (P0,T0)

F-L1 (P1,T1)

L3,T3
L1-L2 (P2,T2)

Fig. 22.1: Triple Effect Evaporator with Forward feed Arrangement

Equations of energy balance and heat transfer rate for all three effects are-:
Effect -1:
Enthalpy balance:
F(hf-h1)+V00=(F-L1)1

22.1

Rate equation:
U1A(T0-T1)=V0
Effect -2:

22.2

Enthalpy balance:
L1(h1-h2)+(F-L1)=(L1-L2)2

22.3

Rate equation:
U2A(T1-T2)=(F-L1)1

22.4

Effect -3 :
Enthalpy balance:
L2(h2-h3)+(L1-L2)2=(L2-L3)3

22.5

Rate equation:
U3A(T2-T3) = (L1-L2)2
The above equation describing the triple effect evaporator system constitute a set of nonlinear algebraic
equations that may be solved in variety of way, one of which is Newton raphson method which we will
learn in next lecture.
Design of multiple effect evaporator with boiling point rise for forward feed:
Most solutions that are concentrated in evaporators are mixtures of water and non-volatile salt. The
boiling temperature of the solution some times depend on the salt concentration. The difference between
the temperature(T) of boiling solution and the temperature of boiling water (pure) at the same pressure is
known as boiling point rise(BPR).thus
T = TW + BPR
BPR is a function of the solute concentration. A graph called the Dhuring chart is commomnly used to
determine BPR.
Calculation procedure:
Specifications:
F,XF,TF,thick liquid composition(XP),P,saturated steam pressure(P0),heat capacity(CP), overall heat
transfer coefficient(U).
To find: Heat transfer area(A)
Step-1:
Corresponding to evaporator temperature find out the boiling point of pure water,TW.
Step-2:
From the figures/empirical correlation determine the solution temperature/BPR at TW,XP. This
temperature is also the temperature of the superheated water vapour leaving the evaporator.

Step-3:
A total mass balance and a component material balance are used to calculate the flow rates.
F=V+L
XF F = LXP
Step-4:
Calculate /determine the enthalpies of three process streams
1) Enthalpy of water vapour from the superheated steam tables by solution temperature and pressure
of the evaporator.
2) The enthalpy of the solution can be calculated from there heat capacities .
H = CP(T-TR)
It can also be calculated by the plots or empirical correlations.
3) Latent heat of vaporization taken from steam table at corresponding steam pressure.
Step-5:
Now , write the enthalpy balance
F(hf-hp) + V0 0 = (F-L)(HV hp)
Or, V0 = [ (F-L)(HV hp) F(hf-hp)]/0
Step-6:
From the rate equation
V0 0 = UA(T0 1)
A = V00/U(T0-1)

Equations for all three effects are-:

Effect -1:
F[h(Tf,Xf)-h(1,X1)]+ V0 0 -(F-L1)[H(1) h(1,X1)]=0
U1A1 (T0 1) - V0 0 =0
XF F L1 X1 = 0

Effect -2:
L1[h(1,X1)-h(2)]+ (F-L1)[H(1)h(T1)]- (L1-L2)[H(2)-h(2)]=0
U2 A2 (T1 2)- (F-L1)[H(1)-h(T1)] = 0
L1 X1 L2 X2 = 0

Effect -3:
L2[h(2,X2)-h(3, X3)] + (L1-L2)[H(2) h(T2)] (L2-L3)[H(3)-h (3, X3)]= 0
U3 A3 (T2 3)-(L1-L2)[H-h]=0
L2 X2 L3 X3 = 0