Search NASA⌕ Search

SEARCH · Search NASA

Results for “LOADS, STEADY - WINGS”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Active Site Dynamics in Molybdenum-Based Silica-Supported Olefin Metathesis Catalysts: Site Renewal and Decay Beyond the Chauvin Cycle

Heterogeneous olefin metathesis catalysts exhibit low active site densities and unpredictable kinetics due to dynamic active site formation and decay processes. Here, in this study, we establish a quantitative framework that captures active site generation, renewal, and decay in olefin metathesis over silica-supported molybdenum oxide catalysts, enabling a mechanistic explanation of catalytic behavior and strategies to achieve high, stable activity. Steady-state active site titrations reveal that 2,3-dimethyl-butene isomers (4MEs) cofeeding increases active site density by up to 4.3-fold, directly correlating with enhanced metathesis rates. Spectroscopic studies demonstrate that 4MEs facilitate Mo(VI) reduction to Mo(IV) and interact strongly with surface Si–OH groups, generating labile protons that promote active site formation via a 1,2-proton shift mechanism. Kinetic modeling indicates that ethylene acts as a decay promoter, shifting kinetic control away from the Chauvin cycle and suppressing metathesis activity. Comparative studies on catalysts with varying Mo loading reveal that promotion is most effective for dispersed molybdate species, with a decline at higher Mo loadings. These findings provide a unified mechanistic framework for heterogeneous olefin metathesis, offering new strategies to enhance active site accessibility, mitigate deactivation, and optimize catalyst design.

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

Analyzing the Effects of Atmospheric Turbulent Fluctuations on the Wake Structure of Wind Turbines and Their Blade Vibrational Dynamics

In recent trends, a rising demand for renewable energy has driven wind turbines to larger proportions, where lighter blade designs are often adopted to reduce the costs associated with logistics and production. This causes modern utility-scale wind turbine blades to be inherently more flexible, and their amplified aeroelastic sensitivity results in complex multi-physics reactions to variant atmospheric conditions, including dynamic patterns of aerodynamic loading at the rotor and vortex structure evolutions within the wake. In this paper, we analyze the influence of inflow variance for wind turbines with large, flexible rotors through simulations of the National Rotor Testbed (NRT) turbine, located at Sandia National Labs’ Scaled Wind Farm Technology (SWiFT) facility in Lubbock, Texas. The Common Ordinary Differential Equation Framework (CODEF) modeling suite is used to simulate wind turbine aeroelastic oscillatory behavior and wind farm vortex wake interactions for a range of flexible NRT blade variations, operating in differing conditions of variant atmospheric flow. CODEF solutions of turbine operation in Steady-In-The-Average (SITA) wind conditions are compared to SITA wind conditions featuring a controlled gust-like pulse overimposed, to isolate the effects of typical wind fluctuations. Finally, simulations of realistic time-varying wind conditions from SWiFT meteorological tower measurements are compared to the solutions of SITA wind conditions. These increasingly complex atmospheric inflow variations are tested to show the differing effects evoked by various patterns of spatiotemporal atmospheric flow fluctuations. An analysis is presented for solutions of wind turbine aeroelastic response and vortex wake evolution, to elucidate the consequences of variant inflow, which pertain to wind turbine dynamics at an individual and farm-collective scale. The comparisons of simulated farm flow for SITA and measured fluctuating wind conditions show that certain regions of the wake contain up to a 12% difference in normalized axial velocity, due to the introduction of wind fluctuations. The findings of this study prove valuable for practical applications in wind farm control and optimization strategies, with particular significance for modern utility-scale wind power plants operating in variant atmospheric conditions.

Farrell, Alayna (ORCID:000000023555720X)↗
DOE OSTI2020

Hydrodynamics of countercurrent flows in a structured packed column: effects of initial wetting and dynamic contact angle

