Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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 361 records · Page 20

A Lattice Boltzmann Method for Predicting Porous Transport Layer Performance During Electrolysis

Electrolysis, the splitting of water into oxygen and hydrogen using electricity, is a sustainable way to produce green hydrogen for energy storage. In polymer electrolyte membrane (PEM) water electrolysis, water is brought into contact with charged catalyst layers and electrochemically separated into oxygen and hydrogen. The hydrogen product formed at the cathode is carried through the catalyst layer for eventual collection, while the oxygen by-product formed at the anode is removed from the surface via a multiphase interaction with circulating water and a solid porous transport layer (PTL). The design of this PTL aids in the detachment and advection of the oxygen by-product and thereby plays a role in the overall efficacy of the catalyst. In this presentation, we present our initial results modeling this multiphase system using a single-component, multiphase lattice Boltzmann method. We use the Shan-Chen model describing inter-particle forces to capture both the cohesion of the water (liquid) and oxygen (gas) phases and their interaction with the PTL (solid) (Shan and Chen, 1993). We use a Carnahan-Starling equation of state to model the effective density governing these inter-particle interactions which allows us to model this relatively high density ratio system (Carnahan and Starling, 1969). With these simulations, we show that the geometry and heterogeneity of the PTL geometry plays a large role in its ability to move oxygen away from the catalyst layer and the resulting bubble structures that are formed within the PTL. The current work demonstrates these effects using synthesized PTL geometries and 2D physics, which will be extended to experimentally-imaged PTL sections and 3D algorithms in the near future.

Boltzmann↗

Proposed Analytical Methods for Determining Filter Media Properties

High Efficiency Particulate Air (HEPA) filters, commonly used in nuclear filtration applications, are an integral part of the waste management processes in nuclear plants. HEPA filters are 99.97% efficient filtration devices characterized by their high resistance to air flow, or pressure drop. From theoretical models, the initial pressure drop across a clean filter is proven to be a function of the filter media fiber diameters and porosity of the media. These filter properties are relatively difficult to obtain with traditional manual methods; therefore, there is a need to develop analytical methods to find the fiber diameter and porosity to determine the pressure drop. Typically, these values can be found by substituting a calculated representative value based upon measured values, such as equivalent fiber diameter based upon a clean pressure drop. However, these values may also be directly recorded and measured by incorporating Scanning Electron Microscope (SEM) image analysis and other measurement methods to analyze the filter media and determine its physical properties without the need for media testing. DiameterJ, an open source Java plug-in, used with ImageJ, can process images of filter media taken by an SEM to find statistical data such as mean fiber diameter and porosity. To produce the raw data, SEM images of two filter media type samples are taken and segmented in DiameterJ using the traditional and statistical region merging segmentation algorithms. Manual segmentation is necessary after the initial segmentation by the algorithms as the images tend to be too complex for the algorithms to output with the necessary accuracy. However, complications exist in the manual segmentation process as these methods can be time intensive and prone to the individual bias of the user. This in turn can skew the final mean fiber diameter result, and lead to either an over or under prediction of the pressure drop. It was also discovered that the porosity data produced by DiameterJ is inaccurate, as the SEM analyzes the three-dimensional filter media by projecting its geometry onto a plane and analyzing it as a two-dimensional binary image. Thus, the porosity is artificially inflated through the segmentation process, rendering this result incorrect. Alternatively, density determination, gravimetric analysis, and thickness testing of the filter media is collectively used to determine the filter fiber porosity. Together, the SEM image analysis and analytical lab methods produce results through direct measurements which allow for the prediction of the initial pressure drop from clean filter media without the need for prior media testing to collect pressure data. By improving upon this proposed analytical method in the future, there is potential to streamline the process of finding these filter properties into a more direct methodology for determining the pressure drop across HEPA filters.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Loss Informed Dispatcher for Battery Systems (LIDBS)

LIDBS takes battery-system operational data as an input. It then cleans the data using density-based clustering methods, forms an initial loss model, and then embeds the loss model into the optimization formulation. LIDBS provides a much more informed dispatch strategy over state-of-the-art approaches.

Latif, Aadil↗

SubTSBR to tackle high noise and outliers for data-driven discovery of differential equations

