Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hard constraint”

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 91 records · Page 5

Microstructure Investigation and Mechanical Properties of Coated Zircaloy Cladding

Accident Tolerant Fuel concepts have been developed and tested in diverse research programs presently around the world. Industry teams have developed proprietary coating compositions and geometries, but significant knowledge gaps exist in terms of the resulting properties of the coating and the cladding and how those components impact each other in reactor conditions. Furthermore, the national laboratory and university teams currently studying the industry concepts are limited in their ability to freely disseminate and discuss results obtained on commercial coatings. The Advanced Fuels Campaign is therefore working to develop an understanding of ‘generic’ coatings applied to commercial nuclear cladding alloys free of these constraints. This milestone report provides insights in the microstructure of a 7-micrometer Cr coating on Zircaloy-2 and a 5- micrometer coating on Zry-4 as reference materials. SEM and TEM methods were utilized to describe and identify the interface, while nanoindentation was used to determine hardness and E-Moduli of coating and cladding material. This milestone summarizes the acquired knowledge to build a necessary verified foundation for future simulation and qualification efforts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Excited-state uncertainties in lattice-QCD calculations of multi-hadron systems

Excited-state effects lead to hard-to-quantify systematic uncertainties in lattice quantum chromodynamics (LQCD) spectroscopy calculations when computationally accessible imaginary times are smaller than inverse excitation gaps, as often arises for multi-hadron systems with signal-to-noise problems. Lanczos residual bounds address this by providing two-sided constraints on energies that do not require assumptions beyond Hermiticity, but often give very conservative systematic uncertainty estimates. Here, a more-constraining set of gap bounds is introduced for hadron spectroscopy. These bounds provide tighter constraints whose validity requires an explicit assumption about an energy gap. Exactly solvable lattice field theory correlators are used to test the utility of residual and gap bounds at finite and infinite statistics. Two-sided bounds and other analysis methods are then applied to a high-statistics LQCD calculation of nucleon-nucleon scattering at $m_π\sim 800$ MeV. Generalized eigenvalue problem (GEVP) and Lanczos energy estimators are compatible when applied to the same correlator data, but analyses including different interpolating operators show statistically significant inconsistencies. However, two-sided bounds from all operators are consistent. Under the assumption that the number of energy levels below $NΔ$ and $ΔΔ$ thresholds is the same as for non-interacting nucleons, gap bounds are sufficient to constrain nucleon-nucleon scattering amplitudes at phenomenologically relevant precision. Lanczos methods further reveal that energy-eigenstate estimates from previously studied asymmetric correlators have not converged over accessible imaginary times. Nevertheless, data-driven examples demonstrate why assumptions are required to draw conclusions about the natures of two-nucleon ground states at these masses.

Detmold, William [MIT, Cambridge, CTP]↗

On the gamma-ray emission from the core of the Sagittarius dwarf galaxy

ABSTRACT We use Fermi-LAT data to analyse the faint gamma-ray source located at the centre of the Sagittarius (Sgr) dwarf spheroidal galaxy. In the 4FGL-DR3 catalogue, this source is associated with the globular cluster, M54. We investigate the spectral energy distribution and spatial extension of this source, with the goal of testing two hypotheses: (1) the emission is due to millisecond pulsars within M54, or (2) the emission is due to annihilating dark matter from the Sgr halo. For the pulsar interpretation, we consider a two-component model which describes both the lower-energy magnetospheric emission and possible high-energy emission arising from inverse Compton scattering. We find that this source has a point-like morphology at low energies, consistent with magnetospheric emission, and find no evidence for a higher-energy component. For the dark matter interpretation, we find the signal favours a dark matter mass of mχ = 29.6 ± 5.8 GeV and an annihilation cross section of $\sigma v = (2.1 \pm 0.59) \times 10^{-26} \, \text{cm}^3$ s−1 for the $b \bar{b}$ channel (or mχ = 8.3 ± 3.8 GeV and $\sigma v = (0.90 \pm 0.25) \times 10^{-26} \, \text{cm}^3$ s−1 for the τ+τ− channel), when adopting a J-factor of $J=10^{19.6} \, \text{GeV}^2 \, \text{cm}^{-5}$. This parameter space is consistent with gamma-ray constraints from other dwarf galaxies and with dark matter interpretations of the Galactic Centre Gamma-Ray Excess.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimization-Driven Scenario Grouping

