Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Surrogates”

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,387 records · Page 77

Polynomial Chaos–Based Bayesian Inference of K -Profile Parameterization in a General Circulation Model of the Tropical Pacific

The authors present a polynomial chaos (PC)–based Bayesian inference method for quantifying the uncertainties of the K-profile parameterization (KPP) within the MIT general circulation model (MITgcm) of the tropical Pacific. The inference of the uncertain parameters is based on a Markov chain Monte Carlo (MCMC) scheme that utilizes a newly formulated test statistic taking into account the different components representing the structures of turbulent mixing on both daily and seasonal time scales in addition to the data quality, and filters for the effects of parameter perturbations over those as a result of changes in the wind. To avoid the prohibitive computational cost of integrating the MITgcm model at each MCMC iteration, a surrogate model for the test statistic using the PC method is built. Because of the noise in the model predictions, a basis-pursuit-denoising (BPDN) compressed sensing approach is employed to determine the PC coefficients of a representative surrogate model. The PC surrogate is then used to evaluate the test statistic in the MCMC step for sampling the posterior of the uncertain parameters. Results of the posteriors indicate good agreement with the default values for two parameters of the KPP model, namely the critical bulk and gradient Richardson numbers; while the posteriors of the remaining parameters were barely informative.

54 ENVIRONMENTAL SCIENCES↗

Uncertainty Propagation in Coupled Atmosphere–Wave–Ocean Prediction System: A Study of Hurricane Earl (2010)

This study focuses on understanding the evolution of Hurricane Earl (2010) with respect to random perturbations in the storm’s initial strength, size, and asymmetry in wind distribution. We rely on the Unified Wave Interface-Coupled Model (UWIN-CM), a fully coupled atmosphere–wave–ocean system to generate a storm realization ensemble, and use polynomial chaos (PC) expansions to build surrogate models for time evolution of both the maximum wind speed and minimum sea level pressure in Earl. The resulting PC surrogate models provide statistical insights on probability distributions of model responses throughout the simulation time span. Statistical analysis of rapid intensification (RI) suggests that initial perturbations having intensified and counterclockwise-rotated winds are more likely to undergo RI. In addition, for the range of initial conditions considered RI seems mostly sensitive to azimuthally averaged maximum wind speed and asymmetry orientation, rather than storm size and asymmetry magnitude; this is consistent with global sensitivity analysis of PC surrogate models. Finally, we combine initial condition perturbations with a stochastic kinetic energy backscatter scheme (SKEBS) forcing in the UWIN-CM simulations and conclude that the storm tracks are substantially influenced by the SKEBS forcing perturbations, whereas the perturbations in initial conditions alone had only limited impact on the storm-track forecast.

54 ENVIRONMENTAL SCIENCES↗

Application of an automated machine learning-genetic algorithm (AutoML-GA) coupled with computational fluid dynamics simulations for rapid engine design optimization

In recent years, the use of machine learning-based surrogate models for computational fluid dynamics (CFD) simulations has emerged as a promising technique for reducing the computational cost associated with engine design optimization. However, such methods still suffer from drawbacks. One main disadvantage is that the default machine learning (ML) hyperparameters are often severely suboptimal for a given problem. This has often been addressed by manually trying out different hyperparameter settings, but this solution is ineffective in case of a high-dimensional hyperparameter space. Besides this problem, the amount of data needed for training is also not known a priori. In response to these issues that need to be addressed, the present work describes and validates an automated active learning approach, AutoML-GA, for surrogate-based optimization of internal combustion engines. In this approach, a Bayesian optimization technique is used to find the best machine learning hyperparameters based on an initial dataset obtained from a small number of CFD simulations. Subsequently, a genetic algorithm is employed to locate the design optimum on the ML surrogate surface. In the vicinity of the design optimum, the solution is refined by repeatedly running CFD simulations at the projected optima and adding the newly obtained data to the training dataset. It is demonstrated that AutoML-GA leads to a better optimum with a lower number of CFD simulations, compared to the use of default hyperparameters. The proposed framework offers the advantage of being a more hands-off approach that can be readily utilized by researchers and engineers in industry who do not have extensive machine learning expertise.

