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 847 records · Page 47

Paradigm for global gyrokinetic turbulence

Using a novel wave-number-advection algorithm, we describe a method to add global profile curvature terms to the local gyrokinetic equations. This approach enables a high-precision global simulation capability without sacrificing the efficiency and spectral accuracy of local simulations. Preliminary numerical results suggest that there are two distinct types of global correction to the local result: (1) a trivial shift of the local result in radius, (2) a non-trivial decorrelation that is roughly equivalent to the well-known E×B shear effect. This interpretation provides a new, intuitive framework for understanding global gyrokinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

thornado+FLASH-X: A Hybrid Discontinuous Galerkin–Implicit-explicit and Finite-volume Framework for Neutrino-radiation Hydrodynamics in Core-collapse Supernovae

We present neutrino-transport algorithms implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado) and their coupling to self-gravitating hydrodynamics within the adaptive mesh refinement–based multiphysics simulation framework FLASH-X. thornado, developed primarily for simulations of core-collapse supernovae (CCSNe), employs a spectral, six-species two-moment formulation with algebraic closure and special-relativistic observer corrections accurate to $\mathcal{O}(v/c)$, and uses discontinuous Galerkin (DG) methods for phase-space discretization combined with implicit-explicit time stepping. A key development is a nonlinear neutrino–matter coupling algorithm based on nested fixed-point iteration with Anderson acceleration, enabling fully implicit treatment of collisional processes, including energy-coupling interactions such as neutrino–electron scattering and pair production. Coupling to finite-volume (FV) hydrodynamics is achieved through a hybrid DG-FV representation of the fluid variables and operator-split evolution within FLASH-X. The implementation is verified using basic transport tests with idealized opacities and relaxation and deleptonization problems with tabulated microphysics. Spherically symmetric CCSN simulations demonstrate accuracy and robustness of the coupled scheme, including close agreement with the CCSN simulation code Chimera. An axisymmetric CCSN simulation further demonstrates the viability of DG-based neutrino transport for multidimensional supernova modeling within FLASH-X. thornado’s neutrino-transport solver is GPU-enabled using OpenMP offloading or OpenACC, and all CCSN applications included in this work use the GPU implementation. Together, these results establish a foundation for future enhancements in physics fidelity, numerical algorithms, and computational performance, for increasingly realistic large-scale CCSN simulations.

