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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Spectral element methods - Algorithms and architectures

Spectral element methods are high-order weighted residual techniques for partial differential equations that combine the geometric flexibility of finite element methods with the rapid convergence of spectral techniques. Spectral element methods are described for the simulation of incompressible fluid flows, with special emphasis on implementation of spectral element techniques on medium-grained parallel processors. Two parallel architectures are considered; the first, a commercially available message-passing hypercube system; the second, a developmental reconfigurable architecture based on Geometry-Defining Processors. High parallel efficiency is obtained in hypercube spectral element computations, indicating that load balancing and communication issues can be successfully addressed by a high-order technique/medium-grained processor algorithm-architecture coupling.

Fischer, Paul↗
NASA NTRS2006

Self-Consistent Model of Magnetospheric Ring Current and Propagating Electromagnetic Ion Cyclotron Waves: Waves in Multi-Ion Magnetosphere

The further development of a self-consistent theoretical model of interacting ring current ions and electromagnetic ion cyclotron waves (Khazanov et al., 2003) is presented In order to adequately take into account wave propagation and refraction in a multi-ion magnetosphere, we explicitly include the ray tracing equations in our previous self-consistent model and use the general form of the wave kinetic equation. This is a major new feature of the present model and, to the best of our knowledge, the ray tracing equations for the first time are explicitly employed on a global magnetospheric scale in order to self-consistently simulate the spatial, temporal, and spectral evolution of the ring current and of electromagnetic ion cyclotron waves To demonstrate the effects of EMIC wave propagation and refraction on the wave energy distribution and evolution, we simulate the May 1998 storm. The main findings of our simulation can be summarized as follows. First, owing to the density gradient at the plasmapause, the net wave refraction is suppressed, and He+-mode grows preferably at the plasmapause. This result is in total agreement with previous ray tracing studies and is very clearly found in presented B field spectrograms. Second, comparison of global wave distributions with the results from another ring current model (Kozyra et al., 1997) reveals that this new model provides more intense and more highly plasmapause-organized wave distributions during the May 1998 storm period Finally, it is found that He(+)-mode energy distributions are not Gaussian distributions and most important that wave energy can occupy not only the region of generation, i.e., the region of small wave normal angles, but all wave normal angles, including those to near 90 . The latter is extremely crucial for energy transfer to thermal plasmaspheric electrons by resonant Landau damping and subsequent downward heat transport and excitation of stable auroral red arcs.

Khazanov, G. V.↗
DOE OSTI2025

Infrared spectral signatures of interfacial water at TiO 2 –electrolyte interfaces from deep potential molecular dynamics

Vibrational spectroscopy is a powerful tool for probing water at oxide–electrolyte interfaces, but its molecular interpretation can be challenging. Here, we employ deep potential long-range molecular dynamics simulations with layer-resolved spectral analysis to investigate the microscopic origins of the infrared (IR) response of water at the interface with anatase TiO 2 (101), a prototypical oxide surface. The calculated interfacial spectra exhibit characteristic modifications compared to bulk water IR spectra, including enhanced intensities, a red shifted and broadened stretching band, and a higher-frequency shoulder, in qualitative agreement with experiments. Spectral decomposition shows that these signatures originate mainly from the first interfacial water layer, dominated by surface-bound H 2 O at Ti 5C sites, with secondary contributions from the second layer. A moderate salt concentration (0.4 M NaCl) leaves both the interfacial structure and the spectra essentially unchanged, while tuning the pH strongly modulates the spectral intensity. We establish a scaling relation linking the spectral intensity to the surface water dissociation fraction and the dipole moment, both governed by interfacial electric fields. These findings provide a microscopic framework for interpreting IR spectra of oxide–electrolyte interfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1984

An analysis of the effect of biological and physical parameters of a wetlands grass biome on the spectral modeling of phytomass and primary productivity

Aircraft simulated thematic mapper data and field data were acquired in the fall and spring to analyze the relationship of spectral response and biomass for the marsh grass Spartina patens. Regression results indicate no simple relationship exists for TMS spectral response and biomass with a high R sq. However, results show a consistent relationship between spectral response and the percent live vegetation (by weight) and percent interstitial standing surface water (by area) as independent variables. It is suggested that the reflected energy of a pixel represents a mixture of surface constituents. It is recommended that alternative remote sensors be employed to account for the pixel constituents of live and dead vegetation, litter, and standing water.

