Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure modeling”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Summary of Developing a Hanford Grout Modeling Framework Initial Workshop

The “Development of the Hanford Grout Modeling Framework” project as part of the Network of National Laboratories for Environmental Management and Stewardship (NNLEMS) funding program is underway. An initial project team workshop was held in Richland, WA on June 4 and 5, 2024 with participants from PNNL, SRNL, SNL, WRPS and CRESP (VU and UCF). The workshop covered an overview of the current state of the art in grout modeling tools provided by the team. Based on a focus on liquid waste immobilization grout and encapsulation grout for solid wastes and tank closure grout at Hanford, the team identified preliminary conceptual models for these three applications which will be used to form the basis of the Hanford Grout Modeling Framework. The key processes that can drive mechanisms impacting the grout performance were identified for each test case. The programmatic plan was also identified with five core Themes being selected that will form the focus of the project. Within each theme, individual tasks for Year 1 were planned in support of the development of an updated Hanford Grout Modeling Framework. This report presents a summary of the information presented during the workshop.

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

Theory for transport in magnetized plasmas

We present a two fluid model for anomalous thermal, particle and momentum transport which has been derived from kinetic theory using a small parameter ε~10-2. This small parameter can come from either the ratio of wave frequency to cyclotron frequency or the ratio of potential to thermal energies. This fluid model is strongly nonlinear and has an exact closure within the accuracy ε. It can be used to simulate the full radius of a tokamak including internal and edge transport barriers at the same time without introducing any local assumptions for physics model or grid size. It involves a full transport matrix, thus giving possibility for both density and heat pinches. The model can be derived from the Vlasov equation, adding also close collisions and includes particle pinches of the type found in the magnetosphere and in the Massachusetts Institute of Technology elevated dipole machine. Pinches are also in agreement for the heat pinch experiments on Doublet III D including simultaneous particle and heat pinches. For tokamaks it gives the right scaling of confinement time versus heating power. Finally, we have also recovered the Turbulent Equipartition result that the density scales as inverse safety factor q in typical simulations. Thus this model is both fundamental and of high practical use.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Digital Infrastructure Industry Engagement

The commercial nuclear sector faces unprecedented financial challenges driven by low natural gas prices and subsidized renewables in a market that does not reward carbon-free baseload capacity. These circumstances, along with increasingly antiquated labor-centric operating models and analog technology, have forced the early closure of multiple nuclear facilities and placed a much larger population of nuclear stations at risk. Nuclear plant economic survival in current and forecasted market conditions requires an efficient and technology-centric operating model that harvests the native efficiencies of advanced technology. This is analogous to transformations that have occurred in other industries.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Understanding the Spatiotemporal Variability of Tropical Orographic Rainfall Using Convective Plume Buoyancy

Abstract Mechanical forcing by orography affects precipitating convection across many tropical regions, but controls on the intensity and horizontal extent of the orographic precipitation peak and rain shadow remain poorly understood. A recent theory explains this control of precipitation as arising from modulation of lower-tropospheric temperature and moisture by orographic mechanical forcing, setting the distribution of convective rainfall by controlling parcel buoyancy. Using satellite and reanalysis data, we evaluate this theory by investigating spatiotemporal precipitation variations in six mountainous tropical regions spanning South and Southeast Asia, and the Maritime Continent. We show that a strong relationship holds in these regions between daily precipitation and a measure of convective plume buoyancy. This measure depends on boundary layer thermodynamic properties and lower-free-tropospheric moisture and temperature. Consistent with the theory, temporal variations in lower-free-tropospheric temperature are primarily modulated by orographic mechanical lifting through changes in cross-slope wind speed. However, winds directed along background horizontal moisture gradients also influence lower-tropospheric moisture variations in some regions. The buoyancy measure is also shown to explain many aspects of the spatial patterns of precipitation. Finally, we present a linear model with two horizontal dimensions that combines mountain wave dynamics with a linearized closure exploiting the relationship between precipitation and plume buoyancy. In some regions, this model skillfully captures the spatial structure and intensity of rainfall; it underestimates rainfall in regions where time-mean ascent in large-scale convergence zones shapes lower-tropospheric humidity. Overall, these results provide new understanding of fundamental processes controlling subseasonal and spatial variations in tropical orographic precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Business Case Analysis for Digital Safety-Related Instrumentation & Control System Modernizations

