Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Measurements and kinetic modeling of O 2 vibrational kinetics in O 2 –Ar mixtures partially dissociated by a Ns pulse discharge

Vibrational kinetics of O 2 is studied during the O atom recombination in an O 2 –Ar mixture, partially dissociated by a burst of ns discharge pulses in a heated plasma flow reactor. The time-resolved temperature in the discharge afterglow is determined by Rayleigh scattering. Time-resolved O atom number density is measured by ps Two-Photon absorption Laser Induced Fluorescence, calibrated in xenon. Time-resolved vibrational level populations of molecular oxygen, O 2 (v= 8–20), are measured by ps Laser Induced Fluorescence (LIF), with the absolute calibration by NO LIF. Time-resolved ozone number density is monitored by broadband UV absorption. The results are compared with the predictions of a state-specific kinetic model. The experimental data indicate a rapid initial decay of O 2 (v) populations generated by electron impact in the discharge, due to the vibration-translation (V–T) relaxation by O atoms. This is followed by a slower population reduction, on the time scale much longer compared to that for V–T relaxation or vibration-vibration (V–V) exchange. Both O atoms and the O 2 (v) populations decay on the same time scale, indicating that chemical reactions initiated by the O atom recombination result in the generation of vibrationally excited O 2 molecules. These trends are reproduced by the kinetic model, which shows that the reaction of O atoms with ozone is the dominant pathway of O 2 (v) generation at the present conditions. The predicted relative O 2 (v) populations are close to the experimental results, but absolute number densities differ from the experimental data. This is likely due to uncertainties in the absolute calibration of LIF measurements and in the spectroscopic model used in the data reduction. The present work demonstrates the capability for the absolute, time-resolved measurements of vibrationally excited O 2 in recombining gas flows, to quantify the energy partition in the recombination reactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Erosion and Mobilisation of Highly Active Simulant Suspensions using Impinging Vertical Jets - 19378

Based on the nuclear industry in the U.K., storage of nuclear waste poses a particularly complex challenge. Specifically, this research is concerned with post-PUREX (Plutonium Uranium Redox Extraction) high level waste, known as HAL (Highly Active Liquor), containing nitric acid and a host of fission products. HAL is concentrated, and stored in HASTs (Highly Active Storage Tanks) at the HALES (Highly Active Liquor Evaporation and Storage) facility at the Sellafield site, UK. HASTs are large concrete containers designed to keep the HAL cool and continually circulated prior to vitrification for permanent storage and disposal. Without comprehensive understanding of how this HAL is behaving within the HASTs, a multitude of issues could arise, and understanding such issues is the basis of this research. Components of the HAL begin to form solid particles in the HASTS, and impinging vertical liquid jets are used to mobilize the sedimented particle bed to reduce hot-spot formation. Experimental validation of the jets and associated modelling is required to assess the jets efficacy, and by considering the erosion based on dimensionless parameters, will allow the results to be scaled. This is of particular interest due to the wide reaching applications to all industries that require utilisation of impinging jets, for mass or heat transfer. An experimental rig of 0.5 m{sup 3} has been used to experimentally model the central jet of the HASTs. A spherical glass simulant (SiO{sub 2}, spherical) of various sizes has been used for the initial non-cohesive particle testing, and modelling. Each experiment used 8 different jet diameters. Dimensionless analysis takes the form of the erosion parameter, Ec, which allows comparison of systems of different scales and a means of correlating multiple simulants against one another, so patterns are easily identifiable. Part of the dimensionless analysis was also to deduce the modified Shields parameter from vane viscometry to allow for a comparison of erosion behaviour taking into account particle properties, flow properties and the overall influence these have on the erosion profile. Particle size distribution was found to have a direct effect on yield stress and modified Shields parameter, with larger particles being dominated by inertial forces and past a certain size threshold surface forces appear to dominate. Ec data found the larger glass simulants to behave as expected for a non-cohesive particles, however surface forces also appear to be present for the smaller size distributions. (authors)

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

Mesoscale to Microscale Coupling for a Wind Ramp Case over Complex Terrain

