Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “scale coupling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,279 records · Page 71
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Photoinduced Dirac semimetal in ZrTe5

Abstract Novel phases of matter with unique properties that emerge from quantum and topological protection present an important thrust of modern research. Of particular interest is to engineer these phases on demand using ultrafast external stimuli, such as photoexcitation, which offers prospects of their integration into future devices compatible with optical communication and information technology. Here, we use MeV Ultrafast Electron Diffraction (UED) to show how a transient three-dimensional (3D) Dirac semimetal state can be induced by a femtosecond laser pulse in a topological insulator ZrTe 5 . We observe marked changes in Bragg diffraction, which are characteristic of bond distortions in the photoinduced state. Using the atomic positions refined from the UED, we perform density functional theory (DFT) analysis of the electronic band structure. Our results reveal that the equilibrium state of ZrTe 5 is a topological insulator with a small band gap of ~ 25 meV, consistent with angle-resolved photoemission (ARPES) experiments. However, the gap is closed in the presence of strong spin-orbit coupling (SOC) in the photoinduced transient state, where massless Dirac fermions emerge in the chiral band structure. The time scale of the relaxation dynamics to the transient Dirac semimetal state is remarkably long, τ ~ 160 ps, which is two orders of magnitude longer than the conventional phonon-driven structural relaxation. The long relaxation is consistent with the vanishing density of states in Dirac spectrum and slow spin-repolarization of the SOC-controlled band structure accompanying the emergence of Dirac fermions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

RF sheath induced sputtering on Proto-MPEX. I. Sheath equivalent dielectric layer for modeling the RF sheath

The pulsed linear plasma device Prototype Material Plasma Exposure eXperiment (Proto-MPEX) uses a radio frequency (RF) helicon antenna with an aluminum nitride ceramic window for plasma production. The RF sheath created under the helicon antenna is sufficient to cause ion impact energies to be greater than the sputtering threshold of the AlN helicon window material and for impurities to be created. In this study, we investigate the RF sheath on the inner diameter of the helicon window and its impact on the impurity production rates in Proto-MPEX. Three models—a 3D COMSOL finite element RF model of the Proto-MPEX helicon region, a rectified DC sheath potential model, and the 3D Global Impurity TRansport code—are coupled together to study impurity production and transportation. This novel method of impurity generation and transport modeling spans length scales ranging from the sheath (millimeters or less) up to the full device (meters) and can be applied to other radio frequency sources and antennas in a wide range of plasma physics studies, including basic plasmas, low-temperature processing plasmas, plasma thrusters, and fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

From Spin to Pseudospin Symmetry: The Origin of Magic Numbers in Nuclear Structure

Magic numbers lie at the heart of nuclear structure, reflecting enhanced stability in nuclei with closed shells. While the emergence of magic numbers beyond 20 is commonly attributed to strong spin-orbit coupling, the microscopic origin of the spin-orbit potential remains elusive, owing to its dependence on the resolution scale and renormalization scheme of nuclear forces. Here, we investigate the evolution of nuclear shell structure with varying momentum resolution in nuclear interactions derived from chiral effective field theory, using the similarity renormalization group to link different scales. We uncover a novel transition from spin symmetry to pseudospin symmetry as the resolution scale decreases, during which magic numbers emerge naturally. A similar pattern is found in calculations using relativistic one-boson-exchange potentials, underscoring the robustness of the phenomenon. This establishes a direct connection between realistic nuclear forces with a high resolution scale and effective nuclear forces at coarse-grained scales, offering a first-principles explanation for the origin of magic numbers and pseudospin symmetry in nuclear shell structure and new insights into the structure of exotic nuclei far from stability

Energy levels↗
DOE OSTI2025

Time-delayed gamma-ray signatures of heavy axions from core-collapse supernovae

Heavy axions that couple to both quantum electrodynamics and quantum chromodynamics with masses on the order of MeV–GeV and high-scale decay constants in excess of ∼ 10 8 GeV may arise generically in, e.g., axiverse constructions. In this work we provide the most sensitive search to date for the existence of such heavy axions using Fermi-LAT data toward four recent supernovae (SN): Cassiopea A, SN1987A, SN2023ixf, and SN2024ggi. We account for heavy axion production in the proto-neutron-star cores through nuclear and electromagnetic processes and then the subsequent decay of the axions into photons. While previous works have searched for gamma rays from SN1987A using the Solar Maximum Mission that observed SN1987A during the SN itself, we show that using Fermi Large Area Telescope data provides an approximately 5 orders of magnitude improvement in flux sensitivity for axions with lifetimes larger than around 10 yrs. We find no evidence for heavy axions and exclude large regions of previously unexplored parameter space.

