Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiscale systems”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Multiscale Interaction of Turbulence, Temperature, and Soot Formation: Measurements for Critical Assessments of Chemical Kinetics and Mechanisms

This program targeted an improved understanding of how turbulence affects soot formation chemistry (directly and indirectly) with a particular focus on detailing the multiscale coupling between flow turbulence, mixing, thermal transport, and soot formation kinetics in gas-phase reacting systems. In this project, simultaneous, quantitative multi-dimensional velocity, soot volume fraction, and gas-phase temperature measurements in turbulent non-premixed sooting flames were developed. These new diagnostics tools can be used to elucidate the relative effects of fluid kinematics, mixing, thermal processes, and turbulence/scalar time/length scales on soot formation chemistry and topology. Simultaneous multi-parameter measurements were acquired within a series of flames with different Reynolds numbers and spatial locations to investigate the inter-dependence of soot, temperature, and turbulence within the soot inception, growth, and oxidation region. Preliminary assessment of the results has been completed at the project end and detailed analysis is planned following the completed program.

02 PETROLEUM↗
DOE OSTI2023

AMORE-Isoprene v1.0: a new reduced mechanism for gas-phase isoprene oxidation

Abstract. Gas-phase oxidation of isoprene by ozone (O3) and the hydroxyl (OH) and nitrate (NO3) radicals significantly impacts tropospheric oxidant levels and secondary organic aerosol formation. The most comprehensive and up-to-date chemical mechanism for isoprene oxidation consists of several hundred species and over 800 reactions. Therefore, the computational expense of including the entire mechanism in large-scale atmospheric chemical transport models is usually prohibitive, and most models employ reduced isoprene mechanisms ranging in size from ∼ 10 to ∼ 200 species. We have developed a new reduced isoprene oxidation mechanism using a directed-graph path-based automated model reduction approach, with minimal manual adjustment of the output mechanism. The approach takes as inputs a full isoprene oxidation mechanism, the environmental parameter space, and a list of priority species which are protected from elimination during the reduction process. Our reduced mechanism, AMORE-Isoprene (where AMORE stands for Automated Model Reduction), consists of 12 species which are unique to the isoprene mechanism as well as 22 reactions. We demonstrate its performance in a box model in comparison with experimental data from the literature and other current isoprene oxidation mechanisms. AMORE-Isoprene's performance with respect to predicting the time evolution of isoprene oxidation products, including isoprene epoxydiols (IEPOX) and formaldehyde, is favorable compared with other similarly sized mechanisms. When AMORE-Isoprene is included in the Community Regional Atmospheric Chemistry Multiphase Mechanism 1.0 (CRACMM1AMORE) in the Community Multiscale Air Quality Model (CMAQ, v5.3.3), O3 and formaldehyde agreement with Environmental Protection Agency (EPA) Air Quality System observations is improved. O3 bias is reduced by 3.4 ppb under daytime conditions for O3 concentrations over 50 ppb. Formaldehyde bias is reduced by 0.26 ppb on average for all formaldehyde measurements compared with the base CRACMM1. There was no significant change in computation time between CRACMM1AMORE and the base CRACMM. AMORE-Isoprene shows a 35 % improvement in agreement between simulated IEPOX concentrations and chamber data over the base CRACMM1 mechanism when compared in the Framework for 0-D Atmospheric Modeling (F0AM) box model framework. This work demonstrates a new highly reduced isoprene mechanism and shows the potential value of automated model reduction for complex reaction systems.

58 GEOSCIENCES↗
DOE OSTI2024

Multiscale Dynamics of Reactive Fronts in the Subsurface

Understanding and predicting flow and reactive transport in rocks (i.e. geologic porous media) is critical to many technologies at the heart of the energy transition, including CO 2 sequestration and H2 storage. However, accurate modeling and prediction of these systems is very complex because physico-chemical processes that occur at very small spatial scales, i.e. in the pores of the rocks, can dramatically control the system performance at the field scale (km). For example, precipitation reactions at the pore-scale can lead to large permeability changes at the field scale and dramatically alter the migration of stored gases. Properly accounting for multi-scale coupling effects is critical to achieve predictivity and confidence in model outputs, which then can guide design and optimization at the system-scale. This can be achieved through the development of rigorous mathematical models that can appropriately account for fine-scale effects at the large scale. The final report of the Early Career award DE-SC0019075 “Multiscale dynamics of reactive fronts in the subsurface” summarizes the mathematical, numerical and experimental advancements to study and predict reactive transport in geologic porous media across scales.

