Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameterization”

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,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Iodine Waste Form Corrosion: PNNL Summary Testing Report

A test method to correlate the laboratory corrosion behavior of iodine-containing waste forms (IWFs) to a conceptual long-term corrosion models is required. Several candidate methods have been evaluated previously, and each method presents unique challenges and limitations. Dynamic leach tests have emerged as the preferred candidate test method to evaluate IWFs and parameterize an eventual long-term corrosion model. This report was prepared in support of this effort to define the test suite to evaluate the corrosion of IWFs, with the goal of defining and parameterizing a conceptual long-term corrosion model. A review was presented highlighting the high degree of variability between the IWF corrosion data that is available in the literature and the need for a standardized test suite. The targeted test suite is composed of a single interval screening test, a dynamic leach test, and post-test solid characterization of the monolithic IWF samples. The data presented serves two purposes: to expand the range of IWFs evaluated with the test suite to assess method variability and to further understand the corrosion of IWFs to best inform conceptual long-term corrosion model development. A series of IWF types were included in the testing, including hot isostatically pressed Ag-mordenite, silver functionalized silica-based aerogels, iodosodalites (hot isostatic pressing, sintered, glass bonded ceramic), iodoapatite, tellurite glass and glass composite materials.

42 ENGINEERING↗
DOE OSTI2021

Facilitating better and faster simulations of aerosol-cloud interactions in Earth system models

Focal Area(s): 1. Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising a hierarchy of models. 2. Insight gleaned from complex data (both observed and simulated) using AI, big data analytics, and other advanced methods, including explainable AI and physics- or knowledge-guided AI. Science Challenge: One major challenge that Earth system models (ESMs) face in providing credible prediction of the Earth system and its water cycle characteristics (e.g., mean state, variability, and extreme events) is to accurately simulate aerosol-cloud interactions (ACI). The physical, chemical, and dynamical processes affecting ACI are extremely complex and they range from nanoscale to planetary scale. In each model development cycle, scientists spend significant efforts investigating model deficiencies and uncertainties associated with aerosols (e.g., emissions, chemical processes, aerosol microphysics, and transport) and clouds (e.g., macrophysics, microphysics, turbulence, and large-scale circulation) in order to develop improved treatments. However, despite decades of active research, ACI is still a major source of uncertainty in climate projections, even though great progress has been made. Specific scientific challenges include: (i) Parameterizations are developed based on limited data; (ii) The complexity of a parameterization required for accurate predictions is not understood; (iii) Incomplete and unknown physics leads to errors in the fully coupled Earth system; and (iv) Complex physics is computationally too expensive to employ in ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Structurally flexible cloud microphysics, observationally constrained at all scales via ML-accelerated Bayesian inference

We discuss the challenge of developing observationally informed parameterizations of microphysics for use at a hierarchy of modeling scales. Our proposed approach is applicable to any domain that suffers from a two-fold parameterization problem, where physical processes are not resolved at the model scale (the first problem), and those processes are uncertain at any scale (the second problem). For such problems, a physical approach facilitates modeling across scales, as well as systematic observational inference accelerated by machine learning (ML) surrogate models, and quantification of physical uncertainties.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Wind Energy Forecasting with the Weather Research and Forecasting Model

This was a collaborative effort between Lawrence Livermore National Security, LLC as manager and operator of Lawrence Livermore National Laboratory (LLNL) and Siemens Energy, Inc. (Siemens) to develop a wind resource forecasting tool. LLNL was to develop an independent high-resolution mesoscale modeling capability forecasting tool that could be implemented in conjunction with existing wind farm control and monitoring software to provide forecasting of wind resources using local observations of winds and temperature. Research with LLNL’s state-of-the-art large-eddy simulation meteorological prediction model, based on the community WRF model and innovative turbulence parameterizations, would improve that model’s applicability to large wind farms offshore and in complex terrain. The modeling capability would include uncertainty quantification. Finally, the application of the modeling tool and existing global climate change predictions would enable the delineation of the likely effects of climate change on wind resources. Siemens was to provide high time resolution hub-height wind speed and other meteorological data streams, including temperature profiles from wind farms, for LLNL to incorporate into the modeling system, to validate and tune this forecasting model for their locations of interest. These data streams would also be used for longer-term studies of correlations of wind resources to climate oscillations to indicate how long-term climate change trends may affect the available wind resource. Siemens would also provide information and observations of turbine wakes for incorporation into the modeling tool. By implementing state-of-the-art turbulence parameterizations into a simulation model and/or ensembles of simulation models, and by integrating real-time hub height wind speed and other meteorological datastreams from wind farms into that model or ensemble of models, LLNL would develop a forecasting tool that could be implemented by Siemens as an add-on to existing wind farm control and monitoring software to provide owners with useful resource forecasting. The desired outcome was that the accuracy level of the output would be sufficient to substantiate power output commitments. The final deliverable for this work would consist of a document outlining the algorithms and software tools that could be integrated into Siemens Wind Park Supervisor.

