How To Model Phase Equilibrium For Stripping Weak Electrolytes? Guide For Wastewater Pilot Plants

Modeling the phase equilibrium for stripping weak volatile electrolytes is not a standard vapor-liquid equilibrium (VLE) problem. It demands a coupled, multi‑constraint thermodynamic framework that simultaneously resolves chemical dissociation, ionic activity, physical solubility, mass balances, and electroneutrality. You cannot simply apply Raoult’s or Henry’s Law in isolation—the system’s true driving force depends on the concentration of the neutral molecular species, which is controlled by pH‑dependent chemical reactions.

Stripping ammonia, hydrogen sulfide, or sulfur dioxide from wastewater requires modeling the liquid phase as a reactive electrolyte solution. The phase equilibrium is built by combining chemical dissociation constants (for liquid‑phase ionic equilibrium) with Henry’s constants (for molecular solute vapor-liquid equilibrium), all expressed through unsymmetrically normalized activity coefficients that approach unity at pure water, and constrained by elemental mass balances and electroneutrality.

The Dual‑Equilibrium Challenge

Weak volatile electrolytes like ammonia ($NH_3$), hydrogen sulfide ($H_2S$), and sulfur dioxide ($SO_2$) partially dissociate in water. The resulting ionic species ($NH_4^+$, $HS^-$, $S^{2-}$, $HSO_3^-$, $SO_3^{2-}$) are essentially nonvolatile. Only the neutral molecular form can transfer to the vapor phase.

This creates a fundamental modeling complexity. You must first compute how much of the total analyte exists as the free molecule, then apply the appropriate vapor‑phase relationship to that fraction.

Why a Simple Henry’s Law is Not Enough

A direct Henry’s Law – $p_i = H_i x_i$ – fails because the composition variable $x_i$ typically represents total analytical concentration. That total includes ionized forms that do not vaporize.

Using the apparent mole fraction of “ammonia” without correcting for dissociation will grossly over‑predict stripping rates. The model must isolate the molecular molality, $m_{i,mol}$. That quantity is hidden behind a dissociation equilibrium.

The Unusual Reference State in Dilute Electrolyte Solutions

In non‑electrolyte mixtures, activity coefficients are often symmetric, reducing to unity for pure component $i$. But here the pure volatile electrolyte liquid is not a realistic reference for ions in water.

The proper approach is an unsymmetric normalization. Activity coefficients $\gamma_i^$ are defined so that $\gamma_i^ \rightarrow 1$ as the solution becomes pure solvent (water). This convention decouples the model from the unphysical pure‑solute state and correctly describes the highly dilute ionic environment.

Building the Thermodynamic Framework

A rigorous model rests on four simultaneous conditions. Each must be enforced at every stage of a column or pilot‑plant simulation.

1. Mass Balances for Every Reactive Element

Total concentrations of elements (nitrogen, sulfur, sodium, etc.) are conserved across molecular and ionic forms.

For a single weak electrolyte like ammonia: $$ m_{NH_3,total} = m_{NH_3,mol} + m_{NH_4^+} $$ For systems with multiple acids and bases, you write balances linking all species containing sulfur, carbon, or nitrogen. These balances are the backbone of the computation.

2. Chemical Dissociation Equilibrium in the Liquid

Dissociation is treated as a reversible reaction with a thermodynamic equilibrium constant $K$. For example: $$ K_{NH_3} = \frac{a_{NH_4^+} , a_{OH^-}}{a_{NH_3}} $$ Here, $a_i$ denotes activity, typically expressed as the product of molality and the unsymmetric activity coefficient: $a_i = \gamma_i^* , m_i$.

Because $K$ depends strongly on temperature, pilot plants must record temperature profiles to match the correct dissociation constant. pH control directly manipulates these equilibria, shifting the molecular fraction and therefore the volatility.

3. Electroneutrality Constraint

The liquid phase carries no net charge. The sum of positive ionic charges must equal the sum of negative charges: $$ \sum z_+ m_+ = \sum z_- m_- $$ This condition locks the model. Any change in pH or added salt shifts the distribution of all ions, altering the molecular fraction available for stripping. It closes the system of equations and makes the problem deterministic.

4. Vapor‑Liquid Equilibrium for the Molecular Species

Once the molecular molality $m_{i,mol}$ is known from the equilibria above, physical phase equilibrium is described by an extended Henry’s Law: $$ p_i = y_i P = H_i , \gamma_i^* , m_{i,mol} , / , \phi_i^V $$

$H_i$ is the Henry’s constant (temperature‑dependent).
$\gamma_i^*$ is the unsymmetrically normalized activity coefficient of the molecular species.
$\phi_i^V$ is the vapor‑phase fugacity coefficient, often near unity at low pressure but included for completeness.

Pilot‑plant data for off‑gas composition directly validate this equation.

Capturing Non‑Idealities: Activity Coefficient Models

When ionic strength rises—due to added salts, multiple weak electrolytes, or pH adjusters—the molecular activity coefficient $\gamma_i^*$ deviates from unity. Simply ignoring that deviation introduces error into stripping efficiency predictions.

Pitzer‑Based Models for Ionic Interactions

For concentrated or high‑ionic‑strength solutions, Pitzer’s theory offers a robust semi‑empirical framework. It calculates activity coefficients for both molecular and ionic species from binary and ternary interaction parameters.

