Search NASA⌕ Search

SEARCH · Search NASA

Results for “Energy usage prediction”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Battery asset management with cycle life prognosis

We report Battery Asset Management problem determines the minimum cost replacement schedules for each individual asset in a group of battery assets that operate in parallel. Battery cycle life varies under different operating conditions including temperature, depth of discharge (DOD), charge rate, etc., and a battery deteriorates due to usage, which cannot be handled by current asset management models. This paper presents a new battery asset management methodology where battery cycle life prognosis is integrated with parallel asset management to reduce lifecycle cost of the Battery Energy Storage Systems (BESS). For the battery failure time prognosis, a nonlinear physics-based battery capacity fade model is developed and incorporated in parallel asset management model to update battery capacity over time. Experiment results have shown that the developed battery asset management methodology can be conveniently used to facilitate BESS asset management decision making thereby decreasing asset lifecycle costs.

25 ENERGY STORAGE↗
DOE OSTI2024

Investigating explainable transfer learning for battery lifetime prediction under state transitions

Battery lifetime prediction at early cycles is crucial for researchers and manufacturers to examine product quality and promote technology development. Machine learning has been widely utilized to construct data-driven solutions for high-accuracy predictions. However, the internal mechanisms of batteries are sensitive to many factors, such as charging/discharging protocols, manufacturing/storage conditions, and usage patterns. These factors will induce state transitions, thereby decreasing the prediction accuracy of data-driven approaches. Transfer learning is a promising technique that overcomes this difficulty and achieves accurate predictions by jointly utilizing information from various sources. Hence, we develop two transfer learning methods, Bayesian Model Fusion and Weighted Orthogonal Matching Pursuit, to strategically combine prior knowledge with limited information from the target dataset to achieve superior prediction performance. From our results, our transfer learning methods reduce root-mean-squared error by 41% through adapting to the target domain. Furthermore, the transfer learning strategies identify the variations of impactful features across different sets of batteries and therefore disentangle the battery degradation mechanisms and the root cause of state transitions from the perspective of data mining. These findings suggest that the transfer learning strategies proposed in our work are capable of acquiring knowledge across multiple data sources for solving specialized issues.

25 ENERGY STORAGE↗
DOE OSTI2024

Association Kinetics for Perfluorinated n -Alkyl Radicals

Radical-radical reaction channels are important in the pyrolysis and oxidation chemistry of perfluoroalkyl substances (PFAS). In particular, unimolecular dissociation reactions within unbranched n-perfluoroalkyl chains, and their corresponding reverse barrierless association reactions, are expected to be significant contributors to the gas-phase thermal decomposition of families of species such as perfluorinated carboxylic acids and perfluorinated sulfonic acids. Unfortunately, experimental data for these reactions are scarce and uncertain. Furthermore, obtaining reliable theoretical predictions for such reactions is a laborious and computationally intensive task. Here, in this work, the chemical kinetics of the various association/decomposition reactions producing/decomposing the C 2 -C 4 series of unbranched n-perfluoroalkanes (C 2 F 6 , C 3 F 8 , and C 4 F 10 ) are examined using state-of-the-art ab initio transition-state-theory-based master-equation calculations. The variable-reaction-coordinate transition-state theory (VRC-TST) formalism is employed in computing the microcanonical and canonical rates for the association reactions. Reaction thermochemistry is obtained via composite quantum chemistry calculations and the laddering of error-canceling reaction schemes via a connectivity-based hierarchy approach employing ANL1/ANL0-style reference energies. Lennard-Jones collision model parameters for the considered systems were estimated by a direct dynamics approach, and collisional energy transfer parameters were obtained from analogies to systems of similar size and heavy-atom connectivity. A one-dimensional master equation approach was used to convert the microcanonical rate coefficients from the VRC-TST analysis into temperature- and pressure-dependent rate constants for the association reactions and the reverse dissociation reactions. The data are reported in standardized formats for usage in comprehensive chemical kinetic models for PFAS thermal destruction.

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

Porosity, swelling, and composition evolution in high-burnup monolithic U-Mo fuel