17 WIND ENERGY↗
DOE OSTI2021

Aerosol Chemical Speciation Monitor (ACSM) Composition-Dependent Collection Efficiency (CDCE) Value-Added Product Report

Aerosol particles influence the Earth’s radiation balance directly by absorbing and scattering light and indirectly by influencing cloud formation, properties, and lifetimes. Measurements of aerosol particle optical properties, mass loading, size distributions, microphysical properties, cloud formation properties, and chemical composition are important for understanding the aerosol life cycle and for validating earth system models that predict these quantities. The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) user facility’s Aerosol Observing System (AOS) is a highly instrumented platform that measures many of these aerosol properties in situ. The Aerodyne aerosol chemical speciation monitor (ACSM) is a baseline instrument deployed in the AOS. The ACSM provides a quantitative measurement of aerosol particle chemical composition for non-refractory (operationally defined as components that evaporate on the 600ºC vaporizer) aerosol components in real time. Standard output is the mass concentration of particulate organics, nitrate, sulfate, chloride, and ammonium. One well-known limitation to the accuracy of the ACSM data is in evaluating the fraction of the ambient aerosol particles that are detected by the instrument. This quantity is referred to as the collection efficiency (CE) and is often less than unity. This is attributed to particles rebounding after impaction onto the heated vaporizer rather than being trapped, volatilized rapidly, and detected. Other factors, such as divergence of the aerosol particle beam, may also impact CE, but the particle rebound effect is the dominant factor. Scientists will sometimes assume CE = 0.5 (i.e., one half of sampled particles are detected) for ambient particles collected during field missions. However, parameterizations have been developed that express CE as a function of the measured chemical composition, referred to as the composition-dependent collection efficiency (CDCE). The physical explanation for a CDCE, supported by laboratory studies, is that rebound from the vaporizer is a function of the particle phase, with liquid-like particles “sticking” to the vaporizer and solid-like or crystalline particles “bouncing” from the vaporizer. Thus, particles with compositions are liquids-like under the measurement conditions (e.g., certain organics, acidic particles, particles enriched in nitrate) have CE close to unity while crystalline or solid particles (e.g., deliquesced ammonium sulfate) have lower CE. This value-added product (VAP) implements the procedure described by Middlebrook et al. (2012) to correct the ACSM data for the composition-dependent collection efficiency. Applying this parameterization improves the accuracy of the ACSM data and brings them into better agreement with other co-located aerosol measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assessing the drivers of Isoprene SOA: laboratory studies, field observations and modeling (Final Technical Report)

This report describes the research results and products stemming from the funding award DE-SC0018221. The main goals of this project were to develop a quantitative understanding of the factors driving secondary organic aerosol (SOA) formation from biogenic volatile organic compounds, namely isoprene and monoterpenes, using a combination of laboratory measurements, field data analysis, and modeling. Specifically, we aimed to: 1) develop an observationally constrained volatility distribution of trace gases produced from the oxidation of biogenic hydrocarbons that can explain the formation of SOA from in situ observations made during DOE ASR/ARM field campaigns such as HI-SCALE and BAECC 2) use new and existing (e.g. DOE ASR funded SOAFFEE experiments) laboratory chamber experiments to develop detailed parameterizations of isoprene-derived epxoy diol (IEPOX) reactive uptake, including product branching, product volatility, and solubility constants, as well as SOA formation generally from the formation of highly oxygenated organic molecules (HOM) formed from isoprene and monoterpene photo-oxidation that can be directly incorporated into models. 3) collaborate with PNNL modeling teams to incorporate new parameterizations into detailed box models, such as MOSAIC, and regional or Earth System models such as WRF.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ellipsoidal Fitting Methodology for Defect Clusters in Gallium Arsenide