Owoyele, Opeoluwa↗

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗

Polynomial Chaos Characterization of Uncertainty in Multiscale Models and Behavior of Carbon Reinforced Composites

Design of non-crimp fabric (NCF) composites entails major challenges pertaining to (1) the complex fine-scale morphology of the constituents, (2) the manufacturing-produced inconsistency of this morphology spatially, and thus (3) the ability to build reliable, robust, and efficient computational surrogate models to account for this complex nature. Traditional approaches to construct computational surrogate models have been to average over the fluctuations of the material properties at different scale lengths. This fails to account for the fine-scale features and fluctuations in morphology, material properties of the constituents, as well as fine-scale phenomena such as damage and cracks. In addition, it fails to accurately predict the scatter in macroscopic properties, which is vital to the design process and behavior prediction. In this work, funded in part by the Department of Energy, we present an approach for addressing these challenges by relying on polynomial chaos representations of both input parameters and material properties at different scales. Moreover, we emphasize the efficiency and robustness of integrating the polynomial chaos expansion with multiscale tools to perform multiscale assimilation, characterization, propagation, and prediction, all of which are necessary to construct the data-driven surrogate models required to design under the uncertainty of composites. These data-driven constructions provide an accurate map from parameters (and their uncertainties) at all scales and the system-level behavior relevant for design. While this perspective is quite general and applicable to all multiscale systems, NCF composites present a particular hierarchy of scales that permits the efficient implementation of these concepts.

42 ENGINEERING↗

Characterization of swift heavy ion irradiation damage in ceria

We have examined microstructural evolution in irradiated ceria (CeO2) using swift heavy ion irradiation, electron microscopy, and atomistic simulation. CeO2, a UO2 fuel surrogate, was irradiated with gold ions at an energy of 1 GeV to fluences up to 1x1014 ions/cm2. Transmission electron microscopy accompanied by electron energy loss spectroscopy showed that the ion tracks were of similar size at all fluences, and that there was no chemical change in the ion track core. Classical molecular dynamics simulations of thermal spikes in CeO2 with energy deposition of 12 and 36 keV/nm show damage consisting of isolated point defects at the lower energy and defect clusters at 36 keV/nm, with no amorphization at either energy. Inferences are drawn from modeling about density changes in the ion track and the formation of interstitial loops that shed light on features observed by electron microscopy of swift heavy ion irradiated ceria.

Yablinsky, Clarissa↗

IDAES-PSE 2.6.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications. IDAES-PSE 2.6.0 Release Highlights Upcoming Changes IDAES will be switching to the new Pyomo solver interface in the next release. Whilst this will hopefully be a smooth transition for most users, there are a few important changes to be aware of. The new solver interface uses a different version of the IPOPT writer (“ipopt_v2”) and thus any custom configuration options you might have set for IPOPT will not carry over and will need to be reset. By default, the new Pyomo linear presolver will be activated with ipopt_v2. Whilst are working to identify any bugs in the presolver, it is possible that some edge cases will remain. IDAES will begin deploying a new set of scaling tools and APIs over the next few releases that make use of the new solver writers. The old scaling tools and APIs will remain for backward compatibility but will begin to be deprecated. New Models, Tools and Features New Intersphinx extension automatically linking Jupyter notebook examples to project documentation New end-to-end diagnostics example demonstrated on a real problem New complementarity formulation for VLE with cubic equations of state, backward compatibility for old formulation New solver interface with presolve (ipopt_v2) in support of upcoming changes to the initialization and APIs methods, with default set to ipopt to maintain backwards compatibility; this will deprecate once all examples have been updated New forecaster and parameterized bidder methods within grid integration library Updated surrogates API and examples to support Keras 3, with backwards compatibility for older formats such as TensorFlow SavedModel (TFSM) Updated costing base dictionary to include the 2023 cost year index value Updated ProcessBlock to include information on the constructing block class Updated Flowsheet Visualizer to allow visualize() method to return value and functions Bug Fixes Fixed bug in the Modular Property Framework that would cause errors when trying to use phase-based material balances with phase equilibria. Fixed bug in Modular Properties Framework that caused errors when initializing models with non-vapor-liquid phase equilibria. Fixed typos flagged by June update to crate-ci/typos and removed DMF-related exceptions Minor corrections of units of measurement handling in power plant waste/transport costing expressions, control volume material holdup expressions, and BTX property package parameters Fixed throwing >7500 numpy deprecation warnings by replacing scalar value assignment with element extraction and item iteration calls Testing and Robustness Migrated slow tests (>10s) to integration, impacting test coverage but also yielding a nearly 30% decrease in local test runtime Pinned pint to avoid issues with older supported Python versions Pinned codecov versions to avoid tokenless upload behavior with latest version Bumped extensions to version 3.4.2 to allow pointing to non-standard install location Deprecations and Removals Python 3.8 is no longer supported. The supported Python versions are 3.9 through 3.12 The Data Management Framework (DMF) is no longer supported. Importing idaes.core.dmf will cause a deprecation warning to be displayed until the next release The SOFC Keras surrogates have been removed. The current version of the SOFC surrogate model in the examples repository is a PySMO Kriging model.

