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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Improving the Physical Basis for Updraft Dynamics in Deep Convection Parameterizations

This article presents a new deep convective parameterization that determines cloud characteristics based on a specified cloud size distribution. The vertical profiles of cloud properties are determined by analytic equations, which formulate entrainment with an inverse relationship to cloud width. In line with recent studies of large eddy simulations (LES), cloud widths are assumed to be constant with height and vertical mass flux (M) characteristics of the clouds are therefore regulated by the vertical velocity profile. The parameterization is configured to work with existing cloud base M closure formulations, with the closure predicting the total cloud area rather than the cloud base M directly. Analytic formula are also used to connect the vertical wind shear magnitude to the cloud size distribution, wherein larger shear magnitudes result in more numerous large updrafts than weaker shear magnitudes, which is in line with recent research results. The parameterization is compared against 10 deep convective LES with varying thermodynamic and vertical wind shear profiles. Results show dramatic improvements in the prediction of normalized M, detrainment, and the properties of detrained air over the existing Zhang and McFarlane (1995) scheme. In particular, the new model is able to correctly portray the transition from a bottom-heavy M profile in weakly sheared environments, to a top-heavy M profile in strongly sheared environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Linking Deep and Shallow Convective Mass Fluxes via an Assumed Entrainment Distribution in CAM5‐CLUBB: Parameterization and Simulated Precipitation Variability

Abstract We modify the Zhang‐McFarlane (ZM) deep convection scheme in the Community Atmosphere Model version 5 to couple it with a unified parameterization for boundary‐layer turbulence and shallow convection, that is, Cloud Layers Unified by Binormals (CLUBB). By assuming a lognormal distribution of entrainment rate across the entire moist convective regimes, we link mass fluxes between shallow and deep convection, which are partitioned by the entrainment rate of the shallowest deep convective plume. Hence, a new deep convective closure is established which is coupled to the sub‐grid vertical motion variability in CLUBB. The convection feedback (or memory) effects are also considered to decrease the entrainment spectrum width and enhance the vertical velocity variability that further affect deep convection. Results show that the revised scheme improves the precipitation simulations in terms of the mean state and variability at various timescales, such as the alleviated double‐intertropical convergence zone and more realistic simulations of the seasonal variation of monsoon precipitation over East Asia, Madden‐Julian Oscillation, and precipitation diurnal phase propagations downstream of large terrains. The improvements are still seen in many aspects such as the mean‐state precipitation when turning off the convection feedback impacts in the revised scheme, emphasizing the benefits of using the modified mass‐flux closure. However, the convection feedbacks have considerable effects on the precipitation diurnal cycle simulations over regions with late‐afternoon precipitation peaks. Overall, the revised scheme provides a unified treatment for sub‐grid vertical motions across regimes of boundary‐layer turbulence, shallow convection, and deep convection, leading to better‐simulated precipitation at various timescales.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Automated Symbolic Upscaling: 1. Model Generation for Extended Applicability Regimes

Abstract In porous media theory, upscaling techniques are fundamental to deriving rigorous Darcy‐scale models for flow and reactive transport in subsurface systems. Due to limitations in classical techniques, a number of ad hoc approaches have been proposed to address physical regimes in which reactive time scales are similar to, or faster than, diffusive time scales. In Part 1 of this two part series, we present a strategy for expanding the applicability of classical homogenization theory by generalizing the assumed closure form. We detail the implementation of this strategy on two reactive mass transport problems with moderately reactive physics. The strategy produces nontrivial homogenized models with emergent terms and effective parameters that couple reactive, diffusive, and advective transport. The differences in equation forms between the macroscopic and pore‐scale descriptions advise caution to further studies where the forms of macroscopic equations are assumed, as opposed to rigorously derived. Numerical validation is provided for each problem to show that the error estimates of homogenization theory are satisfied, and to justify the implemented strategy. In Part 2, the presented strategy is automated using symbolic computing to expedite its implementation.

