Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “balance function”

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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Artificial Intelligence for Digital Security and Protections

Proper functioning of nuclear power plants relies on a mix of well-regulated human and machine-driven workflows. This regulation supports nuclear safety through a series of processes and many of the tasks that support these processes have a repetitive nature that make artificial intelligence (AI) informed by machine learning (ML) a potential aid in a variety of tasks. AI is being evaluated for activities that include inspections, fuel processing, monitoring, and other activities. The introduction of any new technology presents a potential new attack vector. In the case of AI/ML, there are many attacks that have already been discovered and over time the attacks can be expected to follow the growth pattern observed in cyber security. While future planning is necessary, current efforts need to be established now to predict the threat emergence over the next year 10 years and mitigate potential threats. Based on these observations, AI/ML will need to become trustworthy, which corresponds to techniques and procedures that emphasize AI explainability along with resilience techniques to data, algorithms, models, and systems. This kind of system robustness is the foundation for defenses against AI/ML-specific attacks. Attempting to look forward and take a broad view of capabilities provides input to research roadmaps and the ability to distill vulnerabilities into specific use cases may provide greater assistance in understanding the technology benefits while introducing new risks. The impact of current and future AI in three areas—capabilities, challenges, and recovery strategies—represents an initial attempt at balancing both.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2025

Biogas sequestration to carbon nanofibers via tandem catalytic strategies

Upgrading decentralized biogas represents a sustainable route to produce valuable products while mitigating two potent greenhouse gases, namely, methane (CH 4 ) and carbon dioxide (CO 2 ). Conventional dry reforming of CH 4 with CO 2 yields syngas with low H 2 /CO ratios (≤1) and requires high temperatures (>800 °C) to overcome equilibrium constraints and abate coke deposition, which limits commercial implementation. Here we demonstrate the conversion of biogas into value-added carbon nanofibers via reaction integration in tandem reactors, while reducing the reaction temperature, shifting equilibrium limits and yielding H 2 -enriched syngas (H 2 /CO = 2–3) as a byproduct. Experimental and theoretical insights reveal that potassium (K) modification enhances carbon nanofiber formation due to synergistic effects via a balanced interplay between KO x -induced cobalt facets and cobalt carbide species. In conclusion, the energy cost and CO 2 footprint analyses highlight the potential advantages of tandem processes for the sustainable upgrading of biogas into valuable solid carbon products.

09 BIOMASS FUELS↗
DOE OSTI2022

A Two-Level Model Predictive Control-Based Approach for Building Energy Management including Photovoltaics, Energy Storage, Solar Forecasting and Building Loads

This paper uses a two-level model predictive control-based approach for the coordinated control and energy management of an integrated system that includes photovoltaic (PV) generation, energy storage, and building loads. Novel features of the proposed local controller include (1) the ability to simultaneously manage building loads and energy storage to achieve different operational objectives such as energy efficiency, economic cost efficiency, demand response and grid optimization through the design of specific power trajectory tracking performance functionals, (2) an energy trim function that minimizes the impact of solar forecasting errors on system performance, and (3) the design of a state of charge controller that uses day-ahead forecast of solar power and building loads to intialize energy storage at the start of each day. The local controller is tested in simulation using an exemplary system with PV generation, energy storage and dispatchable building loads. Two sample days with different PV forecasts and multiple case scenarios are considered, and the performance of the algorithm in managing the real and reactive net building load trajectories and the ramp rate of PV injections into the utility network are evaluated. The simulations are based on actual forecasted and measured PV data, and the results show that the local controller meets the tracking requirements for real and reactive power within the operating constraints of the building.

14 SOLAR ENERGY↗
DOE OSTI2025

Surrogate-driven design optimization with uncertainty constraints in Monte Carlo simulations

In multi-objective design tasks, the computational cost increases rapidly when high-fidelity simulations are used to evaluate objective functions. Surrogate models help mitigate this cost by approximating the simulation output, simplifying the design process. However, under high uncertainty, surrogate models trained on noisy data can produce inaccurate predictions, as their performance depends heavily on the quality of training data. This study investigates the impact of data uncertainty on two multi-objective design problems modelled using Monte Carlo transport simulations: a neutron moderator and an ion-to-neutron converter. For each, a grid search was performed using five different tally uncertainty levels to generate training data for neural network surrogate models. These models were then optimized using NSGA-III. The recovered Pareto-fronts were analyzed across uncertainty levels: in the moderator problem, normalized hypervolume dropped from 0.886 at 1.0% uncertainty to 0.748 at 10% uncertainty, while in the converter problem it remained near 0.50 for all cases. Average simulation times were also compared to evaluate the trade-off between accuracy and computational cost. Results show that the influence of simulation uncertainty is strongly problem-dependent. In the neutron moderator case, higher uncertainties led to exaggerated objective sensitivities and distorted Pareto-fronts, reducing normalized hypervolume. In contrast, the ion-to-neutron converter task was less affected—low-fidelity simulations produced results similar to those from high-fidelity data. These findings suggest that a fixed-fidelity approach is not optimal. Surrogate models can recover the Pareto-front under noisy conditions, and multi-fidelity studies help identify suitable uncertainty levels for each problem to balance efficiency and accuracy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

Aerosol Chemical Speciation Monitor (ACSM) Composition-Dependent Collection Efficiency (CDCE) Value-Added Product Report

