Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Improved treatment of multi-material cells in thermal radiation transport codes

High-energy-density physics simulations with non-conformal meshes of materials require a multi-material (MM) closure that affects the thermal radiation transport (TRT). We propose a set of novel closures that work for an arbitrary number of materials, both grey and multigroup energy discretizations, and any angular discretization (such as Sn, IMC, or diffusion). For each spatial cell, our closures let each species (ion and electron) of each material have its own temperature, density, and internal energy, but use a single radiation distribution that interacts with all materials within the cell. Our closures maintain energy conservation, do not incur increased computational cost in the TRT solve itself, are compatible with single-material TRT solvers, and do not make any temperature-equilibrium assumptions. Here we test our closures on a wide range of increasingly realistic problems and find them to be robust.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Spatiotemporally dynamic implicit large eddy simulation using machine learning classifiers

In this article, we utilize machine learning to dynamically determine if a point on the computational grid requires implicit numerical dissipation for large eddy simulation (LES). The decision making process is learnt through a priori training on quantities derived from direct numerical simulation (DNS) data. In particular, we compute eddy-viscosities obtained through the coarse-graining of DNS quantities and utilize their projection onto a Gaussian distribution to categorize areas that may require dissipation. If our learning determines that closure is necessary, an upwinded scheme is utilized for computing the non-linear Jacobian. In contrast, if it is determined that closure is unnecessary, a symmetric and second-order accurate energy and enstrophy preserving Arakawa scheme is utilized instead. This results in a closure framework that precludes the specification of any model-form for the small scale contributions of turbulence but deploys an appropriate numerical dissipation from explicit closure driven hypotheses. Here, this methodology is deployed for the Kraichnan turbulence test-case and assessed through various statistical quantities such as angle-averaged kinetic energy spectra and vorticity structure functions. Our framework thus establishes a link between the use of explicit LES ideologies for closure and numerical dissipation-based modeling of turbulence leading to improved statistical fidelity of a posteriori simulations.

42 ENGINEERING↗
DOE OSTI2023

The starch-deficient plastidic PHOSPHOGLUCOMUTASE mutant of the constitutive crassulacean acid metabolism (CAM) species Kalanchoë fedtschenkoi impacts diel regulation and timing of stomatal CO2 responsiveness

Abstract Background and Aims Crassulacean acid metabolism (CAM) is a specialized type of photosynthesis characterized by a diel pattern of stomatal opening at night and closure during the day, which increases water-use efficiency. Starch degradation is a key regulator of CAM, providing phosphoenolpyruvate as a substrate in the mesophyll for nocturnal assimilation of CO2. Growing recognition of a key role for starch degradation in C3 photosynthesis guard cells for mediating daytime stomatal opening presents the possibility that starch degradation might also impact CAM by regulating the provision of energy and osmolytes to increase guard cell turgor and drive stomatal opening at night. In this study, we tested the hypothesis that the timing of diel starch turnover in CAM guard cells has been reprogrammed during evolution to enable nocturnal stomatal opening and daytime closure. Methods Biochemical and genetic characterization of wild-type and starch-deficient RNAi lines of Kalanchoë fedtschenkoi with reduced activity of plastidic phosphoglucomutase (PGM) constituted a preliminary approach for the understanding of starch metabolism and its implications for stomatal regulation in CAM plants. Key Results Starch deficiency reduced nocturnal net CO2 uptake but had negligible impact on nocturnal stomatal opening. In contrast, daytime stomatal closure was reduced in magnitude and duration in the starch-deficient rPGM RNAi lines, and their stomata were unable to remain closed in response to elevated concentrations of atmospheric CO2 administered during the day. Curtailed daytime stomatal closure was linked to higher soluble sugar contents in the epidermis and mesophyll. Conclusions Nocturnal stomatal opening is not reliant upon starch degradation, but starch biosynthesis is an important sink for carbohydrates, ensuring daytime stomatal closure in this CAM species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Charge radii in covariant density functional theory: A global view

