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 937 records · Page 52

Validation of the EDGES Low-band Antenna Beam Model

The response of the antenna is a source of uncertainty in measurements with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES). We aim to validate the electromagnetic beam model of the low-band (50–100 MHz) dipole antenna with comparisons between models and against data. We find that simulations of a simplified model of the antenna over an infinite perfectly conducting ground plane are, with one exception, robust to changes in the numerical electromagnetic solver code or algorithm. For simulations of the antenna with the actual finite ground plane and realistic soil properties, we find that two out of three numerical solvers agree well. Applying our analysis pipeline to a simulated drift-scan observation from an early EDGES low-band instrument that had a 10 m × 10 m ground plane, we find residual levels after fitting and removing a five-term foreground model from the simulated data binned in local sidereal time (LST) average about 250 mK with ±40 mK variation between numerical solvers. A similar analysis of the primary 30 m × 30 m sawtooth ground plane reduced the LST-averaged residuals to about 90 mK with ±10 mK between the two viable solvers. More broadly we show that larger ground planes generally perform better than smaller ground planes. Simulated data have a power that is within 4% of real observations, a limitation of net accuracy of the sky and beam models. We observe that residual spectral structures after foreground model fits match qualitatively between simulated data and observations, suggesting that the frequency dependence of the beam is reasonably represented by the models. We find that a soil conductivity of 0.02 S m{sup -1} and relative permittivity of 3.5 yield good agreement between simulated spectra and observations. This is consistent with the soil properties reported by Sutinjo et al. for the Murchison Radio-astronomy Observatory, where EDGES is located.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Stratospheric limb radiance simulations for the HCl channel of CLAES

The earthlimb radiance of the HCl detector array of the CLAES instrument is simulated in order to specify and optimize the free spectral range for the etalon used in this spectral region and to select sampling positions which facilitate the recovery of HCl profiles. Considerations entering into the limb radiance computation of each contributing species are discussed. Simulations are then presented for tangent heights of 15 to 5 km for daytime and nighttime conditions, and their significance for HCl retrievals is addressed. The simulation results indicate that the HCl (P2) feature should be observable in the earth limb emission spectrum, and that an etalon with an FSR of 12.6/cm optimally avoids interference from OH emission lines blended with the HCl feature by adjacent transmission peaks.

Mergenthaler, J. L.↗

Conforming Chebyshev spectral collocation methods for the solution of laminar flow in a constricted channel

The numerical simulation of steady planar two-dimensional, laminar flow of an incompressible fluid through an abruptly contracting channel using spectral domain decomposition methods is described. The key features of the method are the decomposition of the flow region into a number of rectangular subregions and spectral approximations which are pointwise C(1) continuous across subregion interfaces. Spectral approximations to the solution are obtained for Reynolds numbers in the range 0 to 500. The size of the salient corner vortex decreases as the Reynolds number increases from 0 to around 45. As the Reynolds number is increased further the vortex grows slowly. A vortex is detected downstream of the contraction at a Reynolds number of around 175 that continues to grow as the Reynolds number is increased further.

Karageorghis, Andreas↗

Effect of sun and sensor geometry, canopy structure and density, and atmospheric condition on the spectral response of vegetation, with particular emphasis on across-track pointing

A computer modeling and simulation study carried out to assess the effects of various sun and sensor geometries and atmospheric conditions on the directional reflected radiance of several vegetated targets is described. Spectral responses at two wavelengths, 0.68 micron and 0.80 micron, are simulated at nine sensor zenith angles, five sensor azimuths, and nine solar zenith angles for six vegetation canopies under three atmspheric conditions. The six canopies comprise two different geometries of grass canopies at low, medium, and high leaf density. The results suggest that off-nadir viewing effects are more pronounced in the red than in the IR. However, the use of such transformations as the normalized difference index is found to reduce much of the variability seen in the bands. The magnitude of off-nadir viewing effects is found to be a function of canopy geometry.

Schnetzler, C. C.↗

Short-range order and dynamics of atoms in liquid gallium