Within the wind energy community, it has become increasingly important to simulate realistic weather events at high fidelity in order to improve wind turbine design, wind farm efficiency, and power forecasting. One of the major difficulties in simulating atmospheric events with numerical models at high fidelity is in accurately transitioning from the mesoscale to the microscale. The mesoscale flow, when input into the large-eddy simulation (LES) domain, will be free of resolved turbulence; needing sufficient time and space (known as fetch) for resolved-scale turbulence to develop. In order to reduce this fetch, several methods exist in which the temperature field, wind field, or both, are stochastically perturbed in order to accelerate the generation of turbulence. However, for cases involving complex terrain, it remains unclear if these perturbation processes are necessary, and if the interaction of the flow over the terrain, alone, is sufficient to generate adequate turbulence. In this study, a wind ramp event observed at the National Wind Technology Center (NWTC) outside of Boulder, CO is simulated first with a mesoscale model and then, using the simulated mesoscale solution as initial and boundary conditions, within a microscale LES solver. The mesoscale solution generated by the Weather Research and Forecasting (WRF) model accurately captures the wind ramp event with only a slight temporal shift and a small temperature bias when compared to the meteorological towers in the NWTC site. This information is then taken as the initial and boundary conditions within the Simulator for Wind Farm Applications (SOWFA) microscale LES model with varying terrain inputs. These inputs include the high-resolution (roughly 30 m grid spacing) terrain data from the Shuttle Radar Topography Mission (SRTM) (Figure 1a), as well as, three intermediate terrain resolution data sets produced by applying a uniform convolution window to the SRTM data with varying window sizes (Figure 1b-d). These are compared to the WRF terrain data interpolated to the high-resolution microscale grid (Figure 1e) in order to determine the impact of low-resolution (about 1 km grid spacing) terrain data on turbulence generation. The wind field is compared at several heights and distances along the dominant wind flow axis between each of the terrain cases. It is found that the development of turbulence across the domain is accelerated as terrain information increases in fidelity. Further, because this ramp event at the location of interest is largely influenced by the interaction of the short-wave trough with the terrain, as the topographic information decreases in fidelity the accuracy of the ramp event is decreased. To make generalizations and describe the benefits and shortcomings of this approach, several of the existing perturbation techniques are employed and compared to that of turbulence generation through topographic features.

17 WIND ENERGY↗
DOE OSTI2020

Extreme Scale Infrasound Inversion and Prediction for Weather Characterization and Acute Event Detection

Accurate and timely weather predictions are critical to many aspects of society with a profound impact on our economy, general well-being, and national security. In particular, our ability to forecast severe weather systems is necessary to avoid injuries and fatalities, but also important to minimize infrastructure damage and maximize mitigation strategies. The weather community has developed a range of sophisticated numerical models that are executed at various spatial and temporal scales in an attempt to issue global, regional, and local forecasts in pseudo real time. The accuracy however depends on the time period of the forecast, the nonlinearities of the dynamics, and the target spatial resolution. Significant uncertainties plague these predictions including errors in initial conditions, material properties, data, and model approximations. To address these shortcomings, a continuous data collection occurs at an effort level that is even larger than the modeling process. It has been demonstrated that the accuracy of the predictions depends on the quality of the data and is independent to a certain extent on the sophistication of the numerical models. Data assimilation has become one of the more critical steps in the overall weather prediction business and consequently substantial improvements in the quality of the data would have transformational benefits. This paper describes the use of infrasound inversion technology, enabled through exascale computing, that could potentially achieve orders of magnitude improvement in data quality and therefore transform weather predictions with significant impact on many aspects of our society.

58 GEOSCIENCES↗
DOE OSTI2023

Disruption event characterization and forecasting in tokamaks

