Jan 312012
*Note: This article assumes you have rudimentary knowledge of how a simple distillation column operates.  Though much of the basics of the simple model for a distillation column will be covered below, it should not be considered as a stand alone reference.


While taking my separations and mass transfer operations course at UF, I took a particularly interesting exam where I was tasked with predicting the performance of a distillation column with two inputs.  This may be getting ahead of myself, as the problem's information was written into the proof, but the problem was as follows:

The Test Problem

Given a distillation column separating n-pentane and n-heptane, with a feed rate of 200 kgmol/hr of a 40 mol% n-pentane liquid at bubble point, a 95 mol% distillate stream, a 5 mol% bottoms stream, a 30 mol% side stream with a flow rate out equal to that of the bottoms leaving, a reflux ratio equal to twice that of the minimum, and a 50% average tray efficiency, use the McCabe-Thiele graphical method and the provided vapor-liquid equilibrium data to determine a) the flow rate and composition of all streams, b) the minimum reflux ratio, c) the number of theoretical plates required, and d) the optimum placement of the feed stream and side stream. So, the information given so far is
  • F=200 kgmol/hr
  • x_F=0.40
  • q_F=1
  • x_D=0.95
  • x_W=0.05
  • G=-W (negative sign due to opposite orientation)
  • x_G=0.30
  • q_G=1
  • R_D=2R_{min}
  • \epsilon_a=0.50
Vapor-Liquid Equilibrium Data for the n-Pentane/n-Heptane System
x y
1.000 1.000
0.867 0.984
0.594 0.925
0.398 0.836
0.254 0.701
0.145 0.521
0.059 0.271
0.000 0.000

Two Solution Methods

There are two solution methods that may be used when using the McCabe-Thiele graphical method.  The classical method uses the assumptions placed on the system explicitly to determine an analytical solution for either the stripping operating line or the enriching operating line, then the operating line between the two feeds, after which the remaining operating line is found trivially by applying the rules implicit to Cartesian geometry (two points define a single line in a plane).  The other method uses the lever rule to combine the feeds and allows for the trivial solution of a single feed column to be used.  Beyond the use of the lever rule, all that is required is on of the effluent operating lines.  After determining the combined system, all that is left is to connect the remaining points with their respective operating lines.  The basis of the proof is that the method employing the lever rule is simply a mathematical manipulation of the classical method.  As it is a simple transform rule, the two methods are analytically identical.

How to Prove the Methods are the Same

To prove the two methods are identical, the intersections of the various lines were determined and the resulting coordinate pairs were then set as ratios of like coordinates and then reduced to eliminate variables.  The resulting ratios were then inspected to determine whether or not the coordinates would be the same (a ratio of one) in all cases.  Needless to say, all inspected intersections produced ratios of one, but the proof did lend itself to generalization, which will be discussed after the proof is applied to the initial problem.

System Definition

The system is a distillation column with two feed streams located by their respective optimum placements, a distillate effluent stream, and a bottoms effluent stream.  Each respective theoretical plate in the column is designated by its index i relative to the streams crossing the boundary between the top of the column and the top tray region and is associated with liquid streams entering from above (L_{i-1} ) and exiting below (L_i ) and vapor streams entering from below (V_{i+1} ) and exiting above (V_i ).  The primary feed stream (F ) is designated as entering the column at tray n and the secondary feed stream (G ) at tray m .  The vapor stream leaving the top of the column (V_1 ) is condensed with an ideal heat exchanger with thermal duty Q_c and then splits into the distillate stream (D ) and the reflux stream (L_0 ).  The liquid stream leaving the bottom (L_N ) is sent to an ideal reboiler, with thermal duty Q_W , which splits the flow between the liquid bottoms stream (W ) and the vaporized boilup stream (V_{N+1} ).  All liquid compositions are designated as associated mole fractions x_i and all vapor compositions are similarly designated as y_i .  All side streams are also associated with their fractional quality (q_i ) which is the degree to which the stream has passed the VLE transition and is defined by the difference of the saturated vapor enthalpy and the stream enthalpy scaled by the difference of enthalpies for each respectively saturated fluid phase (eg. q_F=\frac{H_{V,F}-H_F}{H_{V,F}-H_{L,F}} ).

System Assumptions

These are intentionally written as braodly as possible as they will be used for the generalization of the proof to columns of arbitrary design.
  • Equimolar overflow (L_0=L_1=L_2=\cdots =L_i and V_1=V_2=V_3=\cdots =V_i )
  • No sub-cooling of the reflux entering at the top or super-heating of the boilup at  the bottom
  • Net positive feed into the column from side streams (\sum_{i}{F_i}\geq 0 )
  • Side streams are at bubble point (q_i=1 )

