Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “WAKEFULNESS”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

On Rank Selection for Nonnegative Matrix Factorization

Rank selection, i.e. the choice of factorization rank, is the first step in constructing Nonnegative Matrix Factorization (NMF) models. It is a long-standing problem which is not unique to NMF, but arises in most models which attempt to decompose data into its underlying components. Since these models are often used in the unsupervised setting, the rank selection problem is further complicated by the lack of ground truth labels. In this paper, we review and empirically evaluate the most commonly used schemes for NMF rank selection.

Eswar, Srinivas [Argonne National Laboratory]↗
DOE OSTI2025

An infrasound source analysis of the OSIRIS-REx sample return capsule hypersonic re-entry

The OSIRIS-REx sample return capsule's hypersonic re-entry into the atmosphere is a rare opportunity to test a variety of sonic boom source models since the projectile dimensions are well characterized. While the as-flown flight path is unknown, the predicted flight path enables a rough approximation of the source Mach number and location. Six infrasound microphones deployed in the boom carpet along the predicted flight path recorded impulsive signals from the OSIRIS-REx re-entry. Using a suite of atmosphere profiles and the geometric acoustics approximation, we estimate locations with uncertainty estimates along the flight path from which the signals were emitted. Acoustic overpressure and signal duration predictions from Whitham's far field theory, Carlson's simplified sonic boom prediction method, and a drag-dominated hypersonic model are analyzed with uncertainty estimates from the location estimate. While the Carlson simplified sonic boom prediction method could be accurate, our preference is for the drag-dominated source model. Using this source model with an inviscid Burgers's equation solver for propagation, we obtained an excellent match to the recorded data. In conclusion, these results will help better inform future sample return capsule re-entry observation campaigns as well as contribute to a better understanding of high altitude infrasonic sources.

58 GEOSCIENCES↗
DOE OSTI2025

Randomized Algorithms for Symmetric Nonnegative Matrix Factorization

Symmetric Nonnegative Matrix Factorization (SymNMF) is a technique in data analysis and machine learning that approximates a matrix with a product of a nonnegative, low-rank matrix and it transpose. To design faster and more scalable algorithms for SymNMF we develop two randomized algorithms for its computation. The first method uses randomized matrix sketching to compute an initial low-rank approximation to the input matrix and proceeds to uses this as a low-rank input to rapidly compute a SymNMF. The second methods uses randomized leverage score sampling to approximately solve constrained least squares problems. Many successful methods for SymNMF rely on (approximately) solving sequences of constrained least squares problems. Here, we prove theoretically that leverage score sampling can approximately solve constrained least squares problems to e-accuracy. Finally we demonstrate both methods work in practice by applying them to graph clustering tasks on large real world data sets. These experiments show that our methods approximately maintain solution quality and achieve significant speed ups for both large dense and large sparse problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Communication Lower Bounds and Optimal Algorithms for Symmetric Matrix Computations

In this article, we focus on the communication costs of three symmetric matrix computations: (i) multiplying a matrix with its transpose, known as a symmetric rank-k update (SYRK) (ii) adding the result of the multiplication of a matrix with the transpose of another matrix and the transpose of that result, known as a symmetric rank-2k update (SYR2K) (iii) performing matrix multiplication with a symmetric input matrix (SYMM). All three computations appear in the Level 3 Basic Linear Algebra Subroutines (BLAS) and have wide use in applications involving symmetric matrices. We establish communication lower bounds for these kernels using sequential and distributed-memory parallel computational models, and we show that our bounds are tight by presenting communication-optimal algorithms for each setting. Our lower bound proofs rely on applying a geometric inequality for symmetric computations and analytically solving constrained nonlinear optimization problems. As a result, the symmetric matrix and its corresponding computations are accessed and performed according to a triangular block partitioning scheme in the optimal algorithms.