In assessing the initial spatial distribution of defects from neutron or heavy ion irradiation, it is useful to have a reliable, automated, and fast-running tool to evaluate characteristic metrics such as the number of sub-clusters or the overall cluster volume. The latter metric, for instance, can be utilized to estimate a reference neutron fluence level at which inter-cluster interaction effects begin to become significant. This paper details a methodology to fit an arbitrarily complex defect map with a set of ellipsoids (one per identified sub-cluster) in which the constituent defects of a sub-cluster are determined using fuzzy degree-of-membership analysis. Specifically, a parameterized model is developed for point defects in gallium arsenide. Cluster volume calculations based on the model are compared against convex hull and single- ellipsoid representations. Results show that the parameterized sub-cluster model begins to deviate from the two reference models at a recoil energy of about 100 keV in GaAs, with the convex hull and single-ellipsoid representations increasingly overestimating the volume thereafter.

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

Small field campaign: aerosol – ice formation closure pilot study. Final Report

Prediction of atmospheric ice formation from aerosol particles by heterogeneous nucleation represents one of the grand challenges in atmospheric science. Our insufficient predictive understanding of primary ice formation is the reason that climate models typically do not include heterogeneous ice nucleation with subsequent effects on climate uncertainty. Mixed-phase clouds, where supercooled water droplets and ice crystals coexist play globally an important role regulating climate. This is especially the case for the Arctic region that experiences the greatest warming due to climate change compared to other regions in the world. Immersion freezing initiated by ice-nucleating particles (INPs) in supercooled water droplets is recognized as the dominant primary ice formation pathway in mixed-phase cloud regimes. For this reason, it is crucial to evaluate our capability to predict immersion freezing for a given ambient aerosol population. The goal of this project is conducting a field-based pilot study at the U.S. DOE Atmospheric Radiation Measurement (ARM) user facility at Southern Great Plains (SGP) to evaluate our capability to predict the number concentration of aerosol particles that serve as INPs in the immersion freezing mode. Successful prediction of INP number concentrations is also termed “closure”. This field-observational approach represents a first-of-its kind attempt of an aerosol–ice formation closure study (AEROICESTUDY). Very few closure studies related to INPs have been conducted, and to our knowledge, none using robust size-resolved ambient aerosol composition measurements as a starting point. Achievement of aerosol–ice formation closure relies on our ability to characterize the ambient aerosol population with respect to particles size and composition and to determine INP number concentrations for specified freezing temperatures. This requires numerous online and offline instrumentation resulting in this pilot field campaign being a multi-institutional and community-collaborative effort. We chose the ARM SGP megasite for this first aerosol-ice formation closure pilot study due to its significant measurement capabilities available to obtain detailed physical characterization of the local aerosol population including size distribution, mass loading, and chemical composition of non-refractory aerosol particles. The overall objective of this project is to identify ice nucleation parameterizations that produce the most robust predictions of INP numbers and thus are best suited to be included in cloud and climate models. This objective includes the following goals for this field and laboratory-based project: i) What are the crucial aerosol physicochemical property measurements needed to accurately guide ice nucleation representations in models and long-term INP measurements? ii) What level of parameter details needs to be known to achieve aerosol–ice formation closure? iii) What are the leading causes for climate model bias in INP predictions? We found that the advances in our understanding of immersion freezing garnered over the last 20 years allowed us to yield partial and full closures of atmospheric immersion freezing from ambient aerosol particles. When the aerosol population is physicochemically complex and parameterizations for representative INP types are not yet available, we still struggle to accurately predict INP number concentrations. This project clearly demonstrates that with more laboratory and field measurements that are accompanied by particle composition analysis, the necessary datasets to achieve aerosol–ice formation closure for various locations will emerge, thus providing a robust foundation for guiding the representation of INPs in cloud and climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

FY2022 Status Update: A Probabilistic Model for Stress Corrosion Cracking of SNF Dry Storage Canisters

Understanding the potential risk of stress corrosion cracking of spent nuclear fuel dry storage canisters has been identified as a knowledge gap for determining the safety of long-term interim storage of spent nuclear fuel. To address this, the DOE is funding a multi-lab DOE effort to understand the timing, occurrence, and consequences of potential canister SCC. Sandia National Laboratories has developed a probabilistic model for canister penetration by SCC. This model has been continuously updated at SNL since 2014. Model uncertainties are treated using a nested loop structure, where the outer loop accounts for uncertainties due to lack of data and the inner aleatoric loop accounts for uncertainties due to variation in nature. By separating uncertainties into these categories, it is possible to focus future work on reducing the most influential epistemic uncertainties. Several experimental studies have already been performed to improve the modeling approach through expanded process understanding and improved model parameterization. The resulting code is physics-based and intended to inform future work by identifying (1) important modeling assumptions, (2) experimental data needs, and (3) necessary model developments. In this document, several of the sub-models in the probabilistic SCC model have been exercised, and the intermediate results, as the model progresses from one sub-model to the next, are presented. Evaluating the sub-models in this manner provides a better understanding of sub-model outputs and has identified several unintended consequences of model assumptions or parameterizations, requiring updates to the modeling approach. The following updates have been made, and future updates have been identified.

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

