Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parametric analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Benchmark Exercise Report for Experimental Study of Bubble Scrubbing in Sodium Pool

Mechanistic source term (MST) analyses are likely to be an important part of advanced reactor licensing applications. For the purpose, an MST analysis code applicable to newly introduced advanced reactors, called SRT (Simplified Radionuclide Transport) code, has been developed by Argonne National Laboratory. SRT can track overall behaviors of radionuclides especially in metal fuel-based sodium fast reactors (SFRs) and microreactors. Throughout the simulation, migration inside fuel pins before failure, interaction with coolant (for SFR), removal/leakage in cover gas and containment (or confinement), and environmental dose impacts are considered alongside radioactive decay for short-lived nuclides. Among the postulated process, pool scrubbing phenomenon, especially under sodium pool condition, has been identified as high importance with limited supportive data. The phenomenon plays a crucial role in assessing the degree of radiological impacts as radioactive aerosols or vapors are efficiently and effectively removed during the process. To provide validation basis for SRT in assessing pool scrubbing performance inside sodium pools, the University of Wisconsin-Madison performed tests including extensive parametric effects. Separate effect tests were conducted to directly evaluate the SRT models and to estimate degree of contribution by each contributing factor. Specifically, bubble size, aerosol size, aerosol density, aerosol concentration, pool depth, system temperature, and bubble swarm effects were considered. According to the parametric effects, decontamination performance enhances with decreasing bubble size, large density, and deeper pool height. Aerosol concentration provides no effect for the whole range of interest, and pool temperature variation shows minor effects under the considered temperature condition. When multiple bubbles are injected generating a bubble swarm condition, DF performance further enhances by bubble interactions and turbulence characteristics. The measurement shows the exceptional importance of aerosol size range considered, with the lowest decontamination, where most radionuclides are expected to escape.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Heat Pump Retrofits for Central Plant Hydronic Heating Systems: A Software Toolkit for Screening and Design

Retrofitting existing central plants with high-efficiency heat pump technologies can play a crucial role in achieving long-term planning goals. Modern heat pump technologies are able to use waste heat recovery to meet a building's heating demand, but there is a lack of accessible tools designed for non-HVAC experts, such as building owners, to quickly and easily conduct what-if analysis, e.g., estimating retrofit costs and payback period for their partial or full equipment replacement. This paper introduces an open-source software toolkit designed to facilitate the initial screening and decision-making of heat pump retrofits in existing central plants using a building's yearly load profile from metered or utility bill data. The toolkit evaluates the technical and economic viability of replacing traditional central plant equipment with various options including water-to-water or air-to-water heat pumps, which can provide efficient and lower-cost heating and cooling. It allows users to compare current central plant configurations with retrofit scenarios, assessing energy consumption, life-cycle costs, and environmental impact. The toolkit offers (1) a web-based tool designed for user-friendly access by a broad audience and (2) Python-based source code for researchers and engineers conducting parametric studies and design parameter optimization. The toolkit compares a typical central plant configuration to a configuration that uses a heat pump to supply hydronic heating and cooling. The output metrics include energy consumption and output of each equipment, life-cycle cost analyses and metrics, and environmental impact of the system.

Excell, L↗
DOE OSTI2022

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: measurement of the growth rate of structure from the small-scale clustering of the luminous red galaxy sample