AS↗

CCSI Toolset 3.19 Release

CCSI Toolset 3.19 Release Highlights A gradient generation tool was developed to support GENN models in FOQUS. Certain machine learning tools train gradient-enhanced neural network (GENN) models which can be more accurate for complex datasets given a priori knowledge of model derivatives. However, the derivatives must be known beforehand and are not often available for process data. This tool automatically predicts the gradients for a training dataset in a form usable by common GENN trainers, such as Surrogate Modeling Toolbox. Support was added for Surrogate Modeling Toolbox GENN models in FOQUS, including updates to the run methods, node properties, test framework, documentation and optional dependencies list. Users can train/save Surrogate Modeling Toolbox gradient-enhanced neural network (GENN) models with custom objects and produce .pkl files compatible with the Machine Learning/Artificial Intelligence Plugin in FOQUS. A simpler implementation of the ordering algorithm in the Sequential Design of Experiments (SDOE) module was included. The SDOE examples documentation was updated. The Optimality-Based Design of Experiments was updated to improve the error handling when the results are None.

AS↗

Deep Learning–Assisted Multiobjective Optimization of Geological CO 2 Storage Performance under Geomechanical Risks

In geological CO 2 storage, designing the optimal well control strategy for CO 2 injection to maximize CO 2 storage while minimizing the associated geomechanical risks is not trivial. This challenge arises due to pressure buildup, CO 2 plume migration, the highly nonlinear nature of geomechanical responses to rock-fluid interaction, and the high computational cost associated with coupled flow and geomechanics simulations. In this paper, we introduce a novel optimization framework to address these challenges. The optimization problem is formulated as follows: maximize total CO 2 storage while minimizing geomechanical risks by adjusting the injection schedules within bounded constraints. The geomechanical risks are primarily driven by injection-induced pressure build-up, which is characterized by ground displacement and the induced microseismicity. We used the Fourier neural operator (FNO)-based deep learning model to construct surrogate models, replacing the time-consuming coupled flow and geomechanics simulations for evaluating the aforementioned objective functions. The developed surrogate models have been incorporated into a multiobjective optimization framework through a genetic algorithm to reduce the computational burden. The proposed optimization framework reduces the computational cost from approximately 2,400 hours, when using objective function evaluations based on physics-based simulations, to around 20 minutes. A set of Pareto-optimal solutions of the proposed workflow yields nontrivial optimal decisions, reducing the microseismicity potential and the vertical displacement. This Pareto front highlights the optimal trade-offs between CO 2 storage amount, safety, and ground displacement, emphasizing the need for careful optimization and management of injection strategies to achieve a balanced outcome. The novelty of this work is twofold. First, we demonstrate the importance of incorporating the minimization of the geomechanical risks as objective functions into the CO 2 storage optimization workflow to mitigate the potential risk of induced microseismicity and ground displacement. Second, we leverage the FNO-based surrogate models to optimize a real-field CO 2 storage operation.

