Some books and learning resources use steady-state distillation columns as an example problem to introduce linear systems of equations in linear algebra and numerical computing; for instance:
- https://iastate.pressbooks.pub/chemicalengineeringseparations/chapter/distillation/
- Chapter 2.1 in Mostoufi, Constantinides - Applied Numerical Methods for Chemical Engineers (2022)
- https://www.chegg.com/homework-help/questions-and-answers/question-2-20-marks-two-stage-distillation-column-shown-figure-2-used-separate-mixture-met-q95394507#question-transcript
- https://cse.sc.edu/~gatzke/550/linalg (probably a simpler example, since it does not seem to involve intermediate stages with condensers but just a balance of mass; this one is clear to me).
I am an applied mathematician, and I have difficulties understanding the modelling part of these examples. I would like to get an explanation of the general idea of the chemical process behind these examples, assuming little previous knowledge of chemical engineering. From what I understand, these columns contain a sequence of "stages" that separate two different chemicals with different boiling points; but it is unclear to me how the 'balance' associated to these stages is computed. For instance, from the first reference:
$L$ mass or molar flow rate of the liquid reflux returned to the column from the condenser (mass time$^-1$ or mole time$^-1$); also generic flow rate of the liquid phase in the rectifying section
$\overline L$ = mass or molar flow rate of the liquid leaving the bottom of the column and entering the reboiler (mass time-1 or mole time-1); also generic flow rate of the liquid phase in the stripping section
$n$ = generic stage number, stage 1 is at the top of the column
$V$ = mass or molar flow rate of vapor leaving the top of the column and entering the condenser (mass time-1 or mole time-1); also generic flow rate of the vapor phase in the rectifying section
$\overline V$ = mass or molar flow rate of the gaseous boilup returned to the column from the reboiler (mass time-1 or mole time-1); also generic flow rate of the vapor phase in the stripping section
$x_n$ = mass or mole fraction of the light key in the liquid leaving stage $n$
$y_n$ = mass or mole fraction of the light key in the vapor leaving stage $n$
[...]
Material balance on stages $1$-$n$, the rectifying section of the column
$$ y_{n+1} = \frac{L}{V} x_n + y_1 - \frac{L}{V}x_0 $$
What are L and V? How does one obtain the equation? Which quantities are known and which are unknown?
You can assume a simplified model for clarity: e.g., two chemical species, ignoring the units at the top/bottom named "reboiler/condenser" and focus on the column itself.