Scenario decomposition algorithms for stochastic programs compute bounds by dualizing all nonanticipativity constraints and solving individual scenario problems independently. Here, we develop an approach that improves on these bounds by reinforcing a carefully chosen subset of nonanticipativity constraints, effectively placing scenarios into groups. Specifically, we formulate an optimization problem for grouping scenarios that aims to improve the bound by optimizing a proxy metric based on information obtained from evaluating a subset of candidate feasible solutions. We show that the proposed grouping problem is NP-hard in general, identify a polynomially solvable case, and present two formulations for solving the problem: a matching formulation for a special case and a mixed-integer programming formulation for the general case. We use the proposed grouping scheme as a preprocessing step for a particular scenario decomposition algorithm and demonstrate its effectiveness in solving standard test instances of two-stage 0–1 stochastic programs. Using this approach, we are able to prove optimality for all previously unsolved instances of a standard test set. Additionally, we implement this scheme as a preprocessing step for PySP, a publicly available and widely used implementation of progressive hedging, and compare this grouping approach with standard grouping approaches on large-scale stochastic unit commitment instances. Finally, the idea is extended to propose a finitely convergent algorithm for two-stage stochastic programs with a finite feasible region.

97 MATHEMATICS AND COMPUTING↗

Probing the axion-electron coupling with NuSTAR observations of galaxies

We search for the existence of ultralight axions coupling to electrons and photons using data from the NuSTAR telescope directed toward the galaxies M82, M87, and M31. We focus on electron bremsstrahlung and Compton scattering for axion production in stars, summing over the stellar populations found in the target galaxies when computing the axion luminosity. We then compute the hard x-ray signal that arises from the conversion of these axions to photons in each galaxy’s magnetic fields, inferred from analog galaxies in cosmological magnetohydrodynamic simulations. Analyzing NuSTAR data toward these galaxies between roughly 20 to 70 keV, we find no evidence for axions and set leading constraints on the combined axion-electron and axion-photon coupling at the level of |𝑔 𝑎⁢𝑒⁢𝑒 ×𝑔 𝑎⁢𝛾⁢𝛾 | ≲ 8.3 ×10 −27 GeV −1 for 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence, with M82 providing the most stringent constraints.

Axion-like particles↗

Chalcogenide phase-change material advances programmable terahertz metamaterials: a non-volatile perspective for reconfigurable intelligent surfaces

Terahertz (THz) waves have gained considerable attention in the rising 6G communication due to their large bandwidth. However, the cost and power consumption become the major constraints for the commercialization of 6G THz systems as the frequency increases. Reconfigurable intelligent surface (RIS) comprising active metasurfaces and digital controllers has been proposed for beamforming in the 6G multiple-input-multiple-output systems, showing good potential to suppress the system size, weight, and power consumption (SWaP). Currently, their controlling diodes can hardly work up to THz frequencies. Therefore, several active stimuli have been investigated as alternatives. Among them, chalcogenide phase-change material Ge 2 Sb 2 Te 5 (GST) addresses large modulation depth, picosecond switching speed, and non-volatile properties. Notably, the non-volatile GST may enable RIS systems with memory and low control power. This work briefly reviews the advances of GST-tuned THz metamaterials (MTMs), discusses the current obstacles to overcome, and gives a perspective of GST applications in the rising 6G communications.

6G↗

Optimising the processing and storage of visibilities using lossy compression

The next-generation radio astronomy instruments are providing a massive increase in sensitivity and coverage, largely through increasing the number of stations in the array and the frequency span sampled. The two primary problems encountered when processing the resultant avalanche of data are the need for abundant storage and the constraints imposed by I/O, as I/O bandwidths drop significantly on cold storage. An example of this is the data deluge expected from the SKA Telescopes of more than 60 PB per day, all to be stored on the buffer filesystem. While compressing the data is an obvious solution, the impacts on the final data products are hard to predict. In this paper, we chose an error-controlled compressor – MGARD – and applied it to simulated SKA-Mid and real pathfinder visibility data, in noise-free and noise-dominated regimes. As the data have an implicit error level in the system temperature, using an error bound in compression provides a natural metric for compression. MGARD ensures the compression incurred errors adhere to the user-prescribed tolerance. To measure the degradation of images reconstructed using the lossy compressed data, we proposed a list of diagnostic measures, exploring the trade-off between these error bounds and the corresponding compression ratios, as well as the impact on science quality derived from the lossy compressed data products through a series of experiments. We studied the global and local impacts on the output images for continuum and spectral line examples. We found relative error bounds of as much as 10%, which provide compression ratios of about 20, have a limited impact on the continuum imaging as the increased noise is less than the image RMS, whereas a 1% error bound (compression ratio of 8) introduces an increase in noise of about an order of magnitude less than the image RMS. For extremely sensitive observations and for very precious data, we would recommend a 0.1% error bound with compression ratios of about 4. These have noise impacts two orders of magnitude less than the image RMS levels. At these levels, the limits are due to instabilities in the deconvolution methods. We compared the results to the alternative compression tool DYSCO, in both the impacts on the images and in the relative flexibility. MGARD provides better compression for similar error bounds and has a host of potentially powerful additional features.