42 ENGINEERING↗

Sulfur Polymer Stabilization/Solidification Treatability Study of Mercury Contaminated Soil from the Y-12 Site

As a result of past operations, the Department of Energy’s (DOE) Oak Ridge Y-12 National Security Complex (Y-12 Plant) has extensive mercury-contamination in building structures, soils, storm sewer sediments, and stream sediments, which are a source of pollution to the local ecosystem. Because of mercury’s toxicity and potential impacts on human health and the environment, DOE continues to investigate and implement projects to support the remediation of the Y-12 site. URS and #9122;CH2M Oak Ridge LLC (UCOR) under its prime contract with DOE has cleanup responsibilities on the DOE Oak Ridge Reservation and is investigating potential mercury-contaminated soil treatment technologies through an agreement with Babcock and Wilcox (B and W) Y-12, the Y-12 operating contractor to DOE. As part of its investigations, UCOR has subcontracted with Brookhaven National Laboratory (BNL) to conduct laboratory-scale studies evaluating the applicability of the Sulfur Polymer Stabilization/Solidification (SPSS) process using surrogate and actual mixed waste Y-12 soils containing mercury (Hg) at 135, 2,000, and 10,000 ppm. SPSS uses a thermoplastic sulfur binder to convert Hg to stable mercury sulfide (HgS) and solidifies the chemically stable product in a monolithic solid final waste form to reduce dispersion and permeability. Formulations containing 40 – 60 dry wt% Y-12 soil were fabricated and samples were prepared in triplicate for Environmental Protection Agency Toxicity Characteristic Leaching Procedure (TCLP) testing by an independent laboratory. Those containing 50 and 60 wt% soil easily met the study criteria for maximum allowable Hg concentrations (47 and 1 ppb, respectively compared with the TCLP limit of 200 ppb Hg). The lowest waste loading of 40 wt% yielded TCLP Hg concentrations slightly higher (240 ppb) than the allowable limit. Since the Y-12 soil tended to form clumps, the improved leaching at higher waste loadings was probably due to reduction in particle size from friction of the soil mixing, which creates more surface area for chemical conversion. This was corroborated by the fact that the same waste loading pre-treated by ball milling to reduce particle size prior to SPSS processing yielded TCLP concentrations almost 30 times lower, and at 8.5 ppb Hg was well below EPA limits. Pre-treatment by ball milling also allowed a reduction in the time required for stabilization, thus potentially reducing total process times by 30%.Additional performance testing was conducted including measurement of compressive strength to confirm mechanical integrity and immersion testing to determine the potential impacts of storage or disposal under saturated conditions. For both surrogate and actual Y-12 treated soils, waste form compressive strengths ranged between 2,300 and 6,500 psi, indicating very strong mechanical integrity (a minimum of greater than 40 times greater than the NRC guidance for low-level radioactive waste). In general, compressive strength increases with waste loading as the soil acts as an aggregate in the sulfur concrete waste forms. No statistically significant loss in strength was recorded for the 30 and 40 wt% surrogate waste samples and only a minor reduction in strength was measured for the 43 wt% waste forms. The 30 wt% Y-12 soil did not show a significant loss in strength but the 50 wt% samples were severely degraded in immersion due to swelling of the clay soil. The impact on Hg leaching, if any, was not determined.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Decomposition Analysis of Neodymium Loaded Resin

