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,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

DIFFUSE PeV NEUTRINOS FROM GAMMA-RAY BURSTS

The IceCube Collaboration recently reported the potential detection of two cascade neutrino events in the energy range 1-10 PeV. We study the possibility that these PeV neutrinos are produced by gamma-ray bursts (GRBs), paying special attention to the contribution by untriggered GRBs that elude detection due to their low photon flux. Based on the luminosity function, rate distribution with redshift and spectral properties of GRBs, we generate, using a Monte Carlo simulation, a GRB sample that reproduces the observed fluence distribution of Fermi/GBM GRBs and an accompanying sample of untriggered GRBs simultaneously. The neutrino flux of every individual GRB is calculated in the standard internal shock scenario, so that the accumulative flux of the whole samples can be obtained. We find that the neutrino flux in PeV energies produced by untriggered GRBs is about two times higher than that produced by the triggered ones. Considering the existing IceCube limit on the neutrino flux of triggered GRBs, we find that the total flux of triggered and untriggered GRBs can reach at most a level of {approx}10{sup -9} GeV cm{sup -2} s{sup -1} sr{sup -1}, which is insufficient to account for the reported two PeV neutrinos. Possible contributions to diffuse neutrinos by low-luminosity GRBs and the earliest population of GRBs are also discussed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

A CR-hydro-NEI model of the structure and broadband emission from Tycho's supernova remnant

Tycho's supernova remnant (SNR) is well-established as a source of particle acceleration to very high energies. Constraints from numerous studies indicate that the observed γ-ray emission results primarily from hadronic processes, providing direct evidence of highly relativistic ions that have been accelerated by the SNR. Here we present an investigation of the dynamical and spectral evolution of Tycho's SNR by carrying out hydrodynamical simulations that include diffusive shock acceleration of particles in the amplified magnetic field at the forward shock of the SNR. Our simulations provide a consistent view of the shock positions, the nonthermal emission, the thermal X-ray emission from the forward shock, and the brightness profiles of the radio and X-ray emission. We compare these with the observed properties of Tycho to determine the density of the ambient material, the particle acceleration efficiency and maximum energy, the accelerated electron-to-proton ratio, and the properties of the shocked gas downstream of the expanding SNR shell. We find that evolution of a typical Type Ia supernova in a low ambient density (n {sub 0} ∼ 0.3 cm{sup –3}), with an upstream magnetic field of ∼5 μG, and with ∼16% of the SNR kinetic energy being converted into relativistic electrons and ions through diffusive shock acceleration, reproduces the observed properties of Tycho. Under such a scenario, the bulk of observed γ-ray emission at high energies is produced by π{sup 0}-decay resulting from the collisions of energetic hadrons, while inverse-Compton emission is significant at lower energies, comprising roughly half of the flux between 1 and 10 GeV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

STUDYING THE INTERSTELLAR MAGNETIC FIELD FROM ANISOTROPIES IN VELOCITY CHANNELS

Turbulence in the interstellar medium is anisotropic due to the ubiquitous magnetic fields. This anisotropy depends on the strength of the magnetic field and leaves an imprint on observations of spectral line maps. We use a grid of ideal magnetohydrodynamic simulations of driven turbulence and produce synthetic position–position–velocity maps to study the turbulence anisotropy in velocity channels of various resolutions. We found that the average structure function of velocity channels is aligned with the projection of the magnetic field on the plane of the sky. We also found that the degree of such anisotropy increases with the magnitude of the magnetic field. For thick velocity channels (low velocity resolution), the anisotropy is dominated by density, and the degree of anisotropy in these maps allows one to distinguish sub-Alfvénic and super-Alfvénic turbulence regimes, but it also depends strongly on the sonic Mach number. For thin channels (high velocity resolution), we find that the anisotropy depends less on the sonic Mach number. An important limitation of this technique is that it only gives a lower limit on the magnetic field strength because the anisotropy is related only to the magnetic field component on the plane of the sky. It can, and should, be used in combination with other techniques to estimate the magnetic field, such as the Fermi-Chandrasekhar method, anisotropies in centroids, Faraday rotation measurements, or direct line-of-sight determinations of the field from Zeeman effect observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2011

ABSORPTION CROSS SECTION OF GASEOUS ACETYLENE AT 85 K IN THE WAVELENGTH RANGE 110-155 nm

Absorption spectra and absorption cross sections of gaseous acetylene, C{sub 2}H{sub 2}, at 298 and 85 K were measured in the wavelength range 110-155 nm with a slit-jet system coupled to a synchrotron as a source of vacuum ultraviolet light. Using published spectral parameters of C{sub 2}H{sub 2}, we simulated the absorption profile for the Rydberg transition to state 4R{sub 0} in the range 124.6-125.1 nm, according to which the temperature of the jet-expanded sample at stagnation pressure 200 Torr is 85 {+-} 5 K. Our cross sections of C{sub 2}H{sub 2} are applicable for determining properties sensitive to temperature for diagnostic work on Saturn and Titan.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Extending Supernova Spectral Templates for Next-generation Space Telescope Observations

