Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Building performance simulation”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Human Factors and Technologies Design to Improve User Acceptance of Pooled Rideshare for Increasing Transportation System Energy Efficiency

This multi-year project delivered a comprehensive, human-factors-driven framework to understand, model, and improve pooled rideshare (PR) adoption in the United States. Through three large-scale national survey studies involving more than 16,000 participants across multiple cities and demographic groups, the research established one of the most extensive datasets to date on user perceptions, behavioral barriers, and service expectations related to pooled rideshare. These data revealed key human factors barriers of user acceptance of PR and suggested potential actionable experience optimizations that could lead to increased PR usage. This foundational knowledge guided the development of novel human-factors models and behavioral choice models that quantify how psychological, demographic, and trip-level factors influence willingness to pool. Building on these empirical insights, the project developed advanced behavioral modeling tools, including mixed logit and integrated choice and latent variable models, to capture both observable and latent influences on PR adoption. These models significantly improved the ability to predict riders’ acceptance of pooled trips, explaining choice heterogeneity through latent constructs such as safety, service experience, privacy concerns, time sensitivity, and environmental attitudes. Together, these models provide a robust analytical foundation for designing PR systems that more effectively meet user needs. The project translated human-factors insights and behavioral models into actionable technology innovations by extending POLARIS—an agent-based, activity-based travel simulation platform—into a fully functional pooled rideshare simulation environment. New PR modules, acceptance models, and regional scenarios were implemented for Greenville, SC and Austin, TX, enabling high-fidelity validation of algorithmic strategies under realistic demand and traffic conditions. The simulation platform supported the development and evaluation of adaptive discount-based assignment algorithms, enhanced willingness-to-pay formulations, demographic-aware incentive mechanisms, and a proactive joint assignment and repositioning strategy. Simulation results demonstrated substantial gains in pooling uptake, average vehicle occupancy, energy efficiency, and fleet profitability. In Greenville, pooling adoption more than doubled, while reductions in vehicle-miles traveled and energy consumption were significant. In Austin, pooling improvements were achieved with minimal service-quality trade-offs, and profitability increased across all fleet sizes. Through this research, we developed a comprehensive understanding of the human factors barriers that limit user acceptance of pooled rideshare services. These insights enabled the design of human-factors-aware pooled rideshare technologies that more effectively address user concerns and improve adoption rates. By integrating these models into an advanced agent-based simulation framework, we demonstrated that higher adoption of pooled rideshare can lead to measurable improvements in energy efficiency and system performance. Together, these contributions establish a validated pathway from human-centered analysis to technology development and energy-saving outcomes, supporting national goals for more sustainable and efficient mobility systems.

Jia, Yunyi↗
DOE OSTI2022

Development and Demonstration of a Fuel-Efficient, Class 8 Tractor & Trailer Engine System (SuperTruck II)

Navistar presents the SuperTruck II (ST II) Final Report to the Unites States Department of Energy (US DOE), which covers the five Budget Periods (BPs) from 10-1-2016 through 6-30-2022. For ST II, Navistar built on the achievements of the SuperTruck I (ST I) Program as a catalyst to continue critical research, design and development, testing, and operations to reach the ambitious goals of the ST II project. This approach allowed Navistar to continue contributing to the essential needs of our nation for safe, efficient, and cost-effective delivery of goods and services, as we reduced negative environmental effects and improved operational productivity. This document contains information specified in DOE F 4600.2, Final Scientific/Technical Report DOE F 241.3, B. SCIENTIFIC/TECHNICAL REPORTS, explaining how we met and exceeded program requirements. Throughout this Final Report, Navistar extracted information from documents prepared during the project that represent our management, design and development, building, and testing efforts to meet and exceed SuperTruck II project goals. Navistar followed Plan requirements to achieve / exceed Project Objectives: a) >100% improvement in vehicle freight efficiency (FE) (on ton-MPG basis) relative to 2009 baseline with stretch goal of 140% improvement [actual: 170%); b) >55% engine brake thermal efficiency (BTE) demonstrated in operational engine at a 65-mph cruise point on a dynamometer – ≥31% increase from 2009 baseline [actual: 55.20% of combined BTE) ; and c) development and implementation of commercially cost effective technologies (in terms of a simple payback). Technology selection / development path focused on developing technologies applicable for production within 3-year approach, while ensuring technology readiness and cost of ownership for end users. The Program was organized into five budget periods: Requirements / Technology Assessment and Initial Hardware Testing; Technology Development and Concept Readiness Demonstration; Technology Finalization and Validation Tractor / Trailer Fabrication, Integration and Commissioning Demonstration; and Fuel Economy (FE) and Brake Thermal Efficiency (BTE) and Program Completion. Leadership was provided by DOE, with tasks performed by laboratories (Argonne National Laboratory, Lawrence Livermore National Laboratory); partners at Bosch, TPI, Dana, and J.B. Hunt; , and support from University of Michigan and Clemson University. Navistar lead this team with Principal Investigator / Contracting Officer; Project Manager (PM); Vehicle, Engine, and Aftertreatment Engineers; Finance Manager, Technical Program Leads, and Legal/IP; and other key personnel. Work also included personnel in risk management; funding / budget / finance. Work involved analysis, development, testing, and down selection of individual/system engine, aftertreatment, and vehicle technologies, with integration of selected technologies into a prototype vehicle for demonstration of fuel-efficiency gain. Work also included component/integrated system level development of truck and trailer aerodynamics, base engine efficiency, advanced aftertreatment, combustion efficiency, waste heat recovery, hybrid powertrain, reduced rolling resistance, weight reduction, idle reduction, and driver feedback. As ST II progressed, Navistar performed computer-based modeling / simulations of technologies focused on the primary operational areas: Engine, Aftertreatment, and Vehicle. During the ST II Program, the COVID Virus outbreak unexpectedly challenged by the effects of, which affected staffing, scheduling, design, supplies, availability of materials, production procedures, and testing. The DOE responded by extending the program by three quarters to ensure that project tasks were completed for this vital project. Focus continued on analyzing, developing, testing, and down selecting individual-/system-level engine and vehicle technologies for integration of the final selected technologies into a prototype vehicle that would demonstrate fuel-efficiency gains made possible through these technologies. This included component/integrated system-level development of truck and trailer aerodynamics, base engine efficiency, advanced aftertreatment, combustion efficiency, waste heat recovery, solar power, distributed and intelligent vehicle power, hybrid powertrain, reduced rolling resistance, weight reduction, idle reduction, and driver feedback. Throughout the program, function, reliability, and performance at all levels were ensured through testing. Proof of this approach was demonstrated in multiple, on-road demonstrations: Scenario A (Flatland) Fuel Economy, Scenario B (Hilly) Fuel Economy, and City Cycle Tests. Other benefits derived from ST II included new/improved products, publications, patents, and next-step capabilities related to electric/hydrogen vehicles and autonomous driving.

