Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation”

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 253 records · Page 14

Integrity Monitoring, Prediction and Assessment, Corrosion Control, and Repair of the Hanford Storage Tanks – 25213

The Hanford Nuclear Reservation site contains approximately 211 million liters of radioactive and chemically hazardous waste arising from nuclear weapons production, beginning with World War II, and continuing through the Cold War [1]. The waste is stored in 177 carbon-steel underground storage tanks, of which 149 are single-shell tanks (SSTs) and the remaining are double-shell tanks (DSTs). The mission of an ongoing River Protection Project is to retrieve the waste from the underground storage tanks and then treat and immobilize (i.e., vitrify) it for disposal. Waste from the older SSTs is being progressively retrieved into the newer DSTs for storage pending treatment, immobilization, and disposal. Figure 1 depicts a typical DST design [2]. The tank is approximately 23 m in diameter and 9 m high and has a domed structure and has a capacity of 4000 m3. The tank wall and floor vary in thickness between approximately 10 mm and 25 mm depending on location. The thicker wall sections are near the curved transition between the tank wall and the floor, while the thinner sections are located near the top of the tank wall. The tank floor thickness varies from 25 mm at the tank center to 10 mm near the tank wall. The tanks were constructed of either ASTM A516 Grade 65 or ASTM A537 Class 1 carbon steel and were post-weld heat treated to reduce the risk of SCC.