Techniques: interferometric↗

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↗

Pulse and Polarization Structures in Axion-Converted X-Rays from Pulsars

Neutron stars (NSs) with their strong magnetic fields and hot dense cores could be powerful probes of axions, a classic benchmark of feebly coupled new particles, through abundant production of axions with the axion-nucleon coupling and subsequent conversion into x-rays due to the axion-photon coupling. In this Letter, we point out that the pulsation structures in both the intensity and polarization of x-rays from NSs could provide us additional information about axions and their couplings. We develop new analytical formalisms of pulsation-polarization structure applicable to a wide range of NSs in the axion scenario and argue that they hold in complicated astrophysical environments. As a case study, we apply our formalism to a representative X-ray Dim Isolated Neutron Star, RX J1856.6-3754, with an unexpected hard x-ray excess which might be axion induced. We show with an updated fit that the axion explanation is compatible with both the intensity and pulsation data available, and combining the pulsation data does not shift the posterior by more than 1⁢𝜎. Yet, the preferred parameter space is close to being excluded by other astrophysical constraints. With a 75% reduction of the uncertainties in the pulsation data, we could potentially draw a definite conclusion on the axion-induced x-rays at more than 3⁢𝜎 level.

Axions↗

Lyman continuum leaker candidates among highly ionised, low-redshift dwarf galaxies selected from He II

Contemporary research suggests that the reionisation of the intergalactic medium (IGM) in the early Universe was predominantly realised by star-forming (proto-)galaxies (SFGs). Due to observational constraints, our knowledge on the origins of sufficient amounts of ionising Lyman continuum (LyC) photons and the mechanisms facilitating their transport into the IGM remains sparse. Recent efforts have thus focussed on the study of local analogues to these high-redshift objects. We aim to acquire a set of very low-redshift SFGs that exhibit signs of a hard radiation field being present. A subsequent analysis of their emission line properties is intended to shed light on how the conditions prevalent in these objects compare to those predicted to be present in early SFGs that are thought to be LyC emitters (LCEs). We used archival spectroscopic SDSS DR12 data to select a sample of low-redshift He II 4686 emitters and restricted it to a set of SFGs with an emission line diagnostic sensitive to the presence of an active galactic nucleus, which serves as our only selection criterion. We performed a population spectral synthesis with FADO to reconstruct these galaxies’ star-formation histories (SFHs). Utilising the spectroscopic information at hand, we constrained the predominant ionisation mechanisms in these galaxies and inferred information on interstellar medium (ISM) conditions relevant for the escape of LyC radiation. Our final sample consists of eighteen ionised, metal-poor galaxies (IMPs). These low-mass (6.2 ≤ log(M * /M ⊙ ) ≤ 8.8), low-metallicity (7.54 ≤log(O/H) + 12 ≤ 8.13) dwarf galaxies appear to be predominantly ionised by stellar sources. We find large [OIII] 5007/[OII] 3727 ratios and [SII] 6717,6731/Hα deficiencies, which provide strong indications for these galaxies to be LCEs. At least 40% of these objects are candidates for featuring cosmologically significant LyC escape fractions ≳10%. The IMPs’ SFHs exhibit strong similarities and almost all galaxies appear to contain an old (> 1 Gyr) stellar component, while also harbouring a young, two-stage (~10 Myr and < 1 Myr) starburst, which we speculate might be related to LyC escape. The properties of the compact emission line galaxies presented here align well with those observed in many local LCEs. In fact, our sample may prove as an extension to the rather small catalogue of local LCEs, as the extreme ISM conditions we find are assumed to facilitate LyC leakage. Notably, all of our eighteen candidates are significantly closer (z < 0.1) than most established LCEs. If the inferred LyC photon loss is genuine, this demonstrates that selecting SFGs from He II 4686 is a powerful selection criterion in the search for LCEs.

