Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interface coupling”

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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular model of TFIIH recruitment to the transcription-coupled repair machinery

Transcription-coupled repair (TCR) is a vital nucleotide excision repair sub-pathway that removes DNA lesions from actively transcribed DNA strands. Binding of CSB to lesion-stalled RNA Polymerase II (Pol II) initiates TCR by triggering the recruitment of downstream repair factors. Yet it remains unknown how transcription factor IIH (TFIIH) is recruited to the intact TCR complex. Combining existing structural data with AlphaFold predictions, we build an integrative model of the initial TFIIH-bound TCR complex. We show how TFIIH can be first recruited in an open repair-inhibited conformation, which requires subsequent CAK module removal and conformational closure to process damaged DNA. In our model, CSB, CSA, UVSSA, elongation factor 1 (ELOF1), and specific Pol II and UVSSA-bound ubiquitin moieties come together to provide interaction interfaces needed for TFIIH recruitment. STK19 acts as a linchpin of the assembly, orienting the incoming TFIIH and bridging Pol II to core TCR factors and DNA. Molecular simulations of the TCR-associated CRL4CSA ubiquitin ligase complex unveil the interplay of segmental DDB1 flexibility, continuous Cullin4A flexibility, and the key role of ELOF1 for Pol II ubiquitination that enables TCR. Collectively, these findings elucidate the coordinated assembly of repair proteins in early TCR.

Paul, Tanmoy↗
DOE OSTI2024

Controlled patterning of crystalline domains by frontal polymerization

Materials with hierarchical architectures that combine soft and hard material domains with coalesced interfaces possess superior properties compared with their homogeneous counterparts. These architectures in synthetic materials have been achieved through deterministic manufacturing strategies such as 3D printing, which require an a priori design and active intervention throughout the process to achieve architectures spanning multiple length scales. Here we harness frontal polymerization spin mode dynamics to autonomously fabricate patterned crystalline domains in poly(cyclooctadiene) with multiscale organization. This rapid, dissipative processing method leads to the formation of amorphous and semi-crystalline domains emerging from the internal interfaces generated between the solid polymer and the propagating cure front. The size, spacing and arrangement of the domains are controlled by the interplay between the reaction kinetics, thermochemistry and boundary conditions. Small perturbations in the fabrication conditions reproducibly lead to remarkable changes in the patterned microstructure and the resulting strength, elastic modulus and toughness of the polymer. Furthermore, this ability to control mechanical properties and performance solely through the initial conditions and the mode of front propagation represents a marked advancement in the design and manufacturing of advanced multiscale materials. Drawing inspiration from biological systems in which structural complexity develops through dissipative reaction–diffusion processes, this study explores a transformative synthetic manufacturing strategy aimed at harnessing the principles underpinning morphogenic growth, unlocking new avenues for advanced materials design and fabrication. Synthetic coupled reaction-transport processes offer a versatile yet relatively underexplored method to manipulate the spatial attributes of synthetic materials10. Here we introduce an innovative manufacturing approach based on frontal ring-opening metathesis polymerization (FROMP) that draws parallels with morphogenic growth and development, enabling the formation of patterned microstructures within polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing Electrolyte Influence on CO 2 Reduction in Aprotic Solvents

Selective CO 2 capture and electrochemical conversion are important tools in the fight against climate change. Industrially, CO 2 is captured using a variety of aprotic solvents due to their high CO 2 solubility. However, most research efforts on electrochemical CO 2 conversion use aqueous media and are plagued by competing hydrogen evolution reaction (HER) from water breakdown. Fortunately, aprotic solvents can circumvent HER, making it important to develop strategies that enable integrated CO 2 capture and conversion. However, the influence of ion solvation and solvent selection within nonaqueous electrolytes for efficient and selective CO 2 reduction is unclear. In this work, we show that the bulk solvation behavior within the nonaqueous electrolyte can control the CO2 reduction reaction and product distribution occurring at the catalyst-electrolyte interface. We study different tetrabutylammonium (TBA) salts in two electrolyte systems with glyme ethers (e.g., 1,2 dimethoxyethane or DME) and dimethyl sulfoxide (DMSO) as a low and high dielectric constant medium, respectively. Using spectroscopic tools, we quantify the fraction of ion pairs that forms within the electrolyte. Also, we show how ion pair formation is prevalent in DME and is dependent on the anion type. More importantly, we show that as ion pair formation decreases within the electrolyte, CO 2 current densities increase, and a higher CO Faradaic efficiency is observed at low overpotentials. Meanwhile, in an electrolyte medium where the ion pair fraction does not change with the anion type (such as in DMSO), a smaller influence of solvation is observed on CO 2 current densities and product distribution. By directly coupling bulk solvation to interfacial reactions and product distribution, we showcase the importance and utility of controlling the reaction microenvironment in tuning the electrocatalytic reaction pathways. Insights gained from this work will enable novel electrolyte designs for efficient and selective CO 2 conversion to desired fuels and chemicals.

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

