Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “performance bound”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Operational space for lower hybrid heating scenarios in the full tungsten environment of WEST

In tungsten—W—Environment in Steady-state Tokamak (WEST), the lower hybrid current drive (LHCD) system is key for achieving long pulse operation by providing most of the non-inductive plasma current, as well as a crucial source of electron heating. Therefore, determining the operational space for its application is fundamental. In the present study, the LHCD operational space is deeply analyzed for 0.5 MA pulses. This space is bounded by three limits: (i) the ratio of the LHCD power over density must be above a threshold to compensate tungsten radiation with enough core heating, (ii) the line-averaged density must be high enough to allow good coupling of the hybrid wave with the plasma, and (iii) fast electron ripple losses must be below a limit to avoid reaching a thermal threshold on plasma-facing components. If the tungsten radiation peak or burn-through phase is not safely overcome, a maximum electron temperature of 1.5 keV is obtained, confinement is degraded, and magnetohydrodynamic activity is frequently triggered, potentially causing a disruption. From experimental measurements and interpretative simulations, we highlight the main mechanisms that prevent the plasma from heating up during LHCD power ramp-up. Three parameters play a major role: plasma density, tungsten concentration and LHCD power deposition. A strategy to overcome this limitation is found: a precise density ramp-up performed simultaneously with the increase in LHCD power. Additionally, we show that boronization greatly facilitates the burn-through of tungsten by lowering its content during the heating phase. Finally, taking into account the three constraints given above, the LHCD operational space is determined at power ramp-up and during constant heating phases.

lower hybrid heating and current drive↗
DOE OSTI2025

Convergence Analysis of the Alternating Anderson–Picard Method for Nonlinear Fixed-Point Problems

Anderson acceleration (AA) has been widely used to solve nonlinear fixed-point problems due to its rapid convergence. This work focuses on a variant of AA in which multiple Picard iterations are performed between each AA step, referred to as the Alternating Anderson–Picard (AAP) method. Furthermore, despite introducing more “slow” Picard iterations, this method has been shown to be efficient and even more robust in both linear and nonlinear cases. However, there is a lack of theoretical analysis for AAP in the nonlinear case. In this paper, we address this gap by establishing the equivalence between AAP and a multisecant-GMRES method that uses GMRES to solve a multisecant linear system at each iteration. From this perspective, we show that AAP “converges” to the Newton-GMRES method. Specifically, as the residual approaches zero, the multisecant matrix, the approximate Jacobian inverse, the search direction, and the optimization gain of AAP converge to their counterparts in the Newton-GMRES method. These connections provide insights for analyzing the asymptotic convergence properties of AAP. Consequently, we show that AAP is locally 𝑞-linear convergent and provide an upper bound for the convergence factor of AAP. To validate the theoretical results, numerical examples are provided.

Anderson acceleration↗
DOE OSTI2023

Estimate of Gamma Dose Rates from Arrays of Fermi-1 Blanket Elements During the MEDE Process

The Enrico Fermi Atomic Power Plant Unit 1 (Fermi-1) was a sodium-cooled fast breeder reactor located in Monroe County, Michigan. The reactor was powered by a core of enriched uranium metal alloy driver fuel, which was enveloped by an axial and radial blanket material consisting of depleted uranium metal alloyed with 3 wt.% molybdenum. There are 406 axial and 559 radial irradiated sodium-bonded Fermi-1 blanket assemblies in storage at INL, totaling 34 metric tons of heavy metal. Disposal of the Fermi-1 blanket material directly into a geological repository is prohibited due to the reactive characteristic of its bond sodium. A melt drain evaporate (MEDE) process can effectively remove bond sodium from Fermi-1 blanket material. Consequently, planning is underway to apply a MEDE process to treat the 34 metric tons of heavy metal of irradiated Fermi-1 blanket material. Given the irradiated Fermi-1 blanket material’s relatively low power history and 50+ years of decay, it can be treated in a shielded glovebox. To assess the requisite shielding, the dose rates of the Fermi-1 blanket elements in various process configurations are needed. The method to perform the dose rate calculations is to first generate an average source term and associated photon source spectra for the Fermi-1 blanket material in SCALE and then use the associated photon spectra to calculate dose rates using MCNP6.2 in seven representative geometries. Dose rates in rem/h were calculated on contact (1 cm from outer geometry surface), 30 cm away from outer geometry surface, and 1 m from outer geometry surface at axial heights spanning the length of the blanket material within the geometry. The maximum average dose rate for the single Fermi-1 radial blanket element is ~0.60 rem/h, for the Fermi-1 radial blanket assembly ~1.10 rem/h, for the single Fermi-1 axial blanket element ~0.5 rem/h, for the Fermi-1 axial blanket assembly ~1.01 rem/h, for the MEDE can ~ 0.98 rem/h, for the Fermi Storage Canister ~ 0.83 rem/h, and for the MEDE cans in the DOE Standard Canister ~0.51 rem/h. Based on these maximum average values, the bounding dose rates are assumed to be 3× the average dose rates calculated for the axial blanket material and 7.6× the average dose rates calculated for the radial blanket material, based on the distribution of Fermi-1 blanket material assembly burnup.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Search for resonances in events with photon and jet final states in proton-proton collisions at $ \sqrt{s} $ = 13 TeV

