SEARCH · Engineering Papers

Results for “SCALE”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

All sources NASA NTRS DOE OSTI

DOE OSTI2015

Effects of forcing time scale on the simulated turbulent flows and turbulent collision statistics of inertial particles

In this paper, we study systematically the effects of forcing time scale in the large-scale stochastic forcing scheme of Eswaran and Pope [“An examination of forcing in direct numerical simulations of turbulence,” Comput. Fluids 16, 257 (1988)] on the simulated flow structures and statistics of forced turbulence. Using direct numerical simulations, we find that the forcing time scale affects the flow dissipation rate and flow Reynolds number. Other flow statistics can be predicted using the altered flow dissipation rate and flow Reynolds number, except when the forcing time scale is made unrealistically large to yield a Taylor microscale flow Reynolds number of 30 and less. We then study the effects of forcing time scale on the kinematic collision statistics of inertial particles. We show that the radial distribution function and the radial relative velocity may depend on the forcing time scale when it becomes comparable to the eddy turnover time. This dependence, however, can be largely explained in terms of altered flow Reynolds number and the changing range of flow length scales present in the turbulent flow. We argue that removing this dependence is important when studying the Reynolds number dependence of the turbulent collision statistics. The results are also compared to those based on a deterministic forcing scheme to better understand the role of large-scale forcing, relative to that of the small-scale turbulence, on turbulent collision of inertial particles. To further elucidate the correlation between the altered flow structures and dynamics of inertial particles, a conditional analysis has been performed, showing that the regions of higher collision rate of inertial particles are well correlated with the regions of lower vorticity. Regions of higher concentration of pairs at contact are found to be highly correlated with the region of high energy dissipation rate.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

Evidence for Bolgiano-Obukhov scaling in rotating stratified turbulence using high-resolution direct numerical simulations

We report results on rotating stratified turbulence in the absence of forcing, with large-scale isotropic initial conditions, using direct numerical simulations computed on grids of up to $4096^3$ points. The Reynolds and Froude numbers are respectively equal to $Re=5.4\times 10^4$ and $Fr=0.0242$. The ratio of the Brunt-V\"ais\"al\"a to the inertial wave frequency, $N/f$, is taken to be equal to 5, a choice appropriate to model the dynamics of the southern abyssal ocean at mid latitudes. This gives a global buoyancy Reynolds number $R_B=ReFr^2=32$, a value sufficient for some isotropy to be recovered in the small scales beyond the Ozmidov scale, but still moderate enough that the intermediate scales where waves are prevalent are well resolved. We concentrate on the large-scale dynamics and confirm that the Froude number based on a typical vertical length scale is of order unity, with strong gradients in the vertical. Two characteristic scales emerge from this computation, and are identified from sharp variations in the spectral distribution of either total energy or helicity. A spectral break is also observed at a scale at which the partition of energy between the kinetic and potential modes changes abruptly, and beyond which a Kolmogorov-like spectrum recovers. Large slanted layers are ubiquitous in the flow in the velocity and temperature fields, and a large-scale enhancement of energy is also observed, directly attributable to the effect of rotation.

Bolgiano-Obukhov↗

DOE OSTI2016

The role of zonal flows in the saturation of multi-scale gyrokinetic turbulence

The 2D spectrum of the saturated electric potential from gyrokinetic turbulence simulations that include both ion and electron scales (multi-scale) in axisymmetric tokamak geometry is analyzed. The paradigm that the turbulence is saturated when the zonal (axisymmetic) ExB flow shearing rate competes with linear growth is shown to not apply to the electron scale turbulence. Instead, it is the mixing rate by the zonal ExB velocity spectrum with the turbulent distribution function that competes with linear growth. A model of this mechanism is shown to be able to capture the suppression of electron-scale turbulence by ion-scale turbulence and the threshold for the increase in electron scale turbulence when the ion-scale turbulence is reduced. The model computes the strength of the zonal flow velocity and the saturated potential spectrum from the linear growth rate spectrum. The model for the saturated electric potential spectrum is applied to a quasilinear transport model and shown to accurately reproduce the electron and ion energy fluxes of the non-linear gyrokinetic multi-scale simulations. The zonal flow mixing saturation model is also shown to reproduce the non-linear upshift in the critical temperature gradient caused by zonal flows in ion-scale gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Assessment of dynamic-screw-pinch-driven, current-scaled MagLIF target implosion performance using 3D magnetohydrodynamic simulations

