The research, led by Keju An and colleagues at the National Laboratory of the Rockies in Golden, Colorado, tackles a stubborn computational bottleneck. Detailed models of adsorption systems—devices where solid sorbent materials grab CO2 molecules from air and then release them when heated and placed under vacuum—are so computationally expensive that simulating a full year of hour-by-hour operation was previously impractical. Without year-long simulations, engineers could not properly explore how machines should behave when weather and grid conditions change around them.

The team’s solution was a physics-informed reduced-order model: a stripped-down but physically faithful mathematical description of a temperature-vacuum swing adsorption (TVSA) cycle, the six-stage process of vacuum extraction, heating, desorption, cooling, pressurization, and adsorption. Crucially, the model retains the phenomena that matter most. It tracks CO2 and water co-adsorption on the sorbent, uses a two-pool formulation in which some CO2 binds quickly to accessible sites while a slower pool represents stubbornly retained molecules, and captures incomplete regeneration—the fact that each cycle never fully clears the sorbent, leaving residual loading that carries over as a memory into the next cycle.

That cycle-to-cycle memory is the heart of the study’s argument. Because the sorbent’s state at the start of each cycle depends on everything that happened before, performance cannot be predicted from any single steady-state snapshot. The loading states act as memory variables that determine the working capacity available in the next cycle, the humidity-dependent adsorption behavior, and the energy required for regeneration. The model tracks six coupled states—fast and slow CO2 pools, free and bound water, gas-phase CO2 concentration, and bed temperature—propagating them forward step by step through the year.

To make the framework fast enough for serious optimization work, the researchers generated roughly 120,000 cycle-resolved simulations using Latin hypercube sampling across environmental conditions, operating strategies, system designs, and initial sorbent states. Each simulated cycle recorded its ambient conditions, control settings, internal sorbent states, and resulting performance. This dataset then trained a compact neural network—a single hidden layer of 39 neurons—that reproduces the physics model’s state transitions at a tiny fraction of the computational cost, creating a surrogate fast enough to embed inside a model predictive control framework.

Before trusting the surrogate, the team benchmarked the physics model against published experimental data from Sendi and colleagues. Across twelve benchmark cases spanning typical temperature and humidity conditions, the reduced-order model reproduced productivity with a mean absolute percentage error of 10.2 percent and electricity demand with an error of 11.4 percent. Deviations grew largest under extreme conditions—high temperature with low humidity, or very cold conditions—reflecting the simplifications inherent in a lumped, zero-dimensional formulation. But the model faithfully captured the direction and relative sensitivity of performance to environmental changes, which matters more for large-scale screening than exact point accuracy.

Sensitivity analysis revealed a starkly asymmetric performance landscape. Regeneration temperature dominated everything: raising it from 70 to 130 degrees Celsius cut annual productivity by roughly 10 tonnes of CO2 per cubic meter of bed while adding 0.7 to 0.8 megawatt-hours of electricity demand per tonne captured—a punishing double penalty from hotter, slower cycles. Humidity was the second most influential factor for energy, with relative humidity climbing from 0.1 to 0.9 adding roughly 0.5 megawatt-hours per tonne, a direct consequence of water co-adsorbing alongside CO2 and demanding extra thermal energy during regeneration.

The real payoff came when the neural-network surrogate was handed to the optimization controller and run against hourly 2025 weather, electricity price, and grid carbon intensity data for two contrasting American locations: humid Atoka, Texas and dry, cooler Madison, Wyoming. A baseline machine running fixed settings completed about 1,094 cycles per year. The adaptive controller pushed that to roughly 1,449 cycles in Atoka and 1,379 in Madison by shifting energy-hungry regeneration steps into hours when electricity was cheap and clean, and exploiting favorable humidity windows for capture.

The results were striking. In Atoka, electricity demand fell 18.8 percent while annual CO2 productivity rose 25.5 percent. In Madison, energy dropped 8.9 percent and productivity climbed 22.4 percent. The most dramatic gains appeared in carbon-adjusted net productivity—capture minus the emissions embedded in the electricity consumed—which nearly doubled at both sites, rising 95.2 percent in Texas and 94.4 percent in Wyoming. Notably, the economics did not always align with the climate math: electricity cost fell slightly in Atoka but rose 3.7 percent in Madison, showing that carbon-aware control sometimes accepts higher operating costs to harvest cleaner power.

