McCabeThiele.jl

Installing and Loading McCabeThiele

McCabeThiele can be installed and loaded either from the JuliaHub repository (last released version) or from the maintainer's repository.

Last Released Version

The last version of McCabeThiele can be installed from JuliaHub repository:

using Pkg
Pkg.add("McCabeThiele")
using McCabeThiele

If McCabeThiele is already installed, it can be updated:

using Pkg
Pkg.update("McCabeThiele")
using McCabeThiele

Pre-Release (Under Construction) Version

The pre-release (under construction) version of McCabeThiele can be installed from the maintainer's repository.

using Pkg
Pkg.add(path="https://github.com/aumpierre-unb/McCabeThiele.jl")
using McCabeThiele

Citation of McCabeThiele

You can cite all versions (both released and pre-released), by using 10.5281/zenodo.7126164.

This DOI represents all versions, and will always resolve to the latest one.

For citation of the last released version of InternalFluidFlow, please check CITATION file at the maintainer's repository.

The McCabeThiele Module for Julia

McCabeThiele provides the following functions:

stages
refmin
qR2S

stages computes the number of theoretical stages of a distillation column using the method of McCabe-Thiele, given a function y = y(x) that relates the liquid fraction x and the vapor fraction y, or a x-y matrix of the liquid and the vapor fractions, the vector of the fractions of the products and the feed, the feed quality, and the reflux ratio at the top of the column.

If feed is a saturated liquid, feed quality q = 1, feed quality is reset to q = 1 - 1e-10.

By default, theoretical stages are computed from the stripping section to the rectifying section, updown = true.

If updown = false is given, theoretical stages are computed from the rectifying section to the stripping section.

By default, stages plots a schematic diagram of the solution, fig = true.

If fig = false is given, no plot is shown.

Syntax:

N=stages(y,X,q,R[,updown[,fig]])

Examples:

Compute the number of theoretical stages of a distillation column from the bottom of the column, given a matrix that relates the liquid fraction and the vapor fraction, the composition of the column's bottom product is 11 %, the composition of the distillate is 88 %, the composition of the feed is 46 %, the feed quality is 54 %, and the reflux ratio at the top of the column 70 % higher that the minimum reflux ratio:

data=[0.  0.;
      0.1 0.212;
      0.2 0.384;
      0.3 0.529;
      0.4 0.651;
      0.5 0.752;
      0.6 0.833;
      0.7 0.895;
      0.8 0.942;
      0.9 0.974;
      1.  1.];
x=[0.88 0.46 0.11];
q=0.54;
r=refmin(data,x,q)
R=1.70*r;
N=stages(data,x,q,R,false,false)

Compute the number of theoretical stages of a distillation column from the top of the column, given the function that compute the vapor fraction given the liquid fraction, the composition of the column's bottom product is 11 %, the composition of the distillate is 88 %, the composition of the feed is 46 %, the feed quality is 54 %, and the reflux ratio at the top of the column 70 % higher that the minimum reflux ratio, and plot a schematic diagram of the solution:

y(x)=x.^1.11 .* (1-x).^1.09 + x;
x=[0.88 0.46 0.11];
q=0.54;
r=refmin(y,x,q)
R=1.70*r;
N=stages(y,x,q,R)

refmin

refmin computes the minimum value of the reflux ratio of a distillation column, given a function y = y(x) that relates the liquid fraction x and the vapor fraction y, or a x-y matrix of the liquid and the vapor fractions, the vector of the fractions of the distillate and the feed, and the feed quality.

If feed is a saturated liquid, feed quality q = 1, feed quality is reset to q = 1 - 1e-10.

Syntax:

r=refmin(y,X,q)

Examples:

Compute the minimum value of the reflux ratio of a distillation column, given a matrix that relates the liquid fraction and the vapor fraction, the composition of the column's bottom is 11 %, the composition of the distillate is 88 %, and the feed quality is 54 %:

data=[0.  0.;
      0.1 0.212;
      0.2 0.384;
      0.3 0.529;
      0.4 0.651;
      0.5 0.752;
      0.6 0.833;
      0.7 0.895;
      0.8 0.942;
      0.9 0.974;
      1.  1.];
x=[0.88 0.46];
q=0.54;
r=refmin(data,x,q)

y(x)=x.^1.11 .* (1-x).^1.09 + x;
x=[0.88 0.46];
q=0.54;
r=refmin(y,x,q)

qR2S

qR2S computes the reflux ratio at the bottom of the column, given the reflux ratio at the top of the column, the vector of the fractions of the products and the feed, and the feed quality.

If feed is a saturated liquid, feed quality q = 1, feed quality is reset to q = 1 - 1e-10.

Syntax:

S=qR2S(R,X,q)

Examples:

Compute the reflux ratio at the bottom of the column, given the reflux ratio R = 2 at the top of the column, the composition xB = 11 % of the column's bottom, the composition xD = 88 % of the distillate, the composition xF = 46 % of the feed, and the feed quality q = 54 %:

R=2;
x=[0.88 0.46 0.11];
q=0.54;
S=qR2S(R,x,q)

email: aumpierre@gmail.com