The commercial nuclear sector faces unprecedented financial challenges driven by low natural gas prices and subsidized renewables in a marketplace that does not reward carbon-free baseload capacity. These circumstances, along with increasingly antiquated labor-centric operating models and analog technology, have forced the premature closure of multiple nuclear facilities and placed a much larger population of nuclear power stations at risk. Nuclear plant economic survival in current and forecasted market conditions requires an efficient and technology-centric operating model that harvests the native efficiencies of advanced technology. This is analogous to transformations in nearly every other industry. In light of previous industry experience in modernizing safety Instrumentation and Controls (I&C) systems, nuclear utilities are reluctant to pursue these upgrades due to uncertainty in licensing and cost.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Synchrotron self-Compton in a radiative-adiabatic fireball scenario: modelling the multiwavelength observations in some Fermi /LAT bursts

Energetic GeV photons expected from the closest and the most energetic Gamma-ray bursts (GRBs) provide a unique opportunity to study the very-high-energy emission as well as the possible correlations with lower energy bands in realistic GRB afterglow models. In the standard GRB afterglow model, the relativistic homogeneous shock is usually considered to be fully adiabatic, however, it could be partially radiative. Based on the external forward-shock scenario in both stellar wind and constant-density medium, we present a radiative-adiabatic analytical model of the synchrotron self-Compton (SSC) and synchrotron processes considering an electron energy distribution with a power-law index of $1\lt p\lt 2$ and $2\le p$. We show that the SSC scenario plays a relevant role in the radiative parameter $\epsilon$, leading to a prolonged evolution during the slow cooling regime. In a particular case, we derive the Fermi/LAT light curves together with the photons with energies $\ge 100$ MeV in a sample of nine bursts from the second Fermi/LAT GRB catalogue that exhibited temporal and spectral indices with $\gtrsim 1.5$ and $\approx 2$, respectively. These events can hardly be described with closure relations of the standard synchrotron afterglow model, and also exhibit energetic photons above the synchrotron limit. We have modelled the multiwavelength observations of our sample to constrain the microphysical parameters, the circumburst density, the bulk Lorentz factor, and the mechanism responsible for explaining the energetic GeV photons.

acceleration of particles↗
DOE OSTI2021

Factors Governing Cloud Growth and Entrainment Rates in Shallow Cumulus and Cumulus Congestus During GoAmazon2014/5

Shallow cumulus and cumulus congestus clouds play an important role in the large-scale tropical circulation by mixing heat and moisture vertically and preconditioning the environment for deeper convection. Different representations of these shallow clouds account for much of the spread in General Circulation Model (GCMs) climate sensitivity, potentially because of how entrainment is represented in GCM parameterizations. This study uses observations from the Department of Energy's Atmospheric Radiation Measurement (ARM) mobile facility deployed at Manacapuru, Brazil, during the Green Ocean Amazon (GoAmazon2014/5) Campaign. Environmental thermodynamic profiles and observations of cloud top height (CTH) are used to constrain an entraining plume model to estimate bulk entrainment rates (ERs). Estimates of CTH are obtained from a combination of vertically pointing W-band ARM cloud radar and 1,290 MHz Radar Wind Profiler observations. A combination of radiosonde, microwave radiometer profiler, and microwave radiometer observations provides new best estimates of the environmental thermodynamic state. We quantify uncertainty in ERs considering uncertainties in estimated CTH, environmental thermodynamic properties, and assumed initial parcel characteristics. We find ERs ranging from 0.16 to 2.8 km -1 with an average of 0.58 ± 0.10 km -1 over a selected population of 469 shallow cumulus and cumulus congestus clouds. Using the retrieved estimates of ER, we evaluate several entrainment closures that are currently used in atmospheric models or have been proposed based on theory or large eddy simulation. Finally, entrainment rates in cumulus clouds are weakly correlated with low-level buoyancy, cloud depth, and cloud size.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Precision mass measurement of lightweight self-conjugate nucleus 80 Zr

Protons and neutrons in the atomic nucleus move in shells analogous to the electronic shell structures of atoms. The nuclear shell structure varies due to changes of the nuclear mean field with the number of neutrons N and protons Z, and these variations can be probed by measuring mass differences between nuclei. The N=Z=40 self-conjugate nucleus 80 Zr is of particular interest as its proton and neutron shell structures are expected to be very similar, and its ground state is highly deformed. Here, we provide evidence for the existence of a deformed double shell closure in 80Zr through high precision Penning trap mass measurements of 80-83 Zr. Our mass values show that 80Zr is significantly lighter, and thus more strongly bound than predicted. This can be attributed to the deformed shell closure at N=Z=40 and the large Wigner energy. A statistical Bayesian model mixing analysis employing several global nuclear mass models demonstrates difficulties with reproducing the observed mass anomaly using current theory.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Uncertainties in estimating global potential yields and their impacts for long-term modeling