The microstructural progression of very high-burnup (>8 × 10 21 fissions/cm 3 ) monolithic uranium-molybdenum (U-10wt.%Mo) was analyzed, providing crucial insights into the behavior of post-recrystallized nuclear fuel, where scant data exists. Three focused ion beam cuboids sourced from a fuel plate with varying local burnups of 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 were characterized. The porosity and composition of the samples were evaluated to characterize the evolution of the microstructure as a function of fission density and different locations on the fuel plate, while simultaneously isolating plate-specific parameters such as Zr diffusion barrier thickness, hot-isostatic-pressing conditions, enrichment, and reactor conditions. The porosity was segmented, and the three-dimensional distribution of the porosity was extrapolated from the two-dimensional segmentation. The composition was assessed and quantified using energy-dispersive X-ray spectroscopy areal mapping. The porosity fraction increased as a function of the burnup from 27.77±0.51, 35.12±1.54, and 37.71±0.44 % for 8.86 × 10 21 , 9.05 × 10 21 , and 9.36 × 10 21 fissions/cm 3 , respectively. When compared to literature, the porosity volume fraction plateaus at burnups greater than 6 × 10 21 fissions/cm 3 , while the pore size grows linearly as a function of fission density. The number of large pores increased in number density as a function of burnup, while the smallest pores (<0.3 µm) increased up to 9.05 × 10 21 fissions/cm 3 , followed by a decrease at 9.36 × 10 21 fissions/cm 3 . The delamination and cracking in the fuel plate propagated through an interconnected porosity sublayer identified ∼5 µm from the diffusion barrier. The local swelling of the specimens was within or near the prediction bounds of the Robinson-Williams model for local swelling. The fission products, strontium, barium, cerium, and cesium, precipitated into the pores, while neodymium accumulated adjacent to the pores. Furthermore, these findings have direct implications for the development of fuel performance codes and the accurate documentation of the microstructure evolution in high burnup U-Mo, thus enhancing the safety and efficiency of nuclear fuel usage.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Regional Medium-Term Hourly Electricity Demand Forecasting Based on LSTM

This paper aims to forecast high-resolution (hourly) aggregated load for a certain region in the medium term (a few days to over a year). One region is defined as some places with similar climate characteristics because the climate influences people's daily lifestyles and hence the electric usage. We decompose the electric usage records into two parts: base load and seasonal load. Considering both temperature and time factors, different deep-learning methods are adopted to characterize them. The first goal of our approach is to predict the peak load which is critical for power system planning. Furthermore, our proposed forecast method can provide the depiction of the hourly load profile to provide customized load curves for high-level real-time applications. The proposed method is tested on real-world historical data collected by CAISO, BPA, and PACW. The experimental results show that trained by three years of data, our method could reduce the prediction error for a one-year lead hourly load below $5\%$ MAPE, and predict the occurrence of the peak load for next year in CAISO with an error within three days. Furthermore, as a byproduct, an interesting observation on the impact of COVID-19 on human life was made and discussed based on these case studies.

deep learning↗
DOE OSTI2021

Timed to save: the added value of accounting for hourly incidence of electricity savings from residential space-conditioning measures

Abstract Previous research has recognized that the value of measures that reduce electricity usage depends upon the timing of the savings generated, but the lack of hourly savings shapes has limited the demonstration of this concept. We develop empirical hourly savings shapes for residential space-conditioning measures from nearly 18,000 efficiency projects in California and show how they combine with the diurnal and seasonal variation in electricity system costs. We find that these measures (cooling replacements; windows, doors, and skylights; and other envelope measures) tend to save electricity when system costs are highest and that the hourly savings account for 1.4–1.5 times as much value as non-time-sensitive estimates of efficiency would predict. We present these impact multipliers for each measure to quantify the additional value revealed by the time-sensitive approach. We show that this additional value is similar in an evolving electricity grid with storage, rather than natural gas generation, as the marginal resource.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
NASA NTRSDate not supplied

Runway Configuration Management with Offline Reinforcement Learning