A search for resonances in events with the γ+jet final state has been performed using proton-proton collision data collected at $\sqrt s$ = 13 TeV by the CMS experiment at the LHC. The total data analyzed correspond to an integrated luminosity of 138 fb –1 . Models of excited quarks and quantum black holes are considered. Using a wide-jet reconstruction for the candidate jet, the γ+jet invariant mass spectrum measured in data is examined for the presence of resonances over the standard model continuum background. The background is estimated by fitting this mass distribution with a functional form. The data exhibit no statistically significant deviations from the expected standard model background. Exclusion limits at 95% confidence level on the resonance mass and other parameters are set. Excited light-flavor quarks (excited bottom quarks) are excluded up to a mass of 6.0 (3.8) TeV. Quantum black hole production is excluded for masses up to 7.5 (5.2) TeV in the Arkani-Hamed-Dimopoulos-Dvali (Randall-Sundrum) model. These lower mass bounds are the most stringent to date among those obtained in the γ+jet final state.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

RESOLVING THE ELECTROCHEMICAL EQUATIONS OF A SOLID OXIDE FUEL CELL FOR USE IN TRANSIENT SIMULATION AND INTEGRATION INTO CYBER-PHYSICAL SYSTEMS

A major challenge with complex cyber-physical systems stems from long model computational time that creates a mismatch between the model system and the physical system. The numerical modeling of solid oxide fuel cells (SOFCs) presents particular challenges due to the highly coupled nature of the underlying equations and the multiphysics needed to fully resolve their behavior during a transient event. To this end current approaches revolve around splitting the computational efforts into resolving temperature effects and resolving electrochemical effects. Current methods employed for the transient simulation of an SOFC for implementation in the Hybrid Performance (HyPer) facility cyberphysical plant at the National Energy Technology Laboratory reveal a distinct need for accelerated results with a high degree of stability. To this aim, an investigation into the computational time for the code reveals that the underlying electrochemical algorithm takes an order of magnitude more time than its thermal counterpart and has a tendency to vary in terms of iteration time and as such a rework of the underlying system is proposed. The primary method for accelerated electrochemical algorithm solutions is to employ higher order root finding recipes for the resolution of the highly coupled electrochemical equations. This is done with the intention to reduce the overall number of subiterations necessary for resolving voltage, current density, and species concentration, properties of the fuel cell that are all directly coupled and require nested iterative approaches. The overall objective of this approach is an order of magnitude reduction in calculation time without sacrificing stability and increasing accuracy. Specific approaches involve using both bounded and unbounded techniques, such as the False Position method and the Secant method (or if applicable Newton-Raphson) respectively, the drawbacks being slower convergence for False Position and instability for the Secant or Newton-Raphson methods. Current preliminary results on simplified versions of the parent functions involved for electrochemical calculations indicate a reduction in computational steps by a factor of two for the secant method and a factor of three for Newton-Raphson. When implemented into new modified electrochemical algorithms, the results indicate a possible order of magnitude reduction in calculation time.

Arias, Jesus↗
DOE OSTI2025

Radiation-induced bowing of SiC/SiC composites under neutron flux gradients—integral experimental data for model validation

