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 811 records · Page 45

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↗

Affinity propagation clustering of full-field, high-spatial-dimensional measurements for robust output-only modal identification: A proof-of-concept study

Determination of the model order is a challenging problem in system identification, especially in output-only or operational modal identification where some modes are weakly excited. While existing methods such as the stabilization diagram method (spectral information) are effective, they do not scale to high-dimensional data, which is usually needed for high-fidelity characterization of structural dynamics and has been made available in the emerging full-field measurement techniques using optical methods such as photogrammetry and laser vibrometers. In this proof-of-concept study we present a new non-parametric, data-driven approach for robust output-only identification of high-spatial-dimensional modal parameters of basic structures by efficiently processing and interactively exploiting the full-field measurement (i.e., very dense spatial measurement points). Specifically, we first over-estimate the system model once, producing a pool of candidate modes associated with their modal frequencies and full-field, high-spatial-dimensional mode shapes. This is accomplished by a data-driven method termed affinity propagation clustering (APC), where the active clusters, which are the active modes in our formulations, emerge from the "message-passing" procedure and does not require a pre-determination of the cluster number (mode or model order). Next, rather than using the spectral information to distinguish the physical and spurious modes in the stabilization diagram method, we exploit and visualize the spatial, full-field mode shape associated with each candidate mode to do so. We conduct extensive experiments on basic structural models with comparisons to a few existing methods. The results indicate that the new method is computationally efficient for identifying high-spatial-dimensional modal parameters, and robust to identify weak modes by exploiting the full-field measurement. We also discuss its applicability and limitations for structures with complex geometry (shapes).

42 ENGINEERING↗

Unified Wind-Wave Growth and Spectrum Functions for All Water Depths: Field Observations and Model Results

Abstract Wind-wave development is governed by the fetch- or duration-limited growth principle that is expressed as a pair of similarity functions relating the dimensionless elevation variance (wave energy) and spectral peak frequency to fetch or duration. Combining the pair of similarity functions, the fetch or duration variable can be removed to form a dimensionless function of elevation variance and spectral peak frequency, which is interpreted as the wave energy evolution with wave age. The relationship is initially developed for quasi-neural stability and quasi-steady wind forcing conditions. Further analyses show that the same fetch, duration, and wave-age similarity functions are applicable to unsteady wind forcing conditions, including rapidly accelerating and decelerating mountain gap wind episodes and tropical cyclone (TC) wind fields. Here it is shown that with the dimensionless frequency converted to dimensionless wavenumber using the surface wave dispersion relationship, the same similarity function is applicable in all water depths. Field data collected in shallow to deep waters and mild to TC wind conditions and synthetic data generated by spectrum model computations are assembled to illustrate the applicability. For the simulation work, the finite-depth wind-wave spectrum model and its shoaling function are formulated for variable spectral slopes. Given wind speed, wave age, and water depth, the measured and spectrum-computed significant wave heights and the associated growth parameters are in good agreement in forcing conditions from mild to TC winds and in all depths from deep ocean to shallow lake. Significance Statement This paper presents a growth function and spectrum model to describe wind-wave development in all water depths. Their applicability covers a wide range of wind forcing conditions including steady, accelerating, decelerating, and tropical cyclone events. Support for the unified spectrum model and growth function is presented with field observations and numerical computations.

Hwang, Paul A.↗

On reduced modelling of the modulational dynamics in magnetohydrodynamics

