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

Level scheme study of Mo 91 : Weak-coupling approximation in the N = 50 region

The level scheme of 91 Mo has been studied using the GALILEO γ-ray spectrometer at Laboratori Nazionali di Legnaro, INFN. Five transitions together with four levels are identified for the first time by the γ-γ coincidence measurement, considerably enriching the negative band in 91 Mo. Furthermore, the discrepancy in the previous parity assignments for the states higher than 3807.5 keV has been clarified with the help of new experimental results. A shell-model calculation employing the JUN45 interaction is performed in order to get a further understanding of the level structure of 91 Mo. Being close to N=50 shell closure, the weak-coupling approximation between the valence g 9/2 neutron-hole with the even-even core ( 92 Mo) is discussed in 91Mo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Observation of collective modes of excitations in 59 Co, 59 Ni, and 61 Co and the influence of the $g_{9/2}$ orbital

High-spin states in 59 Co (Z = 27), 59 Ni (Z = 28), and 61 Co have been populated by the fusion evaporation reactions, 48 Ti( 14 C, p2n ) 59 Co, 48 Ti( 14 C, 3 n ) 59 Ni, and 50 Ti( 14 C, p2n ) 61 Co. The 9 MV tandem accelerator at the John D. Fox Laboratory, Florida State University (FSU) was used to accelerate the 14 C beam and the deexciting γ rays were detected by the FSU detector array consisting of six high-purity germanium (HPGe) clover detectors and three single crystals. Directional correlation of the $γ$ rays deexciting oriented states (DCO ratios) and polarization asymmetry measurements helped to establish spin and parities of the excited states whenever possible. The level scheme of 59 Co has been expanded with the inclusion of positive-parity states up to 31/2 + at around 11 MeV. The 59 Ni positive-parity states known from previous study were verified with modifications to some of the spins and parities. On the other hand, the negative-parity states were extended to 31/2 at an excitation energy of 12 MeV. No new transition was observed for 61 Co, but one of the major bands has been reassigned as consisting of positive-parity states by reason of this study. Importantly, excitations observed within the f 7/2 , p 3/2 , f 5/2 , and p 1/2 orbitals, and also across the N = 40 subshell closure into the g 9/2 orbital was established by comparison with large-scale shell-model calculations for the three nuclei studied.

59 ≤ A ≤ 89↗
DOE OSTI2024

Direct cross-section measurement of the weak r -process Sr 88 ( α , n ) Zr 91 reaction in ν -driven winds of core-collapse supernovae

About half of the heavy elements beyond iron are known to be produced by the rapid neutron capture process, known as the r process. However, the astrophysical site producing the r process is still uncertain. Chemical abundances observed in several cosmic sites indicate that different mechanisms should be at play. For instance, the abundances around silver measured in a subset of metal-poor stars indicate the presence of a weak r process. This process may be active in neutrino-driven winds of core collapse supernovae where (α,n) reactions dominate the synthesis of Ζ ≈ 40 elements in the expelled materials. Scarcely measured, the rates of (α,n) reactions are determined from statistical Hauser-Feshbach calculations with α-optical-model potentials, which are still poorly constrained. Further, the uncertainties of the (α,n) reaction rates therefore make a significant contribution to the uncertainties of the abundances determined from stellar modeling. In this work, the 88 Sr (α,n)⁢ 91 Zr reaction which impacts the weak r-process abundances has been probed at astrophysics energy for the first time; directly measuring the total cross sections at astrophysical energies of 8.37–13.09 MeV in the center of mass (3.8–7.5 GK). Two measurements were performed at ATLAS with the electrically segmented ionization chamber MUSIC, in inverse kinematics, while following the active target technique. The cross sections of this α-induced reaction on 88 Sr, located at the shell closure N = 50, have been found to be lower than expected, by a factor of 3, despite recent statistical calculations validated by measurements on neighboring nuclei. This result encourages more experimental investigations of (α,n) reactions, at Ν = 50 and towards the neutron-rich side, to further test the predictive power and reliability of such calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