Two Dimensional Topology Optimization of Heat Exchangers with the Density and Level-Set Methods

We design heat exchangers using two topology optimization approaches: the density, i.e. volume fraction and level set methods. Our goal is to maximize the heat exchange between two fluids in separate channels while constraining the pressure drop across each channel. The heat exchanger is modeled with a coupled thermal-flow formulation. The flow is governed by an isothermal and incompressible Stokes-Brinkman equation and the heat transfer is governed by a convection-diffusion equation with high Peclet number. We solve one set of Stokes-Brinkman equations per fluid. Each Brinkman term in the flow equation serves to model the other phase as a solid, thereby preventing mixing. We first represent the solid and fluid phases using a volume fraction variable and apply a SIMP-like penalization in the Brinkman term to drive the optimization to a discrete design. The cost and constraint function derivatives are automatically calculated with the library pyadjoint and the optimization is performed by the Method of Moving Asymptotes. In a second optimization formulation, we use the level set approach to define the interface that separates the two fluids. Pyadjoint calculates the shape derivatives of the cost and constraint functions and the Hamilton-Jacobi advects the interface, allowing for topological changes. We present results in two dimensions and discuss the advantages and disadvantages of each approach.

42 ENGINEERING↗
DOE OSTI2025

GR-Athena++: General-relativistic Magnetohydrodynamics Simulations of Neutron Star Spacetimes

We present the extension of GR-Athena++ to general-relativistic magnetohydrodynamics (GRMHD) for applications to neutron star spacetimes. The new solver couples the constrained transport implementation of Athena++ to the Z4c formulation of the Einstein equations to simulate dynamical spacetimes with GRMHD using oct-tree adaptive mesh refinement. We consider benchmark problems for isolated and binary neutron star spacetimes demonstrating stable and convergent results at relatively low resolutions and without grid symmetries imposed. The code correctly captures magnetic field instabilities in nonrotating stars with total relative violation of the divergence-free constraint of 10 –16 . It handles evolutions with a microphysical equation of state and black hole formation in the gravitational collapse of a rapidly rotating star. For binaries, we demonstrate correctness of the evolution under the gravitational radiation reaction and show convergence of gravitational waveforms. We showcase the use of adaptive mesh refinement to resolve the Kelvin–Helmholtz instability at the collisional interface in a merger of magnetised binary neutron stars. GR-Athena++ shows strong scaling efficiencies above 80% in excess of 10 5 CPU cores and excellent weak scaling is shown up to ~5 × 10 5 CPU cores in a realistic production setup. GR-Athena++ allows for the robust simulation of GRMHD flows in strong and dynamical gravity with exa-scale computers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Terrain-Influenced Winds and Fire-Fire Interactions in Wildland Fire Simulations [Dissertation]

