Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spectral simulations”

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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
NASA NTRS2005

The QBO as Potential Amplifier and Conduit to Lower Altitudes of Solar Cycle Influence?

The solar cycle (SC) effect in the lower atmosphere has been linked observationally to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves (GW). Following up on a 2D study with our Numerical Spectral Model (NSM), we discuss here a 3D study in which we simulated the QBO under the influence of the SC. For a SC period of 10 years, the relative amplitude of radiative forcing is taken to vary exponentially with height, i.e., 0.2% at the surface, 2% at 50 km, 20% at 100 km and above. Applying spectral analysis to filter out and identify the SC signature, this model produces a relatively large modulation in the QBO of the lower stratospheric circulation, which is in qualitative agreement with the results obtained by Salby and Callaghan (2000) who analyzed zonal wind observations covering more than 40 years. The modeled SC modulation of the QBO extends to high latitudes where it produces temperature variations of about 1 K in the troposphere below 10 km. We report that the SC also generates in the model a hemispherically symmetric Equatorial Annual Oscillation (EAO, with 12-month period) that is largely confined to low latitudes. Under the influence of the GWs, the SC modulated EA0 propagates down into the lower stratosphere like the QBO. As is the case for the QBO, the energy of this EA0 is partially redistributed by the meridional circulation and planetary waves, presumably, to generate measurable SC signatures in the tropospheric temperature of the polar regions, which may be related to the so called Arctic Oscillation (Thompson and Wallace, 1998). The larger SC influence at higher altitudes is apparently transferred to the lower and denser regions of the atmosphere by tapping the momentum from the upward propagating GWs that drive the oscillations. The SC modulation of the QBO period could prove to be very effective in this process, as our earlier 2D study indicated. Further studies are needed, (1) to make sure that the SC effects are real and the numerical results are robust, and (2) to explore more fully the mechanism(s) that may conspire to amplify the SC effect. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with, or aid, the SC influence.

Mayr, Hans G.↗
DOE OSTI2018

Precise Time Delays from Strongly Gravitationally Lensed Type Ia Supernovae with Chromatically Microlensed Images

Time delays between the multiple images of strongly gravitationally lensed Type Ia supernovae (glSNe Ia) have the potential to deliver precise cosmological constraints, but the effects of microlensing on time delay extraction have not been studied in detail. Here we quantify the effect of microlensing on the glSN Ia yield of the Large Synoptic Survey Telescope (LSST) and the effect of microlensing on the precision and accuracy of time delays that can be extracted from LSST glSNe Ia. Microlensing has a negligible effect on the LSST glSN Ia yield, but it can be increased by a factor of ~2 over previous predictions to 930 systems using a novel photometric identification technique based on spectral template fitting. Crucially, the microlensing of glSNe Ia is achromatic until three rest-frame weeks after the explosion, making the early-time color curves microlensing-insensitive time delay indicators. By fitting simulated flux and color observations of microlensed glSNe Ia with their underlying, unlensed spectral templates, we forecast the distribution of absolute time delay error due to microlensing for LSST, which is unbiased at the sub-percent level and peaked at 1% for color curve observations in the achromatic phase, while for light-curve observations it is comparable to state-of-the-art mass modeling uncertainties (4%). About 70% of LSST glSN Ia images should be discovered during the achromatic phase, indicating that microlensing time delay uncertainties can be minimized if prompt multicolor follow-up observations are obtained. Lastly, accounting for microlensing, the 1-2 day time delay on the recently discovered glSN Ia iPTF16geu can be measured to 40% precision, limiting its cosmological utility.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Improving Spectral Resolution from Real-time Evolution for Correlated Systems

Abstract The quality of numerically simulated spectra using real-time evolution methods for strongly correlated systems is affected by both the length of simulation time and the system size, limiting resolution in both frequency and momentum. In this work, we propose a computationally cheap, linear autoregressive machine learning-based framework to extend short-time and short-distance results over a wider range. We use the proposed method to extend the lesser Green’s function for both the Hubbard model and the much more computationally challenging Hubbard-extended Holstein model. This technique significantly improves both the frequency and momentum resolution of the single-particle removal spectrum $${\mathcal{A}}(k,\omega )$$ A ( k , ω ) , allowing the observation of otherwise obscured spectral features due to electron-phonon coupling.