Empirical models of supernova (SN) spectral energy distributions (SEDs) are widely used for SN survey simulations and photometric classifications. The existing library of SED models has excellent optical templates, but limited, poorly constrained coverage of ultraviolet (UV) and infrared (IR) wavelengths. However, both regimes are critical for the design and operation of future SN surveys, particularly at IR wavelengths that will be accessible with the James Webb Space Telescope and the Wide-Field Infrared Survey Telescope. We create a public repository of improved empirical SED templates using a sampling of Type Ia and core-collapse (CC) photometric light curves to extend the Type Ia parameterized SALT2 model and a set of SN Ib, SN Ic, and SN II SED templates into the UV and near-IR. We apply this new repository of extrapolated SN SED models to examine how future surveys can discriminate between CC and Type Ia SNe at UV and IR wavelengths, and present an open-source software package written in Python, SNSEDextend, that enables users to generate their own extrapolated SEDs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A novel approach for calculating galaxy rotation curves using spaxel cross-correlation and iterative smoothing

ABSTRACT Precise measurements of the internal dynamics of galaxies have proven of great importance for understanding the internal dark matter distribution of galaxies. We present a novel method for measuring the line-of-sight (LOS) velocities across the face of galaxies by cross-correlation of spectral pixels (spaxels) and an iterative method of smoothing. On simulated data the method can accurately recover the input LOS velocities for different types of spectra (absorption-line dominated, emission-line dominated, and differing shapes of the continuum), and can handle stellar population radial gradients. Most important of all, it continues to provide reliable measurements of LOS velocities with reasonable uncertainties even when the spectra are very low signal-to-noise ratio (approaching ∼1), which is a challenge for traditional template-fitting approaches. We apply our method to data from a real MaNGA galaxy as a demonstration and find promising results with good precision. This novel approach can be complementary to existing methods primarily based on template fitting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Intergalactic magnetic field studies by means of γ -ray emission from GRB 190114C

The presence of delayed GeV emission after a strong transient, such as a gamma-ray burst (GRB), in the very-high energy (VHE) E > 100 GeV band can be the signature of a nonzero magnetic field in the intergalactic medium. Here we used a synchrotron self-Compton multiwavelength model to infer an analytical description of the intrinsic VHE spectrum (corrected for absorption by the extragalactic background light) of GRB 190114C to predict the light curves and spectral energy distributions of the delayed emission with Monte Carlo simulations for different intergalactic magnetic field (IGMF) configurations (strengths B = 8 × 10 -21 G, 10 -20 G, 3 × 10 -20 G and correlation length $\lambda$ > 1 Mpc), and compared them with the Fermi Large Area Telescope (Fermi LAT) limits computed for several exposure times. We found that the Fermi LAT is not sensitive enough to constrain any IGMF strengths using GRB 190114C.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Dynamic Phase Alignment in Navier-Stokes Turbulence

In Navier-Stokes turbulence, energy and helicity injected at large scales are subject to a joint direct cascade, with both quantities exhibiting a spectral scaling ∝k -5/3 . We demonstrate via direct numerical simulations that the two cascades are compatible due to the existence of a strong scale-dependent phase alignment between velocity and vorticity fluctuations, with the phase alignment angle scaling as cosα k ∝k -1 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Spectral Rear Irradiance Testing and Modeling for Degradation and Performance of Solar Fields

This work investigates how the spectrum of irradiance incident on the rear of solar modules impacts the degradation and performance of backsheets. We model the spectral irradiance incident on the rear of modules through raytrace simulations and validate with measured field data collected from a 75kW single axis-tracked bifacial test site. A generic equation to estimate relative degradation is proposed, and we show that current acceleration factors for UV damage in chambers can be sub-estimate up to 4.5% absolute from the usually assumed 10% dosage on the rear surfaces.

bifacial performance↗
DOE OSTI2016

Formation and evolution of flapping and ballooning waves in magnetospheric plasma sheet

By adopting Lembége & Pellat’s 2D plasma-sheet model, we investigate the flankward flapping motion and Sunward ballooning propagation driven by an external source (e.g., magnetic reconnection) produced initially at the sheet center. Within the ideal MHD framework, we adopt the WKB approximation to obtain the Taylor–Goldstein equation of magnetic perturbations. Fourier spectral method and Runge–Kutta method are employed in numerical simulations, respectively, under the flapping and ballooning conditions. Studies expose that the magnetic shears in the sheet are responsible for the flapping waves, while the magnetic curvature and the plasma gradient are responsible for the ballooning waves. In addition, the flapping motion has three phases in its temporal development: fast damping phase, slow recovery phase, and quasi-stabilized phase; it is also characterized by two patterns in space: propagating wave pattern and standing wave pattern. Moreover, the ballooning modes are gradually damped toward the Earth, with a wavelength in a scale size of magnetic curvature or plasma inhomogeneity, only 1–7% of the flapping one; the envelops of the ballooning waves are similar to that of the observed bursty bulk flows moving toward the Earth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Conservative Multimoment Transport along Characteristics for Discontinuous Galerkin Methods