Toward Better Understanding of Microphysical Processes and Resulting Precipitation Physics: A Merger Of Observations and Cloud Models (Final Report)

A large database of global disdrometer observations and a diverse set of simulations from the Regional Atmospheric Modeling System (RAMS) were compared using principal component analysis (PCA) in order to better understand warm and ice-based precipitation processes. The analysis demonstrated that six distinct precipitation groups (PGs) with common characteristics were revealed in both the observations and model simulations. These PGs were defined by similar co-variability of rain parameters. However, the model showed large concentrations of drops with a given size compared to the observations. Detailed investigations determined the parameterization of drop breakup was forcing the drops to an equilibrium size with a stronger constraint than was indicated by the observations. A series of sensitivity studies permutating the rain shape parameter in the RAMS 2-moment bulk scheme demonstrated how the assumed shape of the rain distribution influences the microphysics as well as precipitation characteristics. Results were compared to the same simulations performed with a bin microphysics scheme where the rain shape parameter is freely evolving. In the bin scheme, a wide range of shape parameters were found, with horizontal and vertical variability, in addition to being a function of storm lifetime. The results of this study highlight the limitations of a fixed assumed rain distribution in 2-moment microphysics schemes. Model simulations were used to probe the microphysical origins of the six PGs. Rain budgets from a supercell case showed that several of the groups hypothesized to be associated with strong ice microphysical processes had significant contributions from melting hail, but the complexity of ice and warm-rain processes in the supercell precluded definitive determination of microphysical origins for many of the PGs. A case study from the Mid-latitude Continental Convection and Clouds Experiment (MC3E) compared the model PGs to disdrometer observations, and demonstrated the spatial and temporal variability of the PGs which is not possible from disdrometer point measurements alone. The analyses performed during this project demonstrated how statistical analysis can provide a robust method for comparing model simulations and observations, which ultimately resulted identification of limitations in some current microphysics parameterizations as well as insights into precipitation variability which ultimately benefits both observations and modeling efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

HPC for Fission Modeling in Support of Nuclear Data [Slides]

The presentation discusses the background of the High-Performance Computing (HPC) division at the Los Alamos National Laboratory (LANL), with emphasis on its computational needs and computational techniques. It also discusses the highlights of the HPC and the results of calculations, including energy sharing parameterized from DFT calculations, parameterization used in CGMF calculations, angular momentum from microscopic calculations, and correlations between FF angular momenta. The presentation concludes with an outlook including what remains to be done in the HPC division and how it will benefit nuclear data, specifically how it will provide a better understanding of the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Illuminating the pathways to carbon liberation: a systems approach to characterizing the consequential unknowns of carbon transformation and loss from thawing permafrost peatlands (Final Report)