ABSTRACT We measure the small-scale clustering of the Data Release 16 extended Baryon Oscillation Spectroscopic Survey Luminous Red Galaxy sample, corrected for fibre-collisions using Pairwise Inverse Probability weights, which give unbiased clustering measurements on all scales. We fit to the monopole and quadrupole moments and to the projected correlation function over the separation range $7-60\, h^{-1}{\rm Mpc}$ with a model based on the aemulus cosmological emulator to measure the growth rate of cosmic structure, parametrized by fσ8. We obtain a measurement of fσ8(z = 0.737) = 0.408 ± 0.038, which is 1.4σ lower than the value expected from 2018 Planck data for a flat ΛCDM model, and is more consistent with recent weak-lensing measurements. The level of precision achieved is 1.7 times better than more standard measurements made using only the large-scale modes of the same sample. We also fit to the data using the full range of scales $0.1\text{--}60\, h^{-1}{\rm Mpc}$ modelled by the aemulus cosmological emulator and find a 4.5σ tension in the amplitude of the halo velocity field with the Planck + ΛCDM model, driven by a mismatch on the non-linear scales. This may not be cosmological in origin, and could be due to a breakdown in the Halo Occupation Distribution model used in the emulator. Finally, we perform a robust analysis of possible sources of systematics, including the effects of redshift uncertainty and incompleteness due to target selection that were not included in previous analyses fitting to clustering measurements on small scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Nature of LoBAL QSOs. II. HST/WFC3 Observations Reveal Host Galaxies Dominated by Mergers

Low-ionization broad absorption line QSOs (LoBALs) are suspected to be merging systems in which extreme, active galactic nucleus-driven outflows have been triggered. Whether or not LoBALs are uniquely associated with mergers, however, has yet to be established. To characterize the morphologies of LoBALs, we present the first high-resolution morphological analysis of a volume-limited sample of 22 Sloan Digital Sky Survey (SDSS)-selected LoBALs at 0.5 < z < 0.6 from Hubble Space Telescope Wide Field Camera 3 observations. Host galaxies are resolved in 86% of the systems in F125W, which is sensitive to old stellar populations, while only 18% are detected in F475W, which traces young, unobscured stellar populations. Signs of recent or ongoing tidal interaction are present in 45%–64% of the hosts, including double nuclei, tidal tails, bridges, plumes, shells, and extended debris. Ongoing interaction with a companion is apparent in 27%-41% of the LoBALs, with as much as 1/3 of the sample representing late-stage mergers at projected nuclear separations <10 kpc. Detailed surface brightness modeling indicates that 41% of the hosts are bulge dominated while only 18% are disks. We discuss trends in various properties as a function of merger stage and parametric morphology. Notably, mergers are associated with slower, dustier winds than those seen in undisturbed/unresolved hosts. Our results favor an evolutionary scenario in which quasar-level accretion during various merger stages is associated with the observed outflows in low-z LoBALs. We discuss differences between LoBALs and FeLoBALs and show that selection via the traditional balnicity index would have excluded all but one of the mergers.

Quasars↗
DOE OSTI2024

Parametric description of intermittent probability distribution functions in solar wind and magnetohydrodynamic turbulence

ABSTRACT In this work, we find empirical evidence that the scale-dependent statistical properties of solar wind and magnetohydrodynamic (MHD) turbulence can be described in terms of a family of parametric probability distribution functions (PDFs) known as Normal Inverse Gaussian (NIG). Understanding these PDFs is one of the most important goals in turbulence theory, as they are inherently connected to the intermittent properties of solar wind turbulence. We investigate the properties of PDFs of Elsasser increments based on a large statistical sample from solar wind observations and high-resolution numerical simulations of MHD turbulence. In order to measure the PDFs and their corresponding properties, three experiments are presented: fast and slow solar wind for experimental data and a simulation of reduced MHD (RMHD) turbulence. Conditional statistics on a 23-yr-long sample of WIND data near 1 au and high-resolution pseudo-spectral simulation of steadily driven RMHD turbulence on a $2048^3$ mesh are used to construct scale-dependent PDFs. The empirical PDFs are fitted to NIG distributions, which depend on four free parameters. Our analysis shows that NIG distributions accurately capture the evolution of the PDFs, with scale-dependent parameters, from large scales characterized by a Gaussian distribution, turning to exponential tails within the inertial range and stretched exponentials at dissipative scales. We also show that empirically-measured NIG parameters exhibit well-defined scaling properties that are similar across the three empirical data sets, which may be indicative of universal behaviour.

Astronomy & Astrophysics↗
DOE OSTI2022

Sensitivity Analysis of Wind and Turbulence Predictions With Mesoscale‐Coupled Large Eddy Simulations Using Ensemble Machine Learning

