Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PARTITION”

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 973 records · Page 54

Molten-salt-mediated carbon dioxide capture and superequilibrium utilization with ethane oxidative dehydrogenation

Existing CO 2 -mediated oxidative dehydrogenation (CO 2 -ODH) of ethane has yet to demonstrate >60% single-pass CO yield due to the intrinsic equilibrium limitations. We report a unique approach with mixed molten carbonates as a reaction medium for CO 2 -ODH, which strategically partitions the CO 2 -ODH reactions into gas and molten-salt phases and facilitates integrated CO 2 capture from power plant flue gases. An 89% CO yield was achieved at 770°C, doubling the equilibrium limitation of conventional CO 2 -ODH. The high CO yield in turn enhances ethylene formation. Further characterizations confirmed that molten-salt mediated ODH (MM-ODH) proceeds through a gas-phase cracking and molten-salt mediated reverse water-gas-shift reaction pathway. Based on this understanding, thermodynamic analysis and ab initio molecular dynamics simulations were conducted to develop general principles to optimize the molten-salt reaction medium. Process analyses confirm that MM-ODH has the potential to be significantly more efficient for CO 2 capture and utilization than conventional CO 2 -ODH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Proteome reallocation enables the selective de novo biosynthesis of non-linear, branched-chain acetate esters

The one-carbon recursive ketoacid elongation pathway is responsible for making various branched-chain amino acids, aldehydes, alcohols, and acetate esters in living cells. Controlling selective microbial biosynthesis of these target molecules at high efficiency is challenging due to enzyme promiscuity, regulation, and metabolic burden. In this study, we present a systematic modular design approach to control proteome reallocation for selective microbial biosynthesis of branched-chain acetate esters. Through pathway modularization, we partitioned the branched-chain ester pathways into four submodules including keto-isovalerate submodule for converting pyruvate to keto-isovalerate, ketoacid elongation submodule for producing longer carbon-chain keto-acids, ketoacid decarboxylase submodule for converting ketoacids to alcohols, and alcohol acyltransferase submodule for producing branched-chain acetate esters by condensing alcohols and acetyl-CoA. By systematic manipulation of pathway gene replication and transcription, enzyme specificity of the first committed steps of these submodules, and downstream competing pathways, we demonstrated selective microbial production of isoamyl acetate over isobutyl acetate. In this work, we found that the optimized isoamyl acetate pathway globally redistributed the amino acid fractions in the proteomes and required up to 23-31% proteome reallocation at the expense of other cellular resources, such as those required to generate precursor metabolites and energy for growth and amino acid biosynthesis. From glucose fed-batch fermentation, the engineered strains produced isoamyl acetate up to a titer of 8.8 g/L (> 0.25 g/L toxicity limit), a yield of 0.22 g/g (61% of maximal theoretical value), and 86% selectivity, achieving the highest titers, yields and selectivity of isoamyl acetate reported to date.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Controlling circuitry underlies the growth optimization of Saccharomyces cerevisiae

Microbial growth emerges from coordinated synthesis of various cellular components from limited resources. In Saccharomyces cerevisiae, cyclic AMP (cAMP)-mediated signaling is shown to orchestrate cellular metabolism; however, it remains unclear quantitatively how the controlling circuit drives resource partition and subsequently shapes biomass growth. Here we combined experiment with mathematical modeling to dissect the signaling-mediated growth optimization of S. cerevisiae. We showed that, through cAMP-mediated control, the organism achieves maximal or nearly maximal steady-state growth during the utilization of multiple tested substrates as well as under perturbations impairing glucose uptake. However, the optimal cAMP concentration varies across cases, suggesting that different modes of resource allocation are adopted for varied conditions. Under settings with nutrient alterations, S. cerevisiae tunes its cAMP level to dynamically reprogram itself to realize rapid adaptation. Moreover, to achieve growth maximization, cells employ additional regulatory systems such as the GCN2-mediated amino acid control. This study establishes a systematic understanding of global resource allocation in S. cerevisiae, providing insights into quantitative yeast physiology as well as metabolic strain engineering for biotechnological applications.

59 BASIC BIOLOGICAL SCIENCES↗