A systematic global investigation of differential charge radii has been performed within the CDFT framework for the first time. Theoretical results obtained with conventional covariant energy density functionals and the separable pairing interaction of Tian et al. [Phys. Lett. B 676, 44 (2009)] are compared with experimental differential charge radii in the regions of the nuclear chart in which available experimental data crosses the neutron shell closures at N=28,50,82, and 126. The analysis of absolute differential radii of different isotopic chains and their relative properties indicate clearly that such properties are reasonably well described in model calculations in the cases when the mean-field approximation is justified. However, while the observed clusterization of differential charge radii of different isotopic chains is well described above the N=50 and N=126 shell closures, it is more difficult to reproduce it above the N=28 and N=82 shell closures because of possible deficiencies in the underlying single-particle structure. The impact of the latter has been evaluated for spherical shapes and it was shown that the relative energies of the single-particle states and the patterns of their occupation with increasing neutron number have an appreciable impact on the evolution of the δ$\langle$r2$\rangle$N,N' values. These factors also limit the predictive power of model calculations in the regions of high densities of the single-particle states of different origin. It is shown that the kinks in the charge radii at neutron shell closures are due to the underlying single-particle structure and due to weakening or collapse of pairing at these closures. The regions of the nuclear chart in which the correlations beyond mean field are expected to have an impact on charge radii are indicated; the analysis shows that the assignment of a calculated excited prolate minimum to the experimental ground state allows us to understand the trends of the evolution of differential charge radii with neutron number in many cases of shape coexistence even at the mean-field level. It is usually assumed that pairing is a dominant contributor to odd-even staggering (OES) in charge radii. Our analysis paints a more complicated picture. It suggests a new mechanism in which the fragmentation of the single-particle content of the ground state in odd-mass nuclei due to particle-vibration coupling provides a significant contribution to OES in charge radii.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Machine Learning Assisted Safety Modeling and Analysis of Advanced Reactors

With the advances in computational power and numerical methods, analysts can now rely on first-principle simulations to predict ultra-fine details in a variety of applications. Advances in machine learning (ML) have produced algorithms that can now learn high-level abstractions via hierarchical models. This project aims to leverage advances in ML techniques and the available high-resolution simulation data to develop a novel modeling and simulation (M\&S) methodology for reactor safety analysis. While application-agnostic ML techniques are available, complex physics constraints need to be incorporated into ML techniques to build ML-based closures for computationally efficient predictive simulations. This project intends to develop a physics-guided data-driven multi-scale methodology for M\&S of advanced reactors. The project focuses on thermal fluid (T/F) phenomena, which play major roles in advanced reactor safety. Specifically, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor (SFR). The model has a coarse-mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data with Reynolds-averaged Navier-Stokes (RANS) turbulence model to ensure accuracy. Three different neural networks are developed and tested for loss-of-flow transients in the hot pool of SFR, i.e. the densely connected convolutional neural network (DCNN), long-short-term-memory network based on proper orthogonal decomposition (POD-LSTM), and the DCNN informed by LSTM (DCNN-LSTM). The performances of these three neural networks are evaluated based on baseline models. The DCNN-LSTM model has been chosen for further hyperparameter optimization. Furthermore, based on a simplified two-dimensional case, uncertainty quantification (UQ) of the developed ML-based closure are investigated with three methods, i.e. Monte Carlo dropout, deep ensemble, and Bayesian neural network. The developed ML-based turbulent viscosity closure relation based on deep ensemble is then integrated into the system analysis module SAM and serves as a term in the conservation equations. Such a SAM-ML based procedure guarantees that the obtained results are consistent with the physical constraints of the thermal-fluid system. The SAM-ML simulation on the same loss-of-flow transient showed comparable accuracy with the CFD simulation but with a much coarser mesh setup. Last but not least, the ML-based closure improvement with the support of higher-fidelity data from large eddy simulation (LES) is discussed. As a first step towards this direction, a baseline LES simulation is performed to obtain comparable data with RANS results. Based on the early results, future investigation on further improving the ML-based closure is discussed. We believe the developed approach that combines scientific machine learning with nuclear system analysis code can benefit the advanced reactor community as more accurate safety analyses will better characterize reactor safety margins and reduce licensing efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Composite Analysis for Low-Level Waste Disposal in the Hanford Site Central Plateau (FY 2020) (Rev. 2)