In an effort to understand and quantify the carbon and sulfur contents of a curium oxide production process performed in the hot cells at Oak Ridge National Laboratory’s Radiochemical Engineering Development Center, a non-radioactive surrogate was used in a similar oxide formation process. Neodymium (Nd), which was chosen as a substitute for curium, was loaded onto a Dowex 50W-X8 resin by use of ion exchange column techniques and then treated under a variety of processing conditions. The surrogate product was analyzed for carbon and sulfur impurities, using inductively coupled plasma–mass spectrometry for sulfur assessment, a LECO combustion furnace/analyzer for carbon quantification, and X-ray diffraction for compound identification. The results indicate that the carbon and sulfur contents in the Nd oxide were similar to current estimations for curium oxide using the reference hot cell process. It was also seen that the major surrogate product of the hot cell process contained sulfur in a Nd oxysulfate (Nd2O2SO4) compound. Parametric analysis of the resin decomposition process revealed that increasing the temperature to 1400°C yielded considerably better removal of carbon and sulfur.

Curium Targets↗

Initiating the Validation of CCIM Processability for Multi-phase all Ceramic (SYNROC) HLW Form: Plan for Test BFY14CCIM-C

This plan covers test BFY14CCIM-C which will be a first–of–its-kind demonstration for the complete non-radioactive surrogate production of multi-phase ceramic (SYNROC) High Level Waste Forms (HLW) using Cold Crucible Induction Melting (CCIM) Technology. The test will occur in the Idaho National Laboratory’s (INL) CCIM Pilot Plant and is tentatively scheduled for the week of September 15, 2014. The purpose of the test is to begin collecting qualitative data for validating the ceramic HLW form processability advantages using CCIM technology- as opposed to existing ceramic–lined Joule Heated Melters (JHM) currently producing BSG HLW forms. The major objectives of BFY14CCIM-C are to complete crystalline melt initiation with a new joule-heated resistive starter ring, sustain inductive melting at temperatures between 1600 to 1700°C for two different relatively high conductive materials representative of the SYNROC ceramic formation inclusive of a HLW surrogate, complete melter tapping and pouring of molten ceramic material in to a preheated 4 inch graphite canister and a similar canister at room temperature. Other goals include assessing the performance of a new crucible specially designed to accommodate the tapping and pouring of pure crystalline forms in contrast to less recalcitrant amorphous glass, assessing the overall operational effectiveness of melt initiation using a resistive starter ring with a dedicated power source, and observing the tapped molten flow and subsequent relatively quick crystallization behavior in pans with areas identical to standard HLW disposal canisters. Surrogate waste compositions with ceramic SYNROC forming additives and their measured properties for inductive melting, testing parameters, pre-test conditions and modifications, data collection requirements, and sampling/post-demonstration analysis requirements for the produced forms are provided and defined.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Determining the Release of Radionuclides from Tank Waste Residual Solids: FY2015 Report

Methodology development for pore water leaching studies has been continued to support Savannah River Site High Level Waste tank closure efforts. For FY2015, the primary goal of this testing was the achievement of target pH and Eh values for pore water solutions representative of local groundwater in the presence of grout or grout-representative (CaCO 3 or FeS) solids as well as waste surrogate solids representative of residual solids expected to be present in a closed tank. For oxidizing conditions representative of a closed tank after aging, a focus was placed on using solid phases believed to be controlling pH and E h at equilibrium conditions. For three pore water conditions (shown below), the target pH values were achieved to within 0.5 pH units. Tank 18 residual surrogate solids leaching studies were conducted over an E h range of approximately 630 mV. Significantly higher Eh values were achieved for the oxidizing conditions (ORII and ORIII) than were previously observed. For the ORII condition, the target Eh value was nearly achieved (within 50 mV). However, E h values observed for the ORIII condition were approximately 160 mV less positive than the target. E h values observed for the RRII condition were approximately 370 mV less negative than the target. Achievement of more positive and more negative E h values is believed to require the addition of non-representative oxidants and reductants, respectively. Plutonium and uranium concentrations measured during Tank 18 residual surrogate solids leaching studies under these conditions (shown below) followed the general trends predicted for plutonium and uranium oxide phases, assuming equilibrium with dissolved oxygen. The highest plutonium and uranium concentrations were observed for the ORIII condition and the lowest concentrations were observed for the RRII condition. Based on these results, it is recommended that these test methodologies be used to conduct leaching studies with actual Tank 18 residual solids material. Actual waste testing will include leaching evaluations of technetium and neptunium, as well as plutonium and uranium.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Reduced Order Model Implementation in the Risk-Informed Safety Margin Characterization Toolkit