1.0 Solution Method 1: The Classical McCabe-Thiele Graphical Method

1.1.0 Mass Balances

1.1.1 Overall Mass Balance

First, the sum of the mass entering and the mass exiting must be equal to zero.  There are no reactions or other structural changes to the chemicals, therefore molar flows may be used as a proxy for the associated mass flows.  For sake of convenience, the distillate and bottoms streams are designated as negative flows into the column as well as all explicit side draw streams.  Feed streams are assumed to be positive flows into the column, but may flow in either direction so long as the column receives a net positive flow. Eq. 1-1) \displaystyle 0=F+G-D-W\leadsto F+G=D+W Next, given the assumptions above, the sum of all flows of a given chemical species must also be zero. Eq. 1-2) \displaystyle 0=Fx_F+Gx_G-Dx_D-Wx_W\leadsto Fx_F+Gx_G=Dx_D+Wx_W

1.1.2 Mass Balance on the Top (Enriching) Section

Beginning with a rearrangement of the mass balance on the top plate: Eq. 1-3 and Eq. 1-4) \displaystyle\left\{\begin{matrix}V_1+L_{n-1}=V_n+L_0 \\ V_1y_1+L_{n-1}x_{n-1}=V_ny_n+L_0x_0\end{matrix}\right. Using the mass balance across the condenser: Eq. 1-5 and Eq. 1-6) \displaystyle \hookrightarrow\left\{ \begin{matrix}\underbrace{V_1-L_0}_{D}=V_n-L_{n-1}\leadsto D=V_1-L_0 \\ V_1\underbrace{y_1}_{x_D}+L_{n-1}x_{n-1}-L_0\underbrace{x_0}_{x_D}=V_ny_n\end{matrix}\right. Eq. 1-7) \displaystyle \hookrightarrow\underbrace{\left(V_1-L_0\right)}_{D}x_D+\underbrace{L_{n-1}}_{L_0}x_{n-1}=\underbrace{V_n}_{V_1}y_n\leadsto y_n=\frac{L_0}{V_1}x_{n-1}+\frac{D}{V_1}x_D Introducing a quantity to signify the ratio of flows exiting and recycling at the top of the column: Eq. 1-8) \displaystyle \boxed{R_D=\frac{L_0}{D}} Rearranging to allow substitution and simplification: Eq. 1-9) \displaystyle y_n=\cfrac{\cfrac{1}{D}}{\cfrac{1}{D}}\frac{L_0}{L_0+D}x_{n-1}+\cfrac{\cfrac{1}{D}}{\cfrac{1}{D}}\frac{D}{L_0+D}x_D\leadsto y_n=\cfrac{\cfrac{L_0}{D}}{\cfrac{L_0}{D}+\cfrac{D}{D}}x_{n-1}+\cfrac{\cfrac{D}{D}}{\cfrac{L_0}{D}+\cfrac{D}{D}}x_D Finally, the operating line for the enriching section is produced: Eq. 1-10) \displaystyle \boxed{y=\frac{R_D}{R_D+1}x+\frac{x_D}{R_D+1}}

1.1.3 Mass Balance Across Plate n