Pietrzyk, Kyle↗

DOE OSTI2021

Benchmarking NIMROD continuum kinetic formulations through the steady-state poloidal flow

In this work, continuum kinetic formulations are employed as a mechanism to include closure physics in an extended magnetohydrodynamics model. Two continuum kinetic approaches have been implemented in the plasma fluid code NIMROD [Sovinec et al., “Nonlinear magnetohydrodynamics with high-order finite elements,” J. Comput. Phys. 195, 355 (2004)] including a Chapman–Enskog-like (CEL) formulation and a more conventional df approach. Ion kinetic closure schemes are employed to describe the neoclassical flow properties in axisymmetric toroidal geometry. In particular, predictions for steady-state values of poloidal flow profiles in tokamak geometry are provided using both the df formulation and two different solution techniques for the CEL approach. These results are benchmarked against analytic theory predictions as well as results from the drift kinetic code DK4D. The continuum kinetic formulations employed here show agreement with both the analytic theory and DK4D results, and offer a novel velocity space representation involving higher-order finite elements in pitch angle.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Development of five-moment two-fluid modeling for Z-pinch physics

The Z-pinch m = 0 instability as well as its stabilization by radially sheared axial flow is studied using the nonlinear ideal five-moment twofluid (5M2F) model with an extension of that model to include Braginskii heat and momentum transport. Using the ideal 5M2F model, linear growth rate results are compared with prior work using MHD and Hall MHD. At small normalized wavenumber, 1 < k a < 4 , where a is the effective pinch radius, 5M2F results agree with Hall MHD within 20% in scenarios without radially sheared axial flow. With the sheared flow and focusing on ka = 10/3, agreement with Hall MHD is excellent. In the limit of small ion inertial length, results also match with MHD. A comparison with PIC modeling of shear-free m = 0 stability focuses on a plasma scenario based on recent experimental results. In a scan of mode wavenumber, ideal 5M2F results are qualitatively similar to PIC: the growth rate rises to a peak at a moderate wavenumber and declines at a large wavenumber in contrast to MHD results, which show the saturation of the growth rate with the increasing wavenumber rather than a decline. The peak normalized 5M2F growth rate is γ τ A = 1.5 , where sA is the Alfven transit time across the pinch. The peak occurs at normalized wavenumber ka = 10. For comparison, PIC results have a peak growth of γ τ A = 0.77 at ka = 5. Including Braginskiibased closure of the 5M2F model does not qualitatively change the ideal results in this particular case. Nonlinear saturation is studied using the 5M2F model with the dissipative Braginskii-based closure in cases with pinch-edge sheared-flow speed equal to half the Alfven speed. Nonlinear mixing due to the sheared flow yields a quasi-steady state after modest losses of pinch ion inventory and pinch thermal energy, approximately 30% and 10%, respectively. 5M2F modeling captures the essential physics of m = 0 instability and offers a computationally tractable route to high-fidelity modeling of 3D Z-pinch behavior, including m = 1 instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Compton rocket effect due to the action of radiation reaction force in degenerate plasma

A closed set of fluid equations with radiation reaction force (RRF) are constructed from the moments of the appropriate single particle kinetic equation describing a relativistic degenerate (high density) electron plasma. The closure, in analogy with the Maxwellian closure for non-degenerate plasmas, is affected via a parametrized Fermi-Dirac distribution. It is shown that the degeneracy increases RRF just as will be predicted from the so-called “Compton Rocket” effect.

Physics↗

DOE OSTI2024

A Review of Remote Welding and Nondestructive Examination Technologies for the DOE Standard Canister

