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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Quasilinear theory for inhomogeneous plasma

Here, this paper presents quasilinear theory (QLT) for a classical plasma interacting with inhomogeneous turbulence. The particle Hamiltonian is kept general; for example, relativistic, electromagnetic and gravitational effects are subsumed. A Fokker–Planck equation for the dressed ‘oscillation-centre’ distribution is derived from the Klimontovich equation and captures quasilinear diffusion, interaction with the background fields and ponderomotive effects simultaneously. The local diffusion coefficient is manifestly positive-semidefinite. Waves are allowed to be off-shell (i.e. not constrained by a dispersion relation), and a collision integral of the Balescu–Lenard type emerges in a form that is not restricted to any particular Hamiltonian. This operator conserves particles, momentum and energy, and it also satisfies the H -theorem, as usual. As a spin-off, a general expression for the spectrum of microscopic fluctuations is derived. For on-shell waves, which satisfy a quasilinear wave-kinetic equation, the theory conserves the momentum and energy of the wave–plasma system. The action of non-resonant waves is also conserved, unlike in the standard version of QLT. Dewar's oscillation-centre QLT of electrostatic turbulence is proven formally as a particular case and given a concise formulation. Also discussed as examples are relativistic electromagnetic and gravitational interactions, and QLT for gravitational waves is proposed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ultra-flat and long-lived plasmons in a strongly correlated oxide

Plasmons in strongly correlated systems are attracting considerable attention due to their unconventional behavior caused by electronic correlation effects. Recently, flat plasmons with nearly dispersionless frequency-wave vector relations have drawn significant interest because of their intriguing physical origin and promising applications. However, these flat plasmons exist primarily in low-dimensional materials with limited wave vector magnitudes (q< ~0.7 Å –1 ). Here, we show that long-lived flat plasmons can propagate up to ~1.2 Å –1 in α-Ti 2 O 3 , a strongly correlated three-dimensional Mott-insulator, with an ultra-small energy fluctuation (<40 meV). The strong correlation effect renormalizes the electronic bands near Fermi level with a small bandwidth, which is responsible for the flat plasmons in α-Ti 2 O 3 . Moreover, these flat plasmons are not affected by Landau damping over a wide range of wave vectors (q < ~1.2 Å –1 ) due to symmetry constrains on the electron wavefunctions. Our work provides a strategy for exploring flat plasmons in strongly correlated systems, which in turn may give rise to novel plasmonic devices in which flat and long-lived plasmons are desirable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Neural network reconstruction of the dense matter equation of state from neutron star observables

The Equation of State (EoS) of strongly interacting cold and hot ultra-dense QCD matter remains a major challenge in the field of nuclear astrophysics. With the advancements in measurements of neutron star masses, radii, and tidal deformabilities, from electromagnetic and gravitational wave observations, neutron stars play an important role in constraining the ultra-dense QCD matter EoS. Here, in this work, we present a novel method that exploits deep learning techniques to reconstruct the neutron star EoS from mass-radius (M-R) observations. We employ neural networks (NNs) to represent the EoS in a model-independent way, within the range of ~1-7 times the nuclear saturation density. The unsupervised Automatic Differentiation (AD) framework is implemented to optimize the EoS, so as to yield through TOV equations, an M-R curve that best fits the observations. We demonstrate that this method works by rebuilding the EoS on mock data, i.e., mass-radius pairs derived from a randomly generated polytropic EoS. The reconstructed EoS fits the mock data with reasonable accuracy, using just 11 mock M-R pairs observations, close to the current number of actual observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Characterizing a supernova’s standing accretion shock instability with neutrinos and gravitational waves

Here, we perform a novel multi-messenger analysis for the identification and parameter estimation of the Standing Accretion Shock Instability (SASI) in a core collapse supernova with neutrino and gravitational wave (GW) signals. In the neutrino channel, this method performs a likelihood ratio test for the presence of SASI in the frequency domain. For gravitational wave signals we process an event with a modified constrained likelihood method. Using simulated supernova signals, the properties of the Hyper-Kamiokande neutrino detector, and O3 LIGO Interferometric data, we produce the two-dimensional probability density distribution (PDF) of the SASI activity indicator and calculate the probability of detection P D as well as the false identification probability P FI . We discuss the probability to establish the presence of the SASI as a function of the source distance in each observational channel, as well as jointly. Compared to a single-messenger approach, the joint analysis results in P D (at P FI = 0.1) of SASI activities that is larger by up to ≈ 40% for a distance to the supernova of 5 kpc. We also discuss how accurately the frequency and duration of the SASI activity can be estimated in each channel separately. Our methodology is suitable for implementation in a realistic data analysis and a multi-messenger setting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Bayesian Inference of Hybrid Star Properties from Future High-precision Measurements of Their Radii