Analytic studies and two-dimensional “clean” radiation-magnetohydrodynamic (rad-MHD) simulations employing dynamical similarity driver-target scaling prescriptions [Ruiz et al., Phys. Plasmas 30, 032708 (2023)] suggest that Magnetized Liner Inertial Fusion (MagLIF) target implosions can scale to > 10 MJ DT fusion yields when peak drive current is increased beyond 60 MA. We present results from three-dimensional (3D) rad-MHD simulations of similarity-scaled MagLIF target implosions at peak drive currents ranging from 15 to 40 MA. Simulations in this study suggest that magneto-Rayleigh–Taylor instability (MRTI) growth and feedthrough to the fuel region are more severe at higher drive current scales, which reduces the fusion yield compared to prior analytic and 2D clean simulation predictions. In contrast to standard MagLIF, simulations of current-scaled MagLIF target implosions driven by a dynamic screw pinch (DSP) demonstrate reduced MRTI feedthrough and greater fuel magnetization, resulting in improved thermonuclear performance and enhanced performance scaling with peak drive current. DSP drive enables additional scaling of the liner mass to increase liner radius but maintain implosion time, resulting in higher implosion velocities at the expense of increased susceptibility to MRTI. We present a current- and mass-scaled simulated DSP-MagLIF target implosion at the ∼ 40 MA peak current level that produces ignition scale performance, demonstrating a burn-averaged Lawson ignition parameter above unity and DT fusion yield above 1 MJ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

On large-scale dynamo action at high magnetic Reynolds number

We consider the generation of magnetic activity—dynamo waves—in the astrophysical limit of very large magnetic Reynolds number. We consider kinematic dynamo action for a system consisting of helical flow and large-scale shear. We demonstrate that large-scale dynamo waves persist at high Rm if the helical flow is characterized by a narrow band of spatial scales and the shear is large enough. However, for a wide band of scales the dynamo becomes small scale with a further increase of Rm, with dynamo waves re-emerging only if the shear is then increased. We show that at high Rm, the key effect of the shear is to suppress small-scale dynamo action, allowing large-scale dynamo action to be observed. We conjecture that this supports a general 'suppression principle'—large-scale dynamo action can only be observed if there is a mechanism that suppresses the small-scale fluctuations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2017

Multi-scale gyrokinetic simulations of an Alcator C-Mod, ELM-y H-mode plasma

High fidelity, multi-scale gyrokinetic simulations capable of capturing both ion (${k}_{\theta }{\rho }_{s}\sim { \mathcal O }(1.0)$) and electron-scale (${k}_{\theta }{\rho }_{e}\sim { \mathcal O }(1.0)$) turbulence were performed in the core of an Alcator C-Mod ELM-y H-mode discharge which exhibits reactor-relevant characteristics. These simulations, performed with all experimental inputs and realistic ion to electron mass ratio (${({m}_{i}/{m}_{e})}^{1/2}=60.0$) provide insight into the physics fidelity that may be needed for accurate simulation of the core of fusion reactor discharges. Three multi-scale simulations and series of separate ion and electron-scale simulations performed using the GYRO code are presented. As with earlier multi-scale results in L-mode conditions, both ion and multi-scale simulations results are compared with experimentally inferred ion and electron heat fluxes, as well as the measured values of electron incremental thermal diffusivities—indicative of the experimental electron temperature profile stiffness. Consistent with the L-mode results, cross-scale coupling is found to play an important role in the simulation of these H-mode conditions. Extremely stiff ion-scale transport is observed in these high-performance conditions which is shown to likely play and important role in the reproduction of measurements of perturbative transport. In conclusion, these results provide important insight into the role of multi-scale plasma turbulence in the core of reactor-relevant plasmas and establish important constraints on the fidelity of models needed for predictive simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2015

Is there scale-dependent bias in single-field inflation?