The features of the microscopic structure, as well as one-particle and collective dynamics of liquid gallium in the temperature range from T = 313 to 1273 K, are studied on the p = 1.0 atm isobar. Detailed analysis of the data on diffraction of neutrons and X-rays, as well as the results of atomic dynamics simulation, lead to some conclusions about the structure. In particular, for preset conditions, gallium is in the equilibrium liquid phase showing no features of any stable local crystalline clusters. The pronounced asymmetry of the principle peak of the static structure factor and the characteristic “shoulder” in its right-hand part appearing at temperatures close to the melting point, which are clearly observed in the diffraction data, are due to the fact that the arrangement of the nearest neighbors of an arbitrary atom in the system is estimated statistically from the range of correlation length values and not by a single value as in the case of simple liquids. Compactly located dimers with a very short bond make a significant contribution to the statistics of nearest neighbors. The temperature dependence of the self-diffusion coefficient calculated from atomic dynamics simulation agrees well with the results obtained from experimental spectra of the incoherent scattering function. Interpolation of the temperature dependence of the self-diffusion coefficient on a logarithmic scale reveals two linear regions with a transition temperature of about 600 K. The spectra of the dynamic structure factor and spectral densities of the local current calculated by simulating the atomic dynamics indicate the existence of acoustic vibrations with longitudinal and transverse polarizations in liquid gallium, which is confirmed by experimental data on inelastic scattering of neutrons and X-rays. It is found that the vibrational density of states is completely reproduced by the generalized Debye model, which makes it possible to decompose the total vibrational motion into individual contributions associated with the formation of acoustic waves with longitudinal and transverse polarizations. Comparison of the heights of the low-frequency component and of the high-frequency peak in the spectral density of vibrational states also indicates a temperature of T ≈ 600 K, at which the diffusion type of one-particle dynamics changes to the vibrational type upon a decrease in temperature. It is demonstrated that the modified Einstein–Stokes relation can be derived using the generalized Debye model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Verification of Loop Diagnostics

Many different techniques have been used to characterize the plasma in the solar corona: density-sensitive spectral line ratios are used to infer the density, the evolution of coronal structures in different passbands is used to infer the temperature evolution, and the simultaneous intensities measured in multiple passbands are used to determine the emission measure. All these analysis techniques assume that the intensity of the structures can be isolated through background subtraction. In this paper, we use simulated observations from a 3D hydrodynamic simulation of a coronal active region to verify these diagnostics. The density and temperature from the simulation are used to generate images in several passbands and spectral lines. We identify loop structures in the simulated images and calculate the loop background. We then determine the density, temperature and emission measure distribution as a function of time from the observations and compare with the true temperature and density of the loop. We find that the overall characteristics of the temperature, density, and emission measure are recovered by the analysis methods, but the details of the true temperature and density are not. For instance, the emission measure curves calculated from the simulated observations are much broader than the true emission measure distribution, though the average temperature evolution is similar. These differences are due, in part, to inadequate background subtraction, but also indicate a limitation of the analysis methods.

Winebarger, A.↗

CGYRO

CGYRO is a multi-species spectral gyrokinetic solver with sonic rotation capability. It is an Eulerian GK solver specifically designed and optimized for collisional, electromagnetic, multiscale simulation. A key algorithmic aspect of CGYRO is the radially spectral formulation used to reduce the complicated integral gyroaveraging kernel into a multiplication in wavenumber space, but retaining the ability to treat profile variation important for edge plasmas. This makes it suitable for next-generation computational systems that require high levels of parallel concurrency.

Belli, EmilyA↗

Verification of coronal loop diagnostics using realistic three-dimensional hydrodynamic models

Many different techniques have been used to characterize the plasma in the solar corona: density-sensitive spectral line ratios are used to infer the density, the evolution of coronal structures in different passbands is used to infer the temperature evolution, and the simultaneous intensities measured in multiple passbands are used to determine the emission measure distributions. All these analysis techniques assume that the intensity of the structures can be isolated through background subtraction. In this paper, we use simulated observations from a three-dimensional hydrodynamic simulation of a coronal active region to verify these diagnostics. The density and temperature from the simulation are used to generate images in several passbands and spectral lines. We identify loop structures in the simulated images and calculate the background. We then determine the density, temperature, and emission measure distribution as a function of time from the observations and compare these with the true temperature and density of the loop. We find that the overall characteristics of the temperature, density, and emission measure are recovered by the analysis methods, but the details are not. For instance, the emission measure curves calculated from the simulated observations are much broader than the true emission measure distribution, though the average temperature evolution is similar. These differences are due, in part, to a limitation of the analysis methods, but also to inadequate background subtraction.

79 ASTRONOMY AND ASTROPHYSICS↗

Gyrokinetic particle simulation of a field reversed configuration