Shukla, Pavan K. [Savannah River National Laborato↗

Window Cooling Studies and Disk Vibration Testing on a Subset of Mo-100 Disks

Production of metastable Technetium-99 (Tc-99m), a radioactive tracer that emits gamma rays, is vital to the medical imaging community. Tc-99m is extracted from the decay of Molybdenum-99 (Mo-99) which has a half-life of about 2-3 days. The work presented in this report is part of the NNSA’s mission to produce Mo-99 commercially, within the US, without the use of highly enriched uranium (HEU) in support of nonproliferation and global security. Los Alamos National Laboratory (LANL) is working with NorthStar Medical Radioisotopes (NMR) on their efforts to produce Mo-99 through the irradiation of Mo-100 targets using an electron beam. The NMR target comprises a stack of approximately 60-70 Mo-100 disks with diameter 24 mm and thickness 0.74 mm held in stainless steel laminations, each separated using 0.25 mm thick stainless-steel spacers. The symmetric target stack is housed in an Inconel vessel with two Inconel windows on either side. Two electron accelerators are used to produce 40 MeV, 3.16 µA electron beams each that penetrate the Inconel windows and irradiate the Mo-100 disks. Approximately 90% of the total 250 kW beam power is deposited in the NMR target during the irradiation process, with a smaller percentage adding up to 2.2 kW of heat deposited on the Inconel window. During irradiation, pressurized helium gas flows through thin gaps between the disks cooling the beam window, target disks, disk laminations and spacers. Both NMR and LANL have found during cold testing of the target system (no heat deposition) that the Mo 100 disks undergo significant mass loss and disk breakage due to vibrations induced by the flowing helium gas. The mass loss is not only undesirable due to monetary loss from reduced final quantities of Mo-99, but also due to the hazards associated with radioactive material trapped in the cooling lines and particle filters. The effect of flow rate and target geometry on the flow induced vibrations need to be quantified, and recommendations provided to minimize this mass loss. LANL has previously also tested NMR’s Inconel beam window by heating the window, while flowing pressurized helium, using the average heat deposited on the window. However, the NMR beam is pulsed with a duty cycle of 12.5%, which introduces oscillation in temperature around the nominal 600 °C steady state value with each pulse. Available fatigue curves for Inconel are few, established for room temperature, and they are based on mechanical strain cycles not thermally induced strain as in the NMR target. The effect of pulsed beam heating on the Inconel window therefore needs to be quantified. This report details the experiments conducted to assess the factors that lead to mass loss in the NMR target disks as well as to understand the effect of a pulsed beam on NMR’s Inconel window. This work describes LANL’s experimental characterization of the flow induced vibrations and disk mass loss in a reduced scale set-up containing 5 to 10 Mo-100 disks. We use high speed imaging, displacement measurements and microphone measurements combined with signal processing to estimate the vibration frequency of each disk. The effect of disk thickness, target fit and duration of testing on the mass loss is described. We find that in the current configuration of NMR targets, the vibrations and mass loss on the first disk are minimized, while those in the adjacent disks are highest. The microphone and high-speed image data show that increased flow rates and increased duration of testing increases vibration frequency and mass loss. The mass loss is due to both disk rotation and back and forth motion. There are visible wear marks on the disks with the highest mass loss. We also note that the current NMR window gap reduces flow induced vibrations compared to the previous smaller gaps. Improved target holders significantly reduce disk mass loss to almost negligible quantities. This work finds that the larger window to first disk gap and improved target holder geometry should allow NMR to successfully conduct irradiations with minimal mass loss. The window tests were conducted to understand the effect of a pulsed beam on both the window longevity and to estimate the window temperature and displacement during pulsing. The experiments presented here were performed at significantly low power, due to the limitations of the induction heating system. The window temperature rose to approximately 73 °C with a significantly reduced power of 45 W without beam pulsing. With a 5 Hz pulse rate, 12.5% duty cycle, the window temperature remained constant at 26 °C. These experiments will be repeated with improved coil geometry and reported in upcoming journal papers.

42 ENGINEERING↗

4.3.1.1 TCF: Recovery of Lithium from Geothermal Brine with Lithium-Aluminum Layered Double Hydroxide Chloride Sorbents

Lithium has gained importance as a near-critical element for various energy storage applications, and efforts are ongoing to develop simple and cost-effective methods for lithium recovery. For instance, in a typical 50 MW geothermal power plant in the Salton Sea, United States, approximately 15,000 metric tons of Li 2 CO 3 or LiOH salts can be produced annually by recovering and converting LiCl from geothermal plant waste solutions. Currently, approximately 390 MW of geothermal power is generated in the Salton Sea Known Geothermal Resource Area, which corresponds to potential Li 2 CO 3 production of approximately 120,000 metric tons per year. This makes lithium recovery from geothermal brines a promising, potentially low-cost source of lithium, contributing significantly to the global supply chain.

15 GEOTHERMAL ENERGY↗

Williston Basin Resource Study for Commercial-Scale Subsurface Hydrogen Storage

The Energy & Environmental Research Center (EERC), in partnership with the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL), the EERC’s State Energy Research Center (SERC), MPLX Operations LLC, and TC Energy Development Holdings Inc. (a subsidiary of TC Energy Corporation), studied the potential for subsurface hydrogen storage and recovery in the Williston Basin of western North Dakota. The project’s goal was to evaluate the feasibility of large-scale, secure geologic H 2 storage to support future hydrogen generation, storage, and use. This work included laboratory testing, H 2 –rock–fluid exposure experiments, literature reviews on H 2 embrittlement, and reservoir modeling and simulations. The study included an assessment of storage potential across three types of storage reservoirs using both reservoir simulation and DOE’s web-based tool SHASTA-HELP (Subsurface Hydrogen Assessment, Storage, and Technology Acceleration – Hydrogen Estimator for Logistical Planning), as well as investigation of potential H 2 production and markets for commercial-scale deployment. Building on prior EERC gas storage research, three storage options were selected for detailed evaluation: the Broom Creek Formation (a clastic saline reservoir), the Dickinson Lodgepole Mounds (DLM) complex (carbonate mud mound structures) of the Lodgepole Formation (an active oil and gas producing reservoir), and the Dunham Salt Interval of the Piper Formation (to be used for engineered salt cavern development). These targets were prioritized based on prior EERC research using datasets related to seal capacity, reservoir quality, mechanical integrity, and injectivity. Exposure tests on Broom Creek and DLM samples showed mineral dissolution and precipitation that increased brine salinity and altered reservoir rock surfaces. Although these results provide useful insight, they are limited by small sample sizes and short-term (30-day) exposure, requiring further study to assess long-term storage integrity. Salt formations were not tested because of their known nonreactivity and established mechanical stability. Results of reservoir simulations performed for a single site demonstrated that the Broom Creek Formation may be capable of receiving up to 42,000 tonnes of injected H 2 over 7 months via one well. H 2 recovery took place over 5 months, resulting in approximately 26,000 tonnes (~62% without cushion gas [CG]). This work suggests water production may be important and subsequent cycles of injection and production may perform more efficiently; however, significant site-specific work in the future is needed to assess actual reservoir performance of injection and withdrawal of H 2 storage. For oil reservoir potential, a multiple-well model was used to simulate injection of approximately 32,000 tonnes of H 2 into a single wellbore while simultaneously producing in place reservoir fluids from four offset wells to maintain reservoir pressure. The simulation results suggested a high recovery (~98%); in addition, cost advantages through existing infrastructure could be realized. Challenges in this reservoir include vi managing gas purity and leakage risks. In both scenarios, production of H 2 takes place in a single-well scenario with 10 cycles (7 months of injection and 5 months of production) over 10 years. Finally, the use of engineered caverns in the Dunhan Salt was evaluated, and the results suggest that while they have a smaller capacity (<1000 tonnes per cavern), they exhibit nearly complete gas recovery (>99%), fast response times, and low purity risk. While caverns in North Dakota may be smaller in capacity, fields can be developed in galleries to accommodate the volumetric needs and rapid turnaround times necessary to meet market demands. Geographic limitations and thin salt intervals in North Dakota may represent less total storage potential than salt domes elsewhere, but significant opportunities exist to expand this market for gas storage in North Dakota. A basinwide assessment was performed to estimate a first-of-its-kind value for H 2 storage on a large scale. DOE’s SHASTA-HELP, combined with EERC simulation work, was used to perform the assessment. Estimated H 2 storage potential varied widely for each formation type. The Broom Creek saline formation was estimated to have a storage potential of approximately 1.7–90.5 million tonnes (MMt). The DLM oil reservoirs were estimated to have 0.07–0.19 MMt of capacity. Notably, each of these estimates relies on significant assumptions regarding reservoir thickness, porosity, permeability, and CG needed for operation. Much research is needed to understand the true site-specific storage resource potential of each formation. Using the Dunham Salt Interval for cavern development may result in as much as 4.79 MMt (up to 2.87 MMt working gas) of H 2 storage potential. An important note for consideration is that the values presented here need significantly more geological characterization and engineering assessments prior to gaining confidence in performance. This will be a focal point for future research and development needs. The basinwide evaluation also indicated that North Dakota has significant H 2 generation potential, with estimates up to a possible 13 MMt annually, suggesting a substantial opportunity for H 2 market development and thus the need for commercial-scale H 2 storage to facilitate growth and resilience. Key Recommendations 1. Conduct detailed site characterization (3D seismic, well logs, core sampling) to reduce geologic uncertainty. 2. Perform techno-economic analyses incorporating market, regulatory, and incentive frameworks. 3. Investigate long-term interactions among H 2 , CGs, well materials, and formations to assess risks. 4. Develop pilot- and field-scale demonstrations to validate models and establish best practices.