Abstract Coupling between mesoscale models and large‐eddy simulation (LES) models is increasingly used to more realistically represent the wide range of scales of atmospheric motions affecting boundary layer winds and turbulence that need to be simulated accurately for applications such as wind energy. However, such mesoscale‐to‐microscale coupled modeling frameworks are potentially affected by a large number of uncertain closure parameters. Here, we investigate the sensitivity associated with six closure parameters related to a 1.5‐order subgrid‐scale turbulence closure for an ensemble of mesoscale‐coupled LES. The simulations are performed using the Weather Research and Forecasting model nested from horizontal resolutions of greater than a kilometer down to tens of meters. Closure parameters are varied to generate perturbed parameter ensembles for two case studies of highly sheared, convective boundary layers observed in the Columbia Basin of Oregon and Washington during the Second Wind Forecast Improvement Project. Machine learning algorithms are used to explore the sensitivity of LES predictions, considering the effects of the perturbed physical parameters alongside categorical factors such as the case study identity, measurement location, and LES resolution. For the conditions we examine, a single parameter, the eddy viscosity coefficient, is the dominant source of parametric sensitivity and its importance is comparable to the categorical factors for several of the simulation response variables we examine.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Spectral Mesh-Free Quadrature for Planar Regions Bounded by Rational Parametric Curves

This article presents spectral, mesh-free, Green’s theorem-based numerical quadrature schemes for integrating functions over planar regions bounded by rational parametric curves. Our algorithm proceeds in two steps: (1) We first find intermediate quadrature rules for line integrals along the region’s boundary curves corresponding to Green’s theorem. (2) We then use a high-order quadrature rule to compute the numerical antiderivative of the integrand along a coordinate axis, which is used to evaluate the Green’s theorem line integral. We present two methods to compute the intermediate quadrature rule. The first is spectrally accurate (it converges faster than any algebraic order with respect to number of quadrature points) and is relatively easy to implement, but has no guarantee of polynomial exactness. The second guarantees exactness for polynomial integrands up to a pre-specified degree with an a priori-known number of quadrature points and retains the convergence properties of the first, but is slightly more complicated. The quadrature schemes have applications to computation of geometric moments, immersogeometric analysis, conservative field transfer between high-order meshes, and initialization of multi-material simulations with rational geometry. We compare the quadrature schemes produced using our method to other methods in the literature and show that they are much more efficient both in terms of number of quadrature points and computational time. We provide an open-source implementation of the algorithm in MATLAB.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Benders Decomposition Using Graph Modeling and Multi-Parametric Programming

Benders decomposition is a widely used method for solving large and structured optimization problems, but its performance is affected by the repeated solution of subproblems. We propose a flexible and modular algorithmic framework for accelerating Benders decomposition. Specifically, we express the problem structure by using a graph-theoretic modeling abstraction in which nodes represent optimization subproblems and edges represent connectivity between subproblems. A key innovation of our approach is that we embed multiparametric programming (mp) surrogates for node subproblems, which maps the exact analytical map of the subproblem solution space. The use of mp surrogates allows us to replace subproblem solves with fast look-ups and function evaluations for primal and dual variables during the iterative Benders process. We formally show the equivalence between classical Benders cuts and those derived from the mp solution. We implement our framework in the open-source PlasmoBenders.jl software package. To demonstrate the capabilities of the proposed framework, we apply it to a two-stage stochastic programming problem, which aims to make optimal capacity expansion decisions under market uncertainty. We evaluate both single-cut and multicut variants of Benders decomposition and show that the use of mp surrogates achieves substantial speedups in subproblem solve time, while preserving the convergence guarantees of Benders decomposition. We highlight advantages in solution analysis and interpretability that is enabled by mp critical region tracking; specifically, we show that these reveal how decisions evolve geometrically across the Benders search. Our results aim to demonstrate that combining surrogate modeling with graph modeling offers a promising and extensible foundation for structure-exploiting decomposition. In addition, by decomposing the problem into more tractable subproblems, the proposed approach also aims to overcome scalability issues of mp. Finally, the use of mp surrogates provides a unifying and modular optimization framework that enables the representation of heterogeneous node subproblems as modeling objects with a homogeneous structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Radiative and continuum dampings of reversed shear Alfvén eigenmodes and perturbative analysis limitations for tokamaks

