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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Small-scale properties from exascale computations of turbulence on a $\mathbf{32\,768^3}$ periodic cube

To study the physics of small-scale properties of homogeneous isotropic turbulence at increasingly high Reynolds numbers, direct numerical simulation results have been obtained for forced isotropic turbulence at Taylor-scale Reynolds number R λ = 2500 on a 32 768 3 three-dimensional periodic domain using a GPU pseudo-spectral code on a 1.1 exaflop GPU supercomputer (Frontier). These simulations employ the multi-resolution independent simulation (MRIS) technique (Yeung & Ravikumar 2020, Phys. Rev. Fluids, vol. 5, 110517) where ensemble averaging is performed over multiple short segments initiated from velocity fields at modest resolution, and subsequently taken to higher resolution in both space and time. Reynolds numbers are increased by reducing the viscosity with the large-scale forcing parameters unchanged. Although MRIS segments at the highest resolution for each Reynolds number last for only a few Kolmogorov time scales, small-scale physics in the dissipation range is well captured – for instance, in the probability density functions and higher moments of the dissipation rate and enstrophy density, which appear to show monotonic trends persisting well beyond the Reynolds number range in prior works in the literature. Attainment of range of length and time scales consistent with classical scaling also reinforces the potential utility of the present high-resolution data for studies of short-time-scale turbulence physics at high Reynolds numbers where full-length simulations spanning many large-eddy time scales are still not accessible. A single snapshot of the 32 768 3 data is publicly available for further analyses via the Johns Hopkins Turbulence Database.

intermittency↗
NASA NTRS1985

The utilization of satellite data and dynamics in understanding and predicting global weather phenomena

A two layer spectral quasi-geostrophic model is used to simulate the effects of topography on the equilibria, the stability, and the long term evaluation of incipient unstable waves. The flow is forced by latitudinally dependent radiational heating. The nature of the form drag instability of high index equilibria is investigated. The proximity of the equilibrium shear to a resonant value is essential for the instability, provided the equilibrium occurs at a slightly stronger shear than resonance. The properties of the steady Hadley and Rossby required for a thermally forced rotating fluid on a sphere are further explained. An objective parameterization technique is developed for general nonlinear hydrodynamical systems. The typical structure is one in which the rates of change of the dependent variables depend on homogeneous quadratic and linear forms, as well as on inhomogeneous forcing terms. Also documented is a steady, axisymmetric model of the general circulation developed as a basis for climate stability studies. The model includes the effects of heating, rotation, and internal friction, but neglects topography. Included is further research on cloud street phenomena. Orientation angles and horizontal wavelengths of boundary layer rolls and cloud streets are determined from an analysis of a truncated spectral model of three dimensional shallow moist Boussinesq convection in a shearing environment is further explained. Relatively broadly spaced roll clouds have orientations for which the Fourier component of the roll perpendicular shear is nearly zero, but the second corresponds to narrowly spaced rolls having orientations for which the Fourier coefficients of both the perpendicular and the parallel components of the shear are nearly equal.

Shirer, H. N.↗
NASA NTRS2001

An Analysis of Simulated and Observed Global Mean Near-Surface Air Temperature Anomalies from 1979 to 1999: Trends and Attribution of Causes