Ensemble-based approaches to prescribed fire planning cannot be supported by computational fluid dynamics based models like FIRETEC and the Wildland-Urban Interface Fire Dynamics Simulator (WFDS) because they are too computationally expensive and cannot leverage large eddy simulation approaches like CAWFE and WRF-SFIRE because they have too coarse of resolution. QUIC-Fire was developed to fill this gap but it cannot currently address complex terrain, that is typical for instance in the Western United States. This dissertation describes a variety of improvements made to QUIC-Fire and its various incorporated algorithms in an effort to make it a viable tool in simulating wildland and prescribed fires on terrain. The modifications made to QUIC-Fire are described in three chapters. The first chapter describes the extension of the diagnostic wind model QUIC-URB, the wind engine of QUIC-Fire, to a terrain-following coordinate system. The terraininfluenced winds it generates are analyzed and compared. In particular, this chapter presents the mathematical derivation of the wind solver leading to a linear system of equations that are solved through the successive over-relaxation method. The model is validated against a standard test used in previous works (the Askervein Hill) and against a new dataset from measurements in the Socorro Mountains, New Mexico. The terrain-following implementation captures the correct phenomenology for the isolated Askervein Hill, with a wind speed up at the top of the hill. The model agrees well with measurements on the upwind side of the peak, but overestimates speed-up on the downwind side of the hill. This is due to the inability of the model to generate flow separation and wake-eddy dynamics. On a common laptop, the divergence-free wind field is obtained in 6 s, making the solver appealing for coupled fire-atmosphere simulations. The Socorro Mountain is highly complex, with many cliff faces, peaks, and valleys. Although the model captures the magnitude and direction of inlet and outlet areas of the domain, it performs rather poorly in the valley region and in the regions near the steep cliffs. Hence, the model shows good agreement with data in areas of open sloped terrain but lacks in areas where flow separation and thermally driven effects may be present (neither effect is addressed in this work). In the second chapter the implementation of the terrain-following version of QUIC-URB into QUIC-Fire, and the necessary changes needed to include terrain are described. No changes to the underlying fire spread algorithm are made other than what is required to correctly account for the inclusion of terrain. Previously published FIRETEC results that use five different topographies that share the same centerline profile are compared to simulation results from the modified QUIC-Fire that use the same topographies and fuels. QUIC-Fire results show overall similar behaviors in terms of how the topographies affect fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all topographies are over-predicted when compared to FIRETEC. Lateral fire growth shows similar trends with FIRETEC between topographies but does not capture the increase in spread due to a diagonal interface between grassland and forested fuel region of the domain. These results suggest that there are three algorithms within QUIC-Fire that could use improvement: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB. Lastly, the third chapter describes a modification to the initial guess used for the calculation of the QUIC-URB mass-conserved wind solution during fire simulations. The modification is aimed at improving fire-fire interactions in QUIC-Fire simulations. The modification consists of using the solution from the previous timestep as the starting point for the calculation of the solution for the next timestep. Fire-fire interactions is greatly improved by the change but a new source of error is introduced. Due to how plumes are modelled in QUIC-Fire the new solution contains errors where gaps in the plume structure are present. However, these errors are mostly limited to the upper atmosphere, where they do not affect fire behavior at the surface, and their magnitude isn’t significant enough to discount the amount of new fire phenomenology now captured in QUIC-Fire with the change.

58 GEOSCIENCES↗
DOE OSTI2025

Revealing Coexisting Cu0–Cu+ Sites in Cu3N Nanoensembles for Selective CC Coupling of CO2 Under Low Overpotential

To address a long-existing debate on what copper species are responsible for efficient CC coupling, especially ethanol formation, in electrochemical CO2 reduction reaction, herein, a comprehensive study using Cu3N nanocubes with a uniform size and shape, alongside a single crystalline phase is reported. The Cu3N nanoensemble electrode has a remarkable Faradaic efficiency (FE) of 64% for ethanol production at a relatively low potential of -0.6 V versus reversible hydrogen electrode. Through in-operando X-ray absorption spectroscopy study, a dynamic phase evolution that directly correlates with changes in FE across varying applied potentials is observed. Notably, the nanoensemble with a composition of ≈71% Cu+ and 29% Cu0 is identified as being selective for ethanol formation at the low overpotential. Conversely, a predominantly metallic Cu phase formed at potentials more negative than -0.6 V favors the hydrogen evolution reaction. Density functional theory calculations at the Cu3N-Cu interface substantiate that the coexistence of Cu0-Cu+ not only energetically favors the ethanol reaction pathway but also destabilizes the intermediates for ethylene pathway.