Zukouski, Russ↗
DOE OSTI2020

Data analysis and modeling pipelines for controlled networked social science experiments

There is large interest in networked social science experiments for understanding human behavior at-scale. Significant effort is required to perform data analytics on experimental outputs and for computational modeling of custom experiments. Moreover, experiments and modeling are often performed in a cycle, enabling iterative experimental refinement and data modeling to uncover interesting insights and to generate/refute hypotheses about social behaviors. The current practice for social analysts is to develop tailor-made computer programs and analytical scripts for experiments and modeling. This often leads to inefficiencies and duplication of effort. In this work, we propose a pipeline framework to take a significant step towards overcoming these challenges. Our contribution is to describe the design and implementation of a software system to automate many of the steps involved in analyzing social science experimental data, building models to capture the behavior of human subjects, and providing data to test hypotheses. The proposed pipeline framework consists of formal models, formal algorithms, and theoretical models as the basis for the design and implementation. We propose a formal data model, such that if an experiment can be described in terms of this model, then our pipeline software can be used to analyze data efficiently. The merits of the proposed pipeline framework is elaborated by several case studies of networked social science experiments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Using Dense Array Waveform Correlations to Build a Velocity Model with Stochastic Variability

In an ongoing effort to improve 3D seismic-wave propagation modeling for frequencies up to 10 Hz, we used cross correlations between vertical-component waveforms from an underground chemical explosion to estimate the statistical properties of small-scale velocity heterogeneities. The waveforms were recorded by a dense 2D seismic array deployed during the Source Physics Experiments for event number 5 (SPE-5) in a series of six underground chemical explosions, conducted at the Nevada National Security Site. The array consisted of 996 geophones with a 50–100 m grid spacing, deployed at the SPE site at the north end of the Yucca Flat basin. The SPE were conducted to investigate the generation and propagation of seismic and acoustic waves from underground explosions. Comparisons of decay rates of waveform cross correlations as function of interstation distance, computed for observed and synthetic seismograms from the SPE-5 chemical explosion, were used to constrain statistical properties of correlated stochastic velocity perturbations representing small-scale heterogeneities added to a geology-based velocity model of the Yucca Flat basin. Using comparisons between recorded and simulated waveform cross correlations, we were able to recover sets of statistical properties of small-scale velocity perturbations in the velocity model that produce the best-fit between the recorded and simulated ground motion. The stochastic velocity fluctuations in the velocity model that produced the smallest misfits have a horizontal correlation distance of between 400 and 800 m, a vertical correlation distance between 100 and 200 m, and a standard deviation of 10% from the nominal model velocity in the alluvium basin layers. They also have a horizontal correlation distance of 1000 m, a vertical correlation distance of 250 m, and a standard deviation of 6% in the underlying and consolidated sedimentary layers, up to a depth of 4 km. Comparisons between observed and simulated wavefields were used to assess the proposed small-scale heterogeneity enhancements to the Yucca Flat basin model. We found that adding a depth-resolved stochastic variability to the geology-based velocity model improves the overall performance of ground-motion simulations of an SPE-5 explosion in the modeled frequency range up to 10 Hz. The results may be applicable to other similar basins.

58 GEOSCIENCES↗
DOE OSTI2023