58 GEOSCIENCES↗
DOE OSTI2025

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Analyzing the Multiscale Impacts of Implementing Energy-Efficient HVAC Improvements Through Energy Audits and Economic Input–Output Analysis

Abstract Heating, ventilation, and air-conditioning (HVAC) systems are usually an industry’s highest consumer of energy, most of which goes toward space cooling in buildings. Industrial energy-efficiency audits not only benefit manufacturers but also generate significant economic and environmental benefits to localities, states, and the nation. This article analyzes the micro- and macro scale impacts of implementing energy-efficient HVAC systems by integrating the industrial building energy data with the macroeconomic regional economic flow model. Micro-scale data include 10 years of historical energy, cost, and carbon dioxide savings achieved from energy-efficient HVAC implementation offered to manufacturers through industrial energy audits. The data were integrated into the macroeconomic modeling framework to illuminate the cascading regional economic impacts of implementing energy-efficient HVAC recommendations in manufacturing facilities. Results show that if recommendations had been implemented throughout all manufacturers in the region, $656 M energy costs would have been directly saved, 7.8 million metric tons of carbon dioxide emissions would have been avoided, and 4387 jobs could have been created, resulting in a total annual economic impact of $899 M stemming from direct, indirect, and induced impacts. The results offer insight into how industrial energy systems can be designed and provide models for how communities can accomplish a net-zero society.

Energy & Fuels↗
DOE OSTI2025

Multiscale Thermal-hydraulic analysis of the MARVEL micro-reactor using coupled MOOSE Subchannel (SCM) and SAM

MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory’s Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor’s primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL’s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Fracture‐Resistant and Thermally Insulating Ultrahigh‐Temperature Carbide Foams

Dense ultrahigh‐temperature ceramics (UHTCs) carbides are recognized as potential materials for thermal protection systems (TPS) owing to properties beyond existing structural materials’ capabilities. Recent advances in UHTCs have enabled the development of multiscale porous microstructures. Herein, it is highlighted that the porosity in UHTCs are no longer treated as a defect but as a functional property specifically tailored for thermal insulation. It is a promising solution to design and fabricate bulk UHTC foams via a freeze‐drying (FD) approach followed by calibrated pressureless spark plasma sintering. Herein, monolithic TaC and HfC UHTC foams and their composite show the partial solid–solution formation of (Ta, Hf)C with porosity ≥50%. TaC–HfC foam (≈80–92 N) shows an intermediate load‐bearing capability compared to monolithic TaC (≈120–135 N) and HfC (≈28–35 N) foams, with no evident cracking on the sample surface. The thermal conductivity of partial solid‐solution TaC–HfC foam increases up to fivefold compared to parent UHTC foams. In the results, solid solutions’ efficacy and pores’ unidirectionality in providing thermal insulation to TaC–HfC while maintaining its high‐load bearing capability are illustrated. In conclusion, the developed technique establishes a new paradigm shift in UHTCs, expanding their potential for TPS in extreme environments.

36 MATERIALS SCIENCE↗
DOE OSTI2024

SQMS Quantum R&D in Machine Learning, Optimization and Sensing beyond Fundamental Physics Applications