The mass balance across plate n is given by Eq. 1-11 and Eq. 1-12) \displaystyle \left\{\begin{matrix}V_n=V_{n+1}+F\left(1-q_F\right) \\ L_n=L_{n-1}+Fq_F\end{matrix}\right.\leadsto\left\{\begin{matrix}V_1=L_0+D=V_{n+1}+F\left(1-q_F\right) \\ L_n=L_0+Fq_F\end{matrix}\right.

1.1.4 Mass Balance Between Plate m and Plate n

The mass balance between plate m and plate n is given by Eq. 1-13 and Eq. 1-14\displaystyle \left\{\begin{matrix}V_{n+1}+L_{m-1}=V_m+L_n \\ V_{n+1}y_{n+1}+L_{m-1}x_{m-1}=V_my_m+L_nx_n\end{matrix}\right.

1.1.5 Mass Balance Across Plate m

The mass balance across plate m is given by Eq. 1-15 and Eq. 1-16) \displaystyle \left\{\begin{matrix}V_m=V_{m+1}+G\left(1-q_G\right) \\ L_m=L_{m-1}+Gq_G\end{matrix}\right.\leadsto\left\{\begin{matrix}V_m=V_{N+1}+G\left(1-q_G\right) \\ L_N=L_{m-1}+Gq_G\end{matrix}\right.

1.1.6 Mass Balance on Bottom (Stripping) Section

The mass balance across the bottom (stripping) section is given by Eq. 1-17 and Eq. 1-18) \displaystyle \left\{\begin{matrix}V_{m+1}+L_N=V_{N+1}+L_m \\ V_{m+1}y_{m+1}+L_Nx_N=V_{N+1}y_{N+1}+L_mx_m\end{matrix}\right.\leadsto\left\{\begin{matrix}V_{m+1}-L_m=\boxed{V_{N+1}-L_N=W} \\ V_{m+1}y_{m+1}=L_mx_m+\underbrace{V_{N+1}y_{N+1}-L_Nx_N}_{Wx_W}\end{matrix}\right.

1.2.0 Operation Lines and q-Lines

1.2.1 Enriching Operation Line

This is the result, Eq. 1-10, derived in section 1.1.2: Eq. 1-19) \displaystyle \boxed{y=\frac{R_D}{R_D+1}x+\frac{x_D}{R_D+1}}

1.2.2 Stripping Operation Line (Incomplete Form)

Starting from Eq. 1-18 from section 1.1.6: Eq. 1-20) \displaystyle V_{m+1}y_{m+1}=L_mx_m+Wx_W\leadsto y_{m+1}=\frac{L_N}{V_{m+1}}x_m+\frac{W}{V_{m+1}}x_W Including from section 1.1.5: Eq. 1-21 and Eq. 1-22) \displaystyle \left\{\begin{matrix}L_m=L_N \\ V_{m+1}=V_{N+1}=L_N+W\end{matrix}\right. And then substituting yields: Eq. 1-23) \displaystyle \boxed{y_{m+1}=\frac{L_N}{L_N+W}x_m+\frac{W}{L_N+W}x_W} \,

1.2.3 Operation Line Between Stream F and Stream G

Starting with Eq. 1-13 and Eq. 1-14 from section 1.1.4: Eq. 1-24 and Eq. 1-25) \displaystyle\left\{\begin{matrix}V_{n+1}+L_{m-1}=V_m+L_n \\ V_{n+1}y_{n+1}+L_{m-1}x_{m-1}=V_my_m+L_nx_n\end{matrix}\right.\leadsto\left\{\begin{matrix}V_{n+1}-L_n=V_m-L_{m-1}=\text{const.} \\ V_{n+1}y_{n+1}=L_nx_n+V_my_m-L_{m-1}x_{m-1}\end{matrix}\right. Which rearranges to yield the incomplete form: Eq. 1-26) \displaystyle\hookrightarrow\boxed{y_{n+1}=\frac{L_n}{V_{n+1}}x_n+\frac{V_my_m-L_{m-1}x_{m-1}}{V_{n+1}}} Using substitutions from Eq. 1-11 and Eq. 1-12 in section 1.1.3 and Eq. 1-15 and Eq. 1-16 from section 1.1.5 yields: Eq. 1-27) \displaystyle y_{n+1}=\frac{L_0+Fq_F}{L_0+D-F\left(1-q_F\right)}x_n+\frac{\left(V_{N+1}+G\left(1-q_G\right)\right)y_m-\left(L_N-Gq_G\right)x_{m-1}}{L_0+D-F\left(1-q_F\right)} Further substitution using Eq. 1-17 from section 1.1.6 in the numerator of the second term and subsequent rearrangement of said numerator yields: Eq. 1-28) \displaystyle y_{n+1}=\frac{L_0+Fq_F}{L_0+D-F\left(1-q_F\right)}x_n+\frac{\left(L_N-Gq_G\right)\left(y_m-x_{m-1}\right)+\left(W+G\right)y_m}{L_0+D-F\left(1-q_F\right)} Combining Eq. 1-14, Eq. 1-16, Eq. 1-17, and Eq. 1-18 and rearranging yields the relationship Eq. 1-29) \displaystyle Wx_W-Gx_G=L_{m-1}x_{m-1}-V_my_m Substitution using Eq. 1-29 into Eq. 1-26 yields Eq. 1-30) \displaystyle y_{n+1}=\frac{L_0+Fq_F}{L_0+D-F\left(1-q_F\right)}x_n+\frac{Gx_G-Wx_W}{L_0+D-F\left(1-q_F\right)} Which rearranges to yield Eq. 1-31) \displaystyle\hookrightarrow y_{n+1}=\cfrac{\cfrac{L_0}{D}+\cfrac{Fq_F}{D}}{\cfrac{L_0}{D}+\cfrac{D}{D}-\cfrac{F\left(1-q_F\right)}{D}}x_n+\cfrac{\cfrac{Gx_G}{D}-\cfrac{Wx_W}{D}}{\cfrac{L_0}{D}+\cfrac{D}{D}-\cfrac{F\left(1-q_F\right)}{D}} \displaystyle\leadsto y_{n+1}=\cfrac{R_D+\cfrac{F}{D}q_F}{R_D+1-\cfrac{F}{D}\left(1-q_F\right)}x_n+\cfrac{Gx_G-Wx_W}{\left(R_D+1-\cfrac{F}{D}\left(1-q_F\right)\right)D} Which finally yields the completed operating line Eq. 1-32) \displaystyle\boxed{y=\cfrac{R_D+\cfrac{F}{D}q_F}{R_D+\cfrac{F}{D}\left(q_F-1\right)+1}x+\cfrac{Gx_G-Wx_W}{D\left(R_D+\cfrac{F}{D}\left(q_F-1\right)+1\right)}}