Adroit Materials Final Scientific/Technical Report kV-class GaN-based Junction Barrier Schottky diodes using ion implantation

The primary aim of this research was to develop GaN-based Junction Barrier Schottky (JBS) diodes using an innovative ion implantation process previously established in ARPA-E funded projects. The central focus of our proposed technology revolves around selective area p-type doping, accomplished through the implantation of Mg ions. This approach builds upon our successes in the ARPA-E PNDIODES program, advancing towards commercial device integration. Selective area p-doping plays a pivotal role in realizing the next generation of GaN-based power devices, capable of significantly reducing the carbon footprint in the United States by several million tons. While ion implantation is a well-established technique for achieving selective area doping in SiC and Si materials, its feasibility in GaN had not been demonstrated until now. To fabricate high voltage GaN JBS diodes, we initially created thick n-type drift layers with high carrier concentrations ranging from 5×1015 cm-3 to 2×1016 cm-3 and very high mobilities. Subsequently, Mg ions were selectively implanted to form p-type islands within the n-type drift layer. To reduce electric field crowding at the edge of the diode and to achieve high breakdown voltage, junction edge termination (JTE) and floating field rings (FFRs) were formed using Mg implantation. A high temperature, high-pressure post-implantation annealing process was carried out to activate the implanted Mg ions. As a result, we were able to demonstrate GaN JBS diodes with a breakdown voltage of 915 V and an on-resistance of 0.6 mΩ·cm2. These diodes exhibited a forward bias current density of 1 kA/cm2 at 1.5 V. Subsequently, we achieved GaN JBS diodes with a remarkable breakdown voltage of 1900 V and an on-resistance of 1.9 mΩ·cm2, capable of sustaining a forward bias current density of 0.5 kA/cm2 at 1.5 V. Importantly, the ON and OFF state performance of these GaN JBS diodes surpassed that of Si and SiC-based power diodes reported in existing literature. Lastly, we successfully grew 60 μm thick GaN:Si layers using HVPE with a carrier concentration of approximately 3 to 5×1015 cm-3. Based on simulation and empirical data these devices represent 5 kV GaN JBS power diodes, leveraging the developed processes in this project.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of Mineral Impurities and Heat on Uranium(VI) Sorption onto Bentonite (Final Report)

Most nuclear waste disposal options include compacted bentonite, consisting of montmorillonite clay, as part of a barrier system to minimize contaminant mobility. Uranium (U) is the primary element in spent nuclear fuel, and from an environmental perspective, a potential contaminant of water resources. Furthermore, uranium may control nuclear fuel degradation rates and the consequent release of other radioactive contaminants based on its diffusive transport away from waste canisters. Uranium sorption onto clay and other mineral surfaces is expected to limit U(VI) mobility in these systems. However, at this point a prediction of U(VI) sorption and transport behavior in performance assessment (PA) models is complicated by a series of factors, such as: (1) the presence of various U(VI) solution species with different charges and sorption characteristics; (2) the complex microstructure of montmorillonite clay resulting in two types of clay porosities and multiple sorption sites; (3) the largely unknown effects of bentonite mineral impurities on pore water chemistry, U(VI) sorption and diffusion behavior; and (4) the potential impacts of heat, generated by the decay of spent fuel, on mineralogical and microstructural transformations, and any subsequent effects on radionuclide sorption and mobility. In this project, we investigated the effects of calcite impurities on U(VI) sorption and diffusion onto montmorillonite before and after mineral exposure to heat. Overall, our main goals were to: 1. Gain a mechanistic understanding of the effects of calcite impurities on U(VI) sorption and diffusion in montmorillonite/bentonite at variable chemical conditions, and before and after mineral exposure to heat; and to 2. Work towards the development of new surface complexation and diffusion models to predict U(VI) sorption and mobility in montmorillonite/bentonite while accounting for the impacts of calcite impurities and heat. In order to achieve these goals, we pursued the following research objectives and tasks: 1. Characterize the effects of calcite impurities on U(VI) sorption onto Na-montmorillonite, and evaluate if and how these effects change due to solids’ exposure to heat (Task 1). 2. Characterize U(VI) surface speciation in solid samples from Task 1, with a particular focus on U(VI) precipitates and co-precipitates (Task 2). 3. Test the relevance of observed changes in U(VI) sorption characteristics for U(VI) diffusive transport behavior in bentonite barriers (Task 3). 4. Evaluate the exclusion of anionic U(VI) solution complexes from clay nanopores (Task 4). 5. Building on Tasks 1-4, start the development of new surface complexation and reactive transport models that can simulate and predict U(VI) sorption and diffusion in complex bentonite systems (Task 5). A detailed description of all tasks, and their specific results and conclusions is provided in this report.

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

Tailored Fiber Placement for Complex Preforms