The 1979 - 1999 response of the climate system to variations in solar spectral irradiance is estimated by comparing the global averaged surface temperature anomalies simulated by a 2D (two dimensional) energy balance climate model to observed temperature anomalies. We perform a multiple regression of southern oscillation index and the individual model responses to solar irradiance variations, stratospheric and tropospheric aerosol loading, stratospheric ozone trends, and greenhouse gases onto each of five near-surface temperature anomaly data sets. We estimate the observed difference in global mean near surface air temperature attributable to the solar irradiance difference between solar maximum and solar minimum to be between 0.06 and 0.11 K, and that 1.1 - 3.8% of the total variance in monthly mean near-surface air temperature data is attributable to nations in solar spectral irradiance. For the five temperature data sets used in our analysis, the trends in raw monthly mean temperature anomaly data have a large range, spanning a factor of 3 from 0.06 to 0.17 K/decade. However. our analysis suggests that trends in monthly temperature anomalies attributable to the combination of greenhouse gas, stratospheric ozone, and tropospheric sulfate aerosol variations are much more consistent among data sets, ranging from 0.16 to 0.24 K/decade. Our model results suggest that roughly half of the warming from greenhouse gases is cancelled by the cooling from changes in stratospheric ozone. Tropospheric sulfate aerosol loading in the present day atmospheric contributes significantly to the net radiative forcing of the present day climate system. However, because the change in magnitude and latitudinal distribution of tropospheric sulfate aerosol has been small over the past 20 years, the change in the direct radiative forcing attributable to changes in aerosol loading over this time is also small.

MacKay, R. M.↗
NASA NTRS1993

Simulations of Deep Galaxy Fields I: Monte-Carlo Simulations of Optical and Near Infrared Counts

Monte-Carlo simulations of three dimensional galaxy distributions are performed, following the prescription of Chokishi & Wright (1988), to study the photometric properties of evolving galaxy populations in the optical and near infrared bands to high redshifts. Galaxies are spatially distributed according to the spatial two point correlation function in a range of cluster to field environments andin volume elements appropriate to the cosmological model under consideration. Details of individual galaxy properties, including luminosities, morphologies, disk-to-bulge ratios, and size distributions are simulated to match local observed galaxy properties. Galaxies have evolving gas fraction, metallicity and spectral energy distributions that include both stellar emission and internal dust absorption and re-emission. The simulations result in noiseless two dimensional galaxy distributions on the sky that canbe compared to the observed deep images in the blue and near-infrared bands...

faint↗
DOE OSTI2018

Compression of femtosecond ytterbium fibre laser pulses using nonlinear processes in silica fibre

The temporal compression of 100-fs nanojoule laser pulses using nonlinear processes in silica fibre has been studied by numerical simulation and experimentally. The simulation results demonstrate that, in the case of spectral broadening in a normal dispersion fibre, followed by the temporal compression of the pulse in a compressor with a negative second-order dispersion, up to 90 % of the compressed pulse energy can be concentrated in the main peak, whose duration is a factor of 8 shorter than the initial laser pulse duration. In our experiments, 100-fs ytterbium fibre laser pulses have been compressed to pulses with an ∼13-fs main peak accounting for 87 % of the compressed pulse energy. At an average ytterbium laser output power of 4.2 W, the average power at the compressor output has been 3.1 W. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Verification of the Two-Fluid Model in the Spectral Element Code Nek-2P for Adiabatic Flow Conditions

The simulation of multiphase flow is a difficult physical, computational, and numerical problem. Sophisticated tools are necessary to provide results encompassing a wide range of physical phenomena. The development of these tools must be supported by rigorous analysis and ongoing efforts of verification and validation to ensure the accuracy of simulations. The Nek-2P multiphase code is being developed at Argonne with the general goal of simulating a wide range of multiphase flows and the primary goal of simulating boiling flows in LWR fuel assemblies. As part of this, there is an ongoing effort to verify the code and validate the underlying physical models. The objective of this work is to present recent results in the verification of Nek-2P. As a starting point, this work is focused on adiabatic flows for the sake of simplicity with the immediate next planned step to expand the effort to include heat and mass transfer. Details of the verification of the implemented mass fraction transport, momentum conservation, and pressure equations in Nek-2P are presented. Results of simulations are compared to analytic solutions when possible, however due to the complexity of multiphase flows, these are seldom available. To overcome this difficulty, an independent model is used to provide solutions for comparison.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2014

Simulating signatures of two-dimensional electronic spectra of the Fenna-Matthews-Olson complex: By using a numerical path integral