Axions↗
DOE OSTI2022

Quantifying local rearrangements in three-dimensional granular materials: Rearrangement measures, correlations, and relationship to stresses

Quantifying the ways in which local particle rearrangements contribute to macroscopic plasticity is one of the fundamental pursuits of granular mechanics and soft matter physics. Here we examine local rearrangements that occur naturally during the deformation of three samples of 3D granular materials subjected to distinct boundary conditions by employing in situ x-ray measurements of particle-resolved structure and stress. We focus on five distinct rearrangement measures, their statistics, interrelationships, contributions to macroscopic deformation, repeatability, and dependence on local structure and stress. Our most significant findings are that local rearrangements (1) are correlated on a scale of three to four particle diameters, (2) exhibit volumetric strain-shear strain and nonaffine displacement-rotation coupling, (3) exhibit correlations that suggest either rearrangement repeatability or that rearrangements span multiple steps of incremental sample strain, and (4) show little dependence on local stress but correlate with quantities describing local structure, such as porosity. Our results are presented in the context of relevant plasticity theories and are consistent with recent findings suggesting that local structure may play at least as important of a role as local stress in determining the nature of local rearrangements.

58 GEOSCIENCES↗
DOE OSTI2020

Uncovering Non-Fermi-Liquid Behavior in Hund Metals: Conformal Field Theory Analysis of an SU(2)×SU(3) Spin-Orbital Kondo Model

Hund metals have attracted attention in recent years due to their unconventional superconductivity, which supposedly originates from non-Fermi-liquid (NFL) properties of the normal state. When studying Hund metals using dynamical mean-field theory, one arrives at a self-consistent “Hund impurity problem” involving a multiorbital quantum impurity with nonzero Hund coupling interacting with a metallic bath. If its spin and orbital degrees of freedom are screened at different energy scales, T sp < T orb , the intermediate energy window is governed by a novel NFL fixed point, whose nature had not yet been clarified. We resolve this problem by providing an analytical solution of a paradigmatic example of a Hund impurity problem, involving two spin and three orbital degrees of freedom. To this end, we combine a state-of-the-art implementation of the numerical renormalization group, capable of exploiting non-Abelian symmetries, with a generalization of Affleck and Ludwig’s conformal field theory (CFT) approach for multichannel Kondo models. We characterize the NFL fixed point of Hund metals in detail for a Kondo model with an impurity forming an SU ( 2 ) × SU ( 3 ) spin-orbital multiplet, tuned such that the NFL energy window is very wide. The impurity’s spin and orbital susceptibilities then exhibit striking power-law behavior, which we explain using CFT arguments. We find excellent agreement between CFT predictions and numerical renormalization group results. Our main physical conclusion is that the regime of spin-orbital separation, where orbital degrees of freedom have been screened but spin degrees of freedom have not, features anomalously strong local spin fluctuations: the impurity susceptibility increases as χ sp imp ~ ω – γ , with γ > 1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Urbanization-induced land and aerosol impacts on storm propagation and hail characteristics

Changes in land surface and aerosol characteristics from urbanization can affect dynamic and microphysical properties of severe storms, thus affecting hazardous weather events resulting from such storms such as hail and tornado. We examine the joint and individual effects of urban land and anthropogenic aerosols of Kansas City on a severe convective storm observed during the 2015 Plains Elevated Convection At Night (PECAN) field campaign, focusing on storm evolution, convective intensity, and hail characteristics. The simulations are carried out at the cloud-resolving scale (1 km) using a version of WRF-Chem in which the spectral-bin microphysics (SBM) is coupled with the Model for Simulating Aerosol Interactions and Chemistry (MOSAIC). It is found that the urban land effect of Kansas City initiated a much stronger convective cell and the storm got further intensified when interacting with stronger turbulence induced by the urban land. The urban land effect also changed the storm path by diverting the storm toward the city mainly resulting from enhanced urban land-induced convergence in the urban area and around the urban-rural boundaries. The joint effect of urban land and anthropogenic aerosols enhances occurrences of both severe hail and significant severe hail by ~ 20% by enhancing hail formation and growth from riming. Overall the urban land effect on convective intensity and hail is relatively larger than the anthropogenic aerosol effect, but the joint effect is much notable than either of the individual effects, emphasizing the importance to consider both effects in evaluating urbanization effects.

