Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained wave”

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 235 records · Page 13

Impact of dark matter on tidal signatures in neutron star mergers with the Einstein Telescope

If dark matter (DM) accumulates inside neutron stars (NS), it changes their internal structure and causes a shift of the tidal deformability from the value predicted by the dense-matter equation of state (EOS). In principle, this shift could be observable in the gravitational-wave (GW) signal of binary neutron star (BNS) mergers. We investigate the effect of fermionic, noninteracting DM when observing a large number of GW events from DM-admixed BNSs with the precision of the proposed Einstein telescope (ET). Specifically, we study the impact on the recovery of the baryonic EOS and whether DM properties can be constrained. For this purpose, we create event catalogs of BNS mock events with DM fraction up to 1%, from which we reconstruct the posterior uncertainties with the Fisher matrix approach. Using this data, we perform joint Bayesian inference on the baryonic EOS, DM particle mass, and DM particle fraction in each event. Here, our results reveal that when falsely ignoring DM effects, the EOS posterior is biased toward softer EOSs, though the offset is rather small. Further, we find that within our assumptions of our DM model and population, ET will likely not be able to test the presence of DM in BNSs, even when combining many events and adding Cosmic Explorer (CE) to the next-generation detector network. Likewise, the potential constraints on the DM particle mass will remain weak because of degeneracies with the fraction and EOS.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Experimental shock Hugoniot and melt boundary temperatures of aluminum

Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.

