Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLOSURES”

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 487 records · Page 27

Study of the Alfven eigenmodes stability in CFQS plasma using a Landau closure model

The aim of this study is to analyze the stability of the Alfven eigenmodes (AE) in the Chinese First Quasi-axisymmetric Stellarator (CFQS). The AE stability is calculated using the code FAR3d that solves the reduced MHD equations to describe the linear evolution of the poloidal flux and the toroidal component of the vorticity in a full 3D system, coupled with equations of density and parallel velocity moment for the energetic particles (EP) species including the effect of the helical couplings and acoustic modes. The Landau damping and resonant destabilization effects are added in the model by a given closure relation. The simulation results indicate the destabilization of n = 1 to 4 AEs by EP during the slowing down process, particularly n = 1 and n = 2 toroidal AEs (TAE), n = 3 elliptical AE (EAE) and n = 4 non circular AE (NAE). If the resonance is caused by EPs with an energy above 17 keV (weakly thermalized EP), n = 2 EAEs and n = 3 NAEs are unstable. On the other hand, EPs with an energy below 17 keV (late thermalization stage) lead to the destabilization of n = 3 and n = 4 TAEs. The simulations for an off-axis NBI injection indicate the further destabilization of n = 2 to 4 AEs although the growth rate of the n = 1 AEs slightly decreases, so no clear optimization trend with respect to the NBI deposition region is identified. In addition, n = 2, 4 helical AE (HAE) are unstable above an EP β threshold. Also, if the thermal β of the simulation increases (higher thermal plasma density) the AE stability of the plasma improves. Overall, the simulations including the effect of the finite Larmor radius and electron-ion Landau damping show the stabilization of the n = 1 to 4 EAE/NAEs as well as a decrease of the growth rate and frequency of the n = 1 to 4 BAE/TAEs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear dynamics and transport driven by energetic particle instabilities using a gyro-Landau closure model *

Abstract Energetic particle (EP) destabilized Alfvén eigenmode (AE) instabilities are simulated for a DIII-D experimental case with a pulsed neutral beam using a gyro-Landau moments model which introduces EP phase-mixing effects through closure relations. This provides a computationally efficient reduced model which is applied here in the nonlinear regime over timescales that would be difficult to address with more complete models. The long timescale nonlinear evolution and related collective transport losses are examined including the effects of zonal flow/current generation, nonlinear energy cascades, and EP profile flattening. The model predicts frequencies and mode structures that are consistent with experimental observations. These calculations address issues that have not been considered in previous modelling: the EP critical gradient profile evolution in the presence of zonal flows/currents, and the dynamical nature of the saturated state. A strong level of intermittency is present in the predicted instability-driven transport; this is connected to the zonal flow growth and decay cycles and nonlinear energy transfers. Simulation of intermittent AE-enhanced EP transport will be an important issue for the protection of plasma facing components in the next generation of fusion devices.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Machine-learned closure of URANS for stably stratified turbulence: connecting physical timescales & data hyperparameters of deep time-series models

Stably stratified turbulence (SST), a model that is representative of the turbulence found in the oceans and atmosphere, is strongly affected by fine balances between forces and becomes more anisotropic in time for decaying scenarios. Moreover, there is a limited understanding of the physical phenomena described by some of the terms in the Unsteady Reynolds-Averaged Navier–Stokes (URANS) equations—used to numerically simulate approximate solutions for such turbulent flows. Rather than attempting to model each term in URANS separately, it is attractive to explore the capability of machine learning (ML) to model groups of terms, i.e. to directly model the force balances. We develop deep time-series ML for closure modeling of the URANS equations applied to SST. We consider decaying SST which are homogeneous and stably stratified by a uniform density gradient, enabling dimensionality reduction. We consider two time-series ML models: long short-term memory and neural ordinary differential equation. Both models perform accurately and are numerically stable in a posteriori (online) tests. Furthermore, we explore the data requirements of the time-series ML models by extracting physically relevant timescales of the complex system. We find that the ratio of the timescales of the minimum information required by the ML models to accurately capture the dynamics of the SST corresponds to the Reynolds number of the flow. The current framework provides the backbone to explore the capability of such models to capture the dynamics of high-dimensional complex dynamical system like SST flows.