This paper explores structure formation in two-dimensional magnetohydrodynamic (MHD) turbulence as a modulational instability (MI) of turbulent fluctuations. We focus on the early stages of structure formation and consider simple backgrounds that allow for a tractable model of the MI while retaining the full chain of modulational harmonics. This approach allows us to systematically examine the validity of popular closures such as the quasilinear approximation and other low-order truncations. We find that, although such simple closures can provide quantitatively accurate approximations of the MI growth rates in some regimes, they can fail to capture the modulational dynamics in adjacent regimes even qualitatively, falsely predicting MI when the system is actually stable. We find that this discrepancy is due to the excitation of propagating spectral waves (PSWs) which can ballistically transport energy along the modulational spectrum, unimpeded until dissipative scales, thereby breaking the feedback loops that would otherwise sustain MIs. The PSWs can be self-maintained as global modes with real frequencies and drain energy from the primary structure at a constant rate until the primary structure is depleted. To describe these waves within a reduced model, we propose an approximate spectral closure that captures them and MIs on the same footing. We also find that introducing corrections to ideal MHD, conservative or dissipative, can suppress PSWs and reinstate the accuracy of the quasilinear approximation. In this sense, ideal MHD is a ‘singular’ system that is particularly sensitive to the accuracy of the closure within mean-field models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Landsat near-infrared (NIR) band and ELM-FATES sensitivity to forest disturbances and regrowth in the Central Amazon

Forest disturbance and regrowth are key processes in forest dynamics, but detailed information on these processes is difficult to obtain in remote forests such as the Amazon. We used chronosequences of Landsat satellite imagery (Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper Plus) to determine the sensitivity of surface reflectance from all spectral bands to windthrow, clear-cut, and clear-cut and burned (cut + burn) and their successional pathways of forest regrowth in the Central Amazon. We also assessed whether the forest demography model Functionally Assembled Terrestrial Ecosystem Simulator (FATES) implemented in the Energy Exascale Earth System Model (E3SM) Land Model (ELM), ELM-FATES, accurately represents the changes for windthrow and clear-cut. The results show that all spectral bands from the Landsat satellites were sensitive to the disturbances but after 3 to 6 years only the near-infrared (NIR) band had significant changes associated with the successional pathways of forest regrowth for all the disturbances considered. In general, the NIR values decreased immediately after disturbance, increased to maximum values with the establishment of pioneers and early successional tree species, and then decreased slowly and almost linearly to pre-disturbance conditions with the dynamics of forest succession. Statistical methods predict that NIR values will return to pre-disturbance values in about 39, 36, and 56 years for windthrow, clear-cut, and cut + burn disturbances, respectively. The NIR band captured the observed, and different, successional pathways of forest regrowth after windthrow, clear-cut, and cut + burn. Consistent with inferences from the NIR observations, ELM-FATES predicted higher peaks of biomass and stem density after clear-cuts than after windthrows. ELM-FATES also predicted recovery of forest structure and canopy coverage back to pre-disturbance conditions in 38 years after windthrows and 41 years after clear-cut. The similarity of ELM-FATES predictions of regrowth patterns after windthrow and clear-cut to those of the NIR results suggests the NIR band can be used to benchmark forest regrowth in ecosystem models. Our results show the potential of Landsat imagery data for mapping forest regrowth from different types of disturbances, benchmarking, and the improvement of forest regrowth models.

54 ENVIRONMENTAL SCIENCES↗

Wolf-Rayet Galaxies in SDSS-IV MaNGA. II. Metallicity Dependence of the High-mass Slope of the Stellar Initial Mass Function

As hosts of living high-mass stars, Wolf-Rayet (WR) regions or WR galaxies are ideal objects for constraining the high-mass end of the stellar initial mass function (IMF). We construct a large sample of 910 WR galaxies/regions that cover a wide range of stellar metallicity (from Z ~ 0.001 to 0.03) by combining three catalogs of WR galaxies/regions previously selected from the SDSS and SDSS-IV/MaNGA surveys. We measure the equivalent widths of the WR blue bump at ~4650 Å for each spectrum. They are compared with predictions from stellar evolutionary models Starburst99 and BPASS, with different IMF assumptions (high-mass slope α of the IMF ranging from 1.0 to 3.3). Both singular evolution and binary evolution are considered. We also use a Bayesian inference code to perform full spectral fitting to WR spectra with stellar population spectra from BPASS as fitting templates. We then make a model selection among different α assumptions based on Bayesian evidence. These analyses have consistently led to a positive correlation of the IMF high-mass slope α with stellar metallicity Z, i.e., with a steeper IMF (more bottom-heavy) at higher metallicities. Specifically, an IMF with α = 1.00 is preferred at the lowest metallicity (Z ~ 0.001), and an Salpeter or even steeper IMF is preferred at the highest metallicity (Z ~ 0.03). These conclusions hold even when binary population models are adopted.