Tailored Fiber Placement (TFP) offers a novel approach to optimize fiber architecture for the fabrication of complex, structural parts not traditionally suitable for advanced composites. This technology not only offers new routes for weight reduction via metal substitution, it also offers cost reduction through minimization of material scrap and reduced labor. This reduction in component weight leads to increased fuel efficiency, and reduced production energy consumption, thereby, helping to achieve the stated IACMI technical goals. This technology leverages centuries of manufacturing development in support of the textile and embroidery industry. One major drawback to this technology is the lack of commercial or non- proprietary structural performance data and robust analytical tools used to optimize fiber architecture and predict performance. This project was structured to utilize common sub-element features to validate analytical performance tools, generate performance data, and gather cost and performance data on components of interest. This project was designed to give industry sponsors the confidence and ability to take full advantage of TFP to fabricate primary, highly loaded structure and integrate features such as metallic fasteners. The project focused principally on the use of high strength carbon fiber, such as T700, and the use of aerospace epoxy resin matrix to primarily support development of new composite applications in vehicle, aerospace, and industrial markets. This project applied previously developed analytical tools to predict the performance of TFP produced parts. This work focused on developing the pipeline to characterize material in order to accurately predict component performance when modifying the TFP print paths and stitch density. This focused on experimental characterization via standardized ASTM testing, alongside experimental testing of more representative service components by testing curved beam strength, beam shear performance, a large scale TFP lug, and ultimately designing a fully TFP clip bracket that reduced weight and cost compared to a traditional metallic component. The new knowledge gained from this program included: 1) development and demonstration of novel analytical tools applied to analysis of TFP preforms; 2) development and demonstration of a building block approach using coupons and sub-elements to optimize the design of a more complex component; 3) demonstration that optimized fiber orientation using TFP can exceed performance of conventional textile composite materials and can open new applications currently limited to metallic components; 4) Demonstration of performance and cost benefits of the TFP process as compared to metallic and conventional textile composites. Recommendations for follow-on work include development of design allowables to assess the impact of high temperature/moisture exposure or saturation during loading, tracking the impact of stitching needle wear on the performance of parts and ability to stitch thicker preforms, using TFP preforms as local reinforcement at areas of bearing or complex loading, and topology optimization of components by tow steering. The expertise developed during the course of this project can be leveraged to provide commercial engineering design and fabrication services using TFP. UDRI is in the process of formalizing their partnership with Spintech, who will serve as the commercialization partner for this technology and provide molding services and deliver finished components to the end user. UDRI will continue to produce the preforms until the economics allow Spintech to procure its own TFP equipment or lease UDRI equipment, at which point UDRI will step away from manufacture and serve as the engineering and design lead on product development.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Fission with Exotic Nuclei (Abbreviated Report)

Nuclear fission is a key mechanism involved in the synthesis of heavy elements in the Cosmos and is the primary explanation for the stability of superheavy elements. Nevertheless, our knowledge of fission remains extremely fragmented. Most experiments have been conducted only on a tiny number of stable actinide nuclei and are often incomplete, leading to gaps in our basic understanding of the process. For many radioactive isotopes, basic fission data such as the charge or mass distribution of the fragments is unknown. These gaps cannot always be filled by simulation alone. Common fission models contain too many free parameters and lack predictive power. In contrast, the fundamental theory of fission under development at LLNL is much more predictive, but its current computational cost is too high to be used extensively for data evaluations. A unique window of opportunity to resolve these limitations has recently opened: the U.S. nuclear science community is ramping up major experimental programs at the Facility for Rare Isotope Beams (FRIB, the DOE flagship facility in low-energy nuclear science), and techniques from machine learning have shown great potential to simplify the use of a fundamental, quantum-mechanical theory of fission. This project has two components. On the experimental side, we acquired and deployed at the HIGS facility a new dual Frisch-Grid ionization chamber to measure correlated fragment-mass, kinetic energy, and angular distributions of fission fragments from induced fission. This new device was used to perform measurements of charge, mass and total kinetic energy of fission fragments in the photofission of 238 U and eight gamma-ray beam energies between 6.2 and 13 MeV, which allowed extracting high-precision independent yields for this reaction. The device was also used to perform measurements of the same quantities in the neutron-induced fission of 234 U with monoenergetic beams of energy between 5 and 8 MeV. In parallel, we collaborated with a team at Commissariat à l’énergie atomique et aux énergies alternatives (CEA) to perform a series of measurements of fission yields in inverse kinematics for the two isotopes of 236 U and 240 Pu. The experiment took place at the Grand Accélérateur National d’Ions Lourds in France in June 2023. The deployment of the VAMOS spectrometer with a new array called PISTA allowed determining the excitation energy of the fissioning system within 1 Mega-electronvolts. The second component of the project involved using deep neural networks to build fast and reliable emulators of our current fission models. In an invited paper published in Frontier in Physics, we showed that autoencoders could successfully compress nuclear wavefunctions in nuclear density functional theory. We achieved a dimensionality reduction of the order of two orders of magnitude while keeping the error in the total energy to less than 0.01%. In a second paper submitted to Physical Review Letters in June 2023 with our collaborators at CEA, we showed that variational autoencoders can learn the collective degrees of freedom driving the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Transportation Hub Infrastructure Expansion: Decision Support Under Uncertainty