A careful theoretical analysis of the excitation of Alfvén eigenmodes (AEs), such as TAE (toroidicity-induced AE) and RSAE (reversed shear AE), by superalfvenic energetic particles is required for reliable predictions of energetic ion relaxation in present day fusion experiments. This includes the evaluation of different AE damping mechanisms including radiative and continuum dampings which are the focus of this study. A recent comprehensive benchmark of different eigenmode solvers including gyrokinetic, gyrofluid and hybrid magenetohydrodynamics (MHD) has shown that employed models may have deficiencies when addressing some of them (Taimourzadeh et al., Nucl. Fusion, vol. 59, 2019, 066006). Here, in this paper, we are studying the radiative and continuum dampings of RSAEs in details which were missing in hybrid NOVA/NOVA-C calculations to prepare a NOVA-C package with a substantial upgrade. Both dampings require the finite Larmor radius (FLR) corrections to AE mode structures to be accounted for. Accurately calculating different damping rates and understanding their parametric dependencies, we resolve the limitation coming out of the perturbative approach. In particular, here, the radiative damping is included perturbatively, whereas the continuum damping is computed non-perturbatively. Our comparison leads to the conclusion that the non-perturbative treatment of the unstable RSAE modes is needed to find the agreement with the gyrokinetic calculations. We expect that the RSAE mode structure modification plays a dominant role in determining the RSAE stability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Data-driven recovery of hidden physics in reduced order modeling of fluid flows

In this article, we introduce a modular hybrid analysis and modeling (HAM) approach to account for hidden physics in reduced order modeling (ROM) of parameterized systems relevant to fluid dynamics. The hybrid ROM framework is based on using first principles to model the known physics in conjunction with utilizing the data-driven machine learning tools to model the remaining residual that is hidden in data. This framework employs proper orthogonal decomposition as a compression tool to construct orthonormal bases and Galerkin projection (GP) as a model to build the dynamical core of the system. Our proposed methodology hence compensates structural or epistemic uncertainties in models and utilizes the observed data snapshots to compute true modal coefficients spanned by these bases. The GP model is then corrected at every time step with a data-driven rectification using a long short-term memory (LSTM) neural network architecture to incorporate hidden physics. A Grassmannian manifold approach is also adopted for interpolating basis functions to unseen parametric conditions. The control parameter governing the system's behavior is thus implicitly considered through true modal coefficients as input features to the LSTM network. The effectiveness of the HAM approach is then discussed through illustrative examples that are generated synthetically to take hidden physics into account. Furthermore, our approach thus provides insights addressing a fundamental limitation of the physics-based models when the governing equations are incomplete to represent underlying physical processes.

42 ENGINEERING↗
DOE OSTI2020

Weak structure functions in $\nu_l-N$ and $\nu_l-A$ scattering with nonperturbative and higher order perturbative QCD effects

We study the effect of various perturbative and nonperturbative QCD corrections on the free nucleon structure functions [FiNWI(x,Q2);i=1–3] and their implications in the determination of nuclear structure functions. The evaluation of the nucleon structure functions has been performed by using the MMHT 2014 parton distribution functions (PDFs) parametrization, and the target mass correction (TMC) and higher twist (HT) effects are incorporated following the works of Schienbein et al. and Dasgupta et al., respectively. These nucleon structure functions are taken as input in the determination of nuclear structure functions. The numerical calculations for the νl/ν¯l-A deep inelastic scattering (DIS) process have been performed by incorporating the nuclear medium effects like Fermi motion, binding energy, nucleon correlations, mesonic contributions, shadowing, and antishadowing in several nuclear targets such as carbon, polystyrene scintillator, iron, and lead, which are being used in MINERνA, and in argon nuclei, which is relevant for the ArgoNeuT and DUNE experiments. The differential scattering cross sections d2σAWIdxdy and (dσAWIdx/dσCHWIdx) have also been studied in the kinematic region of the MINERνA experiment. The theoretical results are compared with the recent experimental data of MINERνA and the earlier data of the NuTeV, CCFR, CDHSW, and CHORUS Collaborations. Moreover, a comparative analysis of the present results for the ratio (dσAWIdx/dσCHWIdx), and the results from the Monte Carlo (MC) generator GENIE and other phenomenological models of Bodek and Yang, and Cloet et al., has been performed in the context of the MINERνA experiment. The predictions have also been made for the ν¯l-A cross section relevant for the MINERνA experiment.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Chiral-odd generalized parton distributions in the large-𝑁 𝑐 limit of QCD: Spin-flavor structure, polynomiality, and sum rules