The Pitzer model explicitly accounts for short‑range ion‑ion and ion‑molecule interactions. In ammonia‑$CO_2$‑$H_2S$‑water systems common to sour water strippers, this approach has become a standard because it reliably predicts volatility shifts with salt content and pH.

Regression from Pilot‑Plant Data to Improve Accuracy

Activity coefficient models (NRTL, UNIQUAC, or Pitzer) rely on binary interaction parameters. These parameters are not universal constants—they are regressed from experimental data and lose reliability as mixture complexity grows.

Pilot‑plant runs generate the very data needed to refit those parameters. By measuring temperature, pressure, liquid‑phase composition, and vapor‑phase composition across the column, researchers can feed regression subroutines in process simulators. The outcome is a tuned model that captures the specific waste matrix, dramatically improving design reliability.

Understanding the Trade‑offs

No single model is universally optimal. The choice carries inherent tensions between accuracy, experimental burden, and computational effort.

Detail versus Practicality

A full Pitzer treatment with dozens of interactive parameters can describe almost any weak electrolyte mixture. However, each additional parameter requires more experimental data to fit.

In many dilute municipal wastewaters, the ionic strength is low enough that $\gamma_i^* \approx 1$ is acceptable. Here, the simpler model of “$K + H$” with ideal dilutions still gives acceptable stripping predictions. The engineer must decide if the added accuracy justifies the extensive pilot efforts needed for parameterization.

Handling Gas‑Phase Non‑Ideality

At atmospheric pressure, $\phi_i^V \approx 1$ often suffices. But in pressurized stripping columns or when the gas contains high molecular‑weight compounds, fugacity corrections become significant. Neglecting them introduces bias that grows with pressure. Pilot plants that recycle gas or operate under slight pressure must verify this assumption.

The Hidden Trap of Incomplete Speciation

A common modeling pitfall is to account for only the first dissociation step. For hydrogen sulfide, the presence of $HS^-$, $S^{2-}$, and possibly polysulfides must be considered when pH exceeds ~7. Missing the second dissociation leads to an under‑prediction of the sulfidic load remaining in the liquid and an over‑prediction of stripping performance—an error that can cause an environmental violation.

How to Apply This to Your Pilot Plant

The best approach depends on your primary objective. Use these goal‑oriented strategies to select and implement the right phase‑equilibrium model.

If your primary focus is rapid scoping for dilute, single‑contaminant streams: Start with the simple coupled model: use the dissociation constant $K$ to compute molecular fraction, apply Henry’s constant at infinite dilution, and assume $\gamma_i^* = 1$. Validate with a few pH‑variation runs to confirm sensitivity.
If your primary focus is high‑fidelity design for complex, high‑ionic‑strength mixtures (e.g., sour water strippers): Implement a Pitzer or eNRTL framework. Run the pilot plant to generate a matrix of temperature, pH, and ionic strength data, then regress binary interaction parameters. This investment pays off when scaling up to industrial units.
If your primary focus is real‑time process control and automation: Embed the equilibrated thermodynamic model (electroneutrality + $K$ + Henry) into the control logic. Use pH and temperature as feed‑forward variables to predict off‑gas load and adjust stripping steam or air flow immediately.
If your primary focus is comparing chemical dosing strategies for pH adjustment: Model the full acid‑base system, including the added cation/anion. A weak electrolyte model that ignores the counter‑ion from the neutralization chemical will mispredict the pH response and thus the stripping behavior. Include those ionic species in the electroneutrality equation.

Resolving phase equilibrium for these reactive systems transforms pilot‑plant data from simple performance curves into a predictive, transferable model. That model becomes the foundation for safe, cost‑effective, and regulation‑compliant full‑scale design.

Summary Table:

Modeling Component	Operational Role & Purpose	Key Equation / Constraint
Chemical Dissociation	Determines the fraction of volatile molecular species vs. nonvolatile ions	$K = \frac{a_{ion}^+ \cdot a_{ion}^-}{a_{molecular}}$
Physical Phase Equilibrium	Relates volatile molecular species in liquid to vapor phase	$p_i = H_i \gamma_i^* m_{i,mol}$ (Extended Henry's Law)
Mass Balances	Conserves total analytical concentration across all species	$m_{total} = m_{molecular} + \sum m_{ionic}$
Electroneutrality	Ensures liquid phase carries no net charge, locking the pH equilibrium	$\sum z_+ m_+ = \sum z_- m_-$
Activity Models (Pitzer/eNRTL)	Corrects for ionic strength and non-idealities in concentrated solutions	$\gamma_i^*$ (Unsymmetric activity coefficient)

Maximize the Accuracy of Your Environmental & Water Treatment Pilot Plants

Translating complex thermodynamic phase equilibria into reliable, scalable industrial processes requires high-fidelity experimental data.

LABPARK designs and delivers advanced Educational and Vocational Unit Operations Pilot Plants across chemical engineering, bioprocess & biotech, and environmental & water treatment. Specifically tailored for universities, research institutes, and enterprises, our pilot systems enable precise measurement of temperature profiles, pH adjustments, and off-gas compositions—allowing you to regress binary interaction parameters with confidence and ensure regulation-compliant designs.

Ready to elevate your engineering training and research capabilities? Contact us today to explore our pilot plant solutions!