Tang, Ta↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2015

Enhanced laser-radiation-pressure-driven proton acceleration by moving focusing electric-fields in a foil-in-cone target

A foil-in-cone target is proposed to enhance stable laser-radiation-pressure-driven proton acceleration by avoiding the beam degradation in whole stage of acceleration. Two and three-dimensional particle-in-cell simulations demonstrate that the guiding cone can substantially improve the spectral and spatial properties of the ion beam and lead to better preservation of the beam quality. This can be attributed to the focusing effect of the radial sheath electric fields formed on the inner walls of the cone, which co-move with the accelerated foil and effectively suppress the undesirable transverse explosion of the foil. It is shown that, by using a transversely Gaussian laser pulse with intensity of ∼2.74 × 10{sup 22 }W∕cm{sup 2}, a quasi-monoenergetic proton beam with a peak energy of ∼1.5 GeV/u, density ∼10n{sub c}, and transverse size ∼1λ{sub 0} can be obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Whistler anisotropy instabilities as the source of banded chorus: Van Allen Probes observations and particle-in-cell simulations

Magnetospheric banded chorus is enhanced whistler waves with frequencies ω r < Ω e , where Ω e is the electron cyclotron frequency, and a characteristic spectral gap at ω r ≃ Ω e /2. This paper uses spacecraft observations and two-dimensional particle-in-cell simulations in a magnetized, homogeneous, collisionless plasma to test the hypothesis that banded chorus is due to local linear growth of two branches of the whistler anisotropy instability excited by two distinct, anisotropic electron components of significantly different temperatures. The electron densities and temperatures are derived from Helium, Oxygen, Proton, and Electron instrument measurements on the Van Allen Probes A satellite during a banded chorus event on 1 November 2012. The observations are consistent with a three-component electron model consisting of a cold (a few tens of eV) population, a warm (a few hundred eV) anisotropic population, and a hot (a few keV) anisotropic population. The simulations use plasma and field parameters as measured from the satellite during this event except for two numbers: the anisotropies of the warm and the hot electron components are enhanced over the measured values in order to obtain relatively rapid instability growth. The simulations show that the warm component drives the quasi-electrostatic upper band chorus and that the hot component drives the electromagnetic lower band chorus; the gap at ~Ω e /2 is a natural consequence of the growth of two whistler modes with different properties.

58 GEOSCIENCES↗
NASA NTRS1986

Energetics analysis of the observed and simulated general circulation using three-dimensional normal mode expansions

The energetics characteristics of the observed and simulated general circulation are analyzed using three-dimensional normal mode expansions. The data sets involved are the Goddard Laboratory for Atmospheres (GLA) analysis and simulation data and the Geophysical Fluid Dynamics Laboratory (GFDL) analysis data. The spectral energy properties of the Rossby and gravity modes and energy transformations are presented. Significant influences of model characteristics and the assimilation techniques are observed in the barotropic energy spectrum, particularly for the gravity mode. Energy transformations of the zonal mean field in the GLA analysis and simulation are similar, but distinctly different from that in the GFDL analysis. However, overall, the energy generation in the baroclinic mode is largely balanced by the sink in the barotropic mode. The present study may demonstrate utilities of the three-dimensional normal mode energetics in the analysis of the general circulation.

Tanaka, Hiroshi↗
NASA NTRS2013

A Review of Spectral Methods for Variable Amplitude Fatigue Prediction and New Results

A comprehensive review of the available methods for estimating fatigue damage from variable amplitude loading is presented. The dependence of fatigue damage accumulation on power spectral density (psd) is investigated for random processes relevant to real structures such as in offshore or aerospace applications. Beginning with the Rayleigh (or narrow band) approximation, attempts at improved approximations or corrections to the Rayleigh approximation are examined by comparison to rainflow analysis of time histories simulated from psd functions representative of simple theoretical and real world applications. Spectral methods investigated include corrections by Wirsching and Light, Ortiz and Chen, the Dirlik formula, and the Single-Moment method, among other more recent proposed methods. Good agreement is obtained between the spectral methods and the time-domain rainflow identification for most cases, with some limitations. Guidelines are given for using the several spectral methods to increase confidence in the damage estimate.

