Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spectral modeling”

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

Color–color diagrams as tools for assessment of the variable absorption in high mass X-ray binaries

High mass X-ray binaries hold the promise of allowing us to understand the structure of the winds of their supermassive companion stars by using the emission from the compact object as a backlight to evaluate the variable absorption in the structured stellar wind. The wind along the line of sight can change on timescales as short as minutes and below. However, such short timescales are not available for the direct measurement of absorption through X-ray spectroscopy with the current generation of X-ray telescopes. In this paper, we demonstrate the usability of color–color diagrams for assessing the variable absorption in wind accreting high mass X-ray binary systems. Furthermore, we employ partial covering models to describe the spectral shape of high mass X-ray binaries and assess the implication of different absorbers and their variability on the shape of color–color tracks. We show that in taking into account, the ionization of the absorber, and in particular accounting for the variation of ionization with absorption depth, is crucial to describe the observed behavior well.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A simple generalization of the energy gap law for nonradiative processes

For more than 50 years, an elegant energy gap (EG) law developed by Englman and Jortner [Mol. Phys. 18, 145 (1970)] has served as a key theory to understand and model the nearly exponential dependence of nonradiative transition rates on the difference of energy between the initial and final states. This work revisits the theory, clarifies the key assumptions involved in the rate expression, and provides a generalization for the cases where the effects of temperature dependence and low-frequency modes cannot be ignored. For a specific example where the low-frequency vibrational and/or solvation responses can be modeled as an Ohmic spectral density, a simple generalization of the EG law is provided. Test calculations demonstrate that this generalized EG law brings significant improvement over the original EG law. Both the original and generalized EG laws are also compared with the stationary phase approximations developed for electron transfer theory, which suggests the possibility of a simple interpolation formula valid for any value of EG.

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

Constraining chaos: Enforcing dynamical invariants in the training of reservoir computers

Drawing on ergodic theory, we introduce a novel training method for machine learning based forecasting methods for chaotic dynamical systems. The training enforces dynamical invariants—such as the Lyapunov exponent spectrum and the fractal dimension—in the systems of interest, enabling longer and more stable forecasts when operating with limited data. The technique is demonstrated in detail using reservoir computing, a specific kind of recurrent neural network. Finally, results are given for the Lorenz 1996 chaotic dynamical system and a spectral quasi-geostrophic model of the atmosphere, both typical test cases for numerical weather prediction.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Diagnosing up-scattered deuterium–tritium fusion neutrons produced in burning plasmas at the National Ignition Facility (invited)

In the push to higher performance fusion plasmas, two critical quantities to diagnose are α-heat deposition that can improve and impurities mixed into the plasma that can limit performance. In high-density, highly collisional inertial confinement fusion burning plasmas, there is a significant probability that deuterium–tritium (DT) fusion products, 14.1 MeV neutrons and 3.5 MeV α-particles, will collide with and deposit energy onto (“up-scatter”) surrounding deuterium and tritium fuel ions. These up-scattered D and T ions can then undergo fusion while in-flight and produce an up-scattered neutron (15–30 MeV). These reaction-in-flight (RIF) neutrons can then be uniquely identified in the measured neutron energy spectrum. Further, the magnitude, shape, and relative size of this spectral feature can inform models of stopping-power in the DT plasma and hence is directly proportional to α-heat deposition. In addition, the RIF spectrum can be related to mix into the burning fuel, particularly relevant for high-Z shell and other emerging National Ignition Facility platforms. The neutron time-of-flight diagnostic upgrades needed to obtain this small signal, ~10 –5 times the primary DT neutron peak, will be discussed. Results from several gain > 1 implosions will be shown and compared to previous RIF spectra. Finally, comparisons of experimental data to a simplified computational model will be made.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Introduction to spectral line shape theory

Spectral line-shape models are an important part of understanding high-energy-density (HED) plasmas. Models are needed for calculating opacity of materials and can serve as diagnostics for astrophysical and laboratory plasmas. However, much of the literature on line shapes is directed toward specialists. This perspective makes it difficult for non-specialists to enter the field. We have two broad goals with this topical review. First, we aim to give information so that others in HED physics may better understand the current field. Additionally, this first goal may help guide future experiments to test different aspects of the theory. Second, we provide an introduction for those who might be interested in line-shape theory, and enough materials to be able to navigate the field and the literature. We give a high-level overview of line broadening process, as well as dive into the formalism, available methods, and approximations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Polariton creation in coupled cavity arrays with spectrally disordered emitters