Estimating realistic potential yields by crop type and region is challenging; such yields depend on both biophysical characteristics (e.g., soil characteristics, climate, etc.), and the crop management practices available in any site or region (e.g., mechanization, irrigation, crop cultivars). A broad body of literature has assessed potential yields for selected crops and regions, using several strategies. In this study we first analyze future potential yields of major crop types globally by two different estimation methods, one of which is based on historical observed yields (“Empirical”), while the other is based on biophysical conditions (“Simulated”). Potential yields by major crop and region are quite different between the two methods; in particular, Simulated potential yields are typically 200% higher than Empirical potential yields in tropical regions for major crops. Applying both of these potential yields in yield gap closure scenarios in a global agro-economic model, GCAM, the two estimates of future potential yields lead to very different outcomes for the agricultural sector globally. In the Simulated potential yield closure scenario, Africa, Asia, and South America see comparatively favorable outcomes for agricultural sustainability over time: low land use change emissions, low crop prices, and high levels of self-sufficiency. In contrast, the Empirical potential yield scenario is characterized by a heavy reliance on production and exports in temperate regions that currently practice industrial agriculture. At the global level, this scenario has comparatively high crop commodity prices, and more land allocated to crop production (and associated land use change emissions) than either the baseline or Simulated potential yield scenarios. This study highlights the importance of the choice of methods of estimating potential yields for agro-economic modeling.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Ion-Exchange Modeling of Crystalline Silicotitanate for Cesium Removal - 20283

The Tank Closure Cesium Removal (TCCR) system is a Savannah River Site (SRS) demonstration 'at-tank' process designed to remove {sup 137}Cs from the high-level aqueous tank waste so that the decontaminated solution can be disposed as low-level waste. Cesium is removed by ion exchange (IX) columns using engineered IONSIV{sup TM} R9120-B form of the Crystalline Silicotitanate (CST) media. The TCCR system is deployed at Tanks 10 and Tank 11 in the SRS H Tank Farm. Water is added to the salt-cake in Tank 10 H to dissolve it. The dissolved salt solution waste is pumped out of Tank 10H (feed tank), through filters and IX columns. The decontaminated salt solution is transferred to Tank 11 (receipt tank), and on to Tank 50H for final disposal in the Saltstone Production Facility. The current TCCR can accommodate lead-lag (two-column) or lead-lag-guard (three-column) configurations to optimize media utilization and achieve the target decontamination. To assist the TCCR operations, a parametric study was conducted to evaluate the impact of different parameters (e.g., column configurations (single column, two columns or three columns in series), waste characteristics, operating temperature, process flow rate, CST average particle size) on the IX column performance including CST bed utilization. The initial results indicate that the IX column performance is improved at slower process flow rate, at lower operating temperature, and with smaller CST average particle size. Multi-column configurations are recommended, because the single-column configuration does not utilize CST bed effectively. This paper demonstrates the versatility of the ion exchange modeling to evaluate the effects of CST characteristics and operational parameters on IX column performances. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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

Aircraft Observations of Turbulence in Cloudy and Cloud-Free Boundary Layers Over the Western North Atlantic Ocean From ACTIVATE and Implications for the Earth System Model Evaluation and Development

This study examines boundary layer turbulence derived from high temporal resolution meteorological measurements from 40 research flights over the western North Atlantic Ocean during the 2020 deployments of ACTIVATE. Frequency distributions of various turbulent quantities reveal stronger turbulence during the winter deployment than in summer and for cloud‐topped than in cloud‐free boundary layers during the summer deployment. Maximum turbulence kinetic energy (TKE) is most often within cloud from observations in winter and summer, whereas it is mostly below cloud in both seasons by a global model turbulence parameterization. Further, bivariate frequency distributions are consistent with the bivariate Gaussian probability distribution functions assumed for the closure of higher‐order turbulence/shallow convection parameterizations used by some global models. Turbulence simulated by the Community Atmosphere Model version 6 and the Energy Exascale Earth System Model Atmosphere Model version 2 using such parameterizations is not as strong as observed, with more TKE going into vertical wind perturbations rather than into zonal wind perturbations as observed, suggesting that the treatment of turbulence in Earth system models still needs to be further improved.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Better practices for inferring ecosystem water use strategy from eddy covariance data