Metabolic flux optimization of iterative pathways through orthogonal gene expression control: Application to the β-oxidation reversal

Balancing relative expression of pathway genes to minimize flux bottlenecks and metabolic burden is one of the key challenges in metabolic engineering. This is especially relevant for iterative pathways, such as reverse β-oxidation (rBOX) pathway, which require control of flux partition at multiple nodes to achieve efficient synthesis of target products. Here, we develop a plasmid-based inducible system for orthogonal control of gene expression (referred to as the TriO system) and demonstrate its utility in the rBOX pathway. Leveraging effortless construction of TriO vectors in a plug-and-play manner, we simultaneously explored the solution space for enzyme choice and relative expression levels. Remarkably, varying individual expression levels led to substantial change in product specificity ranging from no production to optimal performance of about 90% of the theoretical yield of the desired products. We obtained titers of 6.3 g/L butyrate, 2.2 g/L butanol and 4.0 g/L hexanoate from glycerol in E. coli, which exceed the best titers previously reported using equivalent enzyme combinations. Since a similar system behavior was observed with alternative termination routes and higher-order iterations, we envision our approach to be broadly applicable to other iterative pathways besides the rBOX. Here, considering that high throughput, automated strain construction using combinatorial promoter and RBS libraries remain out of reach for many researchers, especially in academia, tools like the TriO system could democratize the testing and evaluation of pathway designs by reducing cost, time and infrastructure requirements.

59 BASIC BIOLOGICAL SCIENCES↗

Phase transition to blob-hole coherent structure in the Hasegawa–Mima model for plasmas

An equilibrium statistical mechanics theory for the Hasegawa–Mima equations of toroidal plasmas, with canonical constraint on energy and microcanonical constraint on potential enstrophy, is solved exactly as a spherical model. The use of a canonical energy constraint instead of a fixed-energy microcanonical approach is justified by the preference for viewing real plasmas as an open system. A significant consequence of the results obtained from the partition function, free energy and critical temperature, is the condensation into a ground state exhibiting a blob-hole-like structure observed in real plasmas.

Physics↗

Modelling the nonlinear plasma response to externally applied three-dimensional fields with the Stepped Pressure Equilibrium Code

Small-amplitude, symmetry-breaking magnetic field perturbations, including resonant magnetic perturbations (RMPs) and error fields, can profoundly impact plasma properties in both tokamaks and stellarators. In this work, we perform the first comparison between the Stepped Pressure Equilibrium Code (SPEC) (a comparatively fast and efficient equilibrium code based on energy-minimisation principles) and M3D-C 1 (a high-fidelity albeit computationally expensive initial-value extended-magnetohydro-dynamic (MHD) code) to assess the conditions under which SPEC can be used to model the nonlinear, non-ideal plasma response to an externally applied (m = 2,n = 1) RMP field in an experimentally relevant geometry. We find that SPEC is able to capture the plasma response in the weakly nonlinear regime – meaning perturbation amplitudes below the threshold for break up of the separatrix and onset of secondary magnetic island formation – when around half of the total toroidal flux is enclosed in the volume containing the q = 2 resonant surface. Further, the observed dependence of SPEC solutions on input parameters, including toroidal flux and the number of volumes into which the plasma is partitioned, indicates that additional exploration of the underlying Multi-Region Relaxed MHD physics model is needed to constrain the choice of parameters. Nonetheless, this work suggests promising applications of SPEC to optimisation and fusion plasma design.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temperature separation under compression of moderately coupled plasma

In moderately coupled plasmas, a significant fraction of the internal energy resides in electric fields. As these plasmas are heated or compressed, the shifting partition of energy between particles and fields leads to surprising effects, particularly when ions and electrons have different temperatures. In this work, quasi-equations of state (quasi-EOS) are derived for two-temperature moderately coupled plasma in a thermodynamic framework and expressed in a simple form. These quasi-EOS readily yield expressions for correlation heating, in which heating of the electrons causes a rapid increase in ion temperature even in the absence of collisional energy exchange between species. It is also shown that, remarkably, compression of moderately coupled plasma drives a temperature difference between electrons and ions, even when the species start at equal temperatures. These additional channels for ion heating may be relevant in designing ignition schemes for inertial confinement fusion.