Runway configuration management (RCM) is a challenging task, and it affects the efficiency of the National Airspace System (NAS) and airport surface operations significantly. Each airport, depending on the geometry, capacity, local climate patterns, etc. has multiple configurations for the runway usage for arriving and departing flights. Many factors such as the incoming/outgoing traffic load, wind direction and speed, convective weather, cloud ceiling and other environmental factors might affect the choice of a runway configuration at any point in time. However, other factors such as safety measures and regulations, noise abatement, capacity of each configuration, and preference of the air traffic controllers (ATCs) can also play a significant role in selecting the configuration. A sub-optimal selection of the runway configuration, or delay in making configuration changes might result in significant increase in taxi times for aircraft on the surface of the airport, fuel and energy use of the aircraft, and maintenance costs. It can also lead to safety concerns, such as an aircraft performing one or more go-arounds before being able to land. All these factors make RCM an extremely important and challenging decision-making process for the ATCs. The current state of practice sets the runway configuration by the ATCs based on relevant information available at the time including weather, traffic, noise abatement, safety bounds, etc. This makes the decision-making process subjective based on the accuracy of the available information and the bias in human decision making. Unfortunately, this approach yields poor results (e.g., significant delays) if the predicted outcomes are uncertain and their relative impact is not well understood. This is especially evident when the uncertainty increases the size of possible predicted outcomes (combinatorial explosion in possible scenarios) that cannot be handled by human reasoning. On the other hand, an automated approach based on machine intelligence can make use of historical data and search through all (or significant amount of) possible scenarios under uncertainty and make well-informed decisions.

Milad Memarzadeh↗
DOE OSTI2021

Bamboo Bio-composite Truck/Trailer Decking

Bamboo is one of the fastest growing plants in the world. Massively productive, bamboo will maintain that productivity with limited inputs and minimal management resulting in predictable volume and operating margins and with low overhead. The goal of this project was to develop a bamboo bio-composite trailer decking product that replaces apitong, is lighter weight, stronger and luminescent. The collaboration partners, Fontaine, Resource Fiber and IACMI, worked closely to develop the bamboo bio-composite decking to Fontaine Trailer’s (end-user) specifications. The various technological objectives were (a) To conduct a comprehensive Design of Experiments study for bamboo composites to understand the structure-property relationships for different resin composites and bamboo forms (strips, woven, bulk etc); (b) To establish the nail pull out characteristics of trailer deck geometry bamboo composite form. The nail pull out is a critical test in trailer decks; (c) To design, process and prototype select number of trailer decking planks for Fontaine evaluation; and (d) To conduct life cycle analysis to conduct energy calculations from the various conversion steps of the bamboo from crop to product. The project meets DOE/IACMI metrics of reduced embodied energy, lightweighting and lowering the cost of the end product. It also offers a green solution to a value-added application, i.e. trailer decking. The development goals and benefits of bamboo composite decking included: A. Decreased weight when compared to Apitong and Gen 1, thereby reducing operating costs and petroleum usage due to increased miles per gallon. B. Flexibility to allow camber design into trailers. C. Luminescence to increase safety during low light and dark conditions. D. Composite material to maximize use of bamboo bio-composites as substitute for petroleum-based composites. E. Embedded layer(s) of conditioned bamboo for added strength and stability. F. Improved safety over Apitong by temporary cargo indentation in planks for better stability. G. Cradle-to-cradle design so end-of-life becomes beginning-of-life for other products. Resource Fiber’s bamboo biocomposite trailer decking passed nail pullout tests as compared to apitong (incumbent), was lighter weight than apitong, and used less embodied energy particularly when bamboo is sourced domestically. The planks were optimized in the lab setting to the extent possible and successfully installed on the trailer at Fontaine for field testing. Despite the process not being fully optimized due to through-heat and tooling limitations within the budget constraints, the prototype decking withstood 400 cycles of reverse fatigue loading of a 107,000 lb Caterpillar 349F which was extreme conditions of field testing by Fontaine Trailer. Resource Fiber plans to outsource production of a next round of planks with a commercial pultrusion processor, then to do a re-test with Fontaine Trailers. Tooling specific to the part is required. Long-term commercial plans are to continue outsourcing production while supplying bamboo fiber and mats. Commercial markets include heavy haul trailers, military trailers, and decking for marine, industrial and residential use.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Efficient light-trapping in ultrathin GaAs solar cells using quasi-random photonic crystals