This document provides an updated Hanford Site composite analysis (CA). A Hanford Site CA was prepared and issued in 1998 (PNNL-11800) with an addendum provided in 2002 (PNNL-11800 Addendum 1). The CA was approved in 2002 (Frei, 2002) and has been maintained (DOE/RL-2000-29 and subsequent revisions) to support low-level waste disposal performance assessments (PAs) and disposal authorizations for facilities at the Hanford Site, including the following: Continued operation of the Environmental Restoration Disposal Facility (ERDF) and the 200 East and 200 West Low-Level Burial Grounds; Construction of the Integrated Disposal Facility (IDF); Forthcoming closure of tank residual waste systems such as Waste Management Area (WMA) C. The CA maintenance program resulted in a determination in 2015 (DOE/RL-2015-66) that the Hanford Site CA needed an update for the following reasons: While the initial Hanford Site CA has been maintained since 2001, the accumulation of basis changes reported in the annual summary reports over the succeeding 14 years merit evaluation in an updated analysis; The U.S. Department of Energy (DOE) Headquarters requested in a memorandum in 2015 (Gilbertson and Marcinowski, 2015) that “as soon as the relevant PAs are complete, the CA will be revised to account for all of the new information.” This updated Hanford Site CA accounts for the following new information: 1. Inclusion of a detailed Hanford Site baseline disposition that projects remedial activities through site closure. There have been significant changes through decision making in the Comprehensive Environmental Response, Compensation, and Liability Act of 1980 (CERCLA) process that were not available when the original CA was produced; 2. Inclusion of an updated inventory basis, new modeling capabilities, and new decisions reached in the associated record of decision (ROD) that was provided by issuance of DOE/EIS-0391 in fiscal year (FY) 2013. Development of a Hanford Site groundwater model from the baseline provided in a technical transfer of models for the Final Tank Closure and Waste Management Environmental Impact Statement (DOE/EIS-0391) commenced in FY 2014 and resulted in the plateau-to-river (P2R) groundwater flow model that is used in this updated Hanford Site CA; 3. Accounting for pump and treat systems, which were not evaluated in the initial CA. Pump and treat systems have had significant impact on groundwater flow system behavior, contaminant transport, and contaminant removal from Hanford Site groundwater. This process is accounted for in the use of the P2R model for this updated Hanford Site CA; 4. Inclusion of water-level data collected since the initial CA was completed. Data collected as the unconfined aquifer water levels continue to recede since the cessation of large liquid discharges in the late 1990s have led to marked improvement in understanding of the flow system for future conditions, particularly regarding northward flow potential in the critical Gable Gap area. The calibration in the latest version of the P2R model accounts for this information; 5. Inclusion of results from the revised ERDF PA completed in 2013 (WCH-5209) that account for updated inventory and expansion of the ERDF facility to about twice the size that was evaluated in the original CA; 6. Use of updated geoframeworks to provide the structural basis for numerical fate and transport models in the groundwater pathway of this updated CA. The geologic basis for groundwater models has continued to improve with additional data collection and interpretation with the creation and maintenance of the Hanford South Geoframework and the Central Plateau Vadose Zone Geoframework tools; 7. Incorporation of updated tank residual inventory estimates. Tank residual inventory estimates have improved with the incorporation of tank retrieval inventory data for those tanks that have completed retrieval. The CA inventory data package includes this updated information; 8. Incorporation of WMA C PA results. Two additional PAs for tank farm closure decisions are in preparation during the period required to prepare an updated Hanford Site CA: WMA C (FY 2016) and WMA A-AX (in preparation). The WMA C PA results are incorporated into the updated Hanford Site CA, and its grouted residuals model is used as the basis for a release model to account for the other tank farm systems modeling in this updated CA; 9. An update to the IDF PA was submitted in FY 2017 and has been reviewed and approved. The results of the IDF PA are incorporated into this updated CA; 10. Updated risk assessment scenarios. The risk assessment scenarios currently in use for Hanford Site CERCLA and Resource Conservation and Recovery Act of 1976 (RCRA) analyses differ from those evaluated in the initial Hanford Site CA. The representative person exposure scenario evaluated in this updated Hanford Site CA is consistent with recent PAs and CERCLA and RCRA analyses as the Hanford Site. This updated Hanford Site CA provides the following: A comparison of the updated Hanford Site CA all-pathway dose results with the performance measures during the compliance period, which is assumed to begin with site closure in calendar 2070 with the last scheduled disposal action and continue for 1,000 years postclosure (to calendar year 3069); A comparison of Hanford Site CA all-pathway dose results with the performance measures during the postcompliance period to address potential peaks beyond the compliance period. This is accomplished by evaluating dose in the period for 9,000 years following the compliance period (i.e., from calendar year 3070 to 12070).

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