Lin, Yun↗
DOE OSTI2021

Conversion efficiency of soliton Kerr combs

We investigate the conversion efficiency (CE) of soliton modelocked Kerr frequency combs. Our analysis reveals three distinct scaling regimes of CE with the cavity free spectral range (FSR), which depends on the relative contributions of the coupling and propagation loss to the total cavity loss. Our measurements, for the case of critical coupling, verify our theoretical prediction over a range of FSRs and pump powers. Our numerical simulations also indicate that mode crossings have an adverse effect on the achievable CE. Our results indicate that microresonator combs operating with spacings in the electronically detectable regime are highly inefficient, which could have implications for integrated Kerr comb devices.

Jang, Jae K. (ORCID:0000000187786058)↗
DOE OSTI2021

Frontier nonequilibrium materials science enabled by ultrafast electron methods

Most innovation in materials science and engineering resides in our ability to understand and control the intimate relationship between the structure of materials and their properties. Conventionally, the only route to discovering innovative new material properties has been to explore the structural and compositional phase space that is accessible at (or near) thermodynamic equilibrium. The characterization, manipulation and, ultimately, control of material properties far from equilibrium offers almost completely untapped possibilities for uncovering novel states and phases of materials. Investigations of materials far from equilibrium require the development of new techniques, which are capable of following dynamic processes in materials with extreme spatiotemporal resolution. Thus, ultrafast electron-based methods have become a major new frontier in materials science due to the capability of following dynamics on time scales as short as femtoseconds with the high spatial resolution and sensitivity afforded by electrons. The articles in this issue provide an exciting cross section of the novel research directions that have been enabled at this frontier. Imaging on ultrafast time scales carries critical information on phase transitions, fundamental processes involved in light-harvesting, magnetic, and plasmonic dynamics, the coupling of electronic and nuclear degrees of freedom in materials and has revealed previously hidden, nonequilibrium metastable states of matter that have no equilibrium analog.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A tail of eternal inflation

Non-trivial inflaton self-interactions can yield calculable signatures of primordial non-Gaussianity that are measurable in cosmic surveys. We calculate the non-Gaussian corrections to Stochastic Inflation within the framework of Soft de Sitter Effective Theory, from which we derive the associated probability distribution for the scalar fluctuations. As a consequence of this new result, we show that the phase transition to slow-roll eternal inflation is often incalculable in these models. Instead, this transition is sensitive to the non-Gaussian tail of the distribution of scalar fluctuations, which probes physics inside the horizon, potentially beyond the cutoff scale of the Effective Field Theory of Inflation. We delineate the parameter space consistent with current observations and weak coupling at horizon crossing in which the large fluctuations relevant for eternal inflation can only be determined by appealing to a UV completion. We also argue that this breakdown of the perturbative description is required for the de Sitter entropy to reflect the number of de Sitter microstates.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Domestication of Algae for Increasing Biomass Productivity

Microalgae cultivation processes have been developed for the production of a variety of bioproducts, however currently only a few species are used in commercial applications. Their domestication, that is strain improvements, is still in its infancy, with major advances required, specifically to maximize biomass productivity a limiting factor in microalgae production. This requires a deep understanding of algal biology, in particular to develop superior strains without the need of genetic technologies that would require lengthy regulatory permits, and often limit consumer acceptance. Adaptive Laboratory Evolution techniques, alone or in conjunction with sexual recombination, can allow for rapid develop of improved strains and their industrial production. Light harvesting antenna reduction has been a major approach to achieve increased photon utilization efficiency by cultures operating under full sunlight conditions due to higher light saturation levels, allowing for higher productivities under outdoor conditions. Decades of research yielded some promising results under controlled conditions with a few specific mutant strains. However, these failed to achieve the anticipated higher productivities in actual algal mass cultures, in part due to the inability of single mutations to overcome photoinhibition, reactive oxygen species, and other pleiotropic impacts on the complex metabolic processes of photosynthesis. Higher productivity strains will require multiple genetic improvements. We report on recent Adaptive Laboratory Evolution with the green alga Scenedesmus obliquus resulting in higher biomass productivity in open pond cultivation. Coupling our approach with sexual recombination and genome sequencing provides a path to algal domestication suitable for large-scale, low-cost biomass production.