03 NATURAL GAS↗

Mathematical methods for optimal polynomial recovery of high-dimensional systems from noisy data

The goal of our Early Career Research Project (ECRP) is to establish a modern mathematical foundation that will enable next-generation computational methods for polynomial approximation of high-dimensional systems, having a certain set of constraints, from a limited amount of noisy data. Such a foundation is critical to realizing the future potential of the DOE user facilities, and will ultimately empower scientists to address a fundamental question, namely, “how many realizations of a nonlinear manifold are required to recover the entire high-dimensional solution map, with optimal approximation guarantees and minimal computational cost?” The central theme of this effort aims to conquer this challenge by pioneering the development of extraordinarily innovative theoretical analysis and transformational non-intrusive computational methodologies. Such approaches will enable the reconstruction of the entire high-dimensional solution map, with accuracy comparable to the best approximation, while utilizing an optimal number of samples. During this reporting period we have made significant progress on four thrusts.

97 MATHEMATICS AND COMPUTING↗

Dual-loop Solvent-based CCS for Net Negative CO 2 Emissions with Lower Cost

This final technical report details the successful design, construction, and operational validation of an innovative dual-loop CO 2 capture technology designed to achieve deep decarbonization (>99%) from Natural Gas Combined Cycle (NGCC) power plants that results in the electricity with carbon intensity of approximate 42 kg CO 2 -eq/MWh, less than the electricity produced by solar PV. The integrated process couples a primary aqueous solvent (in this project, a water lean solvent – WLS) absorption loop for bulk CO 2 removal with a secondary potassium hydroxide (KOH) polishing loop featuring electrochemical regeneration. This architecture leverages the higher exergy efficiency of the primary loop while utilizing the fastest kinetic of the secondary loop to capture dilute residual CO 2 , achieving an overall capture rate of 99.9% and co-producing pure hydrogen after moisture being condensed and dehydrated. Technical feasibility was established through a comprehensive 2,000-hour experimental campaign on a bench-scale fully-integrated unit (using 4” absorber and 4” stripper) with the feeding flue gas flowrate in the range of 8-20 cfm, confirming the attainment of Technology Readiness Level (TRL) 4. The project executed extensive parametric testing followed by 1,000 hours of continuous steady-state testing, demonstrating exceptional process stability with the electrochemical regenerator exhibiting less than a 10% reduction in electrical conductivity over the duration of the campaign. Operational characterization gathering on the bench unit revealed distinct energy profiles for the hybrid system: the primary loop required approximately 280 kJ mol -1 for bulk removal, while the polishing loop required approximately 1,600 kJ mol -1 specifically when reducing dilute CO 2 concentrations from ~740 ppm down to <50 ppm. (Please note those energy values/numbers can only be viewed as relative relationship and should not be extrapolated as absolute values required for CO 2 capture). Furthermore, dynamic testing validated the system’s flexibility for utility applications, demonstrating a rapid process response time of <30 minutes to changes in flue gas flowrate. Emission monitoring confirmed that the dual-loop architecture effectively mitigates solvent losses, utilizing a water wash to remove entrained aerosols to <1 ppm. The Techno-Economic Analysis (TEA) indicates a cost of capture of $\$$59.3/tonne and a Levelized Cost of Electricity (LCOE) of 71.7 $\$$/MWh at the overall capture efficiency of 99.8% of total carbon in the flue gas stream, with sensitivity analysis identifying an economic optimum when the primary loop captures 97% of the total CO 2 . A Life Cycle Assessment (LCA) confirms the technology’s potential for net-negative emissions, determining a Global Warming Potential (GWP) of 52 kg CO 2 -eq/MWh—significantly lower than the baseline—which further decreases to 42 kg CO 2 -eq/MWh when crediting the displacement of conventional hydrogen production.