Future high-precision X-ray and gravitational-wave observations of neutron stars (NSs) are expected to constrain NS radii with uncertainties as small as σ ≃ 0.1 km. Such unprecedented precision offers a unique opportunity to extract new information about the nature and equation of state (EOS) of supradense matter in NS cores. Using mock radius data with uncertainties ranging from σ = 1.0 to 0.1 km, together with a flexible meta-model NS EOS that allows for a first-order hadron–quark phase transition, we perform a Bayesian statistical analysis to assess the impact of radius measurements on EOS constraints. We find that high-precision radius measurements, particularly for massive NSs, significantly tighten constraints on the hadron–quark transition density ρ t , the quark matter mass fraction in NS cores, and several parameters characterizing the EOS of supranuclear hadronic matter, although the degree of improvement depends on the assumed prior range of ρ t . In contrast, even with the highest precision considered, NS radii—including those of massive stars—remain largely insensitive to the stiffness of quark matter, independent of the measurement accuracy or the prior range adopted for ρ t .

Neutron stars↗
DOE OSTI2023

Constraint of pionless EFT using two-nucleon spectra from lattice QCD

Finite-volume pionless effective field theory (FVEFT π/ ) at next-to-leading order (NLO) is used to analyze the two-nucleon lattice QCD spectrum of Ref. [1], performed at quark masses corresponding to a pion mass of approximately 800 MeV. Specifically, the effective theory is formulated in finite volume, and variational sets of wave functions are optimized using differential programming. Using these wave functions projected to the appropriate finite-volume symmetry group, variational bounds from FVEFT π/ are obtained for the ground state, as well as excited states. By comparison with the lattice QCD GEVP spectrum, different low energy constants (LECs) are constrained. Furthermore, relativistic corrections are incorporated, allowing for the s-d-wave mixing term in the deuteron channel.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Isolating the Source Region of Infrasound Travel Time Variability Using Acoustic Sensors on High-Altitude Balloons

High-altitude balloons carrying infrasound sensor payloads can be leveraged toward monitoring efforts to provide some advantages over other sensing modalities. On 10 July 2020, three sets of controlled surface explosions generated infrasound waves detected by a high-altitude floating sensor. One of the signal arrivals, detected when the balloon was in the acoustic shadow zone, could not be predicted via propagation modeling using a model atmosphere. Considering that the balloon’s horizontal motion showed direct evidence of gravity waves, we examined their role in infrasound propagation. Implementation of gravity wave perturbations to the wind field explained the signal detection and aided in correctly predicting infrasound travel times. Our results show that the impact of gravity waves is negligible below 20 km altitude; however, their effect is important above that height. The results presented here demonstrate the utility of balloon-borne acoustic sensing toward constraining the source region of variability, as well as the relevance of complexities surrounding infrasound wave propagation at short ranges for elevated sensing platforms.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

A search for transit timing variations in the HATS-18 planetary system

ABSTRACT HATS-18 b is a transiting planet with a large mass and a short orbital period, and is one of the best candidates for the detection of orbital decay induced by tidal effects. We present extensive photometry of HATS-18 from which we measure 27 times of mid-transit. Two further transit times were measured from data from the Transiting Exoplanet Survey Satellite (TESS) and three more taken from the literature. The transit timings were fitted with linear and quadratic ephemerides and an upper limit on orbital decay was determined. This corresponds to a lower limit on the modified stellar tidal quality factor of $Q_\star ^{\, \prime } \gt 10^{5.11 \pm 0.04}$. This is at the cusp of constraining the presence of enhanced tidal dissipation due to internal gravity waves. We also refine the measured physical properties of the HATS-18 system, place upper limits on the masses of third bodies, and compare the relative performance of TESS and the 1.54 m Danish Telescope in measuring transit times for this system.