Scale-dependent halo bias due to local primordial non-Gaussianity provides a strong test of single-field inflation. While it is universally understood that single-field inflation predicts negligible scale-dependent bias compared to current observational uncertainties, there is still disagreement on the exact level of scale-dependent bias at a level that could strongly impact inferences made from future surveys. In this paper, we clarify this confusion and derive in various ways that there is exactly zero scale-dependent bias in single-field inflation. Much of the current confusion follows from the fact that single-field inflation does predict a mode coupling of matter perturbations at the level of f{sub NL}{sup local}; ≈ −5/3, which naively would lead to scale-dependent bias. However, we show explicitly that this mode coupling cancels out when perturbations are evaluated at a fixed physical scale rather than fixed coordinate scale. Furthermore, we show how the absence of scale-dependent bias can be derived easily in any gauge. This result can then be incorporated into a complete description of the observed galaxy clustering, including the previously studied general relativistic terms, which are important at the same level as scale-dependent bias of order f{sub NL}{sup local} ∼ 1. This description will allow us to draw unbiased conclusions about inflation from future galaxy clustering data.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Magnetorotational dynamo can generate large-scale vertical magnetic fields in 3D GRMHD simulations of accreting black holes

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $t_{\rm gen}$, and one for the large-scale scale advection, $t_{\rm adv}$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $\alpha$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $\alpha$-dynamos.

Astronomy & Astrophysics↗

DOE OSTI2018

Scale-dependent bias and bispectrum in neutrino separate universe simulations

Cosmic background neutrinos have a large velocity dispersion, which causes the evolution of long-wavelength density perturbations to depend on scale. As such, this scale-dependent growth leads to the well-known suppression in the linear theory matter power spectrum that is used to probe neutrino mass. In this paper, we study the impact of long-wavelength density perturbations on small-scale structure formation. By performing separate universe simulations where the long-wavelength mode is absorbed into the local expansion, we measure the responses of the cold dark matter (CDM) power spectrum and halo mass function, which correspond to the squeezed-limit bispectrum and halo bias. We find that the scale-dependent evolution of the long-wavelength modes causes these quantities to depend on scale and provide simple expressions to model them in terms of scale and the amount of massive neutrinos. Importantly, this scale-dependent bias reduces the suppression in the linear halo power spectrum due to massive neutrinos by 13% and 26% for objects of bias b ¯ = 2 and b ¯ &#187 1 , respectively. We demonstrate with high statistical significance that the scale-dependent halo bias cannot be modeled by the CDM and neutrino density transfer functions at the time when the halos are identified. This reinforces the importance of the temporal nonlocality of structure formation, especially when the growth is scale dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Non-Boussinesq subgrid-scale model with dynamic tensorial coefficients

A major drawback of Boussinesq-type subgrid-scale stress models used in large-eddy simulations is the inherent assumption of alignment between large-scale strain rates and filtered subgrid-stresses. A priori analyses using direct numerical simulation (DNS) data have shown that this assumption is invalid locally as subgrid-scale stresses are poorly correlated with the large-scale strain rates [J. Bardina, J. Ferziger, and W. Reynolds, Improved subgrid-scale models for large-eddy simulation, in Proceedings of the 13th Fluid and Plasmadynamics Conference, AIAA (1980); C. Meneveau and K. Katz, Scale-invariance and turbulence models for large-eddy simulation, Ann. Rev. Fluid Mech. 32, 1 (2000)]. In the present work, a new, non-Boussinesq subgrid-scale model is presented where the model coefficients are computed dynamically. Some previous non-Boussinesq models have observed issues in providing adequate dissipation of turbulent kinetic energy [e.g., Bardina et al., Proceedings of the 13th Fluid and Plasmadynamics Conference (1980); R. A. Clark, J. Ferziger, and W.C. Reynolds. Evaluation of subgrid-scale models using an accurately simulated turbulent flow, J. Fluid Mech. 91, 1 (1979); S. Stolz and N. A. Adams, An approximate deconvolution procedure for large-eddy simulation, Phys. Fluids 11, 1699 (1999)]; however, the present model is shown to provide sufficient dissipation using dynamic coefficients. Modeled subgrid-scale Reynolds stresses satisfy the consistency requirements of the governing equations for large-eddy simulation (LES), vanish in laminar flow and at solid boundaries, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The new model, referred to as the dynamic tensor-coefficient Smagorinsky model (DTCSM), has been tested in simulations of canonical flows: decaying and forced homogeneous isotropic turbulence, and wall-modeled turbulent channel flow at high Reynolds numbers. The results show favorable agreement with DNS data. It has been shown that DTCSM offers similar predictive capabilities as the dynamic Smagorinsky model for canonical flows. In order to assess the performance of DTCSM in more complex flows, wall-modeled simulations of high Reynolds number flow over a Gaussian bump (Boeing speed bump) exhibiting smooth-body flow separation are performed. Predictions of surface pressure and skin friction, compared against DNS and experimental data, show improved accuracy from DTCSM in comparison to existing static coefficient (Vreman) and dynamic Smagorinsky model. The computational cost of performing LES with this model is up to 15% higher than the dynamic Smagorinsky model.