Eddy covariance data are critical for inferring ecosystem water use strategies. Yet, such inferences are sensitive to a range of assumptions applied across studies, hindering our understanding of water use strategies within and across eddy covariance sites. A recent analysis across 151 FLUXNET2015 and AmeriFlux-FLUXNET datasets found that poor model performance was the key driver of non-robust inferences of ecosystem water use strategies. Here, we leverage this previous analysis to (i) identify the specific assumptions that improve inference model performance across most sites, (ii) explain the mechanisms behind the performance improvements, and (iii) check whether better performance improves water use inference. We find that the common practice of fitting a model to canopy conductance (G c ) derived from the evapotranspiration (ET) observations, rather than to observed ET itself, artificially amplifies data errors and degrades the model performance. Next, accounting for vegetation dynamics by applying a growing season filter or incorporating satellite LAI data improves performance, but the former practice may remove soil water stress periods. Lastly, using the leaf-to-air vapor pressure deficit (VPD l ) derived from ET observations as a model input may artificially inflate performance. Based on these results, we recommend selecting observed ET (rather than derived G c ) as the response variable, carefully accounting for vegetation dynamics, and avoiding derived VPD l as a model input; these best practices improve model performance by c. 20% and robustness by c. 80% across all eddy covariance sites. Nevertheless, the performance improvements do not always correspond to more robust inference of water use strategies, as model parameter selection and surface energy budget closure corrections still strongly influence the ecosystem water use parameter estimation in a site-specific manner.

AmeriFlux↗
DOE OSTI2021

The role of micro-inertia on the shock structure in porous metals

The behavior of porous materials under shock loading is a multi-scale problem bridging orders of magnitude across the macroscale geometry and the microscale pores. Under static loading, this problem is well understood, relating mechanisms of pore closure and crushing to the equivalent macroscale models. The dynamic response of porous solids under shock loading is related to the effects of viscoplasticity and micro-acceleration fields around the void boundaries. The significance of the micro-inertia effects in modeling the dynamic behavior of porous materials remains an open question. In this work, an experimental investigation on closed-cell porous aluminum with small porosity provides the evidence for the first time of micro-inertia’s fundamental role in describing the shock structure in these materials. Materials with different levels of porosity were manufactured using a modified process of additive manufacturing to achieve a mean pore size below 50$\mu$m. Plate impact experiments on porous aluminum samples were conducted at pressures in the range of 2 to 11 GPa. The structure of the steady shock was characterized as a function of porosity and shown to validate behavior revealed by an analytical approach (Czarnota et al. [J. Mech. Phys. Solids 107 (2017)]), highlighting the fundamental role of micro-inertia effects in such cases.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Checking nonflow assumptions and results via PHENIX published correlations in p+p, p+Au, d+Au, and 3 He+Au at √ s NN=200 GeV

Recently the PHENIX Collaboration has made available two-particle correlation Fourier coefficients for multiple detector combinations in minimum bias p+p and 0–5% central p+Au, d+Au, and 3 He+Au collisions at √ s NN = 200 GeV [Phys. Rev. C 105, 024901 (2022)]. Using these coefficients for three sets of two-particle correlations, azimuthal anisotropy coefficients v 2 and v 3 are extracted for midrapidity charged hadrons as a function of transverse momentum. In this paper, we use the available coefficients to explore various nonflow hypotheses as well as to compare the results with theoretical model calculations. The nonflow methods fail basic closure tests with ampt and pythia/angantyr, particularly when including correlations with particles in the low multiplicity light-projectile going direction. In data, the nonflow adjusted v 2 results are modestly lower in p+Au and the adjusted v 3 results are more significantly higher in p+Au and d+Au. However, the resulting higher values for the ratio v 3 /v 2 in p+Au at RHIC compared to p+Pb at the LHC is additional evidence for a significant overcorrection. Incorporating these additional checks, the conclusion that these flow coefficients are dominated by initial geometry coupled with final-state interactions (e.g., hydrodynamic expansion of quark-gluon plasma) remains true, and explanations based solely on initial-state glasma are ruled out. The detailed balance between intrinsic and fluctuation-driven geometry and the exact role of weakly versus strongly coupled prehydrodynamic evolution remains an open question for triangular flow, requiring further theoretical and experimental investigation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

GDSA - Modeling and Integration – ORNL: PFLOTRAN modeling