A framework for simulating electronic spectra from photon-echo experiments is constructed by using a numerical path integral technique. This method is non-Markovian and nonperturbative and, more importantly, is not limited by a fixed form of the spectral density functions of the environment. Next, a two-dimensional (2D) third-order electronic spectrum of a dimer system is simulated. The spectrum is in agreement with the experimental and theoretical results previously reported [for example, M. Khalil, N. Demirdöven, and A. Tokmakoff, Phys. Rev. Lett. 90, 047401 (2003)]. Finally, a 2D third-order electronic spectrum of the Fenna-Matthews-Olson (FMO) complex is simulated by using the Debye, Ohmic, and Adolphs and Renger spectral density functions. It is shown that this method can clearly produce the spectral signatures of the FMO complex by using only the Adolphs and Renger spectral density function. Plots of the evolution of the diagonal and cross-peaks show that they are oscillating with the population time.

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

Noise Reduction for Solar-Induced Fluorescence Retrievals Using Machine Learning and Principal Component Analysis: Simulations and Applications to GOME-2 Satellite Retrievals

We use a spectral-based approach that employs principal component analysis along with a relatively shallow artificial neural network (NN) to substantially reduce noise and other artifacts in terrestrial chlorophyll solar-induced fluorescence (SIF) retrievals. SIF is a very small emission at red and far-red wavelengths that is difficult to measure and is highly sensitive to random errors and systematic artifacts. Our approach relies upon an assumption that a trained NN can effectively reconstruct the total SIF signal from a relatively small number of leading principal components of the satellite-observed far-red radiance spectra without using information from the trailing modes that contain most of the random errors. We test the approach with simulated reflectance spectra produced with a full atmospheric and surface radiative transfer model using different observing and geophysical parameters and various noise levels. Resulting noisy and noise-reduced retrieved SIF values are compared with true values to assess performance. We then apply our noise reduction approach to SIF derived from two different satellite spectrometers. For evaluation, since the truth in this case is unknown, we compare SIF retrievals from two independent sensors with each other. We also compare the noise-reduced SIF temporal variations with those from an independent gross primary product (GPP) product that should display similar variations. Results show that our noise reduction approach improves capture of SIF seasonal and inter-annual variability. Our approach should be applicable to many noisy data products derived from spectral measurements. Our methodology does not replace the original retrieval algorithms; rather the original noisy retrievals are needed as the target for the NN training process.

noise reduction, solar-induced fluorescence, SIF, ↗
NASA NTRS2004

Simulating Remote Sensing Systems

The Application Research Toolbox (ART) is a collection of computer programs that implement algorithms and mathematical models for simulating remote sensing systems. The ART is intended to be especially useful for performing design-tradeoff studies and statistical analyses to support the rational development of design requirements for multispectral imaging systems. Among other things, the ART affords a capability to synthesize coarser-spatial-resolution image-data products. The ART also provides for simulations of image-degradation effects, including point-spread functions, misregistration of spectral images, and noise. The ART can utilize real or synthetic data sets, along with sensor specifications, to create simulated data sets. In one example of a particular application, simulated imagery of a coarse resolution system was created using high-resolution imagery from another system in order to perform a radiometric cross-comparison. In the case of a proposed sensor system, the simulated data can be used to conduct trade studies and statistical analyses to ensure that the sensor system will satisfy the requirements of potential scientific, academic, and commercial user communities.

Zanoni, Vicki↗
NASA NTRS2016

DNS of Low-Pressure Turbine Cascade Flows with Elevated Inflow Turbulence Using a Discontinuous-Galerkin Spectral-Element Method