42 ENGINEERING↗

DOE OSTI2023

Measuring scale-dependent shape anisotropy by coarse-graining: Application to inhomogeneous Rayleigh-Taylor turbulence

Here, we generalize the “filtering spectrum” to probe scales along different directions by spatial coarse-graining. This multidimensional filtering spectrum quantifies the spectral content of flows that are not necessarily homogeneous. From multidimensional spectral information, we propose a simple metric for shape anisotropy at various scales. The method is applied to simulations of two-dimensional (2D) and 3D Rayleigh-Taylor (RT) turbulence, which is inhomogeneous and anisotropic. We show that 3D RT has clear shape anisotropy at large scales with approximately 4:3 vertical to horizontal aspect ratio but tends toward isotropy at small scales as expected. In sharp contrast, we find that RT in 2D simulations, which are still the main modeling framework for many applications, is isotropic at large scales and its shape anisotropy increases at smaller scales where structures tend to be horizontally elongated. While this may be surprising, it is consistent with recent results that large-scale isotropy in 2D RT is due to the generation of a large-scale overturning circulation via an upscale cascade, while small-scale anisotropy is due to the stable stratification resultant from such overturning and the inefficient mixing in 2D.

42 ENGINEERING↗

DOE OSTI2025

Scaling Laws of Graph Neural Networks for Atomistic Materials Modeling

Atomistic materials modeling is a critical task with wide-ranging applications, from drug discovery to materials science, where accurate predictions of the target material property can lead to significant advancements in scientific discovery. Graph Neural Networks (GNNs) represent the state-of-the-art approach for modeling atomistic material data thanks to their capacity to capture complex relational structures. While machine learning performance has historically improved with larger models and datasets, GNNs for atomistic materials modeling remain relatively small compared to large language models (LLMs), which leverage billions of parameters and terabyte-scale datasets to achieve remarkable performance in their respective domains. To address this gap, we explore the scaling limits of GNNs for atomistic materials modeling by developing a foundational model with billions of parameters, trained on extensive datasets in terabytescale. Our approach incorporates techniques from LLM libraries to efficiently manage large-scale data and models, enabling both effective training and deployment of these large-scale GNN models. This work addresses three fundamental questions in scaling GNNs: the potential for scaling GNN model architectures, the effect of dataset size on model accuracy, and the applicability of LLM-inspired techniques to GNN architectures. Specifically, the outcomes of this study include (1) insights into the scaling laws for GNNs, highlighting the relationship between model size, dataset volume, and accuracy, (2) a foundational GNN model optimized for atomistic materials modeling, and (3) a GNN codebase enhanced with advanced LLM-based training techniques. Our findings lay the groundwork for large-scale GNNs with billions of parameters and terabyte-scale datasets, establishing a scalable pathway for future advancements in atomistic materials modeling.

Li, Chaojian [ORNL] (ORCID:0000000340309777)↗

DOE OSTI2021

Changes in Future Precipitation Mean and Variability across Scales