79 ASTRONOMY AND ASTROPHYSICS↗

Analog and symbolic computation through the Koopman framework

We develop a Koopman operator framework for studying the computational structure of dynamical systems. Specifically, we show that the resolvent of the Koopman operator provides a natural abstraction of halting, yielding a ‘Koopman halting problem’ that is recursively enumerable in general. For symbolic systems, such as those defined on Cantor space, this operator formulation captures reachability between clopen sets, while for equicontinuous systems we prove that the Koopman halting problem is decidable. Our framework demonstrates that absorbing (halting) states in coarse-grained finite automata correspond to Koopman eigenfunctions with eigenvalue one, while cycles in the transition graph impose spectral constraints associated with periodic dynamics. These results provide a unifying perspective on computation in symbolic and analog systems, showing how computational universality is reflected in operator spectra, invariant subspaces, and algebraic structures. Beyond symbolic dynamics, this operator-theoretic lens opens pathways to analyze the computational properties of a broader class of dynamical systems, including polynomial and analog models, and suggests that computational hardness may admit dynamical signatures in terms of Koopman spectral structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

First detection of VHE gamma-ray emission from TXS 1515–273, study of its X-ray variability and spectral energy distribution

ABSTRACT We report here on the first multiwavelength (MWL) campaign on the blazar TXS 1515–273, undertaken in 2019 and extending from radio to very-high-energy gamma-rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotron-peaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6σ. The average integral VHE flux of the source is 6 ± 1 per cent of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMM–Newton, and NuSTAR. The long continuous X-ray observations were separated by ∼9 h, both showing clear hour scale flares. In the XMM–Newton data, both the rise and decay time-scales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability time-scales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme synchrotron-peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band.

Acciari, V. A.↗

Efficient Topology Design Algorithms for Power Grid Stability

The dynamic response of power grids to small disturbances influences their overall stability. This letter examines the effect of network topology on the linearized time-invariant dynamics of electric power systems. The proposed framework utilizes H 2 -norm based stability metrics to study the optimal placement of lines on existing networks as well as the topology design of new networks. The design task is first posed as an NP-hard mixed-integer nonlinear program (MINLP) that is exactly reformulated as a mixed-integer linear program (MILP) using McCormick linearization. To improve computation time, graph-theoretic properties are exploited to derive valid inequalities (cuts) and tighten bounds on the continuous optimization variables. Moreover, a cutting plane generation procedure is put forth that is able to interject the MILP solver and augment additional constraints to the problem on-the-fly. Finally, the efficacy of our approach in designing optimal grid topologies is demonstrated through numerical tests on the IEEE 39-bus network.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Axion lines from nuclear de-excitations in galactic stellar populations

We show that mono-energetic axions are produced in abundance through nuclear de-excitations in nearby galaxies such as M87, which is the central galaxy of the Virgo cluster, and the starburst galaxy M82. If the axion couples to both nucleons and photons and is ultralight, then monochromatic hard x-ray signatures are induced by the subsequent axion-to-photon conversion in the magnetic fields permeating these systems. We search for evidence of such signals using NuSTAR data, focusing specifically on the 57 Fe de-excitation line at 14.4 keV, and we catalog other potentially relevant nuclear lines. We find no evidence for axions from M87 or M82 and set stringent constraints on the combined axion-nucleon and axion-photon coupling at the level of |𝑔 𝑎⁢𝑛⁢𝑛 × 𝑔 𝑎⁢𝛾⁢𝛾 | ≲ 1.1 × 10 −22 GeV −1 in the limit 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence.