The authors are careful about scope. Their model is a control-oriented screening tool, not a plant-specific predictor, and because the surrogate learns from the physics model, its uncertainties propagate into the optimization results—meaning the reported gains are relative, model-based estimates rather than guarantees. Still, the framework offers something the field has lacked: a scalable path from adsorption physics through fast dynamic simulation to real-world control, demonstrating that the biggest untapped efficiency reserve in direct air capture may not be new materials or bigger machines, but smarter timing.

Subject of Research: Dynamic modeling and neural-network-based optimization of solid-sorbent direct air capture systems under varying weather and electricity conditions

Article Title: Cycle-resolved dynamic modeling of atmospheric CO2 separation systems under time-varying ambient and electricity conditions

Article References: An, K., Wald, D., Chen, K., & King, R. (2026). Cycle-resolved dynamic modeling of atmospheric CO2 separation systems under time-varying ambient and electricity conditions. Frontiers in Environmental Engineering, 5, Article 1881554. https://doi.org/10.3389/fenve.2026.1881554

Image Credits: AI Generated

DOI: 10.3389/fenve.2026.1881554

Keywords: direct air capture, carbon dioxide removal, temperature-vacuum swing adsorption, reduced-order model, neural network surrogate, model predictive control, humidity, electricity carbon intensity, energy optimization, solid sorbents, climate engineering, Cycle-resolved

The starting point for the new work is a Bernoulli-based description of the discharge along a representative streamline that runs from the gas bulk inside the vessel, through the outlet pipe, and out into the undisturbed ambient atmosphere. Because the vessel volume is far larger than the outlet cross section, the gas velocity in the bulk is negligible, and the unsteady acceleration along the streamline can be treated as quasi-stationary. What the formulation does keep, and what distinguishes it from many purely thermodynamic blowdown models, is the explicit treatment of irreversible pressure losses: local inertial losses at the vessel-to-pipe transition and at the discharge into ambient air, plus distributed pipe friction described by the Darcy-Weisbach equation with a friction factor obtained iteratively from the Colebrook relation. These losses are folded directly into an expression for the area-averaged outlet velocity, which in turn feeds a mass balance that yields the time-resolved pressure decay inside the vessel.

The first major extension concerns choked flow. When the ratio of ambient pressure to vessel pressure falls below the critical pressure ratio, which depends on the heat capacity ratio of the gas, the flow at the narrowest section reaches sonic conditions and the mass flow rate can no longer increase no matter how much the downstream pressure drops. The team uses the classical ideal-gas critical-pressure-ratio criterion as a regime-switching rule: at every time step of the numerical integration, the instantaneous vessel pressure is compared with the critical pressure, and the choked-flow pressure-decay equation is applied as long as the vessel pressure remains above that threshold. Because friction and inertia reduce the effective mass flow, the model introduces a loss-corrected equivalent outlet velocity, obtained by dividing the local speed of sound by the effective loss coefficient. The authors are careful to stress that this is an engineering approximation: the model does not resolve the spatially varying compressible flow field and does not claim to locate the sonic section precisely, which in a real vessel-pipe combination may sit near a vena contracta just downstream of the inlet or at the pipe exit depending on geometry and friction.

The second extension brings thermal physics into the picture. As gas expands during discharge, it cools, and in applications such as hydrogen storage that cooling can push temperatures below limits tolerated by polymer liners and seals. The Dortmund team couples the discharge model to an energy balance in which heat flows from the vessel wall to the gas by convection, with the heat-transfer coefficient derived from empirical Nusselt correlations based on the Rayleigh number for natural convection. Natural convection is the appropriate regime here because numerical simulations of the same geometry show that gas velocities inside the vessel are very low except near the outlet, so buoyancy driven by temperature-induced density gradients dominates the internal heat exchange. A conservative sensitivity analysis confirmed that the stainless-steel vessel behaves as a thermal reservoir on the timescale of the experiments: even under deliberately pessimistic assumptions, the predicted wall-temperature change during a single discharge remained close to one kelvin, justifying a constant wall-temperature boundary condition.

Numerically, the friction-and-heat-transfer formulation is implemented as a sequential two-pass procedure rather than a fully implicit coupled solve. In the first pass, a preliminary pressure history is computed without pipe friction, and at each time step the Rayleigh number, Nusselt number, and heat-transfer coefficient are evaluated to determine the heat transferred from the wall to the gas. From the ratio of diabatic to isentropic enthalpy changes, a global heat-transfer polytropic exponent is extracted. In the second pass, the pressure history is recalculated with the full frictional losses, using the global thermal exponent in place of the isentropic one, while the Darcy friction factor is iterated to convergence at every step. A local effective polytropic exponent is then matched to the loss-corrected outlet velocity for post-processing the gas temperature. The whole scheme runs with an explicit fourth-order Runge-Kutta method at a time step of one ten-thousandth of a second, and a Richardson extrapolation showed the numerical error to be negligible compared with the experimental deviations.