Abstract Integrated photonics has been a promising platform for analog quantum simulation of condensed matter phenomena in strongly correlated systems. To that end, we explore the implementation of all-photonic quantum simulators in coupled cavity arrays with integrated ensembles of spectrally disordered emitters. Our model is reflective of color center ensembles integrated into photonic crystal cavity arrays. Using the Quantum Master equation and the Effective Hamiltonian approaches, we study energy band formation and wavefunction properties in the open quantum Tavis–Cummings–Hubbard framework. We find conditions for polariton creation and (de)localization under experimentally relevant values of disorder in emitter frequencies, cavity resonance frequencies, and emitter-cavity coupling rates. To quantify these properties, we introduce two metrics, the polaritonic and nodal participation ratios, that characterize the light-matter hybridization and the node delocalization of the wavefunction, respectively. These new metrics combined with the Effective Hamiltonian approach prove to be a powerful toolbox for cavity quantum electrodynamical engineering of solid-state systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Dust from evolved stars: a pilot analysis of the AGB to PN transition

ABSTRACT We present a novel approach to address dust production by low- and intermediate-mass stars. We study the asymptotic giant branch (AGB) phase, during which the formation of dust takes place, from the perspective of post-AGB and planetary nebula (PN) evolutionary stage. Using results from stellar evolution and dust formation modelling, we interpret the spectral energy distribution of carbon-dust-rich sources currently evolving through different evolutionary phases, believed to descend from progenitors of similar mass and chemical composition. Comparing the results of different stages along the AGB to PNe transition, we can provide distinct insights on the amount of dust and gas released during the very late AGB phases. While the post-AGB traces the history of dust production back to the tip of the AGB phase, investigating the PNe is important to reconstruct the mass-loss process experienced after the last thermal pulse. The dust surrounding the post-AGB was formed soon after the tip of the AGB. The PNe dust-to-gas ratio is ∼10−3, 2.5 times smaller than what expected for the same initial mass star during the last AGB interpulse, possibly suggesting that dust might be destroyed during the PN phase. Measuring the amount of dust present in the nebula can constrain the capacity of the dust to survive the central star heating.

Dell’Agli, F. (ORCID:0000000324426981)↗
DOE OSTI2025

An obscured quasar census with the 4MOST IR AGN survey: design, predicted properties, and scientific goals

ABSTRACT We present the 4MOST (4-metre Multi-Object Spectroscopic Telescope) infrared (IR) AGN survey, the first large-scale optical spectroscopic survey characterizing mid-infrared (MIR) selected obscured active galactic nuclei (AGNs). The survey targets $\approx 212\,000$ obscured IR AGN candidates over $\approx 10\,000 \rm \: deg^2$ down to a magnitude limit of $r_{\rm AB}=22.8 \, \rm mag$ and will be $\approx 100 \times$ larger than any existing obscured IR AGN spectroscopic sample. We select the targets using an MIR colour criterion applied to the unWISE catalogue from the WISE (Wide-field Infrared Survey Explorer) all-sky survey, and then apply a $r-W2\ge 5.9 \rm \: mag$ cut; we demonstrate that this selection will mostly identify sources obscured by $N_{\rm H}>10^{22} \rm \: cm^{-2}$. The survey complements the 4MOST X-ray survey, which will follow up $\sim 1\,\rm M$ eROSITA (extended ROentgen Survey with an Imaging Telescope Array)-selected (typically unobscured) AGN. We perform simulations to predict the quality of the spectra that we will obtain and validate our MIR–optical colour-selection method using X-ray spectral constraints and UV-to-far-IR spectral energy distribution (SED) modelling in four well-observed deep-sky fields. We find that: (1) $\approx 80-87{{\ \rm per\ cent}}$ of the WISE-selected targets are AGN down to $r_{\rm AB}=22.1-22.8 \: \rm mag$ of which $\approx 70{{\ \rm per\ cent}}$ are obscured by $N_{\rm H}>10^{22} \: \rm cm^{-2}$, and (2) $\approx 80{{\ \rm per\ cent}}$ of the 4MOST IR AGN sample will remain undetected by the deepest eROSITA observations due to extreme absorption. Our SED-fitting results show that the 4MOST IR AGN survey will primarily identify obscured AGN and quasars ($\approx 55{{\ \rm per\ cent}}$ of the sample is expected to have $L_{\rm AGN,IR}>10^{45} \rm \: erg \: s^{-1}$) residing in massive galaxies ($M_{\star }\approx 10^{10}-10^{12} \rm \: M_{\odot }$) at $z\approx 0.5-3.5$ with $\approx 33{{\ \rm per\ cent}}$ expected to be hosted by starburst galaxies.