Fetsch, H. (ORCID:0000000256133975)↗

Simulations of radiatively cooled magnetic reconnection driven by pulsed power

Magnetic reconnection is an important process in astrophysical environments, as it reconfigures magnetic field topology and converts magnetic energy into thermal and kinetic energy. In extreme astrophysical systems, such as black hole coronae and pulsar magnetospheres, radiative cooling modifies the energy partition by radiating away internal energy, which can lead to the radiative collapse of the reconnection layer. In this paper, we perform two- and three-dimensional simulations to model the MARZ (Magnetic Reconnection on Z) experiments, which are designed to access cooling rates in the laboratory necessary to investigate reconnection in a previously unexplored radiatively cooled regime. These simulations are performed in GORGON, an Eulerian two-temperature resistive magnetohydrodynamic code, which models the experimental geometry comprising two exploding wire arrays driven by 20 MA of current on the Z machine (Sandia National Laboratories). Radiative losses are implemented using non-local thermodynamic equilibrium tables computed using the atomic code Spk, and we probe the effects of radiation transport by implementing both a local radiation loss model and P 1/3 multi-group radiation transport. The load produces highly collisional, super-Alfvénic (Alfvén Mach number M A ≈ 1.5), supersonic (Sonic Mach number M S ≈ 4–5) strongly driven plasma flows which generate an elongated reconnection layer (Aspect Ratio L/δ≈100 , Lundquist number S L ≈ 400). The reconnection layer undergoes radiative collapse when the radiative losses exceed the rates of ohmic and compressional heating (cooling rate/hydrodynamic transit rate = τ$^{–1}_{cool}$/τ$^{–1}_{H}$ ≈100); this generates a cold strongly compressed current sheet, leading to an accelerated reconnection rate, consistent with theoretical predictions. Lastly, the current sheet is also unstable to the plasmoid instability, but the magnetic islands are extinguished by strong radiative cooling before ejection from the layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Aerosol Optical Tweezers Elucidate the Chemistry, Acidity, Phase Separations, and Morphology of Atmospheric Microdroplets

Aerosol particles represent unique chemical environments because of their high surface area-to-volume ratio that promotes the effects of interfacial chemistry in confined environments. Properties such as viscosity, diffusivity, water content, pH, and morphology—following liquid–liquid phase separation—can strongly alter how a particle interacts with condensable vapors and reactive trace gases, thus modifying its continual evolution and environmental effects. Our understanding of this chemical evolution of atmospheric particulate matter and its environmental impacts is largely limited by our ability to directly observe how these critical particle properties respond to the addition or reactive uptake of new chemical components. Aerosol optical tweezers (AOT) stably trap particles in focused laser beams, providing positional control and the retrieval of many of these critical properties required to understand and predict the chemistry of aerosolized microdroplets. The analytical power of the AOT stems from the retrieval of the cavity-enhanced Raman spectrum induced by the trapping laser. Analysis of the whispering gallery modes (WGMs) that resonate as a standing wave around the droplet’s interface, provide high accuracy measurements of the droplet’s size, refractive index (and thus a measurement of composition), and can distinguish between core–shell, partially engulfed, and homogeneous morphologies. In this study, we have advanced the ability to determine the properties of the core and shell phases in biphasic droplets, including obtaining high-accuracy pH measurements. These capabilities were applied to perform AOT physical chemistry experiments on authentic secondary organic aerosol (SOA) produced directly in the AOT chamber by ozonolysis of terpene vapors. The propensity of the SOA to phase separate as a shell from a wide range of nonpolar to polar core phases was observed, along with the discovery of a stable emulsified state of SOA particles in an aqueous salt droplet. Micron-thick SOA shells did not impede the gain or loss of water or squalane from the core to the surrounding air, indicating no significant diffusional limitations to condensational growth or partitioning even under dry conditions. These experiments formed the foundation of a new framework that predicts how the phase-separated morphology of complex aerosols containing organic carbon evolves during continual atmospheric oxidation processes. Increases in oxidation state will quickly drive conversion from a partially engulfed to core–shell morphology that has dramatically different chemical reactivity since the core phase is completely concealed by the shell. The recent advances in the experimental capabilities of the AOT technique such as presented here enable novel experimental methodologies that provide insights into the chemistry and multidimensional properties of aerosol microdroplets, and how these coevolve and respond to continual chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab Initio Direct Dynamics