Validation rested on 108 discharge experiments: air, carbon dioxide, and helium released from a 2.085 cubic meter stainless-steel vessel through a 25 millimeter outlet pipe, at twelve initial overpressures from 100 to 1500 millibar, with three repetitions per condition. A quick-opening valve, reproducibly actuated within 34 milliseconds, released the gas, while a phase Doppler anemometer measured outlet velocities on the jet axis one pipe diameter downstream. Fog droplets of ethylene glycol and water served as tracers, and only particle-size classes up to 11.5 micrometers were retained, because exit Stokes numbers for those classes were far below unity, confirming that the particles faithfully followed the rapidly accelerating gas. Statistical analysis with linear mixed models, including adjusted intraclass correlation coefficients and conditional coefficients of determination, showed that repetition effects were negligible and justified combining the three repetitions at each condition.

The agreement between model and experiment is striking for a reduced-order formulation. Across all gases and pressures, the mean relative deviation in discharge duration was just 3.60 percent, with experimental scatter between repetitions ranging from 0.1 to 2.6 percent. The measured discharge durations fell between the idealized isentropic and isothermal limits, exactly as the polytropic formulation with friction and heat transfer predicts. Maximum outlet velocities agreed even more closely for the heavier gases: deviations ranged from minus 1.7 to plus 0.9 percent for air and from minus 1.0 to plus 1.5 percent for carbon dioxide, comparable to the sub-one-percent nominal uncertainty of the optical technique itself. Notably, the maximum velocity occurs at the very start of the discharge, before wall heat transfer has had time to influence the gas state, which explains why varying the empirical Nusselt correlation parameters by plus or minus ten percent barely moved the predictions.

Helium exposed both the power and the limits of the approach. Its low molar mass means the vessel empties in seconds rather than tens of seconds, and its high speed of sound drives outlet velocities beyond the 433 meters per second ceiling of the phase Doppler anemometer at almost all tested pressures. Direct velocity validation was therefore confined to the lowest overpressure, where the model deviated by minus 6.5 percent. Helium discharge durations, however, were predicted over the entire pressure range, with the model tending to slightly underpredict, a trend the authors attribute to fixed-timescale effects such as valve opening and sensor response becoming relatively more important in short discharges, and to the sensitivity of light-gas discharge to the effective loss coefficient and the choked-to-subcritical transition. The team is transparent that helium velocities at higher pressures remain model predictions rather than validated measurements.

The practical payoff is a tool that delivers, in seconds on an ordinary computer, what would otherwise demand computationally expensive three-dimensional simulations: time-resolved vessel pressure, discharge duration, mass flow rate, and maximum outlet velocity, with choking, friction, and thermal interaction all included. That makes the model well suited to sizing pressure relief and venting systems, assessing mechanical loads and jet momentum, running parameter studies over vent-line geometries, and supplying boundary conditions for downstream gas dispersion and risk analyses. It is particularly timely for alternative energy carriers such as hydrogen, where low molar mass produces exactly the high-velocity, choked conditions the model now handles. The authors caution that for safety-critical final designs the model should be paired with appropriate safety factors and, where necessary, detailed simulation, and they identify improved natural-convection correlations for enclosures and spatially resolved temperature measurements as the key next steps for extending the thermal validation to larger vessels, longer discharge times, and different wall materials.

Subject of Research: Transient gas discharge modeling from pressurized vessels including choked flow and heat transfer

Article Title: Modeling and validation of transient gas outflows from pressurized vessels under choked and thermal effects

Article References: Fischer, M.-D., Ryll, F., Baier, S., Boettcher, K. E., & Dieguez-Alonso, A. (2026). Modeling and validation of transient gas outflows from pressurized vessels under choked and thermal effects. Results in Engineering, 32, Article 113271. https://doi.org/10.1016/j.rineng.2026.113271

Image Credits: AI Generated

DOI: 10.1016/j.rineng.2026.113271

Keywords: pressurized vessels, gas discharge, choked flow, heat transfer, blowdown modeling, phase Doppler anemometry, hydrogen safety, pressure relief, reduced-order model, natural convection, fluid dynamics, process safety