Southworth, John (ORCID:0000000238073198)↗
DOE OSTI2023

Laser Characterization of Optically Smooth BREAD parts

The Broadband Reflector Experiment for Axion Detection is a dish antenna experiment probing the axion and dark photon parameter space in the [1 μeV, 1 eV] mass range, or [.02, 200] THz frequency range. In the pilot design for the InfraBREAD model, the dish antenna includes an optically smooth coaxial parabolic reflector to focus incoming photons in the near Infrared spectrum, created by the spontaneous photoconversion of axions and dark photons, onto a Superconducting Nanowire Single Photon Detector at the focus point of the reflector. As the pilot begins as a dark photon detector, sensitivity estimates hypothesize that it could cover unprobed parameter space in the search for the dark photon. Because of this, it is of critical importance that the reflector and SNSPD be placed at the correct points in the dish antenna, as previous work has shown that in addition to the focal spot of the photons being smeared to an area larger than the SNSPD, detection efficiency is s ensitive to shifts on the millimeter scale. In order to constrain this barrier to efficiency, the reflector had its focal spot, roughness, and length of surface waves measured by lasers in an optical bench setup. Results show the focal point diverging slightly from its theorized position in the antenna, while uncertainty in the height at which the detector should be placed was constrained to 10s of microns, the uncertainty in radial position could not be constrained to smaller than the detector given the methods used. Results for the roughness and surface wavelength show that the reflector was uniformly machined in a diamond-turning process used to fabricate optically smooth parts, as well as support a roughly 90% detection efficiency, reinforcing the sensitivity estimates presented.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Arctic Deployment of a Fully Integrated Self-Powered Drifting Buoy Harvesting Wave Energy via a Triboelectric Nanogenerator

The Arctic Ocean remains one of the most poorly sampled regions on Earth, where improved in situ environmental monitoring is vital for advancing oceanographic and atmospheric studies. However, data collection efforts are constrained by the short operational lifespans and high costs of conventional systems. Drifting buoys powered by pendulum-driven wave energy harvesters offer a cost-effective alternative, yet earlier designs have neither been optimized for real-world wave conditions nor validated in the Arctic. In this study, we develop a self-powered drifting buoy that integrates a pendulum-driven triboelectric nanogenerator (TENG) system with a mechanical motion rectifier, a high-gear-ratio transmission, and power management circuits. Through coupled buoy–pendulum dynamic simulations and laboratory testing using a motion simulator, we identify an optimal pendulum mass of 1.6 kg (12.7% of total buoy weight) that maximizes energy output while maintaining buoy stability. Laboratory experiments achieved average power outputs of 12.7 mW under Arctic-like wave and temperature conditions. The system was successfully deployed in the Bering Sea, where it generated 11 J of energy in 3.1 m waves, marking the first Arctic deployment of a TENG-based drifting buoy for sea surface temperature monitoring. This work establishes a cost-effective framework for designing self-powered Arctic monitoring platforms and advances the feasibility of long-term environmental observations in real Arctic waters.

marine enerby↗
DOE OSTI2025

Shear-Wave Radiation Patterns from Explosive and Earthquake Sources in Scattering, Heterogeneous Media

Distinguishing whether a seismic event is an earthquake or an explosion is a core problem in explosion monitoring. For simplistic models, earthquakes produce a predictable S-wavefield based on their radiation pattern, while an explosion produces no S-wavefield. However, observations from nuclear tests show that explosions can produce significant S-wave energy on both horizontal components. We perform numerical experiments using SPECFEM to constrain when an S-wavefield generated from an explosion might differ from one produced by an earthquake. We generate the S-wavefield for our explosions by placing the source location within a small region consisting of large velocity heterogeneities. For the 2D cases investigated here, a crucial condition is that the size of the region of heterogeneities must be much smaller (approximately 1/12th) than the minimum wavelength of the wavefield for our simulation setups. We investigate both isotropic heterogeneities and anisotropic heterogeneities. Initial results demonstrate that the heterogeneous region lowers the P/S amplitude ratios of an explosion. The presence of strong anisotropic heterogeneities near the source can also produce an S-wavefield with similar amplitudes to the P-wavefield. Lastly, we show that S-wave polarization angles vary as a function of azimuth and distance for explosive sources in the presence of small-scale heterogeneities.