Continued Development and Advanced Testing of DPC Filler Cements (on FY22 R&D and Demonstration Activities) (Progress Report)

Commercial generation of energy by nuclear power plants in the United States (U.S.) has produced thousands of metric tons of spent nuclear fuel (SNF), the disposal of which is the responsibility of the U.S. Department of Energy (DOE). Utilities typically utilize the practice of storing this SNF in dual-purpose canisters (DPCs). DPCs were designed, licensed, and loaded to meet Nuclear Regulatory Commission (NRC) requirements that preclude the possibility of a criticality event during SNF storage and transport, but were not designed or loaded to preclude the possibility of a criticality event during the regulated post-closure period following disposal, which could be up to 1,000,000 years (Price, 2019). There are several options being investigated that could facilitate the disposal of SNF stored in DPCs in a geologic repository (Hardin et al., 2015; SNL 2020b; SNL 2021b). These include: (1) repackage the SNF into canisters that are designed to prevent criticality during the regulated post-closure period following disposal, but with an increased disposal cost estimated at approximately $\$$20B in United States dollars (USD) (Freeze et al., 2019); (2) analysis of the probability and consequences of criticality from the direct disposal of DPCs during a 1,000,000-year post-closure period in several geologic disposal media (Price, 2019); and (3) filling the void space of a DPC with a material before its disposal that significantly limits the potential for criticality over the post-closure regulatory period. This report further investigates the third option, filling DPC already containing SNF with a material to limit the potential for criticality over the post-closure regulatory period.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Streamlined Approach for Environmental Restoration (SAFER) Plan for Corrective Action Unit 114: Area 25 EMAD Facility Nevada National Security Site, Nevada

This Streamlined Approach for Environmental Restoration (SAFER) Plan addresses the actions needed to achieve closure for Corrective Action Unit (CAU) 114, Area 25 EMAD Facility, identified in the Federal Facility Agreement and Consent Order (FFACO). CAU 114 comprises the following corrective action sites (CASs) located in Area 25 of the Nevada National Security Site: • 25-41-03, EMAD Facility (Building 3900) • 25-99-23, Manned Control Car (MCC) and Engine Installation Vehicle (EIV) • 25-33-05, Building 3901, Engine Transport System Maintenance Building (Train Shed) This plan provides the methodology for field activities needed to gather the necessary information for closing CAU 114. There is sufficient information and process knowledge from historical documentation and investigations of similar sites regarding the expected nature and extent of potential contaminants to recommend closure of CAU 114 using the SAFER process. Additional information will be obtained by conducting a field investigation before selecting the appropriate corrective actions for CAU 114. It is anticipated that the results of the field investigation and implementation of corrective actions will support a defensible recommendation that no further corrective action is necessary. The purpose of the corrective action investigation will be to document and verify the adequacy of existing information; to affirm the decision for either clean closure, closure in place, or no further action; and to provide sufficient data to implement the corrective action. The actual corrective action selected will be based on characterization activities implemented under this SAFER Plan. If specific conditions or findings fall outside the bounds of the conceptual site model, such as an unanticipated release, the Nevada Division of Environmental Protection (NDEP) will be consulted to determine the path forward before proceeding. Upon completion of SAFER activities, a closure report (CR) will be prepared and submitted to NDEP for review and approval. The schedule for completion of the CR will be established in consultation with NDEP.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