03 NATURAL GAS↗

Annual Summary Report (FY 2025) Performance Assessment for the Environmental Restoration Disposal Facility

Beginning in 1996, the Environmental Restoration Disposal Facility (ERDF) started accepting low-level radioactive, hazardous, and mixed wastes that were generated during cleanup activities at the Hanford Site. ERDF is composed of a series of cells or disposal areas and can accommodate future design expansions as needed. Currently, there are eight cells and two supercells in ERDF. Each supercell is the equivalent of two cells. During this reporting period, fiscal year (FY) 2025, extending from October 1, 2024, through September 30, 2025, approximately 50,653 metric tons (55,835 U.S. tons) of waste was disposed at ERDF. From ERDF inception through September 30, 2025, approximately 17.3 million metric tons (19.1 million U.S. tons) of waste has been disposed at ERDF, which equates to consumption of approximately 91.0% of the currently constructed disposal volume. According to the design of ERDF, the facility has the ability to be expanded as needed. In 2023, the Department of Energy approved the construction of supercell 11 to provide additional disposal capacity to meet the cleanup mission at the Hanford Site. The construction of supercell 11 began in FY 2025, and it is expected to take up to two years. Construction of the new cell will not affect ongoing waste disposal activities at ERDF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Charge And Dynamic Current On Tubular Antennas For Various Drive Conditions

The mixed boundary value problem of a tubular conductor is solved using an approximate representation of its Fourier coefficients. A two term solution is derived, which represents the solution over an extremely broad range of aspect ratios. This representation is used to find the electrostatic solution and capacitance of a charged tube as well as the solution of a tube in a uniform field and its dipole moment. This second case is directly useful as a model for a monopole electric field probe. This approximation is a special case of a representation using a combination of Chebyshev and Legendre polynomials. Combining the charged tube and tube in a uniform field allows the solution of voltage driven tubular antennas. Comparisons are made with numerical solutions using piecewise sinusoidal representations of the current. The results are also generalized to the dynamic case and up to and beyond the first resonance. Simple corrections for finite gap and delta gap drives to magnetic frill drives are examined using infinite tube integral transform representations. Corrections between magnetic frill drives and coaxial drives are also given. Approximate drive corrections using conformal mapping and an effective radius are also discussed. Finally, this efficient current representation is applied to the magnetic problem involving simple tubular solenoids.