Here, the radiation-induced swelling of SiC and its composites, including strong dependencies on temperature and dose, can drive significant lateral bowing in the presence of temperature and/or dose gradients. In recent years, simulations have been performed to assess the extent of bowing in SiC composite light-water reactor (LWR) fuel cladding and boiling water reactor (BWR) channel boxes. However, to date, no integral experimental data exist to validate these models. This work provides the first experimental bowing evaluation of three ∼380 mm long SiC composite specimens irradiated under varying neutron dose gradients (∼50°C–60°C, 0.03–0.06 dpa): two tubes (∼9.8 mm diameter) and a miniature BWR channel box (∼30 mm square). The measured radiation-induced length swelling (∼0.3%–0.7% linear) was consistently 10%–21% higher than values obtained from 3D finite element structural analyses with inputs from 3D radiation transport calculations. This discrepancy could be at least partially explained by differences in dose rate (∼10 -8 dpa/s) compared to the literature data (∼10-6 dpa/s) used to establish the dose-to-swelling correlations in the model. Nevertheless, the modeled bowing magnitudes (<2 mm) obtained from finite element analyses and simple analytical equations were within the bounds of the experimental measurements for all specimens. With improved confidence in the ability to predict the structural response and measure the macroscopic deformations, future experiments will target transient bowing under neutron flux gradients at representative LWR temperatures and assess whether grid spacers can mitigate the tens of millimeters of bowing that would otherwise be expected in ∼4 m long LWR components.

bowing↗
DOE OSTI2026

Basin-scale analysis of Mokelumne River Formation for multi-well CO 2 injection

Large-scale carbon sequestration will likely require multiple projects injecting CO 2 into the same subsurface formation, raising concerns about safe operation and efficient use of storage capacity. This study evaluates the long-term response of the Mokelumne River Formation in California’s Sacramento Basin to multi-megaton CO 2 injection using three geologic models of the formation and the open-source simulator GEOS. The analysis focuses on three aspects of reservoir performance: (1) average pressure increase and dissolved CO 2 mass after 30 years for varying well counts and injection rates, (2) pressure interference in a multi-well configuration, and (3) dynamic storage capacity with identification of overpressure-prone regions. The results show that average formation pressure increases linearly with injected mass, while CO 2 dissolution exhibits mixed scaling: approximately linear with the number of wells but sublinear with injection rate, indicating that distributing injection across more wells enhances dissolution more effectively than increasing per-well rates. Pressure-interference effects are significant, with lower-permeability conditions delaying their onset but amplifying their magnitude at later times. Dynamic capacity, defined by the first occurrence of pressure exceeding the local overburden-based limit anywhere in the formation, varies across geologic models and assumed overburden pressure gradients. A lower fidelity geologic model predicts nearly twice the storage capacity of the two higher fidelity models, which consistently estimate approximately 1 Gt under the upper-bound overburden pressure gradient considered for the Sacramento Basin. In all model scenarios, overpressure develops away from injection wells, particularly in higher-elevation regions, highlighting the importance of basin-scale modelling for identifying risks beyond the immediate well vicinity.

Basin-scale↗
DOE OSTI2024

3-D Geological Modeling for Numerical Flow Simulation Studies of Gas Hydrate Reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. Here, this study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-stepping grid were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a three-dimensional (3-D) trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

58 GEOSCIENCES↗
DOE OSTI2026

Broadband hard X-ray attosecond pulses from extremely chirped electron beams

Attosecond pulses from free-electron lasers have opened the doors to atomic site-specific studies of bound electronic dynamics on their natural, sub-femtosecond timescales. Key to their success has been electron beam shaping techniques enabling the generation of sub-femtosecond current spikes with peak currents on the order of 10 kA. We demonstrate in an RF linac the generation of current spikes with extreme chirps on the order of 350 MeV/micron, directly competitive with the chirps expected from beam-driven plasma wakefield accelerators. Leveraging chirp-taper compensation, we use these highly chirped beams to generate hard X-ray attosecond pulses with bandwidths exceeding 30 eV, a factor of two beyond previous demonstrations. We simultaneously present the first explicit experimental evidence of chirp-taper compensation in an attosecond XFEL, finding that optimal tapering improves the bandwidth and pulse energy by factors of two and five, respectively, for our conditions. In addition to the immediate utility of such broadband hard X-ray pulses, electron beams with such extreme chirps can be utilized for unique new experimental modalities by performing further compression after the undulators. Such post-lasing compression can enable subsequent superradiant light emission at longer wavelengths, or direct excitation of quantum systems with the beam's intense space-charge field for unique attosecond pump-probe possibilities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Computational studies reveal Fluorine based quinolines to be potent inhibitors for proteins involved in SARS-CoV-2 assembly