Hartsfield, Thomas Murray [Sandia National Laborat↗

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗

Quantum Annealing with Inequality Constraints: The Set Cover Problem

Abstract Quantum annealing is a promising method for solving hard optimization problems by transforming them into quadratic unconstrained binary optimization (QUBO) problems. However, when constraints are involved, particularly multiple inequality constraints, incorporating them into the objective function poses challenges. In this paper, the authors present two novel approaches for solving problems with multiple inequality constraints on a quantum annealer and apply them to the set cover problem (SCP). The first approach uses the augmented Lagrangian method to represent the constraints, while the second approach employs a higher‐order binary optimization (HUBO) formulation. The experiments show that both approaches outperform the standard approach for solving the SCP on the D‐Wave Advantage quantum annealer. The HUBO formulation performs slightly better than the augmented Lagrangian method in solving the SCP, but its scalability in terms of embeddability in the quantum chip is worse. The results demonstrate that the proposed augmented Lagrangian and HUBO methods can successfully implement a large number of inequality constraints, making them applicable to a broad range of constrained problems beyond the SCP.

Djidjev, Hristo N.↗

On the Scaling and Growth Limit of Fluvial Dunes

Abstract Bedforms emerge in a variety of shapes and sizes when a granular erodible surface is subject to a strong enough shear flow, as observed in topographic data from submarine canyons, rivers, deserts, and planetary bodies. The two salient features of bedforms are the ability to collectively transport particles, and to generate form drag thus increasing flow resistance. These two mechanisms are in competition and contribute to force bedforms of increasing size migrating more slowly. In a dedicated large‐scale open channel flow facility, hierarchies of fluvial bedforms were generated and measured in equilibrium conditions. The corresponding scale‐dependent migration velocity and mass flux contributions were quantified in the wave number and frequency domains. Experimental results are paired with a validated set of theoretical models to demonstrate that ripples or dunes reach an equilibrium state when drag partitioning ensures high enough frictional drag to sustain the bedload transport of sediment, and low enough form drag to enable the migration of the largest bedform size. This mechanism is inferred to constrain the growth of bedforms when sediment supply is not the limiting factor.

58 GEOSCIENCES↗

A finite-volume scheme for modeling compressible magnetohydrodynamic flows at low Mach numbers in stellar interiors

Fully compressible magnetohydrodynamic (MHD) simulations are a fundamental tool for investigating the role of dynamo amplification in the generation of magnetic fields in deep convective layers of stars. The flows that arise in such environments are characterized by low (sonic) Mach numbers ($\mathscr{M}$ son ≲ 10 –2 ). In these regimes, conventional MHD codes typically show excessive dissipation and tend to be inefficient as the Courant–Friedrichs–Lewy (CFL) constraint on the time step becomes too strict. In this work we present a new method for efficiently simulating MHD flows at low Mach numbers in a space-dependent gravitational potential while still retaining all effects of compressibility. The proposed scheme is implemented in the finite-volume SEVEN-LEAGUE HYDRO (SLH) code, and it makes use of a low-Mach version of the five-wave Harten–Lax–van Leer discontinuities (HLLD) solver to reduce numerical dissipation, an implicit–explicit time discretization technique based on Strang splitting to overcome the overly strict CFL constraint, and a well-balancing method that dramatically reduces the magnitude of spatial discretization errors in strongly stratified setups. The solenoidal constraint on the magnetic field is enforced by using a constrained transport method on a staggered grid. We carry out five verification tests, including the simulation of a small-scale dynamo in a star-like environment at $\mathscr{M}$ son ~ 10 –3 . We demonstrate that the proposed scheme can be used to accurately simulate compressible MHD flows in regimes of low Mach numbers and strongly stratified setups even with moderately coarse grids.

79 ASTRONOMY AND ASTROPHYSICS↗

Machine learning framework for quantum sampling of highly constrained, continuous optimization problems

In recent years, there is growing interest in using quantum computers for solving combinatorial optimization problems. In this work, we developed a generic, machine learning-based framework for mapping continuous-space inverse design problems into surrogate quadratic unconstrained binary optimization (QUBO) problems by employing a binary variational autoencoder and a factorization machine. The factorization machine is trained as a low-dimensional, binary surrogate model for the continuous design space and sampled using various QUBO samplers. Using the D-Wave Advantage hybrid sampler and simulated annealing, we demonstrate that by repeated resampling and retraining of the factorization machine, our framework finds designs that exhibit figures of merit exceeding those of its training set. We showcase the framework’s performance on two inverse design problems by optimizing (i) thermal emitter topologies for thermophotovoltaic applications and (ii) diffractive meta-gratings for highly efficient beam steering. This technique can be further scaled to leverage future developments in quantum optimization to solve advanced inverse design problems for science and engineering applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The luminous, slow-rising orphan afterglow AT2019pim as a candidate moderately relativistic outflow

ABSTRACT Classical gamma-ray bursts (GRBs) have two distinct emission episodes: prompt emission from ultrarelativistic ejecta and afterglow from shocked circumstellar material. While both components are extremely luminous in known GRBs, a variety of scenarios predict the existence of luminous afterglow emission with little or no associated high-energy prompt emission. We present AT 2019pim, the first spectroscopically confirmed afterglow with no observed high-energy emission to be identified. Serendipitously discovered during follow-up observations of a gravitational-wave trigger and located in a contemporaneous TESS sector, it is hallmarked by a fast-rising ($t \approx 2$ h), luminous ($M_{\rm UV,peak} \approx -24.4$ mag) optical transient with accompanying luminous X-ray and radio emission. No gamma-ray emission consistent with the time and location of the transient was detected by Fermi-GBM or by Konus, placing constraining limits on an accompanying GRB. We investigate several independent observational aspects of the afterglow in the context of constraints on relativistic motion and find all of them are consistent with an initial Lorentz factor of $\Gamma _0 \approx$ 10–30 for the on-axis material, significantly lower than in any well-observed GRB and consistent with the theoretically predicted ‘dirty fireball’ scenario in which the high-energy prompt emission is stifled by pair production. However, we cannot rule out a structured jet model in which only the line-of-sight material was ejected at low-$\Gamma$, off-axis from a classical high-$\Gamma$ jet core, and an on-axis GRB with below-average gamma-ray efficiency also remains a possibility. This event represents a milestone in orphan afterglow searches, demonstrating that luminous optical afterglows lacking detected GRB counterparts can be identified and spectroscopically confirmed in real time.

Perley, Daniel A. (ORCID:0000000184721996)↗

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Bridging the gap between molecules and materials in quantum chemistry with localized active spaces

The number of materials that “bridge the gap” between single molecules and extended solids, such as metal-organic frameworks and organic semiconductors, has been increasing. Consequently, there is a growing need for modeling approaches that effectively integrate the real-space molecular perspective employed by computational chemists and the reciprocal-space dispersive perspective employed by computational physicists. Here, we propose the localized active space (LAS) approach as a promising method to successfully bridge this gap. The LAS approach extends the active space concept from multiconfigurational methods such as complete active space self-consistent field theory to multiple molecular fragments via a product-form wave function ansatz. Here, we apply this method to solid state phenomena by treating each unit cell as a fragment with different sets of local quantum numbers (e.g., charge and excitation number). State interaction between these LAS states (LASSI) thus provides a comprehensive basis for the study of charge and energy transfer, meeting and surpassing the capabilities of single-reference fragmentation approaches such as constrained density functional theory (cDFT). Most centrally, we show how combining this LASSI approach with multiconfigurational pair-density functional theory (MC-PDFT) provides an elegant and efficient method to compute band structures that capture multiconfigurational character. We apply the LASSI band structure approach to the computation of band gaps in stretched hydrogen chain, polyacetylene, and bulk nickel oxide (NiO), finding good or excellent quantitative agreement with reference values in all cases. Additionally, we use the LAS basis in one-dimensional model systems to demonstrate its ability to treat difficult solid-state phenomena such as exciton transfer and excitation at p-n junctions.

method development↗

pygwb: A Python-based Library for Gravitational-wave Background Searches

The collection of gravitational waves (GWs) that are either too weak or too numerous to be individually resolved is commonly referred to as the gravitational-wave background (GWB). A confident detection and model-driven characterization of such a signal will provide invaluable information about the evolution of the universe and the population of GW sources within it. We present a new, user-friendly, Python-based package for GW data analysis to search for an isotropic GWB in ground-based interferometer data. We employ cross-correlation spectra of GW detector pairs to construct an optimal estimator of the Gaussian and isotropic GWB, and Bayesian parameter estimation to constrain GWB models. The modularity and clarity of the code allow for both a shallow learning curve and flexibility in adjusting the analysis to one’s own needs. We describe the individual modules that make up pygwb, following the traditional steps of stochastic analyses carried out within the LIGO, Virgo, and KAGRA Collaboration. We then describe the built-in pipeline that combines the different modules and validate it with both mock data and real GW data from the O3 Advanced LIGO and Virgo observing run. We successfully recover all mock data injections and reproduce published results.

79 ASTRONOMY AND ASTROPHYSICS↗

Gravitational wave measurement in the mid-band with atom interferometers

Abstract Gravitational Waves (GWs) have been detected in the ∼ 100 Hz and nHz bands, but most of the gravitational spectrum remains unobserved. A variety of detector concepts have been proposed to expand the range of observable frequencies. In this work, we study the capability of GW detectors in the “mid-band”, the ∼ 30 mHz– 10 Hz range between LISA and LIGO, to measure the signals from and constrain the properties of ∼ 1 – 100M ⊙ compact binaries. We focus on atom-interferometer-based detectors. We describe a Fisher matrix code,AIMforGW, which we created to evaluate their capabilities, and present numerical results for two benchmarks: terrestrial km-scale detectors, and satellite-borne detectors in medium Earth orbit. Mid-band GW detectors are particularly well-suited to pinpointing the location of GW sources on the sky. We demonstrate that a satellite-borne detector could achieve sub-degree sky localization for any detectable source with chirp mass ℳ c ≲ 50M ⊙ . We also compare different detector configurations, including different locations of terrestrial detectors and various choices of the orbit of a satellite-borne detector. As we show, a network of only two terrestrial single-baseline detectors or one single-baseline satellite-borne detector would each provide close-to-uniform sky-coverage, with signal-to-noise ratios varying by less than a factor of two across the entire sky. We hope that this work contributes to the efforts of the GW community to assess the merits of different detector proposals.

Astronomy & Astrophysics↗

Symmetry Breaking with the SCAN Density Functional Describes Strong Correlation in the Singlet Carbon Dimer

The SCAN (strongly constrained and appropriately normed) meta-generalized gradient approximation (meta-GGA), which satisfies all 17 exact constraints that a meta-GGA can satisfy, accurately describes equilibrium bonds that are normally correlated. With symmetry breaking, it also accurately describes some sd equilibrium bonds that are strongly correlated. While sp equilibrium bonds are nearly always normally correlated, the C 2 singlet ground state is known from correlated wave function theory to be a rare case of strong correlation in an sp equilibrium bond. Earlier work that calculated atomization energies of the molecular sequence B 2 , C 2 , O 2 , and F 2 in the local spin density approximation (LSDA), the Perdew–Burke–Ernzerhof (PBE) GGA, and the SCAN meta-GGA, without symmetry breaking in the molecule, found that only SCAN was accurate enough to reveal an anomalous under-binding for C 2 . Here, this work shows that spin symmetry breaking in singlet C 2 , which involves the appearance of net up- and down-spin densities on opposite sides (not ends) of the bond, corrects that underbinding, with a small SCAN atomization-energy error more like that of the other three molecules, suggesting that symmetry breaking with an advanced density functional might reliably describe strong correlation. This article also discusses some general aspects of symmetry breaking and the insights into strong correlation that symmetry breaking can bring. The normally correlated low-lying triplet excited state has the right vertical excitation energy in SCAN but not in LSDA or PBE, where the triplet is a false ground state. Fractional occupation numbers are found only for the symmetry-unbroken singlet and only in LSDA and PBE GGA.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Searching for multimessenger transients and lensed supernovae with the dark energy survey

Advances in data collection technology for optical astronomy, such as the Dark Energy Camera, have facilitated the collection of some of the largest astronomical datasets in terms of sky area, depth, number of objects, number of observations, and frequency of observations. Consequentially, physical systems that have been historically rare become represented in small, but discoverable, numbers. I have developed techniques to identify rare physical systems (counterparts to high-energy neutrinos, counterparts to gravitational waves, and gravitationally lensed supernovae) that can shed new light on the physical processes that govern the Universe. For each system of interest, I collect new or search archival observations with the Dark Energy Camera, develop selection methodologies to find the system of interest in the observations, and quantify the efficiency of the search using simulations. Each of the three analyses has its own scientific implications: I constrain the observing program required to determine the contribution of supernovae to the diffuse high-energy neutrino flux; I set constraints on the physical properties of an optical counterpart to a compact object merger; and I identify candidate lensed supernovae using a novel machine learning technique. Looking forward, the work presented in this dissertation informs future rare object searches and facilitates measurements of the underlying properties of the Universe during the next era of optical survey astronomy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Gravitational Wave Eigenfrequencies from Neutrino-driven Core-collapse Supernovae

Abstract Core-collapse supernovae (CCSNe) are predicted to produce gravitational waves (GWs) that may be detectable by Advanced LIGO/Virgo. These GW signals carry information from the heart of these cataclysmic events, where matter reaches nuclear densities. Recent studies have shown that it may be possible to infer the properties of the proto-neutron star (PNS) via GWs generated by hydrodynamic perturbations of the PNS. However, we lack a comprehensive understanding of how these relationships may change with the properties of CCSNe. In this work, we build a self-consistent suite of over 1000 exploding CCSNe from a grid of progenitor masses and metallicities combined with six different nuclear equations of state (EOS). Performing a linear perturbation analysis on each model, we compute the resonant GW frequencies of the PNS, and we motivate a time-agnostic method for identifying characteristic frequencies of the dominant GW emission. From this, we identify two characteristic frequencies, of the early- and late-time signal, that measure the surface gravity of the cold remnant neutron star, and simultaneously constrain the hot nuclear EOS. However, we find that the details of the CCSN model, such as the treatment of gravity or the neutrino transport, and whether it explodes, noticeably change the magnitude and evolution of the PNS eigenfrequencies.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding carbon sourcing and transport originating from the helicon antenna surfaces during high-power helicon discharge in DIII-D Tokamak

The high-power helicon wave system in the DIII-D tokamak could potentially introduce new plasma–material interaction (PMI) challenges owing to rectified RF sheath potentials that develop near the antenna and surrounding plasma-facing components. We present the first application of the STRIPE (Simulated Transport of RF Impurity Production and Emission) framework to helicon-induced PMIs, extending previous STRIPE studies of ICRH antennas by incorporating net erosion, local re-deposition, and three-dimensional global impurity transport. The integrated workflow couples SOLPS, COMSOL, RustBCA, GITR, and GITRm to simulate carbon erosion, re-deposition, and impurity transport for two experimentally constrained DIII-D H-mode helicon operating scenarios with different antenna–plasma gaps, coupled RF powers, and edge plasma conditions. COMSOL predicts rectified RF sheath potentials of 1–5 kV localized near the lower portion of the antenna, where the magnetic field intersects the surface at grazing incidence. Carbon self-sputtering dominates the erosion source, whereas RF-accelerated D + ions contribute approximately 1% of the total gross erosion. The smaller-gap operating scenario exhibits substantially stronger gross erosion, enhanced local re-deposition (∼12%), and a larger confined carbon inventory owing to increased plasma accessibility and broader RF sheath coverage. Comparison with available DIII-D measurements shows no distinct change in the global carbon signal that correlates with the helicon RF pulse, consistent with the simulations indicating that the helicon-generated carbon source remains small compared with the existing background carbon inventory under the present graphite-wall operating conditions. These results demonstrate the capability of STRIPE to integrate RF sheath modeling, plasma transport, surface interaction physics, and three-dimensional impurity transport for the interpretation of helicon-induced PMIs. The study further identifies the principal sources of modeling uncertainty, including grazing-angle RF sheath physics, slow-wave resolution, plasma-background extrapolation, and trace-impurity transport assumptions, providing a framework for future validation and model development.

Kumar, Atul [Oak Ridge National Laboratory (ORNL),↗

Fast Seismic Anomalies Under Continents Explained by the Delaminated Lower Continental Crust—Implications From High Pressure‐Temperature Elasticity of Jadeite

Seismic tomography has shown that the shear wave velocities (Vs) under continents, especially under cratons, are extremely fast at 100–200 km depth, which is difficult to explain by low temperatures or high Mg#. Alternatively, delaminated eclogitic lower continental crust has been proposed to account for these fast seismic anomalies. However, the thermoelastic properties of jadeite which constitutes up to 60–80 mol% of clinopyroxene in the potentially delaminated lower continental crust are not well constrained. In this study, we measured the single-crystal elasticity of jadeite by Brillouin spectroscopy under simultaneous high pressure and temperature conditions for the first time. We found that the temperature dependence of Vs of jadeite is extremely small if not negligible. The seismic velocities of the potentially delaminated lower continental crusts were subsequently modeled and found to match the widely observed fast seismic anomalies under cratons between 100 and 200 km depth.

Geology↗

Near-canopy horizontal concentration heterogeneity of semivolatile oxygenated organic compounds and implications for 2-methyltetrols primary emissions

Semivolatile oxygenated organic compounds (SV-OVOCs) are important atmospheric species, in particular for the production and chemistry of atmospheric particulate matter and related impacts on air quality and climate. In this study, SV-OVOCs were collected in the horizontal plane of the roughness layer over the tropical forest in the central Amazon during the wet season of 2018. A sampler mounted to a copter-type, hovering unmanned aerial vehicle was used. Underlying the collection region, a plateau forest transitioned into a slope forest across several hundred meters. The concentrations of pinonic and pinic acids, which are monoterpene oxidation products, had no statistical difference over the two forests. By comparison, across the study period, differences in the concentration of 2-methyltetrols, which are products of isoprene oxidation, ranged from -70% to +480% over the two forests. The chemical lifetime of 2-methyltetrols in the atmosphere is sufficiently long that heterogeneity in the isoprene emission rate from the two forests followed by atmospheric oxidation does not explain the concentration heterogeneity of 2-methyltetrols. Standing waves and local meteorology also cannot account for the heterogeneity. Of the possibilities considered, the most plausible explanation is the direct emission from the forest of 2-methyltetrols produced through biological processes within the plants. Under this explanation, the rate of direct SV-OVOC emissions should be modulated by forest type and related environmental stressors. Direct emissions of SV-OVOCs should be more broadly considered for constraining and improving models of atmospheric composition, transport, and chemistry over tropical forests.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