Larsen, Curtis E.↗
DOE OSTI2015

Pushing core-collapse supernovae to explosions in spherical symmetry. I. The model and the case of SN 1987A

We report on a method, PUSH, for artificially triggering core-collapse supernova explosions of massive stars in spherical symmetry. We explore basic explosion properties and calibrate PUSH to reproduce SN 1987A observables. Our simulations are based on the GR hydrodynamics code AGILE combined with the neutrino transport scheme isotropic diffusion source approximation for electron neutrinos and advanced spectral leakage for the heavy flavor neutrinos. To trigger explosions in the otherwise non-exploding simulations, the PUSH method increases the energy deposition in the gain region proportionally to the heavy flavor neutrino fluxes. We explore the progenitor range 18–21 M{sub ⊙}. Our studies reveal a distinction between high compactness (HC; compactness parameter ξ{sub 1.75}>0.45) and low compactness (LC; ξ{sub 1.75}<0.45) progenitor models, where LC models tend to explode earlier, with a lower explosion energy, and with a lower remnant mass. HC models are needed to obtain explosion energies around 1 Bethe, as observed for SN 1987A. However, all the models with sufficiently high explosion energy overproduce {sup 56}Ni and fallback is needed to reproduce the observed nucleosynthesis yields. {sup 57–58}Ni yields depend sensitively on the electron fraction and on the location of the mass cut with respect to the shell structure of the progenitor. We identify a progenitor and a suitable set of parameters that fit the explosion properties of SN 1987A assuming 0.1 M{sub ⊙} of fallback. We predict a neutron star with a gravitational mass of 1.50 M{sub ⊙}. We find correlations between explosion properties and the compactness of the progenitor model in the explored mass range. However, a more complete analysis will require exploring of a larger set of progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Stochastic equilibrium Raman spectroscopy (STERS)

In this manuscript, we propose a new method for cavity- and surface-enhanced Raman spectroscopy (SERS) with improved temporal resolution in the measurement of stochastic Raman spectral fluctuations. Our approach combines Fourier spectroscopy and photon correlation to decouple the integration time from the temporal resolution. Using statistical optics Monte Carlo simulations, we establish the relationship between time resolution and Raman signal strength, revealing that typical Raman spectral fluctuations, commensurate with molecular conformational dynamics, can theoretically be resolved on micro- to millisecond timescales. The method can further extract average single-molecule dynamics from small sub-ensembles, thereby potentially mitigating challenges in achieving strictly single-molecule isolation on SERS substrates.