*CH-CHO as an intermediate for ethanol production↗
DOE OSTI2021

Observation of magnetic adatom-induced Majorana vortex and its hybridization with field-induced Majorana vortex in an iron-based superconductor

Braiding Majorana zero modes is essential for fault-tolerant topological quantum computing. Iron-based superconductors with nontrivial band topology have recently emerged as a surprisingly promising platform for creating distinct Majorana zero modes in magnetic vortices in a single material and at relatively high temperatures. The magnetic field-induced Abrikosov vortex lattice makes it difficult to braid a set of Majorana zero modes or to study the coupling of a Majorana doublet due to overlapping wave functions. Here we report the observation of the proposed quantum anomalous vortex with integer quantized vortex core states and the Majorana zero mode induced by magnetic Fe adatoms deposited on the surface. We observe its hybridization with a nearby field-induced Majorana vortex in iron-based superconductor FeTe 0.55 Se 0.45 . We also observe vortex-free Yu-Shiba-Rusinov bound states at the Fe adatoms with a weaker coupling to the substrate, and discover a reversible transition between Yu-Shiba-Rusinov states and Majorana zero mode by manipulating the exchange coupling strength. The dual origin of the Majorana zero modes, from magnetic adatoms and external magnetic field, provides a new single-material platform for studying their interactions and braiding in superconductors bearing topological band structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Meshfree Multiphysics Damage Modeling of Li-ion Battery Materials

Through repeated charging and discharging cycles, the electrodes of a Li-ion battery experience significant swelling and contraction due to the movement of lithium, also known as intercalation. This research focuses specifically on the chemo-mechanical cracking in the cathode. Because cathode particles are comprised of many randomly-oriented grains, which have highly anisotropic material properties, the expansion and contraction is very non-uniform. As a result, stress concentrations tend to form between grains, which necessitates the modeling of crack propagation largely along grain boundaries and material interfaces. Chemo-mechanical damage models are generally simulated by one of two methods: the cohesive zone model (CZM) and the continuous damage model (CDM). The CZM is more accurate at capturing the sharp discontinuities of a crack but is very computationally expensive and intractable for large-scale models as a result. Conversely, the CDM is easily computed but not well-suited to easily allow for discontinuous field variables, which are inherent across a crack. This study aims to improve the CDM's ability to capture discontinuous cracks. Current versions of the CDM use the finite element method (FEM), which is one of the most widely used approaches for spatial discretization. This research investigates the use of the Reproducing Kernel Particle Method (RKPM), a meshfree method, for spatial discretization and aims to achieve a coupled chemo-mechanical crack propagation model that enhances accuracy while maintaining high computational efficiency. This model features a fully coupled, iterative electrochemistry solution, which informs the meshfree damage model of impending crack formation. To gain a better understanding of how the meshfree model compares to the current FEM model in capturing cathode crack propagation, further investigation is needed.

damage modeling↗
DOE OSTI2023

Single-cell multimodal imaging uncovers energy conversion pathways in biohybrids

Microbe–semiconductor biohybrids, which integrate microbial enzymatic synthesis with the light-harvesting capabilities of inorganic semiconductors, have emerged as promising solar-to-chemical conversion systems. Improving the electron transport at the nano–bio interface and inside cells is important for boosting conversion efficiencies, yet the underlying mechanism is challenging to study by bulk measurements owing to the heterogeneities of both constituents. Here we develop a generalizable, quantitative multimodal microscopy platform that combines multi-channel optical imaging and photocurrent mapping to probe such biohybrids down to single- to sub-cell/particle levels. We uncover and differentiate the critical roles of different hydrogenases in the lithoautotrophic bacterium Ralstonia eutropha for bioplastic formation, discover this bacterium’s surprisingly large nanoampere-level electron-uptake capability, and dissect the cross-membrane electron-transport pathways. Furthermore, this imaging platform, and the associated analytical framework, can uncover electron-transport mechanisms in various types of biohybrid, and potentially offers a means to use and engineer R. eutropha for efficient chemical production coupled with photocatalytic materials.

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

Ultrasonic waves from radial mode excitation of a piezoelectric disc on the surface of an elastic solid