Ultrathin solar cells reduce material usage and allow the use of lower-quality materials thanks to their one order of magnitude smaller thickness than their conventional counterparts. However, efficient photonic light-trapping is required to harvest the incident light efficiently for an otherwise insufficient absorber thickness. Quasi-random photonic crystals are predicted to have high efficient light-trapping while being more robust under angle and thickness variations than simple photonic crystals. In this study, we experimentally demonstrate a light-trapping solution based on quasi-random photonic crystals fabricated by polymer blend lithography. We control the average lattice parameter by modifying the spin-coating speed. We demonstrate an ultrathin GaAs cell of 260 nm with a rear quasi-random pattern with submicron features, and a J sc = 26.4 mA/cm 2 and an efficiency of 22.35% under the global solar spectrum.

14 SOLAR ENERGY↗
NASA NTRS1999

A Study of Direct and Cloud-Mediated Radiative Forcing of Climate Due to Aerosols

The Intergovernmental Panel on Climate Change (IPCC) has reported that in the southeastern US and eastern China, the general greenhouse warming due to anthropogenic gaseous emissions is dominated by the cooling effect of anthropogenic aerosols. To verify this model prediction in eastern China and southeastern US, we analyzed regional patterns of climate changes at 72 stations in eastern China during 1951- 94 (44 years), and at 52 stations in the southeastern US during 1949-94 (46 years) to detect the fingerprint of aerosol radiative forcing. It was found that the mean rates of change of annual mean daily, maximum, minimum temperatures and diurnal temperature range (DTR) in eastern China were 0.8, -0.2, 1.8, and -2.0 C/100 years respectively, while the mean rates of change of annual mean daily, maximum, minimum temperatures and DTR in the southeastern US were -0.2, -0.6, 0.2, and -0.8 C/100 years, respectively. This indicates that the high rate of increase in annual mean minimum temperature in eastern China results in a slightly warming trend of daily temperature, while the high rate of decrease in annual mean maximum temperature and low rate of increase in annual mean minimum temperature lead to the cooling trend of daily temperature in the southeastern US. We found that the warming from the longwave forcing due to both greenhouse gases and aerosols was completely counteracted by the shortwave aerosol forcing in the southeastern US in the past 46 years. A slightly overall warming trend in eastern China is evident; winters have become milder. This finding is explained by hypothesizing that increasing energy usage during the past 44 years has resulted in more coal and biomass burning, thus increasing the emission of absorbing soot and organic aerosols in eastern China. Such emissions, in addition to well-known Asia dust and greenhouse gases, may be responsible for the winter warming trend in eastern China that we have reported here. The sensitivity of aerosol radiative properties to aerosol composition, size distribution, relative humidity (RH) is examined for the following aerosol systems: inorganic and organic ions (Cl-, Br-, NO3 -, SO4 2-, Na+, NH4 +, K+, Ca2+, Mg2+, HCOO-, CH3COO-, CH3CH2COO-, CH3COCOO-, OOCCOO2-, MSA-1); water-insoluble inorganic and organic compounds (elemental carbon, n-alkanes, SiO2, Al2O3, Fe2O3 and other organic compounds). The partial molar refraction method was used to calculate the real part of the refractive index. It was found that the asymmetry factor increased by approximately 48% with the real part varying from 1.40 to 1.65, and the single scattering albedo decreased by 24% with the imaginary part varying from -0.005 to -0.1. The asymmetry factor increased by 5.4 times with the geometric standard deviation varying from 1.2 to 3.0. The radiation transmission is very sensitive to the change in size distribution; other factors are not as significant. To determine the aerosol direct radiative forcing (ADRF), the aerosol optical depth (AOD) values at the three operational wavelengths (415, 500 and 673 nm) were determined at a regionally representative site, namely, Mt. Gibbs (35.78 deg N, 82.29 deg W, elevation 2006 m) in Mt. Mitchell State Park, NC, and a site located in an adjacent valley (Black Mountain, 35.66 deg N, 82.38 deg W, elevation 951 m) in the southeastern US. The two sites are separated horizontally by 10 km and vertically by 1 km. It was found that the representative total AOD values at 500 nm at the valley site for highly polluted (HP), marine (M) and continental (C) air masses were 0.68 +/- 0.33, 0.29 +/- 0.19 and 0.10 +/- 0.04, respectively. A search-graph method was used to retrieve the columnar size distribution (number concentration N, effective radius reff and geometric standard deviation=?g) from the optical depth observations at three operational wavelengths. The ground albedo, single scattering albedo and imaginary part of the refractive index were calculated using a mathematically unique procedure involving a Mie code and a radiative transfer code in conjunction with the retrieved aerosol size distribution, AOD, and diffuse-direct irradiance ratio. It was found that N, r(eff) and sigma(g) were in the ranges of 10 to 1.7 x 10(exp 4)/cubic cm, 0.09 to 0.68 micrometers and 1.12 to 2.95, respectively. The asymmetry factor and single scattering albedo were in the ranges of 0.63 to 0.75 and 0.74 to 0.97 respectively. The ground albedo for the forested terrain and imaginary part of refractive index were found to be in the ranges of 0.06 to 0.29 and 0.005 to 0.051 respectively. On the basis of these aerosol radiative properties obtained at the research sites and computations using the Column Radiation Model (CRM) of National Center of Atmospheric Research (NCAR) Community Climate Model (CCM3), it was found that the average cloud-free 24-hour ADRF values were -13 +/- 8, -8 +/- 3, -33 +/- 16 W/square m for marine, continental, and polluted air masses, respectively. On the assumption that the fractional coverage of clouds is 0.61, it was estimated that the annual mean ADRF was 7 +/- 2 W/square m in the southeastern US. The review with respect to the current knowledge of organic acids shows that aerosol formate and acetate concentrations range from 0.02 to 5.3 nmol/cubic m and from 0.03 to 12.4 nmol/cubic m respectively, and that between 34% to 77% of formate and between 21% to 66% of acetate are present in the fine fraction of aerosols. It was found that although most (98-99%) of these volatile organic acids were present in the gas phase, their concentrations in the aerosol particles were sufficient to make them a good candidate for cloud condensation nuclei (CCN). It is hypothesized that organic acids are at least one of the primary sources of CCN in the atmosphere due to their ubiquitous presence in the troposphere, especially over the continental forested areas. The results of our measurements at Palmer Station, Antarctica show that the daily average CCN concentrations at 0.3% and 1% supersaturations ranged from 0.3 to 160/cubic cm and from 4 to 168/cubic cm, respectively, during the period from 17 January to 26 February, 1994. New evidence for substantial and definitive CCN enhancement near and within cloud has been observed at Mt. Mitchell, North Carolina. The results show that the average monthly CCN concentrations were 460 +/- 217, 386 +/- 286, 429 +/- 228 and 238 +/- 134/cubic cm for in-cloud, overcast, clear and rainy conditions, respectively. The typical CCN spectra show that there were a lot of small CCN produced and the ion concentrations (especially H+ and SO4 2-) were very high during the CCN enhancement period. The significantly positive correlation between black carbon (BC) and CCN at 1% supersaturation indicates that a percentage of the BC measured at the site may be in the form of an internal mixture and participated in the formation of CCN.