97 MATHEMATICS AND COMPUTING↗

Storage Field Development Plan: One Earth Energy

This Storage Field Development plan presents the Storage Complex characterization results, construction, monitoring, and operational plans, and costs associated with the proposed One Earth Sequestration Carbon Capture and Storage (OES-CCS) site in McLean County, Illinois, near Gibson City. The proposed storage complex, known as the Mt. Simon Storage Complex, comprises the Cambrian Mt. Simon Sandstone reservoir and the primary seal, the Cambrian Eau Claire Formation. The lowermost Underground Source of Drinking Water (USDW) identified for the site is the Ordovician St. Peter Sandstone. Geologic characterization of the Mt. Simon Storage Complex at the OES-CCS site was performed by the Illinois Storage Corridor CarbonSAFE Phase III project, which also prepared and submitted three UIC Class VI applications to construct three injection wells; the permit applications were submitted and are in the federal EPA review process. A characterization well, OEE #1, was drilled to collect site-specific data. These data were analyzed and used to develop the UIC Class VI applications. The OEE #1 well will be converted to an in-zone monitoring (IZM) well for the injection phase. The proposed buildout for the OES-CCS site includes (1) three injection wells (OES #1, OES #2, and OES #3), (2) two IZM wells, (3) two above confining zone (ACZ) monitoring wells, one of which will be used to monitor the lowermost USDW, (4) capture and compression facilities, and (5) transportation facilities, i.e., pipelines. A pre-operational testing program was proposed in the Class VI permit application and will be employed at the site pending approval. Additional pre-injection (baseline), syn-injection, and post-injection monitoring and site care procedures will be followed by OES to ensure that injection activities are protective of human health and the environment. Injection is scheduled to begin in 2025, distributed across the three injection wells in accordance with the permit operating conditions. One Earth Sequestration intends to inject up to 90 million tonnes of CO 2 over a period of approximately 20 years. Injection will begin at approximately 0.5 million tonnes of CO 2 annually and ramp up to a maximum of 4.5 million tonnes annually. Daily injection rates are expected to range from 1,400 to 1,500 tonnes per day initially and reach a maximum of approximately 4,225 tonnes per day, depending on site geology and injectivity at each injection well location, and CO 2 availability. The costs associated with the OES-CCS project include pre-operational costs (e. g. additional seismic data acquisition and well drilling), capture and transportation facility and equipment costs, predicted field operating expenditures (OpEx), and decommissioning and post-injection site care (PISC) costs. The risks associated with project activities, such as site construction, injection operations, and verification of secure storage were evaluated, and mitigation strategies proposed to alleviate those risks.

09 BIOMASS FUELS↗

Non-Neutron Transmutation of Used Nuclear Fuel (Final Report)