Changes in precipitation variability can have large societal consequences, whether at the short time scales of flash floods or the longer time scales of multiyear droughts. Recent studies have suggested that in future climate projections, precipitation variability rises more steeply than does its mean, leading to concerns about societal impacts. Here we evaluate changes in mean precipitation over a broad range of spatial and temporal scales using a range of models from high-resolution regional simulations to millennial-scale global simulations. Results show that changes depend on the scale of aggregation and involve strong regional differences. On local scales that resolve individual rainfall events (hours and tens of kilometers), changes in precipitation distributions are complex and variances rise substantially more than means, as is required given the well-known disproportionate rise in precipitation intensity. On scales that aggregate across many events, distributional changes become simpler and variability changes smaller. At regional scale, future precipitation distributions can be largely reproduced by a simple transformation of present-day precipitation involving a multiplicative shift and a small additive term. The “extra” broadening is negatively correlated with changes in mean precipitation: in strongly “wetting” areas, distributions broaden less than expected from a simple multiplicative mean change; in “drying” areas, distributions narrow less. Precipitation variability changes are therefore of especial concern in the subtropics, which tend to dry under climate change. Outside the tropics, variability changes are similar on time scales from days to decades (i.e., show little frequency dependence). This behavior is highly robust across models, suggesting it may stem from some fundamental constraint.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

SCALE Code System 6.2.2

The SCALE Code System is a widely used modeling and simulation suite for nuclear safety analysis and design that is developed, maintained, tested, and managed by the Reactor and Nuclear Systems Division (RNSD) of Oak Ridge National Laboratory (ORNL). SCALE provides a comprehensive, verified and validated, user-friendly tool set for criticality safety, reactor physics, radiation shielding, radioactive source term characterization, and sensitivity and uncertainty analysis. Since 1980, regulators, licensees, and research institutions around the world have used SCALE for safety analysis and design. SCALE provides an integrated framework with dozens of computational modules including 3 deterministic and 3 Monte Carlo radiation transport solvers that are selected based on the desired solution strategy. SCALE includes current nuclear data libraries and problem-dependent processing tools for continuous-energy (CE) and multigroup (MG) neutronics and coupled neutron-gamma calculations, as well as activation, depletion, and decay calculations. SCALE includes unique capabilities for automated variance reduction for shielding calculations, as well as sensitivity and uncertainty analysis. SCALE’s graphical user interfaces assist with accurate system modeling, visualization of nuclear data, and convenient access to desired results. SCALE 6.2 represents one of the most comprehensive revisions in the history of SCALE, providing several new capabilities and significant improvements in many existing features.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2018

Extreme-Scale Bayesian Inference for Uncertainty Quantification of Complex Simulations

Uncertainty quantification (UQ)—that is, quantifying uncertainties in complex mathematical models and their large-scale computational implementations—is widely viewed as one of the outstanding challenges facing the field of CS&E over the coming decade. The EUREKA project set to address the most difficult class of UQ problems: those for which both the underlying PDE model as well as the uncertain parameters are of extreme scale. In the project we worked on these extreme-scale challenges in the following four areas: 1. Scalable parallel algorithms for sampling and characterizing the posterior distribution that exploit the structure of the underlying PDEs and parameter-to-observable map. These include structure-exploiting versions of the randomized maximum likelihood method, which aims to overcome the intractability of employing conventional MCMC methods for solving extreme-scale Bayesian inversion problems by appealing to and adapting ideas from large-scale PDE-constrained optimization, which have been very successful at exploring high-dimensional spaces. 2. Scalable parallel algorithms for construction of prior and likelihood functions based on learning methods and non-parametric density estimation. Constructing problem-specific priors remains a critical challenge in Bayesian inference, and more so in high dimensions. Another challenge is construction of likelihood functions that capture unmodeled couplings between observations and parameters. We will create parallel algorithms for non-parametric density estimation using high dimensional N-body methods and combine them with supervised learning techniques for the construction of priors and likelihood functions. 3. Bayesian inadequacy models, which augment physics models with stochastic models that represent their imperfections. The success of the Bayesian inference framework depends on the ability to represent the uncertainty due to imperfections of the mathematical model of the phenomena of interest. This is a central challenge in UQ, especially for large-scale models. We propose to develop the mathematical tools to address these challenges in the context of extreme-scale problems. 4. Parallel scalable algorithms for Bayesian optimal experimental design (OED). Bayesian inversion yields quantified uncertainties in the model parameters, which can be propagated forward through the model to yield uncertainty in outputs of interest. This opens the way for designing new experiments to reduce the uncertainties in the model parameters and model predictions. Such experimental design problems have been intractable for large-scale problems using conventional methods; we will create OED algorithms that exploit the structure of the PDE model and the parameter-to-output map to overcome these challenges. Parallel algorithms for these four problems were created, analyzed, prototyped, implemented, tuned, and scaled up for leading-edge supercomputers, including UT-Austin’s own 10 petaflops Stampede system, ANL’s Mira system, and ORNL’s Titan system. While our focus is on fundamental mathematical/computational methods and algorithms, we will assess our methods on model problems derived from several DOE mission applications, including multiscale mechanics and ice sheet dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Final Report, University of California Merced: Uranium and strontium fate in waste-weathered sediments: Scaling of molecular processes to predict reactive transport

