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 199 records · Page 11

Solar Cycle Variations in Mesospheric Carbon Monoxide

As an extension of Lee et al. (2013), solar cycle variation of carbon monoxide (CO) is analyzed with MLS observation, which covers more than thirteen years (2004-2017) including maximum of solar cycle 24. Being produced primarily by the carbon dioxide (CO2) photolysis in the lower thermosphere, the variations of the mesospheric CO concentration are largely driven by the solar cycle modulated ultraviolet (UV) variation. This solar signal extends down to the lower altitudes by the dynamical descent in the winter polar vortex, showing a time lag that is consistent with the average descent velocity. To characterize a global distribution of the solar impact, MLS CO is correlated with the SORCE measured total solar irradiance (TSI) and UV. As high as 0.8 in most of the polar mesosphere, the linear correlation coefficients between CO and UV/TSI are more robust than those found in the previous work. The photochemical contribution explains most (68%) of the total variance of CO while the dynamical contribution accounts for 21% of the total variance at upper mesosphere. The photochemistry driven CO anomaly signal is extended in the tropics by vertical mixing. The solar cycle signal in CO is further examined with the Whole Atmosphere Community Climate Model (WACCM) 3.5 simulation by implementing two different modeled Spectral Solar Irradiances (SSIs): SRPM 2012 and NRLSSI. The model simulations underestimate the mean CO amount and solar cycle variations of CO, by a factor of 3, compared to those obtained from MLS observation. Different inputs of the solar spectrum have small impacts on CO variation.

MLS↗

Broad-band Spectroscopy of Hercules X-1 with Suzaku

Hercules X-1 was observed with Suzaku in the main-on state from 2005 to 2010. The 0.4–100 keV wide-band spectra obtained in four observations showed a broad hump around 4–9 keV in addition to narrow Fe lines at 6.4 and 6.7 keV. The hump was seen in all the four observations regardless of the selection of the continuum models. Thus it is considered a stable and intrinsic spectral feature in Her X-1. The broad hump lacked a sharp structure like an absorption edge. Thus it was represented by two different spectral models: an ionized partial covering or an additional broad line at 6.5 keV. The former required a persistently existing ionized absorber, whose origin was unclear. In the latter case, the Gaussian fitting of the 6.5-keV line needs a large width of σ = 1.0–1.5 keV and a large equivalent width of 400–900 eV. If the broad line originates from Fe fluorescence of accreting matter, its large width may be explained by the Doppler broadening in the accretion flow. However, the large equivalent width may be inconsistent with a simple accretion geometry.

Asami, Fumi↗

FISM 2.0: Improved Spectral Range, Resolution, and Accuracy

The Flare Irradiance Spectral Model (FISM) was first released in 2005 to provide accurate estimates of the solar VUV (0.1-190 nm) irradiance to the Space Weather community. This model was based on TIMED SEE as well as UARS and SORCE SOLSTICE measurements, and was the first model to include a 60 second temporal variation to estimate the variations due to solar flares. Along with flares, FISM also estimates the tradition solar cycle and solar rotational variations over months and decades back to 1947. This model has been highly successful in providing driving inputs to study the affect of solar irradiance variations on the Earth's ionosphere and thermosphere, lunar dust charging, as well as the Martian ionosphere. The second version of FISM, FISM2, is currently being updated to be based on the more accurate SDO/EVE data, which will provide much more accurate estimations in the 0.1-105 nm range, as well as extending the 'daily' model variation up to 300 nm based on the SOLSTICE measurements. with the spectral resolution of SDO/EVE along with SOLSTICE and the TIMED and SORCE XPS 'model' products, the entire range from 0.1-300 nm will also be available at 0.1 nm, allowing FISM2 to be improved a similar 0.1nm spectral bins. FISM also will have a TSI component that will estimate the total radiated energy during flares based on the few TSI flares observed to date. Presented here will be initial results of the FISM2 modeling efforts, as well as some challenges that will need to be overcome in order for FISM2 to accurately model the solar variations on time scales of seconds to decades.