The structured packed column has gained considerable attention as a prospective technology to mitigate carbon dioxide emissions from thermal power plants. The computational flow investigation of such systems is a multiscale problem, and local hydrodynamics plays a key role in overall column efficiency. Accordingly, multiphase flow studies in the representative elementary unit (REU) of a packed column can offer insights into the local flow behavior such as flow pattern, interfacial area, etc. The interfacial area plays a key role in the overall interphase mass transfer and thereby in separation efficiency. The predicted effective area matches well with corresponding ones in experiments for Mellapak 250.Y packings. The contact angle (i.e., solid substrate behavior) is one of the critical factors that dictates wettings, and thereby the interfacial area. The dynamic contact angle (DCA) was considered in order to explore the effects of contact angle hysteresis on the interfacial area. DCA has more pronounced impact on the interfacial area for a solvent possessing a higher value of surface tension and equilibrium contact angle. The interfacial area shows temporal undulation and does not achieve a pseudo steady state due to the capillary force. In contrary, the interfacial area gets a net value for a solvent having low surface tension value. Further, specification of the static contact angle at the solid substrate shows the least value of the interfacial area. We also extensively studied the effect of the initial sheet condition (dry vs wet) on the interfacial area. The initially wetted sheets show slightly higher value of the interfacial area as compared to the initially dry sheets at a fixed flow rate for a solvent having high surface tension value. Wetting hysteresis decreases with increasing liquid loads. On the other hand, wetting hysteresis does not occur for a solvent having value of low surface tension and equilibrium contact angle.

structured packings, dynamic contact angle, interf↗
DOE OSTI2024

THERMAL STRUCTURAL AND SHOCK EVENT EVALUATIONS OF THE FUELING PELLET INJECTION SYSTEM FOR ITER

The US is among seven partner nations in a collaborative effort to design, build, and demonstrate fusion’s ability to be a large-scale carbon free energy source. Each country has its own Domestic Agencies (DA) that contribute directly to the ITER project. US ITER, which is a DOE Office of Science project managed by Oak Ridge National Laboratory, is developing world class engineering solutions to the design, construction, and assembly of the burning plasma experiment that can demonstrate the scientific and technological feasibility of fusion. US ITER’s scope includes completing the preliminary and final design, qualifying materials and processes for manufacture and testing, executing manufacturing, and delivering the Fueling Pellet Injection System (FPIS) to the ITER site for assembly. The US-ITER FPIS system is designed to inject cryogenically frozen pellets of deuterium-tritium (D-T) into the plasma. The FPIS has two main functions: 1. Provide a steady supply of deuterium and tritium fuel, 2. Mitigating the impact of edge localized modes on the plasma facing components. The FPIS will have the capability to reside in three port cells within a tritium second barrier containment cask. Each pellet cask contains three flight tubes linking with penetrations on the torus cryopump housing and vacuum vessel (VV): two for the magnetic high field side (HFS) pellet injection and one for the magnetic low field side (LFS) pellet injection. This paper presents results of the thermal-structural analyses of the FPIS flight tube structural components when subjected to various loads such as electromagnetic (EM) loads, nuclear heating, seismic, and operational, and dynamic shock events. The resulting temperatures and stresses under combined conditions have been found to satisfy the design criteria to ensure safe and reliable operation of the FPIS flight tubes within the vacuum vessel.

Martinez, Oscar↗
DOE OSTI2023

$\mathrm{RADICAL}$-Pilot and $\mathrm{PMIx}$/$\mathrm{PRRTE}$: Executing Heterogeneous Workloads at Large Scale on Partitioned $\mathrm{HPC}$ Resources