This newly formed team at SQMS under the Ecosystem Thrust is looking to develop capabilities impacting societal advances outside the core domain of HEP and condensed matter physics. We explicitly leverage the experimental and algorithmic innovations developed across all groups as well as connect to broad-scope external projects of the diverse team of PIs. As the inaugural set of projects, we are studying numerically quantum machine learning models inspired by efficiently trainable echo-state and orthogonal neural networks and developing designs for related experiments to be performed on quantum processors based on SQMS SRF cQED technology and Rigetti s transmon arrays. Investigated models exploit ideas and lessons learned from multiple prior work by SQMS team members in a variety of internal and external activities [R1]. Target initial applications include noisy signal processing, potentially captured by quantum sensors or noisy QPUs, as well as simulation and classification of healthcare data. For instance, image reconstruction of the brain s electrical properties by solving the inverse Maxwell equation problem with uncertainty [R2] through a hybrid quantum-classical physics-informed architecture for time-dependent processes [R3]. The group is also investigating the application and development of novel quantum sensors based on magnetic levitation of a superconducting sphere coupled to a superconducting qubit. This coupling enables high-precision measurements of the position of the sphere, which can be used for sensitive detection of forces, enabling practical applications such as gravimetry for geophysics analysis, or accelerometry for GPS-denied navigation [R4] [R1] Rieffel, Eleanor G., Ata Akbari Asanjan, M. Sohaib Alam, Namit Anand, David E. Bernal Neira, Sophie Block, Lucas T. Brady et al. "Assessing and advancing the potential of quantum computing: A NASA case study." Future Generation Computer Systems (2024). [R2] Yu, X., Serrall s, J.E., Giannakopoulos, I.I., Liu, Z., Daniel, L., Lattanzi, R. and Zhang, Z., 2023. Pifon-ept: Mr-based electrical property tomography using physics-informed fourier networks. IEEE Journal on Multiscale and Multiphysics Computational Techniques. [R3] Wudarski, Filip, Daniel OConnor, Shaun Geaney, Ata Akbari Asanjan, Max Wilson, Elena Strbac, P. Aaron Lott, and Davide Venturelli. "Hybrid quantum-classical reservoir computing for simulating chaotic systems." arXiv preprint arXiv:2311.14105 (2023). [R4] Higgins, Gerard, Saarik Kalia, and Zhen Liu. "Maglev for dark matter: Dark-photon and axion dark matter sensing with levitated superconductors." Physical Review D 109.5 (2024): 055024.

Venturelli, Davide↗
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2022

A distributed temperature profiling system for vertically and laterally dense acquisition of soil and snow temperature

Abstract. Measuring soil and snow temperature with high vertical and lateral resolution is critical for advancing the predictive understanding of thermal and hydro-biogeochemical processes that govern the behavior of environmental systems. Vertically resolved soil temperature measurements enable the estimation of soil thermal regimes, frozen-/thawed-layer thickness, thermal parameters, and heat and/or water fluxes. Similarly, they can be used to capture the snow depth and the snowpack thermal parameters and fluxes. However, these measurements are challenging to acquire using conventional approaches due to their total cost, their limited vertical resolution, and their large installation footprint. This study presents the development and validation of a novel distributed temperature profiling (DTP) system that addresses these challenges. The system leverages digital temperature sensors to provide unprecedented, finely resolved depth profiles of temperature measurements with flexibility in system geometry and vertical resolution. The integrated miniaturized logger enables automated data acquisition, management, and wireless transfer. A novel calibration approach adapted to the DTP system confirms the factory-assured sensor accuracy of ±0.1 ∘C and enables improving it to ±0.015 ∘C. Numerical experiments indicate that, under normal environmental conditions, an additional error of 0.01 % in amplitude and 70 s time delay in amplitude for a diurnal period can be expected, owing to the DTP housing. We demonstrate the DTP systems capability at two field sites, one focused on understanding how snow dynamics influence mountainous water resources and the other focused on understanding how soil properties influence carbon cycling. Results indicate that the DTP system reliably captures the dynamics in snow depth and soil freezing and thawing depth, enabling advances in understanding the intensity and timing in surface processes and their impact on subsurface thermohydrological regimes. Overall, the DTP system fulfills the needs for data accuracy, minimal power consumption, and low total cost, enabling advances in the multiscale understanding of various cryospheric and hydro-biogeochemical processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Generative learning for slow manifolds and bifurcation diagrams