In this work, we discuss the ultrasonic wave propagation characteristics in an isotropic elastic solid due to the radial mode excitation of a piezoelectric disc actuator bonded to its surface. Finite element simulations using coupled electromechanical modeling is employed to investigate the wave propagation behavior. We find that the radial mode vibrations on the surface of the elastic solid generates all the three types of ultrasonic waves: longitudinal, shear, and surface waves. The waves in the solid are comprised of a central lobe and multiple side lobes based on the frequency of the radial mode excitation. The central lobe is predominantly composed of longitudinal waves and the side lobes are composed of shear waves. While the longitudinal waves have a strong central lobe and weak side lobes that are not fully developed, shear waves consists of fully developed side lobes. In addition, we observe that the longitudinal waves have fewer side lobes inside the solid compared to the shear waves. A semi-analytical approach is presented to explain the above observations. The interface traction boundary condition between the piezoelectric disc and the elastic solid is approximated as a truncated 'Bessel' excitation over the area of the piezoelectric disc in contact with the solid. A Fourier–Bessel series expansion technique is used to investigate the wave propagation behavior from the above traction boundary condition. The results obtained explain the observations from the finite element simulations with regard to the number and distribution of side lobes pertaining to the longitudinal and shear waves generated from the radial modes of the piezoelectric disc. The proposed semi-analytical approach is general and can be applied to any arbitrary axisymmetric excitations on the surface of an elastic solid.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Proximity-enhanced valley Zeeman splitting at the WS 2 /graphene interface

The valley Zeeman physics of excitons in monolayer transition metal dichalcogenides provides valuable insight into the spin and orbital degrees of freedom inherent to these materials. Being atomically-thin materials, these degrees of freedom can be influenced by the presence of adjacent layers, due to proximity interactions that arise from wave function overlap across the 2D interface. Here, we report 60 T magnetoreflection spectroscopy of the A- and B- excitons in monolayer WS 2 , systematically encapsulated in monolayer graphene. While the observed variations of the valley Zeeman effect for the A- exciton are qualitatively in accord with expectations from the bandgap reduction and modification of the exciton binding energy due to the graphene-induced dielectric screening, the valley Zeeman effect for the B- exciton behaves markedly different. We investigate prototypical WS 2 /graphene stacks employing first-principles calculations and find that the lower conduction band of WS 2 at the <?CDATA $K/K^{^{\prime}}$?> K / K ′ valleys (the <:tex-math><?CDATA $\mathrm{CB}^-$?> C B − band) is strongly influenced by the graphene layer on the orbital level. Specifically, our detailed microscopic analysis reveals that the conduction band at the Q point of WS 2 mediates the coupling between <?CDATA $\mathrm{CB}^-$?> C B − and graphene due to resonant energy conditions and strong coupling to the Dirac cone. This leads to variations in the valley Zeeman physics of the B- exciton, consistent with the experimental observations. Our results therefore expand the consequences of proximity effects in multilayer semiconductor stacks, showing that wave function hybridization can be a multi-step energetically resonant process, with different bands mediating the interlayer interactions. Such effects can be further exploited to resonantly engineer the spin-valley degrees of freedom in van der Waals and moiré heterostructures.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polar coupling enabled nonlinear optical filtering at MoS 2 /ferroelectric heterointerfaces

Complex oxide heterointerfaces and van der Waals heterostructures present two versatile but intrinsically different platforms for exploring emergent quantum phenomena and designing new functionalities. The rich opportunity offered by the synergy between these two classes of materials, however, is yet to be charted. Here, we report an unconventional nonlinear optical filtering effect resulting from the interfacial polar alignment between monolayer MoS 2 and a neighboring ferroelectric oxide thin film. The second harmonic generation response at the heterointerface is either substantially enhanced or almost entirely quenched by an underlying ferroelectric domain wall depending on its chirality, and can be further tailored by the polar domains. Unlike the extensively studied coupling mechanisms driven by charge, spin, and lattice, the interfacial tailoring effect is solely mediated by the polar symmetry, as well explained via our density functional theory calculations, pointing to a new material strategy for the functional design of nanoscale reconfigurable optical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Zero-Strain Metal-Insulator Transition by the Local Fluctuation of Cation Dimerization