Aerosol particles influence the Earth’s radiation balance directly by absorbing and scattering light and indirectly by influencing cloud formation, properties, and lifetimes. Measurements of aerosol particle optical properties, mass loading, size distributions, microphysical properties, cloud formation properties, and chemical composition are important for understanding the aerosol life cycle and for validating earth system models that predict these quantities. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Aerosol Observing System (AOS) is a highly instrumented platform that measures many of these aerosol properties in situ. The Aerodyne aerosol chemical speciation monitor (ACSM) is a baseline instrument deployed in the AOS. The ACSM provides a quantitative measurement of aerosol particle chemical composition for non-refractory (operationally defined as components that evaporate on the 600ºC vaporizer) aerosol components in real time. Standard output is the mass concentration of particulate organics, nitrate, sulfate, chloride, and ammonium. One well-known limitation to the accuracy of the ACSM data is in evaluating the fraction of the ambient aerosol particles that are detected by the instrument. This quantity is referred to as the collection efficiency (CE) and is often less than unity. This is attributed to particles rebounding after impaction onto the heated vaporizer rather than being trapped, volatilized rapidly, and detected. Other factors, such as divergence of the aerosol particle beam, may also impact CE, but the particle rebound effect is the dominant factor. Scientists will sometimes assume CE = 0.5 (i.e., one half of sampled particles are detected) for ambient particles collected during field missions. However, parameterizations have been developed that express CE as a function of the measured chemical composition, referred to as the composition-dependent collection efficiency (CDCE). The physical explanation for a CDCE, supported by laboratory studies, is that rebound from the vaporizer is a function of the particle phase, with liquid-like particles “sticking” to the vaporizer and solid-like or crystalline particles “bouncing” from the vaporizer. Thus, particles with compositions are liquids-like under the measurement conditions (e.g., certain organics, acidic particles, particles enriched in nitrate) have CE close to unity while crystalline or solid particles (e.g., deliquesced ammonium sulfate) have lower CE. This value-added product (VAP) implements the procedure described by Middlebrook et al. (2012) to correct the ACSM data for the composition-dependent collection efficiency. Applying this parameterization improves the accuracy of the ACSM data and brings them into better agreement with other co-located aerosol measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

When more data hurts: Optimizing data coverage while mitigating diversity-induced underfitting in an ultrafast machine-learned potential

Machine-learned interatomic potentials (MLIPs) are becoming an essential tool in materials modeling. However, optimizing the generation of training data used to parametrize the MLIPs remains a significant challenge. This is because MLIPs can fail when encountering local environments too different from those present in the training data. The difficulty of determining a priori the environments that will be encountered during molecular dynamics simulation necessitates diverse, high-quality training data. Here, this study investigates how training data diversity affects the performance of MLIPs using the Ultra-Fast force field (UF 3 ) to model amorphous silicon nitride. We employ expert and autonomously generated data to create the training data and fit four force field variants to subsets of the data. Our findings reveal a critical balance in training data diversity: insufficient diversity hinders generalization, while excessive diversity can exceed the MLIP's learning capacity, reducing simulation accuracy. Specifically, we found that the UF 3 variant trained on a subset of the training data, in which nitrogen-rich structures were removed, offered vastly better prediction and simulation accuracy than any other variant. By comparing these UF 3 variants, we highlight the nuanced requirements for creating accurate MLIPs, emphasizing the importance of application-specific training data to achieve optimal performance in modeling complex material behaviors.

ab initio molecular dynamics↗
DOE OSTI2020

Roles of Alkali Metals and Ionic Networks in Directing the Formation of Anionic Metal–Organic Frameworks

Integrating different types of fundamental chemical interactions to direct the self-assembly processes and to control the dimension of the network formed by each type of interactions is a promising strategy to develop advanced functional materials. Here, the ionic interaction of alkali metals (Na, K, Cs) towards multi-functional 2,2’-disulfo-4,4’-oxydibenzoic acid is utilized in two solvent systems, aqueous-amide and uro-amide, to synthesize six novel anionic framework materials with dramatically different features and dimensionalities: (Na 3 [Zn 4 O(L 4- ) 3 (e-murea)(H 2 O) 3 ](H 2 N(CH 3 ) 2 +) 3 (CPM-s1), [Zn 4 O(L 4- ) 3 (DMF)]((NH 2 (CH 3 ) 2 +) 6 ) (CPM-s2), Na 6 [(Zn 2 OH 2 ) 3 (L 3- ) 6 ](H 2 O) 3 (CPM-s3), K 3 [Zn 2 (L 4- ) 2 ](H 2 N(CH 3 ) 2 +)(e-murea) 2 (H 2 O) 3 (CPM-s4), K 4 [Zn(L 3- ) 2 ](H 2 O) 2 (CPM-s5), and Cs 6 [Zn(L 4- ) 2 ](H 2 O) 2 (CPM-s6), L 4- = 2,2’-disulfo-4,4’-oxydibenzoate, DMF = N,N-dimethylformamide, e-murea = 1,3-dimethyl-1-2-imidazolidinone). In CPM-s1 to CPM-s6, the alkali-sulfonate ionic bonding networks transformed from 0D clusters with Na + to 1D chains and 2D sheets with K + , and 3D network with Cs + . These vastly differing ionic domains influenced the coordination modes and spatial arrangements of the metal-carboxylate linkers, which consequently, directed the dimensionality of the coordination networks from 3D with Na + to 1D with Cs + . CPM-s1 and CPM-s2, made from tetrameric Zn 4 O clusters are rare examples of anionic MOF-5-type structures. CPM-s3 features an unprecedented hexameric Zn 6 cluster. Even though a large portion of the charge-balancing cations in these materials would hamper hydrogen bonding networks, some of them still exhibit very high ionic conductivity (e.g. 1.25 x 10 -3 S cm -1 in CPM-s3), suggesting the potential of these materials as solid-state electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The role of thermal stratification on the co‐spectral properties of momentum transport above an Amazonian forest