Gyrokinetic particle simulation of the field-reversed configuration (FRC) has been developed using the gyrokinetic toroidal code (GTC). The magnetohydrodynamic equilibrium is mapped from cylindrical coordinates to Boozer coordinates for the FRC core and scrape-off layer (SOL), respectively. A field-aligned mesh is constructed for solving self-consistent electric fields using a semi-spectral solver in a partial torus FRC geometry. This new simulation capability has been successfully verified and driftwave instability in the FRC has been studied using the gyrokinetic simulation for the first time. Initial GTC simulations find that in the FRC core, the ion-scale driftwave is stabilized by the large ion gyroradius. In the SOL, the driftwave is unstable on both ion and electron scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of one-dimensional spectral dynamic equations of the Reynolds stresses in homogeneous anisotropic turbulence: Application to split-spectrum modeling

The CTR numerical data base generated by direct simulation of homogeneous anisotropic turbulence was used to calculate all of the terms in the spectral balance equations for the turbulent Reynolds stresses. The aim in not only to test the main closure assumptions used in the split-spectrum models, but also to try to devise improved hypotheses deduced from the statistical information. Numerical simulations of turbulent flows provide a large amount of data, a thought provoking wealth of information. The main advantage of this type of comparison is that a great variety of flows can be considered, and this is necessary to test closure hypotheses. Moreover various initial conditions can be introduced in the calculation, even if they are not experimentally feasible. All the terms in the spectral equations can be calculated. The limited Reynolds numbers of the simulations and the statistical noise caused by a small sample, particularly at the large scales, causes some difficulty in the interpretation of the results, but the method of approach proved to be a powerful tool for testing and improving spectral closures.

Schiestel, R.↗

A perturbative analysis of synchrotron spectral index variation over the microwave sky

In this paper, we implement a perturbative approach, first proposed by Bouchet and Gispert, to estimate the variation of the spectral index of galactic polarized synchrotron emission, using a linear combination of simulated Stokes Q polarization maps of selected frequency bands from WMAP and Planck observations on a region of sky dominated by the synchrotron Stokes Q signal. We find that a first order perturbative analysis recovers the input spectral index map well. Along with the spectral index variation map, our method provides a fixed reference index, β-hat {sub 0s}, over the sky portion being analyzed. Using Monte-Carlo simulations, we find that 〈 β-hat {sub 0s}〉=−2.84±0.01, which matches very closely with the position of a peak at β{sub s}(p)=−2.85 of the empirical probability density function of input synchrotron indices obtained from the same sky region. For thermal dust, the mean recovered spectral index 〈 β-hat {sub d}〉=2.00±0.004 from simulations, matches very well with the spatially fixed input thermal dust spectral index β{sub d}=2.00. As accompanying results of the method, we reconstruct cosmic microwave background, thermal dust, and a synchrotron template components with fixed spectral indices over the entire sky region. We use, in our analysis, full pixel–pixel noise covariance matrices of all frequency bands, estimated from the sky region being analyzed. The perturbative technique of this work (1) can build a model with an arbitrary but sufficient degree of accuracy (and precession) as allowed by the data and (2) can produce maximum likelihood estimators for reference indices and templates asymptotically.

79 ASTRONOMY AND ASTROPHYSICS↗

High-order accurate solutions of the point kinetics equations with the spectral deferred correction method

Solving initial value problems with higher-order methods can improve the accuracy of the simulation results or the efficiency of the calculation. In this paper, we apply the spectral deferred correction (SDC) method to solve the initial value problem of the point kinetics equations (PKE). SDC is a stable, robust, and efficient high-order time-integration method capable of an arbitrary order of accuracy. For our implementation we show that it is A-stable for orders up to 8 and the order of accuracy is verified for PKE problems with a range of different reactivities. A 5.-order SDC method was then implemented to solve the exact PKE (EPKE) in the Transient Multilevel (TML) method of MPACT. The error from solutions of the EPKE is shown to be negligible. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Many-Body Level Statistics of Single-Particle Quantum Chaos

We consider a noninteracting many-fermion system populating levels of a unitary random matrix ensemble (equivalent to the q = 2 complex Sachdev-Ye-Kitaev model)—a generic model of single-particle quantum chaos. We study the corresponding many-particle level statistics by calculating the spectral form factor analytically using algebraic methods of random matrix theory, and match it with an exact numerical simulation. Despite the integrability of the theory, the many-body spectral rigidity is found to have a surprisingly rich landscape. In particular, we find a residual repulsion of distant many-body levels stemming from single-particle chaos, together with islands of level attraction. These results are encoded in an exponential ramp in the spectral form factor, which we show to be a universal feature of nonergodic many-fermion systems embedded in a chaotic medium.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Studies on spectral analysis of randomly sampled signals: Application to laser velocimetry data