Butera, M. K.↗
NASA NTRS2019

Scale–Resolving Simulations of a Fundamental Trailing–Edge Cooling Slot Using a Discontinuous–Galerkin Spectral–Element Method

The accurate prediction of turbulent mixing in high-pressure turbines that incorporate various airfoil surface-cooling strategies is becoming increasing critical to the design of modern gas turbine engines where the quest for improved efficiency is driving compressor overall pressure ratios and turbine inlet temperatures to much higher levels than ever before. In the present paper, a recently developed computational capability for accurate and efficient scale–resolving simulations of turbomachinery is extended to study the turbulent mixing mechanism of a simplified abstraction of an airfoil trailing-edge cooling slot - a plane wall jet with finite lip thickness discharging into an ambient flow. The computational capability is based on an entropy– stable, discontinuous–Galerkin approach that extends to arbitrarily high orders of spatial and temporal accuracy. The numerical results show that the present simulations capture the trends observed in the experiments. Discrepancies between the simulations and experiments are believed to be due to differences in the inflow profiles and tunnel side–wall effects. The thick lip configuration leads to a thicker wake and higher unsteadiness in the wall jet compared to the thin lip. A detailed comparison of the turbulent flowfields is presented to highlight differences arising due to lip thickness variations.

Scale Resolving↗
NASA NTRS2018

Scale-Resolving Simulations of Low-Pressure Turbine Cascades with Wall Roughness Using A Spectral-Element Method

The accurate prediction of wall-roughness effects in turbomachinery is becoming critical as turbine designers address airfoil surface quality and degradation concerns arising from the shift to advanced ceramic matrix composite (CMC) or additively-manufactured airfoils operating in higher temperature environments. In this paper, a recently developed computational capability for accurate and efficient scale-resolving simulations of turbomachinery is extended to analyze the boundary- layer separation and transition characteristics in a rough-wall low-pressure turbine (LPT) cascade. The computational capability is based on an entropy-stable discontinuous-Galerkin spectral-element approach that extends to arbitrarily high orders of spatial and temporal accuracy, and is implemented in an efficient manner for a modern high performance computer architecture. Results from the scale-resolving simulations of both smooth and rough airfoil cascades are presented and compared to previous experiments and numerical simulations. The results show that the suction surface boundary layer undergoes laminar separation, transition, and turbulent reattachment for the smooth airfoil cascade, while in the presence of roughness the separation and transition behavior of the suction surface boundary layer is substantially modified. The differences between the smooth and rough airfoil cascades are then highlighted by a detailed analysis of their respective turbulent flow fields.

Spectral-Element↗
NASA NTRS1981

The rigid shell component for superrotation in planetary atmospheres: Angular momentum budget, mechanical analog and simulation of the spin up process

An analysis of superrotation in the atmosphere of planets, with rotation axis perpendicular to the orbital plane is presented. As the atmosphere expands, Hadley cells develop producing a redistribution of mass and angular momentum. A three dimensional thermally driven zonally symmetric spectral model and Laplace transformation simulate the time evolution of a fluid leading from corotation under globally uniform heating to superrotation under globally nonuniform heating. For high viscosities the rigid shell component of atmospheric superrotation can be understood in analogy with a pirouette. During spin up angular momentum is transferred to the planet. For low iscosities, the process is reversed. A tendency toward geostrophy, combined with increase of surface pressure toward the poles (due to meridional mass transport), induces the atmosphere to subrotate temporarily at lower altitudes. Resultant viscous shear near the surface permits angular momentum to flow from the planet into the atmosphere propagating upwards to produce high altitude superrotation rates.

Mayr, H. G.↗
NASA NTRS1994

Mercury: Mid-infrared (7.3 - 13.5 microns) spectroscopic observations showing features characteristic of plagioclase