The RISMC project aims to develop new advanced simulation-based tools to perform Probabilistic Risk Analysis (PRA) for the existing fleet of U.S. nuclear power plants (NPPs). These tools numerically model not only the thermo-hydraulic behavior of the reactor primary and secondary systems but also external events temporal evolution and components/system ageing. Thus, this is not only a multi-physics problem but also a multi-scale problem (both spatial, µm-mm-m, and temporal, ms-s-minutes-years). As part of the RISMC PRA approach, a large amount of computationally expensive simulation runs are required. An important aspect is that even though computational power is regularly growing, the overall computational cost of a RISMC analysis may be not viable for certain cases. A solution that is being evaluated is the use of reduce order modeling techniques. During the FY2015, we investigated and applied reduced order modeling techniques to decrease the RICM analysis computational cost by decreasing the number of simulations runs to perform and employ surrogate models instead of the actual simulation codes. This report focuses on the use of reduced order modeling techniques that can be applied to any RISMC analysis to generate, analyze and visualize data. In particular, we focus on surrogate models that approximate the simulation results but in a much faster time (µs instead of hours/days). We apply reduced order and surrogate modeling techniques to several RISMC types of analyses using RAVEN and RELAP-7 and show the advantages that can be gained.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Collaborative Project: Building improved optimized parameter estimation algorithms to improve methane and nitrogen fluxes in a climate model

Soils in natural and managed ecosystems and wetlands are well known sources of methane, nitrous oxides, and reactive nitrogen gases, but the magnitudes of gas flux to the atmosphere are still poorly constrained. Thus, the reasons for the large increases in atmospheric concentrations of methane and nitrous oxide since the preindustrial time period are not well understood. The low atmospheric concentrations of methane and nitrous oxide, despite being more potent greenhouse gases than carbon dioxide, complicate empirical studies to provide explanations. In addition to climate concerns, the emissions of reactive nitrogen gases from soils are important to the changing nitrogen balance in the earth system, subject to human management, and may change substantially in the future. Thus improved modeling of the emission fluxes of these species from the land surface is important. Currently, there are emission modules for methane and some nitrogen species in the Community Earth System Model’s Community Land Model (CLM-ME/N); however, there are large uncertainties and problems in the simulations, resulting in coarse estimates. In this proposal, we seek to improve these emission modules by combining state-of-the-art process modules for emissions, available data, and new optimization methods. In earth science problems, we often have substantial data and knowledge of processes in disparate systems, and thus we need to combine data and a general process level understanding into a model for projections of future climate that are as accurate as possible. The best methodologies for optimization of parameters in earth system models are still being developed. In this proposal we will develop and apply surrogate algorithms that a) were especially developed for computationally expensive simulations like CLM-ME/N models; b) were (in the earlier surrogate optimization Stochastic RBF) demonstrated to perform very well on computationally expensive complex partial differential equations in earth science with limited numbers of simulations; and, c) will be (as part of the proposed research) significantly improved both by adding asynchronous parallelism, early truncation of unsuccessful simulations, and the improvement of both serial and parallel performance by the use of derivative and sensitivity information from global and local surrogate approximations S(x) . The algorithm development and testing will be focused on the CLM-ME/N model application, but the methods are general and are expected to also perform well on optimization for parameter estimation of other climate models and other classes of continuous multimodal optimization problems arising from complex simulation models. In addition, this proposal will compile available datasets of emissions of methane, nitrous oxides and reactive nitrogen species and develop protocols for site level comparisons with the CLM-ME/N. Once the model parameters are optimized against site level data, the model will be simulated at the global level and compared to atmospheric concentration measurements for the current climate, and future emissions will be estimated using climate change as simulated by the CESM. This proposal combines experts in earth system modeling, optimization, computer science, and process level understanding of soil gas emissions in an interdisciplinary team in order to improve the modeling of methane and nitrogen gas emissions. This proposal thus meets the requirements of the SciDAC RFP, by integrating state-of-the-art computer science and earth system to build an improved earth system model.