The influence of thermal stratification on the turbulent kinetic energy balance has been widely studied; however, its influence on the turbulent stress remains less explored in the presence of tall vegetated canopies and less ideal micrometeorological conditions. Here, the impact of thermal stratification on turbulent momentum flux is considered in the roughness sublayer (RSL) and the atmospheric surface layer (ASL) using the Amazon Tall Tower Observatory (ATTO) in Brazil. A scalewise co‐spectral budget (CSB) model is developed using standard closure schemes for the pressure–velocity decorrelation. The CSB revealed that the co‐spectrum ${F}_{wu}\left({k}_x\right)$ between longitudinal (u') and vertical (w') velocity fluctuations is impacted by the energy spectrum of the vertical velocity ${E}_{ww}\left({k}_x\right)$ and the much less studied longitudinal heat‐flux co‐spectrum ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$, where ${\theta}_{\mathrm{v}}^{\prime }$ are temperature fluctuations and ${k}_x$ is the longitudinal wavenumber. Under stable, very stable, and dynamic–convective conditions, the scaling exponent ${F}_{wu}\left({k}_x\right)$ in for the inertial subrange (ISR) scales is dominated by ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ instead of ${E}_{ww}\left({k}_x\right)$. A near ${k}_x^{-7/3}$scaling in ${F}_{u{\theta}_{\mathrm{v}}}\left({k}_x\right)$ robust to large variations in thermal stratification is found, whereas the Kolmogorov ISR scaling for ${E}_{ww}\left({k}_x\right)\sim {k}_x^{-5/3}$ is not found. The scale‐dependent decorrelation time between u' and w' is dominated by ${\epsilon}^{-1/3}{k}_x^{-2/3}$ in the ISR, but is nearly constant for eddies larger than the vertical velocity integral scale, regardless of stability. Implications of these findings for generalized stability correction functions that are based on the turbulent stress budget instead of the turbulent kinetic energy budget are discussed.

canopy turbulence↗
DOE OSTI2023

Biologically Cooling the Planet: A BRISK Approach that Directly Connects the Ground to Deep Space

Climate change is a grave threat to our ability to thrive and prosper on Earth. Accumulation of greenhouse gases in the atmosphere stifle the Earth’s ability to cool itself, which leads to a rise of the average temperature – Global warming. The planetary cooling is exclusively radiative – thermal radiation escapes into cold outer space via a flux of photons at the low-infrared range. The penetration of these photons through Earth’s atmosphere is dictated by atmospheric transparency, which, serendipitously possesses a window centered around the peak of the Plank spectrum for a black body at about room temperature – the average temperature of Earth. This window spanning about 8 --13 μm wavelengths allows most of the photons in this range to escape into the Cosmos carrying our planet’s heat away. The efficiency of this outer heat flux – and thus planet’s cooling ability – can be controlled via emissivity – a function of wavelength which multiplies the fundamental Plank spectrum of a black body at a given temperature. In the presence of non-trivial atmospheric transparency it can be shown that radiative cooling can be enhanced with properly “designed” emissivity function. For the last 100 years, or so, people have been looking into ways to engineer materials to take advantage of such enhancement for various cooling applications, from making ice to cool clothes, and substantial progress has been achieved, including cooling in direct sunlight. The beauty if this approach is that it is entirely passive – no energy input is required. The objective of this project was to investigate whether biological systems possess such advanced cooling capacity, developed via millennia of evolution, to enhance plant’s adaptation to a hot environment, to increase cooling capability, and possibly modulate spectral emissivities on short or long time scales. We wanted to know if BRISK (Biological Radiative Sky Cooling) occurs and to what extent. It is unknown whether biological systems can perform this "spectral redistribution." It is reasonable to expect, however, that they might: in particular, desert plants under high heat stress and arid conditions might evolve this mechanism in order to maintain a tolerable internal temperature. Our intended approach was twofold: extensive literature search and supporting experimental investigations on select plants. We obtained the answers from the former, while the latter proved to be more complicated to perform with required accuracy in the timespan/budget of a small seedling project. After one year of research, we concluded that among the plant species studied there seem to be no evidence of BRISK, but it’s not the end of the story. Our planet’s biosphere is immense, and the evolution is continuing. With our preliminary results and a better understanding of what to look for further discoveries await, considering the accelerated pace of infrared imaging and spectroscopy technology development, from microscopic to satellite-based studies of plants and canopies. In the end, Earth’s energy balance is ultimately tied to radiative cooling so that further investigation into this problem is directly tied to solving the Climate Change challenge.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Marine Hydrokinetic Tidal Turbine

