Feb 292012
 February 29, 2012  Posted by at 4:00 pm Chemical Engineering, Education, Projects Tagged with: , , , , , , , ,

*Note: 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 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.

4.0 Generalization to Columns of Arbitrary Design

4.1.0 Motivation for Generalizing to an Arbitrary Design

The motivation to generalize to a system with an arbitrary topology is simply to take the next logical step and expand on what has already been accomplished. The case for two feed/side streams was itself an extension of the simple case of a single feed stream.

4.2.0 Physical Description for a Generalized Column Model

To generalize to an arbitrary topology, some choices as far as the idealized general topology have to be made. This leads to the least complex yet most regularly structured limit of including a single feed stream and a single side stream placed optimally at each theoretical plate. So, the column is now comprised of \displaystyle n plates, each of which is associated with a feed stream (\displaystyle F_i , \displaystyle x_{F_i} , and \displaystyle q_{F_i} ), a side draw stream (\displaystyle S_i , \displaystyle x_{S_i} , and \displaystyle q_{S_i} ), a liquid overflow entering from above (\displaystyle L_{i-1} and \displaystyle x_{i-1} ) and underflow leaving below (\displaystyle L_i and \displaystyle x_i ), a vapor overflow leaving above (\displaystyle V_i and \displaystyle y_i ) and underflow entering from below (\displaystyle V_{i+1} and \displaystyle y_{i+1} ), and the distillate and waste streams with their associated values as described in part 1.

4.3.0 Method of Model Generalization

4.3.1 Use and Extension of the Base Model

Now that there is ample precedent for the theoretical extension to multiple feed and side streams, as demonstrated in part 1, the same extension may be made utilizing the same framing of the transport balances about the plate where the feed enters from the base model. This will be further extended initially by the addition of a side draw stream to the mass balance.

4.3.2 Refactoring of the Approach

The prior method of allowing the balances about a side attachment to be located at an arbitrary position will be left alone as it is the basis for the balances being derived. However, the exact definition of the flow rates between adjacent plates will be refactored to include the contributions of the preceding tray’s feed and side streams as well as the compounded contributions before it on either side. This leads to the slight complication of equimolar overflow being a relevant assumption, but only between locations where there is fluid entering or exiting the column. Therefore, this fundamental assumption is no longer a strictly binding assumption as discrete differentials in holdup between plates may later be rectified using quasi-streams in place of the newly provided feed and side streams.

The other refactoring that will be done is the removal of the assumption that the streams crossing the boundary of the overall system are at saturated liquids at bubble point. This is a relatively minor change that will likely be reinstated for example purposes as was effectively done during the derivation for the original case (see part 1 and part 2).

4.3.3 Inherent Assumptions of the Derivation

In light of the above restructuring and enhancement of the forthcoming derivation’s approach and the prior assumptions listed in part 1, the revised list of assumptions is somewhat simpler.

  • Equimolar overflow across individual plates
  • 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 across all trays (\displaystyle \sum_{i}{F_i-S_i}\geq 0 )

4.4.0 Derivation of the Generalized Model

4.4.1 Overall Mass Balance

Beginning with the overall mass balance at steady state and rearranging to indicate the idealized direction of material flows yields

Eq. 4-1 and Eq. 4-2) \displaystyle\left\{\begin{matrix}\displaystyle\sum^n_{i=1}{F_i}&=&D+W+\displaystyle\sum^n_{i=1}{S_i} \\ \displaystyle\sum^n_{i=1}{F_ix_{F_i}}&=&Dx_D+Wx_W+\displaystyle\sum^n_{i=1}{S_ix_{S_i}}\end{matrix}\right.

4.4.2 Mass Balance for an Arbitrary Plate

For an arbitrary plate, the equimolar flow provision is modified to include contributions from the feed and side streams as

Eq. 4-3 and Eq. 4-4) \displaystyle\left\{\begin{matrix}L_{i-1}+F_iq_{F_i}&=&L_i+S_iq_{S_i} \\ V_{i+1}+\left(1-q_{F_i}\right)F_i&=&V_i+\left(1-q_{S_i}\right)S_i\end{matrix}\right.

Combining these two equations yields

Eq. 4-5) \displaystyle L_{i-1}+F_iq_{F_i}+V_{i+1}+\left(1-q_{F_i}\right)F_i=L_i+S_iq_{S_i}+V_i+\left(1-q_{S_i}\right)S_i