E-Area Low-Level Waste Facility Inadvertent Human Intruder Limits and Doses in Support of the PA2022

This report documents the inadvertent human intruder (IHI) analysis for the E-Area Low-Level Waste Facility (ELLWF) at the Savannah River Site (SRS), near Aiken, South Carolina. This analysis supports the revised ELLWF Performance Assessment (PA), complying with the Department of Energy standard for operation of low-level waste disposal facilities (USDOE, 2017). The ELLWF is an operating waste disposal facility and is scheduled to continue accepting waste to 2065. One task of the revised PA is to establish waste inventory limits for the various disposal units at ELLWF. This is done by modeling future contaminant release and transport through applicable pathways to human receptors, comparing predicted doses per disposed curie with applicable performance measures, to obtain inventory limits which will assure that doses to receptors do not exceed performance measures. This report documents results of modeling future doses to one class of receptor, the inadvertent human intruder. It is assumed that after site closure, public knowledge of the site is lost, and IHIs will engage in activities on the ELLWF that will disrupt the closure cap, causing dose to the IHI. Following USDOE (2017), six different stylized exposure scenarios are considered, simulating activities by an IHI which could result in a radiological dose. The six scenarios are: • Acute – Basement Construction: IHI constructs a basement and encounters waste during excavation which is inadvertently mixed with clean soil and diluted. • Acute – Well Drilling: IHI drills a water well through waste and is exposed to drill cuttings mixed with clean soil that are brought to the surface. • Acute – Discovery: IHI begins constructing a basement but stops when encountering the riprap in the final closure cap and is exposed to photon radiation from unexcavated material residing in the undisturbed waste zone. • Chronic – Agriculture: Resident IHI is exposed to waste that was excavated for basement construction and mixed with native soil in the intruder’s vegetable garden. • Chronic – Post-Drilling: Resident IHI is exposed to waste from drill cuttings mixed with native soil and scattered in the garden area. • Chronic – Residential: Resident IHI is exposed to external radiation while in home located above waste with shielding provided by the concrete basement floor and any soil or engineered material remaining between the basement and waste. Dose calculations are performed using the SRNL Dose Toolkit (Aleman, 2023), following the approach of Smith et al (2019). Calculations are performed separately for 27 of the 33 disposal units (DUs) at ELLWF and are radionuclide specific. The results of the IHI analysis include: • Dose Factors: mrem per disposed curie (acute) and mrem/yr per disposed curie (chronic) for each parent radionuclide, for each DU. • Inventory Limits: in curies, for each parent radionuclide, for each DU. • Estimated Dose to IHI: mrem (acute) and mrem/yr (chronic), for each DU, given its projected closure inventory without inventory biases applied. Most DU-specific IHI inventory limits are in the range of 10 3 to 10 7 curies per nuclide. The lowest inventory limits are associated with gamma-emitters such as Sn-126, Ra-226, Th-232, and Cm-248. Radionuclides with short half-lives such as Pu-241, and nuclides which are pure beta emitters or which decay by electron capture, such as Ni-59 and Ni-63, have the highest limits. For the 27 evaluated DUs, predicted IHI doses are shown in Table ES-1. The maximum acute dose is 1.18 mrem, at ST23, much less than the DOE performance measure of 500 mrem (USDOE, 2017). The highest chronic dose is 37.2 mrem/yr at ST02, below the DOE performance measure of 100 mrem/yr. Also shown are estimated inventory sums of fractions (SOFs) at closure in 2065, for groundwater (GW) and IHI pathways. For each DU, the inventory is constrained by the GW pathway. For most DUs, the IHI SOFs are approximately 1000 times lower than the GW SOF values, and the IHI pathway does not drive risk for any disposal unit.

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