The U.S. Department of Energy (DOE) manages a wide variety of spent nuclear fuel (SNF) that poses a unique management challenge. To help address this challenge, the DOE Standard Canister (DOESC), designed to remain sealed during handling, storage, transportation, and disposal, was conceptualized as a standardized containment vessel to accommodate DOE-managed SNF. Since 1999, several welding and examination processes have been independently developed for the DOESC’s closure welds. However, neither the DOESC nor these processes have been realized in an operational capacity. This review paper seeks to present and compare previously developed DOESC closure weld, nondestructive examination, and repair processes and technologies. Specific processes developed for the Idaho Spent Fuel Facility, in preparation for the Yucca Mountain geological repository, and the recent Road-Ready Demonstration Project are discussed. Further, specific focus is given to how different operating constraints and the American Society of Mechanical Engineers Boiler and Pressure Vessel Code (BPVC) have driven certain welding and nondestructive examination requirements. Historical DOESC welding and examination strategies are assessed against current regulatory and BPVC requirements. The comparison of welding processes, technologies, and DOESC designs presented in this review paper will inform further construction efforts for other commercial and DOE-managed SNF containments, including the DOESC.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2020

Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach

Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.

42 ENGINEERING↗

DOE OSTI2023

Bayesian Parameter Estimation of the k-ω Shear Stress Transport Model for Accurate Simulations of Impinging-Jet Heat Transfer

Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

High temperature superconductivity in the candidate phases of solid hydrogen

As the simplest element in nature, unraveling the phase diagram of hydrogen is a primary task for condensed matter physics. As conjectured many decades ago, in the low-temperature and high-pressure part of the phase diagram, solid hydrogen is expected to become metallic with a high superconducting transition temperature. The metallization may occur via band gap closure in the molecular solid or via a transition to the atomic solid. Recently, a few experimental studies pushed the achievable pressures into the 400–500 GPa range. There are strong indications that at some pressure in this range metallization via either of these mechanisms occurs, although there are disagreements between experimental reports. Furthermore, there are multiple good candidate crystal phases that have emerged from recent computational and experimental studies which may be realized in upcoming experiments. Therefore, it is crucial to determine the superconducting properties of these candidate phases. In a recent study, we reported the superconducting properties of the C2/c-24 phase, which we believe to be a strong candidate for metallization via band gap closure. Here, we report the superconducting properties of the Cmca-12, Cmca-4 and I4 1 /amd-2 phases including the anharmonic effects using a Wannier function-based dense k-point and q-point sampling. Here we find that the Cmca-12 phase has a superconducting transition temperature that rises from 86 K at 400 GPa to 212 K at 500 GPa, whereas the Cmca-4 and I4 1 /amd-2 phases show a less pressure-dependent behavior with their T c in the 74–94 K and 307–343 K ranges, respectively. These properties can be used to distinguish between crystal phases in future experiments. Understanding superconductivity in pure hydrogen is also important in the study of high-T c hydrides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2020

Large-Eddy Simulations of Idealized Atmospheric Boundary Layers Using Nalu-Wind

Accurate prediction of wind-plant performance relies, in part, on properly characterizing the turbulent atmospheric boundary layer (ABL) flow in which wind turbines operate. Large-eddy simulation (LES) is a powerful tool for simulating ABLs because it resolves the largest, most energetic scales of three-dimensional turbulent motions. Yet LES predictions are well known to depend on modeling choices such as grid resolution, numerical discretization schemes, and closures for unresolved scales of turbulence. Here, we evaluate how these choices influence predictions of ABL winds using Nalu-Wind, a wind-specific fork of the open-source, generalized, unstructured, massively parallel flow solver NaluCFD/Nalu.

17 WIND ENERGY↗

DOE OSTI2026

Emergence of vorticity and viscous stress in finite-scale quantum hydrodynamics