The coupled electronic and structural transitions in metal-insulator transition (MIT) hinder ultrafast switching and ultimate endurance. Decoupling these transitions and achieving a zero-strain electronic MIT can overcome the fundamental limitations of MIT in solid materials. Here, this study demonstrates that iso-valent Ti dopants in supercooled VO 2 epitaxial films cause MIT with minimal hysteresis without changing unit-cell volume and crystal symmetry. The Ti dopants in the VO 2 lattice locally alter the configuration of V-V pairs, where the long-range ordering in V-V pairs is disrupted, and the nano-domains of V-V dimers are formed. Strikingly, these local V-V dimers persist even above the electronic transition temperature (T MI ), facilitating the zero-strain electronic MIT with nanoscale structural heterogeneity. The geometrically compatible interface between insulating and metallic phases drastically enhances switching speed and endurance during electrically and optically driven zero-strain MIT. In conclusion, this discovery offers a fresh perspective on the scientific understanding of MIT and the improved functionality in terms of device speed and reliability by decoupling electronic and structural transitions.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Geometric frustration of Jahn–Teller order in the infinite-layer lattice

The Jahn–Teller effect, in which electronic configurations with energetically degenerate orbitals induce lattice distortions to lift this degeneracy, has a key role in many symmetry-lowering crystal deformations. Lattices of Jahn–Teller ions can induce a cooperative distortion, as exemplified by LaMnO 3 . Although many examples occur in octahedrally or tetrahedrally coordinated transition metal oxides due to their high orbital degeneracy, this effect has yet to be manifested for square-planar anion coordination, as found in infinite-layer copper, nickel, iron and manganese oxides. Here, for this work, we synthesize single-crystal CaCoO 2 thin films by topotactic reduction of the brownmillerite CaCoO 2.5 phase. We observe a markedly distorted infinite-layer structure, with ångström-scale displacements of the cations from their high-symmetry positions. This can be understood to originate from the Jahn–Teller degeneracy of the d xz and d yz orbitals in the d 7 electronic configuration along with substantial ligand–transition metal mixing. A complex pattern of distortions arises in a $2\sqrt2$ x $2\sqrt2$ x $1$ tetragonal supercell, reflecting the competition between an ordered Jahn–Teller effect on the CoO 2 sublattice and the geometric frustration of the associated displacements of the Ca sublattice, which are strongly coupled in the absence of apical oxygen. As a result of this competition, the CaCoO 2 structure forms an extended two-in–two-out type of Co distortion following ‘ice rules’.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Morphological and crystallographic controls in the replacement of calcite and aragonite by cerussite and otavite

Mineral replacement reactions are essential solid-fluid interactions in natural and industrial processes including metasomatism, diagenesis, and metamorphism, and sequestration of metal toxins from polluted water. In this paper, we explore the morphological evolution of aragonite and calcite (polymorphs of CaCO 3 ) during their replacement by cerussite (PbCO 3 , an isomorph of aragonite) and otavite (CdCO 3 , an isomorph of calcite) in acidic aqueous solutions (with initial pH 3 to 4). Observations using scanning electron and synchrotron-based transmission X-ray microscopies reveal the formation of pseudomorphic shells (~5 to 10 µm thick) of cerussite and otavite when Pb 2+ is reacted with either calcite or aragonite and when Cd 2+ is reacted with aragonite. The formation of pore space at the substrate-precipitate interface and within the precipitate is found to be key in promoting chemical exchange between the dissolving phase and the solution. In contrast, minimal reactivity is observed when Cd 2+ is reacted with calcite. These results demonstrate that the lack of substrate-precipitate epitaxy (either because of inconsistent symmetries or due to large lattice constant mismatch) is a critical factor that enables replacement of primary CaCO 3 minerals by Pb- or Cd-bearing secondary minerals while the presence of epitaxial relationships inhibit those reactions. When calcite and aragonite coexist, the replacement of aragonite by cerussite and otavite occurs preferentially over calcite. We postulate that the relative reactivity of calcite and aragonite to Pb 2+ and Cd 2+ can be determined by interplay between morphological and crystallographic properties of CaCO 3 materials as well as solution chemistry. Our results demonstrate that coupled dissolution-precipitation processes of carbonate minerals play an important role in controlling mobility and bioavailability of lead and cadmium in environments where acidic metal-contaminated water interacts with carbonate-rich rocks and soils.