The IsoGenie3 Project delivered new systems-level insights into carbon cycling in thawing permafrost landscapes, with an emphasis on methane and carbon dioxide emissions. From >200 samples from the site collected over a decade, co-analyzed for geochemistry and microbiology, the team recovered ~1,500 assembled microbial genomes and ~1,900 viral population genomes, revealing appreciable genetic novelty - from a new highly abundant bacterial phylum, to novel methane consumers and their activities, to rampant viral novelty. IsoGenie3 linked these organisms to carbon compound transformations (which define the cycling of organic matter in soils, and the loss of the greenhouse gases carbon dioxide and methane), and saw that the microbes at each stage of permafrost thaw had different genetic potential to degrade categories of compounds, expressed that genetic potential differently, and actually transformed carbon compounds into greenhouse gases in different ways. IsoGenie 3 identified that some of the thaw-stage differences were due to plant-microbiome relationships; the plant species across the thaw gradient contributed different carbon compounds into the soil, and hosted distinct microbiota (differing among parts of plants as well as species). Lastly, microbes in the saturated post-thaw conditions appeared likely to contribute to the mobilization and toxification of mercury released during thaw. In parallel with ongoing field sampling and analysis, hypotheses arising from field observations were tested via lab incubation experiments. When communities are taken out of their native habitats, they behave differently, and the team first rigorously quantified the magnitude of this effect on microbiome composition and functional capacity, organic matter composition, and gas production; overall the main system processes were maintained in the lab incubations under the conditions tested. Further, the microbial data could inform geochemical reaction network models of those processes. Then, the team ran experiments with additions of compounds, varying temperature, and “live” vs. “dead” peat (the latter having been gamma irradiated, with a few additional variants to control for methodological artifacts). From these, we (a) determined the importance of plant-derived soluble phenolic compounds in bogs’ extraordinary recalcitrance of organic matter, and carbon gas emissions skewed to carbon dioxide; (b) proposed an abiotic ‘tanning’ mechanism, which could contribute to Sphagnum’s inhibitory effect on anaerobic decomposition through alteration of N availability. IsoGenie3 illuminated longer-term and landscape-scale interactions of permafrost thaw and carbon cycling, advancing knowledge of the drivers of methane dynamics not only across in the permafrost-associated peatland (where hydrology and plant communities dictate microbiomes) but also their interconnected lakes (where sediment carbon quality and resident microbiota are determined by position within lake, and lake features). By leveraging observations of site methane dynamics extending well before this project, the team was able to construct a 44-year portrait of the interplay of permafrost thaw, hydrology, vegetation dynamics, and carbon gas emissions, and the doubling of the fully-thawed fens over this time. From the detailed study of this focal site, IsoGenie3 also aimed to improve model representation of these kinds of sites and processes. To improve predictions of methane transformations, we incorporated acetate and isotope dynamics into the ‘DNDC’ biogeochemistry model. In addition, recovered genomes were grouped into ‘functional groups’, i.e. the genomes that perform a specific function of interest, then used to parameterize maximum growth rate and optimum growth temperature (via signatures in their sequence composition) for the BioCrunch model. The BioCrunch model was then in turn used to test the impact of increasing functional resolution of the microbes, on the carbon gas emissions. Lastly for modeling, the ecosys model was parameterized from the microbial and other data, and used to evaluate drivers of e.g. change in methane emissions. Finally, this project also led to the development of a range of new methods and tools, a new metric of organic matter decomposability, as well as a graph-database solution to multidisciplinary data storage and querying. This project’s ongoing analyses at our focal site also contributed to broader advancements in understanding elements of genetic plasticity and methane metabolism, climate change microbiology and community assembly, global peatland geochemistry and Arctic lakes’ roles in climate feedbacks.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assimilating Scanning Radar Data into Cloud-Scale Models

The aim of this project is to improve the initial conditions of large eddy simulations (LES) of deep convective storms by assimilating scanning radar reflectivity and radial velocity data as well as local environmental observations. Accomplishing this goal often results in more accurate numerical forecasts of clouds and precipitation than if the background initial meteorological conditions are solely relied upon for convection initiation and intensification, as is typical of most operational and research models. Conducted in a storm-scale LES framework, these more accurate numerical analyses can serve as tool to better understand and parameterize land-atmosphere coupling, boundary layer turbulence, and cloud processes, each used for model parameterization development of convection. Our work examined the sensitivity of the numerical representation of cloud-scale wind and microphysical structures to: i) the assimilation of radiosonde observations characterizing the near-storm inflow, ii) the assimilation of radar data of varying resolution for model simulations with grid spacings of 1-km and 0.25-km, and iii) to the assimilation of three-dimensional plan-projection-indicator (PPI) and range-height indicator (RHI) radar volumes. Few studies have explored deep convective storm analyses and forecasts produced by assimilation of radar observations to cloud-scale LES. Furthermore, to our knowledge, no storm-scale DA study has explored the value of assimilating high-resolution sector RHI volumes to determine if they can improve the retrieval of the vertical kinematic and microphysical storm structure aloft over the use of traditionally utilized PPI scans.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Size-resolved particle and black carbon deposition over the cryosphere (Final Report)

Our project focused on providing observational constraints and investigation of aerosol deposition by performing the first unambiguous direct eddy flux covariance measurements of aerosol dry deposition over the cryosphere. The project aimed to make eddy covariance flux measurements of dry deposition and off-line wet deposition measurements of black carbon containing aerosol plus eddy covariance measurements of size-resolved accumulation mode scattering aerosol over the DOE ARM AMF3 Oliktok Point site in Fall 2020. We aimed to compare our measurements to model parameterizations to constrain uncertainties and systematic bias, and improve deposition parameterizations in models. Due to disruptions in field work from COVID-19, our work was delayed and moved to an alternate site on the North Slope of Alaska. This report summarizes the key findings from the project, including data analysis from previous field projects and from the North Slope of Alaska.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Using ARM Data to Retrieve Entrainment Rates in Stratocumulus Cloud Systems (Final Report)