97 MATHEMATICS AND COMPUTING↗

Level structure of light neutron-rich La isotopes beyond the 𝑁 = 82 shell closure

Here, the high spin excited states of Lanthanum isotopes 140–143 La, above the 𝑁 = 82 closed shell, were populated in fission reactions. The prompt 𝛾-ray transitions were measured using two complementary methods: (a) in coincidence with the isotopically identified fragments produced in the fission of the 238 U + 9 Be system using the Variable Mode Spectrometer (VAMOS++) and the Advanced Gamma Tracking Array (AGATA) spectrometer, and (b) high statistics threefold 𝛾−𝛾−𝛾 and fourfold 𝛾−𝛾−𝛾−𝛾 coincidence data from the spontaneous fission of 252 Cf using the Gammasphere. This work reports the first identification of a pair of parity doublet structures in 143 La and the new high spin level structure in 140–142 La from prompt 𝛾-ray spectroscopy. The level structures are interpreted in terms of the systematics of neighboring odd-𝑍 nuclei above the 𝑍 = 50 shell closure and large-scale shell model calculations. The present results indicate the presence of stable octupole deformation in 143 La. The excitation energy pattern and their comparison with neighboring isotones, moving away from the 𝑁 = 82 closed shell, point towards a transition from single-particle structures to an alternating parity rotational band structure in the La isotopic chain.

