SEARCH · Search NASA
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.
The presentation highlights the market analysis capabilities, both price-taker analysis and the novel market surrogate models approach, applied to optimal design of electricity and hydrogen co-production systems. The research findings were presented at the 2023 AIChE Annual Meeting.
The Deep Underground Neutrino Experiment (DUNE) is a next-generation neutrino experiment that will generate an unprecedented volume of heterogeneous information-from documentation and technical notes to experimental data and reconstruction pipelines. Efficient knowledge retrieval and contextual understanding are increasingly critical for collaboration-wide productivity and onboarding. In this work, we present DUNE-GPT, a prototype framework that leverages large language models (LLMs) and retrieval-augmented generation (RAG) to enable natural-language querying of DUNE's internal documentation and technical resources. The system provides an intelligent interface for DUNE collaborators to interact with experiment-specific knowledge while maintaining data privacy and infrastructure compliance within Fermilab computing resources.
Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).
Mining large-scale data to discover biologically relevant information remains a challenge despite the rapid development of bioinformatics tools. Here, we have developed a new tool, PathTracer, to identify biologically relevant information flows by mining genome-wide protein–protein interaction networks following integration of gene expression data. PathTracer successfully mines interactions between genes and traces the most perturbed paths of perceived activities under the conditions of the study. Here, we further demonstrated the utility of this tool by identifying adaptation mechanisms of hypoxia-induced dormancy in Mycobacterium tuberculosis (Mtb).
Abstract This paper provides a summary of the work done within the OC6 Phase II project, which was focused on the implementation and verification of an advanced soil–structure interaction model for offshore wind system design and analysis. The soil–structure interaction model comes from the REDWIN project and uses an elastoplastic, macroelement model with kinematic hardening, which captures the stiffness and damping characteristics of offshore wind foundations more accurately than more traditional and simplified soil–structure interaction modeling approaches. Participants in the OC6 project integrated this macroelement capability to coupled aero‐hydro‐servo‐elastic offshore wind turbine modeling tools and verified the implementation by comparing simulation results across the modeling tools for an example monopile design. The simulation results were also compared to more traditional soil–structure interaction modeling approaches like apparent fixity, coupled springs, and distributed springs models. The macroelement approach resulted in smaller overall loading in the system due to both shifts in the system frequencies and increased energy dissipation. No validation work was performed, but the macroelement approach has shown increased accuracy within the REDWIN project, resulting in decreased uncertainty in the design. For the monopile design investigated here, that implies a less conservative and thus more cost‐effective offshore wind design.
Satellite quantification of aerosol effects on clouds relies on aerosol optical depth (AOD) as a proxy for aerosol concentration or cloud condensation nuclei (CCN). However, the lack of error characterization of satellite-based results hampers their use for the evaluation and improvement of global climate models. We show that the use of AOD for assessing aerosol–cloud interactions (ACIs) is inadequate over vast oceanic areas in the subtropics. Instead, we postulate that a more physical approach that consists of matching vertically resolved aerosol data from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite at the cloud-layer height with Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua cloud retrievals reduces uncertainties in satellite-based ACI estimates. Combined aerosol extinction coefficients (σ) below cloud top (σBC) from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and cloud droplet number concentrations (N d ) from MODIS Aqua yield high correlations across a broad range of σBC values, with σBC quartile correlations ≥0.78. In contrast, CALIOP-based AOD yields correlations with MODIS Nd of 0.54–0.62 for the two lower AOD quartiles. Moreover, σBC explains 41% of the spatial variance in MODIS N d , whereas AOD only explains 17%, primarily caused by the lack of spatial covariability in the eastern Pacific. Compared with σBC, near-surface σ weakly correlates in space with MODIS N d , accounting for a 16% variance. It is concluded that the linear regression calculated from ln(N d )–ln(σBC) (the standard method for quantifying ACIs) is more physically meaningful than that derived from the Nd–AOD pair.
Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.
Electromagnetic waves propagating through vacuum can polarize virtual electron–positron pairs; this polarization, in turn, nonlinearly modifies their propagation. A semi-classical nonlinear wave equation describing the propagation is derived from the Euler–Heisenberg Lagrangian density, which captures vacuum polarization effects up to the one-loop level. In this article, we present a reduced-actionintegral approach that enables rapid modeling of nonlinear phenomena arising from the Euler– Heisenberg Lagrangian. Application of the variational principle to the reduced action provides equations of motion for familiar light-pulse parameters, such as spot size, phase, polarization, and phase-front curvature, without requiring full-field simulations. Three examples demonstrate the utility of the approach: phase modulation, birefringence, and frequency mixing.
Here, we use nonequilibrium steady states to study the effect of dissipation-assisted operator evolution (DAOE) on the scaling behavior of transport in one-dimensional spin chains. We consider three models in the XXZ family: the XXZ model with staggered anisotropy, which is chaotic; XXZ model with no external field and tunable interaction, which is Bethe-ansatz integrable and (in the zero interaction limit) free-fermion integrable; and the disordered XY model, which is free-fermion integrable and Anderson localized. We find evidence that DAOE's effect on transport is controlled by its effect on the system's conserved quantities. To the extent that DAOE preserves those symmetries, it preserves the scaling of the system's transport properties; to the extent it breaks those conserved quantities, it pushes the system towards diffusive scaling of transport.
Traditional building energy modelling workflows remain labor-intensive and error-prone, requiring specialized expertise that limits broader adoption. This paper introduces a novel Model Context Protocol (MCP)-enabled framework that connects AI assistants to EnergyPlus through MCP, a standardized interface for tool invocation and context management. Two complementary integration paradigms are presented and compared: conversational integration, where users interact through natural language while an AI assistant orchestrates MCP tools on demand, and agentic workflow integration, where specialized agents coordinate autonomously to complete multi-step tasks. Using an experimental testbed for residential buildings, the end-to-end workflows are demonstrated. The conversational approach reduced typical inspection and modification tasks from 1-2 h to under 15 min, while maintaining full transparency through visible tool invocations. The agentic approach automated parametric analysis. These demonstrations establish MCP as a foundational layer for AI-assisted building energy modelling, enabling natural language interactions with simulation tools while preserving professional oversight and decision-making authority.
This project advances a unified, open-source framework for the optimal co-design of thermal, electrical, and control systems in grid-interactive efficient districts (GEDs). As communities integrate growing levels of distributed energy resources, traditional approaches that model thermal and electrical networks independently lead to reduced efficiency, limited flexibility, and missed opportunities for coordinated operation. To address these challenges, the research team developed a comprehensive suite of physics-based models, control algorithms, and software tools that enable holistic simulation, optimization, and demonstration of district-scale energy systems.
Residential buildings’ heating and cooling loads are associated with significant carbon emissions and peak electricity demand. Phase change material (PCM) based thermal energy storage (TES) can be used for space heating and cooling by embedding into the heat pump equipment. Past work on TES integrated with heat pumps (HP) has demonstrated significant load shifting and economic benefits. However, the potential for TES to reduce carbon emissions has not been widely explored. This study evaluates carbon mitigation potential of an ice-based TES coupled to HP (HP-TES) based on a simple rule-based control strategy accounting for electric grid emissions data. A vapor compression HP model using Engineering Equation Solver (EES). The modeled HP-TES system for single-family residential building demonstrates decreased carbon emissions with reduced peak utility cost.
Isotropic integrable spin chains such as the Heisenberg model feature superdiffusive spin transport belonging to an as-yet-unidentified dynamical universality class closely related to that of Kardar, Parisi, and Zhang (KPZ). To determine whether these results extend to more generic one-dimensional models, particularly those realizable in quantum simulators, we investigate spin and energy transport in non-integrable, long-range Heisenberg models using state-of-the-art tensor network methods. Despite the lack of integrability and the asymptotic expectation of diffusion, for power-law models (with exponent 2<α<∞) we observe long-lived z=3/2 superdiffusive spin transport and two-point correlators consistent with KPZ scaling functions, up to times t∼103/J. We conjecture that this KPZ-like transport is due to the proximity of such power-law-interacting models to the integrable family of Inozemtsev models, which we show to also exhibit KPZ-like spin transport across all interaction ranges. Finally, we consider anisotropic spin models naturally realized in Rydberg atom arrays and ultracold polar molecules, demonstrating that a wide range of long-lived, non-diffusive transport can be observed in experimental settings.
A paramount goal in nuclear physics is to unify ab initio treatments of bound and unbound states. The position-space quantum Monte Carlo (QMC) methods have a long history of successful bound-state calculations in light systems but have seen minimal implementation in unbound systems. Here, we introduce a numerical method to improve the efficiency and accuracy of unbound-state calculations in QMC. As an initial application, we compute scattering observables for the smallest system available to probe three-body forces, the neutron-triton system, using variational Monte Carlo (VMC) wave functions. The method involves inferring long-range amplitudes in the wave function from integrals over the short-range region where all the particles interact. This approach using integral relations is well established in the literature; here, we develop it for the QMC framework. We validate our approach with a consistency check between short-range spectroscopic overlap functions computed from direct evaluation and from the integral relations; scattering amplitudes are long-range asymptotics of those overlaps. Comparison against published benchmark calculations using the same potential demonstrates that when applied to the current VMC wave functions, the integral method produces more accurate scattering observables than direct evaluation from the same variational wave function. However, it still differs noticeably from the exact results. Using additional interactions, we then present phase shifts and mixing parameters for the n + 3 H system. In particular, we present one of the first applications of the Norfolk family of local coordinate-space chiral potentials in unbound systems of A > 2. The Norfolk results accurately describe s-wave scattering but predict p-wave cross sections too large. Compared with previous QMC scattering calculations, the integral method avoids difficulties associated with the precise computation of energy differences and with convergence outside the interaction region, which is particularly severe in the variational calculation. Application of the integral method here paves the way for its use in Green's function Monte Carlo (GFMC) calculations. In GFMC, the wave functions are more accurate, but the high-precision convergence of their tails is slow, and there are additional difficulties in reading out amplitudes. Here, the integral methods will address both of those remaining problems.
Nanomaterials offer exceptional properties for photonic, electronic, quantum, magnetic, structural, mechanical, chemical, nuclear, and biological functionality. However, real access to enhanced functionality remains limited without connecting the nanoscale across the mesoscale to macroscale assemblies and the ability to define and control nanomaterials organization, interactions, and interfaces. Integration is essential for accessing and controlling functionality to harness nanomaterials properties while also generating new behaviors and properties. Integration thus provides a route to discovering, generating, and using intrinsic and emergent nanomaterials properties, which define the two focus areas of this leadership theme.
Surface flow and subsurface flow constitute a naturally linked hydrologic continuum that has not traditionally been simulated in an integrated fashion. Recognizing the interactions between these systems has encouraged the development of integrated hydrologic models (IHMs) capable of treating surface and subsurface systems as a single integrated resource. IHMs are dynamically evolving with improvements in technology, and the extent of their current capabilities are often only known to the developers and not general users. This article provides an overview of the core functionality, capability, applications, and ongoing development of one open-source IHM, ParFlow. ParFlow is a parallel, integrated, hydrologic model that simulates surface and subsurface flows. ParFlow solves the Richards equation for three-dimensional variably saturated groundwater flow and the two-dimensional kinematic wave approximation of the shallow water equations for overland flow. The model employs a conservative centered finite-difference scheme and a conservative finite-volume method for subsurface flow and transport, respectively. ParFlow uses multigrid-preconditioned Krylov and Newton–Krylov methods to solve the linear and nonlinear systems within each time step of the flow simulations. The code has demonstrated very efficient parallel solution capabilities. ParFlow has been coupled to geochemical reaction, land surface (e.g., the Common Land Model), and atmospheric models to study the interactions among the subsurface, land surface, and atmosphere systems across different spatial scales. This overview focuses on the current capabilities of the code, the core simulation engine, and the primary couplings of the subsurface model to other codes, taking a high-level perspective.
MLMOD is a software package for incorporating machine learning approaches and models into simulations of microscale mechanics and molecular dynamics in LAMMPS. Recent machine learning approaches provide promising data-driven approaches for learning representations for system behaviors from experimental data and high fidelity simulations. The package facilitates learning and using data-driven models for (i) dynamics of the system at larger spatial-temporal scales (ii) interactions between system components, (iii) features yielding coarser degrees of freedom, and (iv) features for new quantities of interest characterizing system behaviors. MLMOD provides hooks in LAMMPS for (i) modeling dynamics and time-step integration, (ii) modeling interactions, and (iii) computing quantities of interest characterizing system states. The package allows for use of machine learning methods with general model classes including Neural Networks, Gaussian Process Regression, Kernel Models, and other approaches. Here we discuss our prototype C++/Python package, aims, and example usage. For related papers, examples, updates, and additional information see https://github.com/atzberg/mlmod and http://atzberger.org/.
Developing an atomistic understanding of ionizing radiation induced changes to organic materials is necessary for intentional design of greener and more sustainable materials for radiation shielding and detection. Cocrystals are promising for these purposes, but a detailed understanding of how the specific intermolecular interactions within the lattice upon exposure to radiation affect the structural stability of the organic crystalline material is unknown. This study evaluates atomistic-level effects of γ radiation on both single- and multicomponent organic crystalline materials and how specific noncovalent interactions and packing within the crystalline lattice enhance structural stability. Dose studies were performed on all crystalline systems and evaluated via experimental and computational methods. Changes in crystallinity were evaluated by p-XRD and free radical formation was analyzed via EPR spectroscopy. Type of intermolecular interactions and packing within the crystal lattice was delineated and related to the specific free radical species formed and the structural integrity of each material. Periodic DFT and HOMO-LUMO surface mapping calculations provided atomistic-level identifications of the most probable sites for the radicals formed upon exposure to γ radiation and relate intermolecular interactions and molecular packing within the crystalline lattice to experimental results.