World is witnessing one of the worst pandemics of this century caused by SARS-CoV-2 virus which has affected millions of individuals. Despite rapid efforts to develop vaccines and drugs for COVID-19, the disease is still not under control. Chloroquine (CQ) and Hydroxychloroquine (HCQ) are two very promising inhibitors which have shown positive effect in combating the disease in preliminary experimental studies, but their use was reduced due to severe side-effects. Here, we performed a theoretical investigation of the same by studying the binding of the molecules with SARS-COV-2 Spike protein, the complex formed by Spike and ACE2 human receptor and a human serine protease TMPRSS2 which aids in cleavage of the Spike protein to initiate the viral activation in the body. Both the molecules had shown very good docking energies in the range of -6kcal/mol. Subsequently, we did a high throughput screening for other potential quinoline candidates which could be used as inhibitors. From the large pool of ligand candidates, we shortlisted the top three ligands (binding energy -8kcal/mol). We tested the stability of the docked complexes by running Molecular Dynamics (MD) simulations where we observed the stability of the quinoline analogues with the Spike-ACE2 and TMPRSS2 nevertheless the quinolines were not stable with the Spike protein alone. Thus, although the inhibitors bond very well with the protein molecules their intrinsic binding affinity depends on the protein dynamics. Moreover, the quinolines were stable when bound to electronegative pockets of Spike-ACE2 or TMPRSS2 but not with Viral Spike protein. We also observed that a Fluoride based compound: 3-[3-(Trifluoromethyl)phenyl]quinoline helps the inhibitor to bind with both Spike-ACE2 and TMPRSS2 with equal probability. The molecular details presented in this study would be very useful for developing quinoline based drugs for COVID-19 treatment.

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

Thermal Modeling of Advanced Test Reactor Fuel in a Generalized Dry Storage System Under Hypothetical Accident Conditions

Work was performed to investigate the thermal behavior of a sealed aluminum-clad spent nuclear fuel (ASNF) dry storage configuration under hypothetical accident conditions during long-term storage. The considered system consists of a concrete dry storage overpack; a welded over canister, backfilled with argon; and nine Department of Energy (DOE) standardized canisters (DOESCs) loaded with ASNF, backfilled with helium gas. One technical concern associated with the long-term storage of ASNF includes radiolytic hydrogen generation. The yield of large quantities of hydrogen could lead to DOESC over pressurization or a flammable internal atmosphere. While flammability concerns can be resolved by preventing the ingress of oxygen, the canister pressure is partially controlled by the atmosphere temperature. This memorandum summarizes the results of modeling and simulation (M&S) work completed with Star CCM+ (a computational fluid dynamics [CFD] software) considering the thermal effects of a fully engulfing flame and another scenario with the vent ports of the concrete dry storage overpack completely blocked. The goal was to determine the bounding maximum temperature of the DOESC internal atmosphere to evaluate a worst-case pressure scenario, as well as the critical vent port blockage durations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Heterogeneous Supramolecular Catalysis through Immobilization of Anionic M 4 L 6 Assemblies on Cationic Polymers

Although most of the currently developed supramolecular catalysts that emulate enzymatic reactivity with unique selectivity and activity through specific host-guest interactions work under homogeneous conditions, enzymes in nature can operate under heterogeneous conditions as membrane-bound enzymes. In order to develop such a heterogeneous system, an immobilized chiral supramolecular cluster Ga 4 1 6 ( 2 ) was introduced into cross-linked polymers with cationic functionalities. These heterogeneous supramolecular catalysts were used in aza-Prins and aza-Cope reactions and successfully applied to continuous-flow reactions. Additionally, they showed high durability and maintained high turnovers for long periods of time. In sharp contrast to the majority of examples of heterogenized homogeneous catalysts, the newly developed catalysts showed enhanced activity and robustness compared to those exhibited by the corresponding soluble cluster catalyst. An enantioenriched cluster was also immobilized to enable asymmetric catalysis, and activity and enantioselectivity of the supported chiral catalyst were maintained during recovery and reuse experiments and under a continuous-flow process. Significantly, the structure of the ammonium cations in the polymers affected stability, reactivity, and enantioselectivity, which is consistent with the hypothesis that the cationic moieties in the polymer support interact with cluster as an exohedral protecting shell, thereby influencing their catalytic performance.

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