New isomeric transition in Mg 36 : Bridging the N = 20 and N = 28 islands of inversion

Here, we observed a new isomeric γ transition at 168 keV in 36 Mg, with a half-life of T1/2=90⁢($^{+410}_{-50}$) ns. We propose that the observed transition deexcites a new 0 + isomeric state at 833 keV and populates the previously known first 2 + state. The existence of this isomer is consistent with the predictions of the large-scale shell-model calculations of 36 Mg using the sdpf-u-mix interaction. The observed excitation energy of the second 0 + state is caused by the small energy separation between two prolate-deformed configurations where the intruder configuration corresponds to two-neutron excitations from the sd to the pf shell. Within this interpretation, 36 Mg becomes the crossing point between nuclei in which ground state deformed/superdeformed configurations are caused by the dominance of N=20 intruders ( 32,34 Mg) and nuclei where deformed configurations are associated with the breaking of the N=28 closure and a large occupancy of the 1⁢p 3/2 neutron orbit ( 38 Mg and beyond). We found the lack of three-body monopole corrections in other effective interactions results in a predominance of N=20 intruder configurations past 38 Mg incompatible with our observation. We conclude that 36 Mg bridges the N=20 and N=28 islands of inversion, forming the so-called big island of deformation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Yrast structures of 98 Nb and 99 Mo

Neutron-rich nuclei in the A ≈ 100 mass region are interesting due to a rapid shape transition, especially pronounced in the Zr isotopes, and more recently observed in Nb isotopes. 98 Nb, with only one proton and one neutron outside the subshell closure nucleus of 96 Zr, is amenable to shell model calculations. To further examine the rapid shape transition the yrast structure of 98 Nb was established in this work. This was the only yrast structure missing from all immediate neighbors to 96 Zr. The yrast structure of 98 Nb was studied in the fission of the compound systems formed in three heavy-ion induced reactions, 24 Mg (134.5 MeV) + 173 Yb, 23 Na (129 MeV) + 176 Yb, and 18 O (91 MeV) + 208 Pb. Prompt γ-ray spectroscopy was accomplished using the Gammasphere array. Excitation energies up to 3 MeV were observed for the first time in 98 Nb. The yrast structure above the previously known (5) + isomer was established. In the process of studying 98 Nb the yrast structure of positive-parity states in 99 Mo was extended to 3.7 MeV excitation energy, the previously-known 99 Nb level scheme was enriched, and two new levels were added in the level scheme of 97 Zr. The coupling of the odd proton occupying the g 9/2 orbital to the yrast states in the core nucleus of 97 Zr can account for all observed states in 98 Nb. The yrast structure for the positive-parity states of 99 Mo is compared to the deformed ground-state bands of the 101 Ru isotone and of 98,100 Mo.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Improved structure of calcium isotopes from ab initio calculations

The in-medium similarity renormalization group (IMSRG) is a powerful and flexible many-body method to compute the structure of nuclei starting from nuclear forces. Recent developments have extended the IMSRG from its standard truncation at the normal-ordered two-body level, the IMSRG(2), to a precision approximation including normal-ordered three-body operators, the IMSRG(3)-N 7 . This improvement provides a more precise solution to the many-body problem and makes it possible to quantify many-body uncertainties in IMSRG calculations. We explore the structure of 44,48,52 Ca using the IMSRG(3)-N 7 , focusing on understanding existing discrepancies of the IMSRG(2) to experimental results. We find a significantly better description of the first 2 + excitation energy of 48 Ca, improving the description of the shell closure at N=28. At the same time, we find that the IMSRG(3)-N 7 corrections to charge radii do not resolve the systematic underprediction of the puzzling large charge radius difference between 52 Ca and 48 Ca. We present estimates of many-body uncertainties of IMSRG(2) calculations applicable also to other systems based on the size extensivity of the method.