79 ASTRONOMY AND ASTROPHYSICS↗

Simulating a Mesoscale Convective System Using WRF with a New Spectral Bin Microphysics: Part 1: Hail vs Graupel

A modified Fast Spectral Bin Microphysics scheme (FSBM-2) embedded into the Weather Research and Forecasting Model (WRF) is used to simulate a mesoscale deep convective system observed during the Midlatitude Continental Convective Clouds Experiment (MC3E). FSBM-2 uses modified source codes as compared to the current FSBM (FSBM-1). In contrast to FSBM-1, FSBM-2 can simulate hail of several cm in diameter and includes additional processes such as spontaneous breakup of raindrops and aerosol regeneration by drop evaporation. It is shown that allowing large hail particles of diameters exceeding about 1 cm substantially increases the agreement between the simulated and observed squall-line structures both in the convective and stratiform regions. In contrast, if graupel particles are used to represent high-density hydrometeors in convective areas, the ratio of convective-to-stratiform areas diverges from the ratio seen in observations and maximum radar reflectivities are substantially underestimated. Analysis of snow size distributions in the stratiform area shows an important link between the ice crystals formed by homogeneous freezing in the convective area to ice particle number concentration in the stratiform region. Simulated raindrop size distributions in the stratiform area below the melting level from FSBM-2 show a good agreement with observations. The regeneration of CCN by droplet evaporation and detrainment of these CCN to the stratiform region increases the concentration of CCN there up to 8-9 km altitude; the additional CCN penetrate clouds and produce new droplets, leading to some convection intensification. Key words: Mesoscale convective systems (MCS), squall lines, spectral bin microphysics, large hail, in-cloud nucleation, size distributions of hydrometeors in MCS.

Shpund, Kobby↗

Enhanced quantum efficiency from optical interference in alkali antimonide photocathodes: Modeling and experimental results

We present measurements of enhanced quantum efficiency (QE) in thin film alkali antimonide photocathodes from optical interference in the cathode-substrate multilayer. Modulations in the spectral response are observed over a range of visible wavelengths and are shown to increase the QE by more than a factor of two at specific wavelengths. We present a model describing the QE modulations based on the three step photoemission process incorporating cases of both constant density of states and density functional theory-derived density of states and show that the calculated results are in good agreement with the measurements. Model predictions demonstrate that QE can be enhanced by more than a factor of 5 by optimization of cathode and substrate layer thicknesses. Additionally, these calculations reveal that optical interference can yield higher quantum efficiencies in thin films compared to thick, optically dense films. We model the QE vs excitation wavelength of multiple alkali antimonide compounds at different thicknesses. We then discuss the advantages of this interference effect for electron accelerators.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Gamma-ray Spectrum Explanations (GRSE) v0.1

This code repository provides tools for producing various types of machine learning model explanations trained on gamma-ray spectra. It includes implementations of saliency mapping, Grad-CAM, LIME, and Kernel SHAP, including improvements and clarifications to some of the approaches that are specific to their use with gamma-ray spectral data. These tools are in support of the manuscript M.S. Bandstra et al., "Explaining machine-learning models for gamma-ray detection and identification," under review at PLOS ONE.

Bandstra, Mark↗

Anapole moment of neutrinos and radioactive sources near liquid xenon detectors

We show that placing a radioactive source such as 51 Cr near a liquid xenon detector may allow us to detect the contribution induced by the anapole moment to neutrino-electron scattering in the Standard Model (SM) at the 1 − 2⁢𝜎 level. Although the anapole moment of neutrinos induces a scattering rate with the same spectral shape as the neutral and charged current contributions, exposures of ∼ 60 ton × source run at XENONnT or XLZD may be enough to accumulate sufficient statistics for a detection. We also discuss a simple model where the anapole moment of neutrinos is enhanced or decreased with respect to the SM expectation, further demonstrating how a potential measurement of the anapole moment of neutrinos would allow us to constrain new physics.

dark matter direct detection↗