1.2.4 q-Line for Stream F

Given by definition the q-line for stream F is Eq. 1-33) \displaystyle\boxed{y=\frac{q_F}{q_F-1}x-\frac{x_F}{q_F-1}} Combining Eq. 1-10 in section 1.1.2 and Eq. 1-28 in section 1.2.3 and rearranging then canceling out L_0 , D , and F yields Eq. 1-34) \displaystyle y=\frac{q_F}{q_F-1}x+\frac{Dx_D-\left(L_N-Gq_G\right)\left(y_m-x_{m-1}\right)-\left(W+G\right)y_m}{F\left(1-q_F\right)} Which may be equated with Eq. 1-33 and rearranged to yield Eq. 1-35) \displaystyle Fx_F=Dx_D-\left(L_N-Gq_G\right)\left(y_m-x_{m-1}\right)-\left(W+G\right)y_m Combining this with Eq. 1-1 and Eq. 1-2 from section 1.1.1 and rearranging yields Eq. 1-36) \displaystyle Wx_W-Gx_G=\left(L_N-Gq_G\right)x_{m-1}-\left(V_{N+1}+G\left(1-q_G\right)\right)y_m Which also rearranges to yield Eq. 1-29 from section 1.2.3.

1.2.5 q-Line for Stream G

The q-line for stream G is trivially given, as was seen with stream F , by Eq. 1-37) \displaystyle\boxed{y=\frac{q_G}{q_G-1}x-\frac{x_G}{q_G-1}}

 1.2.6 Point of Intersection of Enriching Operation Line and Operation Line Between Stream F and Stream G

Equating Eq. 1-10 from section 1.1.2 and Eq. 1-32 from section 1.2.3 at the intersection \left\{x_{FD},y_{FD}\right\} yields Eq. 1-38) \displaystyle\cfrac{R_D}{R_D+1}x_{FD}+\cfrac{x_D}{R_D+1}=\cfrac{R_D+\cfrac{F}{D}q_F}{R_D+\cfrac{F}{D}\left(q_F-1\right)+1}x_{FD}+\cfrac{Gx_G-Wx_W}{D\left(R_D+\cfrac{F}{D}\left(q_F-1\right)+1\right)} Solving for x_{FD} and simplifying yields Eq. 1-39) \displaystyle\boxed{x_{FD}=\frac{x_D\left(F\left(q_F-1\right)+D\left(R_D+1\right)\right)+\left(R_D+1\right)\left(Wx_W-Gx_G\right)}{F\left(q_F+R_D\right)}} Plugging Eq. 1-39 back into Eq. 1-10 and simplifying yields Eq. 1-40) \displaystyle\boxed{y_{FD}=\frac{x_D\left(Fq_F+DR_D\right)+R_D\left(Wx_W-Gx_G\right)}{F\left(q_F+R_D\right)}}

1.2.7 Point of Intersection of Operation Line Between Stream F and Stream G and q-Line  for Stream G