09 BIOMASS FUELS↗
DOE OSTI2025

Representing lateral groundwater flow from land to river in Earth system models

Lateral groundwater flow (LGF) is an important hydrologic process in controlling water table dynamics. Due to the relatively coarse spatial resolutions of land surface models, the representation of this process is often overlooked or overly simplified. In this study, we developed a hillslope-based lateral groundwater flow model. Specifically, we first developed a hillslope definition model based on an existing watershed delineation model to represent the subgrid spatial variability in topography. Building upon this hillslope definition, we then developed a physical-based lateral groundwater flow using Darcy’s equation. This model explicitly considers the relationships between the groundwater table along the hillslope and the river water table levels. We coupled this intra-grid model to the land component (E3SM Land Model: ELM) and river component (MOdel for Scale Adaptive River Transport: MOSART) of the Energy Exascale Earth System Model (E3SM). We tested both the hillslope definition model and the lateral groundwater flow model and performed sensitivity experiments using different configurations. Simulations for a single grid cell at 0.5°×0.5° within the Amazon basin show that the definition of hillslope is the key to modeling lateral flow processes and the runoff partition between surface and subsurface can be dramatically changed using the hillslope approach. Although our method provides a pathway to improve the lateral flow process, future improvements are needed to better capture the subgrid structure to account for the spatial variability in hillslopes within the simulated grid of land surface models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Graph-Based Representations and Applications to Process Simulation

Rapid and robust convergence of a process flowsheet is critical to enable large-scale simulations that address core scientific questions related to process design, optimization, and sustainability. However, due to the highly coupled and nonlinear nature of chemical processes, efficiently solving a flowsheet remains a challenge. In this work, we show that graph representations of the underlying physical phenomena in unit operations may help identify potential avenues to systematically reformulate the network of equations and enable more robust topology-based convergence of flowsheets. To this end, we developed graph abstractions of the governing equations of vapor-liquid and liquid-liquid equilibrium separation equipment. These graph abstractions consist of a mesh of interconnected variable nodes and equation nodes that are systematically generated through PhenomeNode, a new open-source library in Python developed in this study. We show that partitioning the graph into separate mass, energy, and equilibrium subgraphs can help decouple nonlinearities and guide decomposition algorithms. By employing the graph abstraction on an industrial separation process for separating glacial acetic acid from water, we implemented a new block decomposition scheme in BioSTEAM and demonstrated that this can accelerate convergence over a traditional sequential modular approach.

Distillation↗
DOE OSTI2022

Development Of Thermodynamic and Kinetic Simulation Tools and Testing Procedures for Enhanced Durability of Concrete Containing Industrial By-Products

This project developed screening tools that enable evaluation of alternative cementitious binders that create concretes to significantly reduce energy and emissions while remaining cost competitive on both initial and long-term costs. The team began with the viewpoint that acceptance of new cementitious binder products has a substantially greater chance of successful implementation when capital investment is not excessive and the end product has customers that have experience and resources to use this. As such, a binder system that is based on portland cement with blended using industrial by-products (alternative cementitious materials) has the potential for dramatic and meaningful impact. The team has focused on developing implementable solutions in specifications and current practice. This however requires three main factors: 1) ability to screen byproducts and alternative materials for success, 2) ability to ‘treat’ materials chemically to enhance kinetics, and 3) ability to provide predictions of performance of both binders and concrete from first principles. The project developed/refined a state of the art and scientifically based screening test for SCM called the pozzolanic reactivity test. The team developed kinetic models to simulate these materials as well as experimental approaches to alter selected reactions. Simulation tools were developed that enable the performance of concrete to be predicted based on the chemistry and reactivity of the cement and alternative SCM. Specifically, this project: • enhanced the kinetic reactivity models for use in multi-scale computational programs that use thermodynamics to predict reaction products. 2 • developed scaling models to extend thermodynamic modeling to link these models with pore structure. This enables strength, transport property, and coupled transport prediction. • developed tools to predict performance of cementitious materials using the pozzolanic reactivity test and chemical composition. The predicted properties are consistent with AASHTO R101 and the CEB-FIP model code and can be measured using associated test procedures. • demonstrated mechanical and fracture based modeling tools that thermodynamic predictions and inputs to predict concrete service life. These results have been used to demonstrate the value of enabling specifications to include ASTM C 595 cement as well ASTM C150 cement. In addition, these products are being used to expedite the evaluation of alternative SCM to aide in determining which materials have potential value and what ‘compositions’ of blended cements merit further investment.