Which serves as the basis for the mass balances across an arbitrary plate

Eq. 4-6 and Eq. 4-7) \displaystyle\left\{\begin{matrix}L_{i-1}+V_{i+1}+F_i&=&L_i+V_i+S_i \\ L_{i-1}x_{i-1}+V_{i+1}y_{i+1}+F_ix_{F_i}&=&L_ix_i+V_iy_i+S_ix_{S_i}\end{matrix}\right.

Rearranging Eq. 4-5 in light of the overall mass balance given by Eq. 4-1 yields

Eq. 4-8) \displaystyle\begin{matrix}D&=&V_1&-&L_0 \\ &=&\left(V_{i+1}+\left(1-q_{F_i}\right)F_i\right)&-&\left(L_i+S_iq_{S_i}\right) \\ &=&\left(V_i+\left(1-q_{S_i}\right)S_i\right)&-&\left(L_{i-1}+F_iq_{F_i}\right)\end{matrix}

4.4.3 Descending Liquid Flow for an Arbitrary Plate

Starting with Eq. 4-3 and rearranging yields

Eq. 4-9) \displaystyle L_i=L_{i-1}+F_iq_{F_i}-S_iq_{S_i}

Which, when shifted back one tray and substituted back into itself, yields

Eq. 4-10) \displaystyle L_i=L_{i-2}+F_{i-1}q_{F_{i-1}}-S_{i-1}q_{S_{i-1}}+F_iq_{F_i}-S_iq_{S_i}

Following this pattern through the top plate yields the liquid underflow of any given plate

Eq. 4-11) \displaystyle\boxed{L_i=L_0+\sum^i_{j=1}{\left(F_jq_{F_j}-S_jq_{S_j}\right)}}

4.4.4 Ascending Vapor Flow for an Arbitrary Plate

Starting with Eq. 4-4 and rearranging yields

Eq. 4-12) \displaystyle V_{i+1}=V_i+\left(1-q_{S_i}\right)S_i-\left(1-q_{F_i}\right)F_i

Which, when shifted back one tray and substituted back into itself, yields

Eq. 4-13) \displaystyle V_{i+1}=V_{i-1}+\left(1-q_{S_{i-1}}\right)S_{i-1}-\left(1-q_{F_{i-1}}\right)F_{i-1}+\left(1-q_{S_i}\right)S_i-\left(1-q_{F_i}\right)F_i

Following this procedure one further time to demonstrate the established pattern yields

Eq. 4-14) \displaystyle V_{i+1}=V_{i-2}+\left(1-q_{S_{i-2}}\right)S_{i-2}-\left(1-q_{F_{i-2}}\right)F_{i-2}+\left(1-q_{S_{i-1}}\right)S_{i-1}-\left(1-q_{F_{i-1}}\right)F_{i-1}+\left(1-q_{S_i}\right)S_i-\left(1-q_{F_i}\right)F_i

Following this patterns through the top plate yields the vapor overflow of any given plate

Eq. 4-15) \displaystyle V_{i+1}=\underbrace{V_1}_{D+L_0}+\sum^i_{j=1}{\left(\left(1-q_{S_j}\right)S_j-\left(1-q_{F_j}\right)F_j\right)}

Finally, substitution using Eq. 1-8 yields the useful form

Eq. 4-16) \displaystyle\boxed{V_{i+1}=D\left(1+R_D\right)+\sum^i_{j=1}{\left(\left(1-q_{S_j}\right)S_j-\left(1-q_{F_j}\right)F_j\right)}}