Yu, Shao-Cai↗
DOE OSTI2022

Development of simplified model for injection rate prediction of diesel injectors during transient and steady operation

Determining the injection rate of the diesel injector during the transient short injection operations is a critical matter to enhance the combustion process by improving air/fuel mixture formation. In this regard, numerous injection models have been suggested in the literature to predict the injection rate accurately. However, such algorithms have not precisely predicted the rate of injection (ROI) in transient short injection events. In addition, most of those models require many input parameters to characterize hydraulic and mechanical subsystems in the injector, which can make them inconvenient to tune for different injector specifications, in other words not suitable for universal usage. Therefore, our study aims to develop a simple injection model that can predict injection rate precisely in the short and long injection durations considering relatively fewer input parameters. The model traces the pressure variation inside the sac based on the bulk modulus theory to predict nozzle exit velocity, thus the injection rate of the injector. In the algorithm, the bulk modulus of fuel (E) and discharge coefficients of the nozzle hole and needle seat (C d,H and C d,N ) have been considered to vary rather than set constant during the injection process since such flow parameters would have a critical impact on the flow development at the transient stage, thus affecting the injection rate of diesel injectors. Here with this modelling strategy, the accuracy of the new model on ROI prediction was confirmed with acceptable model error for both short and long injection durations by validating against experimental results. To understand the effect of the new model scheme on ROI prediction, ROI results with the new model have been discussed by comparing it with the ROI results predicted with constant flow parameters.