We study the nonperturbative properties of the nucleon’s chiral-odd generalized parton distributions (transversity GPDs) in the large-𝑁 𝑐 limit of QCD. This includes the parametric ordering of the spin-flavor components, the polynomiality property of the moments, and the sum rules connecting the GPDs with the tensor form factors. A multipole expansion in the transverse momentum transfer is used to enumerate and interpret the structures in the nucleon matrix element of the chiral-odd partonic operator, including monopole, dipole and quadrupole terms. The 1/𝑁 𝑐 expansion of the GPDs is performed using the abstract mean-field picture of baryons in the large-𝑁 𝑐 limit and its symmetries. We derive a large-𝑁 𝑐 relation between the flavor-nonsinglet GPDs 𝐸$^{𝑢−𝑑}_𝑇$ and $\tilde{𝐻}^{𝑢−𝑑}_𝑇$ and test it with recent lattice QCD results. We show that the polynomiality property and sum rules of the GPDs are fulfilled with the restricted realization of translational and rotational invariance in the mean-field picture. The results provide a basis for the phenomenological analysis of chiral-odd GPDs and hard exclusive processes in the large-𝑁 𝑐 limit, and for calculations in specific dynamical models.

generalized parton distributions↗
DOE OSTI2024

Effects of wave damping and finite perpendicular scale on three-dimensional Alfvén wave parametric decay in low-beta plasmas

Shear Alfvén wave parametric decay instability (PDI) provides a potential path toward significant wave dissipation and plasma heating. However, fundamental questions regarding how PDI is excited in a realistic three-dimensional (3D) open system and how the finite perpendicular wave scale—as found in both laboratory and space plasmas—affects the excitation remain poorly understood. Here, we present the first 3D, open-boundary, hybrid kinetic-fluid simulations of kinetic Alfvén wave PDI in low-beta plasmas. Key findings are that the PDI excitation is strongly limited by the wave damping present, including electron–ion collisional damping (represented by a constant resistivity) and geometrical attenuation associated with the finite-scale Alfvén wave, and ion Landau damping of the child acoustic wave. The perpendicular wave scale alone, however, plays no discernible role: waves of different perpendicular scales exhibit similar instability excitation as long as the magnitude of the parallel ponderomotive force remains unchanged. These findings are corroborated by theoretical analysis and estimates. This new understanding of 3D kinetic Alfvén wave PDI physics is essential for laboratory study of the basic plasma process and may also aid future evaluation of the relevance/role of PDI in low-beta space plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Optical Characterization of the Keck Array and BICEP3 CMB Polarimeters from 2016 to 2019