Endeve, Eirik [Oak Ridge National Laboratory (ORNL↗

Repository scale classification and decomposition of tandem mass spectral data

Various studies have shown associations between molecular features and phenotypes of biological samples. These studies, however, focus on a single phenotype per study and are not applicable to repository scale metabolomics data. Here we report MetSummarizer, a method for predicting (i) the biological phenotypes of environmental and host-oriented samples, and (ii) the raw ingredient composition of complex mixtures. We show that the aggregation of various metabolomic datasets can improve the accuracy of predictions. Since these datasets have been collected using different standards at various laboratories, in order to get unbiased results it is crucial to detect and discard standard-specific features during the classification step. We further report high accuracy in prediction of the raw ingredient composition of complex foods from the Global Foodomics Project.

59 BASIC BIOLOGICAL SCIENCES↗

Hyperspectral radiance mapping and chromatic correction for temperature measurement in laser-heated diamond anvil cells

In laser-heated diamond anvil cell (DAC) experiments, the effective heated region typically decreases in size with increasing pressure, leading to steeper thermal gradients. Under these conditions, chromatic aberration in the optical path from sample to detector can significantly create bias in spectro-radiometric temperature measurement. We present a radiance-mapping approach using a hyperspectral camera that records 25 spectral channels spanning 605–875 nm at each pixel in a single exposure, providing spatially and spectrally resolved radiance in each frame. This enables chromatic effects to be recorded and corrected in data processing. We developed a procedure for hyperspectral mapping, involving per-camera calibration, crosstalk removal, measured spectral throughput functions, and optional sub-pixel co-registration to minimize chromatic distortion. The calibrated radiance maps are then used to derive temperature maps of the laser-heated hotspots. For smaller heating spots, the radiance mapping approach reveals chromatic shifts that conventional spectro-radiometric methods cannot quantify. Ambient-pressure heating experiments confirm accurate temperature retrieval. At high pressure, application of the hyperspectral system to a platinum-heating experiment at 12 GPa demonstrates stable temperature reconstruction under steep thermal gradients. Beyond mitigating chromatic aberrations, the ability to diagnose optical artifacts separately from emissivity variations during controlled test experiments or in situ suggests a path toward more rigorous spectral emissivity analysis and improved modeling of thermal transport in laser-heated DAC experiments.

47 OTHER INSTRUMENTATION↗

Multiclustering needlet ILC for CMB B -mode component separation

ABSTRACT The Cosmic Microwave Background (CMB) primordial B-mode signal is predicted to be much lower than the polarized Galactic emission (foregrounds) in any region of the sky pointing to the need for sophisticated component separation methods. Among them, the blind Needlet Internal Linear Combination (NILC) has great relevance given our current poor knowledge of the B-mode foregrounds. However, the expected level of spatial variability of the foreground spectral properties complicates the NILC subtraction of the Galactic contamination. We therefore propose a novel extension of the NILC approach, the Multiclustering NILC (MC-NILC), which performs NILC variance minimization on separate regions of the sky (clusters) properly chosen to have similar spectral properties of the B-mode Galactic emission within them. Clusters are identified thresholding either the ratio of simulated foregrounds-only B modes (ideal case) or the one of cleaned templates of Galactic emission obtained from realistic simulations. In this work we present an application of MC-NILC to the future LiteBIRD satellite, which targets the observation of both reionization and recombination peaks of the primordial B-mode angular power spectrum with a total error on the tensor-to-scalar ratio δr < 0.001. We show that MC-NILC provides a CMB solution with residual foreground and noise contamination that is significantly lower than the NILC one and the primordial signal targeted by LiteBIRD at all angular scales for the ideal case and at the reionization peak for a realistic ratio. Thus, MC-NILC will represent a powerful method to mitigate B-mode foregrounds for future CMB polarization experiments.

Astronomy & Astrophysics↗

Coupled Cluster Green's function formulations based on the effective Hamiltonians

In this work, we demonstrate that the effective Hamiltonians obtained with the downfolding procedure based on double unitary coupled cluster (DUCC) ansatz can be used in the context of Green’s function coupled cluster (GFCC) formalism to calculate spectral functions of molecular systems. This combined approach (DUCC-GFCC) provides a significant reduction of numerical effort and good agreement with the corresponding all-orbital GFCC methods in energy windows that are consistent with the choice of active space. These features are demonstrated on the example of two benchmark systems: H 2 O and N 2 , where DUCC-GFCC calculations were performed for active spaces of various sizes.

74 ATOMIC AND MOLECULAR PHYSICS↗

Quantitative analysis method for multi-color microscopy

We developed an analysis framework for quantitatively interpreting results from multi-color fluorescence microscopy data. When the emission signal from fluorescent probes (i.e. from biological microscopy or in single-particle microscopy) is spectrally separated before detection, the resulting signals contain information about the strengths of each probe. We built a model to deconvolve the raw data into quantitative descriptions of either of the wavelength of a single probe or the relative strengths among multiple probes. This approach requires only spectrally separating the image data into two or more channels, as opposed to fully resolving the spectral profile, thus allowing other microscopy techniques to be used simultaneously. Our analysis framework incorporates the response features of the optical and imaging elements used to produce the spectral separation, yielding highly quantitative interpretations of the underlying behaviors that give rise to the data.

Ryan, Duncan↗

Ultrabroadband flying-focus using an axiparabola-echelon pair

Flying-focus pulses promise to revolutionize laser-driven secondary sources by decoupling the trajectory of the peak intensity from the native group velocity of the medium over distances much longer than a Rayleigh range. Previous demonstrations of the flying focus have either produced an uncontrolled trajectory or a trajectory that is engineered using chromatic methods that limit the duration of the peak intensity to picosecond scales. Here we demonstrate a controllable ultrabroadband flying focus using a nearly achromatic axiparabola-echelon pair. Spectral interferometry using an ultrabroadband superluminescent diode was used to measure designed super- and subluminal flying-focus trajectories and the effective temporal pulse duration as inferred from the measured spectral phase. The measurements demonstrate that a nearly transform- and diffraction-limited moving focus can be created over a centimeter-scale—an extended focal region more than 50 Rayleigh ranges in length. This ultrabroadband flying-focus and the novel axiparabola-echelon configuration used to produce it are ideally suited for applications and scalable to >100 TW peak powers.

47 OTHER INSTRUMENTATION↗

XANES spectroscopy of sulfides stable under reducing conditions

X-ray absorption near-edge structure (XANES) spectroscopy is a powerful technique to quantitatively investigate sulfur speciation in geologically complex materials such as minerals, glasses, soils, organic compounds, industrial slags, and extraterrestrial materials. This technique allows non-destructive investigation of the coordination chemistry and oxidation state of sulfur species ranging from sulfide (2–oxidation state) to sulfate (6+ oxidation state). Each sulfur species has a unique spectral shape with a characteristic K-edge representing the s → p and d hybridization photoelectron transitions. As such, sulfur speciation is used to measure the oxidation state of samples by comparing the overall XANES spectra to that of reference compounds. Although many S XANES spectral standards exist for terrestrial applications under oxidized conditions, new sulfide standards are needed to investigate reduced (oxygen fugacity, f O 2 , below IW) silicate systems relevant for studies of extraterrestrial materials and systems. Sulfides found in certain meteorites (e.g., enstatite chondrites and aubrites) and predicted to exist on Mercury, such as CaS (oldhamite), MgS (niningerite), and FeCr 2 S 4 (daubréelite), are stable at f O 2 below IW-3 but rapidly oxidize to sulfate and/or produce sulfurous gases under terrestrial surface conditions. XANES spectra of these compounds collected to date have been of variable quality, possibly due to the unstable nature of certain sulfides under typical (e.g., oxidizing) laboratory conditions. A new set of compounds was prepared for this study and their XANES spectra are analyzed for comparison with potential extraterrestrial analogs. S K -edge XANES spectra were collected at Argonne National Lab for FeS (troilite), MnS (alabandite), CaS (oldhamite), MgS (niningerite), Ni 1–x S, NiS 2 , CaSO 4 (anhydrite), MgSO 4 , FeSO 4 , Fe 2 (SO 4 ) 3 , FeCr 2 S 4 (daubréelite), Na 2 S, Al 2 S 3 , Ni 7 S 6 , and Ni 3 S 2 ; the latter five were analyzed for the first time using XANES. These standards expand upon the existing S XANES end-member libraries at a higher spectral resolution (0.25 eV steps) near the S K -edge. Processed spectra, those that have been normalized and “flattened,” are compared to quantify uncertainties due to data processing methods. Future investigations that require well-characterized sulfide standards, such as the ones presented here, may have important implications for understanding sulfur speciation in reduced silicate glasses and minerals with applications for the early Earth, Moon, Mercury, and enstatite chondrites.

Anzures, Brendan A.↗

Metrics for Intercomparison of Remapping Algorithms (MIRA) protocol applied to Earth system models

Abstract. Strongly coupled nonlinear phenomena such as those described by Earth system models (ESMs) are composed of multiple component models with independent mesh topologies and scalable numerical solvers. A common operation in ESMs is to remap or interpolate component solution fields defined on their computational mesh to another mesh with a different combinatorial structure and decomposition, e.g., from the atmosphere to the ocean, during the temporal integration of the coupled system. Several remapping schemes are currently in use or available for ESMs. However, a unified approach to compare the properties of these different schemes has not been attempted previously. We present a rigorous methodology for the evaluation and intercomparison of remapping methods through an independently implemented suite of metrics that measure the ability of a method to adhere to constraints such as grid independence, monotonicity, global conservation, and local extrema or feature preservation. A comprehensive set of numerical evaluations is conducted based on a progression of scalar fields from idealized and smooth to more general climate data with strong discontinuities and strict bounds. We examine four remapping algorithms with distinct design approaches, namely ESMF Regrid (Hill et al., 2004), TempestRemap (Ullrich and Taylor, 2015), generalized moving least squares (GMLS) (Trask and Kuberry, 2020) with post-processing filters, and WLS-ENOR (Li et al., 2020). By repeated iterative application of the high-order remapping methods to the test fields, we verify the accuracy of each scheme in terms of their observed convergence order for smooth data and determine the bounded error propagation using challenging, realistic field data on both uniform and regionally refined mesh cases. In addition to retaining high-order accuracy under idealized conditions, the methods also demonstrate robust remapping performance when dealing with non-smooth data. There is a failure to maintain monotonicity in the traditional L2-minimization approaches used in ESMF and TempestRemap, in contrast to stable recovery through nonlinear filters used in both meshless GMLS and hybrid mesh-based WLS-ENOR schemes. Local feature preservation analysis indicates that high-order methods perform better than low-order dissipative schemes for all test cases. The behavior of these remappers remains consistent when applied on regionally refined meshes, indicating mesh-invariant implementations. The MIRA intercomparison protocol proposed in this paper and the detailed comparison of the four algorithms demonstrate that the new schemes, namely GMLS and WLS-ENOR, are competitive compared to standard conservative minimization methods requiring computation of mesh intersections. The work presented in this paper provides a foundation that can be extended to include complex field definitions, realistic mesh topologies, and spectral element discretizations, thereby allowing for a more complete analysis of production-ready remapping packages.

58 GEOSCIENCES↗

Modeling thin layers of analytes on substrates for spectral analysis: use of solid/liquid n and k values to model reflectance spectra

Since solids are only sometimes seen en masse in a pure bulk form, and for liquids other than water almost never, a capability to model reflectance spectra from analytes deposited on various substrates would be highly advantageous. If available, the real, n ( ν ) , and imaginary, k ( ν ) , components of the complex refractive index, n n + i k , can be used to simulate infrared spectra, accounting for reflection, refraction, and absorption phenomena as a function of wavelength. We focus on using the Pacific Northwest National Laboratory (PNNL) derived n / k vectors for solid and liquid analytes deposited as thin layers on different types of substrates including conductors, such as aluminum, and inorganic dielectrics, such as glass. The model is an adaptation of the Monte Carlo ray trace modeling program, TracePro, extended through the use of its macrolanguage. The model is tested using thin films of organic liquids including silicone oil and no. 2 diesel fuel, as well as organic solids such as caffeine and acetaminophen on aluminum and glass. The predicted spectra for the solid films were compared to experimental hemispherical reflectance data measured using a Fourier transform spectrometer with an integrating sphere. The thickness of the calculated layer is a parameter for predicting the (transflectance) spectra and is obtained using the areal density measured from gravimetric methods to generate the thin-layer samples. Comparison of the calculated spectra with experimental hemispherical reflectance data shows excellent agreement, indicating promise for the use of measured n / k data to synthesize reference spectral data. In particular, accounting for the inhomogeneity of the deposits greatly improved the match with experimental data. Finally, the theoretical modeling shows that for thicker layers (ca. 20 to 100 μm) of typical organics possessing moderately strong k values, the longwave infrared features are often saturated and better spectral contrast is obtained from the overtone/combination bands in the shortwave infrared.

42 ENGINEERING↗

Predicting Switchgrass Biomass Yields Using a Spectral Vegetation Index Derived from Multispectral Satellite Imagery

Successful scaling of perennial bioenergy crop production requires a landscape design that optimizes the benefits of finite lands for people, communities, and environments. Utilizing marginal areas is the key to sustainable bioenergy crop production (Ssegane et al., 2015, 2016). Marginal areas are often small-sized lands and unevenly distributed across the agricultural landscape (Ssegane et al., 2016); thus, a systematic, semi-automated remote sensing method is needed as an effective means of estimating bioenergy crop yields across landscapes. Argonne National Laboratory (Argonne) is currently developing a tool, Scaling Up Perennial Bioenergy Economics and Ecosystem Services Tool (SUPERBEEST), to identify marginal agricultural lands and quantify environmental and economic effects of perennial bioenergy crop production systems. The tool aims to provide users a path to foster the sustainable and productive integration of bioenergy crops in the Midwestern agricultural landscape. Reliable, cost-effective, and timely estimation of bioenergy crop yields using remote sensing would help calculate and track the success of integrated bioenergy crops in the landscape for those communities. Argonne previously conducted feasibility studies for estimating biomass yields for bioenergy feedstock, corn and perennial grass using spectral vegetation indices (SVIs)1 derived from optical imagery (Hamada et al., 2015, 2021). In both studies, SVIs, more specifically those sensitive to plant chlorophyll or nitrogen contents, showed potential for estimating or predicting biomass yields with a correlation of determination (R2) ranging from 0.54 to 0.96, indicating a value for further investigation as a viable means of quantifying bioenergy feedstock production across large landscapes. Thus, the goal of this study is to evaluate the feasibility of use of SVIs as a means of estimating or predicting switchgrass biomass yields at harvest using publicly available multispectral satellite imagery. The feasibility analysis was performed using four study areas of mature switchgrass located in Virginia. Objectives are to (1) examine the SVIs and establish their relationships with switchgrass biomass yields at harvest, (2) develop a parsimonious image processing model for predicting at-harvest yields by applying the relationships with the most promising spectral index and (3) map switchgrass yields predicted by the image processing model across the study sites. The calibration to field data will rely on switchgrass biomass yields determined by the baling method, representing a potential challenge to the analysis but an important practical aspect for future applications. With this research design, the study aimed to gain insights into enabling remote sensing-based estimation of bioenergy crop yields in a reliable, cost- effective, and timely manner across large, heterogeneous landscapes.

09 BIOMASS FUELS↗

Iterative self-organizing SCEne-LEvel sampling (ISOSCELES) for large-scale building extraction

Convolutional neural networks (CNN) provide state-of-the-art performance in many computer vision tasks, including those related to remote-sensing image analysis. Successfully training a CNN to generalize well to unseen data, however, requires training on samples that represent the full distribution of variation of both the target classes and their surrounding contexts. With remote sensing data, acquiring a sufficiently representative training set is a challenge due to both the inherent multi-modal variability of satellite or aerial imagery and the general high cost of labeling data. To address this challenge, we have developed ISOSCELES, an Iterative Self-Organizing SCEne LEvel Sampling method for hierarchical sampling of large image sets. Using affinity propagation, ISOSCELES automates the selection of highly representative training images. Compared to random sampling or using available reference data, the distribution of the training is principally data driven, reducing the chance of oversampling uninformative areas or undersampling informative ones. In comparison to manual sample selection by an analyst, ISOSCELES exploits descriptive features, spectral and/or textural, and eliminates human bias in sample selection. Using a hierarchical sampling approach, ISOSCELES can obtain a training set that reflects both between-scene variability, such as in viewing angle and time of day, and within-scene variability at the level of individual training samples. We verify the method by demonstrating its superiority to stratified random sampling in the challenging task of adapting a pre-trained model to a new image and spatial domain for country-scale building extraction. Using a pair of hand-labeled training sets comprising 1,987 sample image chips, a total of 496,000,000 individually labeled pixels, we show, across three distinct model architectures, an increase in accuracy, as measured by F1-score, of 2.2–4.2%.

42 ENGINEERING↗

On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra

The simplest response theory methods for computing vertical excitation spectra in condensed-phase are configuration interaction with single excitations (CIS) and linear-response time-dependent density functional theory (TDDFT) within the Tamm–Dancoff approximation. In applications to X-ray absorption spectroscopy (XAS), methods like CIS and TDDFT that codify only single excitations into the wave function are prone to catastrophic errors in main-edge and post-edge features whose shapes act as a crucial fingerprint in structural analyses of liquids. Here, we show that these errors manifest primarily due to a lack of orbital relaxation in conventional linear-response theories and that core-ionized (n - 1-electron) references, like those of electron-affinity TDDFT, can eliminate the errors in the spectral profile, even in the highest-energy parts of the post-edge. Crucially, we find that single excitations atop core-ionized references are sufficient to elucidate liquid-phase XAS spectra with semi-quantitative accuracy, opening the door for methods like electron-affinity CIS/TDDFT to be used as efficient alternatives to higher-order wave function approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Intermediate polaronic charge transport in organic crystals from a many-body first-principles approach

Abstract Charge transport in organic molecular crystals (OMCs) is conventionally categorized into two limiting regimes − band transport, characterized by weak electron-phonon (e-ph) interactions, and charge hopping due to localized polarons formed by strong e-ph interactions. However, between these two limiting cases there is a less well understood intermediate regime where polarons are present but transport does not occur via hopping. Here we show a many-body first-principles approach that can accurately predict the carrier mobility in this intermediate regime and shed light on its microscopic origin. Our approach combines a finite-temperature cumulant method to describe strong e-ph interactions with Green-Kubo transport calculations. We apply this parameter-free framework to naphthalene crystal, demonstrating electron mobility predictions within a factor of 1.5−2 of experiment between 100 and 300 K. Our analysis reveals the formation of a broad polaron satellite peak in the electron spectral function and the failure of the Boltzmann equation in the intermediate regime.

36 MATERIALS SCIENCE↗

In Situ ATR-FTIR Study of the Cathode–Electrolyte Interphase: Electrolyte Solution Structure, Transition Metal Redox, and Surface Layer Evolution

We present a study of the lithium nickel manganese cobalt oxide (LiNi 0.6 Mn 0.2 Co 0.2 O 2 , NMC622) cathode-electrolyte interphase (CEI) during galvanostatic charging and discharging using in situ attenuated total reflectance Fourier transform infrared (ATR-FTIR) methods to investigate the voltage dependent electrolyte solution structure changes at the interface, transition metal (TM) redox chemistry, and cathode/electrolyte interfacial layer evolution. Furthermore, our in situ cell design provides both reliable electrochemical device testing and strong FTIR vibrational absorption signals near the cathode surface. Specifically, advanced spectral analysis elucidates changes of near-surface Li+ ion (de)solvation by solvent molecules during galvanostatic cycling. Moreover, cathode metal-oxygen vibrational absorptions, sensitive to TM redox behaviors and subsequent local structural variations, were correlated to cathode de-lithiation (and lithiation) and electrolyte solution structure changes. In addition, we have detected the formation and evolution of a CEI surface layer on the NMC622 cathode that contributes to the cell's capacity fade.

25 ENERGY STORAGE↗

Data-Driven Operator Theoretic Methods for Phase Space Learning and Analysis

This paper uses data-driven operator theoretic approaches to explore the global phase space of a dynamical system. In this work, we defined conditions for discovering new invariant subspaces in the state space of a dynamical system starting from an invariant subspace based on the spectral properties of the Koopman operator. When the system evolution is known locally in several invariant subspaces in the state space of a dynamical system, a phase space stitching result is derived that yields the global Koopman operator. Additionally, in the case of equivariant systems, a phase space stitching result is developed to identify the global Koopman operator using the symmetry properties between the invariant subspaces of the dynamical system and time-series data from any one of the invariant subspaces. Finally, these results are extended to topologically conjugate dynamical systems; in particular, the relation between the Koopman tuple of topologically conjugate systems is established. The proposed results are demonstrated on several second-order nonlinear dynamical systems including a bistable toggle switch. Our method elucidates a strategy for designing discovery experiments: experiment execution can be done in many steps, and models from different invariant subspaces can be combined to approximate the global Koopman operator.

42 ENGINEERING↗

Charge self-consistent density functional theory plus ghost rotationally invariant slave-boson theory for correlated materials

We present a charge self-consistent density functional theory combined with the ghost rotationally invariant slave-boson (DFT+gRISB) formalism for studying correlated materials. Here, this method is applied to SrVO 3 and NiO, representing prototypical correlated metals and charge-transfer insulators. For SrVO 3 , we demonstrate that DFT+gRISB yields an accurate equilibrium volume and effective mass close to experimentally observed values. Regarding NiO, DFT+gRISB enables the simultaneous description of charge-transfer and Mott-Hubbard bands, significantly enhancing the accuracy of the original DFT+RISB approach. Furthermore, the calculated equilibrium volume and spectral function reasonably agree with experimental observations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