The U.S. Department of Energy and National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset is part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with other energy technologies. This dataset contains inputs and outputs from simulations of a floating marine hydrokinetic turbine over approximately half a tidal cycle (~6.6 hours). Inflow conditions were derived from field measurements in Alaska’s Cook Inlet and represent a tidal environment in which the current speed ramps from near 0 m/s to a peak of 3 m/s and back. The original acoustic doppler current profiler dataset is publicly available on the Marine and Hydrokinetic Data Repository. In a full tidal cycle, the flow reverses and the rotor would reorient; this reversal was not modeled. In the Cook Inlet campaign , turbulence intensity was similar in both directions. Two inflow cases are included. In the first case, labeled “raw” in the files, the measured current time series was used directly in the InflowWind module of OpenFAST. Speed and direction were applied as a function of time and elevation, uniformly in the horizontal direction. With full spatial coherence, this approach captures high turbulent variability and results in pronounced power fluctuations, so it is considered a conservative, near-worst-case representation of loading. In the second case, labeled “average” in the files, a 30-minute moving average was applied to extract the slowly varying mean speed. The residual fluctuations about this mean were used to generate spatially varying, full-field turbulence inputs with TurbSim, giving a more physically realistic representation of the inflow across the rotor disk. Two random realizations were used to produce distinct inflow conditions for two OpenFAST simulations representing a two-turbine array. The same turbulence intensity is applied across the full time series, producing larger fluctuations at the start and end, where the mean speed is low. The second case is the more appropriate framework for performance and power assessment but overpredicts turbulence at lower flow speeds and underpredicts it at higher speeds. As the floating platform moves and the rotor changes its x-position, Taylor’s frozen turbulence hypothesis used by InflowWind assumes a constant rather than a time-varying mean velocity, introducing some inaccuracy in the velocity plane sampling. The turbine modeled is the 500-kW Reference Model 1, a horizontal-axis two-bladed hydrokinetic turbine on a four-column floating semisubmersible substructure . Simulations were performed using OpenFAST v4.1 with the Reference Open Source Controller (ROSCO) v2.10. All input files required to reproduce the simulations are included. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel . This unit supports up to 2.5 MW, but NLR has only a single 1.25-MW stack. The datasets report hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. The system controls hydrogen production by varying direct current applied to the stack, from a maximum of 3,000 A to a minimum safe operating current of 300 A, or 10%. Because the current–voltage characteristic changes as the stack ages and efficiency degrades, the actual minimum safe operating power changes over time. The simulated tidal turbine time series data was translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz. Each zip file represents a single tidal electrolysis experiment and is named: {technology}_{inflow method}_{number of 500 kW tidal turbines connected} For instance, “tidal-500kW-RM1_average_2.zip” is a 6-hour experiment using the 500-kW tidal reference model, scaled by 2x (1-MW) to better match the electrolyzer maximum of 1.25MW, fed with the 30-minute moving average current case. Each zip folder contains the following files: A .csv file of raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. A .csv file combines all tidal profiles as "combined_tidal_experiments.csv." A separate experiment, “characterization_200.zip,” shows the MC250 electrolyzer steady-state response with 30-minute load steps over 5 hours and is accessible with this entry.

08 HYDROGEN↗
DOE OSTI2024

Modeling scale-up of particle coating by atomic layer deposition

Atomic layer deposition (ALD) is a promising technique to functionalize particle surfaces for energy applications including energy storage, catalysis, and decarbonization. In this work, we present a set of models of ALD particle coating to explore the transition from lab scale to manufacturing. Our models encompass the main particle coating manufacturing approaches including rotary bed, fluidized bed, and continuously vibrating reactors. These models provide key metrics, such as throughput and precursor utilization, required to evaluate the scalability of ALD manufacturing approaches and their feasibility in the context of energy applications. Our results show that designs that force the precursor to flow through fluidized particles transition faster to a transport-limited regime where throughput is maximized. They also exhibit higher precursor utilization. In the context of continuous processes, our models indicate that it is possible to achieve self-extinguishing processes with almost 100% precursor utilization. A comparison with past experimental results of ALD in fluidized bed reactors shows excellent qualitative and quantitative agreement.

25 ENERGY STORAGE↗
DOE OSTI2022

Lithological controls on soil geochemistry and clay mineralogy across Spodosols in the coastal temperate rainforest of southeast Alaska

The interaction of similar soil forming factors with different parent materials determines soil chemical weathering, influencing soil processes and properties. In the coastal temperate rainforests of southeast Alaska, Spodosols is the dominant soil order in well-drained parent materials regardless of the lithology. Yet, the role of lithology on chemical weathering, base cation depletion, and mineral transformation during pedogenesis remains elusive. Here, we established a lithosequence comprised of soils derived from tonalite, phyllite, slate, and metavolcanic rocks to test the hypothesis that despite the soils in southeast Alaska presenting similar taxonomy and morphology, the influence of lithology in chemical weathering intensity and mineralogy can be detected. We evaluated physicochemical properties, clay mineralogy, Fe-Al oxides, chemical weathering, and elemental mass balance on eleven Spodosols sampled along the lithosequence. We also propose a new weathering index, the Weathering Index for Spodosols (WISP), suitable to evaluate the weathering degree of soils with Al mobility and base-rich parent materials. This index fills a gap in traditional weathering indices that uses Al as an immobile element and/or does not evaluate the leaching of main base cations. We found the pedons across the lithosequence expressed similar physicochemical properties, with predominantly thick profiles, presenting andic properties, relatively thick spodic horizons, and thin E horizons. Podzolization imposed similar mineral transformation trends, mainly dissolution of chlorite, mica alteration to interstratified mica-vermiculite, and formation of smectite in E horizons. Kaolinite was detected in the soils from all lithologies, but only in trace amounts in phyllite soils. Fe oxides depth distribution was similar but with a higher concentration in the slate (44.5 ± 4.2 g/kg), followed by metavolcanic (31.6 ± 4.0 g/kg), phyllite (27.9 ± 3.8 g/kg), and tonalite soils (19.9 ± 5.7 g/kg) and predominantly as organometallic complexes (extracted by Na-pyrophosphate). The stronger weathering of chlorite and mica in the E horizons was reflected in larger losses of Mg and K, while Ca and Na were more depleted in the spodic horizons, suggesting plagioclase weathering. Slate soils were the most depleted in base cations, followed by phyllite, tonalite, and metavolcanic soils. The more depleted status of slate soils was supported by the higher WISP values (WISP = 62), followed by phyllite (WISP = 55), tonalite (WISP = 53), and metavolcanic soils (WISP = 45). Our results demonstrated that lithology controlled elemental depletion intensities, Fe oxide concentrations, and partially the clay fraction mineralogy in addition to podzolization acting as a dominant pedogenic process across all sites. These findings advance our understanding of the role of lithology on soil mineralogy/geochemistry that impacts critical soil functions, such as the soil carbon cycle, elemental fluxes, organo-mineral interactions, and solid-solution reactions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Toggling Stereochemical Activity through Interstitial Positioning of Cations between 2D V 2 O 5 Double Layers