54 ENVIRONMENTAL SCIENCES↗

Enhancements to the RAVEN code in FY16

The RAVEN code has been under development at the Idaho National Laboratory since 2012. Its main goal is to create a multi-purpose platform for the deploying of all the capabilities needed for Probabilistic Risk Assessment, uncertainty quantification, data mining analysis and optimization studies. RAVEN has currently reached a good level of maturity in terms of deployed state-of-art and advanced capabilities. The main subject of this report is to show the activities that have been recently accomplished: • Implementation of ensemble modeling for time-series, and • initial implementation of model validation for surrogate models, and • advanced visualization capability for topology based data analysis The development of ensemble modeling for time-series has been performed in order to begin tackling the needs of those RISMC applications that need to communicate 1-Dimensional information (e.g. power histories, etc.) among different models. In this document the implementation details and an application example is reported. The second subject of this report is about the initial development of methods, within the RAVEN framework, to assess the validity of the predictive capabilities of surrogate models. Indeed, after the construction of a surrogate tight to a certain physical model, it is crucial to assess the goodness of its representation, in order to be confident with its prediction. In this initial activity, a cross-validation technique has been employed. This report will highlight the implementation details and proof its correct implementation by an application example. The final subject of this report is about the implementation of advanced visualization capability in RAVEN, for interactive data analysis. Indeed, RAVEN offers several post-processing capabilities that can structurally decompose data extracted from experimental results offering both data clustering/partitioning and dimensionality reduction techniques. A disadvantage of the workflow available in RAVEN is that it treats these as black box operations and the user is expected to know specific information about their data including the number of partitions to expect in some cases or the “correct” parameter settings for a particular algorithm. In order to overcome this limitation, it has been added an interactive user interface that can be run in RAVEN to explore different parameter settings on the fly specifically for the topological post-processor and allows the user to explore the data interactively. The design of this user interface is generalizable to other algorithms and can be used as a model to generate more dynamic and user-friendly visualization capabilities within RAVEN.

97 MATHEMATICS AND COMPUTING↗

Dissolution of Used Nuclear Fuel Using a TBP/N-Paraffin Solvent

The dissolution of unirradiated used nuclear fuel (UNF) pellets pretreated for tritium removal was demonstrated using a tributly phosphate (TBP) solvent. Dissolution of pretreated fuel in TBP could potentially combine dissolution with two cycle of solvent extraction required for separating the actinides and lanthanides from other fission products. Dissolutions were performed using UNF surrogates prepared from both uranyl nitrate and uranium trioxide produced from the pretreatment process by adding selected actinide and stable fission product elements. In laboratory-scale experiments, the U dissolution efficiency ranged from 80-99+% for both the nitrate and oxide surrogate fuels. On average, 80% of the Pu and 50% of the Np and Am in the nitrate surrogate dissolved; however, little of the transuranic elements dissolved in the oxide form. The majority of the 3+ lanthanide elements dissolved. Only small amounts of Sr (0-1.6%) and Mo (0.1-1.7%) and essentially no Cs, Ru, Zr, or Pd dissolved.

##CAI suggestions##↗

PAGOSA Simulations for the LANL Blast Tube

A parameter study was performed for the Los Alamos National Laboratory Blast Tube to tune geometry to best replicate peak blast pressures from experiments. Charge placement, reflecting wall thickness, flange gap width, and reflecting wall stand-off were sampled using a Latin-Hypercube design to train a Radial Basis Function Network surrogate model. This surrogate model was then used with a Nelder - Mead Simplex minimization routine to minimize error for the surrogate predictions.

47 OTHER INSTRUMENTATION↗