Data-driven discovery of differential equations has been an emerging research topic. We propose a novel algorithm subsampling-based threshold sparse Bayesian regression (SubTSBR) to tackle high noise and outliers. The subsampling technique is used for improving the accuracy of the Bayesian learning algorithm. It has two parameters: subsampling size and the number of subsamples. When the subsampling size increases with fixed total sample size, the accuracy of our algorithm goes up and then down. When the number of subsamples increases, the accuracy of our algorithm keeps going up. We demonstrate how to use our algorithm step by step and compare our algorithm with threshold sparse Bayesian regression (TSBR) for the discovery of differential equations. We show that our algorithm produces better results. We also discuss the merits of discovering differential equations from data and demonstrate how to discover models with random initial and boundary condition as well as models with bifurcations. The numerical examples are: (1) predator-prey model with noise, (2) shallow water equations with outliers, (3) heat diffusion with random initial and boundary condition, and (4) fish-harvesting problem with bifurcations.

97 MATHEMATICS AND COMPUTING↗

Initial demonstration of automated fuel performance modeling with 1977 EBR-II metallic fuel pins using BISON code with FIPD and IMIS databases

Using the BISON fuel performance code, simulations were conducted using an automated process to read initial and operating conditions from the Fuels Irradiation and Physics Database (FIPD) and Integral Fast Reactor materials information system (IMIS) database, which contains metallic fuel data from the Experimental Breeder Reactor-II (EBR-II). This work demonstrates use of an integrated framework to access the vast majority of EBR-II experimental fuel pin data to support rapid development of fuel performance models for next-generation metallic fuel systems. With this capability, validation for fuel qualification can be performed rapidly. Between IMIS and FIPD, there is enough information to conduct 1977 unique EBR-II metallic fuel pin histories from 24 different experiments, at varying levels of detail between the two databases. Each of these histories includes a high-resolution power history, flux history, coolant channel flow rates, and coolant channel temperatures. Fission gas release (FGR), cumulative damage fraction (CDF), fuel axial swelling, cladding profilometry, and burnup were all simulated in BISON. The results were compared to post-irradiation examination (PIE) results for the initial demonstration of automated BISON modeling. BISON simulations conducted with IMIS and FIPD were in rough agreement with PIE measurements and calculations. Cladding profilometry, FGR, and fuel axial swelling were found to be in rough agreement with PIE measurements, depending on the physics used within the BISON input files. Here, the mechanical contact solver chosen was found to significantly impact axial fuel swelling and cladding strain predictions. CDF values were assessed to see whether pin failure may have been predicted (CDF ≥ 1). This work suggests that continued development of an automated tool for BISON should focus on inclusion of the Fast Flux Test Facility (FFTF) experimental data for a larger database for metallic fuel, improved physical models to better capture fuel performance, such as fuel-cladding interactions, and a more detailed comparison with available PIE data to further the BISON model development.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Implementing a microphysics model in hydrodynamic simulations to study the initial plasma formation in dielectric ablator materials for direct-drive implosions

A microphysics model to describe the photoionization and impact ionization processes in dielectric ablator materials like plastic has been implemented into the 1-D hydrodynamic code LILAC for planar and spherical targets. Currently, the initial plasma formation during the early stages of a laser drive are modeled in an ad hoc manner, until the formation of a critical surface. Implementation of the physics-based models predict higher values of electron density, electron temperature, and pressure than the ad hoc model. Moreover, the numerical predictions are consistent with previous experimental observations of the shinethrough mechanism in plastic ablators. For planar targets, a decompression of the rear end of the target was observed, that is similar to recent experiments. An application of this model is to understand the laser-imprint mechanism that is caused by nonuniform laser irradiation due to single beam speckle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Data for A Hybrid Biophysical-Machine Learning Framework for Diurnal Surface Energy Flux Estimation Using Proximal Sensing

Thermal infrared-based remote sensing of surface energy fluxes has traditionally relied on high spatial resolution satellite data with revisit frequencies on the order of weeks. In this study, we evaluate a biophysics-based analytical surface energy balance model for predicting latent energy (LE) and sensible heat (H) fluxes using proximal sensing observations. The Surface Temperature Initiated Closure (STIC1.2) model has been extensively validated across a wide range of spatial and temporal scales using various satellite-derived thermal infrared data sets. Here we extend this validation by applying STIC at sub-hourly temporal resolution over multiple growing seasons for four distinct agricultural systems. We further develop and evaluate novel STIC variants that incorporate machine learning (ML) techniques to eliminate the need for surface energy balance observations, specifically net radiation and soil heat flux, thereby enhancing model applicability in data-sparse settings. The integration of a ML component to estimate surface available energy is shown to have strong predictive performance for both LE (R2 = 0.81–0.94) and H (R2 = 0.46–0.72) across all agricultural systems examined here, demonstrating the potential of hybrid biophysical-machine learning approaches for surface energy balance modeling with minimal data requirements. This study concludes with a novel application of explainable machine learning (exML) to diagnose sources of model error. This exML framework attributes residual prediction errors to both model input variables and environmental drivers not explicitly included in the simulation experiments. This approach provides a new pathway for improving model design and integrating previously overlooked yet influential variables into future model iterations.