Al Daas, Hussam [Rutherford Appleton Laboratory, D↗
DOE OSTI2022

Tidal Turbine Test, Downeast Turbines, July 12, 2021

Downeast Turbines tested a tidal turbine prototype with novel rotor/channel system and lateral effluent discharge apparatus (LEDA), during five days of testing in the flume at Alden Lab. Three days of testing (July 12-14, 2021) focused on turbine power metrics of torque and rpm, which were low, and then two days of follow-up testing (July 27-28, 2021) focused on LEDA performance metrics of pressure differential and rates of volumetric flow, with encouraging results. Next step is to to characterize, and even optimize, configurations of the LEDA, using 3D-CFD as a helpful tool, to refine its shape and explore its limits of performance as a means of effluent discharge that augments performance of an instream turbine. An improved configuration of the LEDA will be re-combined with the rotor/channel system of the turbine prototype, and ultimately, the rotor will be re-sized (enlarged), to better match the LEDA's performance capabilities in drawing through a rate of volumetric flow. This submission is Downeast Turbines' Post Access Report for the test event. It includes the files described here (next below), and several reference links. "Downeast TEAMER-Post-Access-Report...docx" is a document file containing the report. "Appendices A, B, and C" are included in this file, and so are "Figures #1-7." "Appendix D - Test Data Workbooks.zip" is an archive file containing all post access data (raw data tables, calculating tables, and graphs), presented in fourteen Excel workbooks as described in the report. "Appendix E - Post Access Figures.zip" is an archive file containing "Figures #8-52," of the report.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Monthly averages of ED2 model simulations initialized with airborne lidar structure, Jan 1981-Dec 2018, Brazilian Amazon

Deforestation and forest degradation (selective logging, fires, fragmentation) have impacted nearly 40% of the original extent of the Brazilian Amazon, and have markedly impacted forest structure across the region. To date, few studies analysed how shifts in forest structure from degradation influence the forest sensitivity to climate extremes, because of the complex interactions between forest structure and micro-environmental conditions. To address this knowledge gap, we carried out a series of simulations across the Brazilian Amazon using the Ecosystem Demography Model (ED2), using observed forest structure derived from 541 airborne lidar transects (375 ha each) and two scenarios representing forest recovery and expansion of degradation to investigate how shifts in forest structure impact ecosystem function under near-average and extreme climate conditions, as part of the manuscript Longo et al 2025 "Degradation and Deforestation Increase the Sensitivity of the Amazon Forest to Climate Extremes". This dataset provides the output results from the ED2 model simulations for the three simulations at monthly time scales, in NetCDF format. For all simulations, we used bias-corrected hourly reanalyses (WFDE5) for most meteorological drivers, except for precipitation, which was obtained from CHIRPS. The meteorological drivers used in the study span 38 years (Jan 1981–Dec 2018). The output results correspond to the last 38 years of simulation (one full cycle of meteorological drivers), in which ED2 simulations used static stand structure (i.e., the forest structure was held constant). The following files are provided:ED2_emean_Global_R004_BrAmaz_s1c0t0l0f0.nc. This corresponds to the Control simulation. The forest structure was obtained from the airborne lidar.ED2_emean_Global_R005_BrAmaz_s1c0t1l1f0.nc. This corresponds to the Degraded simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that expanded deforestation and selective logging across the Amazon.ED2_emean_Global_R006_BrAmaz_s1c0t1l0f0.nc. This corresponds to the Recovery simulation. The forest structure was obtained from a spin-up simulation initialized with airborne lidar and a scenario that completely halted deforestation and degradation, allowing degraded forests to recover for 38 years.We also provide file ED2_zones_R004_BrAmaz_s1c0t0l0f0.nc, which classifies each grid cell into zones used in the reference manuscript: 1: Southeast. 2: South. 3: West. 4: Central. 5: Northeast. 6: North. 7: Northwest". Index 0 corresponds to grid cells excluded from sub-region analyses because they were dominated by flooded forests, deforestation, and naturally non-forest vegetation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Southern Company Services: Project 2: Tall Tower Wind Technical Potential for the Southeast (Modifications 4, 5), Project 4: Offshore Wind Potential for Southern Company Services (Modifications 7, 8)

Project 2: The National Laboratory of the Rockies (NLR) and Southern Company Services (SCS) seek to establish an agreement to assess the technical potential for wind energy in the Southeastern U.S. SCS has requested that NLR conduct a detailed assessment of system performance, available capacity, transmission distance, and levelized cost of electricity (LCOE) at turbine hub heights of 100m to 160m within the SCS service territory. NLR will also provide an assessment of utility photovoltaic (PV) technical potential within the SCS service territory. NLR has unique modeling capabilities to assist SCS with high spatial and temporal resolution assessment of wind and PV technical potential, considering geospatial constraints. Project 4: NLR and SCS seek to develop a spatially explicit techno-economic assessment for offshore wind energy in the Southeastern U.S. The analysis aims to characterize the regional sensitivities of offshore wind plant costs and performance, underwater cabling, potential siting constraints, and landfall spur transmission costs.

17 WIND ENERGY↗
DOE OSTI2026

Wind 2H2 Project: CRADA Final Report (Project 2)

The Xcel-NLR Wind Electrolysis R&D Project is designed to gain knowledge and insight about producing hydrogen from electricity for use as a vehicle fuel and as a "stored wind". The project will analyze producing hydrogen directly from wind as well as from the electric grid. In addition, it will look at storing hydrogen and then using it to put electricity back onto the grid.

17 WIND ENERGY↗
DOE OSTI2020

The Great Markarian 421 Flare of 2010 February: Multiwavelength Variability and Correlation Studies

We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of ~27 Crab Units above 1 TeV was measured in very-highenergy (VHE) γ-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE γ-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE γ-ray (VERITAS, MAGIC), high-energy (HE) γ-ray (Fermi-LAT), X-ray (Swift, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Mets¨ahovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare ‘decline’ epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of ~25–55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor (δ & 33) and the size of the emission region (δ -1RB . 3.8 × 10 13 cm) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Multiwavelength Observations of the Blazar VER J0521+211 during an Elevated TeV Gamma-Ray State

We report on a long-lasting, elevated gamma-ray flux state from VER J0521+211 observed by VERITAS, MAGIC, and Fermi-LAT in 2013 and 2014. The peak integral flux above 200 GeV measured with the nightly binned light curve is (8.8 ± 0.4) × 10 –7 photons m –2 s –1 , or ~37% of the Crab Nebula flux. Multiwavelength observations from X-ray, UV, and optical instruments are also presented. A moderate correlation between the X-ray and TeV gamma-ray fluxes was observed, and the X-ray spectrum appeared harder when the flux was higher. Using the gamma-ray spectrum and four models of the extragalactic background light (EBL), a conservative 95% confidence upper limit on the redshift of the source was found to be z ≤ 0.31. Unlike the gamma-ray and X-ray bands, the optical flux did not increase significantly during the studied period compared to the archival low-state flux. The spectral variability from optical to X-ray bands suggests that the synchrotron peak of the spectral energy distribution (SED) may become broader during flaring states, which can be adequately described with a one-zone synchrotron self-Compton model varying the high-energy end of the underlying particle spectrum. The synchrotron peak frequency of the SED and the radio morphology of the jet from the MOJAVE program are consistent with the source being an intermediate-frequency-peaked BL Lac object.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

SDSS-IV MaNGA: How the Stellar Populations of Passive Central Galaxies Depend on Stellar and Halo Mass

We analyze spatially resolved and co-added SDSS-IV MaNGA spectra with signal-to-noise ratio ~100 from 2200 passive central galaxies (z ~ 0.05) to understand how central galaxy assembly depends on stellar mass (M*) and halo mass (M h ). We control for systematic errors in M h by employing a new group catalog from Tinker and the widely used Yang et al. catalog. At fixed M*, the strengths of several stellar absorption features vary systematically with M h . Completely model-free, this is one of the first indications that the stellar populations of centrals with identical M* are affected by the properties of their host halos. To interpret these variations, we applied full spectral fitting with the code alf. At fixed M*, centrals in more massive halos are older, show lower [Fe/H], and have higher [Mg/Fe] with 3.5σ confidence. We conclude that halos not only dictate how much M* galaxies assemble but also modulate their chemical enrichment histories. Turning to our analysis at fixed M h , high-M* centrals are older, show lower [Fe/H], and have higher [Mg/Fe] for M h > 10 12 h –1 M⊙ with confidence >4σ. While massive passive galaxies are thought to form early and rapidly, our results are among the first to distinguish these trends at fixed M h . They suggest that high-M* centrals experienced unique early formation histories, either through enhanced collapse and gas fueling or because their halos were early forming and highly concentrated, a possible signal of galaxy assembly bias.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Baryonic Content of Galaxies Mapped by MaNGA and the Gas Around Them

We analyze the cool gas in and around 14 nearby galaxies (at z < 0.1) mapped with the Sloan Digital Sky Survey IV MaNGA survey by measuring absorption lines produced by gas in spectra of background quasars/active galactic nuclei at impact parameters of 0–25 effective radii from the galactic centers. Using Hubble Space Telescope/Cosmic Origins Spectrograph, we detect absorption at the galactic redshift and measure or constrain column densities of neutral (H I, N I, O I, and Ar I), low-ionization (Si II, S II, C II, N II, and Fe II), and high-ionization (Si III, Fe III, N v, and O vI) species for 11 galaxies. We derive the ionization parameter and ionization-corrected metallicity using cloudy photoionization models. The H i column density ranges from ~1013 to ~1020 cm-2 and decreases with impact parameter for r ≳ Re. Galaxies with higher stellar mass have weaker H I absorption. Comparing absorption velocities with MaNGA radial velocity maps of ionized gas line emissions in galactic disks, we find that the neutral gas seen in absorption corotates with the disk out to ~10 R e . Sight lines with lower elevation angles show lower metallicities, consistent with the metallicity gradient in the disk derived from MaNGA maps. Higher-elevation angle sight lines show higher ionization, lower H I column density, supersolar metallicity, and velocities consistent with the direction of galactic outflow. Our data offer the first detailed comparisons of circumgalactic medium (CGM) properties (kinematics and metallicity) with extrapolations of detailed galaxy maps from integral field spectroscopy; similar studies for larger samples are needed to more fully understand how galaxies interact with their CGM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Jipole: A Differentiable ipole-based Code for Radiative Transfer in Curved Spacetimes

Recent imaging of supermassive black holes by the Event Horizon Telescope has relied on exhaustive parameter-space searches, matching observations to large, precomputed libraries of theoretical models. As observational data become increasingly precise, the limitations of this computationally expensive approach grow more acute, creating a pressing need for more efficient methods. In this work, we present Jipole, an automatically differentiable (AD), ipole-based code for radiative transfer in curved spacetimes, designed to compute image gradients with respect to underlying model parameters. These gradients quantify how parameter changes—such as the black hole’s spin or the observer’s inclination—affect the image, enabling more efficient parameter estimation and reducing the number of required images. We validate Jipole against ipole in two analytical tests and then compare pixelwise intensity derivatives from AD with those from finite-difference methods. We then demonstrate the utility of these gradients by performing parameter recovery for an analytical model in three increasingly complex cases for the injected image: ideal, blurred, and blurred with added noise. In most cases, high-accuracy fits are obtained in only a few optimization steps, failing only in cases with extremely low signal-to-noise ratios. These results highlight the potential of AD-based methods to accelerate robust, high-fidelity model-data comparisons in current and future black hole imaging efforts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Optimizing Optical Searches for Supermassive Black Hole Binaries in Active Galactic Nuclei Light Curves: Fourier versus Bayesian Periodicity Detection

Simulations predict that supermassive black hole binaries (SMBHBs) will exhibit periodic brightness variations that may exceed the stochastic variability intrinsic to active galactic nuclei (AGN). In this paper, we simulate SMBHBs with damped random walk (DRW) AGN variability and an added sinusoidal signal from the orbital motion, and test three methods—a generalized Lomb–Scargle periodogram (GLSP), a nested Bayesian sampler (NBS), and a weighted wavelet z-transform (or WWZ)—to determine which is best at recovering the periodicity. Our simulated light curves follow the properties of the Catalina Real-Time Transient Survey (or CRTS), Legacy Survey of Space and Time (LSST), and Zwicky Transient Facility (ZTF) to best inform current and future SMBHB searches. We map a broad range of parameter space and identify which DRW-only light curves best mimic periodicity and pass each method’s model selection. The NBS performs best at detecting periodicity and filtering out DRW-only light curves. Combined candidate selection with both the NBS and GLSP significantly reduces false-positive rates (FPRs) with marginal impact on true-positive rates (TPRs). With this joint model selection pipeline, we find the lowest FPRs in ZTF-like simulations and the highest detection rates in LSST-like simulations. Using a modified computation of the false-alarm probability with GLSP, we efficiently triage LSST AGN light curves (∼10 7 light curves in ∼10–30 hr) and achieve TPRs and FPRs of ∼40% and ∼0.5%, respectively.

Banaszak, Sebastian M. [Vanderbilt Univ., Nashvill↗