Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “level set”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Emulating the Lyman-Alpha forest 1D power spectrum from cosmological simulations: new models and constraints from the eBOSS measurement

We present the Lyssa suite of high-resolution cosmological simulations of the Lyman-α forest designed for cosmological analyses. These 18 simulations have been run using the Nyx code with 40963 hydrodynamical cells in a 120 Mpc (∼ 81 Mpc/h) comoving box and individually provide sub-percent level convergence of the Lyman-α forest 1d flux power spectrum. We build a Gaussian process emulator for the Lyssa simulations in the lym1d likelihood framework to interpolate the power spectrum at arbitrary parameter values. We validate this emulator based on leave-one-out tests and based on the parameter constraints for simulations outside of the training set. We also perform comparisons with a previous emulator, showing a percent level accuracy and a good recovery of the expected cosmological parameters. Using this emulator we derive constraints on the linear matter power spectrum amplitude and slope parameters A Lyα and n Lyα . While the best-fit Planck ΛCDM model has A Lyα = 8.79 and n Lyα = -2.363, from DR14 eBOSS data we find that A Lyα < 7.6 (95% CI) and n Lyα = -2.369 ± 0.008. The low value of A Lyα , in tension with Planck, is driven by the correlation of this parameter with the mean transmission of the Lyman-α forest. This tension disappears when imposing a well-motivated external prior on this mean transmission, in which case we find A Lyα = 9.8 ± 1.1 in accordance with Planck.

Walther, Michael↗
DOE OSTI2025

Polaris-PARCS Sensitivity Study on LWR Fuel Cycles: Polaris Input Options

This study is the first of a multi-phase effort to assess the sensitivity of light-water reactor (LWR) core-level prediction biases to changes in lattice-level calculation parameters. Prediction bias is the measured-to-predicted difference in a core-level quantity of interest (QOI) which can be estimated by comparing the simulation results with the plant-measured data for key nuclear parameters. The LWR two-step neutronics codes employed herein are the SCALE–Polaris lattice physics code (v6.3.1) and the Purdue Advanced Reactor Core Simulator (PARCS) nodal diffusion simulator (v3.4.2), both funded and used for confirmatory analysis to support licensing by the US Nuclear Regulatory Commission (NRC). Polaris–PARCS is used to model Watts Bar Unit 1 cycles 1–3 and Peach Bottom Unit 2 cycles 1–3. This study focuses on the impact of changes to Polaris input options such as scattering treatment or quadrature settings and how these input options induce changes in core-level quantities of interest (QOIs)bias. The report documents multiple bias assessments for different modeling choices and compares the bias magnitude to the QOI measurement uncertainties. Future companion reports will investigate the sensitivity of core-level LWR prediction bias to Polaris input options and Polaris-computed QOIs such as few-group assembly-homogenized cross sections to gain an understanding of the key drivers of prediction bias at lattice and core levels for application of a two-step LWR neutronics procedure in a licensing scenario.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Multi-component Dark Matter and small scale structure formation

We consider the evolution of non-thermal dark matter perturbations in models which contain both Weakly Interacting Massive Particles (WIMPs) and axions. Using constraints from existing observations we examine the percentage of WIMPs and axions that may comprise the cosmological dark matter budget in models with an Early Matter Dominated Epoch (EMDE) — where entropy production is important. After carefully tracking the thermal evolution of the various species by solving the Boltzmann equations, we consider the enhancement of perturbations that may have led to early structure formation for axions and WIMPs. We investigate the impact of enhanced perturbations on the parameter space of both species, after imposing existing constraints from indirect detection experiments. Given these constraints we establish the feasibility of axions to form miniclusters in the early universe in EMDEs for a given percentage of allowed WIMPs. We find that EMDEs with low reheat temperatures near the BBN limit are preferred for axion minicluster formation. When the EMDE is caused by string moduli, the WIMP contribution to the relic density is set by the moduli branching to dark matter at the level of ≲ $\mathcal{O}$ (10 −4 ).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Search for $ {\textrm{B}}_{\textrm{s}}^0 $ → ℓ ∓ τ ± with the semi-leptonic tagging method at Belle

