A gas processing company needs to upgrade a natural gas stream from 8% CO2 to pipeline spec (< 2% CO2) at a remote wellhead where amine treating is not practical. Membrane separation is the obvious choice: compact, no chemicals, low maintenance. The engineer selects a commercial polyimide membrane with published CO2/CH4 selectivity of 30 and CO2 permeability of 10 Barrer. The process model predicts a 2-stage system with 500 m2 of membrane area will meet spec at 90% methane recovery.
After installation, the system achieves only 82% methane recovery. The CO2 permeation rate is 30% lower than the published value because the feed gas contains 200 ppm H2S (which competes for sorption sites in the membrane), the operating temperature is 45C (not the 35C at which the published data was measured), and plasticization by CO2 at the high partial pressure in the feed has reduced the effective selectivity from 30 to 18.
These are not installation errors or manufacturing defects. They are predictable consequences of molecular-level phenomena in the membrane material that the solution-diffusion model, with its constant permeability and selectivity coefficients, cannot capture. Every membrane process operates in conditions that differ from the conditions where membrane properties were measured. And the deviations are not small. They routinely cause 20-40% discrepancies between design predictions and field performance.
How Membranes Separate Gases
Dense polymer membranes separate gases through the solution-diffusion mechanism. Gas molecules dissolve into the upstream face of the membrane (sorption), diffuse through the polymer matrix driven by the concentration gradient, and desorb from the downstream face. The permeability P is the product of the solubility coefficient S (how much gas the polymer absorbs) and the diffusion coefficient D (how fast the gas moves through the polymer).
Selectivity between two gases (say CO2 and CH4) depends on both the solubility selectivity (S_CO2/S_CH4) and the diffusivity selectivity (D_CO2/D_CH4). CO2 is more soluble than CH4 in most polymers due to its quadrupole moment and condensability. CO2 is also more mobile than CH4 in glassy polymers because its kinetic diameter is smaller (3.3 A vs 3.8 A). Both factors favor CO2 permeation, giving the useful CO2/CH4 selectivity that makes membrane gas processing possible.
But the standard solution-diffusion model treats S and D as constants, properties of the gas-polymer pair measured at one condition and applied everywhere. In reality, both S and D depend on temperature, pressure, gas composition, and the physical aging history of the membrane. These dependencies are governed by molecular-level interactions between gas molecules and the polymer matrix.
Plasticization: The Performance Killer
Plasticization is the most economically significant deviation from ideal membrane behavior. When CO2 dissolves into a glassy polymer at high partial pressure, it increases the polymer chain mobility, effectively reducing the glass transition temperature. This swells the polymer, increases the free volume, and dramatically increases the permeability of all gas species. The effect on selectivity is devastating: the membrane becomes more permeable but less selective, because the larger free volume allows CH4 (the larger molecule) to permeate faster.
Plasticization onset occurs at a critical CO2 partial pressure that depends on the polymer chemistry, crosslink density, and thermal history. For standard polyimide membranes, onset is typically 8-15 bar CO2 partial pressure. At 8% CO2 and 70 bar total feed pressure, the CO2 partial pressure is 5.6 bar, close to the onset threshold. At 15% CO2 (some wells), the partial pressure exceeds onset and selectivity degrades rapidly.
The standard solution-diffusion model with constant coefficients does not predict plasticization. You either need experimental data at the actual operating conditions (which requires fabricating and testing the membrane at high pressure, a non-trivial experiment) or molecular simulation that computes the CO2-polymer interaction from first principles.
Molecular Simulation of Gas Transport in Polymers
Molecular dynamics simulation of gas permeation starts with building a realistic model of the polymer membrane. An amorphous polymer packing (typically 5,000-50,000 atoms) is generated by placing polymer chains in a simulation box and equilibrating at the target temperature. The resulting structure contains the free volume elements (transient voids between polymer chains) through which gas molecules hop.
Our physics-first simulation approach directly computes gas transport by inserting gas molecules (CO2, CH4, N2, H2S) into the polymer matrix and running molecular dynamics. The simulation captures gas dissolution, with CO2 molecules occupying free volume elements and interacting with polar groups in the polymer (carbonyl, imide, ether). It captures gas diffusion, with molecules hopping between free volume elements as polymer chains undergo segmental motion. And it captures the composition dependence: as more CO2 dissolves, it swells the matrix, creating larger free volume elements and faster diffusion for all species.
The simulation predicts the sorption isotherm (solubility vs. pressure, including the concave-to-convex transition that signals plasticization onset), the diffusion coefficient (as a function of gas concentration in the membrane), and the permeability and selectivity at any operating condition (pressure, temperature, gas composition). These predictions emerge from the molecular interactions without fitted parameters.
Mixed-Gas Effects and Competitive Sorption
Membrane properties measured with pure gases do not predict mixed-gas performance. When CO2 and CH4 are present simultaneously, they compete for sorption sites in the polymer. CO2 preferentially occupies the Langmuir sites (microvoids in the glassy polymer), reducing the available sorption capacity for CH4. The result: CH4 permeability in the mixture is lower than the pure-gas value, while CO2 permeability may be similar or slightly lower. The net effect on selectivity can go either way.
H2S presents an even more complex scenario. H2S is highly condensable and strongly interacting with polar polymer groups. Even at 200 ppm in the feed, H2S can occupy a disproportionate fraction of sorption sites, reducing CO2 sorption and permeation. At higher concentrations (1-5%), H2S can plasticize the membrane more aggressively than CO2, causing catastrophic selectivity loss.
Molecular simulation captures mixed-gas effects naturally. The simulation includes all gas species at their actual partial pressures. Competition for sorption sites, plasticization from multiple condensable gases, and the effect of gas-gas interactions within the polymer matrix are all computed from the molecular forces. The result is a mixed-gas permeability that can differ from the pure-gas value by 20-50%, a deviation that changes the required membrane area and the achievable recovery.
Physical Aging: Why Membrane Performance Drifts
Glassy polymer membranes are not in thermodynamic equilibrium. The polymer chains are frozen in a non-equilibrium packing with excess free volume compared to the equilibrium liquid state. Over time (months to years) the polymer slowly relaxes toward equilibrium, reducing free volume, decreasing permeability, and increasing selectivity. This physical aging causes membrane performance to drift over the operational lifetime.
The aging rate depends on the polymer chemistry (how far below Tg the operating temperature is), the membrane thickness (thin-film composite membranes age faster because the surface-to-volume ratio is larger), and the gas environment (condensable gases like CO2 can accelerate or retard aging depending on concentration). Predicting the performance trajectory over 3-5 years of operation requires understanding the molecular mechanism of aging.
Molecular simulation models physical aging by running long-duration simulations of the polymer packing and tracking the free volume evolution. Accelerated molecular dynamics techniques (temperature-accelerated dynamics, metadynamics) extend the accessible timescales from nanoseconds to the microsecond-millisecond range, capturing the slow relaxation processes that drive aging. The simulation predicts the aging rate and its dependence on conditions, enabling projection of membrane performance over the operational lifetime.
Designing New Membrane Materials
The Robeson upper bound is the empirical tradeoff between permeability and selectivity for polymer membranes: more permeable polymers tend to be less selective, and vice versa. Exceeding the upper bound (achieving both high permeability and high selectivity) requires designing polymer structures with molecular-level precision.
Polymers of intrinsic microporosity (PIMs), thermally rearranged (TR) polymers, and mixed-matrix membranes (MMMs) with metal-organic framework (MOF) fillers are current approaches to exceeding the upper bound. Each creates a specific free volume architecture that enables selective gas transport. PIMs have rigid, contorted backbones that prevent efficient packing, creating permanent microporosity. TR polymers have ladder-like structures with well-defined pore sizes. MOF fillers add crystalline micropores with apertures tuned to discriminate between gas molecules by size.
Molecular simulation guides the design of these advanced membranes by predicting the free volume distribution from the polymer or MOF structure, the gas solubility and diffusivity in each material, the plasticization resistance (how much CO2 the material can absorb before swelling), and the aging rate (how quickly the excess free volume relaxes). For MOF-based MMMs, the simulation also predicts the MOF-polymer interfacial compatibility, since poor compatibility creates non-selective voids at the interface that destroy selectivity.
The design loop is rapid: propose a new polymer structure (backbone chemistry, side groups, crosslink density), simulate the gas transport properties, evaluate against target specifications, and iterate. Each candidate takes 2-5 days of computation. Synthesizing and testing the same candidate experimentally takes 2-6 months. The simulation enables exploration of hundreds of candidates in the time it takes to test five experimentally.
Process Simulation with Realistic Membrane Properties
Membrane process design requires solving the mass balance equations across the membrane module, accounting for the pressure profile along the module length, the changing gas composition (and therefore driving force) from feed to retentate, the temperature profile (permeation is exothermic for condensable gases), and the concentration dependence of permeability (plasticization, competitive sorption).
Standard process simulators use constant permeabilities. This is acceptable for low-pressure, low-concentration applications where the membrane operates well below plasticization onset. For high-pressure natural gas processing, landfill gas upgrading, or biogas purification (where CO2 partial pressures routinely exceed plasticization onset), constant permeabilities give wrong answers.
Multi-scale simulation couples the molecular-level transport model (which provides permeability as a function of local conditions) with the continuum-scale process model (which solves the module mass balance). The permeability is updated at each point along the module length based on the local gas composition and pressure. This captures the fact that the membrane operates in different regimes along the module: near the feed end, high CO2 partial pressure may cause plasticization; near the retentate end, low CO2 partial pressure gives ideal separation. A constant-permeability model averages over these regimes and gets the overall performance wrong.
Membrane Module and Process Optimization
With accurate, condition-dependent membrane properties, the process model can optimize the complete membrane system. Key design variables include membrane area and configuration (single-stage, two-stage with recycle, two-stage with different membranes), operating pressure (higher pressure increases driving force but worsens plasticization), temperature (affects both permeability and selectivity), permeate pressure (lower improves recovery but increases recompression cost), and stage cut (fraction of feed that permeates, which determines product purity and recovery tradeoff).
The optimization finds the minimum-cost system that meets the product specification. For natural gas processing, the cost function includes membrane capital cost ($50-200/m2 installed), compressor capital and operating cost, methane loss (methane that permeates with CO2 is lost revenue), and membrane replacement cost (typically every 3-5 years).
A gas processing company used simulation-guided membrane process design for a 50 MMscfd natural gas sweetening application. The simulation predicted that a 2-stage system with interstage heating (to reverse plasticization between stages) would achieve 95% methane recovery at 2% CO2, compared to 88% recovery with the original isothermal 2-stage design. The 7% improvement in methane recovery at this flow rate is worth $2.4M/year. The interstage heater cost $180K installed.
Economics of simulation-guided membrane process design:
- Pilot testing for membrane selection (5 membranes x 10 conditions): $150K-400K, 3-6 months
- Simulation-guided design (50 virtual + 3 pilot conditions): $30K-100K, 4-6 weeks
- Oversized membrane system (20-40% excess area from constant-permeability model): $200K-2M unnecessary capex
- Methane recovery improvement from accurate process model: 3-7%, worth $500K-5M/year
- Membrane lifetime prediction: avoids premature replacement ($100K-500K per replacement event)
For companies deploying membrane systems across multiple sites (where each site has different gas composition, pressure, and temperature), simulation reduces the per-site engineering cost from $200K+ (custom pilot testing) to $30K-50K (simulation with minimal validation). The physics does not change between sites; only the conditions change, and the simulation handles any conditions from first principles. Explore MolSim for membrane and separation process simulation, or discuss your membrane process design challenges with our engineering team.