58 GEOSCIENCES↗
DOE OSTI2026

pyRMG: A framework for high-throughput, large-cell DFT calculations on supercomputers

Exascale computing delivers the raw power to simulate ever larger and more chemically realistic systems, but realizing this potential requires codes that can efficiently use thousands of processors. Our real-space multigrid (RMG) density functional theory (DFT) code’s grid-decomposition approach scales nearly linearly with the number of graphics processing units (GPUs), even for simulations exceeding thousands of atoms. This scalability makes RMG a compelling tool for high-throughput DFT studies of materials that would otherwise be bottlenecked in other codes (for example, by global fast Fourier transforms in plane-wave DFT). However, the limited workflow infrastructure for RMG has thus far constrained its adoption to a small user community. In this work, we present pyRMG, a Python package designed to streamline the setup and execution of RMG DFT calculations. Built on the pymatgen and ASE (Atomic Simulation Environment) computational materials science Python packages, pyRMG automates input generation and convergence checking, and it integrates with modern job schedulers (e.g., Flux) on leadership-class platforms such as Frontier and Perlmutter. Here, we demonstrate pyRMG for a high-throughput study of strain effects in 2D 2L-Bi 2 Se 3 /2L-NbSe 2 heterostructures, which offers chemical insights into this system and shows that RMG-based workflows can converge with limited user intervention.

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

Strong Interaction Dynamics and Fermi $β$ Decay in the Nucleon and the Nucleus

Nuclear super-allowed β decay has been used to obtain tight limits on the value of the CKM matrix element V ud that is important for unitarity tests and, therefore, for tests of the standard model. Current requirements on precision are so intense that effects formerly thought too small to matter are now considered relevant. This article is a brief review of personal efforts to include the effects of strong interactions on Fermi β decay. First, I examine the role of isospin violation in the decay of the neutron. The size of the necessary correction depends upon detailed strong-interaction dynamics. The isospin violating parts of the nucleon wave function, important at the low energy of β decay, can be constrained by data taken at much higher energies, via measurements, for example, of ed → e'π ± + X reactions at Jefferson Laboratory. The next point of focus is on the role of nuclear short-ranged correlations, which affect the value of the correction needed to account for isospin violation in extracting the value of V ud . The net result is that effects previously considered as irrelevant are now considered relevant for both neutron and nuclear β decay.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

A mixed-integer PDE-constrained optimization formulation for constructing electromagnetic cloaks with multiple materials

We study the design of an electromagnetic cloak from multiple materials with an additional constraint on the mass of the cloak. Our problem is an example of a topology optimization problem, and we formulate this problem as a mixed-integer partial-differential equation constrained optimization (MIPDECO) problem, where Maxwell’s equation models the propagation of the wave through the cloak and surrounding medium. We use binary variables to model the assignment of the different materials, and their relevant properties (permittivity and density). The mass constraint adds a nontrivial constraint to this problem. We propose a two-phase strategy to solve this problem. In the first phase, we solve a continuous relaxation, and then propose a new variant of the feasibility pump that exploits the structure of the PDE to obtain an initial integral solution candidate. In the second phase, we use a trust-region approach to improve this incumbent. We also consider a continuation or mesh-sequencing approach to find better solutions faster on consecutively finer meshes. We present detailed numerical results to illustrate the effectiveness of our approaches for constructing multi-material cloaks with a mass constraint.

Calculus of Variations and Optimization↗
DOE OSTI2026

Error quantification and comparison of binary neutron star gravitational waveforms from numerical relativity codes

Future gravitational wave detections of merging binary neutron star systems have the possibility to tightly constrain the equation of state of dense nuclear matter. In order to extract such constraints, gravitational waveform models need to be calibrated to accurate numerical relativity simulations of the late inspiral and merger. In this work, we take an essential step toward classifying the error and potential systematics in current generation numerical relativity simulations of merging binary neutron stars. To this end, we perform a direct comparison of two codes (fil, spec), which differ in many aspects, including the numerical methods and discretizations used and equations solved. We find that despite these different approaches, the codes are—within current numerical resolution bounds—fully consistent, and broadly comparable in cost for a given accuracy level. Our results indicate that the error in the waveforms is primarily dominated by the hydrodynamic evolution, consistent with earlier findings in the literature. Finally, we also discuss current limitations and cost estimates for numerical relativity simulations to reach the accuracies required in the era of next-generation gravitational detectors.