We present a search for the lepton-flavor-violating decays ${B}_s^0$ → ℓ ∓ τ ± , where ℓ = e, μ, using the full data sample of 121 fb -1 collected at the Υ(5S) resonance with the Belle detector at the KEKB asymmetric-energy e + e - collider. We use ${B}_s^0{\overline{B}}_s^0$ events in which one ${B}_s^0$ meson is reconstructed in a semileptonic decay mode and the other in the signal mode. We find no evidence for (${B}_s^0$ → ℓ ∓ τ ± decays and set upper limits on their branching fractions at 90% confidence level as $\mathcal{B}$(${B}_s^0$ → e ∓ τ ± ) < 14 × 10 -4 and $\mathcal{B}$ (${B}_s^0$ → μ ∓ τ ± ) < 7.3 × 10 -4 . Our result represents the first upper limit on the ${B}_s^0$ → e ∓ τ ± decay rate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A sequential Attacker-Defender game for distribution systems resilience enhancement against extreme weather events

Improving distribution system resilience against frequent extreme weather events is important for reliable power system operations. Especially when dealing with events such as hurricanes that have short-term predictions, proactive pre-event preparedness plays a vital role in system resilience performance. In this paper, we propose a novel approach to construct pre-event resource allocation plans for system operators to cope with upcoming threats through a sequential attacker-defender game framework. The sequential attacker-defender game is designed to model the interaction between the extreme weather and the system operator. In each round of the game, the attacker and the defender sequentially update their current strategies by accounting for the opponent’s action set. The attacker model is formulated as a bi-level problem to identify the severe outage scenarios, and the defender model is formulated as a two-stage optimization problem to determine the allocation of restoration resources including mobile responsive resources and repair crews. Two scale-reduction strategies are proposed to ensure the scalability of the game scheme. Finally, case studies on the IEEE 33-bus and a 7149-node practical utility system validate the effectiveness of the proposed sequential game and the efficiency of the scale-reduction strategies.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Slot-die-coating operability windows for polymer electrolyte membrane fuel cell cathode catalyst layers

Roll-to-roll (R2R) slot-die coating of polymer electrolyte membrane fuel cell (PEMFC) catalyst layers represents a scalable deposition method for producing 10–20 m 2 ·min -1 of catalyst-coated gas diffusion layers (GDLs). This high-throughput production technique will help lower the cost of PEMFC catalyst layers. The uniformity of the wet layer applied by slot die deposition is affected by process parameters such as substrate speed, vacuum pressure applied at the upstream meniscus, gap between the slot die lips and substrate, ink rheology, and other ink and substrate properties. The set of conditions for producing a defect-free coating with a dilute ink typically requires little to no upstream vacuum pressure, so suitable operating conditions can be found easily through trial and error and operator intuition. However, the higher viscosity of more concentrated inks dramatically shifts the range of settings that result in a homogeneous coating to higher vacuum levels, which are harder to find through hit or miss. A predictive model showing the range of operable conditions decreases material wastage inherent in experimentally searching for suitable parameters. In this study, the defect-free coating parameter window is explored experimentally and theoretically for two concentrations of PEFC cathode inks. Both a full capillary hydrodynamic model and a computationally cheaper viscocapillary model successfully predict the experimentally determined coating window within the experimental and model uncertainty limits for inks with 5.3 wt% and 12.0 wt% solids ink while maintaining the 0.1 mg Pt ·cm -2 Pt areal loading target. Here, this paper demonstrates a viable pathway for meeting the $30/kW net ultimate cost target of the United States Department of Energy (U.S. DOE) Hydrogen Fuel Cells Technologies Office (HFTO). The concentrated ink lowers the thermal energy and capital expenditure (CapEx) budget of the coating process by decreasing the amount of time, energy, and floorspace required for drying the coating.

25 ENERGY STORAGE↗
DOE OSTI2020

Two-phase flow measurements and observations in a cooling panel of the reactor cavity cooling system

The analysis of the flow behavior within a cooling panel of a scaled water-cooled reactor cavity cooling system has been performed supported by high-resolution measurements of the void fraction at the risers’ outlet. The results have shown a stable, symmetric distribution of the flow through the risers during subcooled boiling conditions, characterized by small, slow bubbles with average void fraction lower than 0.3 for all risers. As the temperature increased to saturation, flow distribution appeared strongly asymmetrical and unstable, with peaks of void fraction as high as 0.9. Large, faster bubbles were observed during these peaks, immediately followed by flow of subcooled liquid water. Flow stagnation and inversion has also been observed. The high resolution void fraction results, complemented by the measurements of the coolant flow rate and temperatures, not only have further described the behavior of the system under hypothetical accident conditions, but represent a unique set of data to support the validation of system-level codes, and more advanced computational tools.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Prediction of Redox Potentials for U, Np, Pu, and Am in Aqueous Solution

The redox properties of the actinides in aqueous solution are important for fuel production/reprocessing and understanding the environmental impact of nuclear waste. The redox potentials for U, Np, Pu, and Am in oxidation states from 0 up to VII (as appropriate) in aqueous solutions have been predicted at the density functional theory level with the B3LYP functional, Stuttgart small core pseudopotential basis sets for the actinides, and explicit (30H 2 O molecules)/implicit treatment of the aqueous solvent using the self-consistent reaction field COSMO and SMD approaches for the implicit solvation. The predictions of the structural parameters of clusters incorporating first and second solvation shells are consistent with the available experimental data., Our results are typically within 0.2 V of the available experimental data using two explicit solvation shells with an implicit solvent model. The use of the PW91 functional substantially improved the prediction of the Pu(VI/V) redox couple. The redox couples for An(VI/IV) and An(V/IV) which involve the addition of protons and removal of the actinyl oxygens led to slightly larger differences from experiment. Here, the An(IV/0) and An(III/0) couples were reliably predicted with our approach. Predictions of the unknown An(II/I) redox potentials were negative, consistent with expectations, and predictions for unknown An(VII/VI), An(III/II), and An(II/0) redox couples improve prior estimates.

Actinides↗
DOE OSTI2022

Elucidating the Atomic Structures of the Gel Layer Formed during Aluminoborosilicate Glass Dissolution: An Integrated Experimental and Simulation Study

The altered glasses produced during aqueous dissolution of silicate and borosilicate glasses are among the most complex structures to understand at the atomic level due to their amorphous nature, random porosity and various levels of hydration. In this study, we gained insights of the complex atomic structures of altered aluminoborosilicate glasses by combining a range of experimental and computational approaches. The altered glasses were prepared by dissolution of three glasses with varying level of alumina in acid for 7 days. A comprehensive set of experimental [elemental analysis, high-energy X-ray diffraction, 29 Si and 27 Al solid-state nuclear magnetic resonance (NMR), and O 1s X-ray photoelectron spectroscopy (XPS)] and modeling (molecular dynamics (MD) simulations using non-reactive and reactive force fields) were used to study the atomic structures of these altered glasses. Elemental analysis showed that most of the B in the pristine glasses was leached into the solution and was not contained in the altered glass. The 29 Si and 27 Al solid-state NMR spectra revealed that the altered glasses have more polymerized silicate network as compared to those in the pristine glasses due to the reformation of linkages among Si and Al oxygen polyhedral in the altered glasses. The bridging and non-bridging (or hydroxyl O) atoms in the altered glasses were also quantified from their O 1s XPS spectra. Atomic structure models of the altered glasses were constructed by using MD simulations using the reactive force field based on the compositional information obtained from experiments. Various pore structures were generated using the charge-scaling (CS) method by using different initial densities and CS temperatures; the best CS parameters of each alternated glass were then determined by comparing with the experimental structure factors obtained from high energy X-ray diffraction. Pore and the atomic structures and vibrational properties around these pore surfaces were analyzed. Furthermore, these results from this comprehensive study thus provides a realistic insight into the pore morphology, atomic structure, and vibrational properties of altered glasses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Guest-Host Interactions in Clathrate Hydrates: Benchmark MP2 and CCSD(T)/CBS Binding Energies of CH4, CO2 and H2S in (H2O)20 Cages

We present benchmark binding energies of naturally occurring gas molecules CH4, CO2, and H2S in the small cage, namely the pentagonal dodecahedron (512) (H2O)20, which is one of the constituent cages of the 3 major lattices (structures I, II and H) of clathrate hydrates. These weak interactions require higher levels of electron correlation and converge slowly with increasing basis set to the Complete Basis Set (CBS) limit, necessitating the use of large basis sets up to the augcc- pV5Z and subsequent correction for Basis Set Superposition Error (BSSE). For the host hollow (H2O)20 cages we have identified a most stable isomer with binding energy of -200.8 ± 2.1 kcal/mol at the CCSD(T)/CBS limit (-199.2 ± 0.5 kcal/mol at the MP2/CBS limit). Additionally, we report converged second order Moller-Plesset (MP2) CBS binding energies for the encapsulation of guests in the (H2O)20 cage of -4.3 ± 0.1 for CH4@(H2O)20, -6.6 ± 0.1 for CO2@(H2O)20 and -8.5 ± 0.1 kcal/mol for H2S@(H2O)20, respectively. For CH4@(H2O)20, exhibiting the weakest encapsulation affinity among the three, we report CCSD(T)/aug-cc-pVTZ binding energies and, based on them, a CCSD(T)/CBS estimate of -4.75 ± 0.1 kcal/mol. To the best of our knowledge, the CCSD(T)/aug-cc-pVTZ calculation for CH4@(H2O)20 is the largest one reported to date (168 valence electrons, 1978 basis functions and the correlation of 84 doubly occupied and 1873 virtual orbitals) and required a scalable implementation of the (T) module on 6144 nodes (350208 cores) of the “Cori” supercomputer at the National Energy Research Supercomputing Center (NERSC) for a total execution time of 195 minutes (for the (T) part). These efficient scalable implementations of highly correlated methods offer the capability to obtain long-lasting benchmarks of intermolecular interactions in complex systems. They also provide a path towards parametrizing classical potentials needed to study the dynamical and transport properties in these complex systems as well as assess the accuracy of lower scaling electronic structure methods such as Density Functional Theory (DFT) and MP2 including its spin-biased variants.

Heindel, Joseph↗
DOE OSTI2021

Experimental and computational investigation of the bond energy of thorium dicarbonyl cation and theoretical elucidation of its isomerization mechanism to the thermodynamically most stable isomer, thorium oxide ketenylidene cation, OTh + CCO

Collision-induced dissociation (CID) of [Th,2C,2O] + with Xe is performed using a guided ion beam tandem mass spectrometer (GIBMS). The only products observed are ThCO + and Th + by sequential loss of CO ligands. The experimental findings and theoretical calculations support that the structure of [Th,2C,2O] + is the bent homoleptic thorium dicarbonyl cation, Th + (CO) 2 , having quartet spin, which is both thermodynamically and kinetically stable enough in the gas phase to be observed in our GIBMS instrument. Analysis of the kinetic energy-dependent cross sections for this CID reaction yields the first experimental determination of the bond dissociation energy (BDE) of (CO)Th + –CO at 0 K as 1.05 ± 0.09 eV. A theoretical BDE calculated at the CCSD(T) level with cc-pVXZ (X = T and Q) basis sets and a complete basis set (CBS) extrapolation is in very good agreement with the experimental result. Although the doublet spin bent thorium oxide ketenylidene cation, OTh + CCO, is calculated to be the most thermodynamically stable structure, it is not observed in our experiment where [Th,2C,2O] + is formed by association of Th + and CO in a direct current discharge flow tube (DC/FT) ion source. Potential energy profiles of both quartet and doublet spin are constructed to elucidate the isomerization mechanism of Th + (CO) 2 to OTh + CCO. The failure to observe OTh + CCO is attributed to a barrier associated with C–C bond formation, which makes OTh + CCO kinetically inaccessible under our experimental conditions. Furthermore, chemical bonding patterns in low-lying states of linear and bent Th + (CO) 2 and OTh + CCO isomers are also investigated.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

On the Impact of Bus Dwelling on Macroscopic Fundamental Diagrams

Network macroscopic fundamental diagrams (MFDs) have recently been shown to exist in real-world urban traffic networks. When present, MFDs can be used to model traffic dynamics within an urban network by dividing the network into a set of spatially compact homogeneous regions and tracking the average level of congestion in each region. Existing analytical methods to estimate MFD mostly focus on the behavior of a single type of vehicle and do not capture the patterns of mixed traffic (e.g., cars and buses). The existence of buses matters since a bus will block the movements of other vehicles when it dwells at the bus stop. This paper proposes an analytical method to estimate the impact of bus dwelling on a network’s MFD based on the network’s geometric features, traffic control strategies, and bus operation parameters, and validates the performance of the proposed method using simulations based on microscopic traffic models. Comparisons of the analytical and simulation results show that the proposed analytical method can generally provide a good estimate of the lower bound and upper bound of the network’s MFD.

Xu, Guanhao↗
DOE OSTI2020

A modular ultra-high vacuum millikelvin scanning tunneling microscope

We describe the design, construction, and performance of an ultra-high vacuum (UHV) scanning tunneling microscope (STM) capable of imaging at dilution-refrigerator temperatures and equipped with a vector magnet. The primary objective of our design is to achieve a high level of modularity by partitioning the STM system into a set of easily separable, interchangeable components. This naturally segregates the UHV needs of STM instrumentation from the typically non-UHV construction of a dilution refrigerator, facilitating the usage of non-UHV materials while maintaining a fully bakeable UHV chamber that houses the STM. Here, the modular design also permits speedy removal of the microscope head from the rest of the system, allowing for repairs, modifications, and even replacement of the entire microscope head to be made at any time without warming the cryostat or compromising the vacuum. Without using cryogenic filters, we measured an electron temperature of 184 mK on a superconducting Al(100) single crystal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Search for decays of the Higgs boson into pair-produced pseudoscalar particles decaying into τ+τ−τ+τ− using pp collisions at s=13 TeV with the ATLAS detector

A search for a pair of low-mass pseudoscalars a that promptly decay into τ-leptons is presented using 140 fb −1 of proton–proton collision data at 13 TeV centre-of-mass energy recorded with the ATLAS detector at the Large Hadron Collider. The result is used to place constraints on exotic decays of the Higgs boson into four τ-leptons, H→aa→τ+τ−τ+τ−. This search focuses on events with either one or two τ-leptons decaying into hadrons and neutrinos, and the remaining three or two τ-leptons decaying into either electron or muon and neutrinos. No significant excess is observed above the expected Standard Model background and upper limits at the 95% confidence level on B(H→aa→τ+τ−τ+τ−) are set ranging from 0.06 to 0.23, depending on the mass ma ranging from 15 to 60 GeV.

Source record↗
DOE OSTI2022

Search for the weak decay ψ (3686) → $Λ^{+}_{c}$ $\overline{Σ}$ - + c.c. *

Using (448.1 ± 2.9) x 10 6 ψ (3686) events collected with the BESIII detector, we perform the first search for the weak baryonic decay ψ (3686) → $Λ^{+}_{c}$ $\overline{Σ}$ - + c.c. . The analysis procedure is optimized using a blinded method. No significant signal is observed, and the upper limit on the branching fraction ($\mathcal{B}$) of ψ (3686) → $Λ^{+}_{c}$ $\overline{Σ}$ - + c.c. is set as 1.4 x 10 -5 at the 90% confidence level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Helical magnetic order and Fermi surface nesting in noncentrosymmetric ScFeGe

Here, an investigation of the structural, magnetic, thermodynamic, and charge transport properties of noncentrosymmetric hexagonal ScFeGe reveals it to be an anisotropic metal with a transition to a weak itinerant incommensurate helimagnetic state below T N =36 K. Neutron diffraction measurements discovered a temperature and field independent helical wave vector k = (0 0 0.193) with magnetic moments of 0.53 μ B per Fe confined to the ab-plane. Density functional theory calculations are consistent with these measurements and find several bands that cross the Fermi level along the c-axis with a nearly degenerate set of flat bands just above the Fermi energy. The anisotropy found in the electrical transport is reflected in the calculated Fermi surface, which consists of several warped flat sheets along the c-axis with two regions of significant nesting, one of which has a wave vector that closely matches that found in the neutron diffraction. The electronic structure calculations, along with a strong anomaly in the c-axis conductivity at TN, signal a Fermi surface driven magnetic transition, similar to that found in spin density wave materials. Magnetic fields applied in the ab-plane result in a metamagnetic transition with a threshold field of ≈ 6.7 T along with a sharp, strongly temperature dependent discontinuity and a change in sign of the magnetoresistance for in-plane currents. Thus, ScFeGe is an ideal system to investigate the effect of in-plane magnetic fields on a helimagnet with a c-axis propagation vector, where the relative strength of the magnetic interactions and anisotropies determine the topology and magnetic structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Bardeen-Cooper-Schrieffer pairing of composite fermions

Topological pairing of composite fermions has led to remarkable ideas, such as excitations obeying non-Abelian braid statistics and topological quantum computation. Here, we construct a p-wave paired Bardeen-Cooper-Schrieffer (BCS) wave function for composite fermions in the torus geometry, which is a convenient geometry for formulating momentum space pairing as well as for revealing the underlying composite-fermion Fermi sea. Following the standard BCS approach, we minimize the Coulomb interaction energy at half filling in the lowest and the second Landau levels, which correspond to filling factors ν = 1/2 and ν = 5/2 in GaAs quantum wells, by optimizing two variational parameters that are analogous to the gap and the Debye cut-off energy of the BCS theory. Our results show no evidence for pairing at ν = 1/2 but a clear evidence for pairing at ν = 5/2. To a good approximation, the highest overlap between the exact Coulomb ground state at ν = 5/2 and the BCS state is obtained for parameters that minimize the energy of the latter, thereby providing support for the physics of composite-fermion pairing as the mechanism for the 5/2 fractional quantum Hall effect. We discuss the issue of modular covariance of the composite-fermion BCS wave function, and calculate its Hall viscosity and pair correlation function. By similar methods, we look for but do not find an instability to s-wave pairing for a spin-singlet composite-fermion Fermi sea at half-filled lowest Landau level in a system where the Zeeman splitting has been set to zero.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Langevin dynamics of generalized spins as SU( N ) coherent states

Classical models of spin systems traditionally retain only the dipole moments, but a quantum spin state will frequently have additional structure. Spins of magnitude S have N = 2⁢S + 1 levels. Alternatively, the spin state is fully characterized by a set of N 2 - 1 local physical observables, which we interpret as generalized spin components. For example, a spin with S = 1 has three dipole components and five quadrupole components. Further, these components evolve under a generalization of the classical Landau-Lifshitz dynamics, which can be extended with noise and damping terms. In this paper, we reformulate the dynamical equations of motion as a Langevin dynamics of SU(N) coherent states in the Schrödinger picture. This viewpoint is especially useful as the basis for an efficient numerical method to sample spin configurations in thermal equilibrium and to simulate the relaxation and driven motion of topological solitons. To illustrate the approach, we simulate a nonequilibrium relaxation process that creates CP 2 skyrmions, which are topological defects with both dipole and quadrupole character.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