Mid-infrared spectroscopic observations of the surface of Mercury are reported for the wavelength range 7.3 to 13.5 microns. The observed spectral radiance emanated from equatorial and low latitude regions between 110-130 deg Mercurian longitude. The area is primarily an intercrater plain. The spectra show distinct and recognizable features, the principal Christiansen emission peak being the most prominent. The Christiansen feature strongly suggests the presence of plagioclase (Ca,Na)(Al,Si)AlSi2O8, (in particular labradorite: Ab(50) - Ab(30)). In addition we have studied the effects of thermal gradients to gain insight into the effects of thermal conditions on the spectral radiance of rock samples. This simulates the thermophysical effects as the rotating surface of Mercury is alternately heated and cooled. The spectral features of the samples are retained; however, the relative and absolute amplitudes vary as illustrated by laboratory reflectance and emittance spectra from quartzite.

Kozlowski, R. W. H.↗
DOE OSTI2022

Quantitative analysis of Cu XANES spectra using linear combination fitting of binary mixtures simulated by FEFF9

We report x-ray absorption near edge structure (XANES) is a powerful tool to probe the fingerprint of local structures, and when coupled with X-ray microscopy, the small spot size enables one to probe very specific regions of interest in a device or material, e.g. interfaces/bulk, different grains, good/bad electrical areas. In this work, we investigate the use of linear combination fitting (LCF) of XANES spectra for the particular case of Cu doping in CdTe. We show that the experimental data seem to be accurately represented by standards of Cu 2 Te and its substoichiometric counterpart, Cu 1.43 Te. We use Cu in CdTe as a case study to evaluate the accuracy of linear combination fitting using simulated standards, given that experimental standards for certain phases (e.g. Cu 1.43 Te) or defect structures (e.g. vacancies) cannot be readily obtained. We discuss how spectral features of the FEFF9-simulated standards, fitting ranges, and noise levels all dictate the accuracy of this type of analysis. We show that the greater the spectral difference between the two standards, the better the LCF is able to differentiate between the two structures and to tolerate experimental noise. Finally, we estimate the error of the fitted weights for different spectral features and noise levels and propose a framework to study local structures semi-quantitatively by using binary mixtures of FEFF9-simulated standards.

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

Distribution of the Spectral Response of Cells in Silicon Modules - Mechanisms and Implications

Measurements of the spectral response of cells in a PV module can provide important insight into the physics of the device and are also needed for the accurate calibration of the performance of these devices using solar simulators. While techniques for measuring the spectral response of PV cells are well established, measurement of the spectral response of all cells in a module is not straightforward since in this case one has electrical access to a string of cells and not to cells individually. Here we present a new approach for testing the spectral response of individual cells in silicon modules that allows for the fast mapping of the response at different locations of the module. We show that in some cases the spectral response of different cells in a module can vary considerably and propose an underlying mechanism for this variation. We also discuss the implications of this observation for the uncertainty of the spectral correction and of the I SC calibration using a solar simulator.

modules↗
DOE OSTI2022

Arctic sea ice albedo: Spectral composition, spatial heterogeneity, and temporal evolution observed during the MOSAiC drift

The magnitude, spectral composition, and variability of the Arctic sea ice surface albedo are key to understanding and numerically simulating Earth’s shortwave energy budget. Spectral and broadband albedos of Arctic sea ice were spatially and temporally sampled by on-ice observers along individual survey lines throughout the sunlit season (April–September, 2020) during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition. The seasonal evolution of albedo for the MOSAiC year was constructed from spatially averaged broadband albedo values for each line. Specific locations were identified as representative of individual ice surface types, including accumulated dry snow, melting snow, bare and melting ice, melting and refreezing ponded ice, and sediment-laden ice. The area-averaged seasonal progression of total albedo recorded during MOSAiC showed remarkable similarity to that recorded 22 years prior on multiyear sea ice during the Surface Heat Budget of the Arctic Ocean (SHEBA) expedition. In accord with these and other previous field efforts, the spectral albedo of relatively thick, snow-free, melting sea ice shows invariance across location, decade, and ice type. In particular, the albedo of snow-free, melting seasonal ice was indistinguishable from that of snow-free, melting second-year ice, suggesting that the highly scattering surface layer that forms on sea ice during the summer is robust and stabilizing. In contrast, the albedo of ponded ice was observed to be highly variable at visible wavelengths. Notable temporal changes in albedo were documented during melt and freeze onset, formation and deepening of melt ponds, and during melt evolution of sediment-laden ice. While model simulations show considerable agreement with the observed seasonal albedo progression, disparities suggest the need to improve how the albedo of both ponded ice and thin, melting ice are simulated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hydrostatic and Non‐Hydrostatic Baroclinic Instability in the Dynamical Core of the DOE Global Climate Model