Cobb-Bruno, Colburn [University of California, Ber↗
DOE OSTI2025

Angular-spatial hp -adaptivity for radiative transfer with discontinuous Galerkin spectral element methods

Radiative transfer is important for many science and engineering applications, and numerical simulations of radiative transfer can be challenging. For instance, the radiation field is seven-dimensional – three spatial, two angular, one wavelength, and one temporal – and often features steep gradients. Therefore, memory usage is a key issue. To reduce memory, some past work has investigated the use of adaptive mesh refinement (AMR), typically for either the spatial or angular coordinate, and typically for only h -adaptivity. Here, we propose the use of AMR for the spatial and angular coordinates together, and the use of h - and p -adaptivity together as hp -AMR for the potential for further memory savings. We implemented the proposed method for several test cases in two spatial and one angular dimension, with the discontinuous Galerkin spectral element method. These test cases featured highly anisotropic angular radiation, with or without steep spatial gradients. Our primary findings from these test cases were: (1) Angular hp -adaptivity can deliver the radiation solution with the same accuracy as, and with much less computational memory than, uniform angular h - or p -refinements, or angular h -adaptivity alone. This is most obvious when the incoming radiation is highly anisotropic, in which case the savings can be orders of magnitude. (2) Full spatial-angular hp -adaptivity is more efficient in solution representation, compared to solely spatial or solely angular -adaptivity. This is most evident when steep gradients are present in both the spatial and angular distribution. These results suggest that adaptive spatial- hp angular-refinement may perform well in large-scale seven-dimensional applications.

Adaptive refinement↗
NASA NTRS2002

Maximum Likelihood Estimation of the Broken Power Law Spectral Parameters with Detector Design Applications

The method of Maximum Likelihood (ML) is used to estimate the spectral parameters of an assumed broken power law energy spectrum from simulated detector responses. This methodology, which requires the complete specificity of all cosmic-ray detector design parameters, is shown to provide approximately unbiased, minimum variance, and normally distributed spectra information for events detected by an instrument having a wide range of commonly used detector response functions. The ML procedure, coupled with the simulated performance of a proposed space-based detector and its planned life cycle, has proved to be of significant value in the design phase of a new science instrument. The procedure helped make important trade studies in design parameters as a function of the science objectives, which is particularly important for space-based detectors where physical parameters, such as dimension and weight, impose rigorous practical limits to the design envelope. This ML methodology is then generalized to estimate broken power law spectral parameters from real cosmic-ray data sets.

Howell, Leonard W.↗
DOE OSTI2014

Stochastic modeling of the Fermi/LAT γ-ray blazar variability

We study the γ-ray variability of 13 blazars observed with the Fermi/Large Area Telescope (LAT). These blazars have the most complete light curves collected during the first four years of the Fermi sky survey. We model them with the Ornstein-Uhlenbeck (OU) process or a mixture of the OU processes. The OU process has power spectral density (PSD) proportional to 1/f {sup α} with α changing at a characteristic timescale, τ{sub 0}, from 0 (τ >> τ{sub 0}) to 2 (τ << τ{sub 0}). The PSD of the mixed OU process has two characteristic timescales and an additional intermediate region with 0 < α < 2. We show that the OU model provides a good description of the Fermi/LAT light curves of three blazars in our sample. For the first time, we constrain a characteristic γ-ray timescale of variability in two BL Lac sources, 3C 66A and PKS 2155-304 (τ{sub 0} ≅ 25 days and τ{sub 0} ≅ 43 days, respectively, in the observer's frame), which are longer than the soft X-ray timescales detected in blazars and Seyfert galaxies. We find that the mixed OU process approximates the light curves of the remaining 10 blazars better than the OU process. We derive limits on their long and short characteristic timescales, and infer that their Fermi/LAT PSD resemble power-law functions. We constrain the PSD slopes for all but one source in the sample. We find hints for sub-hour Fermi/LAT variability in four flat spectrum radio quasars. We discuss the implications of our results for theoretical models of blazar variability.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Sideband characterization and atmospheric observations with various 340 GHz heterodyne receivers

This paper describes sideband response measurements and atmospheric observations with a double sideband and two Single Sideband (SSB) receiver prototypes developed for the multi-beam limb sounder instrument stratosphere-troposphere exchange and climate monitor radiometer. We first show an advanced Fourier-Transform Spectroscopy (FTS) method for sideband response and spurious signal characterization. We then present sideband response measurements of the different prototype receivers and we compare the results of the SSB receivers with sideband measurements by injecting a continuous wave signal into the upper and lower sidebands. The receivers were integrated into a total-power radiometer and atmospheric observations were carried out. The observed spectra were compared to forward model spectra to conclude on the sideband characteristics of the different receivers. The two sideband characterization methods show a high degree of agreement for both SSB receivers with various local oscillator settings. The measured sideband response was used to correct the forward model simulations. This improves the agreement with the atmospheric observations and explains spectral features caused by an unbalanced sideband response. The FTS method also allows to quantify the influence of spurious harmonic responses of the receiver.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2018

Spectroscopic and Electronic Structure Study of ETHE1: Elucidating the Factors Influencing Sulfur Oxidation and Oxygenation in Mononuclear Nonheme Iron Enzymes

ETHE1 is a member of a growing subclass of nonheme Fe enzymes that catalyzes transformations of sulfur-containing substrates without a cofactor. ETHE1 dioxygenates glutathione persulfide (GSSH) to glutathione (GSH) and sulfite in a reaction which is similar to that of cysteine dioxygenase (CDO), but with monodentate (vs bidentate) substrate coordination and a 2-His/1-Asp (vs 3-His) ligand set. In this study, we demonstrate that GSS– binds directly to the iron active site, causing coordination unsaturation to prime the site for O 2 activation. Nitrosyl complexes without and with GSSH were generated and spectroscopically characterized as unreactive analogues for the invoked ferric superoxide intermediate. New spectral features from persulfide binding to the Fe III include the appearance of a low-energy Fe III ligand field transition, an energy shift of a NO– to Fe III CT transition, and two new GSS– to Fe III CT transitions. Time-dependent density functional theory calculations were used to simulate the experimental spectra to determine the persulfide orientation. Correlation of these spectral features with those of monodentate cysteine binding in isopenicillin N synthase (IPNS) shows that the persulfide is a poorer donor but still results in an equivalent frontier molecular orbital for reactivity. The ETHE1 persulfide dioxygenation reaction coordinate was calculated, and while the initial steps are similar to the reaction coordinate of CDO, an additional hydrolysis step is required in ETHE1 to break the S–S bond. Unlike ETHE1 and CDO, which both oxygenate sulfur, IPNS oxidizes sulfur through an initial H atom abstraction. Thus, factors that determine oxygenase vs oxidase reactivity were evaluated. In general, sulfur oxygenation is thermodynamically favored and has a lower barrier for reactivity. Furthermore, in IPNS, second-sphere residues in the active site pocket constrain the substrate, raising the barrier for sulfur oxygenation relative to oxidation via H atom abstraction.

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

The hard X-ray continuum of Cen a observed with INTEGRAL SPI

We revisit the average hard X-ray spectrum from the active galactic nucleus (AGN) of Centaurus A (Cen A) using 10 yr worth of observations with INTEGRAL SPI. This source has the highest flux observed from any AGNs in the SPI bandpass (23 keV-8 MeV). The 10 year light curve of Cen A is presented, and hardness ratios confirm that the spectral shape changes very little despite the luminosity varying by a factor of a few. Primarily, we establish the presence of a reflection component in the average spectrum by demonstrating an excess between 20 and 60 keV, from extending the spectral shape observed at low energy to the SPI regime. The excess in Chandra HETGS and INTEGRAL SPI data is well described by reflection of the dominant power-law spectrum from a neutral, optically thick atmosphere. We find that the reprocessed emission contributes 20%-25% of the 23-100 keV flux. The existence of a cutoff at tens to hundreds of kiloelectron volts remains controversial. Using simulated spectra, we demonstrate that a high energy cutoff reproduces the observed spectral properties of Cen A more readily than a simple power law. However, we also show that such a cutoff is probably underestimated when neglecting (even modest) reflection, and for Cen A would be at energies >700 keV, with a confidence of >95%. This is atypically high for thermal Comptonizing plasmas observed in AGNs, and we propose that we are in fact modeling the more gradual change in spectral shape expected of synchrotron self-Compton spectra.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

ASTROCHEMICAL CORRELATIONS IN MOLECULAR CLOUDS

We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus identify homologous species. In particular, we show that the pairs C and CO, C{sup +} and CN, and NH{sub 3} and H{sub 2}CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10'', placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

`SkyPy`: A package for modelling the Universe

SkyPy is an open-source Python package for simulating the astrophysical sky. It comprises a library of physical and empirical models across a range of observables and a command-line script to run end-to-end simulations. The library provides functions that sample realisations of sources and their associated properties from probability distributions. Simulation pipelines are constructed from these models using a YAML-based configuration syntax, while task scheduling and data dependencies are handled internally and the modular design allows users to interface with external software. SkyPy is developed and maintained by a diverse community of domain experts with a focus on software sustainability and interoperability. By fostering development, it provides a framework for correlated simulations of a range of cosmological probes including galaxy populations, large scale structure, the cosmic microwave background, supernovae and gravitational waves. Version 0.4 implements functions that model various properties of galaxies including luminosity functions, redshift distributions and optical photometry from spectral energy distribution templates. Future releases will provide additional modules, for example, to simulate populations of dark matter halos and model the galaxy-halo connection, making use of existing software packages from the astrophysics community where appropriate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