neqsim-distillation-design

name	neqsim-distillation-design
description	Distillation column design rules for NeqSim. USE WHEN: setting up distillation columns, troubleshooting convergence, selecting internals (trays/packing), sizing columns, or analyzing column performance. Covers DistillationColumn setup, solver selection, feed tray optimization, reflux ratio, and internals selection per industry standards.
last_verified	2026-07-04

Distillation Design Rules

Guide for distillation column modeling and design in NeqSim.

Column Setup Pattern

import neqsim.process.equipment.distillation.DistillationColumn;
import neqsim.process.equipment.stream.Stream;
import neqsim.thermo.system.SystemSrkEos;
import neqsim.thermo.system.SystemInterface;

// Create feed
SystemInterface feed = new SystemSrkEos(273.15 + 50.0, 15.0);
feed.addComponent("methane", 0.10);
feed.addComponent("ethane", 0.25);
feed.addComponent("propane", 0.30);
feed.addComponent("n-butane", 0.20);
feed.addComponent("n-pentane", 0.15);
feed.setMixingRule("classic");

Stream feedStream = new Stream("feed", feed);
feedStream.setFlowRate(10000.0, "kg/hr");
feedStream.setTemperature(50.0, "C");
feedStream.setPressure(15.0, "bara");
feedStream.run();

// Create column (name, stages, hasCondenser, hasReboiler)
DistillationColumn column = new DistillationColumn("Deethanizer", 15, true, true);
column.addFeedStream(feedStream, 8);  // Feed on stage 8

// Specifications
column.setCondenserTemperature(273.15 - 30.0);  // Kelvin
column.getReboiler().setHeatInput(1e6);          // Watts
// OR set reflux ratio:
// column.getCondenser().setRefluxRatio(2.5);

column.run();

Solver Selection

// Direct substitution (default) - robust sequential tray sweeps
column.setSolverType(DistillationColumn.SolverType.DIRECT_SUBSTITUTION);

// Inside-Out — faster for ideal/near-ideal systems
column.setSolverType(DistillationColumn.SolverType.INSIDE_OUT);

// Adaptive matrix inside-out — bypasses matrix setup on small columns,
// tries a component-balance matrix warm start on larger columns, then
// finishes with rigorous inside-out polishing
column.setSolverType(DistillationColumn.SolverType.MATRIX_INSIDE_OUT);

// Damped substitution - for difficult convergence
column.setSolverType(DistillationColumn.SolverType.DAMPED_SUBSTITUTION);

// MESH residual monitor - for residual auditing
column.setSolverType(DistillationColumn.SolverType.MESH_RESIDUAL);

// Naphtali-Sandholm - guarded simultaneous MESH residual Newton solver
column.setSolverType(DistillationColumn.SolverType.NAPHTALI_SANDHOLM);
// Optional: seed stage temperatures as initial guesses only, not fixed specs
column.setSeedTemperature(3, 273.15 + 45.0);
// Adjust max iterations
column.setMaxNumberOfIterations(200);

Solver Selection Guide

System Type	Recommended Solver	Notes
Ideal HC (demethanizer, deethanizer)	INSIDE_OUT	Fast, robust
Larger HC fractionators	MATRIX_INSIDE_OUT	Adaptive: small columns bypass matrix overhead; larger columns try the matrix warm start before rigorous inside-out polishing
Non-ideal (alcohols, water)	DAMPED_SUBSTITUTION or DIRECT_SUBSTITUTION	More conservative for non-ideal K-values
Absorbers (no condenser/reboiler)	SUM_RATES or DIRECT_SUBSTITUTION	Flow-corrected updates can help absorber/stripper cases
Wide-boiling (C1 to C20+)	DAMPED_SUBSTITUTION	Increase iterations and monitor residuals
Cryogenic (< -100°C)	INSIDE_OUT with residual diagnostics	Careful with phase identification
Solver audit / residual convergence	MESH_RESIDUAL or NAPHTALI_SANDHOLM	MESH_RESIDUAL records diagnostics and enforces the MESH/product-draw gate by default; NAPHTALI_SANDHOLM attempts guarded simultaneous MESH correction

Column Specification Combinations

Bottom Spec	Top Spec	Comment
Reboiler duty	Condenser temperature	Most common
Reboiler duty	Reflux ratio	Alternative
Top product purity	Reboiler duty	Product-quality control
Bottom temperature	Reflux ratio	Direct T control

// Common specification patterns

// Pattern 1: Condenser T + Reboiler duty
column.setCondenserTemperature(273.15 - 30.0);
column.getReboiler().setHeatInput(1.5e6);

// Pattern 2: Reflux ratio + Reboiler duty
column.getCondenser().setRefluxRatio(3.0);
column.getReboiler().setHeatInput(2.0e6);

// Pattern 3: Product quality plus boilup ratio
column.setTopProductPurity("ethane", 0.98);
column.setReboilerBoilupRatio(2.0);

Reading Column Results

column.run();

// Condenser and reboiler duties
double condenserDuty = column.getCondenser().getDuty();  // Watts
double reboilerDuty = column.getReboiler().getDuty();    // Watts