39 ≤ A ≤ 58↗
DOE OSTI2021

Local wave-number model for inhomogeneous two-fluid mixing

We analyze the local wave-number (LWN) model, a two-point spectral closure model for turbulence, as applied to the Rayleigh-Taylor (RT) instability, the flow induced by the relaxation of a statically-unstable density stratification. Model outcomes are validated against data from 3D simulations of the RT instability. In the first part of the study we consider the minimal model terms required to capture inhomogeneous mixing and show that this version, with suitable model coefficients, is sufficient to capture the evolution of important mean global quantities including mix-width, turbulent mass flux velocity, and Reynolds stress, if the start time is chosen such that the earliest transitions are avoided. However, this simple model does not permit the expected finite asymptote of the density-specific-volume covariance b. In the second part of the study, we investigate two forms for a source term for the evolution of the spectrum of density-specific-volume covariance for the LWN model. The first includes an empirically motivated calibration of the source to achieve the final asymptotic state of constant b. Here, the second form does not require calibration but, in conjunction with enhanced diffusion and drag captures the full evolution of all the dynamical quantities, namely, the mix-layer growth, turbulent mass-flux velocity, Reynolds stress, as well as the desired behavior of b.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Fully kinetic simulations of strong steady-state collisional planar plasma shocks

We report on simulations of strong, steady-state collisional planar plasma shocks with fully kinetic ions and electrons, independently confirmed by two fully kinetic codes (an Eulerian continuum and a Lagrangian particle-in-cell). While kinetic electrons do not fundamentally change the shock structure as compared with fluid electrons, we find an appreciable rearrangement of the preheat layer, associated with nonlocal electron heat transport effects. The electron heat-flux profile qualitatively agrees between kinetic- and fluid-electron models, suggesting a certain level of “stiffness,” though substantial nonlocality is observed in the kinetic heat flux. We also find good agreement with nonlocal electron heat-flux closures proposed in the literature. Finally, in contrast to the classical hydrodynamic picture, we find a significant collapse in the “precursor” electric-field shock at the preheat layer leading edge, which correlates with the electron-temperature gradient relaxation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Capturing membrane structure and function in lattice Boltzmann models

Here, we develop a mesoscopic approach to model the nonequilibrium behavior of membranes at the cellular scale. Relying on lattice Boltzmann methods, we develop a solution procedure to recover the Nernst-Planck equations and Gauss's law. A general closure rule is developed to describe mass transport across the membrane, which is able to account for protein-mediated diffusion based on a coarse-grained representation. We demonstrate that our model is able to recover the Goldman equation from first principles and show that hyperpolarization occurs when membrane charging dynamics are controlled by multiple relaxation timescales. The approach provides a promising way to characterize non-equilibrium behaviors that arise due to the role of membranes in mediating transport based on realistic three-dimensional cell geometries.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Machine learning for the identification of phase transitions in interacting agent-based systems: A Desai-Zwanzig example

Deriving closed-form analytical expressions for reduced-order models, and judiciously choosing the closures leading to them, has long been the strategy of choice for studying phase- and noise-induced transitions for agent-based models (ABMs). In this paper, we propose a data-driven framework that pinpoints phase transitions for an ABM—the Desai-Zwanzig model—in its mean-field limit, using a smaller number of variables than traditional closed-form models. To this end, we use the manifold learning algorithm Diffusion Maps to identify a parsimonious set of data-driven latent variables, and we show that they are in one-to-one correspondence with the expected theoretical order parameter of the ABM. We then utilize a deep learning framework to obtain a conformal reparametrization of the data-driven coordinates that facilitates, in our example, the identification of a single parameter-dependent ordinary differential equation (ODE) in these coordinates. Additionally, we identify this ODE through a residual neural network inspired by a numerical integration scheme (forward Euler). We then use the identified ODE—enabled through an odd symmetry transformation—to construct the bifurcation diagram exhibiting the phase transition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Direct statistical simulation of the Lorenz96 system in model reduction approaches

