Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
13 views

Lecture2 PDF

This document discusses mathematical modeling of continuous distillation columns. It describes the steady state modeling approach using mesh equations that account for mass balance (M), equilibrium (E), sum (S), and enthalpy balance (H). The mesh equations relate variables such as compositions, temperatures, pressures, and flows across theoretical plates in the column. The modeling requires specifying feed conditions and fixed parameters like heat inputs, side draws, and reflux ratio. The unknowns are compositions, temperatures, and liquid/vapor flows that are solved using methods like partitioning and coupling or global and simultaneous procedures.

Uploaded by

Ajay Sagar
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
13 views

Lecture2 PDF

This document discusses mathematical modeling of continuous distillation columns. It describes the steady state modeling approach using mesh equations that account for mass balance (M), equilibrium (E), sum (S), and enthalpy balance (H). The mesh equations relate variables such as compositions, temperatures, pressures, and flows across theoretical plates in the column. The modeling requires specifying feed conditions and fixed parameters like heat inputs, side draws, and reflux ratio. The unknowns are compositions, temperatures, and liquid/vapor flows that are solved using methods like partitioning and coupling or global and simultaneous procedures.

Uploaded by

Ajay Sagar
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

26

Mathematical modeling of Continuous distillation


Column
Keywords: Distillation Column, Study State Modeling, MESH, Degree of Freedom
The numerical simulation of multi-component separation process requires the following
points: Reliable predictive method for computing phase equilibrium, more particularly for
calculating equilibrium constant & enthalpies.
A realistic mathematical model with good flexibility & generalization properties.
A reliable and efficient numerical method for computing the solution.
The capability of using the method for design purposes.
Steady state modeling:The mathematical relations describing the counter current separation process like distillation
consists of four sets of basic equations, called mesh equations:
Mass balance equation(M)
Equilibrium relations(E)
Sum(or conservation) equation(S)
Enthalpy balance relations(H)
Presentation and arrangement depend on the chosen iterative variable and on the balance
formulation; so various numerical methods have been proposed. An analysis of their
advantages and disadvantages was reviewed by Wang et.al.(1981). The method can be
classified into two main categories: Partitioning and coupling methods.
Global and simultaneous procedures.
In each class the algorithms can differ from one another according to the following criteria:1. The iterative variable chosen and mathematical model used.
2. The way of arranging the equations and the procedure implemented for their solution.

3. The convergence method and the convergence promoters.

Partitioning and coupling methods: -The mesh equations are grouped either by stage or by
type. Initials estimates for some variables are required for starting the calculations. For a
given set of variables the groups of equations are solved in a predetermined order, the other
variables being kept constant. These last variables are updated according to the previous
result by using direct substitution method. This procedure is continued until a convergence is
achieved. Within this category the different approaches can be used.
1. Stage to stage methods: - In this class of methods, the procedure of Lewis and
Matheson, Thiele and Geddes, and the method of Holland are used. The mesh
equations are grouped stage to stage and solved from bottom to top.
2. Partitioned by type of equation: - in this approach temperature and flow rates must
be initialised. The mass balance equations are combined with the equilibrium
equations to form the first set of equations. These equations are linearised by keeping
the constant flow rates and equilibrium constants, and then they are solved component
by component. By using the new values of compositions, the sum and enthalpy
balance equations can be solved either simultaneously or separately in order to obtain
the new temperature & flow rates profiles.
Global and simultaneous procedures: In this class of methods, the mesh equations are linearized and simultaneously solved by
means of Newtons-Raphson Method. A Jacobian matrix containing a large number of partial
derivatives (numerically or analytically) must be generated and the equations are solved with
the initial set of variable.
Mathematical model: Consider the column shown in following figure, involves N plates numbered from top to
bottom, including a condenser (total or partial) numbered 1(one) and a reboiler numbered N.
The liquid (molar flow L ) and vapour (molar flow V ) coming from the stage j from a liquid
vapour equilibrium. On each stage j, a feed stream F and liquid (SL ) and vapour (SV ) side

stream are planned. Heat supplies or losses Q on each plate are taken into account, so nonadiabatic column can be studied.

Q 1 (Q c)
Stage 1

V2

SV2
F2

V1 (D v )

L1

Stage 2: P2T2

SL1 (DL )
Q2
SL2

VJ+1

SVJ-1
FJ-1

(J-1) PJ-1 TJ-1

QJ-1

VJ

LJ-1

SVJ
FJ

SLJ-1
(J) PJ TJ

QJ

LJ

VJ+1

SVJ+1
FJ+1

(J+1) PJ+1 TJ+1

QJ+1

SLJ+1

LJ+1

VN-1

SVN-1
FN+1

SLJ

(N-1) PN-1 TN-1


VN

SLN-1

LN-1

QN (QB)

Stage N

Fig. 26.1: Continuous Distillation Column

LN (B)
(W)

SVj

Fj

Vj

Lj-1

ij

Xi,j-1

Yi j

Ij-1

Qj

Theoretical plate j

Zi j
Vj+1

SLj

Yi, j+1
Jj+1

Lj

Fig. 26.2: Theoretical Plate of a Continuous Column


When the configuration of the column is fixed, the steady state mesh equations can be
expressed as follows for a mixture involving C components:
Mass balance equation(M) ( j=1,N)
a) Overall mass balance:
Lj-1 + Vj+1 +Fj Lj SLj Vj SVj =0

26.1

b) Partial mass balance( i=1,C)


Lj-1 .Xi, j-1 (Lj +SLj ). Xi j (Vj +SVj ).Yi j + Vj+1 .Yi, j+1+Fj .Zi, j =0

26.2

Equilibrium equations(E) ( i=1,C; j=1,N)


Yi j Ki j .Xi j =0

26.3

This equation is based on the assumption that each plate is a theoretical plate, that is
to say that liquid vapour equilibrium exists on each plate between the two phases
coming from it.
Sum equations(S)

)=0

26.4

Enthalpy balance equations (H) ( j=1, N)


Lj-1 .H j-1 + Vj+1 .H j+1+Fj .HF j (Lj +SLj) .H j (Vj +SVj ) .Hv j Qj =0

26.5

This enthalpy balance is based on the assumption of an ideal heat transfer, i.e. the vapour and
liquid stream coming from each plate have the same temperature. In addition to these
equations, and correlations are required to predict the equilibrium constants and the
enthalpies. These correlations are often complex equations describing relations between
thermodynamic properties and state variable.
Degree of freedom:
In order to define the problem, all the feed flows F feed composition Z

(i=1, C), feed

temperatures and pressure are fixed. All the liquid and vapour side stream flows SLj , SVj , as
well as the amounts of heat supplied Q (2 j N 1) and the pressure Pj are also fixed.

In the case of a total condenser V1 =0, the heat QN supplied to the reboiler is fixed and the
reflux ratio, R=L1/SL1, is also fixed when the side stream flow rate is given. In the case of a
partial condenser QN and the heat supplied to the condenser Q1 are fixed. So the unknowns of
the problem are liquid and vapour compositions Xi j , Yi j , the temperature Tj on each plate,
the liquid(L1 to LN) and vapour(V1 to VN) molar flow i.e to.

You might also like