Execution of heterogeneous workflows on high-performance computing (HPC) platforms present unprecedented resource management and execution coordination challenges for runtime systems. Task heterogeneity increases the complexity of resource and execution management, limiting the scalability and efficiency of workflow execution. Re-source partitioning and distribution of tasks execution over portioned re-sources promises to address those problems but we lack an experimental evaluation of its performance at scale. Here this paper provides a performance evaluation of the Process Management Interface for Exascale (PMIx) and its reference implementation PRRTE on the leadership-class HPC plat-form Summit, when integrated into a pilot-based runtime system called RADICAL-Pilot. We partition resources across multiple PRRTE Distributed Virtual Machine (DVM) environments, responsible for launching tasks via the PMIx interface. We experimentally measure the work-load execution performance in terms of task scheduling/launching rate and distribution of DVM task placement times, DVM startup and termination overheads on the Summit leadership-class HPC platform. Integrated solution with PMIx/PRRTE enables using an abstracted, standardized set of interfaces for orchestrating the launch process, dynamic process management and monitoring capabilities. It extends scaling capabilities allowing to overcome a limitation of other launching mechanisms (e.g., JSM/LSF). Explored different DVM setup configurations provide insights on DVM performance and a layout to leverage it. Our experimental results show that heterogeneous workload of 65,500 tasks on 2048 nodes, and partitioned across 32 DVMs, runs steady with resource utilization not lower than 52%. While having less concurrently executed tasks resource utilization is able to reach up to 85%, based on results of heterogeneous workload of 8200 tasks on 256 nodes and 2 DVMs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Isochronous Architecture-Based Voltage-Active Power Droop for Multi-Inverter Systems

This article proposes an isochronous architecture for parallel inverters with only voltage-active power droop (VP-D) control for improving active power sharing as well as plug-and-play of multi-inverter-based distributed energy resources (DERs). The isochronous framework obviates the need for an explicit regulation of the frequency while allowing for sharing of reactive power. The article shows that the detrimental effects of circulating currents between inverters can be addressed in the framework developed. The isochronous architecture is implemented by employing a GPS to disseminate clock timing signals that enable the microgrid to maintain the nominal system frequency. Small signal eigenvalue analysis of a Multi-inverter Microgrid system near the steady-state operating point is presented to evaluate the system stability. Moreover, unlike traditional droop-based methods, the isochronous strategy lends itself to analytical guarantees which are developed in the article. Furthermore, the effect of delays in receiving GPS clock signals by an inverter is experimentally characterized and an inverter isolation mechanism is proposed that disconnects inverters facing large GPS communication delays. Hardware experiments on an 1.2 kVA-prototype and controller-hardware-in-the-loop validation on the CIGRE distribution network are conducted to demonstrate the effectiveness of the proposed architecture towards active and reactive power sharing between inverters with load scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Neutron irradiation & thermomechanical experiment (NITE) - design

For the reliable long-term operation of fusion power plants, it is crucial to understand and predict the lifetime of materials in use. These materials include all structural and functional materials utilized at the first wall, blanket, magnets, and shielding. The key challenge is, that the harsh environment including high heat fluxes, high thermal stress and stress cycling, neutron irradiation, and sputtering on such materials should not be viewed separately. Currently, the synergistic loads cannot be evaluated experimentally because of the lack of adequate facilities. The purpose of that work is to design a synergetic Neutron Irradiation and Thermomechanical Experiment (NITE) for fusion materials. This design will leverage the existing Advanced-Test-Reactor (ATR), a fission reactor at the Idaho National Laboratory. We also acknowledge that with existing fission reactors the exact fusion condition cannot be created, and the limitations are critically discussed. The combination of neutron irradiation with a high heat flux is the focus. This is realized with an irradiation capsule design that includes a TRISO fueled region inside the capsule to enable a steady-state heat flux on one side of the specimen. In conclusion, the experimental design modeling showed that steady-state heat fluxes of 2.4 MW/m 2 with a thermal gradient of above 250°C can be achieved in a 5 mm thick specimen.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Report on Year-4 of Water NSTF Matrix Testing: Facility Maintenance and Accident Testing