AI/ML↗

Forecasting the Decay of Hybrid Perovskite Performance Using Optical Transmittance or Reflected Dark-Field Imaging

The practicality and economic viability of hybrid perovskite solar cells hinge on their operational lifetime, and methods for forecasting the performance of perovskites under different operational stresses are urgently needed. Here, we explore the evolution of material-level optoelectronic properties as MAPbI 3 degrades and discover universal behaviors where the carrier diffusion length (L D ) decays before quasi-Fermi-level splitting (ΔE F ), regardless of the specific stress protocol (oxygen, humidity, thermal stress, or a combination). We employ a machine learning greedy feature selection model that uses initially measured properties to predict the time it takes L D to decrease to 85% of its initial value with a prediction accuracy of 12.8%. This model reveals a strong correlation between the initial rate of transmittance change and the time until loss of transport. We translate this material-level finding to photovoltaic device-level forecasting by demonstrating that the rate of change of transmittance is equivalent to the rate of change of the spatial standard deviation of dark-field image intensity (i.e., scattered light intensity) collected in reflection mode (and thus applicable to devices with opaque contacts). Furthermore, this work demonstrates that transmittance and scattering methods are highly effective for accelerated material (and device) stability evaluation and forecasting.

14 SOLAR ENERGY↗

Experimental protocol for phase 1 of the APARC QUOCA (QUasibiennial oscillation and Ozone Chemistry interactions in the Atmosphere) working group

The quasi-biennial oscillation (QBO) is the main mode of variability in the tropical stratosphere, influencing the predictability of other regions in the atmosphere through its teleconnections to the stratospheric polar vortices and coupling to surface tropical and extratropical variability. However, climate and forecasting models consistently underestimate QBO amplitudes in the lower stratosphere, likely contributing to their failure to simulate these teleconnections. One underexplored contributor to model biases is missing representation of ozone-radiative feedbacks, which enhance temperature variability in the lower stratosphere, particularly at periods at and greater than the QBO (>28 months). While previous studies suggest that ozone-radiative feedbacks can impact QBO periods, amplitudes and the associated secondary circulation in the lower stratosphere, the reported impacts differ widely among models and are hard to interpret due to differences in methodology. To this end, here we propose a coordinated experimental protocol – held joint between the Atmospheric Processes and their Role in Climate (APARC) Quasi-Biennial Oscillation Initiative (QBOi) and Chemistry Climate Modeling Initiative (CCMI) activities – which is aimed at assessing the coupling between stratospheric ozone, temperature and the circulation. We use the proposed experiments to define the ozone feedback on the QBO in both present-day and idealized (abrupt quadrupling of carbon dioxide) climates. While primary focus is on the QBO, the proposed protocol also enables analysis of other aspects of ozone-radiative-dynamical coupling in the atmosphere, including impacts on the Brewer-Dobson Circulation and tropospheric eddy-driven jet responses to future climate change. Here we document the scientific rationale and design of the QUOCA Phase 1 experiments, summarize the data request, and give a brief overview of participating models. Preliminary results using the NASA Goddard Institute for Space Studies E2-2 climate model are used to illustrate sensitivities to certain methodological choices.

Environmental sciences↗

A study of shock initiation experiments for the explosive PBX 9502 using three reactive burn models

Shock to detonation transition (SDT) experiments are essential in calibrating and validating reactive burn models for explosives. This work investigates the large collection of SDT test data for the explosive PBX 9502 at ambient temperature that was presented by Gustavsen, Sheffield, and Alcon [Journal of Applied Physics, 99, 114907 (2006)]. We first analyze the experimental data and compare two different methods of determining the shock transition time/distance (namely, the bilinear method and the single-curve method). This reveals some of the uncertainty in estimating shock transition points, which contributes to scatter in Pop-plot data. Next, we compare the WSD, AWSD, and SURFplus reactive burn models for a collection of approximately 20 experimental shots using the FLAG hydrocode. Error estimates are used to quantify how well each reactive burn model (and their respective parameter calibrations) performs at predicting the SDT process for a range of loading conditions. Additionally, the importance of mesh resolution and numerical dissipation in SDT simulations will be assessed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