The Athena project (www.athena-mobility.org) has worked to investigate the relationship between the Dallas-Fort Worth Airport (DFW) and the greater Dallas area in order to better understand and therefore better inform future decision-making regarding the critical infrastructure that influence mobility between the airport and the city. Through this work, infrastructure related to curbside pickup and drop-off, parking, public transit, and the road network congestion were identified as critical to the operation of the DFW transportation hub. The infrastructure analysis and expansion aspect of the Athena project is focused on the restructuring of the CTA curb as a hierarchical curb and the building or repurposing of parking infrastructure as the interplay between these two areas. Many sources of uncertainty exist that may impact future airport and transportation hub operations, such as passenger volume growth, population demographic changes over time, electric vehicle (EV) adoption rates, and autonomous vehicle (AV) adoption rates. Due to these sources of uncertainty, we have selected for our research a modeling framework that can capture various types of uncertainty and hedge against those uncertainties in the optimization process. We analyze road network and curb congestion, the rise of transportation networking companies, trends in parking usage, existing policies around this infrastructure, airport revenue streams, and other contributing factors to enable infrastructure decision making with less uncertainty. To accomplish this wholistic analysis, we have developed a novel multi-stage, multi-period stochastic optimization model which considers the airport's decisions from 2025-2045 under different possible future macro trajectories and day-to-day variations in operational conditions captured as "annual representation of operations" scenarios with respective probabilities. This model has also been designed to leverage the outputs of various efforts under the Athena project to create a combined decision framework for infrastructure decisions. These various efforts include the route optimization model, the ASPIRES simulation, the mode choice model, and the SUMO traffic simulation. Our computational experiments of this system at scale have resulted in a working version of our infrastructure model which enables the explicit representation and consideration of various sources of uncertainty in the decision process to enable robust, flexible decision-making. This model has been effectively run on NREL's HPC system, Eagle, with large numbers of stochastic scenarios and shows promise as a scalable tool for robust consideration of uncertainties in airport planning. We have tested our model using 30,240 operational circumstances in total, resulting in a problem with more 200 million variables. This model was solved in several different configurations, and a workflow to simulate the performance of the infrastructure model results was developed and deployed. In general, our results indicate that a combination of remote parking, remote curb infrastructure, and dynamic pricing can generate revenue, reduce emissions, accommodate emerging technologies such as AVs and EVs, and manage airport passenger growth over time. We note the success of the proposed strategy depends on the data collection and forecasting abilities of DFW. We have also seen that the AV adoption by TNCs might necessitate larger amounts of remote curb. The results of this work inform strategies for airport infrastructure decision making, as well as demonstrate the value of an adaptable model, but also indicate that there are avenues remaining where further research would be of value.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Rossi-alpha Analysis of Thermal/Epithermal eXperiments Optimized for Polyethylene Thermal Neutron Scattering

Accurate nuclear data is the foundation for predictive simulations and design of new experimental in the nuclear community. The criticality safety community has particular interest in benchmarking assemblies with thermal neutron data. To address such needs the Nuclear Criticality Safety Program (NCSP) funded the Thermal/Epithermal eXperiments (TEX) campaigns that were to be completed between Lawrence Livermore National Laboratory (LLNL) and Los Alamos National Laboratory (LANL). Specifically, analysis of plutonium with a thermal neutron spectrum expands on previous work to help validate thermal scattering law data, which can have larger impacts in thermal applications. The first of these experiments were successfully conducted in 2018, but this paper will focus on the 2021 measurement focused on investigating the the thermal scattering law (TSL) for polyethylene. These experiments were successfully conducted at the National Criticality Experiments Research Center (NCERC) using the Planet vertical lift critical assembly machine. Polethylene plates were layered with trays of Zero Power Physics Reactor (ZPPR) 24 plates in a 12" by 12" square. The ZPPR plates were 2" by 3" by 0.125" bearing weapons grade plutonium. The polyethylene moderator was either 2" or 1.6875" thick. This work builds on the Rossi-alpha calculations done by McKenzie et al. for the same detector-assembly system and will only focus on the Rossi-alpha neutron noise method. This work will aim to further validate the results of the experiment through a novel neutron noise python package. Following similar methodology to the previous analysis, analysis of TEX evaluated the prompt neutron decay constant at delayed critical, $α_{DC}$ using Rossi-alpha for different polyethylene moderator thicknesses. These results will help improve understanding of TSL in critical experiments. Alpha (α), is the prompt neutron decay constant of the measured system and allows for the evaluation of a systems propensity to sustain fission chains via prompt neutrons. The single value description of the assemblies allow for comparison between experiments regardless of composition, geometry, and reflectors/moderators. Rossi-alpha measurements were performed on the polyethylene moderated TEX experiments to estimate $α_{DC}$.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Dynamic Earth Energy Storage: Terawatt-year, Grid-scale Energy Storage Using Planet Earth as a Thermal Battery (GeoTES): Phase I Project (Final Report)