Under support from the Department of Energy (DOE) and the Office of Advanced Reactor Technologies (ART), a large-scale test facility has been constructed at Argonne National Laboratory to generate NQA-1 qualified validation data for passive decay heat removal systems in advanced reactors. The Natural convection Shutdown heat removal Test Facility (NSTF) reflects key features of a ½ scale, water-based, Reactor Cavity Cooling System (RCCS) and is intended to study the behavior, bound performance, and ultimately guide design decisions for passive decay heat removal systems for advanced reactors. In addition to the experimental activities detailed in this report, a supportive computational modeling effort is on-going which has been demonstrated to significantly strengthen the experimental program while also improving accuracy of the computer models. Together these create a mutually beneficial relationship integral to meeting the overall program objective of examining the heat removal performance of the RCCS concept. This report serves as a summary of maintenance and experimental activities during the program’s fourth year of water-based operation. A planned six-month maintenance period began in August 2021, during which major inspections, repairs, cleaning, and installation of new instrumentation and data acquisition hardware were conducted. Most significantly, two heaters that faulted during Year-3 were repaired, allowing the facility to resume use of the full heated section area and full range of available electric power. The remainder of the year consisted of eight months of test operations, during which the facility logged 211 hours of active heating across one bake-out (following the maintenance period) and seven matrix tests; five classifieds as Accepted per NQA-1, one as Trending, and one as Failed. Testing began by performing two repeat cases to confirm expected facility response and behavior during both single- and two-phase flow conditions, ensuring no changes were introduced during the maintenance period that might have altered the thermal-hydraulic characteristics of the facility. In continuation of the power parametric series initiated in previous years, a high-power test case was then performed examining heat removal performance at a decay heat load equivalent to 2.4 MWt, full-scale, a level exceeding maximum design targets. Additional testing then introduced various blockage scenarios along the network piping, examining the effects of partial and complete blockages of the flow paths on the system behavior and heat removal performance. A study of static boiling tests directed at understanding the geysering two-phase instability was also conducted. The loop was filled only to the bottom of the tank outlet, creating an open loop configuration that prevents any natural circulation flow from occurring, and the heaters were powered on until the facility reached saturation conditions. Following, a series of quasi-steady-state conditions were introduced by adjusting the inventory level in the adiabatic chimney piping at decreasingly lower elevations above the heated region. A strong correlation of geysering characteristics and loop level was observed, with flow and temperature excursions decreasing in intensity, but increasing in frequency, as the fill were reduced to lower elevations along the chimney piping. Once the level fell very low in the chimney, at points near the top of the heated section, the system reached a stable state of continuous boiling without any occurrence of geysering eruptions. A final significant testing accomplishment this year was successful completion of an “accident scenario” test, whose operating conditions were based on a prototypic decay heat curve provided by Framatome and scaled for the NSTF. This test began by establishing steady-state, single-phase “normal operation” conditions, before simulating an accident trip where the availability of active cooling systems was lost. Loop temperatures gradually increased until reaching saturation and subsequent two-phase boiling flow. Over the course of an extended operational period along the defined decay heat curve, steam boil-off caused gradual but continued depletion of liquid inventory until reaching a critically low level causing flow stagnation and cessation of natural circulation heat removal. At this point, after nearly 72 hours of continuous operation, a cold refill was performed to replenish the system inventory and allow the facility to re-establish closed loop natural circulation flow and return to a safe operational state.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

The Effect of Spark-Plug Heat Dispersal Range and Exhaust Valve Opening Timing on Cold-Start Emissions and Cycle-to-Cycle Variability

