Skip to main content
Engineering LibreTexts

9: Putting It All Together

  • Page ID
    8147
  • \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

    ( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\id}{\mathrm{id}}\)

    \( \newcommand{\Span}{\mathrm{span}}\)

    \( \newcommand{\kernel}{\mathrm{null}\,}\)

    \( \newcommand{\range}{\mathrm{range}\,}\)

    \( \newcommand{\RealPart}{\mathrm{Re}}\)

    \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

    \( \newcommand{\Argument}{\mathrm{Arg}}\)

    \( \newcommand{\norm}[1]{\| #1 \|}\)

    \( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

    \( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

    \( \newcommand{\vectorA}[1]{\vec{#1}}      % arrow\)

    \( \newcommand{\vectorAt}[1]{\vec{\text{#1}}}      % arrow\)

    \( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vectorC}[1]{\textbf{#1}} \)

    \( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

    \( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

    \( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

    \( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

    \( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

    \(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)

    Distillation Science (a blend of Chemistry and Chemical Engineering)

    This is Part IX, Putting It All Together of a ten-part series of technical articles on Distillation Science, as is currently practiced on an industrial level. See also Part I, Overview for introductory comments, the scope of the article series, and nomenclature.

    Part IX, Putting It All Together uses the information from Parts III through VII, showing how to combine them in a practical example for bulk separation; and how using the techniques detailed in previous articles give answers that differ from the ideals of basic Raoult/Dalton Law application.

    See previous articles per the following detail:

    • Part III for critical properties and acentric factor, especially Table 3-3
    • Part IV for the recommended vapor pressure equation, especially Equation 4-3 and Table 4-1
    • Part V for Equation of State, especially Equations 5-3 and 5-4
    • Part VI, especially Equation 6-4 for partial pressure and Equations 6-7 and 6-8 for fugacity coefficients
    • Part VII for Liquid Activity Coefficients, especially Equation 7-3

    In this article a practical distillation example is developed, showing the Y/X relationship with temperature for the DCS-TCS binary system with pressure is held constant at 11.0 atmospheres, and showing the results as a TXY/P plot as well as in tabular form. From such data, a distillation column designer would consider the best arrangement for column feed, “tops” and “bottoms” recovery, tray count and external reflux ratio ( and from that, the energy required for the reboiler and condenser, as well as hydraulic loading per unit mass flow of feedstock). Two variations of this example are worked-up, showing the application of distillation science: one which only uses ideal Raoult/Dalton, and one that uses the full set of ideality departure coefficients. Table output and plots below are done using MS Excel. Then at the end of the article, the table results are graphically demonstrated using a McCabe-Theile plot to step off the trays for an appropriate reflux rate and show the best feed tray.

    The table of Figure 9-1 is set up to cover the range of DCS liquid mole fraction from unity to zero, in increments of 0.05 mole fraction, with 0.01/0.99 mole fraction added to illustrate asymptotic effects. For each non-zero table row, a temperature is first assumed, and then the vapor mole fractions calculated per Raoult/Dalton, with p=11.0 atmospheres; then that row’s temperature is iteratively adjusted until the total pressure = 11.0 atmospheres.

    In the plot of Figure 9-1 shows the calculation results in TXY/P form, with the value of XDCS=1.0 representing the tightest possible “tops” product at the distillation column condenser; and a value of XDCS=0.0 representing the tightest possible “bottoms” product at column reboiler, all at a total reflux condition. For a finite column, the mole fraction of the “tops” product would be a cropping of the upper XDCS values of the TXY/P curves; and the mole fraction of the “bottoms” product a would be a cropping of the lesser XDCS values of the TXY/P curves. In such instance, the column feed tray would be identified as that point in the column where the tray temperature is identical to the feed mixture’s boiling point. The accompanying plot of solution shows what the TXY/P curves would look like, for such a Raoult/Dalton solution to the DCS-TCS binary at 11.0 atmospheres.

    \( Y_{1}=(X_{1} \times VP_{1})/(X_{1} \times VP_{1} + X_{2} \times VP_{2}) \nonumber\) \( Y_{2}=(X_{2} \times VP_{2})/(X_{1} \times VP_{1} + X_{2} \times VP_{2}) \nonumber\) repeat Equation 6-1

    Screen Shot 2019-03-06 at 10.22.24 PM.png
    Figure 9-1, TXY calculation table and plot, Raoult/Dalton

    If the data plot of Figure 9-1 plot is carefully examined, one can see that the T-X and T-Y curves are exactly symmetrical when “folded” along a line running from abscissa/ordinate points of (0,124.76) to (1,92.40); or “folded” orthogonally to that same line. In other words, there are no departures from ideality in either vapor or liquid phase, either with DCS or TCS.

    Now X and Y are re-calculated including the departures from ideality, both in fugacity (ø) as well as Liquid Activity Coefficients ( \(\gamma \nonumber\) ).

    Screen Shot 2019-03-06 at 10.23.22 PM.png
    Figure 9-2, TXY/P calculation table, all departure factors included

    The table in above Figure 9-2 is likewise set up and calculations performed row-wise for each value of XDCS, but the liquid and vapor-phase fugacity coefficients (øL and øV) and activity coefficients ( \(\gamma_{i} \nonumber\) ) are included in the calculation of vapor mole fractions (\(\gamma_{DCS} \nonumber\) ). For each row, the temperature is iterated until the sum of partial pressures = 11.0 atmospheres. Note that since the activity coefficients are mild functions of temperature, each row has its own Van Laar constants determined and the value of \(\gamma \nonumber\) calculated for that row’s liquid mole fraction. Thus all departure-from-ideality factors are included.

    \[ Y_{i}= (\phi_{i}^L \times \gamma_{i} \times VP_{i} \times X_{i})/(\phi_{i}^V \times \pi) \label{9-1} \]

    which is an algebraic re-arrangement of Equation 6-4 from previous Part VI.

    To better show how the several ideality departure factors impact the TXY/P plot (as compared to Figure 9‑1 above), the table in Figure 9-2 is plotted. The symmetry of the Figure 9-1 Raoult/Dalton plot is no longer there, neither along the “fold” line running from abscissa/ordinate points of (0,124.76) to (1, 92.40); or “folded” orthogonally to that same line. Also note that the T-X and T-Y curves are closer together in Figure 9‑3, as opposed to Figure 9-1.

    Screen Shot 2019-03-06 at 10.25.39 PM.png
    Figure 9-3, TXY/P plot, all departure factors included

    For a very simple approximation of tray count, for given distillation column “tops” product , “bottoms” product and feed mixture mole fractions, and for a given reflux rate, a McCabe-Thiele plot is constructed in Figure 9-4 to graphically step off the number of trays above and below the feed tray. For almost all industrial applications, this tray-by-tray analysis is done via computer modeling, such as ASPEN, VMG, HYSYS, ChemCad, etc. To apply the technique to a computer model, the user would specify the critical constants ( as per Part III), the vapor pressure equation to be used (as per Part IV), the EOS to be used to calculate fugacity coefficients (as per Part V), and the Liquid Activity Coefficient model to be used (as per Part VI).

    To construct a McCabe-Thiele plot, the X-Y equilibrium (in red below) is plotted on a set of axes where X= liquid mole fraction of the most volatile component = XDCS ; and Y= vapor mole fraction of the most volatile component = YDCS . The total reflux line (in black below) connects (X,Y) = (0,0) and (1,1); representing the maximum reflux possible. The lines are added (in dashed blue below) for the feed (XF), and distillate product (XD) and the bottoms product (XB), from the X-axis to the total reflux line. The feed line is extended toward the equilibrium curve, and from that intersection a line is drawn to the end of the distillate product line. The slope of that drawn line represents the minimum internal reflux needed to make the separation (i.e., with an infinite number of trays). The upper column operating line (in magenta below) is then drawn connecting the end of the distillate product line with the feed line, with a slope somewhat greater than the minimum reflux slope. In normal practice the upper column operating line slope is usually about 1.3 times the minimum slope, depending on the relative value of operating energy vs capital cost. Then the lower column operating line is drawn, intersecting the bottoms product line with the upper column operating line (and the feed line). The number of theoretical trays (aka stages) can then be stepped off ( in green below) from XD to XB, and best tray determined for the feed.

    Screen Shot 2019-03-06 at 10.26.21 PM.png
    Figure 9-4, McCabe Thiele solution to example calculations

    To illustrate an example of using a McCabe-Thiele plot, assume the mixture to be separated is a 35% DCS molar liquid, pre-heated to its boiling point of 108.61°C (as shown on the 35% of Figure 9-2). Further, assume that the column “tops” distillate product is to be 98% molar DCS and the column “bottoms” product is to be 98% molar TCS. Per the text above, the minimum reflux (i.e., the slope of the line between XD = XDCS of 0.98, XB = XDCS of 0.02) is determined to be 3.9:1, reflux-to-distillate. To allow a reasonable number of trays, a reflux rate of 5:1 is chosen (about 1.3 times the minimum), so the slope of the upper column operating line is R/(R+1) = 5/(5+1) =0.8333. Stepping off trays, it is determined that the total number of theoretical trays required is 34, with the feed tray at #22. Conventional tray counting starts at the column condenser = tray #1.

    As can be seen by examining Figure 9-4, as the column’s reflux rate is increased, upper column (magenta) operating line’s slope increases, and the gap opens up between the (red) X-Y Equilibrium Line and the (black) Total Reflux Line. So increasing the reflux decreases the number trays required to make the specified separation.

    In a normal industrial process application, all of the calculations are done via computer simulators, allowing for variations in operating line slope; feed quality condition (e.g., saturated vapor, saturated liquid, sub-cooled liquid, etc); changes in the internal liquid/vapor flow ratios due to physical properties; ambient heat effects; tray efficiencies, etc; as well as integrating the distillation column into the rest of the process design. So the McCabe-Thiele method is just used for graphical explanation.

    However, the primary consideration of the distillation column design is the X-Y equilibrium relationship. If this is considerably in error – such as using only Raoult/Dalton – the column design will be a failure. Given the cost of industrial-scale distillation systems being in the millions of dollars (tens of $MM for larger ones), there is every reason to get the X-Y equilibrium relationship correct.

    In the example of this article, the basic Raoult/Dalton model would result in the X-Y equilibrium relationship in the above example having more curvature (than shown in Figure 9-4) and therefore more open gap between the equilibrium curve and operating lines. The resulting reflux determined would be about 2/3 of required, and number of trays required would also be about 2/3 of needs - so a major error in design. Given what is at stake financially, there is no question that the use of proper distillation science is well worth it, even if the concepts are complex and the calculations are tedious.


    This page titled 9: Putting It All Together is shared under a CC BY 4.0 license and was authored, remixed, and/or curated by Larry Coleman.

    • Was this article helpful?