Recent progress towards developing a new computational capability for accurate and efficient high-fidelity direct numerical simulation (DNS) and large-eddy simulation (LES) of turbomachinery is described. This 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 a computationally efficient manner on a modern high performance computer architecture. An inflow turbulence generation procedure based on a linear forcing approach has been incorporated in this framework and DNS conducted to study the effect of inflow turbulence on the suction- side separation bubble in low-pressure turbine (LPT) cascades. The T106 series of airfoil cascades in both lightly (T106A) and highly loaded (T106C) configurations at exit isentropic Reynolds numbers of 60,000 and 80,000, respectively, are considered. The numerical simulations are performed using 8th-order accurate spatial and 4th-order accurate temporal discretization. The changes in separation bubble topology due to elevated inflow turbulence is captured by the present method and the physical mechanisms leading to the changes are explained. The present results are in good agreement with prior numerical simulations but some expected discrepancies with the experimental data for the T106C case are noted and discussed.

Garai, Anirban↗
DOE OSTI2024

Data-Driven Invertible Neural Surrogates of Atmospheric Transmission

We present Data-Driven Invertible Neural Surrogates of Atmospheric transmission, or DINSAT. DINSAT is a novel framework for inferring an atmospheric transmission profile from a spectral scene. This framework leverages a lightweight, physics-based simulator that is automatically tuned -- by virtue of autodifferentiation and differentiable programming -- to construct a surrogate atmospheric profile to model the observed data. The framework has utility in (i) performing atmospheric correction, (ii) recasting spectral data between various modalities (e.g. radiance and reflectance at the surface and at the sensor), and (iii) inferring atmospheric transmission profiles, such as absorbing bands and their relative magnitudes. We demonstrate the utility of these methods by performing a canonical atmospheric correction task for the purposes of further analysis - in this case, target detection within a scene.

Koch, James V.↗
DOE OSTI2011

K+A GALAXIES AS THE AFTERMATH OF GAS-RICH MERGERS: SIMULATING THE EVOLUTION OF GALAXIES AS SEEN BY SPECTROSCOPIC SURVEYS

Models of poststarburst (or 'K+A') galaxies are constructed by combining fully three-dimensional hydrodynamic simulations of galaxy mergers with radiative transfer calculations of dust attenuation. Spectral line catalogs are generated automatically from moderate-resolution optical spectra calculated as a function of merger progress in each of a large suite of simulations. The mass, gas fraction, orbital parameters, and mass ratio of the merging galaxies are varied systematically, showing that the lifetime and properties of the K+A phase are strong functions of the merger scenario. K+A durations are generally {approx}<0.1-0.3 Gyr, significantly shorter than the commonly assumed 1 Gyr, which is obtained only in rare cases, owing to a wide variation in star formation histories resulting from different orbital and progenitor configurations. Combined with empirical merger rates, the model lifetimes predict rapidly rising K+A fractions as a function of redshift that are consistent with results of large spectroscopic surveys, resolving tension between the observed K+A abundance and that predicted when one assumes the K+A duration is the lifetime of A stars ({approx}1 Gyr). These simulated spectra are spatially resolved on scales of about 1 kpc, and indicate that a centrally concentrated starburst causes the Balmer absorption strengths to increase toward the central few kiloparsecs of the remnant. The effects of dust attenuation, viewing angle, and aperture bias on our models are analyzed. In some cases, the K+A features are longer-lived and more pronounced when active galactic nucleus (AGN) feedback removes dust from the center, uncovering the young stars formed during the burst. In this picture, the K+A phase begins during or shortly after the bright starburst/AGN phase in violent mergers, and thus offers a unique opportunity to study the effects of quasar and star formation feedback on the gas reservoir and evolution of the remnant. Analytic fitting formulae are provided for the estimates of K+A incidence as a function of merger scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2013

Damage Detection in Composite Structures with Wavenumber Array Data Processing

