Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral methods”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A data-driven approach to sampling matrix selection for compressive sensing

Sampling is a fundamental aspect of any implementation of compressive sensing. Typically, the choice of sampling method is guided by the reconstruction basis. However, this approach can be problematic with respect to certain hardware constraints and is not responsive to domain-specific context. We propose a method for defining an order for a sampling basis that is optimal with respect to capturing variance in data, thus allowing for meaningful sensing at any desired level of compression. We focus on the Walsh-Hadamard sampling basis for its relevance to hardware constraints, but our approach applies to any sampling basis of interest. We illustrate the effectiveness of our method on the Physical Sciences Inc. Fabry-Perot interferometer sensor multispectral dataset, the Johns Hopkins Applied Physics Lab FTIR-based longwave infrared sensor hyperspectral dataset, and a Colorado State University Swiss Ranger depth image dataset. The spectral datasets consist of simulant experiments, including releases of chemicals such as GAA and SF6. We combine our sampling and reconstruction with the adaptive coherence estimator (ACE) and bulk coherence for chemical detection and we incorporate an algorithmic threshold for ACE values to determine the presence or absence of a chemical. We compare results across sampling methods in this context. We have successful chemical detection at a compression rate of 90%. For all three datasets, we compare our sampling approach to standard orderings of sampling basis such as random, sequency, and an analog of sequency that we term `frequency.' In one instance, the peak signal to noise ratio was improved by over 30% across a test set of depth images.

compressive sensing, hyperspectral imaging, Walsh-↗
DOE OSTI2025

Tailored Frequency Conversion in Nonlinear Subwavelength Grating Metaslabs

The use of metaslab waveguides consisting of subwavelength gratings (SWGs) is proposed as a highly flexible approach to obtain phase‐matched frequency conversion in nonlinear integrated devices. Control over the phase matching (PM) condition is achieved through modification of the linear modal propagation constants by the SWG. The method is experimentally demonstrated by fabrication of 0.4 mm‐long GaAs metaslabs designed for second harmonic generation of near‐infrared light. A controlled spectral shift of up to 100 nm in the PM condition, relative to a uniform slab is observed. Excellent agreement is found between the experimental results and finite‐difference eigenmode simulations, which predict a tunable frequency conversion bandwidth of ∼300 nm for the design parameters. The different regimes of PM control supported by the structure are discussed. This work highlights the potential of SWGs as building blocks for tailored frequency conversion in nonlinear integrated devices.

Form birefringence↗
DOE OSTI2024

Robust Machine Learning Inference from X-ray Absorption Near Edge Spectra through Featurization

X-ray absorption spectroscopy (XAS) is a commonly employed technique for characterizing functional materials. In particular, X-ray absorption near edge spectra (XANES) encode local coordination and electronic information, and machine learning approaches to extract this information are of significant interest. To date, most ML approaches for XANES have primarily focused on using the raw spectral intensities as input, overlooking the potential benefits of incorporating spectral transformations and dimensionality reduction techniques into ML predictions. Here, in this work, we focused on systematically comparing the impact of different featurization methods on the performance of ML models for XAS analysis. We evaluated the classification and regression capabilities of these models on computed data sets and validated their performance on previously unseen experimental data sets. Our analysis revealed an intriguing discovery: the cumulative distribution function feature achieves both high prediction accuracy and exceptional transferability. This remarkably robust performance can be attributed to its tolerance to horizontal shifts in the spectra, which is crucial when validating models using experimental data. While this work exclusively focuses on XANES analysis, we anticipate that the methodology presented here will hold promise as a versatile asset to the broader spectroscopy community.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Probing atomic physics at ultrahigh pressure using laser-driven implosions

Abstract Spectroscopic measurements of dense plasmas at billions of atmospheres provide tests to our fundamental understanding of how matter behaves at extreme conditions. Developing reliable atomic physics models at these conditions, benchmarked by experimental data, is crucial to an improved understanding of radiation transport in both stars and inertial fusion targets. However, detailed spectroscopic measurements at these conditions are rare, and traditional collisional-radiative equilibrium models, based on isolated-atom calculations and ad hoc continuum lowering models, have proved questionable at and beyond solid density. Here we report time-integrated and time-resolved x-ray spectroscopy measurements at several billion atmospheres using laser-driven implosions of Cu-doped targets. We use the imploding shell and its hot core at stagnation to probe the spectral changes of Cu-doped witness layer. These measurements indicate the necessity and viability of modeling dense plasmas with self-consistent methods like density-functional theory, which impact the accuracy of radiation transport simulations used to describe stellar evolution and the design of inertial fusion targets.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Laser activation of single group-IV colour centres in diamond

Abstract Spin-photon interfaces based on group-IV colour centres in diamond offer a promising platform for quantum networks. A key challenge in the field is realising precise single-defect positioning and activation, which is crucial for scalable device fabrication. Here we address this problem by demonstrating a two-step fabrication method for tin vacancy (SnV − ) centres that uses site-controlled ion implantation followed by local femtosecond laser annealing with in-situ spectral monitoring. The ion implantation is performed with sub-50 nm resolution and a dosage that is controlled from hundreds of ions down to single ions per site, limited by Poissonian statistics. Using this approach, we successfully demonstrate site-selective creation and modification of single SnV − centres. Our in-situ spectral monitoring opens a window onto materials tuning at the single defect level, and provides new insight into defect structures and dynamics during the annealing process. While demonstrated for SnV − centres, this versatile approach can be readily generalised to other implanted colour centres in diamond and wide-bandgap materials.

Science & Technology - Other Topics↗
DOE OSTI2025

Multiple Carrier Generation at an Exceptionally Low Energy Threshold

Multiple carrier generation (MCG), a process wherein two or more carriers are generated from a single high-energy absorbed photon, holds immense promise for quantum sensing, metrology, low-threshold lasers, and photovoltaics. Despite its potential, MCG has faced obstacles such as low efficiency and a high threshold photon energy at least twice the band gap (2?Eg) of the semiconductor, limiting its application only to a class of materials with low Eg. Here, we present a new approach that overcomes this limitation by leveraging carrier-donor scattering to excite secondary electrons from donor states strategically positioned below the conduction band. Our method relies on strong Coulomb interaction, reduced dielectric screening, slow hot carrier cooling, and strictly follows the energy conservation rules. We experimentally demonstrated this idea in a model system of monolayer (1L) MoS2 by exploiting electron-donating chalcogen vacancy states. We observed an exceptionally low MCG threshold of ~1.12?Eg for the first time in 1L MoS2. Remarkably, the quantum yield can be further increased to >3 by increasing the photon energy to 1.65?Eg, representing a substantial advancement over existing methods. Our findings extend the horizon of MCG into next-generation high-performance optoelectronic devices with an on-demand operating spectral range spanning from infrared to ultraviolet.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS↗
DOE OSTI2021

Fast multiscale contrast independent preconditioners for linear elastic topology optimization problems

The goal of this work is to present a fast and viable approach for the numerical solution of the high-contrast state problems arising in topology optimization. The optimization process is iterative, and the gradients are obtained by an adjoint analysis, which requires the numerical solution of large high-contrast linear elastic problems with features spanning several length scales. The size of the discretized problems forces the utilization of iterative linear solvers with solution time dependent on the quality of the preconditioner. The lack of clear separation between the scales, as well as the high-contrast, imposes severe challenges on the standard preconditioning techniques. Thus, here we propose new methods for the high-contrast elasticity equation with performance independent of the high-contrast and the multi-scale structure of the elasticity problem. The solvers are based on two-levels domain decomposition techniques with a carefully constructed coarse level to deal with the high-contrast and multi-scale nature of the problem. The construction utilizes spectral equivalence between scalar diffusion and each displacement block of the elasticity problems and, in contrast to previous solutions proposed in the literature, is able to select the appropriate dimension of the coarse space automatically. The new methods inherit the advantages of domain decomposition techniques, such as easy parallelization and scalability. Finally, the presented numerical experiments demonstrate the excellent performance of the proposed methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Improved Gas Plume Identification Using Nearest Neighbor Methods for Background Estimation

Longwave infrared (LWIR) hyperspectral imaging (HSI) can be used for many tasks in remote sensing, including detecting and identifying effluent gases by LWIR sensors on airborne platforms. Identification is used after detection to increase confidence in weakly detected plumes, reduce false positives from detection, and distinguish between similar and confounding material signatures. Background estimation is an important step used to reveal the unique spectral characteristics of the detected gas, allowing the identification model to determine what the gas is specifically. The importance of proper background estimation increases when dealing with weak signals, large libraries of gases of interest, and uncommon or heterogeneous backgrounds. In this article, we propose two methods for background estimation: a novel k-nearest segments (KNS) algorithm and the standard k-nearest neighbors (KNN) algorithm. We test our methods and three existing background estimation methods for comparison against global background estimation to determine which performs best at estimating the true background radiance under a plume and for increasing identification confidence using a neural network classification model. We compare the different methods using 640 simulated weak plumes in an urban environment. For identification, our KNS algorithm improves median neural network identification confidence by 53.2%. For background radiance estimation, the KNN algorithm provides a median of 49 times less RMSE than global background estimation. Furthermore, KNN is the easiest method to tune for different plumes, making it an excellent “out of the box” background estimator.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

End-to-end GPU acceleration of low-order-refined preconditioning for high-order finite element discretizations

In this article, we present algorithms and implementations for the end-to-end GPU acceleration of matrix-free low-order-refined preconditioning of high-order finite element problems. The methods described here allow for the construction of effective preconditioners for high-order problems with optimal memory usage and computational complexity. The preconditioners are based on the construction of a spectrally equivalent low-order discretization on a refined mesh, which is then amenable to, for example, algebraic multigrid preconditioning. The constants of equivalence are independent of mesh size and polynomial degree. For vector finite element problems in H(curl) and H(div) (e.g., for electromagnetic or radiation diffusion problems), a specially constructed interpolation–histopolation basis is used to ensure fast convergence. Detailed performance studies are carried out to analyze the efficiency of the GPU algorithms. The kernel throughput of each of the main algorithmic components is measured, and the strong and weak parallel scalability of the methods is demonstrated. The different relative weighting and significance of the algorithmic components on GPUs and CPUs is discussed. Results on problems involving adaptively refined nonconforming meshes are shown, and the use of the preconditioners on a large-scale magnetic diffusion problem using all spaces of the finite element de Rham complex is illustrated.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A Dirac R-matrix calculation for the electron-impact excitation of W +

Abstract Aims : tungsten has been chosen for use as a plasma facing component in the divertor for the ITER experiment, and is currently being used on existing tokamaks such as JET. W + plays an integral role in assessing the impurity influx from plasma facing component of tokamaks and subsequent redeposition. Together with previously calculated a neutral tungsten electron-impact dataset this study allows us to determine neighbouring spectral lines in the same wavelength window of the spectrometer, and detect if there is strong blending of overlapping lines between these two ion stages as well as providing ionisation per photon ratios for both species. The new data is to be used for tungsten erosion/redeposition diagnostics. Methods : a significantly modified version of the GRASP0 atomic structure code in conjunction with DARC (Dirac Atomic R -matrix Code) are used to calculate the Einstein A coefficients and collisional rates used to generate a synthetic W II spectrum. The W II spectrum is compared against tungsten spectral emission experiments. Results : this study is used to model the spectrum of W ii , providing the predictive capability of identifying spectral lines from recent experiments. These results provide an integral part of impurity influx and redeposition determination, as the ionisation rates may be used to calculate S/XB ratios.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Calibration-free analysis of Li isotope ratios using laser ablation and laser absorption spectroscopy

We introduce a rapid, calibration-free, all-optical method for high-precision lithium isotope ratio measurements in solid materials using laser ablation combined with tunable laser absorption spectroscopy. A new asynchronous acquisition method is used to acquire time-resolved, high-resolution spectra of the 6Li and 7Li D1 and D2 transitions near 671 nm. Isotope ratios and atomic column densities are extracted from measured spectra via a physics-based fitting model including hyperfine structure. Under 1 Torr air, spectra recorded = 0.75 ms after plasma onset exhibit narrow linewidths corresponding to Doppler temperatures = 400 K, enabling resolution of the isotope peaks with high signal-to-noise ratios. Analysis of LiAlO2 samples with varying 6Li:7Li ratios demonstrates isotopic precisions of 0.6–1.7% for spectra acquired in 30 s. Isotope ratios determined from the spectral fits show accuracy within –0.3% to –1.7% of reference ICP-MS measurements without requiring calibration to external standards. By eliminating sample preparation and enabling spatially resolved isotopic mapping, this method offers a rapid analysis approach to lithium isotope determination in solid materials relevant to nuclear energy, safeguards, and geochemistry.

Phillips, Mark C.↗
DOE OSTI2025

Neural chaos: A spectral stochastic neural operator

Building surrogate models for operators with uncertainty quantification capabilities is essential for many engineering applications where randomness–such as variability in material properties, boundary conditions, and initial conditions–is unavoidable. Polynomial Chaos Expansion (PCE) is widely recognized as a go-to method for constructing stochastic surrogates in both intrusive and non-intrusive ways, and it has recently been used in the context of operator learning. However, its application becomes challenging for complex or high-dimensional processes, as achieving accuracy requires higher-order polynomials, which can increase computational demand and/or the risk of overfitting. Furthermore, PCE requires specialized treatments to manage random variables that are not independent, and these treatments may be problem-dependent or may fail with increasing complexity. Here, in this work, we adopt the same formalism as the spectral expansion used in PCE; however, we replace the classical polynomial basis functions with neural network (NN) basis functions to leverage their expressivity. To achieve this, we propose an algorithm that identifies NN-parameterized basis functions in a purely data-driven manner, without any prior assumptions about the joint distribution of the random variables involved, whether independent or dependent, or about their marginal distributions. The proposed algorithm identifies each NN-parameterized basis function sequentially, ensuring they are orthogonal with respect to the data distribution. The basis functions are constructed directly on the joint stochastic variables without requiring a tensor product structure or assuming independence of the random variables. This approach may offer greater flexibility for complex stochastic models, while simplifying implementation compared to the tensor product structures typically used in PCE to handle random vectors. This is particularly advantageous given the current state of open-source packages, where building and training neural networks can be done with just a few lines of code and extensive community support. We demonstrate the effectiveness of the proposed scheme through several numerical examples of varying complexity and provide comparisons with classical PCE.

Polynomial chaos expansion↗
DOE OSTI2025

Optimal Polynomial Smoothers and One‐Sided V‐Cycles for Poisson Problems

The solution to the Poisson equation arising from the spectral element discretization of the incompressible Navier‐Stokes equations needs robust preconditioning strategies. One such strategy is multigrid. To realize the potential of multigrid methods, effective smoothing strategies are needed. Chebyshev polynomial smoothers, in conjunction with pointwise Jacobi or additive Schwarz methods (ASMs), prove to be an effective smoother. Other polynomial smoothers, however, may provide superior convergence to the multigrid preconditioner. The authors compare the standard Chebyshev polynomial smoothers to both the novel fourth‐kind Chebyshev polynomial smoothers proposed by Lottes as well as smoothers based on the polynomial of best uniform approximation to as proposed by Kraus, Vassilevski, and Zikatanov. At the cost of symmetry, further improvements may be made. For example, a order polynomial smoother on both sides of the V‐cycle may be substituted with an order polynomial smoother on one side at no additional cost. The choice of omitting the postsmoother in favor of higher‐order polynomial presmoothing is advantageous in cases where the multigrid approximation property constant is large. The authors consider a 2D model problem based on finite differences to motivate the choice of polynomial smoother, order, and whether to apply postsmoothing for the target application of high‐order ‐geometric multigrid methods for GPU architectures. Results from both domains demonstrate the substantial improvement of these approaches over the standard Chebyshev polynomial smoother with a symmetric V‐cycle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Pulsed contact glow discharge electrification for online material quantification

Real-time elemental analysis of liquids in extreme environments remains a significant challenge for sensor development, particularly in applications such as advanced nuclear reactors. Plasma-based optical emission spectroscopy offers a promising solution. In this study, we investigate pulsed contact glow discharge electrolysis (CGDE) as an alternative plasma generation method for emission-based material quantification. Using aqueous NaCl solutions, we examine the role of discharge polarity on plasma behavior, emission characteristics, and spectral output. Our results demonstrate that cathodic discharges generate more energetic plasmas, enabling the detection of ionic species such as cerium, while anodic discharges favor atomic emission. These findings reveal a potential dual-mode sensing capability controlled by circuit configuration rather than hardware changes. Through synchronized measurements of emission, current, bubble, and spectral dynamics, we present a model in which electron emission and ohmic heating drive plasma formation. This model enhances our understanding of pulsed-CGDE and supports its potential as a versatile, adaptable platform for future deployment in high-temperature liquids such as molten salts.

98 - NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL↗
DOE OSTI2022

A warm super-Neptune around the G-dwarf star TOI-1710 revealed with TESS, SOPHIE, and HARPS-N

We report the discovery and characterization of the transiting extrasolar planet TOI-1710 b. It was first identified as a promising candidate by the Transiting Exoplanet Survey Satellite. Its planetary nature was then established with SOPHIE and HARPS-N spectroscopic observations via the radial-velocity method. The stellar parameters for the host star are derived from the spectra and a joint Markov chain Monte-Carlo adjustment of the spectral energy distribution and evolutionary tracks of TOI-1710. A joint MCMC analysis of the TESS light curve and the radial-velocity evolution allows us to determine the planetary system properties. From our analysis, TOI 1710 b is found to be a massive warm super-Neptune (M p = 28.3 ± 4.7 M ⊕ and R p = 5.34 ± 0.11 R ⊕ ) orbiting a G5V dwarf star (T eff = 5665 ± 55 K) on a nearly circular 24.3-day orbit (e = 0.16 ± 0.08). The orbital period of this planet is close to the estimated rotation period of its host star P rot = 22.5 ± 2.0 days and it has a low Keplerian semi-amplitude K = 6.4 ± 1.0 m s -1 ; we thus performed additional analyses to show the robustness of the retrieved planetary parameters. With a low bulk density of 1.03 ± 0.23 g cm -3 and orbiting a bright host star (J = 8.3, V = 9.6), TOI-1710 b is one of the best targets in this mass-radius range (near the Neptunian desert) for atmospheric characterization via transmission spectroscopy, a key measurement in constraining planet formation and evolutionary models of sub-Jovian planets.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A numerical study of particle jetting in a dense particle bed driven by an air-blast

In this work, the particle jetting behavior in a blast-driven dense particle bed is studied at early times. Four-way coupled Euler–Lagrange simulations are performed using a high-order discontinuous Galerkin spectral element solver coupled with a high-order Lagrangian particle solver, wherein the inter-particle collisions are resolved using a discrete element method collision model. Following the experiments of Rodriguez et al. [“Formation of particle jetting in a cylindrical shock tube,” Shock Waves 23(6), 619–634 (2013)] and the simulations of Osnes et al. [“Numerical simulation of particle jet formation induced by shock wave acceleration in a Hele-Shaw cell,” Shock Waves 28(3), 451–461 (2018)], the simulations are performed in a quasi-two-dimensional cylindrical geometry (Hele-Shaw cell). Parametric studies are carried out to assess the impact of the coefficient of restitution and the strength of the incident shock on the particle jetting behavior. The deposition of vorticity through a multiphase (gas–particle) analog of Richtmyer–Meshkov instability is observed to play a crucial role in channeling the particles into well-defined jets at the outer edge of the particle bed. This is confirmed by the presence of vortex pairs around the outer jets. Furthermore, the effect of the relaxation of the relative velocity between the two phases on the vorticity generation is explored by analyzing the correlation between the radial velocity of particles and the radial velocity of the gas at the particle location.

42 ENGINEERING↗
DOE OSTI2023

Direct reconstruction of isolated XUV or soft x-ray attosecond pulses from high-harmonic generation streaking spectra

Characterization of an isolated attosecond pulse (IAP) in the extreme ultraviolet (XUV) or soft x-ray (SXR) region is essential for its applications. Here we propose to retrieve an IAP in the time domain directly through the modulation of high-harmonic generation (HHG) spectra in the presence of a time-delayed intense few-cycle infrared or mid-infrared laser. The retrieval algorithm is derived based on the strong-field approximation and an extended quantitative rescattering model. We show that both isolated XUV pulses with a narrow spectral bandwidth and isolated SXR pulses with a broad bandwidth can be well characterized through the HHG streaking spectra. Such an all-optical method for characterizing the IAP differs from the commonly used approach based on the streaked photoelectron spectra that would require electron spectrometers. We check the robustness of the retrieval method by changing the dressing laser or by adjusting the steps of time delay. We also show that the XUV pulse can be accurately retrieved by treating the HHG streaking spectra calculated from solving the time-dependent Schrödinger equation for single atoms as the ‘experimental’ data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Constraints on redshifts of blazars from extragalactic background light attenuation using Fermi - LAT data

ABSTRACT The extragalactic high-energy γ-ray sky is dominated by blazars, which are active galactic nuclei with their jets pointing towards us. Distance measurements are of fundamental importance yet for some of these sources are challenging because any spectral signature from the host galaxy may be outshone by the non-thermal emission from the jet. In this paper, we present a method to constrain redshifts for these sources that relies only on data from the Large Area Telescope on board the Fermi Gamma-ray Space Telescope. This method takes advantage of the signatures that the pair-production interaction between photons with energies larger than approximately 10 GeV and the extragalactic background light leaves on γ-ray spectra. We find upper limits for the distances of 303 γ-ray blazars, classified as 157 BL Lacertae objects, 145 of uncertain class, and 1 flat-spectrum-radio quasar, whose redshifts are otherwise unknown. These derivations can be useful for planning observations with imaging atmospheric Cherenkov telescopes and also for testing theories of supermassive black hole evolution. Our results are applied to estimate the detectability of these blazars with the future Cherenkov Telescope Array, finding that at least 21 of them could be studied in a reasonable exposure of 20 h.

Astronomy & Astrophysics↗