Grid-scale energy storage has been identified by the U.S. Department of Energy’s (DOE) Energy Storage Grand Challenge as a necessary technology to support the continued build-out of intermittent renewable energy resources required to attain a carbon-free energy future. To meet this goal, the 2018 Department of Energy Research and Innovation Act mandated the creation of a comprehensive program to accelerate the development and commercialization of next-generation energy storage technologies. One of numerous energy storage technology options is the storage of excess energy as heated geothermal brine in suitable geologic formations. This concept, known as reservoir thermal energy storage (RTES), geologic thermal energy storage (GeoTES), aquifer thermal energy storage (ATES), etc., relies on the storage of thermal energy in geologic formations for recovery and use in large-scale direct use geothermal (e.g., district heating, industrial processes, etc.) and electrical power generation applications. This thermal energy is derived from excess or waste heat from any high-temperature heat source, such as concentrated solar or from conventional thermal/nuclear generation. As such, RTES can potentially play a significant role in meeting the energy storage shortfall in the coming decades. RTES can provide energy arbitrage through both the storage and production of thermal energy stored in geologic formations for direct use applications and can serve as a source of hot fluids that can be used to generate electricity to support peak demand ramping, thus easing stress on transmission and distribution. This energy storage option has geographic benefits in that energy can be stored locally or regionally depending on the various needs/loads. RTES can also be located across an enormous geographic area, without the need for a traditional hydrothermal resource but where thermal gradients and hydrogeology allow economic exploitation of subsurface heat. The work conducted for this project includes (1) a review of lessons learned from past high-temperature RTES international projects; (2) geochemical experimental investigation and numerical simulations of potential domestic sedimentary reservoirs and (3) development of a thermo-hydrological-mechanical (THM) numerical simulation tool for optimizing formation properties and design parameters to maximize thermal energy storage performance.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Differentiable Multiphysics Codes: A Breakthrough Technology for Simulation and Computing

This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Development, Monitoring, and Control of Fracture Thermal Energy Storage in Crystalline Rock Formations (DEMO-FTES) [Abstract]

Approximately half of global energy consumption is used for heating and cooling. Because fossil fuels are used to meet most of this demand, heating and cooling of buildings produces a large portion of global greenhouse gas (GHG) emissions. The proposed project seeks to demonstrate the feasibility of a seasonal thermal energy storage concept called fracture thermal energy storage (FTES), which has the potential to dramatically lower energy requirements for heating and cooling and improve the resilience of building energy systems. FTES is a technique for building a highly efficient heat exchanger by creating a carefully designed set of fractures in the ground below a building. This heat exchanger allows thermal energy to be stored over seasonal timeframes—for example, hot thermal energy that is easily captured in the summer or energy from waste heat sources—to meet heating needs during colder months. FTES offers a solution with a relatively small footprint and lower cost compared to currently operating aquifer thermal energy storage (ATES) systems and could be used in the many areas without aquifers suitable for ATES. If commercialized, FTES systems could exceed the 2.5 TWh of energy storage per year from the more than 2,800 ATES systems in operation worldwide, which range from 0.1 MW to 30 MW. The savings in CO2 emissions per year are also expected to match or exceed individual ATES projects, with the largest, a 30 MW system used to provide heating and cooling to the University of Technology in Eindhoven in the Netherlands, estimated to save 13,000 tons of CO2 emissions per year. FTES utilizes mature drilling and fracturing technology and therefore has the potential to be rapidly commercialized once demonstrated. The ability to construct and establish flow through an FTES heat exchanger has been demonstrated by a previous project and the potential for efficient, large-capacity energy storage has been shown using numerical models. However, no experimental validation of these numerical estimates of thermal energy storage has been made. The first crucial need to advance FTES technology is to identify the sensitivity of key metrics such as thermal energy storage and production rates, capacities, and efficiencies to design parameters such as the number of fractures, depth/temperature of fractures, size of fractures, and circulation rates. The second crucial need is experimental testing of achievable thermal performance with optimized system design parameters. The proposed scope of work seeks to systematically address these two critical needs through a highly complementary international collaboration spanning theory, laboratory, and mesoscale field evaluation. The proposed work plan calls for using dimensional analysis and existing state-of-the-art numerical simulators to design carefully scaled laboratory and 10-meter-scale field tests of the thermal efficiency of FTES. The existing advanced laboratory and intermediate-scale field testbeds that will be used for this project will allow for detailed monitoring of the system performance during the test and of how the performance changes across time and length scales. These results will determine the feasibility of full- scale FTES systems. If the thermal performance is consistent with model predictions, the results will provide a strong economic justification for rapid commercialization of FTES technology in a wide range of geographical areas.

25 ENERGY STORAGE↗
DOE OSTI2023

Adaptive Ensemble Refinement of Protein Structures in High Resolution Electron Microscopy Density Maps with Radical Augmented Molecular Dynamics Flexible Fitting