LOGOS

LOGOS is a software package which contains a set of discrete optimization models that can be employed for capital budgeting optimization problems. More specifically, provided a set of items (characterized by cost and reward values) and constraints, these models select the best combination of items which maximizes overall reward and satisfies the provided constraints. The developed models are based on different versions of the knapsack optimization algorithms. Two main classes of optimization models have been initially developed: deterministic and stochastic. Stochastic optimization models evolve deterministic models by explicitly considering data uncertainties (associated to constraints or item cost and reward). These models can be employed as stand-alone models or interfaced with the INL developed RAVEN code to propagate data uncertainties and analyze the generated data (i.e., sensitivity analysis).

Wang, Congjian↗

A coupled ground heat flux–surface energy balance model of evaporation using thermal remote sensing observations

Abstract. One of the major undetermined problems in evaporation (ET) retrieval using thermal infrared remote sensing is the lack of a physically based ground heat flux (G) model and its integration within the surface energy balance (SEB) equation. Here, we present a novel approach based on coupling a thermal inertia (TI)-based mechanistic G model with an analytical surface energy balance model, Surface Temperature Initiated Closure (STIC, version STIC1.2). The coupled model is named STIC-TI. The model is driven by noon–night (13:30 and 01:30 local time) land surface temperature, surface albedo, and a vegetation index from MODIS Aqua in conjunction with a clear-sky net radiation sub-model and ancillary meteorological information. SEB flux estimates from STIC-TI were evaluated with respect to the in situ fluxes from eddy covariance measurements in diverse ecosystems of contrasting aridity in both the Northern Hemisphere and Southern Hemisphere. Sensitivity analysis revealed substantial sensitivity of STIC-TI-derived fluxes due to the land surface temperature uncertainty. An evaluation of noontime G (Gi) estimates showed 12 %–21 % error across six flux tower sites, and a comparison between STIC-TI versus empirical G models also revealed the substantially better performance of the former. While the instantaneous noontime net radiation (RNi) and latent heat flux (LEi) were overestimated (15 % and 25 %), sensible heat flux (Hi) was underestimated (22 %). Overestimation (underestimation) of LEi (Hi) was associated with the overestimation of net available energy (RNi−Gi) and use of unclosed surface energy balance flux measurements in LEi (Hi) validation. The mean percent deviations in Gi and Hi estimates were found to be strongly correlated with satellite day–night view angle difference in parabolic and linear pattern, and a relatively weak correlation was found between day–night view angle difference versus LEi deviation. Findings from this parameter-sparse coupled G–ET model can make a valuable contribution to mapping and monitoring the spatiotemporal variability of ecosystem water stress and evaporation using noon–night thermal infrared observations from future Earth observation satellite missions such as TRISHNA, LSTM, and SBG.

Bhattacharya, Bimal K.↗

Hierarchical Multi-scale Modeling of Laser Weld Failure: An Initial Investigation

The failure of 304L laser welds is of interest to system and component designers due to nuclear safety requirements for abnormal environments. Accurately modeling laser weld behavior in full system and component models has proven especially challenging due to three factors: the large variability observed in laser weld characterization tests; the difficulty in isolating the weld material for material characterization and modeling the weld material behavior; and the disparate scales associated with modeling laser welds in large systems. Recent work has shown that meso-scale geometric features of laser welds such as pores and weld root tortuosity are critical to accurately predicting the structural performance of welds. The challenge with modeling these welds is that the geometric features driving their structural performance are generally on the order of ten to hundreds of microns, but can affect the responses of interest in systems and components on the order of centimeters to meters.

36 MATERIALS SCIENCE↗

Complete development of LAROMance-type surrogate models for Zircaloy, as well as initial integration and testing in MOOSE, which addresses industry needs for licensing advanced LWR fuels

This milestone report details the progress in the constitutive surrogate modeling framework for Zircaloy’s mechanical behavior under extreme conditions, within the NEAMS program throughout the U.S. fiscal year 2024 (FY24). The presented work expands on the spectral based mechanical solver and the mechanistic material model for Zircaloy, which were introduced in the NEAMS program during FY22 and detailed in the previous fiscal year (FY23). As part of the NEAMS program, the current effort described in this report focused on the development of a novel viscplastic (creep) surrogate for the mechanistic, polycrystal-based model. To this end, the work was subdivided into two interdependent components related to the necessary improvements to the LApx constitutive model and algorithm and, to the development and use of a surrogate model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