The reactivity and dynamics of molecular systems can be explored computationally by classical trajectory calculations. The traditional approach involves fitting a functional form of a potential energy surface (PES) to the energies from a large number of electronic structure calculations and then integrating numerous trajectories on this fitted PES to model the molecular dynamics. The ever-decreasing cost of computing and continuing advances in computational chemistry software have made it possible to use electronic structure calculations directly in molecular dynamics simulations without first having to construct a fitted PES. In this “on-the-fly” approach, every time the energy and its derivatives are needed for the integration of the equations of motion, they are obtained directly from quantum chemical calculations. This approach started to become practical in the mid-1990s as a result of increased availability of inexpensive computer resources and improved computational chemistry software. The application of direct dynamics calculations has grown rapidly over the last 25 years and would require a lengthy review article. The present Account is limited to some of our contributions to methods development and various applications. To improve the efficiency of direct dynamics calculations, we developed a Hessian-based predictor-corrector algorithm for integrating classical trajectories. Hessian updating made this even more efficient. Furthermore, this approach was also used to improve algorithms for following the steepest descent reaction paths. For larger molecular systems, we developed an extended Lagrangian approach in which the electronic structure is propagated along with the molecular structure. Strong field chemistry is a rapidly growing area, and to improve the accuracy of molecular dynamics in intense laser fields, we included the time-varying electric field in a novel predictor-corrector trajectory integration algorithm. Since intense laser fields can excite and ionize molecules, we extended our studies to include electron dynamics. Specifically, we developed code for time-dependent configuration interaction electron dynamics to simulate strong field ionization by intense laser pulses. Our initial application of ab initio direct dynamics in 1994 was to CH 2 O → H 2 + CO; the calculated vibrational distributions in the products were in very good agreement with experiment. In the intervening years, we have used direct dynamics to explore energy partitioning in various dissociation reactions, unimolecular dissociations yielding three fragments, reactions with branching after the transition state, nonstatistical dynamics of chemically activated molecules, dynamics of molecular fragmentation by intense infrared laser pulses, selective activation of specific dissociation channels by aligned intense infrared laser fields, angular dependence of strong field ionization, and simulation of sequential double ionization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Population Balance Models for Catalytic Depolymerization: From Elementary Steps to Multiphase Reactors

Here, the ongoing accumulation of plastic waste in landfills and in the environment is driving research on chemical processes and catalysts to recycle polymers. Traditional modeling strategies are not applicable to these processes because they involve too many reactants and intermediates, one for each molecular weight and each functionalization. To model the kinetics, we have developed population balance models (PBMs) that account for macromolecular reactants in the bulk and macromolecular catalytic intermediates. These PBMs couple to each other through polymer adsorption and desorption models and to traditional rate equations for small molecule products and co-reactants (like hydrogen or ethylene). The models, in combination with experimental data, are being used in many ways: (i) to test mechanistic hypotheses, (ii) to extract rate parameters, (iii) to quantitatively compare catalyst activities, (iv) to account for mass transfer and vapor–liquid partitioning in two-phase reactors, and (v) to design novel support architectures and catalysts that mimic the processive action of natural depolymerization enzymes. Some key theoretical advances allow PBMs to be constructed from elementary rates and mechanisms, as opposed to traditional formulations with pseudoelementary rate parameters invoked as fitting parameters. We discuss ways to build these models “bottom-up” from first-principles calculations and ways to extract model parameters from “top down” analyses of rate data. The combination provides a quantitative bridge between first-principles calculations and the kinetics of complex macromolecular transformations for polymer upcycling and beyond.