Abstract The dynamical core of the Department of Energy global climate model is used to understand the role of non‐hydrostatic dynamics in the simulation of dry and moist baroclinic waves. To reduce computational cost, the Diabatic Acceleration and REscaling approach is adopted. Scale analysis and numerical simulations suggest that the model solution is not distorted by the change of spherical metric terms due to the change of Earth radius, neither are the inter‐scale interactions strongly altered as indicated by the analysis of spectral flux. Compared with hydrostatic simulations, the onset of baroclinic instability is delayed under non‐hydrostatic dynamics as the associated weaker vertical motions tend to increase the critical wavelength and narrow the range of waves that can be baroclinically unstable. During the development of baroclinic waves, non‐hydrostatic dynamics tends to induce vertical motions in the upper troposphere, accelerating the eastward propagation of upper‐level ridges through their impact on local vorticity tendency. These processes reduce the westward tilt of the vertical ridge axes and suppress the conversion of mean flow available potential energy to eddy kinetic energy, leading to weaker baroclinic eddies than in hydrostatic simulations. The contrasts between hydrostatic and non‐hydrostatic settings hold in supplemental experiments that vary the background flow Rossby number, the amount of water vapor content and vertical resolution. We also find that mesoscale and smaller‐scale activities can be considerably under‐represented when the vertical resolution is limited to that typically used in global climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Performance evaluation of the USGS velocity model for the San Francisco Bay Area

In this study, we evaluated the performance of the United States Geological Survey velocity model developed for the San Francisco Bay Area (SFBA), version 21.1. The evaluation was performed through high-resolution three-dimensional physics-based ground motion simulations of seven small-magnitude earthquakes (ranging from magnitude 3.8 to 4.4) that occurred on the eastern side of the San Francisco Bay. The simulations were performed in the frequency range from 0 to 5 Hz with a minimum shear-wave velocity of 250 m/s, which allowed the capture of wave propagation effects of the near-surface soft materials that characterize local basins. Based on the direct comparison of Fourier amplitude spectra between recorded and simulated ground motions for more than 250 stations, we found that the velocity model generally performs well in the frequency range of 0.2–5 Hz. The median value of the Fourier amplitude residuals was found to be near zero for all seven earthquakes. The slight over-prediction of 0.2 log-natural units at frequencies above 3 Hz in our simulations was attributed to the potentially inaccurate representation of the source radiation pattern by a double-couple point source model, and simple representation of shallow small-scale underground structural complexity in the velocity model. Maps of spectral amplitude differences between the simulated and recorded data were used to identify areas responsible for systematic ground motion over-predictions or under-predictions. For example, while some sub-domains over soft sediments show over-prediction patterns, the block east of the Hayward fault is prone to exhibit patterns of under-prediction. These maps can be used to guide future refinements of the SFBA velocity model. Since our simulation methodology allows for the decoupling of the source and wave propagation effects, the ground motion data generated by our simulations can also be used to quantify the epistemic uncertainty due to the velocity model, in empirically based ground motion estimates for the SFBA.

58 GEOSCIENCES↗
DOE OSTI2023

Large Eddy Simulation of Convective Heat Transfer in a Random Pebble Bed Using the Spectral Element Method

The development of fluoride-cooled high-temperature reactors has drastically increased the demand for an in-depth understanding of the heat transfer (HT) in packed beds cooled by liquid salts. The complex flow fields and space-dependent porosity found in a pebble bed require a detailed understanding to ensure the proper cooling of the reactor core during normal and accident conditions. As detailed experimental data are complicated to obtain for these configurations, high-fidelity simulation such as large eddy simulation and direct numerical simulation (DNS) can be employed to create a high-resolution heat transfer numerical database that can assist in addressing industrial-driven issues associated with the heat transfer behavior of fluoride-cooled high-temperature reactors. In this paper, we performed a series of large eddy simulation using computational fluid dynamics (CFD) code NekRS to investigate the heat transfer for a bed of 1741 pebbles. Further, the characteristics of the flow, such as average, rms, and time series of velocity and temperature, have been analyzed. Porous media averages have also been performed. The simulation results show a good agreement between non-conjugate heat transfer and conjugate heat transfer. The generated data will be used to benchmark heat transfer modeling methods and local maxima/minima of heat transfer parameters. It will also be used for supporting convective heat transfer quantification for Kairos Power and benchmarking lower fidelity models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1997