Andonie, Carolina (ORCID:0000000255804298)↗
DOE OSTI2025

Realizability-preserving discontinuous Galerkin method for spectral two-moment radiation transport in special relativity

Here we present a realizability-preserving numerical method for solving a spectral two-moment model to simulate the transport of massless, neutral particles interacting with a steady background material moving with relativistic velocities. The model is obtained as the special relativistic limit of a four-momentum-conservative general relativistic two-moment model. Using a maximum-entropy closure, we solve for the Eulerian-frame energy and momentum. The proposed numerical method is designed to preserve moment realizability, which corresponds to moments defined by a nonnegative phase-space density. The realizability-preserving method is achieved with the following key components: (i) a discontinuous Galerkin phase-space discretization with specially constructed numerical fluxes in the spatial and energy dimensions; (ii) a strong stability-preserving implicit-explicit time-integration method; (iii) a realizability-preserving conserved to primitive moment solver; (iv) a realizability-preserving implicit collision solver; and (v) a realizability-enforcing limiter. Component (iii) is necessitated by the closure procedure, which closes higher order moments nonlinearly in terms of primitive moments. The nonlinear conserved to primitive and the implicit collision solves are formulated as fixed-point problems, which are solved with custom iterative solvers designed to preserve the realizability of each iterate. With a series of numerical tests, we demonstrate the accuracy and robustness of this discontinuous-Galerkin-implicit-explicit method.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Equating Reference Cell and Broadband Pyranometer Data

An adjustment algorithm has been developed for broadband Class-A pyranometers, to modify the pyranometer data such that it more closely matches that of silicon-based reference cells. The adjustment algorithm uses modeled clear-sky spectral data as one of its inputs. The adjustment also uses a reference cell angle of incidence modifier. A comparison between the pyranometer and the reference cells is performed with and without the adjustment applied. The adjustment aligns the pyranometer and reference cell data to within +-2 % under clear skies when the zenith angle is less than 70 degrees .

incident angle modifier↗
DOE OSTI2023

Semi-Implicit Computation of Fast Modes in a Scheme Integrating Slow Modes by a Leapfrog Method Based on a Selective Implicit Time Filter

Abstract A scheme for integration of atmospheric equations containing terms with differing time scales is developed. The method employs a filtered leapfrog scheme utilizing a fourth-order implicit time filter with one function evaluation per time step to compute slow-propagating phenomena such as advection and rotation. The terms involving fast-propagating modes are handled implicitly with an unconditionally stable method that permits application of larger time steps and faster computations compared to fully explicit treatment. Implementation using explicit and recurrent formulation is provided. Stability analysis demonstrates that the method is conditionally stable for any combination of frequencies involved in the slow and fast terms as they approach the origin. The implicit filter used in the method damps the computational modes without noticeably sacrificing the accuracy of the physical mode. TheO[(Δt 4 )] accuracy for amplitude errors achieved by the implicitly filtered leapfrog is preserved in applications where terms responsible for fast propagation are integrated with a semi-implicit method. Detailed formulation of the method for soundproof nonhydrostatic anelastic equations is provided. Procedures for implementation in global spectral shallow-water models are also given. Examples comparing numerical and analytical solutions for linear gravity waves demonstrate the accuracy of the scheme. The performance is also shown in more practical nonlinear applications, where numerical solutions accomplished by the method are evaluated against those computed from a scheme where the slow terms are handled by the third-order Runge–Kutta scheme. It demonstrates that the method is able to accurately resolve fine-scale dynamics of Kelvin–Helmholtz shear instabilities, the evolution of density current, and nonlinear drifts of twin tropical cyclones.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Spring 2022 update on the status of the Local Wavenumber Model (LWN) in xRAGE