42 ENGINEERING↗
DOE OSTI2020

GPS-Based Gamma Survey for Characterizing and Decommissioning NORM Sites - 20389

Gamma survey techniques are an especially powerful decommissioning tool at naturally occurring radioactive material (NORM) sites due to both the low cost to obtain data over a large spatial scale and the abundance of gamma emitters in the uranium and thorium decay series. Gamma surveys are executed by coupling a detector - most often a sodium iodide crystal - to a global positioning system (GPS), then reporting a location and gross gamma reading coincidentally to a data logger. Systems may be carried by workers or mounted to a car, all-terrain vehicle, or unmanned aerial system (UAS). The resulting data set provides a high-resolution but low precision map of the gamma radiation field over the area surveyed. Frequently this map is also correlated to soil concentrations of NORM radionuclides (most often, Ra-226) and/or exposure rate. Gamma survey parameters such as movement speed, transect spacing, and data logging frequency define the spatial resolution of the resulting surface, and can be optimized depending on the desired survey sensitivity. This paper examines gamma survey as a tool for decommissioning NORM sites and provides an overview of current gamma survey technology designed to improve the efficiency and effectiveness of the decommissioning process. Topics to be discussed in the paper include: - An overview of gamma survey systems, and the utility of different delivery vehicles depending on desired cost, desired spatial resolution, and site topography. - The influence of physical detector characteristics on detection sensitivity and survey planning. - The tradeoff between high-resolution and large spatial extent, but inherently uncertain data, and low-resolution, low spatial extent, but highly certain data, as well as the specific utility of each of these types of data during NORM facility decommissioning. - Confounding variables that may limit the utility of gamma survey at some sites (e.g., radon gas and spatial heterogeneity / hot spots), and methods to plan for and control these conditions. Results show that the confounding variables, such as radon and data output can greatly influence the overall data quality associated with the decommissioning process. In addition, the use of real-time and aerial survey platforms provides a method for ensuring proper spatial extent of the data. When applied thoughtfully, gamma survey is a powerful tool for detecting NORM radionuclides in the environment and a cost-effective technique for identifying areas requiring remediation. However, entities performing or using gamma survey as a decommissioning tool must be aware of both its advantages and its limitations before basing remediation or regulatory action on gamma survey results. (authors)

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

Integrated Molten Salt Reactor Modeling Capabilities in NEAMS Thermal Hydraulics Tools

The DOE neams program supports a full range of computational thermal fluids analysis capabilities and code developments for a broad range of advanced reactor concepts. The research and development approach under the thermal fluids technical area synergistically combines three length and time scales in a hierarchical multi-scale approach. To enable multi-scale thermal fluids capability using these codes, a key joint effort has been underway to develop an integrated system- and engineering-scale thermal fluids analysis capability, through integration of SAM and Pronghorn codes, both based on the MOOSE framework. This report summarizes recent advances in developing an integrated system- and engineering-scale modeling capability for the msr concept, which has gained significant interest in recent years. A consistent framework was established by coupling Pronghorn and SAM through the Saline interface, with thermophysical properties provided by the Molten Salt Thermal Property Database (MSTDB-TP). Further improvements were made to the coupling schemes and domain-overlapping strategies, enhancing the stability and robustness of multi-code simulations. Verification and validation efforts demonstrate the accuracy of this integration across a range of benchmark problems, including one-dimensional heated pipe flows, three-dimensional natural convection loops with evolving isotopic compositions, and \gls{msre} demonstration cases. Within Pronghorn, new capabilities were introduced to model corrosion and noble-metal plating phenomena, supported by an extended thermal-hydraulics framework and refined turbulence treatments. To capture two-phase flow behavior, a multiphase Euler–Euler model was implemented in Pronghorn, including advanced closure relations, high-resolution advection techniques, and capillary force reconstruction. Preliminary verification cases confirm the fidelity of the approach, while planned validation efforts target canonical multiphase benchmarks and application to msr components such as the msre pump bowl. Finally, updates to SAM’s msr mass transfer modeling were extended to consider noble gas migration into porous structures like graphite. The point kinetics model was updated to include reactivity feedback contributions from any defined species, such as xenon. The gas transport model was expanded for applicability to gas mixtures, bubble efflux phenomena, and species transport between liquid and gas phases. A selection of multi-scale Sherwood number correlations from MOSCATO/NekRS and multi-phase correlations from literature have been added for improved accuracy in calculating mass transfer coefficients. A companion effort on developing system-level redox corrosion has also been incorporated into SAM. Collectively, these enhancements strengthen the predictive capability of SAM and Pronghorn for simulating MSR thermal-hydraulics, corrosion, multiphase behavior, and fission-product transport, providing a more complete toolset for design, safety analysis, and licensing support of next-generation \gls{msr}s.