The partnership for advancing combustion engines (PACE) is a US Department of Energy consortium involving multiple national laboratories and includes a goal of addressing key efficiency and emission barriers in light-duty engines fueled with a market-representative E10 gasoline. A major pillar of the initiative is the generation of detailed experimental data and modeling capabilities to understand and predict cold-start behavior. Cold-start, as defined by the time between first engine crank and three-way catalyst light-off, is responsible for a large percentage of NOx, unburned hydrocarbon and particulate matter emissions in light-duty engines. Minimizing emissions during cold-start is a trade-off between achieving faster light-off of the three-way catalyst and engine out emissions during that period. In this study, gaseous and soot emissions were measured at a distance representative of the three-way catalyst position downstream of the engine at a 2 bar net indicated mean effective pressure (NIMEP) steady-state operating condition representative of cold-start. The test matrix included sweeps of ignition timing 15 degrees-before to 10 degrees-after top dead center firing (TDCf) across three different spark-plug heat dispersal ranges (HR). Additionally, the effect of varying exhaust valve opening (EVO) timing on combustion stability and emissions was also studied. Results show that the spark plug HR affects the coefficient of variation (COV) of NIMEP under all cold-start conditions, while the impact on emissions was found to be minimal. At very retarded spark timings, colder spark plugs required higher air and fuel flow to maintain the desired 2bar NIMEP load, but the fraction of fuel energy going into the exhaust was found to be similar for all spark plugs. Finally, retarding exhaust valve timings showed a simultaneous reduction in emissions while increasing the fraction of fuel energy being fed into the exhaust. However, engine COV was also observed to increase with retarded exhaust timings.

42 ENGINEERING↗
DOE OSTI2020

Understanding Microstructural Effects on Dynamic Performance Towards the Development of Shock Metamaterials

With the recent advances in additive manufacturing (AM), long-range periodic lattice assemblies are being developed for vibration and shock mitigation components in aerospace and military applications with unique geometric and topological structures. There has been extensive work in understanding the static properties associated with varying topology of these lattice architectures, but there is almost no understanding of microstructural affects in such structures under high-strain rate dynamic loading conditions. Here we report the shock behavior of lattices with varying intrinsic grain structures achieved by post process annealing. High resolution 316L stainless steel lattices were 3D printed by a laser-powder bed fusion machine and characterized by computed tomography. Subsequent annealing resulted in stress-relieved and recrystallized lattices. Overall the lattices had strong cubic texture aligning with the x-, y- and z-directions of the build with a preference outside the build direction (z). The recrystallized sample had more equiaxed polygonal grains and a layer of BCC ferrite at the surface of the structure approximately 1 grain thick. Upon dynamic compression the as-deposited lattice showed steady compaction behavior while the heat-treated lattices exhibit negative velocity behavior indicative of failure. We attribute this to the stiffer BCC ferrite in the annealed lattices becoming damaged and fragmenting during compression.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Effect of hydrogen isotope on plasma impedance and thermal pinching induced by single aerosol droplets injected into an inductively coupled plasma

The mass of the hydrogen isotope in micrometer-sized aqueous droplets introduced into an ICP has been found to affect the degree of resulting plasma perturbation. Studied plasma perturbations include local plasma cooling, plasma shrinkage due to thermal pinching, and changes in plasma impedance. The local cooling effect caused by a vaporizing aerosol droplet was found to be similar between H 2 O and D 2 O. This finding is likely a result of the similar heat capacities and vaporization rates of H 2 O and D 2 O droplets and of the comparable thermal conductivities and dissociation kinetics of H 2 O and D 2 O vapor. In contrast, a clear isotope effect was observed for plasma shrinkage and impedance change, with H 2 O causing a considerably stronger effect. Plasma shrinkage, gauged by the reduction in Ar emission at the edge of the plasma, was detected at all observation heights within the load-coil region, with D 2 O droplets causing as much as 30% less of an effect than H 2 O. Likewise, D 2 O droplets induce significantly less change in plasma impedance than H 2 O droplets. Depending on experimental conditions and the probing method, changes in plasma impedance caused by monodisperse D 2 O droplets were only 52%–80% that of H 2 O. The lower influence on plasma shrinkage and impedance change for D 2 O droplets is believed to be related to the much lower thermal conductivity of atomic deuterium, around 70% that of atomic hydrogen (protium). The time required for the plasma impedance to return to its original steady-state value was comparatively long (~10 ms) and similar for H 2 O and D 2 O droplets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Competition between roughness and strength for scale-dependent surfaces