The BICEP/Keck experiment (BK) is a series of small-aperture refracting telescopes observing degree-scale cosmic microwave background (CMB) polarization from the South Pole in search of a primordial B-mode signature. This B-mode signal arises from primordial gravitational waves interacting with the CMB and has amplitude parametrized by the tensor-to-scalar ratio r. Since 2016, BICEP3 and the Keck Array have been observing with 4800 total antenna-coupled transition-edge sensor detectors, with frequency bands spanning 95, 150, 220, and 270 GHz. Here we present the optical performance of these receivers from 2016 to 2019, including far-field beams measured in situ with an improved chopped thermal source and instrument spectral response measured with a field-deployable Fourier transform spectrometer. As a pair differencing experiment, an important systematic that must be controlled is the differential beam response between the co-located, orthogonally polarized detectors. We generate per-detector far-field beam maps and the corresponding differential beam mismatch that is used to estimate the temperature-to-polarization leakage in our CMB maps and to give feedback on detector and optics fabrication. The differential beam parameters presented here were estimated using improved low-level beam map analysis techniques, including efficient removal of non-Gaussian noise as well as improved spatial masking. These techniques help minimize systematic uncertainty in the beam analysis, with the goal of constraining the bias on r induced by temperature-to-polarization leakage to be subdominant to the statistical uncertainty. This is essential as we progress to higher detector counts in the next generation of CMB experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Physics-informed machine learning with differentiable programming for heterogeneous underground reservoir pressure management

Abstract Avoiding over-pressurization in subsurface reservoirs is critical for applications like CO $$_2$$ 2 sequestration and wastewater injection. Managing the pressures by controlling injection/extraction are challenging because of complex heterogeneity in the subsurface. The heterogeneity typically requires high-fidelity physics-based models to make predictions on CO $$_2$$ 2 fate. Furthermore, characterizing the heterogeneity accurately is fraught with parametric uncertainty. Accounting for both, heterogeneity and uncertainty, makes this a computationally-intensive problem challenging for current reservoir simulators. To tackle this, we use differentiable programming with a full-physics model and machine learning to determine the fluid extraction rates that prevent over-pressurization at critical reservoir locations. We use DPFEHM framework, which has trustworthy physics based on the standard two-point flux finite volume discretization and is also automatically differentiable like machine learning models. Our physics-informed machine learning framework uses convolutional neural networks to learn an appropriate extraction rate based on the permeability field. We also perform a hyperparameter search to improve the model’s accuracy. Training and testing scenarios are executed to evaluate the feasibility of using physics-informed machine learning to manage reservoir pressures. We constructed and tested a sufficiently accurate simulator that is 400 000 times faster than the underlying physics-based simulator, allowing for near real-time analysis and robust uncertainty quantification.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Cosmology with the Roman Space Telescope – Synergies with CMB lensing

ABSTRACT We explore synergies between the Nancy Grace Roman Space Telescope and CMB lensing data to constrain dark energy and modified gravity scenarios. A simulated likelihood analysis of the galaxy clustering and weak lensing data from the Roman Space Telescope High Latitude Survey combined with CMB lensing data from the Simons Observatory is undertaken, marginalizing over important astrophysical effects and calibration uncertainties. Included in the modelling are the effects of baryons on small-scale clustering, scale-dependent growth suppression by neutrinos, as well as uncertainties in the galaxy clustering biases, in the intrinsic alignment contributions to the lensing signal, in the redshift distributions, and in the galaxy shape calibration. The addition of CMB lensing roughly doubles the dark energy figure-of-merit from Roman photometric survey data alone, varying from a factor of 1.7 to 2.4 improvement depending on the particular Roman survey configuration. Alternatively, the inclusion of CMB lensing information can compensate for uncertainties in the Roman galaxy shape calibration if it falls below the design goals. Furthermore, we report the first forecast of Roman constraints on a model-independent structure growth, parametrized by σ8(z), and on the Hu-Sawicki f(R) gravity as well as an improved forecast of the phenomenological (Σ0, μ0) model. We find that CMB lensing plays a crucial role in constraining σ8(z) at z > 2, with percent-level constraints forecasted out to z = 4. CMB lensing information does not improve constraints on the f(R) model substantially. It does, however, increase the (Σ0, μ0) figure-of-merit by a factor of about 1.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Numerical modeling and performance analysis of a membrane-based air dehumidifier using ionic liquid desiccant