Manis, Lela K. [University of Illinois at Urbana-C↗

Probing Nanoscale Chemical Environments of Zinc in Diatoms

Diatoms impact the biogeochemical cycling of zinc (Zn) due to elevated cellular Zn levels over other phytoplankton, and their strong biological uptake of Zn in the Southern Ocean sets the global distribution of Zn in ocean water columns. Past studies have revealed the abundance and spatial distribution of Zn in individual diatoms, while others have shown diverse Zn chemical species in whole-cell aggregates. However, intracellular and intercellular variations in Zn chemical forms remained unknown. For the first time, we applied a synchrotron X-ray nanoprobe to cultured diatoms (Phaeodactylum tricornutum and Chaetoceros muelleri) and resolved spatial heterogeneities of Zn chemistry. The result shows the partitioning of Zn between intracellular contents (soft parts) and frustules (mainly in the intracellular contents for P. tricornutum and similar Zn concentrations in the intracellular contents and frustules for C. muelleri). We found multiple Zn chemical species unevenly distributed within individual cells, where Zn–phosphoryl complexes were the most abundant followed by cysteine, histidine, biogenic silica, and carboxyl complexes. Cellular Zn species varied between individual organisms and were influenced by Zn availability. To our knowledge, this work presents the first X-ray measurements of trace metal speciation at sub-100 nm resolution in biological samples in their natural state. In conclusion, this research improves our ability to examine the biogeochemistry of Zn on the nanoscale and can help us understand the role of Zn in plankton growth and the Zn cycle in ocean waters.

36 MATERIALS SCIENCE↗

Thermal Gradient Effects on Redox Evolution and Volatility-Driven Fractionation in Ternary U/Ce/Cs Condensates

Understanding how thermal history influences redox evolution and chemical fractionation is essential for characterizing high-temperature condensation in complex materials, including nuclear debris. Here, we tested the hypothesis that distinct thermal regimes in a plasma flow reactor influence redox pathways and elemental partitioning in ternary U/Ce/Cs systems. A configurable plasma flow reactor was modified with an external tube furnace to impose two distinct thermal gradients: continuous ambient cooling and a furnace-assisted thermal hold-up near 1400 K followed by rapid cooling. Transmission electron microscopy characterized phase identity, morphology, and nanoscale element distributions, while inductively coupled plasma-mass spectrometry quantified bulk elemental ratios. Across both thermal regimes, uranium and cerium condensed as UO 2 and CeO 2 as dominant refractory oxide products. Uranium partially oxidized to α-UO 3 during extended ambient cooling, while furnace-assisted hold-up preserved UO 2 and produced partial reduction of cerium to Ce 2 O 3 . Cesium remained volatile upstream and condensed later in the reactor, forming Cs 2 O and Cs-uranate phases with the highest incorporation after thermal hold-up. Bulk ICP-MS measurements supported these observations. U/Ce ratios remained comparatively stable and Cs displayed delayed and apparent transient enrichment that matched the nanoscale measurements. This integrated approach provides a quantitative method for linking thermal gradients to redox evolution and volatility-driven fractionation. These results show how the plasma flow reactor can identify where equilibrium descriptions remain adequate and where kinetic effects from residence time and temperature history must be considered when interpreting condensation behavior in multicomponent systems.

and nuclear chemistry↗

Incorporation and Assembly of a Light-Emitting Enzymatic Reaction into Model Protein Condensates

Eukaryotic cells partition enzymes and other cellular components into distinct subcellular compartments to generate specialized biochemical niches. A subclass of these compartments form in the absence of lipid membranes, via liquid-liquid phase separation of proteins to form biomolecular condensates or “membraneless organelles” such as nucleoli, stress granules, and P-bodies. Because of their ability to form compartments from simple starting materials, membraneless organelles are an attractive target for engineering new functionalities in both living cells and protocells. In this work, we demonstrate incorporation of novel enzymatic activity in protein coacervates with a light-generating enzyme, NanoLuc, to produce bioluminescence. Using condensates comprised of the disordered RGG domain of C. elegans LAF-1, we show functionalization of condensates with enzymatic activity in vitro and that localization to protein coacervates enhances the assembly and activity of split enzymes. To build condensates that function as light emitting reactors, we designed a NanoLuc enzyme flanked by RGG domains. Furthermore, the resulting condensates concentrated NanoLuc by 10-fold over bulk solution and display significantly increased net reaction rates. We further show that condensate viscosity impacts light emission due to diffusion-limited behavior. By splitting NanoLuc enzyme into its constituent components, we demonstrate that NanoLuc activity can be reconstituted via co-condensation. Further, we demonstrate control of the spatial localization of enzyme within condensates by targettng NanoLuc to the surface of in vitro condensates. Collectively, this work demonstrates that membraneless organelles can be endowed with localized enzymatic activity, and that this activity can be spatially and temporally controlled via enzyme reconstitution and design of protein surfactants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of Salt on Phase Behavior and Rheological Properties of Alginate–Chitosan Polyelectrolyte Complexes

Oppositely charged polyelectrolytes often form polyelectrolyte complexes (PECs) due to the association through electrostatic interactions. Obtaining PECs using natural, biocompatible polyelectrolytes is of interest in the food, pharmaceutical, and biomedical industries. In this work, PECs were prepared from two biopolymers, positively charged chitosan and negatively charged alginate. We investigate the changes in the structure and properties of PECs by adding sodium chloride (salt doping) to the system. The shear modulus of PECs can be tuned from ~10 to 10 4 Pa by changing the salt concentration. The addition of salt led to a decrease in the water content of the complex phase with increasing shear modulus. However, at a very high salt concentration, the shear modulus of the complex phase decreased but did not lead to the liquid coacervate formation, typical of synthetic polyelectrolytes. This difference in phase behavior has likely been attributed to the hydrophobicity of chitosan and long semiflexible alginate and chitosan chains that restrict the conformational changes. Large amplitude oscillatory shear experiments captured nonlinear responses of PECs. The compositions of the PECs, determined as a function of salt concentration, signify the preferential partitioning of salt into the complex phase. Small-angle X-ray scattering of the salt-doped PECs indicates that the Kuhn length and radius of the alginate–chitosan associated structure qualitatively agree with the captured phase behavior and rheological data. This study provides insights into the structure–property as a function of salt concentration of natural polymer-based PECs necessary for developing functional materials from natural polyelectrolytes.

59 BASIC BIOLOGICAL SCIENCES↗

Multireference Embedding and Fragmentation Methods for Classical and Quantum Computers: From Model Systems to Realistic Applications

One of the primary challenges in quantum chemistry is the accurate modeling of strong electron correlation. While multireference methods effectively capture such correlation, their steep scaling with system size prohibits their application to large molecules and extended materials. Quantum embedding offers a promising solution by partitioning complex systems into manageable subsystems. In this Review, we highlight recent advances in multireference density matrix embedding and localized active space self-consistent field approaches for complex molecules and extended materials. We discuss both classical implementations and the emerging potential of these methods on quantum computers. Here, by extending classical embedding concepts to the quantum landscape, these algorithms have the potential to expand the reach of multireference methods in quantum chemistry and materials.

Algorithms↗

Impact of a Novel Nickel-Based Catalyst and Phenyl-Acrylate-Based Anion-Exchange Membrane in a Direct Urea Fuel Cell

Developing target-specific catalysts and anion-exchange membranes (AEMs) is crucial for direct urea fuel cell (DUFC) performance. To advance the DUFC system, we developed an anode catalyst with a nickel–iron oxyhydroxide/carbon (NiFeOOH/C) nanofibrous structure for the urea oxidation reaction (UOR), where we optimized the Ni/Fe molar ratio as 6:4. The enhanced electrocatalytic activity of the anode (Ni 6 Fe 4 OOH/C) is attributed to the hydroxide group, which responds with urea molecules to enhance the UOR in a pH-neutral system. Here, we employed a recently developed cross-linked phenyl-acrylate-based AEM (PA/M). A DUFC prepared with the anode and PA/M generates a maximum power density of 11.1 mW/cm 2 and 0.92 V open-circuit voltage under 3 M urea as fuel at 25 °C. We further analyzed the applicability of PA/M in a DUFC system by measuring the urea partition coefficients and permeabilities over a range of concentrations.

10 SYNTHETIC FUELS↗

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