The Madelung equations offer a hydrodynamic description of quantum systems, from single particles to quantum fluids. In this formulation, the probability density is mapped onto the fluid density and the phase is treated as a scalar potential generating the velocity field. As examples of potential flows, quantum fluids described in this way are inherently irrotational, but quantum vortices may arise at discrete points where the phase is undefined. In this paper, starting from this irrotational description of a quantum fluid, a coarse-graining procedure is applied to arrive at a macroscopic description of the quantum fluid in terms of a hierarchy of moments in which the role of velocity is played by a Favre average of the microscopic velocity field. This hierarchy is truncated using an explicit closure derived from an expansion in a finite length scale. The resulting coarse-grained fields are shown to allow for finite vorticity at any point in the fluid. Additionally, it is shown that this vorticity obeys a similar equation to the vorticity equation in classical hydrodynamics and includes a vortex-stretching term. The particular closure employed here also gives rise to a novel stress term in the fluid equations, which in the appropriate limit appears analogous to an artificial viscous stress from computational fluid dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2020

Prompt-delayed γ-ray spectroscopy of neutron-rich 119.121 In isotopes

The Z = 50 shell closure, near N = 82, is unique in the sense that it is the only shell closure with the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , enclosing the magic gap. The interaction of the proton hole/particle in the above-mentioned orbitals with neutrons in the νh 11/2 orbital is an important prerequisite to the understanding of the nuclear structure near N = 82 and the νπ interaction. To explore the structural similarity between the high-spin isomeric states in In ( Z = 49), Sn ( Z = 50), and Sb ( Z = 51) isotopes from a microscopic point of view. In addition, to understand the role of a proton hole or particle in the spin-orbit partner orbitals, π g 9 / 2 and π g 7 / 2 , respectively, with neutron holes in the ν h 11 / 2 orbital on these aforementioned isomers. The fusion and transfer induced fission reaction Be 9 ( U 238 , f ) with 6.2 MeV/u beam energy, using a unique setup consisting of AGATA, VAMOS++, and EXOGAM detectors, was used to populate through the fission process and study the neutron-rich In 119 , 121 isotopes. This setup enabled the prompt-delayed γ -ray spectroscopy of isotopes in the time range of 100 ns – 200 μ s . In the odd- A In 119 , 121 isotopes, indications of a short half-life 19/2 - isomeric state, in addition to the previously known 25/2 + isomeric state, were observed from the present data. Further, new prompt transitions above the 25/2 + isomer in In 121 were identified along with reevaluation of its half-life. The experimental data were compared with the theoretical results obtained in the framework of large-scale shell-model calculations in a restricted model space. Furthermore, the ( π g 9 / 2 ν h 11 / 2 ; I | H ^ | π g 9 / 2 ν h 11 / 2 ; I ) two-body matrix elements of residual interaction were modified to explain the excitation energies and the B ( E 2) transition probabilities in the neutron-rich In isotopes. The (i) decreasing trend of E (29/2 + ) - E (25/2 + ) in odd-In (with dominant configuration π g 9 / 2 - 1 ν h 11 / 2 - 2 and maximum aligned spin of 29/2 + ) and (ii) increasing trend of E (27/2 + ) - E (23/2 + ) in odd-Sb (with dominant configuration π g 7 / 2 + 1 ν h 11 / 2 - 2 and maximum aligned spin of 27/2 + ) with increasing neutron number could be understood as a consequence of hole-hole and particle-hole interactions, respectively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

Spectroscopic quadrupole moments in 124 $\mathrm{Xe}$