Disruption prediction and avoidance is a critical need for next-step tokamaks, such as ITER. Disruption Event Characterization and Forecasting (DECAF) research fully automates analysis of tokamak data to determine chains of events that lead to disruptions and to forecast their evolution allowing sufficient time for mitigation or complete avoidance of the disruption. Disruption event chains related to local rotating or global magnetohydrodynamic (MHD) modes and vertical instability are examined with warnings issued for many off-normal physics events, including density limits, plasma dynamics, confinement transitions, and profile variations. Along with Greenwald density limit evaluation, a local radiative island power balance theory is evaluated and compared to the observation of island growth. Automated decomposition and analysis of rotating tearing modes produce physical event chains leading to disruptions. A total MHD state warning model comprised of 15 separate criteria produces a disruption forecast about 180 ms before a standard locked mode detector warning. Single DECAF event analyses have begun on KSTAR, MAST, and NSTX/-U databases with thousands of shot seconds of device operation using from 0.5 to 1 × 10 6 tested sample times per device. An initial multi-device database comparison illustrates a highly important result that plasma disruptivity does not need to increase as βN increases. Global MHD instabilities, such as resistive wall modes (RWMs), can give the briefest time period of warning before disruption compared to other physics events. In an NSTX database with unstable RWMs, the mode onset, loss of boundary and current control, and disruption event warnings are found in all cases and vertical displacement events are found in 91% of cases. An initial time-dependent reduced physics model of kinetic RWM stabilization created to forecast the disruption chain predicts instability 84% of the time for experimentally unstable cases with a relatively low false positive rate. Instances of the disruption event chain analysis illustrate dynamics including H–L back transitions for rotating MHD and global RWM triggering events. Furthermore, disruption warnings are issued with sufficient time before the disruption (on transport timescales) to potentially allow active profile control for disruption avoidance, active mode control, or mitigation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluating the X2 initial core zero power physics tests with Serpent-Ants

The validation of the Ants nodal neutronics code for VVER applications is started by modelling the zero power physics tests for the initial core of the Khmelnitsky 2 nuclear power plant as described in the X2 benchmark using the Serpent-Ants two step neutronics calculation chain. The Ants prediction compare favorably against an earlier continuous energy Monte-Carlo reference solution as well as the measured data, except for discrepancies in the SCRAM worth between predicted and measured values. Similar discrepancies have been previously reported for VVER-1000 reactors, and the Ants predicted SCRAM worths match well with the Serpent results. Additionally, the assembly power distribution and axial power distribution predicted by Ants for the critical HZP state of the reactor are also compared to a reference Serpent solution with very good results. Lastly, an investigation into the effects of the few-group structure used in the nodal calculations and the use or lack of leakage correction for group constants on the results shows that the best accuracy is reached with seven or eight energy groups without leakage correction, although a reasonable accuracy can also be obtained with two or three energy groups using fundamental mode leakage correction. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2016

Advanced 3D Characterization and Reconstruction of Reactor Materials FY16 Final Report

A coordinated effort to link advanced materials characterization methods and computational modeling approaches is critical to future success for understanding and predicting the behavior of reactor materials that operate at extreme conditions. The difficulty and expense of working with nuclear materials have inhibited the use of modern characterization techniques on this class of materials. Likewise, mesoscale simulation efforts have been impeded due to insufficient experimental data necessary for initialization and validation of the computer models. The objective of this research is to develop methods to integrate advanced materials characterization techniques developed for reactor materials with state-of-the-art mesoscale modeling and simulation tools. Research to develop broad-ion beam sample preparation, high-resolution electron backscatter diffraction, and digital microstructure reconstruction techniques; and methods for integration of these techniques into mesoscale modeling tools are detailed. Results for both irradiated and un-irradiated reactor materials are presented for FY14 - FY16 and final remarks are provided.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

A Moist Potential Vorticity Model for Midlatitude Long-Lived Mesoscale Convective Systems over Land

Mesoscale convective systems (MCSs) bring large amounts of rainfall and strong wind gusts to the midlatitude land regions, with significant impacts on local weather and hydrologic cycle. However, weather and climate models face a huge challenge in accurately modeling the MCS life cycle and the associated precipitation, highlighting an urgent need for a better understanding of the underlying mechanisms of MCS initiation and propagation. From a theoretical perspective, a suitable model to capture the realistic properties of MCSs and isolate the bare-bones mechanisms for their initiation, intensification, and eastward propagation is still lacking. To simulate midlatitude MCSs over land, we develop a simple moist potential vorticity (PV) model that readily describes the interactions among PV perturbations, air moisture, and soil moisture. Multiple experiments with or without various environmental factors and external forcing are used to investigate their impacts on MCS dynamics and mesoscale circulation vertical structures. The result shows that mechanical forcing can induce lower-level updraft and cooling, providing favorable conditions for MCS initiation. A positive feedback among surface winds, evaporation rate, and air moisture similar to the wind-induced surface heat exchange over tropical ocean is found to support MCS intensification. Both background surface westerlies and vertical westerly wind shear are shown to provide favorable conditions for the eastward propagation of MCSs. Last, our result highlights the crucial role of stratiform heating in shaping mesoscale circulation response. In conclusion, the model should serve as a useful tool for understanding the fundamental mechanisms of MCS dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Data-Centric AI and the Open Energy Data Initiative (OEDI)