Chamberlin, Phillip C.↗

A Decentralized Approach for Modeling Organized Convection Based on Thermal Populations on Microgrids

Abstract In this study, a spectral model for convective transport is coupled to a thermal population model on a two‐dimensional horizontal “microgrid,” covering the typical gridbox size of general circulation models. The goal is to explore new ways of representing impacts of spatial organization in cumulus cloud fields. The thermals are considered the smallest building block of convection, with thermal life cycle and movement represented through binomial functions. Thermals interact through two simple rules, reflecting pulsating growth and environmental deformation. Long‐lived thermal clusters thus form on the microgrid, exhibiting scale growth and spacing that represent simple forms of spatial organization and memory. Size distributions of cluster number are diagnosed from the microgrid through an online clustering algorithm, and provided as input to a spectral multiplume eddy‐diffusivity mass flux scheme. This yields a decentralized transport system, in that the thermal clusters acting as independent but interacting nodes that carry information about spatial structure. The main objectives of this study are (a) to seek proof of concept of this approach, and (b) to gain insight into impacts of spatial organization on convective transport. Single‐column model experiments demonstrate satisfactory skill in reproducing two observed cases of continental shallow convection. Metrics expressing self‐organization and spatial organization match well with large‐eddy simulation results. We find that in this coupled system, spatial organization impacts convective transport primarily through the scale break in the size distribution of cluster number. The rooting of saturated plumes in the subcloud mixed layer plays a key role in this process.

54 ENVIRONMENTAL SCIENCES↗

Integrated and spectral energetics of the GLAS general circulation model

Integrated and spectral error energetics of the Goddard Laboratory for Atmospheric Sciences (GLAS) general circulation model are compared with observations for periods in January 1975, 1976, and 1977. For two cases the model shows significant skill in predicting integrated energetics quantities out to two weeks, and for all three cases, the integrated monthly mean energetics show qualitative improvements over previous versions of the model in eddy kinetic energy and barotropic conversions. Fundamental difficulties remain with leakage of energy to the stratospheric level. General circulation model spectral energetics predictions are compared with the corresponding observational spectra on a day by day basis. Eddy kinetic energy can be correct while significant errors occur in the kinetic energy of wavenumber three. Single wavenumber dominance in eddy kinetic energy and the correlation of spectral kinetic and potential energy are demonstrated.

Tenenbaum, J.↗

Technical report series on global modeling and data assimilation. Volume 2: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere

High order horizontal diffusion of the form K Delta(exp 2m) is widely used in spectral models as a means of preventing energy accumulation at the shortest resolved scales. In the spectral context, an implicit formation of such diffusion is trivial to implement. The present note describes an efficient method of implementing implicit high order diffusion in global finite difference models. The method expresses the high order diffusion equation as a sequence of equations involving Delta(exp 2). The solution is obtained by combining fast Fourier transforms in longitude with a finite difference solver for the second order ordinary differential equation in latitude. The implicit diffusion routine is suitable for use in any finite difference global model that uses a regular latitude/longitude grid. The absence of a restriction on the timestep makes it particularly suitable for use in semi-Lagrangian models. The scale selectivity of the high order diffusion gives it an advantage over the uncentering method that has been used to control computational noise in two-time-level semi-Lagrangian models.

Max J. Suarez↗

A low dimensional dynamical system for the wall layer

Low dimensional dynamical systems which model a fully developed turbulent wall layer were derived.The model is based on the optimally fast convergent proper orthogonal decomposition, or Karhunen-Loeve expansion. This decomposition provides a set of eigenfunctions which are derived from the autocorrelation tensor at zero time lag. Via Galerkin projection, low dimensional sets of ordinary differential equations in time, for the coefficients of the expansion, were derived from the Navier-Stokes equations. The energy loss to the unresolved modes was modeled by an eddy viscosity representation, analogous to Heisenberg's spectral model. A set of eigenfunctions and eigenvalues were obtained from direct numerical simulation of a plane channel at a Reynolds number of 6600, based on the mean centerline velocity and the channel width flow and compared with previous work done by Herzog. Using the new eigenvalues and eigenfunctions, a new ten dimensional set of ordinary differential equations were derived using five non-zero cross-stream Fourier modes with a periodic length of 377 wall units. The dynamical system was integrated for a range of the eddy viscosity prameter alpha. This work is encouraging.