The 5/6s 2 lone-pair electrons of p-block cations in their lower oxidation states are a versatile electronic and geometric structure motif that can underpin lattice anharmonicity and often engender electronic and structural instabilities that underpin the function of active elements in nonlinear optics, thermochromics, thermoelectrics, neuromorphic computing, and photocatalysis. In contrast to periodic solids where lone-pair-bearing cations are part of the structural framework, installing lone-pair-bearing cations in the interstitial sites of intercalation hosts provides a means of a systematically modulating electronic structure through the choice of the group and the period of the inserted cation while preserving the overall framework connectivity. The extent of stereochemical activity and the energy positioning of lone-pair-derived mid-gap states depend on the cation identity, stoichiometry, and strength of anion hybridization. V 2 O 5 polymorphs are versatile insertion hosts that can accommodate a broad range of s-, p-, and d-block cations. However, the insertion of lone-pair-bearing cations remains largely underexplored. In this article, we examine the implications of varying the 6s 2 cations situated in interlayer sites between condensed [V 4 O 10 ] n double layers. Systematic modulations of lattice distortions, electronic structure, and magnetic ordering are observed with increasing strength of stereochemical activity from group 12 to group 14 cations. We compare and contrast p-block-layered M x V 2 O 5 (M = Hg, Tl, and Pb) compounds and map the significance of local off-centering arising from the stereochemical activity of lone-pair cations to the emergence of filled antibonding lone-pair 6s 2 –O 2p-hybridized mid-gap states mediated by second-order Jahn–Teller distortions. Crystallographic studies of cation coordination environments and the resulting modulation of V–V interactions have been used in conjunction with variable-energy hard X-ray photoelectron spectroscopy measurements, first-principles electronic structure calculations, and crystal orbital Hamilton population analyses to decipher the origins of stereochemical activity. Magnetic susceptibility measurements reveal antiferromagnetic signatures for all the three compounds. However, the differences in V–V interactions significantly affect the energy balance of the superexchange interactions, resulting in an ordering temperature of 160 and 260 K for Hg 0.5 V 2 O 5 and δ-Tl 0.5 V 2 O 5 , respectively, as compared to 7 K for δ-Pb 0.5 V 2 O 5 . In δ-Pb 0.5 V 2 O 5 , the strong stereochemical activity of electron lone pairs and the resulting electrostatic repulsions enforce superlattice ordering, which strongly modifies the electronic localization patterns along the [V 4 O 10 ] slabs, resulting in disrupted magnetic ordering and an anomalously low ordering temperature. The results demonstrate a versatile strategy for toggling the stereochemical activity of electron lone pairs to modify the electronic structure near the Fermi level and to mediate superexchange interactions.

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

Exploring Fluoropyridine Electrolytes in Li–S Batteries: Balancing Performance and Stability across Temperatures

A novel high-donor 3-fluoropyridine (3FPy) electrolyte has been introduced for use in Li-S batteries, demonstrating an inhibition effect on the polysulfide shuttle, even without the addition of LiNO 3 . In this study, fluoropyridine electrolytes, including 2-fluoropyridine (2FPy) and 3FPy electrolytes, are studied using electrochemical analysis, mass spectrometry (MS), and high-performance liquid chromatography (HPLC) methods. Collision-induced dissociation spectra revealed that Li + preferentially solvates with different fluoropyridines, with 2FPy exhibiting a stronger interaction due to ortho-fluorine's influence, compared to 4FPy and 3FPy. However, MS and HPLC analyses showed that 2FPy is reactive with polysulfides, while 3FPy offers high solubility for polysulfides and sulfur without reacting with them at room temperature. Further, despite 3FPy performing well at room temperature, further electrochemistry studies at elevated (60 °C) and reduced (0 °C) temperatures reveal the challenges. At high temperatures, LiNO 3 is essential to suppress the polysulfide shuttle; and at low temperatures, the performance with the 3FPy electrolyte significantly lags behind that of the ether-based electrolyte.

25 ENERGY STORAGE↗
DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Soil water percolation and nutrient fluxes as a function of topographical, seasonal and soil texture variation in Central Amazonia, Brazil