Objectives of the Project: 1. Determine the process coupling that occurs between mineral transformation and contaminant (U and Sr) speciation in acid-uranium waste weathered Hanford sediments. 2. Establish linkages between molecular-scale contaminant speciation and meso-scale contaminant lability, release and reactive transport. 3. Make conjunctive use of molecular- to bench-scale data to constrain the development of a mechanistic, reactive transport model that includes coupling of contaminant sorption-desorption and mineral transformation reactions. Hypotheses Tested: Uranium and strontium speciation in legacy sediments from the U-8 and U-12 Crib sites can be reproduced in bench-scale weathering experiments conducted on unimpacted Hanford sediments from the same formations; Reactive transport modeling of future uranium and strontium releases from the vadose zone of acid-waste weathered sediments can be effectively constrained by combining molecular-scale information on contaminant bonding environment with grain-scale information on contaminant phase partitioning, and meso-scale kinetic data on contaminant release from the waste-weathered porous media; Although field contamination and laboratory experiments differ in their diagenetic time scales (decades for field vs. months to years for lab), sediment dissolution, neophase nucleation, and crystal growth reactions that occur during the initial disequilibrium induced by waste-sediment interaction leave a strong imprint that persists over subsequent longer-term equilibration time scales and, therefore, give rise to long-term memory effects. Enabling Capabilities Developed: Our team developed an iterative measure-model approach that is broadly applicable to elucidate the mechanistic underpinnings of reactive contaminant transport in geomedia subject to active weathering.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2019

Developing a Macro-scale SiC-cladding Behavior Model Based on Localized Mechanical and Thermal Property Evaluation on Pre- and Post-Irradiation SiC-SiC Composites (Final Report)

Silicon carbide is being investigated for accident tolerant fuel cladding applications due to its good neutronic performance, high temperature strength, exceptional stability under irradiation, and reduced oxidation compared to Zircaloy under accident conditions. The development and investigation of these materials is particularly important in the light of the Fukushima event and subsequent emphasis in DOE on accident tolerant fuel (ATF) concepts. In this work, we propose to develop and improve upon existing small-scale mechanical and thermal characterization methods, including micro-cantilever bend tests and fiber push-out tests, and time-domain thermo-reflectance measurements. These techniques will be applied to evaluate micro-scale properties of SiC-SiC composite constituents (matrix, fibers, interphase). The results will be coupled with meso-scale fiber and void structural information and be used as input to develop a comprehensive, finite element-based model based on constituent properties. The goal of this multi-scale model is to be used to predict anisotropic mechanical and thermal cladding behaviors. The model will be benchmarked against measurements of bulk properties and validated for radiation-damaged materials by utilizing irradiated SiC-SiC composites available from General Atomics. The characterization approach developed under this program will have particular application in the property prediction of reactor irradiated materials, as only small volumes need to be tested, minimizing associated hazards and costs. The micro-scale test methods developed here will be used along with meso-scale structural data to inform a finite element model of the full SiC-based cladding system. This model will capture fabrication and irradiation effects and importantly include for the first-time separate material data for each system component. The incorporation of these micro-scale effects will lead to a more accurate model of SiC-SiC composite behavior, enabling further advancements in the development and design of SiC-based materials for improved fuel cladding performance. The aim of this work is to develop advanced localized material characterization techniques to directly measure mechanical and thermal properties of the individual constituents of SiC-based claddings at the relevant micro-scale. SiC-SiC composites will be evaluated before and after irradiation and these results will be coupled with macro-scale properties and microstructural information in order to provide the input parameters for a comprehensive finite element model, which will be developed in this program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