gravitational waves↗
DOE OSTI2025

Pressure suppresses the density wave order in kagome metal LuNb 6 ⁢Sn 6

The density waves that develop in kagome metals ScV 6 ⁢Sn 6 and LuNb 6 ⁢Sn 6 at low temperature appear to arise from underfilled atomic columns within a V-Sn or Nb-Sn scaffolding. Compressing this network with applied pressure in ScV 6 ⁢Sn 6 suppressed the structural transition temperature by constraining atomic rattling and inhibiting the shifts that define the structural modulation. We predicted that the density wave transition in LuNb 6⁢ Sn 6 at 68 K would be suppressed by pressure as well. In this Letter, we examine the pressure dependence of the density wave transition by measuring resistance vs temperature up to 2.26 GPa. We found the transition temperature is smoothly depressed and disappears around 1.9 GPa. In conclusion, this result not only addresses our prediction, but strengthens the rattling chains origin of structural instabilities in the HfFe 6⁢ Ge 6 -type kagome metals.

Charge density waves↗
DOE OSTI2023

Dark matter or regular matter in neutron stars? How to tell the difference from the coalescence of compact objects

The mirror twin Higgs model is a candidate for (strongly-interacting) complex dark matter, which mirrors SM interactions with heavier quark masses. A consequence of this model are mirror neutron stars—exotic stars made entirely of mirror matter, which are significantly smaller than neutron stars and electromagnetically dark. This makes mergers of two mirror neutron stars detectable and distinguishable in gravitational wave observations, but can we observationally distinguish between regular neutron stars and those that may contain some mirror matter? This is the question we study in this paper, focusing on two possible realizations of mirror matter coupled to standard model matter within a compact object: (i) mirror matter captured by a neutron star and (ii) mirror neutron star-neutron star coalescences. Regarding (i), here we find that (nonrotating) mirror-matter-admixed neutron stars no longer have a single mass-radius sequence, but rather exist in a two-dimensional mass-radius plane. Regarding (ii), we find that binary systems with mirror neutron stars would span a much wider range of chirp masses and completely different binary Love relations, allowing merger remnants to be very light black holes. The implications of this are that gravitational wave observations with advanced LIGO and Virgo, and x-ray observations with NICER, could detect or constrain the existence of mirror matter through searches with wider model and parameter priors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The LIGO HET Response (LIGHETR) Project to Discover and Spectroscopically Follow Optical Transients Associated with Neutron Star Mergers

The LIGO HET Response (LIGHETR) project is an enterprise to follow up optical transients (OTs) discovered as gravitational-wave merger sources by the LIGO/Virgo collaboration (LVC). Early spectroscopy has the potential to constrain crucial parameters such as the aspect angle. The LIGHETR collaboration also includes the capacity to model the spectroscopic evolution of mergers to facilitate a real-time direct comparison of models with our data. The principal facility is the Hobby–Eberly Telescope. LIGHETR uses the massively replicated VIRUS array of spectrographs to search for associated OTs and obtain early blue spectra, and in a complementary role, the low-resolution LRS2 spectrograph is used to obtain spectra of viable candidates as well as a densely sampled series of spectra of true counterparts. Once an OT is identified, the anticipated cadence of spectra would match or considerably exceed anything achieved for GW170817 = AT2017gfo for which there were no spectra in the first 12 hr and thereafter only roughly once daily. We describe special HET-specific software written to facilitate the program and attempts to determine the flux limits to undetected sources. We also describe our campaign to follow up OT candidates during the third observational campaign of the LIGO and Virgo Scientific Collaborations. We obtained VIRUS spectroscopy of candidate galaxy hosts for five LVC gravitational-wave events and LRS2 spectra of one candidate for the OT associated with S190901ap. We identified that candidate, ZTF19abvionh = AT2019pip, as a possible Wolf–Rayet star in an otherwise unrecognized nearby dwarf galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