Guided ultrasonic waves (GUW) have the potential to be an efficient and cost-effective method for rapid damage detection and quantification of large structures. Attractive features include sensitivity to a variety of damage types and the capability of traveling relatively long distances. They have proven to be an efficient approach for crack detection and localization in isotropic materials. However, techniques must be pushed beyond isotropic materials in order to be valid for composite aircraft components. This paper presents our study on GUW propagation and interaction with delamination damage in composite structures using wavenumber array data processing, together with advanced wave propagation simulations. Parallel elastodynamic finite integration technique (EFIT) is used for the example simulations. Multi-dimensional Fourier transform is used to convert time-space wavefield data into frequency-wavenumber domain. Wave propagation in the wavenumber-frequency domain shows clear distinction among the guided wave modes that are present. This allows for extracting a guided wave mode through filtering and reconstruction techniques. Presence of delamination causes spectral change accordingly. Results from 3D CFRP guided wave simulations with delamination damage in flat-plate specimens are used for wave interaction with structural defect study.

Tian, Zhenhua↗
NASA NTRS1991

Computational study of three dimensional viscous flow through a turbine cascade using a multi-domain spectral technique

The three dimensional viscous flow through a planar turbine cascade is numerically simulated by direct solution of the incompressible Navier-Stokes equations. Flow dependence in the spanwise direction is represented by direct expansion in Chebyshev polynomials, while the discretization on planes parallel to the endwalls is accomplished using the spectral element method. Elemental mapping from the physical to the computational space uses an algebraic mapping technique. A fractional time stepping method that consists of an explicit nonlinear convective step, an implicit pressure correction step, and an implicit viscous step is used to advance the Navier-Stokes equations forward in time. Results computed at moderate Reynolds numbers show a three dimensional endwall flow separation, a midspan separation of the blade suction surface boundary layer, and other three-dimensional features such as the presence of a saddle point flow in the endwall region. In addition, the computed skin friction lines are shown to be orthogonal to the surface vorticity lines, demonstrating the accuracy achievable in the present method.

Renaud, Earl W.↗
NASA NTRSDate not supplied

Nanoflare Heating of an X-Ray Bright Point

Nanoflares are thought to be one of the prime candidates that can keep the solar corona to its multimillion kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing the nanoflare heated simulated plasma emissions with the observed emission. Here, we present a simulation of emission from an X-ray Bright Point (XBP) that was observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with concurrent observations from SDO/AIA and Hinode/XRT. We use EBTEL hydrodynamic code to simulate the XBP loops. Length and magnetic field strength of these loops are derived from the potential field extrapolation of the observed photospheric magnetogram by HMI/SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the looplength and magnetic field strength. The simulated outputs are used to predict the intensity of spectrally pure map of Fe-18, Fe-17, Ne-9 ,O-8, O-9, Ne-9 etc, which are then compared with the derived intensity from MaGIXS observation. Further we have predicted the intensity map as observed by AIA and XRT and compared them with the observation. We also estimated the temperature distribution of the XBP from the simulation and found a good agreement with the derived distribution from MaGIXS observation.

coronal heating↗
DOE OSTI2013

Direct numerical simulation of electrokinetic instability and transition to chaotic motion

A new type of instability—electrokinetic instability—and an unusual transition to chaotic motion near a charge-selective surface (semiselective electric membrane, electrode, or system of micro-/nanochannels) was studied by the numerical integration of the Nernst-Planck-Poisson-Stokes system and a weakly nonlinear analysis near the threshold of instability. A special finite-difference method was used for the space discretization along with a semi-implicit 31/3 -step Runge-Kutta scheme for the integration in time. Two kinds of initial conditions were considered: (a) white-noise initial conditions to mimic “room disturbances” and subsequent natural evolution of the solution, and (b) an artificial monochromatic ion distribution with a fixed wave number to simulate regular wave patterns. The results were studied from the viewpoint of hydrodynamic stability and bifurcation theory. The threshold of electroconvective movement was found by the linear spectral stability theory, the results of which were confirmed by numerical simulation of the entire system. Our weakly nonlinear analysis and numerical integration of the entire system predict possibility of both kinds of bifurcations at the critical point, supercritical and subcritical, depending on the system parameters. The following regimes, which replace each other as the potential drop between the selective surfaces increases, were obtained: one-dimensional steady solution, two-dimensional steady electroconvective vortices (stationary point in a proper phase space), unsteady vortices aperiodically changing their parameters (homoclinic contour), periodic motion (limit cycle), and chaotic motion. The transition to chaotic motion does not include Hopf bifurcation. The numerical resolution of the thin concentration polarization layer showed spike-like charge profiles along the surface, which could be, depending on the regime, either steady or aperiodically coalescent. The numerical investigation confirmed the experimentally observed absence of regular (near-sinusoidal) oscillations for the overlimiting regimes. There is a qualitative agreement of the experimental and the theoretical values of the threshold of instability, the dominant size of the observed coherent structures, and the experimental and theoretical volt–current characteristics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Characterization of an x-ray phase contrast imaging system based on the miniature synchrotron MIRRORCLE-6X