Aubry, N.↗

Using Beta Cen to Probe the Photosphere-Wind Connection to B Stars

The goal of the project was to probe the connection between the photospheric pulsations and the X-ray production in the stellar winds of B stars. The B1 III star Beta Centauri is a Beta Cephei variable with a high X-ray count rate, and was observed by XMM. Analysis of XMM X-ray spectra showed characteristic plasma temperatures of 0.1, 0.2 and 0.6 keV. Analysis of He-like forbidden-to-intercombination lines was performed to determine the dominant locations of the X-ray emission. At Prism, work focused on developing and upgrading spectral models to calculate forbidden-to-intercombination line ratios using more detailed atomic models than the simple models typically used in the analysis of X-ray spectral lines, thereby allowing us to assess the effects of a broader variety of populating and depopulating channels for energy levels contributing to the observed X-ray line emission. Based on Ne IX lines, the source location was determined to be no farther than 4 stellar radii. Examination of the time variability of the XMM data for Beta Cen showed no significant correlation with the known optical pulsation period for the star, and also did not find any significant periods of variation at X-ray wavelengths.

MacFarlane, Joseph J.↗

Model calculations of spectral transmission for the CLAES etalons

This paper describes models for calculating spectral transmission for the Cryogenic-Limb-Array-Etalon-Spectrometer (CLAES) etalons. These models involve a convolution of the Airy function for a given thickness with the distribution of surface thicknesses, the effect of absorption in the substrate, and the field of view broadening as a function of etalon tilt angle. A comparison of model calculations with experimental transmission data for CLAES etalons centered at 3.52, 5.72, 8.0, and 11.86 microns showed that these models are able to provide a good description of the CLAES etalons.

James, T. C.↗

Pulsar Emission Geometry and Accelerating Field Strength

The high-quality Fermi LAT observations of gamma-ray pulsars have opened a new window to understanding the generation mechanisms of high-energy emission from these systems, The high statistics allow for careful modeling of the light curve features as well as for phase resolved spectral modeling. We modeled the LAT light curves of the Vela and CTA I pulsars with simulated high-energy light curves generated from geometrical representations of the outer gap and slot gap emission models. within the vacuum retarded dipole and force-free fields. A Markov Chain Monte Carlo maximum likelihood method was used to explore the phase space of the magnetic inclination angle, viewing angle. maximum emission radius, and gap width. We also used the measured spectral cutoff energies to estimate the accelerating parallel electric field dependence on radius. under the assumptions that the high-energy emission is dominated by curvature radiation and the geometry (radius of emission and minimum radius of curvature of the magnetic field lines) is determined by the best fitting light curves for each model. We find that light curves from the vacuum field more closely match the observed light curves and multiwavelength constraints, and that the calculated parallel electric field can place additional constraints on the emission geometry

DeCesar, Megan E.↗

ROSAT PSPC observations of two X-ray-faint early-type galaxies: NGC 4365 and NGC 4382