42 ENGINEERING↗
DOE OSTI2025

Integration of a grey-box refrigerated case model in EnergyPlus via Python plugin

Commercial buildings, in particular grocery stores (due mainly to their large refrigeration load), provide opportunities for energy cost reductions. Grocery stores could offer substantial load flexibility to the power grid through participation in demand response programs because of their usage patterns and relatively high energy intensity. This load flexibility could come from modifying the control of heating, ventilation, and air conditioning (HVAC) systems, refrigeration systems, or both. Although estimation of the HVAC system’s load flexibility potential is relatively targeted in the literature, estimating load flexibility of refrigeration systems is nascent and has been a challenge, in part because of the lack of proper simulation tools that capture the dynamics in the refrigeration cases. The existing refrigerated case model within EnergyPlus, a whole building energy simulation program, assumes a constant case temperature throughout the simulation period and does not explicitly model the cycling of the compressor serving the refrigerated case. In addition, it does not encompass modeling of temperatures of the product inside the refrigerated case. This difference between modeled and actual operation can be a barrier to the development of demand control algorithm and accurate analysis of load flexibility potential. In this paper, we present a grey-box model for modeling refrigerated cases in grocery stores, which include medium temperature and low temperature. Four cases are modeled; two are low-temperature closed cases and two are medium-temperature cases with one closed and one open. Data from an experimental facility are used to train and test the models. Results demonstrate the efficacy of the grey-box models in predicting the temperatures. This model is integrated into EnergyPlus to capture the dynamic effects of case temperature on the environment and enhance the calculation of sensible and latent heat exchange with the environment (case credits). These enhancements can be leveraged more broadly to model advanced refrigeration controls such as defrost, develop and test unique algorithms that could affect refrigeration interactions with HVAC, and refine store design for any commercial building with refrigeration.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

Decarbonization of heat pump dual fuel systems using a practical model predictive control: Field demonstration in a small commercial building

In the transition from fossil fuel to electrified heating, a concerning trend is emerging in certain regions of the US. Owners of buildings with gas-based systems leave them in place after adding heat pumps (HPs). Existing control solutions for these hybrid (dual fuel) systems are rudimentary and fall short of realizing the full carbon reduction potential of these systems. Model predictive control (MPC) is often regarded as the benchmark for achieving optimal control in integrated systems. However, in the case of small-medium commercial buildings (SMCBs), the control and communication infrastructure required to facilitate the implementation of such advanced controls is often lacking. This paper presents a field implementation of easy-to-deploy MPC for a dual fuel heating system consisting of HPs and a gas-fired furnace (GF) for SMCBs. The control system is deployed on an open-source middleware platform and utilizes low-cost sensor devices to be used for real SMCBs without major retrofits. Here, we demonstrated this MPC in a real office building with 5 HPs and 1 GF for 2 months. The test results showed that MPC reduced 27% of cost while completely eliminating GF usage by shifting 23% of the thermal load from occupied-peak time to non-occupied-non-peak times.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2024

A multiphysics coupling framework for exascale simulation of fracture evolution in subsurface energy applications