This presentation emphasizes the critical importance of data-centric AI. The limitations of model-centric AI when dealing with poor or insufficient data are highlighted, and it is illustrated how training models on inaccurate or noisy data leads to suboptimal results. This talk advocates for a hybrid approach that combines a focus on data quality and model parameters to achieve optimal results. The Open Energy Data Initiative (OEDI) is introduced as a valuable resource for obtaining high-quality energy-related datasets, hosting nearly 2,000 publicly accessible datasets, including 99 solar-related datasets, totaling over 2.7 petabytes of data. OEDI's data lakes enable users to query and work with data without extensive transfers. In conclusion, the significance of data-centric AI and adherence to data curation best practices is emphasized, positioning OEDI as a prime source of high-quality data for AI and machine learning in the renewable energy sector.

AI↗
DOE OSTI2023

Bringing Low- and Moderate-Income Solar Financing Models to Scale (Final Technical Report (FTR))

In this final report, the Clean Energy States Alliance(CESA) and its project partners describe their work to accelerate the development of solar projects for three distinct subsets of the low- to moderate-income(LMI) solar market: single-family homes, manufactured homes, and community institutions, including multifamily affordable housing. For each market sector, the project team undertook research, outreach, and market-building activities. By advancing promising LMI solar financing innovations for different housing types and spreading them to new locations, the Scaling Up Solar for Under-Resourced Communities project has built momentum to scale up LMI solar. Although the report describes some research pertaining to solar adoption in LMI communities, this project was original in that it focused on two specific models to tackle single-family homes and community institutions. Until this project, very little research dedicated to solar for manufactured homes existed. For each market sector, the project team undertook research, outreach, and market-building activities. For the single-family homes sector, the project encouraged states to adapt a program model based on a successful Connecticut initiative that has brought solar to thousands of LMI homeowners. For manufactured homes, the project analyzed the potential for using solar for that housing sector and worked with states, utilities, and other stakeholders to launch pilot projects or programs. For community institutions, the project team worked with foundations, lenders, and community service organizations to inventory, analyze, and communicate models for philanthropic investment that accelerate the deployment of solar and solar+storage in multifamily affordable housing, health centers, and other LMI-serving community institutions.

14 SOLAR ENERGY↗
DOE OSTI2019

Influence of hydrological, biogeochemical and temperature transients on subsurface carbon fluxes in a flood plain environment, Biogeochemistry: Dataset

Data tables used in Arora et al. 2016b Influence of hydrological, biogeochemical and temperature transients on subsurface carbon fluxes in a flood plain environment. Biogeochemistry, 127(2-3), 367-396. Files include reactive transport model parameters describing soil physical and thermal parameters, geochemical reaction network and rate laws, as well as initial and boundary conditions. The last table documents modeling results with respect to annual carbon fluxes from the site to the river.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Effects of electron-phonon coupling and electronic thermal conductivity in high energy molecular dynamics simulations of irradiation cascades in nickel

The two-temperature model has been applied to investigate the effects of the electronic subsystem on 150 keV Ni ion cascades in nickel using molecular dynamics simulation. We explore the effects of the magnitude of the electron-phonon coupling and the electronic thermal conductivity on defect production and cluster formation. It has been found that stronger electron-phonon coupling allows larger and more rapid energy feedback to the atomic subsystem, leading to reduction of number of point defects and suppression of the formation of larger defect clusters. It was observed that larger electronic thermal conductivity results in slightly increased number of point defects and larger size vacancy clusters. The latter takes place because of suppression of point defects recombination in faster cooling areas of initial damage.

Zarkadoula, Eva↗
DOE OSTI2017

Local and global Λ polarization in a vortical fluid