Rocks famously have scale-dependent strength, yet the actual dependence is notoriously hard to measure or incorporate into any theoretical framework. Natural rough surfaces present an opportunity to solve the problem. Surfaces sliding in shear evolve as protrusions collide. These asperities can deform or break, thus creating a new surface shape. In particular, natural surfaces have roughness at all scales as well as scale-dependent strength. Based on a scaling analysis, we have previously suggested that the scale-dependent aspect ratio of steady-state surfaces should be proportional to the scale-dependent shear strain at yield. If true, scale-dependent strength could easily be inferred from natural surfaces. Thus, moving beyond the scaling argument to a rigorous treatment of scale-dependent strength for multiscale rough surfaces in shear is important. However, analytic frameworks for analyzing multiscale problems are challenging, as conventional continuum mechanics typically involves a single value for a material property across scales. Here, in this work, we build on the formalism of Persson (2001) that presents a method to compute contact area for rough surfaces with a prescribed topographic spectrum using a stochastic differential equation. The Persson formalism allows for plastic yield under normal loading of otherwise elastic materials and leaves open the possibility of scale-dependent yield stress. In this study, we pursue this route to develop a theory and numerical results for the yielding of a rough, elastoplastic surface with scale-dependent yield stress. Here, we examine surfaces for which the power spectrum of the topography 𝐶 and yield stress 𝑌 follow power laws as a function of scale 𝜆, such that 𝐶∼𝜆 −𝑚 and 𝑌∼𝜆 −𝑛 , respectively. In this formal treatment of the problem, we focus on surfaces in contact and the resulting yield and do not impose shear. Numerical solutions show that the deviation from the elastic scaling solution is bounded as expected by the prior 1D heuristic scaling argument that anticipates the Hurst exponent as 1−𝑛. We also show that the plasticity is expected to erode the contacts if 𝑚 is lower than 𝑛−3, which corresponds to a Hurst exponent lower than 1−𝑛/2. This result is rigorously sound for 2D, i.e., realistic surfaces, and quantitatively different than the prior scaling argument. The theory now permits a correspondingly quantitative approach to interpreting natural surfaces.

elasticity↗
DOE OSTI2025

Lower Snake River Dams Contribution to Grid Services

Hydroelectric generation and water storage have long been components of the clean energy mix, providing both reliable steady output and operational flexibility. As variable renewable energy sources such as wind and solar increasingly replace traditional generation, the role of all flexible resources—including hydropower—in balancing supply and demand continues to evolve. This study examined the contribution of the Lower Snake River (LSR) Dam plants to Bonneville Power Administration grid services in maintaining power system reliability within the Western Interconnection. By analyzing publicly available data, the study evaluated various reliability services through performance metrics including energy capacity, balancing and ramping, voltage and reactive power support, frequency response, and transmission impact. Results indicated that the LSR plants deliver services as expected based on their size, contributing to the balancing process, ramping capabilities, and operational reserves, particularly during peak load conditions and weather events, while also providing measurable frequency and voltage support to the grid.

13 HYDRO ENERGY↗
DOE OSTI2026

EXERGETIC: De-Risking Next-Generation Resilient Geothermal Hybrids via At-Scale Evaluation Using Virtual Emulation Digital Twin Environment for Efficient Operation

