Full-Process Design Review of a 300,000 t/y Methanol Distillation Column

W4 Column Design Series: Final Engineering Case Study

This case study connects the principal elements of column design covered in the previous six articles: valve-tray fundamentals, design standards, hydraulic calculations, high-efficiency trays, operating optimization, and simulation tools.

The project involves the main methanol distillation column T-202 at a coal chemical complex with a stated capacity of 300,000 tonnes per year. The review covers the complete design workflow, from process-package data and preliminary hand calculations to rigorous RadFrac modeling, KG-TOWER hydraulic verification, and final construction drawings.

Process Data Analysis

The coal chemical complex operates a methanol synthesis unit with a stated capacity of 300,000 t/y and 8,000 operating hours per year. Crude methanol from the synthesis section first enters the pre-column T-201, where light components such as dimethyl ether and dissolved gases are removed.

The bottoms stream from T-201 then enters the main column T-202 for methanol–water separation.

T-202 Design Conditions

• Feed rate: 38.5 t/h
• Feed composition: 82.0 wt% methanol and 18.0 wt% water
• Feed condition: liquid at its bubble point
• Overhead product: 36.5 t/h at 99.95 wt% methanol
• Bottoms wastewater: 2.0 t/h with less than 100 ppm methanol
• Operating pressure: 0.03 MPaG at the column top
• Property method: NRTL, with vapor–liquid equilibrium data taken from the Aspen Plus APV72 database

The feed mole fractions are calculated using molecular weights of 32.04 for methanol and 18.02 for water.

A feed rate of 38.5 t/h contains:

• Methanol: 31.57 t/h, equivalent to 985.9 kmol/h
• Water: 6.93 t/h, equivalent to 384.6 kmol/h
• Total molar flow: 1,370.5 kmol/h

The corresponding feed composition is:

xF = 0.719 mole fraction methanol

The stated product compositions are:

• Overhead composition: xD = 0.997
• Bottoms composition: xB = 0.00006

Preliminary Hand Calculations

Minimum Number of Theoretical Stages by the Fenske Equation

The average relative volatility of the methanol–water system is assumed to be:

αavg ≈ 4.5

This value represents an integrated average from the NRTL model across the full concentration range. The relative volatility is approximately 7.5 in the low-methanol region and approximately 3.5 in the high-methanol region.

The minimum number of theoretical stages is calculated as:$$N_{\min}= \frac{ \log\left[ \left(\frac{x_D}{1-x_D}\right) \left(\frac{1-x_B}{x_B}\right) \right] }{ \log(\alpha) }$$

Substituting the stated values:$$N_{\min}= \frac{ \log\left[ \left(\frac{0.997}{0.003}\right) \left(\frac{0.99994}{0.00006}\right) \right] }{ \log(4.5) }$$$$N_{\min}= \frac{\log(332.3 \times 16{,}657)} {\log(4.5)}$$$$N_{\min}= \frac{\log(5{,}538{,}000)}{1.504} = \frac{6.743}{1.504} \approx 4.5$$

The preliminary result is rounded to:

Nmin = 5 stages, including the reboiler

Minimum Reflux Ratio

For a bubble-point liquid feed:$$q=1$$

The equilibrium vapor composition corresponding to the feed composition is estimated as:$$y_F^*= \frac{\alpha x_F} {1+(\alpha-1)x_F}$$$$y_F^*= \frac{4.5\times0.719} {1+3.5\times0.719} = \frac{3.236}{3.517} = 0.920$$

The stated minimum reflux ratio is:$$R_{\min}= \frac{x_D-y_F^*} {y_F^*-x_F}$$$$R_{\min}= \frac{0.997-0.920} {0.920-0.719} = \frac{0.077}{0.201} = 0.38$$

Therefore:

Rmin = 0.38

Actual Number of Theoretical Stages Using the Gilliland Correlation

The selected operating reflux ratio is:$$R=1.8$$

This gives:$$\frac{R}{R_{\min}}=4.7$$

Although relatively high, this reflux ratio is intended to provide AA-grade methanol purity and sufficient operating flexibility.

The Gilliland variables are calculated as:$$X= \frac{R-R_{\min}}{R+1} = \frac{1.8-0.38}{2.8} = 0.507$$$$Y= 0.75(1-X^{0.5668})$$$$Y= 0.75(1-0.507^{0.5668}) = 0.75(1-0.658) = 0.256$$

The theoretical stage count is then estimated as:$$N= \frac{N_{\min}+Y}{1-Y}$$$$N= \frac{5+0.256}{0.744} = 7.1$$

This value includes the condenser and reboiler.

However, this result is based on a constant relative-volatility assumption. In the high-purity methanol region, relative volatility decreases to below 3.5. Using a constant value of 4.5 can therefore underestimate the required number of theoretical stages.

The rigorous RadFrac vapor–liquid equilibrium model indicated a requirement of:

24 theoretical stages, including the condenser and reboiler

This value is approximately 3.4 times the preliminary estimate. The result illustrates an important design principle:

Hand calculations establish the preliminary direction, while rigorous simulation provides the required design accuracy.