Spectral analysis is very useful in determining the frequency characteristics of many turbulent flows, for example, vortex flows, tail buffeting, and other pulsating flows. It is also used for obtaining turbulence spectra from which the time and length scales associated with the turbulence structure can be estimated. These estimates, in turn, can be helpful for validation of theoretical/numerical flow turbulence models. Laser velocimetry (LV) is being extensively used in the experimental investigation of different types of flows, because of its inherent advantages; nonintrusive probing, high frequency response, no calibration requirements, etc. Typically, the output of an individual realization laser velocimeter is a set of randomly sampled velocity data. Spectral analysis of such data requires special techniques to obtain reliable estimates of correlation and power spectral density functions that describe the flow characteristics. FORTRAN codes for obtaining the autocorrelation and power spectral density estimates using the correlation-based slotting technique were developed. Extensive studies have been conducted on simulated first-order spectrum and sine signals to improve the spectral estimates. A first-order spectrum was chosen because it represents the characteristics of a typical one-dimensional turbulence spectrum. Digital prefiltering techniques, to improve the spectral estimates from randomly sampled data were applied. Studies show that the spectral estimates can be increased up to about five times the mean sampling rate.

Sree, David↗

Synthesis of Multispectral Bands from Hyperspectral Data: Validation Based on Images Acquired by AVIRIS, Hyperion, ALI, and ETM+

Spectral band synthesis is a key step in the process of creating a simulated multispectral image from hyperspectral data. In this step, narrow hyperspectral bands are combined into broader multispectral bands. Such an approach has been used quite often, but to the best of our knowledge accuracy of the band synthesis simulations has not been evaluated thus far. Therefore, the main goal of this paper is to provide validation of the spectral band synthesis algorithm used in the ART software. The next section contains a description of the algorithm and an example of its application. Using spectral responses of AVIRIS, Hyperion, ALI, and ETM+, the following section shows how the synthesized spectral bands compare with actual bands, and it presents an evaluation of the simulation accuracy based on results of MODTRAN modeling. In the final sections of the paper, simulated images are compared with data acquired by actual satellite sensors. First, a Landsat 7 ETM+ image is simulated using an AVIRIS hyperspectral data cube. Then, two datasets collected with the Hyperion instrument from the EO-1 satellite are used to simulate multispectral images from the ALI and ETM+ sensors.

Blonski, Slawomir↗

Kolmogorov and scalar spectral regimes in numerical turbulence

Velocity and passive-scalar spectra for turbulent fields generated by a forced three-dimensional simulation and Taylormicroscale Reynolds numbers up to 83 are shown to have distinct spectral regimes, including a Kolmogorov inertial subrange. Both one- and three-dimensional spectra are shown for comparison with experiment and theory, respectively. When normalized by the Kolmogorov dissipation scales velocity spectra collapse to a single curve and a high-wavenumber bulge is seen. The bulge leads to an artificially high Kolmogorov constant, but is consistent with recent measurements of the velocity spectrum in the dissipation regime and the velocity-derivative skewness. Scalar spectra, when normalized by the Oboukov-Corrsin scales, collapse to curves which depend only on Prandtl number and show a universal inertial-convective subrange, independent of Prandtl number. When normalized by the Batchelor scales, the scalar spectra show a universal dissipation regime which is independent of Prandtl numbers from 0.1 to 1.0. The time development of velocity spectra is illustrated by energy-transfer spectra in which distinct pulses propagate to high wavenumbers.

Kerr, R. M.↗

Time-resolved High Spectral Resolution Observation of 2MASSW J0746425+200032AB

Many brown dwarfs (BDs) exhibit photometric variability at levels from tenths to tens of percents. The photometric variability is related to magnetic activity or patchy cloud coverage, characteristic of BDs near the L–T transition. Time-resolved spectral monitoring of BDs provides diagnostics of cloud distribution and condensate properties. However, current time-resolved spectral studies of BDs are limited to low spectral resolution ( R ∼ 100) with the exception of the study of Luhman 16 AB at a resolution of 100,000 using the VLT+CRIRES. This work yielded the first map of BD surface inhomogeneity, highlighting the importance and unique contribution of high spectral resolution observations. Here, we report on the time-resolved high spectral resolution observations of a nearby BD binary, 2MASSW J0746425+200032AB. We find no coherent spectral variability that is modulated with rotation. Based on simulations, we conclude that the coverage of a single spot on 2MASSW J0746425+200032AB is smaller than 1% or 6.25% if spot contrast is 50% or 80% of its surrounding flux, respectively. Future high spectral resolution observations aided by adaptive optics systems can put tighter constraints on the spectral variability of 2MASSW J0746425+200032AB and other nearby BDs.

79 ASTRONOMY AND ASTROPHYSICS↗