An updated implementation of the Local Wavenumber Model (LWN) is discussed, primarily differing from recent versions by placing a greater focus on capturing a wide variety of turbulent flows including compressible flows. New models are introduced for spectral backscatter effects, the effect of bulk compression on the spectra, and incorporating the multispecies variables tracked in BHR4. Methods for reducing the compurational expense of tracking spectra for turbulent quantites are also investigated. Like recent versions of BHR, we test LWN in a number of canonical flows using a single set of coefficients, but additional coefficient tuning is likely to be required to improve the agreement with some of these flows.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Bottom-Up Approach to Rational Design of Crystalline Materials: Investigation of Vibronic Coherences Underlying Exciton Dynamics in Semiconductors

In this project we uncovered structure-function relationships of donor-acceptor co-crystals used to develop next-generation optoelectronic devices. Unraveling the photodynamics of molecular crystalline materials poses many challenges for spectroscopy due to broad, overlapping features representing numerous underlying dynamical processes. This leads researchers to make many assumptions about the dynamics of a system in choosing an appropriate kinetic fitting model. Computationally, electronic structure methods are either prohibitively expensive or underdeveloped for computing the excited state structure of molecular materials, especially states that exhibit charge transfer. Researchers must therefore perform calculations of excited electronic states using truncated models of molecular materials. Here we present a joint experimental-theoretical approach to bridging the gap between the photodynamics of a molecular material and its constituent molecules. We focus our efforts on quantifying the timescales and mechanisms of photoexcitation in donor-acceptor co-crystals and donor-acceptor dimers where the lowest-lying excited state is characterized by charge transfer from the donor to the acceptor. We employ ultrafast UV pump, UV-Vis probe transient absorption spectroscopy to unravel the time-resolved spectroscopic signatures of the photodynamics in both the crystalline material and donor-acceptor dimers in solution. We perform electronic structure and excited state dynamics calculations of the dimers to inform kinetic fitting models and assign the spectral features. The photodynamics of the crystal vs. dimer systems have many similarities, enabling unprecedented insights into the formation and evolution of charge transfer excitons in the crystalline systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

New Projections for Dark Matter Searches with Paleo-Detectors

Paleo-detectors are a proposed experimental technique to search for dark matter (DM). In lieu of the conventional approach of operating a tonne-scale real-time detector to search for DM-induced nuclear recoils, paleo-detectors take advantage of small samples of naturally occurring rocks on Earth that have been deep underground (≳5 km), accumulating nuclear damage tracks from recoiling nuclei for O(1)Gyr. Modern microscopy techniques promise the capability to read out nuclear damage tracks with nanometer resolution in macroscopic samples. Thanks to their O(1)Gyr integration times, paleo-detectors could constitute nuclear recoil detectors with keV recoil energy thresholds and 100 kilotonne-yr exposures. This combination would allow paleo-detectors to probe DM-nucleon cross sections orders of magnitude below existing upper limits from conventional direct detection experiments. In this article, we use improved background modeling and a new spectral analysis technique to update the sensitivity forecast for paleo-detectors. We demonstrate the robustness of the sensitivity forecast to the (lack of) ancillary measurements of the age of the samples and the parameters controlling the backgrounds, systematic mismodeling of the spectral shape of the backgrounds, and the radiopurity of the mineral samples. Specifically, we demonstrate that even if the uranium concentration in paleo-detector samples is 10−8 (per weight), many orders of magnitude larger than what we expect in the most radiopure samples obtained from ultra basic rock or marine evaporite deposits, paleo-detectors could still probe DM-nucleon cross sections below current limits. For DM masses ≲ 10 GeV/c2, the sensitivity of paleo-detectors could still reach down all the way to the conventional neutrino floor in a Xe-based direct detection experiment.

Baum, Sebastian (ORCID:0000000167929381)↗
DOE OSTI2020

Effective Opacity of the Intergalactic Medium from Galaxy Spectra Analysis