Ning, Orion [University of California, Berkeley, C↗

Integration of Soft Data Into Geostatistical Simulation of Categorical Variables

Uncertain or indirect “soft” data, such as geologic interpretation, driller’s logs, geophysical logs or imaging, offer potential constraints or “soft conditioning” to stochastic models of discrete categorical subsurface variables in hydrogeology such as hydrofacies. Previous bivariate geostatistical simulation algorithms have not fully addressed the impact of data uncertainty in formulation of the (co) kriging equations and the objective function in simulated annealing (or quenching). This paper introduces the geostatistical simulation code tsim-s, which accounts for categorical data uncertainty through a data “hardness” parameter. In generating geostatistical realizations with tsim-s, the uncertainty inherent to soft conditioning is factored into both 1) the data declustering and spatial correlation functions in cokriging and 2) the acceptance probability for change of category in simulated quenching. The degree or sensitivity to which soft data conditions a realization as a function of hardness can be quantified by mapping category probabilities derived from multiple realizations. In addition to point or borehole data, arrays of data (e.g., as derived from a depth-dependency function, probability map, or “prior realization”) can be used as soft conditioning. The tsim-s algorithm provides a theoretically sound and general framework for integrating datasets of variable location, resolution, and uncertainty into geostatistical simulation of categorical variables. A practical example shows how tsim-s is capable of generating a large-scale three-dimensional simulation including curvilinear features.

54 ENVIRONMENTAL SCIENCES↗

Axion lines from nuclear de-excitations in galactic stellar populations

We show that mono-energetic axions are produced in abundance through nuclear de-excitations in nearby galaxies such as M87, which is the central galaxy of the Virgo cluster, and the starburst galaxy M82. If the axion couples to both nucleons and photons and is ultralight, then monochromatic hard X-ray signatures are induced by the subsequent axion-to-photon conversion in the magnetic fields permeating these systems. We search for evidence of such signals using NuSTAR data, focusing specifically on the $^{57}$Fe de-excitation line at 14.4 keV, and we catalog other potentially relevant nuclear lines. We find no evidence for axions from M87 or M82 and set leading constraints on the combined axion-nucleon and axion-photon coupling at the level of $|g_{ann} \times g_{aγγ}| \lesssim 1.1 \times 10^{-22}$ GeV$^{-1}$ in the limit $m_a \lesssim 10^{-10}$ eV, at 95% confidence.

Astrophysics of Galaxies (astro-ph.GA)↗

DPC Disposal Concepts of Operations. Final Report

The U.S. Department of Energy supports an R&D program for evaluating approaches to direct disposal of commercial spent fuel in dual-purpose canisters (DPCs). The major thrusts include alternative measures for treating the possibility of internal criticality events in DPC-based waste packages after thousands of years in a repository. These measures include: 1) injectable fillers, 2) analysis of the consequences of criticality events in a repository should they occur, and 3) options for modifying fuel assemblies or baskets in DPCs at the time they are loaded. This report presents a snapshot of progress in each of these areas drawing on deliverable reports generated during FY18 through FY20. Another aspect of the R&D program is to develop concepts of operations for repositories that would permanently dispose of DPC-based waste packages, considering different generic host media (not site-specific). The idea is to examine whether the disposal of large, heavy, heat-generating waste packages is technically feasible, and to identify the engineering challenges that would arise during implementation of the different disposal concepts. Descriptions of repository features are presented for repositories in salt media, argillite (clay/shale) media, crystalline (e.g., granitic) media, and unsaturated media (considering either alluvium or hard rock). Thermal management criteria for each concept are presented in terms of the maximum waste package thermal power at emplacement, when the repository could be opened, and the duration of repository emplacement operations. The overall message of this report is that direct disposal of commercial spent fuel is technically feasible in different types of geologic host media, but that thermal management and postclosure criticality impose different constraints on each concept. Engineering challenges are recognized and discussed. Treatment of postclosure criticality is identified as an important technical question that receives the majority of attention in the R&D program.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Leading Axion-Photon Sensitivity with NuSTAR Observations of M82 and M87

Here, we perform one of the most sensitive searches to date for the existence of ultralight axions using data from the NuSTAR telescope. We search for stellar axion production in the M82 starburst galaxy and the M87 central galaxy of the Virgo cluster and then the subsequent conversion into hard x-rays in the surrounding magnetic fields. We sum over the full stellar populations in these galaxies when computing the axion luminosity, and we account for the conversion of axions to photons by using magnetic field profiles in simulated IllustrisTNG analog galaxies. We show that analyzing NuSTAR data toward these targets between roughly 30 and 70 keV shows no evidence for axions and leads to robust constraints on the axion-photon coupling at the level of |𝑔𝑎⁢𝛾⁢𝛾| ≲ 7.8 × 10 −13 GeV −1 for 𝑚 𝑎 ≲ 10 −10 eV at 95% confidence.

Axions↗