Preliminary Column Diameter Using the Souders–Brown Method

The stated overhead conditions are:

• Temperature: 65°C
• Absolute pressure: 0.13 MPaA
• Vapor density: 1.48 kg/m³
• Liquid density: 750 kg/m³

The overhead vapor flow is estimated as:$$V=D(R+1)$$$$V=985.9\times2.8 = 2{,}760.5\ \text{kmol/h}$$

The corresponding stated values are:

• Vapor mass flow: 88,440 kg/h
• Vapor volumetric flow: 19.2 m³/s

The vapor–liquid flow parameter is:$$F_{LV}= \frac{L}{V} \sqrt{\frac{\rho_V}{\rho_L}}$$$$F_{LV}= \frac{1.8}{2.8} \sqrt{\frac{1.48}{750}}$$$$F_{LV}= 0.643\times0.0444 = 0.0286$$

For a tray spacing of 600 mm, the Fair correlation gives:$$C_{SB0}=0.082\ \text{m/s}$$

The flooding velocity is:$$u_{\text{flood}}= 0.082 \sqrt{\frac{750-1.48}{1.48}}$$$$u_{\text{flood}}= 0.082\times22.5 = 1.84\ \text{m/s}$$

Using 80% of the flooding velocity:$$u=1.47\ \text{m/s}$$

The required cross-sectional area is:$$A= \frac{19.2}{1.47} = 13.1\ \text{m}^2$$

The corresponding diameter is:$$D= \sqrt{\frac{4\times13.1}{\pi}} = 4.08\ \text{m}$$

The preliminary nominal diameter is therefore selected as:

DN4000

Figure 1. T-202 column arrangement with 48 F1 valve trays, a nominal diameter of DN4000, and 24 trays each in the rectifying and stripping sections

Rigorous RadFrac Modeling

The hand calculations establish the preliminary design range, while the rigorous model determines the final separation requirements.

The column contains 48 actual trays. Assuming an overall tray efficiency of 50%, these trays are represented by 24 theoretical stages. The total condenser is defined as stage 1, and the kettle reboiler is defined as stage 24.

The feed enters at stage 12.

RadFrac Model Settings

• Property method: NRTL-RK
• Liquid-phase model: NRTL activity-coefficient model
• Vapor-phase model: Redlich–Kwong equation of state
• Number of stages: 24, including the total condenser and kettle reboiler
• Feed stage: Stage 12
• Feed condition: Bubble-point liquid calculated automatically by Aspen Flash
• Operating specifications: Reflux ratio of 1.8 and distillate flow of 985.9 kmol/h
• Convergence method: Newton–Raphson
• Convergence performance: Four iterations with a residual below 10⁻⁵
• Design specification: Overhead methanol concentration of at least 99.95 wt%, with the reflux ratio adjusted automatically

Converged Results

• Overhead product: 36.5 t/h at 99.952 wt% methanol
• Bottoms wastewater: 2.0 t/h containing 85 ppm methanol
• Condenser duty: 12.8 MW
• Cooling-water temperature rise: 15°C
• Stated cooling-water flow: 61.3 t/h
• Reboiler duty: 13.5 MW
• Heating medium: 0.3 MPaG steam
• Stated steam consumption: 21.8 t/h
• Column-top temperature: 64.8°C
• Feed-stage temperature: 78.3°C
• Column-bottom temperature: 99.2°C

Detailed Hydraulic Verification in KG-TOWER

The vapor and liquid flow profiles obtained from RadFrac are imported into KG-TOWER for tray-by-tray hydraulic verification.

The specified tray parameters are:

• Column diameter: DN4000
• Tray spacing: 600 mm
• Weir height: 50 mm
• Weir length: 2.4 m, equivalent to 0.6D
• Number of F1 valves: 850
• Open-area ratio: 12.8%

Flooding Check

The top tray is identified as the controlling hydraulic condition.$$F_{LV}=0.0286$$$$C_{SB}=0.082\ \text{m/s}$$$$u_{\text{flood}}=1.84\ \text{m/s}$$

The actual superficial vapor velocity is stated as:$$u= \frac{19.2}{12.57} = 1.53\ \text{m/s}$$

The flooding percentage is therefore:$$\text{Flooding percentage}= \frac{1.53}{1.84} = 83.2\%$$

Result: 83.2%, within the stated preferred range of 60%–85%

Valve-Hole Kinetic Energy Factor

For a single-valve opening diameter of 49 mm, the total hole area is:$$A_h= 850\times\frac{\pi}{4}\times0.049^2 = 1.61\ \text{m}^2$$

The hole velocity is:$$u_{\text{hole}}= \frac{19.2}{1.61} = 11.9\ \text{m/s}$$

The kinetic energy factor is:$$F_0= u_{\text{hole}}\sqrt{\rho_V}$$$$F_0= 11.9\sqrt{1.48} = 14.5$$

Result: F₀ = 14.5, within the stated preferred range of 8–15

This value indicates adequate resistance to weeping under the design condition.

Downcomer Residence Time

The stated downcomer area on one side is:$$A_d=0.754\ \text{m}^2$$