Navin, A. [Centre National de la Recherche Scienti↗

Electronic energy gap closure and metal-insulator transition in dense liquid hydrogen

Here, using quantum Monte Carlo (QMC) calculations, we investigate the insulator-metal transition observed in liquid hydrogen at high pressure. Below the critical temperature of the transition from the molecular to the atomic liquid, the fundamental electronic gap closure occurs abruptly, with a small discontinuity reflecting the weak first-order transition in the thermodynamic equation of state. Above the critical temperature, molecular dissociation sets in while the gap is still open. When the gap closes, the decay of the off-diagonal reduced density matrix shows that the liquid enters a gapless, but localized, phase: there is a crossover between the insulating and the metallic liquids. Compared to different density functional theory (DFT) functionals, our QMC calculations provide larger values for the fundamental gap and the electronic density of states close to the band edges, indicating that optical properties from DFT potentially benefit from error cancellations.

36 MATERIALS SCIENCE↗

Equation of state, phase transitions, and band-gap closure in Pb Cl 2 and Sn Cl 2

The equations of state and band-gap closures for Pb Cl 2 and Sn Cl 2 were studied using both experimental and theoretical methods. Here we measured the volume of both materials to a maximum pressure of 70 GPa using synchrotron-based angle-dispersive powder x-ray diffraction. The lattice parameters for both compounds showed anomalous changes between 16–32 GPa, providing evidence of a phase transition from the cotunnite structure to the related Co 2 Si structure, in contrast to the postcotunnite structure as previously suggested. First-principles calculations confirm this finding and predict a second phase transition to a Co 2 Si-like structure between 75– 110 GPa in Pb Cl 2 and 60–75 GPa in Sn Cl 2 . Band gaps were measured under compression to ~70 GPa for Pb Cl 2 and ~66 GPa for Sn Cl 2 and calculated up to 200 GPa for Pb CI 2 and 120 GPa for Sn Cl 2 . We find an excellent agreement between our experimental and theoretical results when using the Heyd-Scuseria-Ernzerhof (HSE06) hybrid functional, which suggests that this functional could reliably be used to calculate the band gap of similar AX 2 compounds. Experimental and calculated band-gap results show discontinuous decreases in the band gap corresponding to phase changes to higher-coordinated crystal structures, giving insight into the relationship between interatomic geometry and metallicity.

36 MATERIALS SCIENCE↗

Point-particle drag, lift, and torque closure models using machine learning: Hierarchical approach and interpretability

Developing deterministic neighborhood-informed point-particle closure models using machine learning has garnered interest recently from the dispersed multiphase flow community. The robustness of neural models for this complex multibody problem is hindered by the availability of particle-resolved data. Here, the present work addresses this unavoidable limitation of data paucity by implementing two strategies: (1) by using a rotation and reflection equivariant neural network and (2) by pursuing a physics-based hierarchical machine learning approach. The resulting machine-learned models are observed to achieve a maximum accuracy of 85% and 96% in the prediction of neighbor-induced force and torque fluctuations, respectively, for a wide range of Reynolds number and volume fraction conditions considered. Furthermore, we pursue force and torque network architectures that provide universal prediction spanning a wide range of Reynolds number (0.25 ≤ Re ≤250) and particle volume fraction (0 ≤ φ ≤0.4). The hierarchical nature of the approach enables improved prediction of quantities such as streamwise torque, by going beyond binary interactions to include trinary interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Energy Gap Closure of Crystalline Molecular Hydrogen with Pressure

Here, we study the gap closure with pressure of crystalline molecular hydrogen. The gaps are obtained from grand-canonical quantum Monte Carlo methods properly extended to quantum and thermal crystals, simulated by coupled electron ion Monte Carlo methods. Nuclear zero point effects cause a large reduction in the gap (~2 eV). Depending on the structure, the fundamental indirect gap closes between 380 and 530 GPa for ideal crystals and 330–380 GPa for quantum crystals. Beyond this pressure the system enters into a bad metal phase where the density of states at the Fermi level increases with pressure up to ~450–500 GPa when the direct gap closes. Our work partially supports the interpretation of recent experiments in high pressure hydrogen.

74 ATOMIC AND MOLECULAR PHYSICS↗

Proton Distribution Radii of 16–24 O: Signatures of New Shell Closures and Neutron Skin

The root mean square radii of the proton density distribution in 16–24 O derived from measurements of charge changing cross sections with a carbon target at ~900A MeV together with the matter radii portray thick neutron skin for 22–24 O despite 22,24 O being doubly magic. Imprints of the shell closures at N=14 and 16 are reflected in local minima of their proton radii that provide evidence for the tensor interaction causing them. The radii agree with ab initio calculations employing the chiral NNLO sat interaction, though skin thickness predictions are challenged. Furthermore, shell model predictions agree well with the data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Atwood effects on nonlocality of the scalar transport closure in Rayleigh-Taylor mixing

The importance of nonlocality is assessed in modeling mean scalar transport for turbulent Rayleigh-Taylor (RT) mixing at different Atwood numbers. Building on the two-dimensional incompressible work of Lavacot et al. [J. Fluid Mech. 985, A47 (2024)], the present work extends the macroscopic forcing method to variable density problems in three-dimensional space to measure moments of the generalized eddy diffusivity kernel in RT mixing for increasing Atwood numbers (𝐴 = 0.05, 0.3, 0.5, 0.8). It is found that as 𝐴 increases, (1) the eddy diffusivity moments become asymmetric and (2) the higher-order eddy diffusivity moments become larger relative to the leading-order diffusivity, indicating that nonlocality becomes more important at higher 𝐴. There is a particularly strong temporal nonlocality at higher 𝐴, suggesting stronger history effects. In conclusion, the implications of these findings for closure modeling for finite-Atwood RT are discussed.

general physics↗

Seniority Structure in Neutron-Rich Nucleus 128 Ag : Evidence for Robustness of 𝑁 = 82 Shell Closure in Silver Isotopes

The spectroscopic studies of very neutron-rich nucleus 128 Ag have been performed for the first time at the Radioactive Isotope Beam Factory of RIKEN. A new seniority isomer with a half-life of 1.60(7) μs has been identified and is proposed to have a spin-parity of 16 - with a maximally aligned configuration comprising three proton holes in the g 9/2 orbital and one neutron hole in the h 11/2 orbital. The new level structure in 128 Ag is quite well described by shell model calculations without invoking excitations across the Z = 50 and N = 82 shell gaps, and presents a good case of seniority scheme in odd-odd nuclei in the south vicinity of the double-magic nucleus 132 Sn. With a classification of various components of the proton-neutron interaction, the inversion of lowest-lying 9 - and 10 - states between 128 Ag and its neighboring isotone 130 In is found to be dynamically ascribed to the seniority-nonconserving proton-neutron interaction components. The structure above 10 - up to the 16 - isomer in 128 Ag shows remarkable similarities to seniority structures in the semimagic nuclei 128 Pd and 130 Cd. These spectroscopic features in 128 Ag indicate that the N = 82 shell closure is still robust in silver isotopes.

Luo, D. W. [Peking University, Beijing (China); et↗

Closure Schemes for Stably Stratified Atmospheric Flows without Turbulence Cutoff

Two recently proposed turbulence closure schemes are compared against the conventional Mellor–Yamada (MY) model for stably stratified atmospheric flows. The Energy- and Flux-Budget (EFB) approach solves the budgets of turbulent momentum and heat fluxes and turbulent kinetic and potential energies. The Cospectral Budget (CSB) approach is formulated in wavenumber space and integrated across all turbulent scales to obtain flow variables in physical space. Unlike the MY model, which is subject to a “critical gradient Richardson number,” both EFB and CSB models allow turbulence to exist at any gradient Richardson number R i and predict a saturation of flux Richardson number ( R f → R fm ) at sufficiently large R i . The CSB approach further predicts the value of R fm and reveals a unique expression linking the Rotta and von Kármán constants. Hence, all constants in the CSB model are nontunable and stability independent. All models agree that the dimensionless sensible heat flux decays with increasing R i . However, the decay rate and subsequent cutoff in the MY model appear abrupt. The MY model further exhibits an abrupt cutoff in the turbulent stress normalized by vertical velocity variance, while the CSB and EFB models display increasing trends. The EFB model produces a rapid increase in the ratio of turbulent potential energy and vertical velocity variance as R fm is approached, suggesting a strong self-preservation mechanism. Vertical anisotropy in the turbulent kinetic energy is parameterized in different ways in MY and EFB, but this consideration is not required in CSB. Differences between EFB and CSB model predictions originate from how the vertical anisotropy is specified in the EFB model.

54 ENVIRONMENTAL SCIENCES↗

A Scale-Adaptive Turbulent Kinetic Energy Closure for the Dry Convective Boundary Layer

Abstract A pragmatic scale-adaptive turbulent kinetic energy (TKE) closure is proposed to simulate the dry convective boundary layer for a variety of horizontal grid resolutions: from 50 m, typical of large-eddy simulation models that use three-dimensional turbulence parameterizations/closures, up to 100 km, typical of climate models that use one-dimensional turbulence and convection parameterizations/closures. Since parameterizations/closures using the TKE approach have been frequently used in these two asymptotic limits, a simple method is proposed to merge them with a mixing-length-scale formulation for intermediate resolutions. This new scale-adaptive mixing length naturally increases with increasing grid length until it saturates as the grid length reaches mesoscale-model resolution. The results obtained using this new approach for dry convective boundary layers are promising. The mean vertical profiles of potential temperature and heat flux remain in good agreement for different resolutions. A continuous transition (in terms of resolution) across the gray zone is illustrated through the partitioning between the model-resolved and the subgrid-scale transports as well as by documenting the transition of the subgrid-scale TKE source/sink terms. In summary, a natural and continuous transition across resolutions (from 50 m to 100 km) is obtained, for dry convection, using exactly the same atmospheric model for all resolutions with a simple scale-adaptive mixing-length formulation.

Kurowski, Marcin J.↗

Concrete Testing Case for Closure of Handford’s 241-C Underground Storage Tanks

Sixty percent of the nation's highly toxic and radioactive mixed wastes are stored at Hanford in 177 deteriorating underground storage tanks. To close or remove these storage tanks from service and place them in a condition that is protective of human health and the environment, the tanks must be physically stabilized to prevent subsidence once wastes have been retrieved. Remaining residual liquid waste in the tanks that cannot be removed must be solidified and the solid wastes encapsulated to meet the Nuclear Regulatory Commission, Department of Energy, Environmental Protection Agency, and the State of Washington requirements. The Department of Energy has developed cementitious flowable concretes to restrict access and provide chemical stabilization for radionuclides. Formulation, laboratory, and field testing for application at Hanford began with flowable, self-leveling structural and non-structural fills. A slump flow equal to or greater than 610 mm, 0% bleed water, and 0.1% (by volume) shrinkage measurements were key parameters guiding reformulation efforts that resulted in highly flowable, self-consolidating concretes that met Hanford 241-C Tank closure short- and long-term regulatory and engineering performance requirements.

formula development↗

Fracture Sustainability Pressure, Temperature, Differential Pressure, and Aperture Closure Data

In these data sets, the experiment time, actual date and time, room temperature, sample temperature, upstream and downstream pressures (measured independently), corrected differential pressure (measured independently and corrected for offset and room temperature) indication of aperture closure by linear variable differential transformer are presented. An indication of the sample is in the file name and in the first line of data.

15 GEOTHERMAL ENERGY↗

Validation of CTF Droplet Entrainment and Annular/Mist Closure Models using Riso Steam/Water Experiments

This report summarizes the work done to validate the droplet entrainment and de-entrainment models as well as two-phase closure models in the CTF code by comparison with experimental data obtained at Riso National Laboratory. The Riso data included a series of over 250 steam/water experiments that were performed in both tube and annulus geometries over a range of various pressures and outlet qualities. Experimental conditions were set so that the majority of cases were in the annular/mist ow regime. Measurements included liquid lm ow rate, droplet ow rate, lm thickness, and two-phase pressure drop. CTF was used to model 180 of the tubular geometry cases, matching experimental geometry, outlet pressure, and outlet ow quality to experimental values. CTF results were compared to the experimental data at the outlet of the test section in terms of vapor and entrained liquid ow fractions, pressure drop per unit length, and liquid lm thickness. The entire process of generating CTF input decks, running cases, extracting data, and generating comparison plots was scripted using Python and Matplotlib for a completely automated validation process. All test cases and scripting tools have been committed to the COBRA-TF master repository and selected cases have been added to the continuous testing system to serve as regression tests. The dierences between the CTF- and experimentally-calculated ow fraction values were con- sistent with previous calculations by Wurtz, who applied the same entrainment correlation to the same data. It has been found that CTF's entrainment/de-entrainment predictive capability in the annular/mist ow regime for this particular facility is comparable to the licensed industry code, COBRAG. While lm and droplet predictions are generally good, it has been found that accuracy is diminished at lower ow qualities. This nding is consistent with the noted deciencies in the Wurtz entrainment model employed by CTF. The CTF predicted two-phase pressure drop in the annular/mist ow regime has been found to be highly inaccurate, exhibiting a clear bias with respect to the experimental data. This inaccuracy led to an investigation that revealed deciencies in the implementation of the annular/mist interfacial friction model, which should be investigated further in the future. Looking to published COBRAG results for this same facility reveal it exhibits no bias with regard to experimental pressure drop results. In addition to the problems with pressure drop prediction, the lm thickness was also signicantly under-predicted by CTF compared to both experimental data and Wurtz's analytical calculations. Film thickness is calculated using a simple geometric relationship and lm void fraction in CTF, which is dependent on slip ratio and interfacial friction. It is possible that the issues aecting the pressure drop and lm void prediction are related.

97 MATHEMATICS AND COMPUTING↗