Contaminant Migration Modeling to Support the In Situ Decommissioning of Hardened Nuclear Facilities - 20215

The decommissioning of hardened nuclear facilities provides a unique challenge in balancing current worker dose with exposure to future receptors and media. In situ decommissioning (ISD) is a cost-effective and safe option for the closure of such facilities but requires a detailed understanding of the threat that radioactive and hazardous constituents remaining in the facilities pose to groundwater and surrounding surface water. The Savannah River Site (SRS) extensively employs contaminant migration models within the decommissioning process to develop site-specific removal, grouting, and monitoring strategies in support of safe and effective final end states of nuclear facilities. This work discusses the use of contaminant migration modeling in the ISD process at SRS, including specific examples from the closure of P-Reactor and R-Reactor, Building 235-F, F-Canyon and F/H Laboratory Complex. Migration models developed during the closure of P-Reactor and R-Reactor were divided into four source areas: the reactor vessel, the process area, the disassembly basin, and the purification wing. Each source area was assigned a specific inventory and migration pathway and was then subject to varying hypothetical removal, capping, and grouting schemes. Modeling indicated that groundwater Maximum Concentration Limits (MCLs) might be exceeded if no action was taken for eleven and ten constituents of concern (COCs) at P-Reactor and R-Reactor, respectively. These exceedances could occur in as few as 200 to 500 years. Alternatively, the migration modeling demonstrated that the selected ISD actions reduced contaminant mobility which allowed for significant radioactive decay and resulted in fewer predicted exceedances of groundwater MCLs (five COCs for P-Reactor; eight COCs for R-Reactor). In response to the modeling results for P-Reactor and R-Reactor, effectiveness monitoring programs were developed to target contaminants, specific to each reactor, that may migrate to groundwater. Contaminant migration modeling also revealed that roof collapse was a large factor in the release of COCs to the environment, giving rise to roof improvements, and an inspection and vegetation control program to ensure roof stability over time. Contaminant migration modeling is also aiding in the closure planning for hardened facilities in F Area. Building 235-F housed the Actinide Billet Line, which produced Np-237 billets for irradiation in SRS reactors, and the Plutonium Fuel Form (PuFF) facility that produced Pu-238 heat sources for the space program. As a result of these missions, areas within Building 235-F contain considerable residual amounts of both Pu-238 and Np-237. Contaminant migration modeling of Building 235-F was originally performed in 2012 to identify the feasibility of ISD and the amount of radioactive material removal required to prevent the exceedance of groundwater MCLs. The original model indicated that a 60% reduction in the PuFF facility Pu-238 inventory could keep groundwater concentrations below standards, while Np-237 did not pose a threat to groundwater. However, updates to the model with an emphasis on source impact pathways revealed that, due to the orientation of the source areas relative to groundwater flow, no amount of reasonable removal of Pu-238 would keep groundwater concentrations below MCLs and that Np-237 could be a large contributor to localized MCL exceedances. With this insight, the refined 2019 model is being used to assist in the development of grouting plans specific to each facility source area, where bentonite may be utilized to slow the migration of Pu-238 and its daughter products from the PuFF facility and a reducing grout may decrease Np-237 transport by ensuring the nuclide remains in the less mobile +IV oxidation state. The beginning phases of contaminant migration modeling are underway for F-Canyon and associated facilities using lessons learned from the Reactors and Building 235-F. A contaminant migration pathway, similar to the pathway used for the reactor vessels, is being developed for the hot and warm canyons. The F-Canyon inventories are being spatially refined to identify specific and localized source zones, comparable to Building 235-F, that may impact groundwater. Contaminant migration modeling provides vital input within the ISD process, from facility investigation to post closure effectiveness monitoring, making it a valuable tool in the development of safe and effective end states for hardened nuclear facilities at SRS. (authors)

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

Data-driven variational multiscale reduced order models