// Product streams
Stream overhead = (Stream) column.getGasOutStream();
Stream bottoms = (Stream) column.getLiquidOutStream();

// Stage temperatures and compositions
for (int stage = 0; stage < column.getTrays().size(); stage++) {
    double stageTemp = column.getTray(stage).getTemperature() - 273.15;
    // Composition on each stage
}

// Convergence metrics
int iterations = column.getLastIterationCount();
double massResidual = column.getLastMassResidual();
double energyResidual = column.getLastEnergyResidual();
boolean matrixWarmStartUsed = column.wasMatrixInsideOutWarmStartUsed();
boolean matrixWarmStartBypassed = column.wasMatrixInsideOutWarmStartBypassed();
int matrixIterations = column.getLastMatrixInsideOutIterationCount();

Feed Tray Location Rules

Kirkbride Correlation

For binary or pseudo-binary separations:

$$ \log\left(\frac{N_R}{N_S}\right) = 0.206 \log\left[\left(\frac{B}{D}\right) \left(\frac{x_{HK,F}}{x_{LK,F}}\right)^2 \left(\frac{x_{LK,B}}{x_{HK,D}}\right)^2 \right] $$

Where $N_R$ = rectifying stages, $N_S$ = stripping stages, $B/D$ = bottoms/distillate ratio.

Rules of Thumb

Column Type	Feed Tray (from top)	Notes
Demethanizer	40-60% of stages	Light key is very volatile
Deethanizer	50-70% of stages	Moderate volatility
Depropanizer	40-60% of stages	Balanced separation
Debutanizer	50-60% of stages	Similar to depropanizer
Crude column	60-80% of stages	Flash zone near bottom

Minimum Stages and Reflux

Fenske Equation (Minimum Stages)

$$ N_{min} = \frac{\log\left(\frac{x_{LK,D}}{x_{HK,D}} \cdot \frac{x_{HK,B}}{x_{LK,B}}\right)}{\log(\alpha_{LK/HK})} $$

Underwood Equation (Minimum Reflux)

$$ R_{min} = \frac{1}{\alpha - 1}\left(\frac{x_D}{\alpha - \theta} - \frac{1 - x_D}{1 - \theta}\right) $$

Design Heuristics

• Actual stages ≈ 2 × minimum stages (Gilliland correlation)
• Actual reflux ≈ 1.2-1.5 × minimum reflux
• Stage efficiency: 50-70% for trays, 70-90% HETP/stage for packing

Convergence Troubleshooting

Problem	Solution
Column does not converge	Increase max iterations to 200-500
Oscillating temperature profile	Reduce condenser/reboiler specs, use DAMPED_SUBSTITUTION; audit with MESH_RESIDUAL before trying NAPHTALI_SANDHOLM
Wrong product split	Check feed tray location and specifications
Negative flows on stages	Too many stages or wrong specifications
Condenser too cold	Check if subcooled liquid is physical (binary dewpoint)
Reboiler too hot	May be decomposing — check component stability

Steps to Debug

1. Start with fewer stages (5-8) and get convergence
2. Gradually increase stages
3. Use liberal specifications first (higher reflux), then tighten
4. Check feed condition (vapor fraction) — subcooled feed may need enthalpy adjustment
5. Verify component K-values make physical sense at column conditions

Column Sizing (Diameter)

Souders-Brown Correlation

$$ V_{flood} = K_{SB} \sqrt{\frac{\rho_L - \rho_V}{\rho_V}} $$

Where $K_{SB}$ = 0.03-0.07 m/s for trays, 0.02-0.05 for packing.

Design velocity = 70-85% of flooding.

// After running column, get phase properties for sizing
SystemInterface topFluid = column.getTray(0).getFluid();
topFluid.initProperties();
double rhoV = topFluid.getPhase("gas").getDensity("kg/m3");
double rhoL = topFluid.getPhase("oil").getDensity("kg/m3");

double Ksb = 0.05;  // m/s for sieve trays
double Vflood = Ksb * Math.sqrt((rhoL - rhoV) / rhoV);
double Vdesign = 0.80 * Vflood;

double gasFlow = column.getGasOutStream().getFlowRate("m3/hr") / 3600.0;
double area = gasFlow / Vdesign;
double diameter = Math.sqrt(4.0 * area / Math.PI);

Internals Selection

Internals	When to Use	Typical HETP (m)
Sieve trays	General service, fouling	0.5-0.7
Valve trays	Variable turndown	0.4-0.6
Bubble cap trays	Low liquid rates	0.5-0.8
Random packing (Pall rings)	Low pressure drop, corrosive	0.3-0.6
Structured packing (Mellapak)	Vacuum, low ΔP	0.2-0.5

Common Pitfalls

1. Feed flash: Ensure feed is at correct T/P for column conditions
2. Missing components: All components in feed must be present in EOS
3. Mixing rule: Always set before column construction
4. Heavy key in top / light key in bottom: Small amounts are normal — zero means perfect separation (unrealistic)
5. Column pressure profile: Default is constant — set stage pressures for realistic profile
6. Condenser type: Total vs partial condenser changes mass balance