The DOE-GTO-funded project, award number 5.1.2.12, entitled "EXERGETIC - De-risking Next Generation Resilient Geothermal Hybrids via at-Scale Evaluation Using a Virtual Emulation Digital Twin Environment for Efficient Operation," advances the solution to these challenges by developing and validating a geothermal co-emulation environment implemented at the National Laboratory of the Rockies (NLR)'s Advanced Research on Integrated Energy Systems (ARIES) platform. This framework enables the de-risking of next-generation geothermal and geothermal hybrid systems through high-fidelity modeling, real-time digital emulation, advanced control strategies, and techno-economic assessment. The project focused on geothermal hybrid configurations that integrate geothermal power plants with concentrated solar power and underground thermal energy storage, enabling enhanced efficiency, flexibility, and grid support capabilities. The main goal of this project was the development of a geothermal digital co-emulation environment to demonstrate the technical and economic value of geothermal hybrid systems and their contribution to grid stability and flexibility. The EXERGETIC framework combined physics-based models, controls, and real assets at ARIES, including digital real-time simulators (DRTS), a 20-MW-scale controllable grid interface (CGI), and a 2-MW conventional generator. Detailed transient models were developed for the key subsystems of a hybrid geothermal plant, including parabolic trough solar collectors, reservoir thermal energy storage (RTES), and a binary Organic Rankine Cycle (ORC) power plant. The ORC model explicitly captured thermal inertia and off-design operation and integrated control strategies to dynamically respond to electric load profiles. The models were validated against published experimental and numerical studies, demonstrating strong agreement and confirming the accuracy and robustness of the modeling approach. The resulting digital twin represents geothermal-solar-storage systems at multiple scales (1 MW to 100 MW) and enables realistic emulation of grid-connected operation. The control architecture allows the geothermal resource to provide stable baseload generation, while solar and stored thermal energy supply flexible, dispatchable support during periods of high demand or variable grid conditions. A key contribution of the EXERGETIC project is the demonstration that geothermal hybrid systems can be designed to be active grid assets rather than passive baseload generators. Using the ARIES platform, the digital twin was evaluated under multiple grid scenarios, including load following, voltage support at the distribution level, and frequency response at the transmission level. Results show that hybrid geothermal systems can respond effectively to dynamic grid conditions, providing inertia-like behavior, primary frequency support, and voltage regulation through coordinated control. In addition to the performance and grid services capability analysis of geothermal and hybrid geothermal systems, the EXERGETIC project also focused on scalability and techno-economic analysis of geothermal hybrid plants. In particular, for the scalability analysis, machine-learning (ML)-based surrogate models were trained using data generated from the geothermal digital twin under different grid-connected scenarios and plant capacities. These ML models demonstrated strong interpolation and extrapolation capabilities across plant sizes, accurately reproducing both steady-state and transient responses with very low errors. Regarding the techno-economic analysis, plant performance results were integrated with cost models for hybrid geothermal systems, and the levelized cost of electricity (LCOE) was used as the main economic metric to evaluate system performance across a range of system capacities, solar shares, solar multiples, and storage durations. Results indicate that economies of scale significantly reduce geothermal LCOE as plant capacity increases, with large-scale systems (25-100 MW) achieving substantially lower costs than small plants. Hybridization with solar thermal energy and storage further improves economic performance by increasing capacity utilization and enabling flexible dispatch. In addition, thermal storage plays a critical role in reducing LCOE by maximizing geothermal, solar, and stored energy resources. In summary, the results from this project demonstrate that geothermal hybrid systems represent a promising alternative for increasing the energy conversion efficiency of geothermal technologies, contributing to the preservation of geothermal resources, and supporting the transition of geothermal plants from traditional baseload resources into flexible, resilient, and cost-competitive energy conversion technologies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

A Preliminary Investigation on the Performance of a Dual-Source Heat Pump using both the Air and the Ground

Air source heat pumps (ASHPs) and ground source heat pumps (GSHPs) are the two most common types of electric-driven heat pumps in the marketplace to replace fossil fuel-based heating systems. However, the performance and efficiency of ASHPs depend on the ambient air conditions. Therefore, ASHPs usually are equipped with electric resistance heaters to provide supplemental heating when the ambient temperature is low, and the heating demand is high. The electric resistance heaters could result in high power draws when they are turned on. On the other hand, due to the relatively steady temperature of the ground, GSHPs are more energy efficient than ASHPs when providing space heating and cooling to the buildings. However, the adoption of GSHP is hindered by its high initial cost, mostly due to the cost of drilling boreholes for installing ground heat exchangers (GHE). To solve the above issues, our study investigates the performance of dual-source heat pumps (DSHPs) with respect to that of ASHPs and GSHPs. The DSHP will use ambient air when its temperature is favorable for the efficient operation of the heat pump. When the ambient temperature is too hot or too cold, the ground source will be utilized to retain the high-efficiency operation of the heat pump. Because the cumulative thermal load of the GHE is shared with the ambient air, the size of the GHE could be smaller than those of the conventional GSHPs. The study will model the DSHP and simulate its performance in providing heating and cooling for a typical single-family house in regions with hot, mild, or cold climates. In addition, the required GHE size of the DSHP system will be determined through annual simulations and compared with that of conventional GSHPs.