We propose a new data-driven reduced order model (ROM) framework that centers around the hierarchical structure of the variational multiscale (VMS) methodology and utilizes data to increase the ROM accuracy at a modest computational cost. The VMS methodology is a natural fit for the hierarchical structure of the ROM basis: In the first step, we use the ROM projection to separate the scales into three categories: (i) resolved large scales, (ii) resolved small scales, and (iii) unresolved scales. In the second step, we explicitly identify the VMS–ROM closure terms, i.e., the terms representing the interactions among the three types of scales. In the third step, we use available data to model the VMS–ROM closure terms. Thus, instead of phenomenological models used in VMS for standard numerical discretizations (e.g., eddy viscosity models), we utilize available data to construct new structural VMS–ROM closure models. Specifically, we build ROM operators (vectors, matrices, and tensors) that are closest to the true ROM closure terms evaluated with the available data. We test the new data-driven VMS–ROM in the numerical simulation of four test cases: (i) the 1D Burgers equation with viscosity coefficient v = 10 -3 ; (ii) a 2D flow past a circular cylinder at Reynolds numbers, Re = 100, Re = 500, and Re= 1000; (iii) the quasi-geostrophic equations at Reynolds number Re = 450 and Rossby number Ro = 0.0036; and (iv) a 2D flow over a backward facing step at Reynolds number Re = 1000. The numerical results show that the data-driven VMS–ROM is significantly more accurate than standard ROMs.

42 ENGINEERING↗
DOE OSTI2023

On the design of stable, consistent, and conservative high-order methods for multi-material hydrodynamics

Obtaining stable and high-order numerical solutions for multi-material hydrodynamics is an open challenge. Although slope limiters are widely used to maintain monotonicity near discontinuities, typical limiting procedures violate closure laws at the discrete level when applied to multi-material hydrodynamics equations. Due to this, the high-order expansions of quantities related by the closure laws are no longer consistent. The commonly observed symptom of this consistency-violation is that the numerical method fails to maintain constant pressure and velocity across material interfaces. This leads to sub-optimal convergence rates for smooth multi-material problems as well. Specialized limiting procedures that satisfy consistency while maintaining conservation need to be developed for such equations. A novel procedure that re-instates consistency into slope-limited high-order discretizations applied to the multi-material hydrodynamics equations is presented here. Using simple examples, it is demonstrated that the presented method satisfies closure laws at the discrete level, while maintaining conservative properties of the high-order method. Furthermore, this procedure involves a projection step which relies on the compact basis of the underlying spatial discretization, i.e. for discontinuous schemes (viz. DG and FV) the projection is local, and does not involve global matrix solves. Comparisons with conventional approaches emphasizes the necessity of the consistent closure-law preserving limiting approach, in order to maintain design order of accuracy for smooth multi-material problems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Crystal mechanics-based thermo-elastic constitutive modeling of orthorhombic uranium using generalized spherical harmonics and first-order bounding theories

In earlier works, a mathematical procedure for invertible microstructure-property linkages was developed using computationally efficient spectral methods for polycrystalline cubic and hexagonal metals. This paper formulates such invertible microstructure–property linkages for orthorhombic polycrystalline metals relying on the generalized spherical harmonics (GSH) spectral basis. The procedure is used to compute property closures of orthorhombic polycrystals. The closures represent the complete set of theoretically possible combinations of effective properties for a selected material. The procedure relies on the first-order bounding theories and considers orientation distribution functions (ODFs) as the main microstructural descriptor influencing homogenized properties. Numerous examples of these closures involving second-rank thermal expansion and fourth-rank elastic stiffness tensorial properties over a broad range of temperatures are presented for α-uranium (α-U). In doing so, certain key properties of these closures are exploited to facilitate their computation with drastically reduced computational effort. Along with the recently developed GSH-based interpolation procedure for ODFs from coarsely spaced experimental measurement grids to finely spaced finite element mesh resolution grids presented in Barrett et al., the developed computationally efficient ODF-effective property linkages are used to establish a crystal mechanics-based simulation framework coupled with the finite element method (FEM). The ODF dependent thermal expansion and elastic stiffness tensors are efficiently calculated at every integration point and used by the FEM to predict the overall distortion of a hemispherical part made of α-U during heating. In conclusion, it is shown that the developed framework can be used to simulate microstructurally heterogeneous components under thermo-mechanical loadings in a computationally efficient manner.