Direct statistical simulation (DSS) of nonlinear dynamical systems bypasses the traditional route of accumulating statistics by lengthy direct numerical simulations by solving the equations that govern the statistics themselves. DSS suffers, however, from the curse of dimensionality as the statistics (such as correlations) generally have higher dimensions than the underlying dynamical variables. Here we investigate two approaches to reduce the dimensionality of DSS, illustrating each method with numerical experiments with the Lorenz96 dynamical system. The forms of DSS chosen here involve approximate closures at second and third order in the equal-time cumulants. We demonstrate significant reduction in computational effort that can be achieved without sacrificing the accuracy of DSS. The methods developed here can be applied to turbulent fluid and magnetohydrodynamical systems. Published by the American Physical Society 2025

Li, Kuan↗
DOE OSTI2021

Lagrangian analysis for turbulent transport in variable-density turbulence

In this report Lagrangian analysis of materially conserved scalars is applied to the problem of turbulent transport in variable-density flows. The consequences of an additional material conserved quantity, the density, is generally not acknowledged and leads to significant and meaningfully different expressions for turbulent transport in the moment equations. The formal Lagrangian analysis produces gradient transport expressions substantially different from those obtained by the physically intuitive “argument by analogy” method used in computational models. Various intuitive arguments, in Favre and Reynolds averaged settings, are contrasted to the formal Lagrangian results. Using expressions from the formal analysis, we derive consistent gradient transport closures for the turbulent transport terms that appear in the first- and second-order Favre moment equations. Results for coupled multispecies turbulent transport are given. The analysis is limited to variable-density turbulence in which the dilatation of the fluctuating velocity is small. The results are applicable to turbulent combustion and to stellar convection problems in which the density fluctuations are on the order of the mean density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Internal energy balance and aerodynamic heating predictions for hypersonic turbulent boundary layers

The elemental equation governing heat transfer in aerodynamic flows is the internal energy equation. For a boundary layer flow, a double integration of the Reynolds-averaged form of this equation provides an expression of the wall heat flux in terms of the integrated effects, over the boundary layer, of various physical processes: turbulent dissipation, mean dissipation, turbulent heat flux, etc. Recently available direct numerical simulation data for a Mach 11 cold-wall turbulent boundary layer allows a comparison of the exact contributions of these terms in the energy equation to the wall heat flux with their counterparts modeled in the Reynolds-averaged Navier-Stokes (RANS) framework. Various approximations involved in RANS, both closure models as well as approximations involved in adapting incompressible RANS models to a compressible form, are assessed through examination of the internal energy balance. There are a number of potentially problematic assumptions and terms identified through this analysis. Here, the effect of compressibility corrections of the dilatational dissipation type is explored, as is the role of the modeled turbulent dissipation, in the context of wall heat flux predictions. The results indicate several potential avenues for RANS model improvement for hypersonic cold-wall boundary-layer flows.

42 ENGINEERING↗
DOE OSTI2020

Charge Radius of the Short-Lived Ni 68 and Correlation with the Dipole Polarizability

We present the first laser spectroscopic measurement of the neutron-rich nucleus 68 Ni at the N = 40 subshell closure and extract its nuclear charge radius. Since this is the only short-lived isotope for which the dipole polarizability $α$ D has been measured, the combination of these observables provides a benchmark for nuclear structure theory. We compare them to novel coupled-cluster calculations based on different chiral two- and three-nucleon interactions, for which a strong correlation between the charge radius and dipole polarizability is observed, similar to the stable nucleus 48 Ca. Three-particle–three-hole correlations in coupled-cluster theory substantially improve the description of the experimental data, which allows to constrain the neutron radius and neutron skin of 68 Ni.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Two-Neutron Halo is Unveiled in F 29