Promise of Graph Sparsification and Decomposition for Noise Reduction in QAOA: Analysis for Trapped-Ion Compilations

We develop new approximate compilation schemes that significantly reduce the expense of compiling the Quantum Approximate Optimization Algorithm (QAOA) for solving the Max-Cut problem. Our main focus is on compilation with trapped-ion simulators using Pauli-X operations and all-to-all Ising Hamiltonian HIsing evolution generated by Molmer-Sorensen or optical dipole force interactions, though some of our results also apply to standard gate-based compilations. Our results are based on principles of graph sparsification and decomposition; the former reduces the number of edges in a graph while maintaining its cut structure, while the latter breaks a weighted graph into a small number of unweighted graphs. Though these techniques have been used as heuristics in various hybrid quantum algorithms, there have been no guarantees on their performance, to the best of our knowledge. This work provides the first provable guarantees using sparsification and decomposition to improve quantum noise resilience and reduce quantum circuit complexity. For quantum hardware that uses edge-by-edge QAOA compilations, sparsification leads to a direct reduction in circuit complexity. For trapped-ion quantum simulators implementing all-to-all HIsing pulses, we show that for a (1−ϵ) factor loss in the Max-Cut approximation (ϵ>0), our compilations improve the (worst-case) number of HIsing pulses from O(n2) to O(nlog(n/ϵ)) and the (worst-case) number of Pauli-X bit flips from O(n2) to O(nlog(n/ϵ)ϵ2) for n-node graphs. This is an asymptotic improvement for any constant ϵ>0. We demonstrate that significant improvements to the approximation ratio are obtained using decomposition in simulated trapped-ion experiments with dephasing noise. We further present a generic argument showing that sparsification results in an exponentially improved circuit fidelity lower bound in digital computing schemes based on one- and two-qubit gates, which are relevant to a wide variety of hardwares such as superconducting qubits and certain neutral atom or trapped ion setups, and more sophisticated noise models. We anticipate these approximate compilation techniques will be useful tools in a variety of future quantum computing experiments.

Moondra, Jai [Georgia Institute of Technology]↗
DOE OSTI2020

Thermal conductivity and its relation to atomic structure for symmetrical tilt grain boundaries in silicon

We perform a systematic study of thermal resistance and conductance of tilt grain boundaries (GBs) in Si using classical molecular dynamics. The GBs studied are naturally divided into three groups according to the structural units forming the GB core. In this work we find that, within each group, the GB thermal conductivity strongly correlates with the excess GB energy. All three groups predict nearly the same GB conductivity extrapolated to the high-energy limit. This limiting value is close to the thermal conductivity of amorphous Si, suggesting similar heat transport mechanisms. While the lattice thermal conductivity decreases with temperature, the GB conductivity slightly increases. However, at high temperatures it turns over and starts decreasing if the GB structure undergoes a premelting transformation. Analysis of vibrational spectra of GBs resolved along different directions sheds light on the mechanisms of their thermal resistance. The existence of alternating tensile and compressive atomic environments in the GB core gives rise to localized vibrational modes, frequency gaps creating acoustic mismatch with lattice phonons, and anharmonic vibrations of loosely bound atoms residing in open atomic environments.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Entropic order parameters in weakly coupled gauge theories

The entropic order parameters measure in a universal geometric way the statistics of non-local operators responsible for generalized symmetries. In this article, we compute entropic order parameters in weakly coupled gauge theories. To perform this computation, the natural route of evaluating expectation values of physical (smeared) non-local operators is prevented by known difficulties in constructing suitable smeared Wilson loops. We circumvent this problem by studying the smeared non-local class operators in the enlarged non-gauge invariant Hilbert space. This provides a generic approach for smeared operators in gauge theories and explicit formulas at weak coupling. In this approach, the Wilson and ’t Hooft loops are labeled by the full weight and co-weight lattices respectively. We study generic Lie groups and discuss couplings with matter fields. Smeared magnetic operators, as opposed to the usual infinitely thin ones, have expectation values that approach one at weak coupling. The corresponding entropic order parameter saturates to its maximum topological value, except for an exponentially small correction, which we compute. On the other hand, smeared ’t Hooft loops and their entropic disorder parameter are exponentially small. We verify that both behaviors match the certainty relation for the relative entropies. In particular, we find upper and lower bounds (that differ by a factor of 2) for the exact coefficient of the linear perimeter law for thin loops at weak coupling. This coefficient is unphysical/non-universal for line operators. We end with some comments regarding the RG flows of entropic parameters through perturbative beta functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Silicate minerals as a direct source of limiting nutrients: Siderophore synthesis and uptake promote ferric iron bioavailability from olivine and microbial growth