36 MATERIALS SCIENCE↗
DOE OSTI2025

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Sensitivity Analysis of Wind and Turbulence Predictions With Mesoscale‐Coupled Large Eddy Simulations Using Ensemble Machine Learning

Abstract Coupling between mesoscale models and large‐eddy simulation (LES) models is increasingly used to more realistically represent the wide range of scales of atmospheric motions affecting boundary layer winds and turbulence that need to be simulated accurately for applications such as wind energy. However, such mesoscale‐to‐microscale coupled modeling frameworks are potentially affected by a large number of uncertain closure parameters. Here, we investigate the sensitivity associated with six closure parameters related to a 1.5‐order subgrid‐scale turbulence closure for an ensemble of mesoscale‐coupled LES. The simulations are performed using the Weather Research and Forecasting model nested from horizontal resolutions of greater than a kilometer down to tens of meters. Closure parameters are varied to generate perturbed parameter ensembles for two case studies of highly sheared, convective boundary layers observed in the Columbia Basin of Oregon and Washington during the Second Wind Forecast Improvement Project. Machine learning algorithms are used to explore the sensitivity of LES predictions, considering the effects of the perturbed physical parameters alongside categorical factors such as the case study identity, measurement location, and LES resolution. For the conditions we examine, a single parameter, the eddy viscosity coefficient, is the dominant source of parametric sensitivity and its importance is comparable to the categorical factors for several of the simulation response variables we examine.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reduction of Pacific Double‐ITCZ Bias by Convection Parameterization in NCAR CESM2.2

Abstract The impact of convective closure on the double‐ITCZ bias in the NCAR CESM2.2 is investigated in this study. The standard CESM2.2 simulates a remarkable double‐ITCZ bias in the central and eastern Pacific, especially in boreal winter and spring. Modifications to the closure in convection parameterization scheme greatly reduce the double‐ITCZ bias in all seasons, demonstrating that convection parameterization can substantially influence the double‐ITCZ bias in CESM2.2. Further analyses suggest that convection parameterization can modulate the tropical atmosphere‐ocean feedback processes, through which it influences the SST in the southern ITCZ region and hence the double‐ITCZ bias. The changes in the upper ocean temperature advection induced by modified convective closure plays important roles in reducing the warm SST bias and double‐ITCZ precipitation bias in the southern ITCZ region. The modified convective closure improves the low‐level cloud and shortwave cloud radiative forcing in the southeastern Pacific. However, surface heat flux plays only a limited role in reducing warm SST bias and double ITCZ bias because the impacts of shortwave radiation changes are largely canceled by changes in longwave radiation and latent heat flux.

Song, Xiaoliang [Scripps Institution of Oceanograp↗
DOE OSTI2023

Memory-based parameterization with differentiable solver: Application to Lorenz ’96

Physical parameterizations (or closures) are used as representations of unresolved subgrid processes within weather and global climate models or coarse-scale turbulent models, whose resolutions are too coarse to resolve small-scale processes. These parameterizations are typically grounded on physically based, yet empirical, representations of the underlying small-scale processes. Machine learning-based parameterizations have recently been proposed as an alternative solution and have shown great promise to reduce uncertainties associated with the parameterization of small-scale processes. Yet, those approaches still show some important mismatches that are often attributed to the stochasticity of the considered process. This stochasticity can be due to coarse temporal resolution, unresolved variables, or simply to the inherent chaotic nature of the process. To address these issues, we propose a new type of parameterization (closure), which is built using memory-based neural networks, to account for the non-instantaneous response of the closure and to enhance its stability and prediction accuracy. We apply the proposed memory-based parameterization, with differentiable solver, to the Lorenz ’96 model in the presence of a coarse temporal resolution and show its capacity to predict skillful forecasts over a long time horizon of the resolved variables compared to instantaneous parameterizations. This approach paves the way for the use of memory-based parameterizations for closure problems.

Mathematics↗