Recent advances in cryo-electron microscopy (cryo-EM) have enabled modeling macromolecular complexes that are essential components of the cellular machinery. The density maps derived from cryo-EM experiments are often integrated with manual, knowledge or artificial intelligence driven, and physics-guided computational methods to build, fit, and refine molecular structures. Going beyond a single stationary- structure determination scheme, it is becoming more common to interpret the experimental data with an ensemble of models, which contributes to an average observation. Hence, there is a need to decide on the quality of an ensemble of protein structures on-the-fly, while refining them against the density maps. Here, we introduce such an adaptive decision making scheme during the molecular dynamics flexible fitting (MDFF) of biomolecules. Using RADICAL-Cybertools, and the new RADICAL augmented MDFF implementation (R-MDFF) is examined in high-performance computing environments for refinement of two protein systems, Adenylate Kinase and Carbon Monoxide Dehydrogenase. For the test cases, use of multiple replicas in flexible fitting with adaptive decision making in R-MDFF improves the overall correlation to the density by 40% relative to the refinements of the brute-force MDFF. The improvements are particularly significant at high, 2 - 3 Å, map resolutions. More importantly, the ensemble model captures key features of biologically relevant molecular dynamics that is inaccessible to a single-model interpretation. Finally, the pipeline is applicable to systems of growing sizes, which is demonstrated using ensemble refinement of capsid proteins from Chimpanzee adenovirus. The overhead for decision making remaining low and robust to computing environments. The software is publicly available on GitHub and includes a short user guide to install the R-MDFF on different computing environments, from local Linux based workstations to High Performance Computing (HPC) environments.

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

The eROSITA Final Equatorial-Depth Survey (eFEDS)

In 2019, the eROSITA telescope on board the Russian-German satellite Spectrum-Roentgen-Gamma (SRG) began to perform a deep all-sky X-ray survey with the aim of identifying ~100 000 clusters and groups over the course of four years. As part of its performance verification phase, a ~140 deg 2 survey, called eROSITA Final Equatorial-Depth Survey (eFEDS), was performed. With a depth typical of the all-sky survey after four years, it allows tests of tools and methods as well as improved predictions for the all-sky survey. As part of this effort, a catalog of 542 X-ray selected galaxy group and cluster candidates was compiled. Here, in this paper, we present the optical follow-up, with the aim of providing redshifts and cluster confirmation for the full sample. Furthermore, we aim to provide additional information on the dynamical state, richness, and optical center of the clusters. Finally, we aim to evaluate the impact of optical cluster confirmation on the purity and completeness of the X-ray selected sample. We used optical imaging data from the Hyper Suprime-Cam Subaru Strategic Program and from the Legacy Survey to identify optical counterparts to the X-ray detected cluster candidates. We make use of the multi-component matched filter cluster confirmation tool (MCMF), as well as of the optical cluster finder CAMIRA to derive cluster redshifts and richnesses. MCMF provided the probabilities with which an optical structure would be a chance superposition with the X-ray candidate. These probabilities were used to identify the best optical counterpart as well as to confirm an X-ray candidate as a cluster. The impact of this confirmation process on catalog purity and completeness was estimated using optical to X-ray scaling relations as well as simulations. The resulting catalog was furthermore matched with public group and cluster catalogs. Optical estimators of the cluster dynamical state were constructed based on density maps of the red-sequence galaxies at the cluster redshift. By providing redshift estimates for all 542 candidates, we construct an optically confirmed sample of 477 clusters and groups with a residual contamination of 6%. Of these, 470 (98.5%) are confirmed using MCMF, and 7 systems are added through cross-matching with spectroscopic group catalogs. Using observable-to-observable scaling and the applied confirmation threshold, we predict that 8 ± 2 real systems have been excluded with the MCMF cut required to build this low-contamination sample. This number agrees well with the 7 systems found through cross-matching that were not confirmed with MCMF. The predicted redshift and mass distribution of this catalog agree well with simulations. Thus, we expect that these 477 systems include >99% of all true clusters in the candidate list. Using an MCMF-independent method, we confirm that the catalog contamination of the confirmed subsample is 6 ± 3%. Application of the same method to the full candidate list yields 17 ± 3%, consistent with estimates coming from the fraction of confirmed systems of ~17% and with expectations from simulations of ~20%. We also present a sample of merging cluster candidates based on the derived estimators of the cluster dynamical state.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Performance study of charcoal-based radon reduction systems for ultraclean rare event detectors

The continuous emanation of radon due to trace amounts of uranium and thorium in detector materials introduces radon to the active detection volume of low-background rare event search detectors. 222 Rn produces a particularly problematic background in the physics region of interest by the "naked" beta decay of its 214 Pb daughter nucleus. While charcoal-based adsorption traps are expected to be effective for radon reduction in auxiliary circulation loops that service the warm components of current ton-scale detectors at slow flow rates (0.5-2 SLPM), radon reduction in the entire circulation loop at high flow rates $\mathcal{O}$(100s SLPM) is necessary to reach high sensitivity in future generation experiments. In this article we explore radon dynamics with a charcoal-based radon reduction system in the main circulation loop of time projection chamber detectors. Here, we find that even for perfect radon traps, circulation speeds of 2,000 SLPM are needed to reduce radon concentration in a 10 ton detector by 90%. This is faster by a factor of four than the highest circulation speeds currently achieved in dark matter detectors. We further find that the effectiveness of vacuum swing adsorption systems, which have been employed very successfully at reducing atmospheric radon levels in clean-rooms, is limited by the intrinsic radon activity of the charcoal adsorbent in ultra-low radon environments. Adsorbents with significantly lower intrinsic radon activity than in currently available activated charcoals would be necessary to build effective vacuum swing adsorption systems operated at room temperature for rare event search experiments. If such VSA systems are cooled to about 190 K, this requirement relaxes drastically.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Implementation of a Mesh refinement algorithm into the quasi-static PIC code QuickPIC