In dynamical systems characterized by separation of time scales, the approximation of so called “slow manifolds”, on which the long term dynamics lie, is a useful step for model reduction. Initializing on such slow manifolds is a useful step in modeling, since it circumvents fast transients, and is crucial in multiscale algorithms (like the equation-free approach) alternating between fine scale (fast) and coarser scale (slow) simulations. In a similar spirit, when one studies the infinite time dynamics of systems depending on parameters, the system attractors (e.g., its steady states) lie on bifurcation diagrams (curves for one-parameter continuation, and more generally, on manifolds in state parameter space. Sampling these manifolds gives us representative attractors (here, steady states of ODEs or PDEs) at different parameter values. Algorithms for the systematic construction of these manifolds (slow manifolds, bifurcation diagrams) are required parts of the “traditional” numerical nonlinear dynamics toolkit. In more recent years, as the field of Machine Learning develops, conditional score-based generative models (cSGMs) have been demonstrated to exhibit remarkable capabilities in generating plausible data from target distributions that are conditioned on some given label. It is tempting to exploit such generative models to produce samples of data distributions (points on a slow manifold, steady states on a bifurcation surface) conditioned on (consistent with) some quantity of interest (QoI, observable). In this work, we present a framework for using cSGMs to quickly (a) initialize on a low-dimensional (reduced-order) slow manifold of a multi-time-scale system consistent with desired value(s) of a QoI (a “label”) on the manifold, and (b) approximate steady states in a bifurcation diagram consistent with a (new, out-of-sample) parameter value. This conditional sampling can help uncover the geometry of the reduced slow-manifold and/or approximately “fill in” missing segments of steady states in a bifurcation diagram. Finally, the quantity of interest, which determines how the sampling is conditioned, is either known a priori or identified using manifold learning-based dimensionality reduction techniques applied to the training data.

Dynamical systems↗
DOE OSTI2025

Effective Li-Ion Transport Quantification in Composite Cathodes for All-Solid-State Batteries via Multiscale Modeling and Experiments

The tortuosity factor of composite cathodes significantly affects the rate performance of all-solid-state batteries (ASSBs) and has significant differences from systems with liquid electrolytes. Here, in this work, we report a simulation-experiment combined approach that quantifies the effective Li-ion transport in an ASSB composite cathode, which links tortuosity factor on ∼ μm scale to terminal voltage during cycling at the cell level (on ∼ cm scale). Two independent approaches of tortuosity factor quantification are considered: fitting electrochemical cycling data and verifying at different cycling rates and calculating from segmented tomography images, with the tortuosity factor quantified from both methods reaching self-consistency. The simulated terminal voltage using the quantified tortuosity factor has a small relative error of <3% compared to the experimental measurements. We find a significantly reduced value of the Bruggeman exponent of the catholyte phase (1.75), and using shape analysis, we show that rod-shaped catholyte particles play an important role in lowering the tortuosity factor.

Yao, Archie Mingze [Univ. of Michigan, Ann Arbor, ↗
DOE OSTI2021

Multiscale Land-Atmosphere Interaction in Tropical Ecosystems

A leading source of uncertainty in the magnitude of future climate change arises from feedback between the physical climate system and the carbon cycle. Output from 15 Earth System Models (ESMs) in the Coupled Model Intercomparison Project (CMIP5) shows that when forced by identical emissions, predicted atmospheric CO 2 concentrations in 2100 ranged from 780 ppm to 1120 ppm due to variations in carbon cycle feedback among the models. The resulting difference in radiative forcing of climate is almost 2 W m -2 , comparable to the total uncertainty in the physical climate system for the representative concentration pathway (RCP) 8.5 scenario in 2100. Scaling future carbon sources and sinks to match contemporary atmospheric CO 2 reduced this uncertainty by a factor of five, suggesting that contemporary observations might dramatically constrain future predictions, but their scaling had no mechanistic basis, so it provides little insight on carbon cycle processes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Absence of Barren Plateaus and Scaling of Gradients in the Energy Optimization of Isometric Tensor Network States

Abstract Vanishing gradients can pose substantial obstacles for high-dimensional optimization problems. Here we consider energy minimization problems for quantum many-body systems with extensive Hamiltonians and finite-range interactions, which can be studied on classical computers or in the form of variational quantum eigensolvers on quantum computers. Barren plateaus correspond to scenarios where the average amplitude of the energy gradient decreases exponentially with increasing system size. This occurs, for example, for quantum neural networks and for brickwall quantum circuits when the depth increases polynomially in the system size. Here we prove that the variational optimization problems for matrix product states, tree tensor networks, and the multiscale entanglement renormalization ansatz are free of barren plateaus. The derived scaling properties for the gradient variance provide an analytical guarantee for the trainability of randomly initialized tensor network states (TNS) and motivate certain initialization schemes. In a suitable representation, unitary tensors that parametrize the TNS are sampled according to the uniform Haar measure. We employ a Riemannian formulation of the gradient based optimizations which simplifies the analytical evaluation.

Barthel, Thomas↗
DOE OSTI2023

Machine learning from RANS and LES to inform coarse grid simulations

Nuclear system thermal hydraulic analysis has historically relied on computationally inexpensive 1D codes. However, such tools are unable to capture multiscale multidimensional effects in large nuclear reactor enclosures. On the other hand, simulations with higher fidelity can be too expensive for such purposes. One of the ways to reduce computational cost is to perform simulations on a coarse grid, which, unfortunately, introduces large discretization errors. In this paper, two high-to-low data-driven approaches are investigated: (1) a coarse grid turbulence model to predict eddy viscosity and (2) correction of errors in coarse grid velocity fields. The approaches aim to reduce grid- and turbulence model-induced errors in coarse grid Reynolds-averaged Navier–Stokes (RANS) simulations. Two sources of high-fidelity data, RANS and large eddy simulations (LES), are explored. To extract the eddy viscosity from the LES data, an inverse optimization problem is solved. However, the LES eddy viscosity is shown to be comparable to the RANS eddy viscosity in terms of error reduction. Therefore, the directly available RANS eddy viscosity was used to develop a coarse grid data-driven turbulence model. Additionally, error correction in velocity is used to reduce the remaining uncertainties and bring the results closer to reality. In conclusion, the performance of the frameworks is demonstrated for a scaled upper plenum of a gas-cooled reactor facility.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair P D1 –P D2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state P D1 – –P D2 + stabilizes energetically below the excitonic state. Importantly, this underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

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

A Multiscale Simulation Study of Influence of Morphology on Ion Transport in Block Copolymeric Ionic Liquids

In this work, we present results of a multiscale simulation study involving coarse-graining and reverse-mapping steps probing at an atomistic resolution the influence of morphology on ion transport in block copolymer polymeric ionic liquids. We provide a detailed description of the multiscale simulation methodology and then describe the results in the context of four microphase-separated morphologies: two lamella systems, a cylinder, and a gyroid morphology. For the framework adopted in this study (in which the total number of ions was maintained fixed), morphology had little influence on the fraction of ions segregating to the interface and the coordination and hopping characteristics of such interfacially present ions. Such results manifested in anion mobilities being insensitive to the morphology once the dimensionality of the morphology was accounted for. Overall, our results are consistent with the hypothesis that the mobility of ions in such microphase-separated morphologies, after accounting for the dimensionality effects, can be roughly correlated to a linear superposition based on the fraction of bulk anions, which possess mobilities unperturbed from their values in the conducting homopolymers, and the interfacial anions, which exhibit much lower mobilities.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Molecular insights into CO 2 -to-bicarbonate transformation in functionalized anion exchange ionomers for electrochemical separations

Bipolar membrane (BPM) electrochemical processes are a promising platform for carbon dioxide (CO 2 ) separations, but the molecular level thermodynamic and kinetic understanding of CO 2 -to-bicarbonate (HCO 3 − ) transformation remain poorly understood. This study employs a multiscale computational approach to systematically explore the adsorption and reactive transformation of CO 2 in five anion exchange ionomer systems. Classical molecular dynamics (MD) simulation results demonstrate that polymers with imidazolium groups significantly reduce CO 2 diffusion and enhance (OH − )–CO 2 interactions due to stronger electrostatic and π-interactions. Compared to the commonly used quaternary ammonium ionomers, imidazolium-functionalized ionomers show improved CO 2 proximity and interaction strength. Ab initio MD and density functional theory (DFT) calculations reveal that the benzyl-substituted imidazolium (IM-Ben) substantially reduces the energy barrier for HCO 3 − formation (∼72 meV lower) compared to the alkyl-substituted IM-nBu, while also mitigating imidazolium deprotonation under moderate hydration conditions. Transition state analysis shows IM-Ben forms more extensive hydrogen-bonding networks, which stabilize the transition state structure and contribute to a lower energy barrier for bicarbonate formation. These findings highlight the advantage of the adjacent benzyl moiety in enabling efficient CO 2 -to-bicarbonate transformation via hydrated hydroxide ion counterions, offering mechanistic insights and clear molecular design principles for optimizing anion exchange ionomers at bipolar membrane interfaces for electrochemical CO 2 separation applications.

Bipolar membranes, Reactive transformation of CO2,↗