An optical systems analysis approach to image resampling

All types of image registration require some type of resampling, either during the registration or as a final step in the registration process. Thus the image(s) must be regridded into a spatially uniform, or angularly uniform, coordinate system with some pre-defined resolution. Frequently the ending resolution is not the resolution at which the data was observed with. The registration algorithm designer and end product user are presented with a multitude of possible resampling methods each of which modify the spatial frequency content of the data in some way. The purpose of this paper is threefold: (1) to show how an imaging system modifies the scene from an end to end optical systems analysis approach, (2) to develop a generalized resampling model, and (3) empirically apply the model to simulated radiometric scene data and tabulate the results. A Hanning windowed sinc interpolator method will be developed based upon the optical characterization of the system. It will be discussed in terms of the effects and limitations of sampling, aliasing, spectral leakage, and computational complexity. Simulated radiometric scene data will be used to demonstrate each of the algorithms. A high resolution scene will be "grown" using a fractal growth algorithm based on mid-point recursion techniques. The result scene data will be convolved with a point spread function representing the optical response. The resultant scene will be convolved with the detection systems response and subsampled to the desired resolution. The resultant data product will be subsequently resampled to the correct grid using the Hanning windowed sinc interpolator and the results and errors tabulated and discussed.

Lyon, Richard G.↗
NASA NTRS1996

Atmospheric Pre-Corrected Differential Absorption Techniques to Retrieve Columnar Water Vapor: Application to AVIRIS 91/95 Data

Water vapor is one of the main forces for weather development as well as for mesoscale air transport processes. The monitoring of water vapor is therefore an important aim in remote sensing of the atmosphere. Current operational systems for water vapor detection use primarily the emission in the thermal infrared (AVHRR, GOES, ATSR, Meteosat) or in the microwave radiation bands (DMSP). The disadvantage of current satellite systems is either a coarse spatial (horizontal) resolution ranging from one to tens of kilometers or a limited insight into the lower atmosphere. Imaging spectrometry on the other hand measures total column water vapor contents at a high spatial horizontal resolution and has therefore the potential of filling these gaps. The sensors of the AVIRIS instrument are capable of acquiring hyperspectral data in 224 bands located in the visible and near infrared at 10 nm resolution. This data includes the information on constituents of the earth's surface as well as of the atmosphere. The optical measurement of water vapor can be performed using sensor channels located in bands or lines of the absorption spectrum. The AVIRIS sensor has been used to retrieve water vapor and with less accuracy carbon dioxide, oxygen and ozone. To retrieve the water vapor amount, the so called differential absorption technique has been applied. The goal of this technique is to eliminate background factors by taking a ratio between channels within the absorption band and others besides the band. Various ratioing methods on the basis of different channels and calculation techniques were developed. The influence of a trace gas of interest on the radiance at the sensor level is usually simulated by using radiative transfer codes. In this study, the spectral transmittance and radiance are calculated by MODTRAN3 simulations with the new DISORT option. The objective of this work is to test the best performing differential absorption techniques for imaging spectrometry of tropospheric water vapor.

Schlaepfer, Daniel↗
NASA NTRS1996

Evaluation of Subgrid-Scale Models for Large Eddy Simulation of Compressible Flows

The objective of this project was to evaluate and develop subgrid-scale (SGS) turbulence models for large eddy simulations (LES) of compressible flows. During the first phase of the project results from LES using the dynamic SGS model were compared to those of direct numerical simulations (DNS) of compressible homogeneous turbulence. The second phase of the project involved implementing the dynamic SGS model in a NASA code for simulating supersonic flow over a flat-plate. The model has been successfully coded and a series of simulations has been completed. One of the major findings of the work is that numerical errors associated with the finite differencing scheme used in the code can overwhelm the SGS model and adversely affect the LES results. Attached to this overview are three submitted papers: 'Evaluation of the Dynamic Model for Simulations of Compressible Decaying Isotropic Turbulence'; 'The effect of the formulation of nonlinear terms on aliasing errors in spectral methods'; and 'Large-Eddy Simulation of a Spatially Evolving Compressible Boundary Layer Flow'.

Blaisdell, Gregory A.↗