A stratocumulus-topped boundary layer (STBL) is a shallow convecting layer in which both the updraft and the downdraft branches of the convective circulation are saturated (cloudy) in the upper part of the boundary layer. The convective circulation is vigorous, yet the cloud-top entrainment velocity, $w_e$, is small and difficult to measure. Differences between currently used parameterizations of $w_e$ in GCMs can lead to STBLs that differ by as much as a factor of two in climatologically important properties such as liquid water path and boundary layer depth. We have developed a new method to retrieve entrainment rates in STBLs which uses vertical velocity retrievals from surface-based cloud radar, in conjunction with other surface-based measurements. These new estimates will allow a more comprehensive evaluation of parameterizations of entrainment rate for STBLs. The new method will also provide ground-based estimates of entrainment in STBLs over long time periods. We use retrievals of vertical velocity and liquid water content profiles to estimate the turbulent liquid water flux, which is then used in combination with observed cloud-top jumps of temperature, water vapor, and calculated radiative flux profiles to estimate the entrainment rate. We present results for three days for which observational data was not significantly contaminated by drizzle. The resulting average entrainment rates are 0.58 cm/s, 1.51 cm/s, and 0.31 cm/s.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Validation of High Resolution Atmospheric and Dispersion Simulations

Weather models typically operate on averages of 1 km or more. This averages out fine scale turbulence. This results in smooth laminar flows and rapid vertical motions are usually averaged out. In this study, those fine scaled features were simulated for a 2013 field study of transport and diffusion of a tracer gas. Parameterized subgrid turbulence was required to produce dilution representative of observations. Despite this, the simulation with the finest grid and contained parameterized subgrid turbulence produced dilution that was consistent with observations.

THOMAS, ANDREW↗
DOE OSTI2024

Controlling radioisotope proportions when randomly sampling from Dirichlet distributions in PyRIID

As machine learning models for radioisotope quantification become more powerful, likewise the need for high-quality synthetic training data grows as well. For problem spaces that involve estimating the relative isotopic proportions of various sources in gamma spectra it is necessary to generate training data that accurately represents the variance of proportions encountered. In this report, we aim to provide guidance on how to target a desired variance of proportions which are randomly when using the PyRIID Seed Mixer, which samples from a Dirichlet distribution. We provide a method for properly parameterizing the Dirichlet distribution in order to maintain a constant variance across an arbitrary number of dimensions, where each dimension represents a distinct source template being mixed. We demonstrate that our method successfully parameterizes the Dirichlet distribution to target a specific variance of proportions, provided that several conditions are met. This allows us to follow a principled technique for controlling how random mixture proportions are generated which are then used downstream in the synthesis process to produce the final, noisy gamma spectra.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Quantifying model prediction sensitivity to model-form uncertainty

Computational and mathematical models are essential to understanding complex systems and phenomena. However, when developing such models, limited knowledge and/or resources necessitates the use of simplifying assumptions. It is therefore crucial to quantify the impact of such simplifying assumptions on the reliability and accuracy of resulting model predictions. This work develops a first-of-its-kind approach to quantify the impact of physics modeling assumptions on predictions. Here, we leverage the emerging field of model-form uncertainty (MFU) representations, which are parameterized modifications to modeling assumptions, in combination with grouped Sobol’ indices to quantitatively measure an assumption’s importance. Specifically, we compute the grouped Sobol’ index for the MFU representation’s parameters as a single importance measure of the assumption for which the MFU representation characterizes uncertainty. To ensure this approach is robust to the subjective choice of how to parameterize a MFU representation, we establish bounds for the difference between sensitivity results for two different MFU representations based on differences in model prediction statistics. The capabilities associated with this approach are demonstrated on three exemplar problems: an upscaled subsurface contaminant transport problem, ablation modeling for hypersonic flight, and nuclear waste repository modeling. We found that our grouped approach is able to assess the impact of modeling assumptions on predictions and offers computational advantages over classical Sobol’ index computation while providing more interpretable results.

97 MATHEMATICS AND COMPUTING↗