We measure the effective opacity (τ{sub eff}) of the intergalactic medium from the composite spectra of 281 Lyman-break galaxies in the redshift range 2 ≲ z ≲ 3. Our spectra are taken from the COSMOS Lyα Mapping And Tomographic Observations survey derived from the Low Resolution Imaging Spectrometer on the W.M. Keck I telescope. We generate composite spectra in two redshift intervals and fit them with spectral energy distribution (SED) models composed of simple stellar populations. Extrapolating these SED models into the Lyα forest, we measure the effective Lyα opacity (τ{sub eff}) in the 2.02 ≤ z ≤ 2.44 range. At z = 2.22, we estimate τ{sub eff} =0.159±0.001 from a power-law fit to the data. These measurements are consistent with estimates from quasar analyses at z < 2.5 indicating that the systematic errors associated with normalizing quasar continua are not substantial. We provide a Gaussian processes model of our results and previous τ{sub eff} measurements that describes the steep redshift evolution in τ{sub eff} from z = 1.5–4.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

First Temperature Profile of a Stellar Flare Using Differential Chromatic Refraction

We present the first derivation of a stellar flare temperature profile from single-band photometry. Stellar flare DWF 030225.574−545707.45129 was detected in 2015 by the Dark Energy Camera as part of the Deeper, Wider, Faster program. The brightness (Δm g = −6.12) of this flare, combined with the high air mass (1.45 ≲ X ≲ 1.75) and blue filter (DES g, 398–548 nm) in which it was observed, provided ideal conditions to measure the zenithward apparent motion of the source due to differential chromatic refraction (DCR) and, from that, infer the effective temperature of the event. We model the flare’s spectral energy distribution as a blackbody to produce the constraints on flare temperature and geometric properties derived from single-band photometry. We additionally demonstrate how simplistic assumptions on the flaring spectrum, as well as on the evolution of flare geometry, can result in solutions that overestimate the effective temperature. Exploiting DCR enables studying chromatic phenomena with ground-based astrophysical surveys, and stellar flares on M dwarfs are a particularly enticing target for such studies due to their ubiquity across the sky and the heightened color contrast between their red quiescent photospheres and the blue flare emission. Our novel method will enable similar temperature constraints for a large sample of objects in upcoming photometric surveys like the Vera C. Rubin Legacy Survey of Space and Time.

Clarke, Riley W. [University of Delaware, Newark, ↗
DOE OSTI2022

Islet: interpolation semi-Lagrangian element-based transport

Abstract. Advection of trace species, or tracers, also called tracer transport, in models of the atmosphere and other physical domains is an important and potentially computationally expensive part of a model's dynamical core. Semi-Lagrangian (SL) advection methods are efficient because they permit a time step much larger than the advective stability limit for explicit Eulerian methods without requiring the solution of a globally coupled system of equations as implicit Eulerian methods do. Thus, to reduce the computational expense of tracer transport, dynamical cores often use SL methods to advect tracers. The class of interpolation semi-Lagrangian (ISL) methods contains potentially extremely efficient SL methods. We describe a finite-element ISL transport method that we call the interpolation semi-Lagrangian element-based transport (Islet) method, such as for use with atmosphere models discretized using the spectral element method. The Islet method uses three grids that share an element grid: a dynamics grid supporting, for example, the Gauss–Legendre–Lobatto basis of degree three; a physics parameterizations grid with a configurable number of finite-volume subcells per element; and a tracer grid supporting use of Islet bases with particular basis again configurable. This method provides extremely accurate tracer transport and excellent diagnostic values in a number of verification problems.

58 GEOSCIENCES↗
DOE OSTI2020

thornado-transport: Anderson- and GPU-accelerated nonlinear solvers for neutrino-matter coupling

Algorithms for neutrino-matter coupling in core-collapse supernovae (CCSNe) are investigated in the context of a spectral two-moment model, which is discretized in space with the discontinuous Galerkin method, integrated in time with implicit-explicit (IMEX) methods, and implemented in the toolkit for high-order neutrino-radiation hydrodynamics (thornado). The model considers electron neutrinos and antineutrinos and tabulated opacities from Bruenn (1985), which includes neutrino-electron scattering and pair processes. The nonlinear system arising from implicit time discretization of the equations governing neutrino-matter coupling is iterated to convergence using Anderson-accelerated fixed-point methods, which avoid formation of Jacobians and inversion of dense linear systems. Numerical experiments show that, for a given tolerance, a nested iteration scheme which aims to reduce opacity evaluations can lower the computational cost. Our initial port to GPUs, using both OpenMP and OpenACC, shows an overall speedup of up to ~ 100× when compared to results using a single CPU core. These results indicate that the algorithms implemented in thornado are well-suited to GPU acceleration.

Laiu, Paul↗