58 GEOSCIENCES↗
DOE OSTI2026

Pore-Scale Transport Effects in Electrochemical CO 2 Reduction on Gold via Coupled Microkinetic-Transport Modeling

A pore-resolved modeling framework is developed to quantify how pore-scale transport affects the intrinsic microkinetics of CO 2 -to-CO on Au. A DFT-informed microkinetic model is coupled self-consistently to a Generalized-Modified Poisson–Nernst–Planck (GMPNP) transport description in a single, electrolyte-filled cylindrical pore, allowing local concentrations and potential to feed back into site-specific reaction rates. FIB-SEM is used to determine pore sizes within realistic electrode materials. Across pore diameters, d p = 10–6000 nm, the surface-averaged CO 2 reduction rate is systematically reduced relative to the ideal microkinetic baseline where mass transport is not accounted for; the effectiveness factor 𝜂 𝑠,CO 2 , which quantifies this ratio, decreases rapidly at more negative potentials and is about 1% near −1.0 V vs SHE due to reactant depletion. Spatial maps reveal pore-bulk alkalization that emerges at higher cathodic bias, with a small, near-wall pH dip due to electrostatic repulsion of hydroxide at the cathode interface. For a fixed aspect ratio L p /d p , narrower pores exhibit larger 𝜂 𝑠,CO 2 by shortening diffusion paths, whereas variations in the aspect ratio L p /d p play a secondary role. A dimensionless analysis (surface/bulk Damköhler numbers) delineates operating regimes. In conclusion, this work offers a concept for incorporating microkinetic models into homogenized porous-electrode models through effectiveness factors and pore-size distribution.

Au-catalyst↗
DOE OSTI2020

235-F GoldSim Fate and Transport Model: Uncertainty Quantification

Building 235-F was configured with two missions in mind: Actinide Billet Line (ABL) and the fabrication of Pu-238 oxide for space program applications. ABL produced Np-237 billets for use in SRS reactors, whereas the design process, fabrication, and examination of Pu-238 oxide powder occurred within the following areas, respectively: Plutonium Experimental Facility (PEF), Plutonium Fuel Form (PuFF), and Old Metallography Lab (OML). By 1990 production ceased and by 2006 de-inventory occurred; however, assays have shown significant holdup remains within ABL and PuFF. As a result, 235-F is a Category 2 nuclear facility, with plans to undergo deactivation and decommission (D and D) via In-Situ Disposal (ISD). The purpose of this project is to ensure United States Environmental Protection Agency (USEPA) groundwater radiation maximum contaminant level (MCL) and dosage standards are met during the D and D of 235-F by quantifying uncertainty through probabilistic modeling and evaluation of various ISD alternatives. GoldSim is a dynamic modeling software package with a graphical, object-oriented interface capable of capturing the influence of complex system input variability on probabilistic system outcomes. A GoldSim stochastic fate and transport model for 235-F was developed and matched with a PORFLOW deterministic model to simulate probabilistic release and flow of radionuclides from ABL and PuFF into the vadose zone, the Upper Three Runs (UTR) Aquifer, and UTR Creek. The GoldSim model was used to probabilistically evaluate four ISD scenarios against USEPA groundwater radiation MCLs and dosage standards. The deterministic 235-F GoldSim fate and transport model continues to be refined to match the results of the PORFLOW deterministic model to ensure the model accurately represents radionuclide movement through the groundwater system. The stochastic variables that are utilized within the GoldSim model are founded on the most current data; a conservative perspective is taken where needed. Alignment with the PORFLOW deterministic model, coupled with input stochastic variability, allows the probabilistic 235-F GoldSim model to capture the conservative breadth of possible outcomes for radionuclide fate within this particular system. This ensures that the USEPA MCLs and dosage limits hold even in the worst case scenarios.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