Abstract Understanding soil water dynamics and transport of nutrients is challenging in tropical rainforests due to the uniqueness of several properties related to soils, vegetation and seasonality that make relating patterns found in temperate environments to tropical sites difficult. We address the need for better edaphic characterization in tropical environments by investigating soil water percolation rates and chemistry across topographic, soil texture and seasonal gradients in a mature tropical rainforest in Central Amazonia, Brazil. We utilized a passive wick flux meter (e.g., drainage lysimeter) to directly measure real‐time percolation fluxes at 60‐cm depth, and to sample a suite of chemical species across plateau, slope and valley topographic positions. We found percolation flux volume and chemical exports generally increase with decreasing elevation and clay content, which was lowest in the valley. Daily percolation flux was observed to be 2.39 ± 0.44 in plateau, 3.01 ± 0.50 in slope and 6.16 ± 0.83 mm in valley. Most solutes were present in small amounts of <1 mg L −1 , such as PO₄ 3− , Fe 2+ /Fe 3+ and Mn 2+ ; however, NO 3 − concentrations were >20 mg L −1 , even exceeding 100 mg L −1 in the valley. Based on additional isotopic analysis, we speculate high NO 3 − concentrations are partially an artefact of root decomposition following installation of the flux meters. The empirical relationships we show among percolation volume and nutrient exports under varying topographies and soil textures can improve Earth System Model performance by better constraining ecohydrological relationships to nutrient fluxes, which can in‐turn better illuminate the important factors that govern their behaviour.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Evaluating the Effectiveness of an Ultrasonic Acoustic Deterrent in Reducing Bat Fatalities at Wind Energy Facilities

This project was designed to use thermal video cameras and fatality monitoring to evaluate the effectiveness of an ultrasonic acoustic deterrent (UAD) on bat activity and mortality, respectively. Our goals were to redesign the UAD device and installation infrastructure, determine the placement on wind turbines to optimize safety, compatibility and functionality, and to compare the mortality among the following conditions: Control (deterrents off and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Deterrent (deterrents on and turbines feathered up to the manufacturer’s cut-in speed of 3.5 m/s), Curtailment (deterrents off and turbines feathered up to 5/ m/s), and combination (deterrents on and turbines feathered up to 5 m/s). The project was divided into a Feasibility Study and Comparative Study. The objectives for the Feasibility Study were to develop an installation strategy, redesign a previous iteration of a UAD to improve performance and weatherization, and test the effectiveness of the deterrents on bat activity. The Feasibility Study was intended to work out potential issues using a relatively small number of devices prior to manufacturing and installing numerous devices for a larger-scale comparative study. The division of the project into these separate studies was based on previous experience and the challenges of assessing the capabilities of an untested UAD. During the initial development and manufacturing of the UAD, NRG Systems decided to use a piezoelectric transducer rather than an electrostatic transducer, which was used by a previous device (i.e. Deaton UAD). NRG Systems conducted lab testing (i.e., IP67 or Ingress Protection) to ensure no water or dust ingress. In addition, shock/drop trials and variations in temperature exposure were conducted as part of the reliability testing. NRG Systems also developed a communications system to allow for continuous performance monitoring of the UADs. For the Feasibility Study, we installed 6 UADs on each of 2 Gamesa G90 2-MW wind turbines (Turbine 14 and 8) at the South Chestnut Wind Energy Facility, Pennsylvania and monitored activity under control and treatment (i.e. Deterrent) conditions using thermal video monitoring. There are two major sources of variation in bat activity (beyond the anticipated treatment effect): 1) environment around the turbines might inherently favor more activity at one than the other; 2) weather conditions on any given night or within season difference (e.g., migration later during the study period) might favor more activity on some nights than on other nights. Because we could only monitor two turbines on any night, we sought to control the potential influences of these two sources by alternating the turbine on which deterrents were activated each night. If there were no loss of data due to technical failures, this design would result in an equal number of deterrent and control nights at each turbine through the monitoring period, balancing the effects of both sources of variation. We compared the time bats spent and the number of events that occurred in overlapping cameras FOV as an indicator of risk, since 80% of the overlapping FOV of the cameras was in the RSA. Equipment failures, majority due to lightning, resulted in only 17 nights with useable data, with unbalanced treatment assignment within turbines and uneven distribution of treatment assignment throughout the observational period. This resulted in a confounding of treatment assignment and seasonal change. Deterrent treatment was measured at Turbine 14 on only 2 of the first 8 usable nights (spanning the period from 8/18-9/17), whereas 6 times on Turbine 8. From 9/18-927, deterrent was on at Turbine 14 on 6 of the remaining 10 nights, and 4 on Turbine 8. We recorded a total of 1,057 bats and observed a reduction in number of events and duration of events at the UAD-activated turbine when it was Turbine 14. When Turbine 8 had the UAD activated, we observed no difference in number or duration of events. Variation between turbines is not unusual and can cause issues when study designs have no true replication, i.e., multiple turbines per treatment. Within-turbine differences suggested a trend for reduced activity when UADs were activated, particularly for Turbine 14. These results may be caused by the overall higher bat activity at Turbine 8 and confounding of treatment assignment and seasonal trends. We mapped 58 bat events in 3-dimensional space (3D), 30 and 28 events during control and treatment conditions, respectively. We observed bats crossing the rotor plane under both control and treatment conditions and observed a total of 40 confirmed or near-collisions (i.e., target close to blade but no visual confirmation of a strike) out of a total of 1,491 medium and high confidence bat observations (880 control, 611 at treatment). Twice as many collisions/possible collisions were observed during control conditions. Given the challenges with the equipment and potential confounding of the data (i.e., different activity levels at the two wind turbines), we were unable to determine whether this initial turbine placement and orientation was optimal. Given no new information on how best to install the devices, we elected to use the same placement and orientation for the comparative study. For the Comparative Study, the objectives were to investigate the relative mortality rates among 4 treatments. We searched the area within 90 m of each turbine daily to recover the highest number of fresh fatalities possible. We were unable to detect a clear reduction in mortality from deterrents alone for any individual species. Surprisingly, mortality rate of the eastern red bat (Lasiurus borealis) was estimated to be 1.3–4.2 times as much when turbines were operating normally and UADs were on than when UADs were off. Reduction in mortality of all bat species combined due to curtailment of turbines was estimated to be between 0%–38%. This effect was nullified when, in addition to curtailment, UADs were on, with 95% confidence interval ranging from a 45% reduction to a 36% increase in mortality. This was likely due to the large proportion of eastern red bats in the total carcass population. Mortality of all low-frequency echolocating bats combined (i.e. hoary bat [L. cinereus], big brown bat [Eptesicus fuscus], silver-haired bat [Lasionycteris noctivagans]) relative to control was lower when curtailed (95% CI: 0%–74%), but the addition of UADs had no detectable effect (95%CI: 13%–79%). The combined treatment reduced mortality in silver-haired bats relative to control by 11%–99%, compared to curtailment (81% reduction–67% increase) or deterrent (82% reduction–67% increase) alone. Because silver-haired bats comprised a large proportion of low-frequency calling bats found during this study, a similar effect was seen for that group. The higher mortality observed for eastern red bats at UAD compared to control could have been caused by several factors, such as the effective range of the UAD, particularly at higher frequencies, behavior, positioning of the devices on the nacelle, or a combination of these. We used 3D thermal videography to compare control and UAD bat behavior from two turbines using a total of 203 3D bat-tracks across 34 nights. We recorded a similar number of bat-tracks between treatment groups, with 51% and 49% for control and UAD, respectively. Due to potential differences in bat behavior around spinning vs stationary turbine blades, we examined UAD effectiveness separately for non-operating turbines (feathered below cut-in speed of 3.5 m/s) and operating (normal operation above wind speed of 3.5 m/s). We found a higher proportion of bat-tracks at operating turbines (82%) compared to non-operating turbines (18%), although this does not account for overall time turbines were operating versus not. At non-operating turbines the UAD appears to be effective at reducing the amount of time, flight length, and number of passes through the rotor plane, compared to control. In addition, we found bats approached turbines similarly between control and treatment turbines, with 61% of control and 63% of UAD bat-tracks originating leeward of the hub. In contrast, at operating turbines, we saw little change in bat behavior in response to UADs. For example, we found an increase in the average duration of bat-tracks between non-operating and operating turbines for UAD but at control turbines we found average duration decreased once turbines became operational. Both control and UAD had a high proportion of bat-tracks that crossed the rotor plane (i.e., collision risk) originate from the windward side when turbines were operational 65% to 92%, respectively. Given that the UAD devices closest to the blades were orientated parallel to the blades, its possible bats were not exposed to the signal until they were close to the turbine blades, as suggested by the slightly higher mean duration within 5 meters of the blades for UAD turbines. Future research should consider concentrating UAD intensity on the areas of risk (i.e. blades) with enough buffer to allow bats to react to the sound before entering the rotor-swept area (RSA). In addition, investigating the potential of installing UAD units windward of the turbine blades (e.g. a hub-mounted UAD), particularly since even under control conditions, a relatively high proportion (65%) of crosses through the blade plane originated windward. 3D thermal videography provided valuable information on future testing strategies (e.g. device placement, UAD orientation) to improve UAD effectiveness when bats are at risk (i.e., operating wind turbines). Across the entire project we experienced issues with the operation and communication with the UADs. Most of the issues occurred during the Feasibility Study and were resolved prior to the Comparability Study. Additional challenges surfaced during the Comparability Study but were remedied immediately and are thought to have little impact on the results. We had logistical constraints at the project that limited our ability use traditional methods in our camera calibration. Several calibrations showed inaccurate scales, which may have been related to inadequate spatial coverage of “points” in the camera calibration volume. Because we were using actual video recordings of bats at a wind turbine, the behavior of bats could have concentrated “points” in specific areas of the turbine (i.e. leeward of nacelle), and limited “points” in other areas, resulting in camera calibration issues. We have plans to address these inconsistencies and improving the software and related methodologies by early 2020.