Predicting the evolution of fractured media is challenging due to coupled thermal, hydrological, chemical and mechanical processes that occur over a broad range of spatial scales, from the microscopic pore scale to field scale. We present a software framework and scientific workflow that couples the pore scale flow and reactive transport simulator Chombo-Crunch with the field scale geomechanics solver in GEOS to simulate fracture evolution in subsurface fluid-rock systems. This new multiphysics coupling capability comprises several novel features. An HDF5 data schema for coupling fracture positions between the two codes is employed and leverages the coarse resolution of the GEOS mechanics solver which limits the size of data coupled, and is, thus, not taxed by data resulting from the high resolution pore scale Chombo-Crunch solver. The coupling framework requires tracking of both before and after coarse nodal positions in GEOS as well as the resolved embedded boundary in Chombo-Crunch. We accomplished this by developing an approach to geometry generation that tracks the fracture interface between the two different methodologies. The GEOS quadrilateral mesh is converted to triangles which are organized into bins and an accessible tree structure; the nodes are then mapped to the Chombo representation using a continuous signed distance function that determines locations inside, on and outside of the fracture boundary. The GEOS positions are retained in memory on the Chombo-Crunch side of the coupling. The time stepping cadence for coupled multiphysics processes of flow, transport, reactions and mechanics is stable and demonstrates temporal reach to experimental time scales. The approach is validated by demonstration of 9 days of simulated time of a core flood experiment with fracture aperture evolution due to invasion of carbonated brine in wellbore-cement and sandstone. We also demonstrate usage of exascale computing resources by simulating a high resolution version of the validation problem on OLCF Frontier.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Development of Numerical Model of Metal Foam with PCM for the Estimation of Effective Thermal Conductivity

Global warming due to climate change is a threat to humankind. Nuclear energy is one of the promising solutions to reduce fossil fuel usage. Nuclear energy can handle the base load, compensating for the volatility of renewable energy. If nuclear energy could achieve load following capability, the combination with renewable energy would be more suitable. Thermal energy storage (TES) is one of the options for enabling load following of nuclear reactors. The TES makes it possible to store surplus nuclear thermal energy and release it later as needed. In Idaho National Laboratory (INL), a new concept of latent heat TES integrated with high-temperature heat pipe has been proposed and is under development, which is called Heat pipe-Integrated Thermal Battery (HITB). HITB exchanges thermal energy between the reactor system and TES via heat pipe. The heat transferred to TES medium, made of phase change material (PCM), stores energy as sensible heat and/or latent heat. As PCM typically has poor thermal conductivity, however, various heat transfer enhancement techniques are required to achieve a rapid charging cycle. There are many techniques to enhance the heat transfer ability of TES medium such as disk, fin, and metal foam. Among them, metal foam is an appropriate option to enhance the heat transfer because it maximizes the heat transfer area through metal wicks. Metal foam is a lightweight metal structure that has a high porosity of over 0.9. The typical materials for metal foam are Aluminum, Copper, Nickel, and Silicon Carbide (SiC). Metal foam not only enhances heat transfer via conduction but also increases contact surface area. In the HITB design , the metal foam is being considered as one of the options to enhance the heat transfer of TES medium (PCM) [1]. To predict the enhanced thermal performance of TES, one should properly estimate the effective thermal conductivity of metal foam combined with PCM material or calculate heat transfer in distributed model. There are many experimental works that provides effective thermal conductivity of metal foam with various PCM [2,3]. Also, many theoretical models were developed based on the unit cell model of metal foam [4,5]. With a distributed model, on the other hand, detail heat transfer characteristics between metal foam and PCM material can be analyzed considering the geometry or buoyancy effect. However, due to the complex geometry of metal foam pores, the computational cost for three-dimensional modeling highly increases. Therefore, if metal foam structure can be modeled in simple and repetitive design, the computational cost would decrease Among the various metal foam models [2], lattice model is one of the simple and extendable design. The porosity and pores per inch (PPI) can be characterized by the size and spatial distance of lattice structure. If the three-dimensional metal foam model consists of lattice structure could properly estimate the heat transfer, which is characterized by effective thermal conductivity, it would be a good option to assess the thermal performance of metal foam with PCM. In this study, a three-dimensional numerical model was developed to simulate conductive heat transfer between metal foam and PCM. The three-dimensional lattice structure of square pillars was selected as a basic structure of the metal foam. The calculation result was characterized by the effective thermal conductivity of the whole domain. A sensitivity study was conducted for mesh size, domain size, and PPI to check whether the calculation result gives a converged result or not. Lastly, the effective thermal conductivity from the lattice model was compared with existing experimental data to validate the model result