We present the results of ROSAT Positive Sensitive Proportional Counter (PSPC) observations of the two early-type galaxies NGC 4365 and NGC 4382. These galaxies are among those observed with Einstein to have the lowest X-ray to optical flux ratios of early-type galaxies. The PSCP data show that for radii r greater than 50 arcsec the radial distributions of the X-ray surface brightness are consistent with the optical distributions of King (1978). We also find that these galaxies have X-ray spectra significantly different from those observed in X-ray-bright ellipticals, with a relative excess of counts detected in the softest spectral channels. This confirms earlier Einstein results. The characteristics of the ROSAT PSPC do not allow us to discriminate between possible spectral models. If we adopt a two-component thermal model on the grounds of physical plausibility, we find that the spectral data can be fitted with a very soft optically thin component, with kT approximately 0.2 keV, and a hard component with kT greater than (1.0-1.5) keV. The hard component has a luminosity consistent with that expected from the integrated emission of a population of low mass-X-ray binaries in these galaxies; the nature of the very soft component is more speculative. Candidates include the coronal emission of late-type stars, supersoft X-ray sources, RS CVn, and perhaps a hot Interstellar Medium (ISM). Alternatively, the spectal data may be fitted with a 0.6-1 keV bremsstrahlung spectrum (expontential plus Gaunt), and may suggest the presence of a totally new population of X-ray sources.

Fabbiano, G.↗

Low-Energy Study of Gamma-Ray Bursts Having Spectral Line Features

Gamma-ray bursts (GRBs) are energetic, short-duration emissions of gamma-rays. The Burst and Transient Source Experiment (BATSE) that was onboard NASA s Compton Gamma-Ray Observatory has done much to advance our understanding of GRBs. Perhaps foremost is to establish that GRBs originate from astronomical sources that exist well beyond our galaxy. Another area in which BATSE has been instrumental is to provide high-resolution data that can be used in spectral studies. Before BATSE, there were many reports of GRB spectra containing what appeared to be spectral absorption lines, whereas Briggs, after an extensive computer search of 117 bright BATSE GRBs, reported finding only one case that might be an absorption line and ten cases that might be emission lines. None of the eleven BATSE cases were definitively identified as spectral lines, and Briggs indicated reasons as to why the pre-BATSE reports should not be taken as conclusive. It remains an open question as to what these spectral-like features are, or if they are even real. The purpose of this work is, for the subset of the eleven BATSE GRBs for which low-energy data are available from two BATSE's Spectroscopy Detectors (SDs), to include these data in the spectral analysis. Such a study will provide additional constraints on the model spectral functions to better ascertain the reality of the line features. The spectral analysis program used was RMFIT. Of the six GRBs that met the selection criteria, the analysis was performed on only three of them due to a lack of time.

Pangia, Michael J.↗

Multi-Band Atmospheric Correction Algorithm for Ocean Color Retrievals

NASA's current Atmospheric Correction (AC) algorithm for ocean color utilizes two bands and their ratio in the Near Infrared (NIR) to estimate aerosol reflectance and aerosol type. The algorithm then extrapolates the spectral dependence of aerosol reflectance to the visible wavelengths based on modeled spectral dependence of the identified aerosol type. Future advanced ocean color sensors, such as the Ocean Color Instrument (OCI) that will be carried on the Plankton, Aerosol, Cloud, and ocean Ecosystem (PACE) satellite, will be capable of measuring the hyperspectral radiance from 340 to 890 nm at 5-nm spectral resolution and at 7 discrete Short-wave Infrared (SWIR) channels: 940, 1038,1250, 1378, 1615, 2130, and 2260 nm. To optimally employ this unprecedented instrument capability, we propose an improved AC algorithm that utilizes all atmospheric-window channels in the NIR to SWIR spectral range to reduce the uncertainty in the AC process. A theoretical uncertainty analysis of this, namely Multi-Band AC (MBAC), indicates that the algorithm can reduce the uncertainty in remote sensing reflectance (Rrs) retrievals of the ocean caused by sensor random noise. Furthermore, in optically complex waters, where the NIR signal is affected by contributions from highly-reflective turbid waters, the MBAC algorithm can be adaptively weighted to the strongly-absorbing SWIR channels to enable improved ocean color retrievals in coastal waters. We provide here a description of the algorithm and demonstrate the improved performance in ocean color retrievals, relative to the current NASA standard AC algorithm, through comparison with field measurements and assessment of propagated uncertainties in applying the MBAC algorithm to MODIS and simulated PACE OCI data.

PACE↗