The primary goal of this study is to develop a national facility concept for transmuting long-lived fission products (LLFP) to substantially reduce the disposal impact by minimizing the need for a geologic-timescale repository. As a charter for this study, the national transmutation facility was required to reduce the radiotoxicity and decay heat of LLFP isotopes by at least 90% relative to their values at discharge from a commercial LWR, while consuming less than 10% of the reactor's energy. The identified LLFP isotopes are Se-79, Zr-93, Tc-99, I-129, Sn-126, and Cs-135, whose radiotoxicity is about 99% of the total radiotoxicity of all fission products at 1,000 years. Approximately ~72 kg of LLFPs is discharged every year from a 1,000 MWe commercial or advanced nuclear reactor. First, LLFP transmutation options with non-neutron beams (photons and protons) were explored. The study concluded that LLFP transmutation is feasible with high-energy, high-intensity photons or protons, but impractical on an engineering scale due to low transmutation rates and the high energy requirements to produce the desired photon or proton beams. As alternatives, LLFP transmutation options with neutrons from fission, fusion, and spallation reactions were additionally explored. The transmutation options using advanced critical reactors are attractive only for selective LLFP isotopes because the production rates of several LLFP isotopes (Zr-93, Sn-126, and Cs-135) from fission reactions are larger than the transmutation rates. The transmutation options with only spallation neutrons are favorable to transmute all LLFP isotopes, but as a tradeoff, the net transmutation rates are reduced. The national transmutation facility concept was developed following an exploration of transmutation options using various incident particles. The proposed national LLFP transmutation comprises a dedicated molten-salt reactor (MSR), a proton accelerator, and a spallation neutron-based transmuter. The MSR power was set at 300 MWt and 120 MWe, with the thermal power approximately 10% of that of a commercial 1,000 MWe PWR. The electricity generated by the MSR powers the accelerator and transmuter. The accelerator produces 1 GeV, 30 mA protons, which are introduced into the spallation neutron-based transmuter. The spallation neutron-based transmuter consists of a central spallation target and LLFP target pins merged in a heavy water tank. The six LLFP isotopes are separated into two groups. Tc-99, I-129, and Se-79, having larger neutron cross sections, belong to group A, while Zr-93, Sn126, and Cs-135, having smaller neutron cross sections, belong to group B. Then, for effective transmutation, LLFPs in groups A and B are transmuted in the dedicated MSR and in a spallation neutron-based transmuter, respectively. The estimated capital cost of the national transmutation facility is approximately $\$$3.1B, and its annual O&M cost is expected to be ~$\$$182M. Radiotoxicity and decay heat of LLFPs were calculated and compared with those of the original LLFPs. It was assumed that the targets were made with elementwise LLFP rather than isotopic LLFP, owing to the potentially high cost of isotopic separation from used nuclear fuels. The decay heat of LLFPs can be reduced by more than 90% using a single national transmutation facility. However, radiotoxicity decreases by 79–84%, which does not meet the transmutation performance requirement, primarily because Cs-135 is produced rather than depleted. Thus, to meet the design requirement, Cs-135 should be separated from other Cs isotopes and irradiated in a spallation neutron-based transmuter. Then, radiotoxicity decreases by ~92%.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Surface and Buried Thermal, and RGB Unexploded Ordnance Data Collection

This document provides a description of a data collection campaign of unexploded ordnance (UXOI) set. The dataset captures a controlled UAV imaging campaign designed to support detection of UXO across varied environmental conditions. Data were collected during three campaigns in Norris and Northeast Knoxville, Tennessee, using RGB, and thermal sensors mounted on Parrot UKR. In total, the dataset contains 9925 images, 26 full-motion video, and approximately 81.99 GB of data, collected across late spring/summer conditions, every hour during sunlight, and multiple surface contexts, including tall grass, short grass, gravel, as well as buried in sand, and other gravel mixtures. The collection was designed to capture thermal and visual variability relevant to UXO detection in agricultural land, bare earth, and subsurface. Review of the imagery showed that ordnance was most detectable during periods of changing solar input, especially approximately 10-60 minutes after sunrise, approximately 20-60 minutes after sunset, and 2-3 min after cloud cover interrupted prolonged solar heating. These conditions increased thermal contrast because many ordnance items retained or released heat differently than the surrounding vegetation and ground surface. This dataset provides a useful resource for developing and evaluating airborne UXO detection methods under realistic field conditions. All ordnance used in the study was inert, and thermal behavior may differ from that of live ordnance. In addition, variation in ordnance type, composition, and placement introduced differences in thermal response that should be considered when interpreting results.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Thermal Reservoir Networks for Modularly Expandable Thermal Microgrids