17 WIND ENERGY↗
DOE OSTI2023

Evaluation of Zero-Net-Rate Pumping Tests

Accurately estimating the distribution of aquifer properties is key to understanding contaminant movement in the subsurface. The distribution of aquifer properties is typically addressed using slug or constant-rate well tests, and the pros and cons of these tests are well known. Slug tests are appealing because they avoid removing contaminated water, but their results are affected by well skin and the small volume of displaced water limits the volume of aquifer that can be evaluated. Constant-rate well tests have the disadvantage of requiring disposal of potentially contaminated water, but they can generate properties that are more representative than slug tests, and they can be used to estimate well efficiency and storativity, which are difficult to characterize using slug tests. Periodic pumping tests are appealing because they have many of the advantages and few of the disadvantages of slug and constant-rate well tests. Periodic pumping tests involve cycling the pumping rate with a regular period and measuring the resulting response in monitoring wells. Some of these tests involve imposing a periodic rate on a constant mean rate. However, other tests involve moving water out and back into the well at rates that are balanced so the net rate after each period is zero. Some zero-net-rate (ZNR) well tests use rates that follow a sinusoidal pattern, whereas others use a square wave pattern that switches between a constant rate of pumping and a constant rate of injection to achieve a zero-net rate. ZNR well tests appear to be a useful compromise between slug and constant rate tests, but methods for conducting and analyzing the results from these tests have received limited evaluation. Periodic pumping tests are typically evaluated through recorded pressure responses to a disturbance source. Recent work has demonstrated that the strain field in the vadose zone shows a response to disturbances in the underlying aquifer. Various instrument designs are available that can record vertical and horizontal strain to very precise resolutions (10-9 ε). Measuring the strain in the vadose zone can present lower costs than measuring pressure in a monitoring well, so using strain could improve the resolution of ZNR tests, as well as make these tests potentially cheaper. A ZNR periodic pumping system was constructed and used to perform tests in the Clemson, South Carolina area. The system was designed to pump and inject with a periodic square-wave, so pumping occurs at a constant rate for half the period and is followed by injection at the same rate for half the period resulting in no net gain or loss of water from the aquifer. The system is designed to generate flow rates from 1 to 3.5 gallons per minute (gpm) with a capacity of 875-gallon per half period. Pressure and strain monitoring points are located in the vicinity of the pumping well. A series of 10 ZNR periodic pumping tests were conducted at the field site, with periods ranging from 2 to 540 minutes. Pressure data were collected and recorded in four monitoring wells in various locations around the pumping well. In addition, vertical and horizontal strain and tilt data were recorded in various locations around the pumping well. Traditional aquifer tests, including slug tests and constant-rate pumping tests were conducted at the field site to provide a baseline of aquifer parameter estimates that are compared the results of ZNR tests. Pressure data measured at monitoring wells during ZNR tests were analyzed utilizing an analytical solution (Streltsova, 1988) that assumes a confined aquifer and estimates hydraulic diffusivity and transmissivity using the time lag and amplitude of the pressure signal. The storativity can be separated out from the hydraulic diffusivity using the transmissivity estimate. The time-lag and amplitude of the pressure signals were estimated using a Fourier transform. The time-lag of the primary period of the pressure increases as a linear function of distance from the pumping well, and it increases as the square root of the period of the ZNR test over a range of periods spanning two orders of magnitudes. The first few harmonics follow similar trends, but the higher frequency harmonics diverge from this trend. Strain and tilt data from various instruments in the vadose and saturated zones are periodic during ZNR tests. The time lags of the strain components increase roughly linearly with distance from the well. Many of the time lags of the strain data from the vadose zone are similar to the time lags observed in the pressure data from a similar distance and pumping period. Theoretical experiments were performed to understand how pressure and strain responds to ZNR periodic pumping. The time-lags of the simulated pressure responses match the field data trends. The results were used to validate the Streltsova (1988) solution for periodic pumping to estimate hydraulic diffusivity. Hydraulic diffusivity was estimated to be 1.8 x 10-2 m2s-1 < Dh < 9.0 x 10-2 m2s-1 using pressure data from all the wells during ZNR tests. Transmissivity was estimated to be 0.8 x 10-4 m2s-1 < T < 4.3 x 10-4 m2s-1. Assuming Storativity = T/Dh gives 1.1 x 10-3 < S < 19 x10-3. Transmissivity is used to estimate the hydraulic conductivity, K, by dividing by the assumed aquifer thickness. Two constant-rate pumping tests were conducted and analyzed using two conceptual models: The Hantush (1961) solution was used to analyze data assuming confined conditions and the Neuman (1974) solution was used to analyze data assuming unconfined conditions. The Hantush (1961) solution estimated transmissivity to be 1.6 x 10-4 m2s-1 < T < 1.8 x 10-4 m2s-1 when using data from all the wells. Storativity was estimated to be 4.4 x 10-3 < S < 8.2 x 10-3. Assuming Dh = T/S gives 0.22 x 10-2 m2s-1 < Dh < 3.6 x 10-2 m2s-1. The Neuman (1974) solution estimated transmissivity to be 0.80 x 10-4 m2s-1 < T < 0.84 x 10-4 m2s-1. Total storativity (S + Sy) was estimated to be 22 x 10-3 < S < 110 x 10-3. Storativity excluding the Sy term was estimated to be 2.2 x 10-3 < S < 2.5 x 10-3. Assuming Dh = T/S gives 0.073 x 10-2 m2s-1 < S < 0.36 x 10-2 m2s-1. Two slug tests were performed on the pumping well, PW-2. The Bouwer and Rice (1976) solution for slug tests was used to analyze the data. This solution works for both confined and unconfined conditions, so both conceptual models were covered through the analysis. Hydraulic conductivity was estimated to be 1.95 x 10-6 ms-1 < K < 2.49 x 10-6 ms-1. Statistical comparison show generally no significant difference in the estimates of Dh, K, or S made using ZNR and conventional pumping tests.. However, the total storage estimated by the Neuman (1974) solution for unconfined settings is larger than that estimate using the Streltsove analysis applied to the ZNR tests. This is likely because the Neuman solution considers delayed yield from storage that the water table, whereas the Streltsova (1988) solution assumes confined conditions. The analysis also indicates that there is a statistically significant difference between parameters estimated with conventional slug tests and those measured with either constant-rate pumping tests or ZNR tests. Hydraulic diffusivity is estimated from the time-lag and the distance from the well. The time-lags of the strain data from the vadose zone are similar to the lags of the pressure data from a similar distance. This indicates that hydraulic diffusivities estimated from the strain data measured in the vadose zone would be similar to the estimates from the pressure data measured in the aquifer. There appear to be sign reversals in some of the strain data that were corrected to estimate the time lag. The dissolved oxygen (CO) concentration was measured during several ZNR tests to evaluate the feasibility of using the procedure to increase contaminant degradation kinetics that are related to DO. The results indicate that DO can be increased during the cycling of ZNR tests in some cases.

Smith-Jones, Austin↗