This report summarizes work done at Oak Ridge National Laboratory (ORNL) related to post-closure criticality consequence analysis of dual-purpose canister (DPC) direct disposal in unsaturated alluvium. PFLOTRAN software was used to model the behavior of a breached DPC with low water infiltration rates and unsaturated conditions in alluvial deposits such that a low-power criticality event occurs with a focus on thermal hydrological constraints on criticality timing and power output. This summary describes initial exploration of those thermal hydrological processes and of the long-term average power output from a criticality event that could be sustained without driving water out of the DPC package and thereby terminating the criticality event. The objectives are (1) to identify, for different groundwater infiltration rates, the time at which the DPC fills with water sufficiently to potentially initiate criticality; (2) to bound the power output that could be produced by a criticality event that would thermally drive water out of the package at various percolation rates; and (3) provide an estimate of how the size of a package breach affects the water infiltration and therefore the possibility timing of criticality.

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

Post-closure Nuclear Criticality Safety Evaluations for Disposition of Criticality Control Overpacks at the Waste Isolation Pilot Plant

The Waste Isolation Pilot Plant (WIPP) is a geological repository in southern New Mexico that provides for disposal of transuranic (TRU) wastes from atomic energy defense activities. The Sandia National Laboratories (Sandia) Report, Consideration of Nuclear Criticality When Disposing of Transuranic Waste at the Waste Isolation Pilot Plant, addresses nuclear criticality safety based on the projected inventory characteristics for the initial compliance certification application of WIPP in 1996. As the inventory, waste forms, and disposal package designs change, revised or new analyses are necessary to demonstrate acceptability for these configurations within the WIPP safety basis and compliance with 10,000-year post-closure standards of the US Environmental Protection Agency (EPA). Saylor and Scaglione evaluated criticality control overpacks (CCOs) in 2017 based on conservative assumptions for post-closure repository structural conditions with resulting effects on containers and container spacing, The Saylor and Scaglione evaluation of CCOs addressed a single waste configuration that represents the Surplus Plutonium Disposition Program’s dilute and dispose waste form and composition. This initial CCO study demonstrated that 50 grams of boron carbide (B 4 C) per CCO is sufficient to ensure post-closure criticality safety based on a well-mixed waste composition, and Oak Ridge National Laboratory (ORNL) subsequently determined that this amount of B 4 C does not require constraints on moisture or plastic present as moderator. The Saylor and Scaglione analysis conservatively assumes repository room closure that eliminates all space between fissile gram equivalent (FGE) 239 Pu masses. The close-packed array was selected based on limited availability of repository salt creep modeling results at that time. In 2019, Brickner provided additional evaluations for pipe overpack containers (POCs), building on the conservative basis provided by Saylor and Scaglione. Brickner’s 2019 analysis made use of new geomechanical data for post-closure spacing that rely on advances in repository modeling as documented in the work by Reedlunn and Bean. This current CCO evaluation for generic waste materials expands on earlier work performed at ORNL and includes evaluation of CCOs across a much broader range of possible waste compositions and geometries. This evaluation is intended to provide input for the required feature, event and process (FEP) screening to determine if post-closure criticality must be included as an event in the 10,000-year regulatory evaluation. As such, the approach to modeling post-closure criticality presented in this report has been coordinated with the Sandia team responsible for FEP screening. The resulting analysis supports disposition of fissile materials in the CCO containing up to 380 FGE 239 Pu and expands conditions acceptable for disposal of fissile material in CCOs. This evaluation builds on the methodology of Saylor and Scaglione and Brickner, using the most recently available geomechanical data for CCO spacing under salt creep compaction scenarios provided by Reedlunn and Bean. The broad range of fissile material configurations analyzed in this report are intended to account for configurations that may occur during the post-closure disposal time period, and it also includes waste configurations that are not physically possible to support analysis of conditions that influence neutron fluence.

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

“Scrutinizing entrainment and mass flux closures in shallow cumulus parameterizations using cloud-radar observations and large-eddy simulation” (Final report for DOE grant)

Warm, shallow clouds are responsible for a substantial portion of Earth System Model (ESM) uncertainty and divergence among model climate projections. A substantial amount of this uncertainty is tied to the representation of convection in ESMs, in particular assumptions underlying the convective parameterizations. This study focuses on entrainment and entrainment closures in ESM parameterizations. The completed project was a joint observational and modeling study to observationally derive estimates of entrainment and to use high-resolution modeling to scrutinize entrainment and mass flux closures in shallow convection parameterizations.

58 GEOSCIENCES↗