The Xe isotopic chain with four valence protons above the Z=50 shell closure is an ideal laboratory for the study of the evolution of nuclear deformation. At the N=82 shell closure, 136 Xe presents all characteristics of a doubly closed shell nucleus with a spherical shape. In the very neutron-deficient isotopes close to N=50, the α-decay chain of Xe was investigated to probe the radioactive decay properties near the drip-line and the magicity of 100 Sn. Additionally, the Xe isotopes present higher order symmetries in the nuclear deformation such as the octupole degree of freedom near N=60 and N=90 or O(6) symmetry in stable isotopes. The relevance of the O(6) symmetry has been investigated by measuring the spectroscopic quadrupole moment of the first excited states in 124 Xe. In the O(6) symmetry limit, the spectroscopic quadrupole moment of collective states is expected to be null. A stable 124 Xe beam with energies of 4.03A MeV and 4.11 A MeV was used to bombard a nat W target at the GANIL facility. Excited states were populated via the safe Coulomb excitation reaction. The collision of the heavy ions with a large Z at low energy make this reaction sensitive to the diagonal E2 matrix element of the excited states. The recoils were detected in the VAMOS++ magnetic spectrometer and the γ rays in the AGATA tracking array. The least squares fitting code gosia was used for the analysis to extract both E2 and M1 transitional and E2 diagonal matrix elements. The rotational ground state band was populated up to the 8$^{+}_{1}$ state as well as the 2$^{+}_{2}$ and 4$^{+}_{2}$ states. Using high precision spectroscopic data to constrain the gosia fit, the spectroscopic quadrupole moments of the 2$^{+}_{1}$, 4$^{+}_{1}$, and 6$^{+}_{1}$ states were determined for the first time. The spectroscopic quadrupole moments were found to be negative, large, and constant in the ground state band underlining the prolate axially deformed ground state band of 124 Xe. Here, the present experimental data confirm that the O(6) symmetry is substantially broken in 124 Xe .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Nucleonic shells and nuclear masses

The binding energy of an isotope is a sensitive indicator of the underlying shell structure as it reflects the net energy content of a nucleus. Since magic nuclei are significantly lighter, or more bound, compared to their neighbors, the presence of nucleonic shell structure makes an imprint on nuclear masses. In this study, using a carefully designed binding-energy indicator, we catalog the appearance of spherical and deformed shell and subshell closures throughout the nuclear landscape. After presenting experimental evidence for shell and subshell closures as seen through the lens of nuclear masses, we study the ability of global nuclear mass models to predict local binding-energy variations related to shell effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Evolution of shell gaps in the neutron-poor calcium region from invariant-mass spectroscopy of 37,38 Sc, 35 Ca, 34 K

A fast secondary beam of 37 Ca impinged on a 9 Be target resulting in a set of reactions populating proton-rich nuclei including 35 Ca and the first observations of 37,38 Sc and 34 K. Invariant-mass spectroscopy, used to reconstruct proton decays for these nuclei, yielded three new ground-state masses and information on their low-lying structures. The newly measured mass excesses are: ΔM( 37 Sc) = 3500(410) keV, ΔM( 38 Sc) = –4656(14) keV, and ΔM( 34 K) = –1487(17) keV. These nuclei straddle the well-known Z = 20 shell closure as well as the N = 16 subshell closure. Furthermore, trends in separation energies help elucidate how nuclear structure evolves showing a fading of the Z = 20 shell gap for N ≥ 18 and indications of a N = 16 subshell gap.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Variations in the charge radii of indium isotopes between 𝑁=52 and 82

Measurements of the 5⁢𝑝 2 𝑃 3/2 → 9⁢𝑠 2 𝑆 1/2 and 5⁢𝑝 2 𝑃 1/2 → 8⁢𝑠 2 𝑆 1/2 transitions in the indium atom, combined with new atomic physics calculations, were used to extract the changes in mean-square nuclear charge radii, 𝛿⁡⟨𝑟 2 ⟩, of the indium (𝑍 = 49) isotopes 101–111 In , 113–123 In , and 125–131 In . With a proton hole in the closed nuclear shell of 𝑍 = 50, indium provides a detailed study of the effect of unpaired nucleons adjacent to the proton-shell closure, allowing investigation into the charge radii for isotopes between the two major neutron-shell closures at 𝑁 = 50 and 𝑁 = 82. A study of the variations in charge radii between neighboring isotopes with neutron number (the ‘odd-even staggering') is presented and provides further insight and challenges for the theoretical description of the size of proton-hole nuclei. Two nuclear theories, density functional theory and the valence-space in-medium similarity renormalization group method, were employed to interpret the data. The new information obtained in this work provides valuable insights into the successes and shortcomings of the theoretical approaches employed.

ab initio calculations↗

DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