Iron is a micronutrient critical to fundamental biological processes including respiration and photosynthesis, and it can therefore impact primary and heterotrophic productivity. Yet in oxic environments, iron is highly insoluble, rendering it, in principle, unavailable as a nutrient for biological growth. Life has “solved” this problem via the invention of iron chelates, known as siderophores, that keep iron available for microbial productivity. In this work, we examined the impact of siderophore synthesis on the speciation, mobility, and bioavailability of iron from rock-forming silicate minerals—shedding new light on the mechanisms by which microbes use mineral substrates to support primary productivity, as well as the consequent effects on silicate dissolution. Growth experiments were performed with Shewanella oneidensis MR-1 in an oxic, iron-depleted minimal medium, amended with olivine minerals as the sole source of iron. Experiments included the wild-type strain MR-1, and a siderophore synthesis gene deletion mutant strain (ΔMR-1). Relative to MR-1, ΔMR-1 exhibited a very pronounced growth penalty and an extended lag phase. However, substantial growth of ΔMR-1, comparable to MR-1 growth, was observed when the mutant strain was provided with siderophores in the form of either filtrate from a well-grown MR-1 culture, or commercially available deferoxamine. These observations suggest that siderophores are critical for S. oneidensis to acquire iron from olivine. Growth-limiting concentrations of deferoxamine amendments were observed to be ≤5–10 µM, concentrations significantly lower than previously recorded as necessary to impact mineral dissolution rates. X-ray photoelectric spectroscopy analyses of the incubated olivine surfaces suggest that siderophores deplete mineral surface layers of ferric iron. Combined, these results demonstrate that low micromolar concentrations of siderophores can effectively mobilize iron bound within silicate minerals, supporting very significant biological growth in limiting environments. The specific mechanism would involve siderophores removing a protective layer of nanometer-thick iron oxides, enhancing silicate dissolution and nutrient bioavailability.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2022

Mercury Removal from Contaminated Water by Wood-Based Biochar Depends on Natural Organic Matter and Ionic Composition

Biochars can remove potentially toxic elements, such as inorganic mercury [Hg(II)] from contaminated waters. However, their performance in complex water matrices is rarely investigated, and the combined roles of natural organic matter (NOM) and ionic composition in the removal of Hg(II) by biochar remain unclear. Here, we investigate the influence of NOM and major ions such as chloride (Cl – ), nitrate (NO 3 – ), calcium (Ca 2+ ), and sodium (Na + ) on Hg(II) removal by a wood-based biochar (SWP700). Multiple sorption sites containing sulfur (S) were located within the porous SWP700. In the absence of NOM, Hg(II) removal was driven by these sites. Ca 2+ bridging was important in enhancing removal of negatively charged Hg(II)-chloro complexes. In the presence of NOM, formation of soluble Hg-NOM complexes (as seen from speciation calculations), which have limited access to biochar pores, suppressed Hg(II) removal, but Cl – and Ca 2+ could still facilitate it. The ability of Ca 2+ to aggregate NOM, including Hg-NOM complexes, promoted Hg(II) removal from the dissolved fraction (<0.45 μm). Hg(II) removal in the presence of Cl – followed a stepwise mechanism. Weakly bound oxygen functional groups in NOM were outcompeted by Cl – , forming smaller-sized Hg(II)-chloro complexes, which could access additional intraparticle sorption sites. Therein, Cl – was outcompeted by S, which finally immobilized Hg(II) in SWP700 as confirmed by extended X-ray absorption fine structure spectroscopy. We conclude that in NOM containing oxic waters, with relatively high molar ratios of Cl – : NOM and Ca 2+ : NOM, Hg(II) removal can still be effective with SWP700.

54 ENVIRONMENTAL SCIENCES↗