The Department of Defense (DoD) faces the substantial challenge of cost-effectively retrofitting one to two installations per month, each comprising approximately 1,000 buildings, to improve resilience, reduce energy consumption, and enhance energy supply security. Achieving these objectives requires optimal system selection and effective risk mitigation during system integration. To address this need, we introduce Platform-Based Design (PBD), a structured, hierarchical methodology adapted from other industrial sectors to the domain of energy system retrofits. We demonstrate the effectiveness of PBD through a techno-economic feasibility study comparing geothermal-coupled thermal energy networks (TENs) with conventional energy systems for heating, cooling, and powering 17 buildings at Joint Base Andrews (JBA) in Maryland. Our analysis illustrates that the PBD approach enables rigorous, data-driven, sequential decision making, resulting in a family of Pareto-optimal systems, among which the TEN emerged as the most promising solution. The selected TEN design integrates geothermal borefields, heat recovery heat pumps, photovoltaic (PV) arrays, and battery storage. Compared to the baseline system – gas heating combined with air-source chillers – the proposed TEN reduces annual imported energy by 74% and peak electricity demand by 45%, achieves a levelized cost of energy of $\$0.210$/kWh, and substantially enhances resilience. Life-cycle costs increase by approximately 6%, and initial investment costs are about 2.5 times higher than the baseline. However, if central plant infrastructure, district loops, and utility-scale PV and battery systems are privately funded and operated, the initial investment would fall below the baseline system cost. Critical to achieving these significant performance improvements were detailed nonlinear dynamic simulations coupling geothermal heat transfer, energy system operation, and realistic feedback control logic. These simulations identified essential design modifications and control strategy refinements that substantially reduced energy use, peak demand, and compressor shortcycling, thereby improving durability and reliability—issues that would have been significantly more expensive to resolve during operation. Additionally, the verification step highlighted sensitivities to key design parameters that could reduce initial investment by approximately $\$2$ million and reduce annual life-cycle costs more than $\$300,000$. We recommend adopting the PBD methodology for future feasibility studies and TEN pilot projects to gain valuable operational experience. Furthermore, we recommend that DoD invest in transferring and scaling the PBD methodology to other installations. This entails developing standardized computational frameworks and component libraries as well as training industry in conducting PBD. Such investments would enable rapid, robust, reliable, and cost-effective retrofits, supporting DoD’s ambitious energy system modernization goals.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Integrative Modeling and Analysis of Fungal Central Carbon Metabolism

Over a thousand fungal genomes have been sequenced, yet manually curated genome-scale metabolic models (GEMs) are available for only a limited number of species. Moreover, these models have often been developed independently, leading to inconsistencies in namespaces, compartment definitions, and pathway representations that hinder comparative analysis, the systematic reuse of prior curation efforts, and the integration of consolidated metabolic knowledge. Here, we present the Consolidated Fungal Core Metabolism Model (CFCMM), constructed by integrating thirteen published fungal models spanning Ascomycota, Mucoromycota, and both Crabtree-positive and Crabtree-negative yeasts. We harmonized metabolites and reactions into a non-redundant shared ModelSEED ontological space, standardized compartmentalization, and refined gene–protein–reaction (GPR) rules. Using pathway-level visualization and systematic gap detection, we further improved the integrated network through literature-guided curation to correct stoichiometry, stereospecificity, and pathway architecture. Orthologous protein family reconstruction and functional annotation workflows were used to validate and inform GPR associations, with particular emphasis on ambiguous enzyme superfamilies and membrane-associated components. Using the resulting CFCMM, we built high-quality central carbon core models for each fungus and performed flux balance analysis to quantify ATP-yield variation under aerobic and anaerobic conditions, explicitly evaluating scenarios driven by differences in electron transport chain (ETC) composition. Simulations reproduced the expected fermentative yield of approximately 2 mmol ATP per mmol glucose under anaerobic conditions and separated the thirteen fungi into two bioenergetic groups under aerobic respiration based on Complex I status, with predicted yields of approximately 30 versus 22 mmol ATP per mmol glucose. Forcing flux through the alternative oxidase bypass further reduced ATP yields to approximately 12 and 4 mmol ATP per mmol glucose in Complex I-containing and Complex I-lacking fungi, respectively. Collectively, this work provides a manually curated, ModelSEED-consistent, and extensible fungal core metabolic template, deployed in DOE KBase as a resource for automated reconstruction of central carbon core models from any sequenced fungal genome. In addition, the CFCMM provides modular components for developing GEMs with more accurate energy predictions and enables robust comparative analyses of fungal bioenergetics and core metabolic diversity

59 BASIC BIOLOGICAL SCIENCES↗

Solid Wastes from Geothermal Energy Production and Implications for Direct Lithium Extraction