Here, we report the measurement of reaction cross sections ( σ R ex ) of F 27 , 29 with a carbon target at RIKEN. The unexpectedly large σ R ex and derived matter radius identify F 29 as the heaviest two-neutron Borromean halo to date. The halo is attributed to neutrons occupying the 2 p 3 / 2 orbital, thereby vanishing the shell closure associated with the neutron number N = 20 . The findings are explained by state-of-the-art shell model calculations. Coupled-cluster computations based on effective field theories of the strong nuclear force describe the matter radius of F 27 but are challenged for F 29 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Solving the Puzzles of the Decay of the Heaviest Known Proton-Emitting Nucleus 185 Bi

Two long-standing puzzles in the decay of 185 Bi, the heaviest known proton-emitting nucleus are revisited. These are the non-observation of the 9/2 – state, which is the ground state of all heavier odd-A Bi isotopes, and the hindered nature of proton and α decays of its presumed 60-μs 1/2 + ground state. The 185 Bi nucleus has now been studied with the 95 Mo( 93 Nb; 3n) reaction in complementary experiments using the Fragment Mass Analyzer and Argonne Gas-Filled Analyzer at Argonne National Laboratory’s ATLAS facility. The experiments have established the existence of two states in 185 Bi; the short-lived T 1/2 = $2.8$$^{+2.3}_{–1.0}$ μs, proton- and α-decaying ground state, and a 58(2)–μs γ-decaying isomer, the half-life of which was previously attributed to the ground state. The reassignment of the ground-state lifetime results in a proton-decay spectroscopic factor close to unity and represents the only known example of a ground-state proton decay to a daughter nucleus ( 184 Pb) with a major shell closure. Furthermore, the data also demonstrate that the ordering of low- and high-spin states in 185 Bi is reversed relative to the heavier odd-A Bi isotopes, with the intruder-based 1/2 + configuration becoming the ground state, similar to the lightest At nuclides.

150 ≤ A ≤ 189↗
DOE OSTI2021

First Observation of the Four-Proton Unbound Nucleus 18 Mg

18 Mg was observed, for the first time, by the invariant-mass reconstruction of 14 O + 4p events. The ground-state decay energy and width are E T = 4.865(34) MeV and Γ = 115(100) keV, respectively. The observed momentum correlations between the five particles are consistent with two sequential steps of prompt 2p decay passing through the ground state of 16 Ne. The invariant-mass spectrum also provides evidence for an excited state at an excitation energy of 1.84(14) MeV, which is likely the first excited 2 + state. As this energy exceeds that for the 2 + state in 20 Mg, this observation provides an argument for the demise of the N = 8 shell closure in nuclei far from stability. Furthermore, in open systems this classical argument for shell strength is compromised by Thomas-Ehrman shifts.

6 ≤ A ≤ 19↗
DOE OSTI2022

β-Delayed One and Two Neutron Emission Probabilities Southeast of 132 Sn and the Odd-Even Systematics in r-Process Nuclide Abundances

Here, the β-delayed one- and two-neutron emission probabilities (P 1n and P 2n ) of 20 neutron-rich nuclei with N≥82 have been measured at the RIBF facility of the RIKEN Nishina Center. P 1n of 130,131 Ag, 133,134 Cd, 135,136 In, and 138,139 Sn were determined for the first time, and stringent upper limits were placed on P2n for nearly all cases. β-delayed two-neutron emission (β2n) was unambiguously identified in 133 Cd and 135,136 In, and their P 2n were measured. Weak β2n was also detected from 137,138 Sn. Our results highlight the effect of the N=82 and Z=50 shell closures on β-delayed neutron emission probability and provide stringent benchmarks for newly developed macroscopic-microscopic and self-consistent global models with the inclusion of a statistical treatment of neutron and γ emission. The impact of our measurements on r-process nucleosynthesis was studied in a neutron star merger scenario. Our P 1n and P 2n have a direct impact on the odd-even staggering of the final abundance, improving the agreement between calculated and observed Solar System abundances. The odd isotope fraction of Ba in r-process-enhanced (r-II) stars is also better reproduced using our new data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