Anees, Fady↗
DOE OSTI2021

Thermal loading analysis of the ring injection dump for the Spallation Neutron Source facility

The ring injection dump (RID) is the largest beam dump in the Spallation Neutron Source facility at Oak Ridge National Laboratory and accepts a fraction of the beam from the ring that is not captured in the ring during injection. Thermo-fluid modeling of the RID components was performed during the initial design by using the software ANSYS CFX to characterize the power rating for the RID. From the calculations, the power rating was lowered from 200 kW to 150 kW, at 1.3 GeV beam energy due to concerns regarding the heating of the concrete structure. The Proton Power Upgrade project scope included reevaluating the power limit to the RID in anticipation that increasing it would provide more operational flexibility. The focus was on the shielding because it was the limiting factor, and data collected from thermocouples installed on the structure were used to benchmark the analysis. A three-stage process was adopted in the latest study. First, the already existing steady-state calculations were validated and extended with the latest Monte Carlo N-Particle source term calculations. Second, transient calculations were conducted to capture the dynamic state of the system in response to the energy dumped to the RID during the following operational period of the beam from 2005–2019 (14 years), and sensitivity analyses of crucial parameters were performed to benchmark the model with reasonable accuracy. Finally, an idealized transient cycle analysis was conducted with realistic duty factors to predict the temperature distribution at higher beam powers. Overall, the results confirmed the accuracy of the original steady state model but allowed the development of a new validated transient model with higher accuracy for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

The physics basis to integrate an MHD stable, high-power hybrid scenario to a cool divertor for steady-state reactor operation

Abstract Coupling a high-performance core to a low heat flux divertor is a crucial step for ITER and a Fusion Pilot Plant or DEMO. Experiments in DIII-D recently expanded the steady-state hybrid scenario to high density and divertor impurity injection to study the feasibility of a radiating mantle solution. This work presents the physics basis for trade-offs between density, current drive and stability to tearing modes (TMs) at high β. EC power is crucial to tailor the plasma profiles into a passively stable state, and to eject impurities from the core. Off-axis EC depositions decrease the heating efficiency, but calculated electron heat transport coefficients show that this effect is partially mitigated by improved confinement inside the EC deposition. Additionally, the reduction in pressure is recovered by increasing the density. This favourable scaling of confinement with density was observed in high power plasmas for years, and this work provides a comprehensive explanation. ELITE predictions indicate that a path in peeling-ballooning stability opens up for certain conditions of density, power, q 95 and shaping, allowing the edge pressure to continue increasing without encountering a limit. In the core, calculated anomalous fast-ion diffusion coefficients are consistent with density fluctuation measurements in the toroidicity-induced Alfvén eigenmode range, showing that smaller fast-ion losses contribute to the enhanced confinement at high density. The edge integration study shows that divertor heat loads can be reduced with Ne and Ar injection, but this eventually triggers a cascade of n = 1, 2, 3 core TMs. We can now show that impurity radiation in the core is small and it is not the cause for the drop in confinement at high Ar and Ne injection rates. The overlap between the core TMs is consistent with the loss of pressure as estimated by the Belt model for the coupled rational surfaces. Optimization of these trade-offs has achieved plasmas with sustained H 98y2 = 1.7, f GW = 0.7 and ∼85% mantle radiation. The scenario and its variations at higher density and on- vs off-axis EC heating has been studied as a candidate for an integrated solution for several reactor designs, such as ITER, ARC, and the ARIES-ACT1 case, showing promising results in terms of fusion power and gain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