Equating Eq. 1-32 from section 1.2.3 and Eq. 1-37 from section 1.2.5 at the intersection \left\{x_{GW},y_{GW}\right\} yields Eq. 1-41) \displaystyle\cfrac{R_D+\cfrac{F}{D}q_F}{R_D+\cfrac{F}{D}\left(q_F-1\right)+1}x+\cfrac{Gx_G-Wx_W}{D\left(R_D+\cfrac{F}{D}\left(q_F-1\right)+1\right)}=\cfrac{q_G}{q_G-1}x-\cfrac{x_G}{q_G-1} Solving for x_{GW} and simplifying yields Eq. 1-42) \displaystyle\boxed{x_{GW}=\frac{x_G\left(F\left(q_F-1\right)+D\left(R_D+1\right)\right)+\left(q_G-1\right)\left(Gx_G-wx_W\right)}{F\left(q_F-q_G\right)+D\left(q_G+R_D\right)}} Plugging Eq. 1-42 back into Eq. 1-37 and simplifying yields Eq. 1-43) \displaystyle\boxed{y_{GW}=\frac{x_G\left(Fq_F+Gq_G\right)+DR_Dx_G-q_GWx_W}{F\left(q_F-q_G\right)+D\left(q_G+R_D\right)}}

1.2.8 Stripping Operation Line (Complete Form)

The final operation line may be constructed from the known points at either endpoint, \left\{x_W,y_W\right\} and \left\{x_{GW},y_{GW}\right\} . The operation line is constructed from the rearranged point-slope form, given by Eq. 1-44) \displaystyle\frac{y_{GW}-x_W}{x_{GW}-x_W}\left(x-x_W\right)+x_W Notice that y_W has been replaced by x_W as they are the same value. Substituting into Eq. 1-44 using Eq. 1-42 and Eq. 1-43 from section 1.2.7 and rearranging yields the complete form of the stripping operation line Eq. 1-45) \displaystyle\boxed{y=\frac{\left(Fq_F+Gq_G+DR_D\right)x_Gx-\left(\left(F\left(q_F-q_G\right)+D\left(q_G+R_D\right)+q_GW\right)x+Wx_G\right)x_W+Wx^2_W}{\left(F\left(q_F-1\right)+G\left(q_G-1\right)+D\left(R_D+1\right)\right)x_G-\left(F\left(q_F-q_G\right)+D\left(q_G+R_D\right)+W\left(q_G-1\right)\right)x_W}}

1.3.0 Solution of the Test Problem

The original problem can now be solved using the data provided and the equations derived.

1.3.1 Plot the VLE Data

First, the equilibrium data provided is plotted.

1.3.2 Plot the VLE Data Scaled for Efficiency

Then the equilibrium line is adjusted to reflect the column's efficiency.

1.3.3 Solve the Overall Mass Balance and Plot Composition Lines

Next, the mass balance is solved to determine the compositions and flow rates of all streams.  Then lines denoting the respective compositions of the individual streams are placed.

1.3.4 Plot q-Lines for All Streams

q-lines are plotted for all streams.  Since the bottoms and distillate streams leave as liquids at bubble point, their respective q-lines are vertical.

1.3.5 Determine the Minimum Reflux Ratio

R_{min} is determined by finding the slope of the line from the distillate composition point to the first pinch point along the bounding VLE curve before or at the feed stream or side stream of interest and then using the slope from Eq. 1-10 in section 1.1.2.

1.3.6 Plot the Operation Lines

Use the equations derived above (Eq. 1-10, Eq. 1-32, and Eq. 1-45) to plot the operation lines.

1.3.7 Plot the Plate Locations

To complete the problem, the plate lines are drawn between the bounding VLE curve and the contiguous operation lines, starting at the distillate composition point and ending at the tray crossing the bottoms composition boundary.  The number of triangles represents the number of trays required and the triangles containing the composition points for the feed streams and side streams are the optimal placement of the respective streams.

1.4.0 Automation Using Mathematica

The plots above were generated using some of the results of this proof as derived in a later section.  This may be jumping the gun a little, but the code for this section's plots is given below.  Reading through the remainder of the proof is recommended before trying to understand the details of the math behind the plot generation. That's it for now. The rest of the proof will be posted in chunks of about the same size as this and the complete proof will be assembled in a single post when all of the parts are completed and reviewed. Thank you for reading my work and please do leave a comment to let me know what you think of it.
Print Friendly
Be Sociable, Share!
Short URL: http://dlvr.it/18Lt90

  4 Responses to “Distillation of a Binary Mixture in a Distillation Column of Arbitrary Design, Part 1”

  1. […] This article is part of a series on a specific research project: Part 1 *Note: This article assumes you have rudimentary knowledge of how a simple distillation column […]

  2. […] This article is part of a series on a specific research project: Part 1, Part 2 *Note: This article assumes you have rudimentary knowledge of how a simple distillation […]

  3. […] This article is part of a series on a specific research project: Part 1, Part 2, Part 3 *Note: This article assumes you have rudimentary knowledge of how a simple […]

  4. Awesome work!

    Do you have anything about distillation design for cyclohexane and n-hexane?


Sorry, the comment form is closed at this time.