We compute the fermion spin distribution in the vortical fluid created in off-central high energy heavy-ion collisions. We employ the event-by-event (3+1)D viscous hydrodynamic model. The spin polarization density is proportional to the local fluid vorticity in quantum kinetic theory. As a result of strong collectivity, the spatial distribution of the local vorticity on the freeze-out hyper-surface strongly correlates to the rapidity and azimuthal angle distribution of fermion spins. We investigate the sensitivity of the local polarization to the initial fluid velocity in the hydrodynamic model and compute the global polarization of Λ hyperons by the AMPT model. The energy dependence of the global polarization agrees with the STAR data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

A REDUCED ORDER MODELING APPROACH TO PROBABILISTIC CREEP-DAMAGE PREDICTIONS IN FINITE ELEMENT ANALYSIS

This paper introduces a computationally efficient Reduced Order Modeling (ROM) approach for the probabilistic prediction of creep-damage failure. Component-level probabilistic simulations are needed to assess the reliability and safety of high-temperature components. Full-scale probabilistic creep-damage modeling in finite element (FE) approach is computationally expensive requiring many hundreds of simulations to replicate the uncertainty of component failure. To that end, ROM is proposed to minimize the elevated computational cost while controlling the loss of accuracy. It is proposed that full-scale probabilistic simulations can be completed in 1D at a reduced cost, the extremum conditions extracted, and those conditions applied for lower-cost 2D/3D probabilistic simulations of components that capture the mean and uncertainty of failure. The probabilistic Sine-hyperbolic (Sinh) model is selected which in previous work was calibrated to alloy 304 stainless steel. The Sinh model includes probability density functions (pdfs) for test condition (stress and temperature), initial damage (i.e. microstructure), and material properties uncertainty. The Sinh model is programmed into ANSYS finite element software using the USERCREEP.F material subroutine. First, the Sinh model and FE code are subject to verification and validation to ensure the accuracy of the simulations. Numerous Monte Carlo simulations are executed in a 1D model to generate probabilistic creep deformation, damage, and rupture data. This data is analyzed and the probabilistic parameters corresponding to extreme creep response are extracted. The ROM concept is applied where only the extreme conditions are applied in the 2D probabilistic prediction of a component. The probabilistic predictions between the 1D and 2D geometry is compared to assess ROM for creep. The accuracy of the probabilistic prediction employing the ROM approach will potentially reduce the time and cost of simulating complex engineering systems. Future studies will introduce multi-stage Sinh, stochasticity, and spatial uncertainty for improved prediction.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Robustness of the Ensemble Score Filter to the Type of Assimilated Observation Networks

Recent advances in data assimilation (DA) have focused on developing more flexible approaches that can better accommodate nonlinearities in models and observations. However, it remains unclear how the performance of these advanced methods depends on the observation network characteristics. In this study, we present initial experiments with the surface quasi‐geostrophic model, in which we compare a recently developed ensemble filter using score‐based diffusion models with the standard Local Ensemble Transform Kalman Filter (LETKF). Our results show that the analysis solutions respond differently to the number, spatial distribution, and nonlinear fraction of assimilated observations. We also find notable changes in the multiscale characteristics of the analysis errors. Given that standard DA techniques will eventually be replaced by more advanced methods, we hope this study sets the ground for future efforts to reassess the value of Earth observing systems in the context of newly emerging algorithms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Verification of Bison fission product species conservation under TRISO reactor conditions

When assessing the reliability and predictive capabilities of a simulation tool, code verification is used to ensure that the implemented numerical algorithm is a faithful representation of its underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. During this process, numerical results in a discrete solution are compared to the analytical solution of the mathematical model. Here, in this paper, the code verification process is applied to one-dimensional spatiotemporal problems that exercise partial differential equation governing the conservation of fission product species (or mass diffusion). Numerical experiments were performed in the Bison fuel performance code to evaluate its predictive capability under various TRISO reactor conditions such as base irradiation and safety heating test conditions for either short- or long-lived fission product species, as well as a case concerning evaporation from the outer surface of a particle. The code predictions were compared with the expected exact results obtained from the analytical expressions, and the fact that they demonstrate the correct analytical behavior provides strong evidence of proper numerical algorithm implementation.

07 ISOTOPE AND RADIATION SOURCES↗