4.4.5 Operation Line for an Arbitrary Plate

Beginning with Eq. 4-7 and solving for the vapor underflow composition yields

Eq. 4-17) \displaystyle y_{i+1}=\frac{L_i}{V_{i+1}}x_i+\frac{V_iy_i+S_ix_{S_i}-\left(L_{i-1}x{i-1}+F_ix_{F_i}\right)}{V_{i+1}}

Substitution using Eq. 4-6 and rearranging yields

Eq. 4-18) \displaystyle y_{i+1}=\frac{L_ix_i+V_iy_i+S_ix_{S_i}-L_{i-1}x_{i-1}-F_ix_{F_i}}{L_i+V_i+S_i-L_{i-1}-F_i}

This form is useful in understanding the context of the final form.

Reversing the substitution in the denominator yields

Eq. 4-19) \displaystyle y_{i+1}=\frac{L_ix_i+V_iy_i+S_ix_{S_i}-L_{i-1}x_{i-1}-F_ix_{F_i}}{V_{i+1}}

Rearranging to segregate the input and output flow differentials on either side of the equation yields

Eq. 4-20) \displaystyle V_iy_i-L_{i-1}x_{i-1}=V_{i+1}y_{i+1}-L_ix_i+F_ix_{F_i}-S_ix_{S_i}

Shifting this equation one plate toward the bottom of the column and substituting back in on the right side of the equation yields

Eq. 4-21) \displaystyle V_iy_i-L_{i-1}x_{i-1}=V_{i+2}y_{i+2}-L_{i+1}x_{i+1}+F_{i+1}x_{F_{i+1}}-S_{i+1}x_{S_{i+1}}+F_ix_{F_i}-S_ix_{S_i}

Following this pattern through the bottom of the column yields

Eq. 4-22) \displaystyle V_iy_i-L_{i-1}x_{i-1}=\underbrace{V_{n+1}y_{n+1}-L_nx_n}_{Wx_W}+\sum^{n}_{j=i}{F_jx_{F_j}-S_jx_{S_j}}

And substituting with Eq. 4-2 yields

Eq. 4-23) \displaystyle V_iy_i-L_{i-1}x_{i-1}=Wx_W+\sum^{n}_{j=i}{F_jx_{F_j}-S_jx_{S_j}}

Substitution using a shift of Eq. 4-20 and both Eq. 4-11 and Eq. 4-16 and again solving for the vapor underflow composition yields the now familiar form for an operation line

Eq. 4-24) \displaystyle y_{i+1}=\frac{L_0+\displaystyle\sum^{i}_{j=1}{\left(F_jq_{F_j}-S_jq_{S_j}\right)}}{D\left(1+R_D\right)+\displaystyle\sum^{i}_{j=1}{\left(\left(1-q_{S_j}\right)S_j-\left(1-q_{F_j}\right)F_j\right)}}x_i+\frac{Wx_W+\displaystyle\sum^{n}_{j=i}{F_jx_{F_j}-S_jx_{S_j}}}{D\left(1+R_D\right)+\displaystyle\sum^{i}_{j=1}{\left(\left(1-q_{S_j}\right)S_j-\left(1-q_{F_j}\right)F_j\right)}}

Substitution using Eq. 1-8 yields the final form

Eq. 4-25) \displaystyle\boxed{y_{i+1}=\cfrac{R_D+\displaystyle\sum^{i}_{j=1}{\left(\cfrac{F_j}{D}q_{F_j}-\cfrac{S_j}{D}q_{S_j}\right)}}{1+R_D+\displaystyle\sum^{i}_{j=1}{\left(\left(1-q_{S_j}\right)\cfrac{S_j}{D}-\left(1-q_{F_j}\right)\cfrac{F_j}{D}\right)}}x_i+\cfrac{\cfrac{W}{D}x_W+\displaystyle\sum^{n}_{j=i}{\cfrac{F_j}{D}x_{F_j}-\cfrac{S_j}{D}x_{S_j}}}{1+R_D+\displaystyle\sum^{i}_{j=1}{\left(\left(1-q_{S_j}\right)\cfrac{S_j}{D}-\left(1-q_{F_j}\right)\cfrac{F_j}{D}\right)}}}