Membrane-based Liquid Desiccant Air Dehumidifier (MLDAD) has the potential to overcome the shortcomings of the conventional open-tower liquid desiccant air dehumidifier. The MLDAD uses membranes with high water vapor permeability to separate air stream from the liquid desiccant and thus eliminates carry-over and the resulting corrosion issues. Recently, the ionic liquid desiccant becomes a promising alternative to conventional liquid desiccant. The ionic liquid desiccant has a large potential of dehumidification, which is also non-corrosive to metals and non-crystallizable. A two-dimensional numerical heat and mass transfer model of the MLDAD using a recently identified ionic liquid desiccant, [EMIM]OAc, is presented here. This model can simulate the performance of the MLDAD with various designs, including different selections of liquid desiccant and membrane materials, dimensions of the MLDAD, and flow patterns. For porous membranes, this model accounts for several micro-scale mass transfer mechanisms of vapor transportation across the membrane. It can also model the mass transfer performance of dehumidifiers using non-porous membranes based on the experimentally measured membrane permeance. The numerical model was validated against the performance data available from literature and experimental tests. The maximum discrepancy of the latent effectiveness between the measured and the model-predicted results was about 6% when porous membranes were used. A parametric study was conducted with the numerical model. Finally, the results indicated that the membrane permeability, air path height, and solution to air mass ratio were the most critical parameters, determining the dehumidification performance of the MLDAD.

42 ENGINEERING↗
DOE OSTI2026

Dual-loop Solvent-based CCS for Net Negative CO 2 Emissions with Lower Cost

This final technical report details the successful design, construction, and operational validation of an innovative dual-loop CO 2 capture technology designed to achieve deep decarbonization (>99%) from Natural Gas Combined Cycle (NGCC) power plants that results in the electricity with carbon intensity of approximate 42 kg CO 2 -eq/MWh, less than the electricity produced by solar PV. The integrated process couples a primary aqueous solvent (in this project, a water lean solvent – WLS) absorption loop for bulk CO 2 removal with a secondary potassium hydroxide (KOH) polishing loop featuring electrochemical regeneration. This architecture leverages the higher exergy efficiency of the primary loop while utilizing the fastest kinetic of the secondary loop to capture dilute residual CO 2 , achieving an overall capture rate of 99.9% and co-producing pure hydrogen after moisture being condensed and dehydrated. Technical feasibility was established through a comprehensive 2,000-hour experimental campaign on a bench-scale fully-integrated unit (using 4” absorber and 4” stripper) with the feeding flue gas flowrate in the range of 8-20 cfm, confirming the attainment of Technology Readiness Level (TRL) 4. The project executed extensive parametric testing followed by 1,000 hours of continuous steady-state testing, demonstrating exceptional process stability with the electrochemical regenerator exhibiting less than a 10% reduction in electrical conductivity over the duration of the campaign. Operational characterization gathering on the bench unit revealed distinct energy profiles for the hybrid system: the primary loop required approximately 280 kJ mol -1 for bulk removal, while the polishing loop required approximately 1,600 kJ mol -1 specifically when reducing dilute CO 2 concentrations from ~740 ppm down to <50 ppm. (Please note those energy values/numbers can only be viewed as relative relationship and should not be extrapolated as absolute values required for CO 2 capture). Furthermore, dynamic testing validated the system’s flexibility for utility applications, demonstrating a rapid process response time of <30 minutes to changes in flue gas flowrate. Emission monitoring confirmed that the dual-loop architecture effectively mitigates solvent losses, utilizing a water wash to remove entrained aerosols to <1 ppm. The Techno-Economic Analysis (TEA) indicates a cost of capture of $\$$59.3/tonne and a Levelized Cost of Electricity (LCOE) of 71.7 $\$$/MWh at the overall capture efficiency of 99.8% of total carbon in the flue gas stream, with sensitivity analysis identifying an economic optimum when the primary loop captures 97% of the total CO 2 . A Life Cycle Assessment (LCA) confirms the technology’s potential for net-negative emissions, determining a Global Warming Potential (GWP) of 52 kg CO 2 -eq/MWh—significantly lower than the baseline—which further decreases to 42 kg CO 2 -eq/MWh when crediting the displacement of conventional hydrogen production.

03 NATURAL GAS↗