Direct lithium extraction (DLE) of brines after geothermal power production offers opportunities to produce environmentally benign “green” lithium; however, some environmental impact is inevitable. We examined solid waste production at geothermal power plants in southern California that are also locations for planned DLE facilities. Currently, the geothermal plants in this region produce approximately 79,800 metric tons (wet weight) per year of solid waste, which represents about 28 metric tons per GWh of net electricity production or approximately 500 mg solids per kg geothermal brine. Approximately 15% of this waste requires management as hazardous waste. Solids produced during power production represent about 0.2% of the total dissolved solids in the brine. Lithium production will require the removal of silica, iron, and other metals as part of the DLE process. Using a mass balance approach, we calculate that precipitation of silica and metals could produce up to an additional 6800 mg solids per kg brine. Calcium occurs at very high concentrations, and the amount of solids disposed in landfills will be dependent on the amount of calcium removed during lithium recovery. Our analysis shows that evaluation of brine chemistry in the context of the DLE process is useful for evaluating the potential solid waste impacts of producing lithium from brines.

42 ENGINEERING↗

The Vera C. Rubin Observatory Data Preview 1

We present Rubin Data Preview 1 (DP1), the first data from the National Science Foundation–Department of Energy Vera C. Rubin Observatory, comprising raw and calibrated single-epoch images, coadds, difference images, detection catalogs, and ancillary data products. DP1 is based on 1792 optical–near-infrared exposures acquired over 48 distinct nights by the Rubin Commissioning Camera (LSSTComCam) on the Simonyi Survey Telescope at the Summit Facility on Cerro Pachón, Chile in late 2024. DP1 covers ∼15 deg 2 distributed across seven roughly equal-sized noncontiguous fields, each independently observed in six broad photometric bands, ugrizy. The median FWHM of the point-spread function across all bands is approximately 1"14, with the sharpest images reaching about 0." 58. The 5σ point-source depths for coadded images in the deepest field, the Extended Chandra Deep Field South, are u = 24.55, g = 26.18, r = 25.96, i = 25.71, z = 25.07, and y = 23.1. Other fields are no more than 2.2 mag shallower in any band, where they have nonzero coverage. DP1 contains approximately 2.3 million distinct astrophysical objects, of which 1.6 million are extended in at least one band in coadds, and 431 solar system objects, of which 93 are new discoveries. DP1 is approximately 3.5 TB in size and is available to Vera C. Rubin Observatory data rights holders via the Rubin Science Platform, a cloud-based environment for the analysis of petascale astronomical data. While small compared to future LSST releases, its high quality and diversity of data support a broad range of early science investigations ahead of full operations in 2026.

Ground-based astronomy↗

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗

Multiqubit entanglement generation with squeezed modes

We present a hybrid continuous variable–discrete variable entanglement generation protocol using linear optics and homodyne measurements, capable of producing multiple high-fidelity Bell pairs per protocol iteration, with an approximate 0.5 success probability. The effectiveness of the protocol is determined by the squeezing strength. To increase the number of Bell pairs, approximately 3 dB of extra squeezing is needed for each additional Bell pair. The protocol also generates single Bell pairs with an approximate 0.75 probability for squeezing strengths ⪅ 15 dB , achievable with current technology.

Macridin, Alexandru [Fermilab] (ORCID:000000022228↗

Algorithm to extract direction in 2D discrete distributions and a continuous Frobenius norm

In this study, we present a novel algorithm for determining directionality in 2D distributions of discrete data. We compare a reference dataset with a known direction to a measured dataset with an unknown direction by the Frobenius norm of the difference (FND) to find the unknown direction. To generalize this concept, we develop a continuous Frobenius norm of the difference (CFND) as a continuous analog of the FND and derive its analytical expression. By relating fitted and normalized 2D Gaussian distributions, we show that the CFND approximates the FND, and we validate this relationship with computer simulations. We find that a first-order approximation of the CFND between two similar Gaussian distributions takes the form of an absolute sine function, offering a simple analytical form with potential for specialized applications in segmented inverse beta decay (IBD) neutrino detectors, astronomy, machine learning, and more. Although this method may easily extend to 3D scalar fields, our focus here is on 2D real-valued fields as it directly applies to directionality. Our methodology consists of modeling a 2D Gaussian distribution, binning the data into a histogram, and encoding it as a square matrix. Rotating this matrix around its geometric center and comparing it to a measured dataset using the FND gives us rotational data that we fit with an absolute sine function. The location of the minimum of this fit is the angle closest to the true angle of the direction in the measured dataset. We present the derivation and discuss initial applications of the CFND in our novel algorithm, demonstrating its success in approximating directionality in 2D distributions.

Data Analysis, Statistics and Probability (physics↗