A set of algorithms based on characteristic discontinuous Galerkin methods is presented for tracer transport on the sphere. The algorithms are designed to reduce message passing interface communication volume per unit of simulated time relative to current methods generally, and to the spectral element scheme employed by the U.S. Department of Energy's Exascale Earth System Model (E3SM) specifically. Two methods are developed to enforce discrete mass conservation when the transport schemes are coupled to a separate dynamics solver; constrained transport and Jacobian-combined transport. A communication-efficient method is introduced to enforce tracer consistency between the transport scheme and dynamics solver; this method also provides the transport scheme's shape preservation capability. A subset of the algorithms derived here is implemented in E3SM and shown to improve transport performance by a factor of 2.2 for the model's standard configuration with 40 tracers at the strong scaling limit of one element per core. Read More: https://epubs.siam.org/doi/10.1137/18M1165943

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗
DOE OSTI2024

Assessing Impacts of Waves on Hub-Height Winds off the U.S. West Coast Using Lidar Buoys and Coupled Modeling Approaches

Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.

17 WIND ENERGY↗
DOE OSTI2024

X-ray Spectral-Timing Pipeline to Investigate an Electron-Scattering Time Delay in Black Hole Accretion Disks [Slides]

The soft lag in black hole X-ray binaries (BHXRBs) refers to the time delay for the soft, thermal disk to respond to hard, variable coronal irradiation. This time lag increases from less than 1 ms to ~10 ms during the hard-to-soft state transition. Interpretations of soft lag trends appeal to changing the light-travel path via an evolving coronal height and/or inner accretion disk radius. Both interpretations neglect a time delay contribution from the reprocessing of irradiation inside the disk, where electron-scattering opacity dominates. A new theory considering a thermalization (electron scattering) time delay in the disk can plausibly produce ~10 ms time delays in the intermediate state. To further investigate the impacts of adding a time delay component from electron scattering in the disk, we are creating a spectral-timing pipeline that can analyze NICER (Neutron Star Interior Composition Explorer) X-ray observations of BHXRBs in outburst. In the near future, we will apply this spectral timing pipeline to develop a reverberation lag model from simulations that include the thermalization time delay from electron scattering in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Studying magnetohydrodynamic turbulence with synchrotron polarization dispersion

We test a new technique for studying magnetohydrodynamic turbulence suggested by Lazarian and Pogosyan, using synthetic observations of synchrotron polarization. This paper focuses on a one-point statistics, which is termed polarization frequency analysis, that is characterized by the variance of polarized emission as a function of the square of the wavelength along a single line of sight. We adopt the ratio η of the standard deviation of the line-of-sight turbulent magnetic field to the line-of-sight mean magnetic field to depict the level of turbulence. When this ratio is large (η≫1), which characterizes a region dominated by turbulent field, or small (η≲0.2), which characterizes a region dominated by the mean field, we obtain the polarization variance 〈P{sup 2}〉∝λ{sup −2} or 〈P{sup 2}〉∝λ{sup −2−2m}, respectively. At small η, i.e., in the region dominated by the mean field, we successfully recover the turbulent spectral index from the polarization variance. We find that our simulations agree well with the theoretical prediction of Lazarian and Pogosyan. With existing and upcoming data cubes from the Low-Frequency Array for Radio Astronomy (LOFAR) and the Square Kilometer Array (SKA), this new technique can be applied to study the magnetic turbulence in the Milky Way and other galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

A New MHD-assisted Stokes Inversion Technique

We present a new method of Stokes inversion of spectropolarimetric data and evaluate it by taking the example of a Sunrise/IMaX observation. An archive of synthetic Stokes profiles is obtained by the spectral synthesis of state-of-the-art magnetohydrodynamics (MHD) simulations and a realistic degradation to the level of the observed data. The definition of a merit function allows the archive to be searched for the synthetic Stokes profiles that best match the observed profiles. In contrast to traditional Stokes inversion codes, which solve the Unno–Rachkovsky equations for the polarized radiative transfer numerically and fit the Stokes profiles iteratively, the new technique provides the full set of atmospheric parameters. This gives us the ability to start an MHD simulation that takes the inversion result as an initial condition. After a relaxation process of half an hour solar time we obtain physically consistent MHD data sets with a target similar to the observation. The new MHD simulation is used to repeat the method in a second iteration, which further improves the match between observation and simulation, resulting in a factor of 2.2 lower mean χ{sup 2} value. One advantage of the new technique is that it provides the physical parameters on a geometrical height scale. It constitutes a first step toward inversions that give results consistent with the MHD equations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Shake-up and shake-off spectra in the electron capture decay of atomic $^7$Be

The most stringent laboratory-based experimental limits on the existence of sub-MeV sterile neutrinos are currently set by decay spectroscopy of radioactive $^7$Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shake-up and shake-off spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wavefunctions ab initio and compute all single and double shake processes in the electron capture decay of atomic $^7$Be. The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the $^7$Be is embedded into. The new models also show that the L/K electron capture ratio of $^7$Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

Atomic Physics (physics.atom-ph)↗