Dust and gas absorption in the High Mass X-ray Binary IGR J16318−4848

Context. With an absorption column density on the order of 10(exp 24)per sq. cm, IGR J16318-4848 is one of the most extreme cases of a highly obscured High Mass X-ray Binary. Besides the overall continuum absorption, the source spectrum exhibits a strong iron and nickel fluorescence line complex at 6.4 keV. Previous empirical modeling of these features and comparison with radiative transfer simulations raised questions about the structure and covering fraction of the absorber and the profile of the fluorescence lines. Aims. We aim at a self-consistent description of the continuum absorption, the absorption edges and the fluorescent lines to constrain properties of the absorbing material, such as ionization structure and geometry. We further investigate the effects of dust absorption on the observed spectra and the possibility of fluorescent emission from dust grains. Methods. We use XMM-Newton and NuSTAR spectra to first constrain empirically the incident continuum and fluorescence lines. Next, we use XSTAR to construct a customized photoionization model where we vary the ionization parameter, column density, and covering fraction. In a third step, we model the absorption and fluorescence in a dusty olivine absorber and employ both, a simple analytical model for the fluorescent line emission and Monte Carlo radiative transfer, spectral shapes and line fluxes that are very close to the data are generated. Results. Our empirical spectral modeling is in agreement with previous works. Our second model, the single gas absorber does not describe the observational data. In particular, irrespective of the ionization state or column density of the absorber, a much higher covering fraction than previously estimated is needed to produce the strong fluorescence lines and the large continuum absorption. A dusty, spherical absorber (modeled as consisting of olivine dust, although the nature of dust cannot be constrained) is well able to produce the observed continuum absorption and edges. Conclusions. A dense, dusty absorber in the direct vicinity of the source consisting of dust offers a consistent description of both the strong continuum absorption and the strong emission features in the X-ray spectrum of IGR J16318􀀀4848. In particular, for low optical depth of individual grains the dust will contribute significantly to the fluorescent emission, which is the case for typical densities and grain size distribution models.

Ralf Ballhausen↗

Discovery of Soft X-ray Pulsations from PSR J1231-1411 using NICER

We report the discovery of soft X-ray pulsations from the nearby millisecond pulsar PSR J1231−1411 using the Neutron Star Interior Composition Explorer (NICER). The pulsed emission is characterized by a broad and asymmetric main pulse and a much fainter secondary interpulse, with a total pulsed count rate of 0.055 c s−1 in the 0.35–1.5 keV band. We analyzed Fermi Large Area Telescope (LAT) data to update the pulse timing model covering 10 yr of data and used that model to coherently combine NICER data over a year of observations. Spectral modeling suggests that the flux is dominated by thermal emission from a hot spot (or spots) on the neutron star surface. The phase relationship between the X-ray pulse and the radio and γ rays provides insight into the geometry of the system.

gamma rays↗

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

Fitting Strategies of Accretion Column Models and Application to the Broadband Spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT data-set of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies

Philipp Thalhammer↗

Reconstructing the Solar VUV Irradiance over the Past 60 Years

The Flare Irradiance Spectral Model (FISM) is an empirical model of the solar irradiance spectrum from 0.1 to 190 nm at I nm spectral resolution and on a 1-minute time cadence. The goal of FISM is to provide accurate solar spectral irradiances over the vacuum ultraviolet (VUV: 0-200 nm) range as input for ionospheric and thermospheric models, as well as climate studies over 60 years. A brief overview of the FISM model will be given, and also discussed is how the Solar Dynamics Observatory (SDO) EUV Variability Experiment (EVE) will contribute to improving FISM and its accuracies. Results will also be shown quantifying the VUV contributions to the total flare energy budget, and more importantly discuss the increased errors associated by not including flares in the solar energy input to the Earth's system. Concluding the talk will be a discussion of the proxies, and their associated uncertainties, used for solar spectral reconstructions prior to 1947 going back hundreds of years.

Chamberlin, Phillip↗