25 ENERGY STORAGE↗
DOE OSTI2022

Resolving Current-Dependent Regimes of Electroplating Mechanisms for Fast Charging Lithium Metal Anodes

Poor fast-charge capabilities limit the usage of rechargeable Li metal anodes. Understanding the connection between charging rate, electroplating mechanism, and Li morphology could enable fast-charging solutions. In this report we develop a combined electroanalytical and nanoscale characterization approach to resolve the current-dependent regimes of Li plating mechanisms and morphology. Measurement of Li + transport through the solid electrolyte interphase (SEI) shows that low currents induce plating at buried Li||SEI interfaces, but high currents initiate SEI-breakdown and plating at fresh Li||electrolyte interfaces. The latter pathway can induce uniform growth of {110}-faceted Li at extremely high currents, suggesting ion-transport limitations alone are insufficient to predict Li morphology. At battery relevant fast-charging rates, SEI-breakdown above a critical current density produces detrimental morphology and poor cyclability. Thus, prevention of both SEI-breakdown and slow ion-transport in the electrolyte is essential. This mechanistic insight can inform further electrolyte engineering and customization of fast-charging protocols for Li metal batteries.

25 ENERGY STORAGE↗
DOE OSTI2022

Manning condensation in ion exchange membranes: A review on ion partitioning and diffusion models

Here, the rational design of ion exchange membranes (IEMs) is becoming more pertinent as their usage becomes broader and as their staple applications (i.e., electrodialysis, flow batteries, and fuel cells) improve in commercial viability. Such efforts would be catalyzed by an improved fundamental understanding of ion transport in IEMs. This review discusses recent progress in modeling ion partitioning and diffusion in IEMs in an effort to relate IEM performance metrics to fundamental membrane properties over which researchers and membrane manufacturers possess direct and sometimes precise control. Central focus is given to the Donnan-Manning model for ion partitioning and the Manning-Meares model for ion diffusion in IEMs. These two frameworks, which are derived from Manning's counter-ion condensation theory for polyelectrolyte solutions, have been widely used within the IEM literature since their recent introduction. To explore this topic, the mathematical derivation of both models is revisited, followed by a survey of experimental and computational discussions of counter-ion condensation in IEMs. Alternative models which fulfill similar roles in predicting IEM transport properties are compared. This review concludes by highlighting the uniquely favorable positions of the Donnan-Manning and Manning-Meares models and discussing their prospects as leading predictors of IEM partitioning and diffusive properties.

25 ENERGY STORAGE↗
DOE OSTI2024

Use of tritium-rich fuel to improve the yield of layered deuterium/tritium inertial fusion capsules

In deuterium–tritium (DT) ice layered implosions, nearly all hot spot mass at peak burn comes from the dense fuel. Accurate prediction of the fuel mass ablation, including the enthalpy associated with mass inflow into the hot spot from the dense fuel, is essential to understanding the energetics and ignition of the hot spot in layered implosions. A recently published boundary layer analysis (Daughton et al., 2023) indicates a faster mass ablation rate than in previous analyses of layered implosions. Inclusion of this effect provides a better match to simulations and leads to a new ignition threshold where the temperature of the dense fuel plays a critical role. This analysis motivates possible new directions for improved capsule performance. Here, the authors present evidence in support of one such approach: the use of tritium-rich ice to decrease 14 MeV neutron scattering and heating of the dense fuel, resulting in less mass ablation and more robust burn of the hot spot. It is found from numerical simulations that despite a less favorable D:T ratio in the ice, the use of a 40:60 D:T ratio leads to an increase in capsule yield of 17% percent compared with that of a 50:50 D:T ratio fuel for capsules resembling those of the recent N210808 ignition experiment on the NIF (Abu-Shawareb et al., 2022) and an increase of 74% compared with that of a 60:40 D:T ratio fuel capsule. Finally, these results are potentially important for modeling all layered implosions, since some degree of DT fractionization may arise naturally during the beta layering process. In addition, this physics is important for the feasibility of high-gain capsule designs that seek to minimize tritium usage, as in some inertial fusion energy concepts.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