Remote sensing of Pu in uranyl nitrate crystals using reflectance spectroscopy and chemometrics

Remote quantification of Pu(VI) (0–5 mol%) co-crystallized with U in uranyl nitrate hexahydrate (UNH) crystals was achieved in a glove box using reflectance spectroscopy coupled with chemometric modeling. Reflectance spectra were also acquired for Pu(IV) and Np(VI) (0–5 mol%) crystallized with UNH; revealing spectral features consistent with their solution-phase analogs. Principal component analysis revealed Pu(IV/VI) and Np(VI) concentrations as the primary source of variation in the data, informing the development of a supervised partial least squares regression model for Pu(VI). The resulting calibration demonstrated robust performance, with replicate root mean square errors near 10% and quantifiable limits near 0.2 mol% Pu(VI) relative to U. The Pu(VI) remained stable in the crystalline UNH matrix for at least one week with minimal reduction to Pu(IV). Notably, Pu(VI) and Np(VI) incorporation in UNH quenched U(VI) fluorescence while Pu(IV) did not. This study presents a noninvasive, spectroscopic approach for solid-state Pu quantification, with direct implications for material accountability and nuclear nonproliferation monitoring.

Sadergaski, Luke R. [Oak Ridge National Laboratory↗

ML-Based Pebble Power Reconstruction for Pebble Bed Reactor Analysis

Pebble power reconstruction has been explored to complement the conventional homogenized modeling approach in pebble bed reactor (PBR) analysis, as detailed heterogeneous geometry calculations are computationally expensive. The random distribution of pebble fuels within the core challenges the application of conventional pin power reconstruction methods. To address this, we introduce a machine learning approach based on the transformer model, composed of encoder and decoder layers, to estimate the flux and power form functions for reconstructing individual pebble neutron fluxes and powers. The homogeneous neutron flux distribution within each spectral zone (SZ) is obtained from finite element solutions of global diffusion or transport calculations. Verification tests demonstrate that the trained transformer model accurately predicts power form functions over a range of conditions, including variations in pebble enrichment, location, type, SZ size, and burnup. In particular, verification using a three-dimensional PBR benchmark with burned pebbles shows good agreement in heterogeneous pebble power distributions between Griffin and Serpent. These results highlight the potential of applying conventional pin power reconstruction approaches to PBR cores with randomly distributed pebbles.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Theory of Photoluminescence Spectral Line Shapes of Semiconductor Nanocrystals

Single-molecule photoluminescence (PL) spectroscopy of semiconductor nanocrystals (NCs) reveals the nature of exciton–phonon interactions in NCs. Understanding the homogeneous spectral line shapes and their temperature dependence remains an open problem. Here, we develop an atomistic model to describe the PL spectrum of NCs, accounting for excitonic effects, phonon dispersion relations, and exciton–phonon couplings. We validate our model using single-NC measurements on CdSe/CdS NCs from T = 4 to 290 K, and we find that the slightly asymmetric main peak at low temperatures is comprised of a narrow zero-phonon line (ZPL) and acoustic phonon sidebands. Furthermore, we identify the specific phonon modes that give rise to the optical phonon sidebands. At temperatures above 200 K, the spectral line width shows a stronger dependence upon the temperature, which we demonstrate to be correlated with higher order exciton–phonon couplings. We also identify the line width dependence upon reorganization energy, NC core sizes, and shell thicknesses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simultaneous single-shot rotation–vibration non-equilibrium thermometry using pure rotational fs/ps CARS coherence beating

Here, we report the development of a simple and sensitive two-beam hybrid femtosecond/picosecond pure rotational coherent anti-Stokes Raman scattering (fs/ps CARS) method to simultaneously measure the rotational and vibrational temperatures of diatomic molecules. Rotation–vibration non-equilibrium plays a key role in the chemistry and thermalization in low-temperature plasmas as well as thermal loading of hypersonic vehicles. This approach uses time-domain interferences between ground state and vibrationally excited N 2 molecules to intentionally induce coherence beating that leads to apparent non-Boltzmann distributions in the pure rotational spectra. These distortions enable simultaneous inference of both the rotational and vibrational temperatures. Coherence beating effects were observed in single-shot fs/ps CARS measurements of a 75 Torr N 2 DC glow discharge and were successfully modeled for rotational and vibrational temperature extraction. We show that this method can be more sensitive than a pure rotational fs/ps CARS approach using a spectrally narrow probe pulse. Lastly, we experimentally measured the beat frequencies via Fourier transform of the time-domain response and obtained excellent agreement with the model.

42 ENGINEERING↗

Realistic inflation in no-scale U(1) R symmetric flipped SU(5)

We have realized non-minimal Higgs inflation and standard hybrid inflation in the supersymmetric flipped SU(5) model with U(1) R symmetry using the no-scale form of the Kähler potential. In non-minimal Higgs inflation the waterfall Higgs field plays the role of inflaton, and in standard hybrid inflation the gauge singlet field S is employed as an inflaton. The predictions of both models are in good agreement with the Planck 2018 data. For numerical calculations we have fixed the gauge symmetry breaking scale, M, around 2 × 10 16 GeV. In both models the inflaton field values are constrained below m P . The tensor to scalar ratio r in non-minimal inflation is of the order of 10 -3 and for standard hybrid inflation r is tiny, of order 10 -15 –10 -4 . The scalar spectral index in both cases lie within the Planck 1-σ bounds, and the running of the scalar spectral index lies in the range, -dn s /d ln k ~ 6 × 10 -4 for non-minimal model and 10 -9 –10 -3 for the standard hybrid model. A realistic scenario of reheating and non-thermal leptogenesis is employed with reheat temperature T r ~ 10 9 GeV for non-minimal model and 10 6 –10 10 GeV for standard hybrid model. The R-symmetry plays a vital role in forbidding rapid proton decay, but at the same time it also suppresses terms responsible for generating right handed neutrino masses. A realistic scenario of right handed neutrino masses is obtained by considering effective R symmetry breaking at the nonrenormalizable level with adequate suppression of rapid proton decay.

79 ASTRONOMY AND ASTROPHYSICS↗

Full radiator-perturber interaction in computer simulations of hydrogenic spectral line broadening by plasmas

Spectral line broadening by plasmas can be computed by solving the equation of motion for the dipole of the radiating system perturbed by a fluctuating potential obtained from computer simulations. Such calculations have relied on the multipole expansion of the radiator-plasma interaction often keeping only the dipole term. With increasing density, however, higher multipoles as well as plasma perturbers overlapping the bound electron wave functions are expected to become important. For hydrogenic systems, the atomic matrix elements of the full Coulomb and screened Coulomb interactions are given by analytical formulas. Using these results, a computer simulation approach that accounts for the full radiator-plasma interaction is developed. One benefit is the removal of inherent strong collision divergences in the multipole expansion approximation. Furthermore, it yields the plasma polarization shift produced by perturbers penetrating the wave function of the radiator bound electrons. Here, the model was applied to hydrogenic argon Ly-α , Ly- β , and Ly-$\gamma$ spectral lines in a dense argon plasma at free electron densities of 10 24 or 10 25 cm ⁻3 and temperature of $800$ $eV$ relevant to plasma diagnostic techniques for inertial confinement fusion implosions.

74 ATOMIC AND MOLECULAR PHYSICS↗

Flexible nonstationary spatiotemporal modeling of high-frequency monitoring data

Many physical datasets are generated by collections of instruments that make measurements at regular time intervals. For such regular monitoring data, we extend the framework of half-spectral covariance functions to the case of nonstationarity in space and time and demonstrate that this method provides a natural and tractable way to incorporate complex behaviors into a covariance model. Further, we use this method with fully time-domain computations to obtain bona fide maximum likelihood estimators—as opposed to using Whittle-type likelihood approximations, for example—that can still be computed conveniently. Additionally, we apply this method to very high-frequency Doppler LIDAR vertical wind velocity measurements, demonstrating that the model can expressively capture the extreme nonstationarity of dynamics above and below the atmospheric boundary layer and, more importantly, the interaction of the process dynamics across it.

54 ENVIRONMENTAL SCIENCES↗

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