The downcomer volume is:$$V_d= 0.754\times0.6 = 0.452\ \text{m}^3$$

The liquid flow rate is:$$87.6\ \text{m}^3/\text{h} = 0.0243\ \text{m}^3/\text{s}$$

The residence time is:$$\tau= \frac{0.452}{0.0243} = 18.6\ \text{s}$$

Result: 18.6 seconds, exceeding the stated minimum requirement of five seconds

The downcomer therefore provides sufficient time for vapor disengagement under the stated design conditions.

Entrainment

Using the Glitsch correlation at a flooding percentage of 83.2%:$$e_v\approx0.038\ \text{kg liquid/kg vapor}$$

Result: 0.038 kg liquid/kg vapor, below the stated limit of 0.1

Operating Flexibility

At 50% load:$$F_0=7.3$$

This remains above the stated minimum value of 5, indicating no significant weeping risk.

At 120% load, the predicted flooding percentage is:$$99.6\%$$

This condition is close to the flooding limit.

The stated turndown ratio is:$$\text{Operating range}= \frac{120\%}{50\%} = 2.4$$

Result: An operating range of 2.4, subject to the limited hydraulic margin at maximum load

Five Critical Design Considerations

Select an Appropriate Thermodynamic Model

The methanol–water system is strongly nonideal and polar. A liquid-phase activity-coefficient model such as NRTL is generally more appropriate than relying solely on cubic equations of state such as Peng–Robinson or SRK.

The original case attributes vapor–liquid equilibrium deviations of up to 30% in the high-methanol region to inappropriate property-model selection. Such deviations can cause substantial errors in the predicted stage requirement.

The NRTL-RK combination uses NRTL for the liquid phase and the Redlich–Kwong equation of state for the vapor phase.

High Purity Should Not Depend Only on Additional Trays

For AA-grade methanol production, purity improvement may involve more than increasing the total tray count.

The proposed arrangement withdraws product from trays 2–4 while operating the top stage as a purification or reflux section. Under suitable design conditions, this configuration may provide higher product purity than direct withdrawal from the extreme top of the column.

However, the terminology and configuration must be defined carefully because a conventional total condenser does not contain an equilibrium tray above the product draw in the same sense as an internal tray section.

Control Entrainment When Bottoms Methanol Must Remain Below 100 ppm

When the methanol concentration in the bottoms must remain below 100 ppm, entrainment in the lower section requires careful evaluation.

The stated concern is that vapor from the lower trays may carry methanol-containing liquid upward. Proposed measures include increasing the spacing of the bottom three to five trays to 800 mm or using directional valve trays to reduce liquid-level gradients.

The actual relationship between entrainment and bottoms methanol concentration should be confirmed using a complete component balance and tray-by-tray simulation.

Verify Feed Flashing

If the T-201 bottoms stream enters T-202 above its bubble-point temperature, flashing may occur on the feed tray and disturb the internal vapor–liquid traffic.

A feed cooler can be used to maintain the feed temperature within approximately ±2°C of its bubble point. In Aspen Plus, specifying a feed vapor fraction of zero allows the corresponding bubble-point condition to be calculated.

Prevent Cold-Column Operation in Winter

A nominally atmospheric column may operate under slight vacuum during winter if heat losses are excessive. A low top temperature can subcool the reflux and reduce the internal vapor flow, increasing the risk of tray weeping.

Potential measures include improving column insulation and installing a reflux heater to maintain the required reflux temperature. The source case proposes a reflux temperature of at least 60°C.

الخاتمة

The design of T-202 follows five principal stages:

Process data analysis → preliminary hand calculations → rigorous simulation → hydraulic verification → construction drawing delivery

Three main engineering lessons are presented.

Hand Calculations Establish the Preliminary Range

Fenske and Underwood calculations can provide an initial design range. However, assuming constant relative volatility may substantially underestimate the number of stages required in the high-purity methanol region.

RadFrac Provides Separation Accuracy

Under the stated NRTL-based vapor–liquid equilibrium model, the column requires 24 theoretical stages rather than the seven stages predicted by the preliminary calculation.

The difference demonstrates the importance of rigorous thermodynamic modeling for high-purity separation.

KG-TOWER Evaluates Hydraulic Feasibility

The stated design produces:

• Flooding percentage: 83.2%
• Valve-hole kinetic energy factor: 14.5
• Downcomer residence time: 18.6 seconds
• Entrainment: 0.038 kg liquid/kg vapor
• Operating range: 2.4

These results suggest that the design is close to the upper end of the preferred hydraulic operating range. Additional verification is required before the column can be considered ready for detailed engineering or construction.

This case concludes the seven-part W4 column design series:

Valve-tray fundamentals → design standards → hydraulic calculations → high-efficiency trays → operating optimization → simulation tools → engineering case study

The next series will focus on reactor design, beginning with the fundamentals of stirred-tank reactors.

Discussion

For comparable methanol distillation columns, key operating data include the selected valve-tray type, measured tray efficiency, flooding margin, turndown performance, and actual separation efficiency. Comparing these values across operating plants can provide useful guidance for future designs.