42 - ENGINEERING↗
DOE OSTI2026

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - NLR Historical Wind

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis from variable sources, hydrogen compression and storage, and 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 represents 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 nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production by conducting a statistical analysis of historical wind data over a five-year period (2020-2025) from a single 1.5MW turbine manufactured by General Electric (GE) located at NLR’s Flatirons Campus, to generate an experimental test profile that was deployed on a 1.25-MW proton exchange membrane type MC250 electrolyzer system manufactured by Nel Hydrogen . [1] While the electrolyzer balance-of-plant supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. The historical wind data provided several metrics, however, the analysis particularly focused on the measured power output by the wind turbine. The power output time series of data for each day was categorized by total energy generation and standard deviation, and the day that represented the highest combination of these two metrics was chosen – December 25th, 2022. This process was then repeated for a moving four-hour window within this day to identify the most statistically variable period. Finally, this four-hour period was scaled by 65% to match the 1.25 MW electrolyzer. The electrolysis system controls hydrogen production by varying DC current applied to the stack, from a maximum of 3000 A to a minimum safe operation 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 historical wind profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1 Hz frequency. For more details on the statistical analysis process, see the presentation labeled “ Public Reference Data for Megawatt-Scale Hydrogen Electrolysis” provided with each data entry. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wind turbine electrolysis experiment and is formatted as follows: {technology}_{scaling factor}-{electrolyzer ramp rate in amperes/second} For instance, “wind-GE1.5MW_0.65-400.zip” represents the hour-long experiment using historical data from the wind-GE1.5MW turbine, scaled to 65%, with the electrolyzer power supply set to a maximum ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production, electrolysis power consumption, and wind power input. A PDF file detailing the historical wind data statistical analysis used to generate the wind profile. An experiment labeled “characterization_200.zip” demonstrates the MC250 electrolyzer steady-state response with 30-minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all simulated wind experiments combined into one dataset labeled "combined_historical_wind_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis [1] nelhydrogen.com/product/mc-series-electrolyser .

08 HYDROGEN↗
DOE OSTI2021

Local-Field Effects in Linear Response Properties within a Polarizable Frozen Density Embedding Method

In this work, we present a polarizable frozen density embedding (FDE) method for calculating polarizabilities of coupled subsystems. The method (FDE-pol) combines a FDE method with an explicit polarization model such that the expensive freeze/thaw cycles can be bypassed, and approximate nonadditive kinetic potentials are avoided by enforcing external orthogonality between the subsystems. To describe the polarization of the frozen environment, we introduce a Hirshfeld partition-based density-dependent method for calculating the atomic polarizabilities of atoms in molecules, which alleviates the need to fit the atomic parameters to a specific system of interest or to a larger general set of molecules. Further, we show that the Hirshfeld partition-based method predicts molecular polarizabilities close to the basis set limit, and thus, a single basis set-dependent scaling parameter can be introduced to improve the agreement against the reference polarizability data. To test the model, we characterized the uncoupled and coupled response of small interacting molecular complexes. Here, the coupled response properties include the perturbation of the frozen system due to the external perturbation which is ignored in the uncoupled response. We show that FDE-pol can accurately reproduce both the exact uncoupled polarizability and the coupled polarizabilities of the supermolecular systems. Using damped response theory, we also demonstrate that the coupled frequency-dependent polarizability can be described by including local field effects. The results emphasize the necessity of including local-field effects for describing the response properties of coupled subsystems, as well as the importance of accurate atomic polarizability models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