Purpose: The implementation of in-line x-ray phase contrast imaging (PCI) for soft-tissue patient imaging is hampered by the lack of a bright and spatially coherent x-ray source that fits into the hospital environment. This article provides a quantitative characterization of the phase-contrast enhancement of a PCI system based on the miniature synchrotron technology MIRRORCLE-6X. Methods: The phase-contrast effect was measured using an edge response of a plexiglass plate as a function of the incident angle of radiation. We have developed a comprehensive x-ray propagation model based on the system's components, properties, and geometry in order to interpret the measurement data. Monte-Carlo simulations are used to estimate the system's spectral properties and resolution. Results: The measured ratio of the detected phase-contrast to the absorption contrast is currently in the range 100% to 200%. Experiments show that with the current implementation of the MIRRORCLE-6X, a target smaller than 30-40 {mu}m does not lead to a larger phase-contrast. The reason for this is that the fraction of x-rays produced by the material (carbon filament and glue) that is used for mounting the target in the electron beam is more than 25% of the total amount of x-rays produced. This increases the apparent source size. The measured phase-contrast is at maximum two times larger than the absorption contrast with the current set-up. Conclusions: Calculations based on our model of the present imaging system predict that the phase-contrast can be up to an order of magnitude larger than the absorption contrast in case the materials used for mounting the target in the electron beam do not (or hardly) produce x-rays. The methods described in this paper provide vital feedback for guiding future modifications to the design of the x-ray target of MIRRORCLE-type system and configuration of the in-line PCI systems in general.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

L 98-59: A Benchmark System of Small Planets for Future Atmospheric Characterization

The M3V dwarf star L 98-59 hosts three small (R < 1.6 R {sub ⊕}) planets. The host star is bright (K = 7.1) and nearby (10.6 pc), making the system a prime target for follow-up characterization with the Hubble Space Telescope (HST) and the upcoming James Webb Space Telescope (JWST). Herein, we use simulated transmission spectroscopy to evaluate the detectability of spectral features with HST and JWST assuming diverse atmospheric scenarios (e.g., atmospheres dominated by H{sub 2}, H{sub 2}O, CO{sub 2}, or O{sub 2}). We find that H{sub 2}O and CH{sub 4} present in a low mean molecular weight atmosphere could be detected with HST in one transit for the two outermost planets, while H{sub 2}O in a clear steam atmosphere could be detected in six transits or fewer with HST for all three planets. We predict that observations using JWST/NIRISS would be capable of detecting a clear steam atmosphere in one transit for each planet and H{sub 2}O absorption in a hazy steam atmosphere in two transits or less. In a clear, desiccated atmosphere, O{sub 2} absorption may be detectable for all three planets with NIRISS. If the L 98-59 planets possess a clear, Venus-like atmosphere, NIRSpec could detect CO{sub 2} within 26 transits for each planet, but the presence of H{sub 2}SO{sub 4} clouds would significantly suppress CO{sub 2} absorption. The L 98-59 system is an excellent laboratory for comparative planetary studies of transiting multiplanet systems, and observations of the system via HST and JWST would present a unique opportunity to test the accuracy of the models presented in this study.

47 OTHER INSTRUMENTATION↗