Plasma-based acceleration (PBA) has emerged as a promising candidate for the accelerator technology used to build a future linear collider and/or an advanced light source. In PBA, a trailing or witness particle beam is accelerated in the plasma wave wakefield (WF) created by a laser or particle beam driver. The WF is often nonlinear and involves the crossing of plasma particle trajectories in real space and thus particle-in-cell methods are used. The distance over which the drive beam evolves is several orders of magnitude larger than the wake wavelength. This large disparity in length scales is amenable to the quasi-static approach. Three-dimensional (3D), quasi-static (QS), particle-in-cell (PIC) codes, e.g., QuickPIC, have been shown to provide high fidelity simulation capability with 2-4 orders of magnitude speedup over 3D fully explicit PIC codes. In PBA, the witness beam needs to be matched to the focusing forces of the WF to reduce the emittance growth. In some linear collider designs, the matched spot size of the witness beam can be 2 to 3 orders of magnitude smaller than the spot size (and wavelength) of the wakefield. Such an additional disparity in length scales is ideal for mesh refinement where the WF within the witness beam is described on a finer mesh than the rest of the WF. A mesh refinement scheme is described that has been implemented into the 3D QS PIC code, QuickPIC. Very fine (high) resolution is used in a small spatial region that includes the witness beam and progressively coarser resolutions in the rest of the simulation domain. A fast multigrid Poisson solver has been implemented for the field solve on the refined meshes and a Fast Fourier Transform (FFT) based Poisson solver is used for the coarse mesh. The code has been parallelized with both MPI and OpenMP, and the parallel scalability has also been improved by using pipelining. A preliminary adaptive mesh refinement technique is described to optimize the computational time for simulations with an evolving witness beam size. Several test problems are used to verify that the mesh refinement algorithm provides accurate results. Additionally, the results are benchmarked against highly resolved simulations exhibiting near-azimuthal symmetry, performed using QPAD—a novel hybrid QS PIC code that uses a PIC description in the coordinates (r, ct – z) and a gridless description in the azimuthal angle, Φ.

Linear collider↗
DOE OSTI2020

Models and Strategies for Optimal Demand Side Management in the Chemical Industries

Deregulation and the increase of renewable electricity generation from wind and solar photovoltaics have transformed the U.S. electricity market. Economic and environmental benefits notwithstanding, the presence of renewables has increased variability and uncertainty on the supply side of the grid. Managing demand, rather than generation – a strategy referred to as “demand response (DR)” – is an attractive approach for mitigating this imbalance. DR efforts aim to reduce electricity usage during peak demand times, lessening stress on the grid. Industrial users are particularly attractive entities for DR participation since they present large, localized loads that can provide significant relief on grid demand and –unlike other large loads, such as buildings – are minimally dependent on human needs and preferences. In this project, we accomplished three main objectives. (1) We developed data-driven low-order DR scheduling-relevant dynamic models of chemical processes. Concurrently, we studied the formulation and solution of the associated optimal DR production scheduling problems. (a) A prototype air separation unit (ASU) model was used to generate simulated operating data for initial modeling efforts, which enabled the later use of industrial data for data-driven modeling. (b) We utilized Hammerstein-Wiener (HW) and Finite Step Response (FSR) models to represent nonlinear plant dynamics. (c) The HW models were linearized using exact linearization so they could potentially be embedded in power system models, which are formulated as mixed integer linear programs (MILPs). (d) We solved DR optimization problems under uncertainty and found that even naïve predictions of electricity price and product demand led to significant cost savings benefits. (2) Our DR scheduling optimization problem formulations are amenable to real-time solution. (a) We utilized Lagrangian Relaxation (LR) to efficiently solve the optimization problem by decoupling subproblems linked by complicating constraints. (b) We have achieved computation times for the 3-day DR scheduling problem of an ASU as low as 1.88 minutes. (3) Our representations of the DR behavior of chemical process as grid-level batteries were embedded in power system models. (a) For a small-scale grid, we found that incorporating the dynamics of the chemical plant in the optimal power flow calculations resulted in better resource management leading to up to 15% and 46% cost reduction for the grid and chemical plant operations, respectively, during periods of power line congestion. We have published several works dedicated to modeling and solving DR optimization problems from the user side. These were published in top peer-reviewed journals and are summarized in this report. The most recent work (and papers in preparation) considers DR scheduling from the grid side. Future efforts will consider networked plants (e.g., air separation units operating on a common pipeline) for DR participation, which is expected to amplify the capabilities of industrial DR participants to perform load-shifting. Our consideration of uncertainty in DR has inspired future directions in this area as well: we plan to develop multistage methods to fully account for the effects of uncertainty in DR scheduling.

24 POWER TRANSMISSION AND DISTRIBUTION↗