4.4.6 q-Line for an Arbitrary Plate

Starting with the general form for a q-line introduced in part 1 by Eq. 1-33 and Eq. 1-37 and using the combined streams concept derived in part 2 yields the condensed form of the q-line for an arbitrary plate

Eq. 4-26) \displaystyle y_{q_i}=\frac{q_{M_i}}{q_{M_i}-1}x_{q_i}-\frac{x_{M_i}}{q_{M_i}-1}

Working backwards from this, the combined composition is defined as

Eq. 4-27) \displaystyle x_{M_i}=\frac{F_ix_{F_i}}{M_i}-\frac{S_ix_{S_i}}{M_i}

The combined fractional quality is given by

Eq. 4-28) \displaystyle x_{M_i}=\frac{F_iq_{F_i}}{M_i}-\frac{S_iq_{S_i}}{M_i}

And the combined flow rate is

Eq. 4-29) \displaystyle M_i=F_i-S_i

Substituting for these combined variables and rearranging yields

Eq. 4-30) \displaystyle\boxed{y_{q_i}=\frac{F_iq_{F_i}-S_iq_{S_i}}{\left({1-q_{S_i}}\right)S_i-\left({1-q_{F_i}}\right)F_i}x_{q_i}-\frac{F_ix_{F_i}-S_ix_{S_i}}{\left({1-q_{S_i}}\right)S_i-\left({1-q_{F_i}}\right)F_i}}

4.5.0 A Simple Example for the Generalized Model

In fact, the best simple example of this generalization in action has already been done. In part 1 the original test problem was solved using the legacy technique of McCabe and Thiele. At the time, I did show the source code for the diagram generation, but cautioned against attempting to understand the math of it without further reading. Going back to this source code with the new generalization in hand should now shed new light on the math involved. There are a few modifications, but they are mostly simplifications rather than additions to the overall form.

4.6.0 Basic Qualitative Properties and Scope of the Generalized Model

This model comes with a few conveniences and a few annoyances, as most generalizations do. First, it is best used for optimal placement. Mal-placed feed modeling will be covered in a later post, but suffice to say, the equations work the exact same way for the general model, but the method for placement is much different on the VLE XY-plot. Second, the explicit designation of independent feed streams and side draw streams was made out of convenience of use in the above derivation. Also, explicitly localizing a distinct topology to each tray is arbitrary in light of the final model. Any number of feed streams and/or side streams may be included on any given position as the actual modeling of optimal positioning is composition and flow rate dependent, not position dependent. It is also possible to simply drop the explicit sign of the flow rate imposed upon the same feed streams and side streams and to cancel one of them out from the final equations.

On the other hand, removing the explicit limits on the topology used to derive the generalized column breaks the equations that explicitly yield respective overflows and underflows. The dependence on the explicit physical topology of the derivation restricts their use to a column with said topology or one with fewer feed streams and/or side streams. One could manufacture a framework agnostic of the explicit placement of streams out of the equations derived, but the burden of overcoming the equimolar overflow criteria using quasi-streams and variable elimination would become excessive for a column of more than a few plates. An alternative would be to refactor the model in terms of the continuum basis afforded by the VLE XY-plot, which could potentially produce an extra degree of freedom to solve for when mapping to a physical system.

The final part will be done soon. It will be a concise synopsis of some analytical and pedagogical conclusions from my work in deriving this model.

Print Friendly
Be Sociable, Share!
Short URL: http://dlvr.it/1GBYfS

Sorry, the comment form is closed at this time.