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,603 records · Page 89
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

Variability simulations with a steady, linearized primitive equations model

Solutions of the steady, primitive equations on a sphere, linearized about a zonally symmetric basic state are computed for the purpose of simulating monthly mean variability in the troposphere. The basic states are observed, winter monthly mean, zonal means of zontal and meridional velocities, temperatures and surface pressures computed from the 15 year NMC time series. A least squares fit to a series of Legendre polynomials is used to compute the basic states between 20 H and the equator, and the hemispheres are assumed symmetric. The model is spectral in the zonal direction, and centered differences are employed in the meridional and vertical directions. Since the model is steady and linear, the solution is obtained by inversion of a block, pente-diagonal matrix. The model simulates the climatology of the GFDL nine level, spectral general circulation model quite closely, particularly in middle latitudes above the boundary layer. This experiment is an extension of that simulation to examine variability of the steady, linear solution.

Kinter, J. L., III↗
DOE OSTI2025

A Data Library of Liquid Clouds Modelled With a Large Eddy Simulation Framework

We describe a library of atmospheric large eddy simulations (LES) of liquid-phase boundary layer clouds constructed to enable aerosol–cloud–turbulence interaction studies, support parameterization evaluation and development, and provide training data for machine learning applications. The simulations use a modern LES framework designed for high numerical accuracy, coupled to a detailed spectral bin microphysical scheme. Case studies are configured to represent observed conditions in four key global cloud regions—the Northeastern Atlantic, Northeastern Pacific, Continental United States and Southern Ocean—following a semi-idealised approach. The library also includes aerosol concentration halving and doubling experiments to expose the sensitivities of the case studies to aerosol perturbations. Simulation results are compared to observations on a case-by-case basis, then the library's coverage is evaluated in terms of spreads in meteorological factors and atmospheric boundary layer attributes.

aerosol↗
DOE OSTI2017

Solar Flare Termination Shock and Synthetic Emission Line Profiles of the Fe xxi 1354.08 Å Line

Solar flares are among the most energetic phenomena that occur in the solar system. In the standard solar flare model, a fast mode shock, often referred to as the flare termination shock (TS), can exist above the loop-top source of hard X-ray emissions. The existence of the TS has been recently related to spectral hardening of a flare’s hard X-ray spectra at energies >300 keV. Observations of the Fe xxi 1354.08 Å line during solar flares by the Interface Region Imaging Spectrograph ( IRIS ) spacecraft have found significant redshifts with >100 km s{sup −1}, which is consistent with a reconnection downflow. The ability to detect such a redshift with IRIS suggests that one may be able to use IRIS observations to identify flare TSs. Using a magnetohydrodynamic simulation to model magnetic reconnection of a solar flare and assuming the existence of a TS in the downflow of the reconnection plasma, we model the synthetic emission of the Fe xxi 1354.08 line in this work. We show that the existence of the TS in the solar flare may manifest itself in the Fe xxi 1354.08 Å line.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Evaluation of the Impact of Horizontal Grid Spacing in Terra Incognita on Coupled Mesoscale–Microscale Simulations Using the WRF Framework

Coupled mesoscale-microscale simulations are required to provide time-varying weather-dependent inflow and forcing for large-eddy simulations under general flow conditions. Such coupling necessarily spans a wide range of spatial scales (i.e., ~10 m to ~10 km). In this paper, we use simulations that involve multiple nested domains with horizontal grid spacings in the terra incognita (i.e., ≲km) that may affect simulated conditions in both the outer and inner domains. We examine the impact on simulated wind speed and turbulence associated with forcing provided by a terrain with grid spacing in the terra incognita. We perform a suite of simulations that use combinations of varying horizontal grid spacings and turbulence parameterization/modeling using the Weather Research and Forecasting (WRF) Model using a combination of planetary boundary layer (PBL) and large-eddy simulation subgrid-scale (LES-SGS) models. The results are analyzed in terms of spectral energy, turbulence kinetic energy, and proper orthogonal decomposition (POD) energy. The results show that the output from the microscale domain depends on the type of turbulence model (e.g., PBL or LES-SGS model) used for a given horizontal grid spacing but is independent of the horizontal grid spacing and turbulence modeling of the parent domain. Simulation using a single domain produced less POD energy in the first few modes compared to a coupled simulation (one-way nesting) for similar horizontal grid spacing, which highlights that coupled simulations are required to accurately pass the mesoscale features into the microscale domain.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1989

Direct simulations of turbulent flow using finite-difference schemes

A high-order accurate finite-difference approach is presented for calculating incompressible turbulent flow. The methods used include a kinetic energy conserving central difference scheme and an upwind difference scheme. The methods are evaluated in test cases for the evolution of small-amplitude disturbances and fully developed turbulent channel flow. It is suggested that the finite-difference approach can be applied to complex geometries more easilty than highly accurate spectral methods. It is concluded that the upwind scheme is a good candidate for direct simulations of turbulent flows over complex geometries.

Rai, Man Mohan↗
NASA NTRS1979

Operator performance and localized muscle fatigue in a simulated space vehicle control task

Fourier transforms in a special purpose computer were utilized to obtain power spectral density functions from electromyograms of the biceps brachii, triceps brachii, brachioradialis, flexor carpi ulnaris, brachialis, and pronator teres in eight subjects performing isometric tracking tasks in two directions utilizing a prototype spacecraft rotational hand controller. Analysis of these spectra in general purpose computers aided in defining muscles involved in performing the task, and yielded a derived measure potentially useful in predicting task termination. The triceps was the only muscle to show significant differences in all possible tests for simple effects in both tasks and, overall, was the most consistently involved of the six muscles. The total power monitored for triceps, biceps, and brachialis dropped to minimal levels across all subjects earlier than for other muscles. However, smaller variances existed for the biceps, brachioradialis, brachialis, and flexor carpi ulnaris muscles and could provide longer predictive times due to smaller standard deviations for a greater population range.

Lewis, J. L., Jr.↗
NASA NTRS1984

CO2 induced climatic change and spectral variations in the outgoing terrestrial infrared radiation

The published temperature changes produced in general circulation model simulations of CO2 induced climate modification are used to compute the top of the atmosphere, clear sky outgoing infrared radiance changes expected for doubled CO2. A significant wavenumber shift is produced, with less radiance emerging in the 500-800 per cm (20.0-12.5 micron) CO2 band and with more emerging in the 800-1200 per cm (12.5-8.3 micron) window. The effect varies greatly with latitude. The radiance shift in the 2300 per cm (4.3 micron) region is of the order of 10-30 percent for doubled CO2. It is suggested that the 2300 per cm region be carefully monitored as an aid in detecting the climatic effects of increasing CO2. The change in the wavenumber-integrated radiant exitance is at most a few percent.

Charlock, T. P.↗
DOE OSTI2018

Bethe-Salpeter Equation calculations of nitrogen-vacancy defects in diamond

Doped diamonds are of significant fundamental, technological, and commercial interest. For example, boron doping leads to a blue coloration at low concentration and superconductivity at high concentration. Similarly nitrogen doping can cause a variety of colors in diamond. Nitrogen-vacancy defects are also promising as quantum-computing qubits. Theoretical Bethe-Salpeter Equation (BSE) simulations of boron-doped diamonds have been found to give optical absorption spectra in good agreement with their observed blue coloration. In contrast, nitrogen defects are more complex, due in part to a plethora of possible N-defect conformations. Here, in this study, we present a theoretical study of several neutral, isolated, N-defect conformations that combines ab initio DFT calculations for relaxed 64-atom unit cells together with first principles GW/BSE calculations of their optical spectra. Spectral results are compared with available experimental data, and the effects of absorption are simulated for realistic rendered diamonds.

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

Numerical investigation of short-pulse laser radiation propagation in a turbulent atmosphere

An algorithm is presented for the numerical simulation of short-pulse optical radiation propagation in a turbulent atmosphere on the basis of the solution to the parabolic wave equation for the complex spectral amplitude of the wave field by the split-step method. We present examples of the use of this algorithm for simulating the propagation of a pulsed coherent spatially limited beam and a plane wave. It is shown that in the regime of strong optical turbulence the relative variance of fluctuations of energy density of femtosecond radiation becomes much smaller than the relative variance of the intensity of cw radiation. (laser radiation scattering)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1992

Research in computational fluid dynamics and analysis of algorithms

Recently, higher-order compact schemes have seen increasing use in the DNS (Direct Numerical Simulations) of the Navier-Stokes equations. Although they do not have the spatial resolution of spectral methods, they offer significant increases in accuracy over conventional second order methods. They can be used on any smooth grid, and do not have an overly restrictive CFL dependence as compared with the O(N(exp -2)) CFL dependence observed in Chebyshev spectral methods on finite domains. In addition, they are generally more robust and less costly than spectral methods. The issue of the relative cost of higher-order schemes (accuracy weighted against physical and numerical cost) is a far more complex issue, depending ultimately on what features of the solution are sought and how accurately they must be resolved. In any event, the further development of the underlying stability theory of these schemes is important. The approach of devising suitable boundary clusters and then testing them with various stability techniques (such as finding the norm) is entirely the wrong approach when dealing with high-order methods. Very seldom are high-order boundary closures stable, making them difficult to isolate. An alternative approach is to begin with a norm which satisfies all the stability criteria for the hyperbolic system, and look for the boundary closure forms which will match the norm exactly. This method was used recently by Strand to isolate stable boundary closure schemes for the explicit central fourth- and sixth-order schemes. The norm used was an energy norm mimicking the norm for the differential equations. Further research should be devoted to BC for high order schemes in order to make sure that the results obtained are reliable. The compact fourth order and sixth order finite difference scheme had been incorporated into a code to simulate flow past circular cylinders. This code will serve as a verification of the full spectral codes. A detailed stability analysis by Carpenter (from the fluid Mechanics Division) and Gottlieb gave analytic conditions for stability as well as asymptotic stability. This had been incorporated in the code in form of stable boundary conditions. Effects of the cylinder rotations had been studied. The results differ from the known theoretical results. We are in the middle of analyzing the results. A detailed analysis of the effects of the heating of the cylinder on the shedding frequency had been studied using the above schemes. It has been found that the shedding frequency decreases when the wire was heated. Experimental work is being carried out to affirm this result.

Gottlieb, David↗
DOE OSTI2025

High-Fidelity CFD Simulation of Mixed Convection and Forced Convection in a Pebble Bed Test Reactor Core

The Hermes low-power [35-MW(thermal)] reactor will be built and operated by Kairos Power LLC (KP) to demonstrate its fluoride salt-cooled high-temperature reactor (FHR) technology. In the KP FHR, the reactor core is composed of randomly packed pebbles with TRISO fuel particles inside with FLiBe flow upward through the core acting as a coolant. Previous numerical and experimental studies have been limited to either a small-size bed or to a lack of detailed measurements for heat transfer. Here, to address the lack of high-fidelity heat transfer data in a real-size FHR core, in this study, we simulated a pebble bed core with 34 374 pebbles randomly packed, similar to the Hermes reactor's size. The core radius was 14 times that of the pebble diameter, while the core height was 45 times. In this work, we were particularly interested in a mixed convection regime, where buoyancy is important. Therefore, we performed several large-eddy simulations at different Reynolds numbers (160 to 1000) with gravitational force included. The spectral element computational fluid dynamics code NekRS with graphics processing unit acceleration was used for this study. The low-Mach number approximation was applied to address property changes in the FLiBe and to account for buoyancy. A pure hexahedral mesh with 60 million elements was generated by the Voronoi cell method. At the polynomial order of 5, the total degrees of freedom was 7.5 billion. The developed case in this work is the first of its kind in terms of size and complexity. The local numerical data across the domain were obtained and compared with empirical correlations. After examining the data, we found the following conclusions. For pressure drop, the Reger correlation predicted less than a 5% error. On the other hand, for heat transfer, the Wakao correlation outperformed the others. Based on our findings, we recommend the use of the Wakao correlation for the Nusselt number calculation, and for pressure drop, the KTA (Kerntechnischer Ausschuss) correclation, among the available experimental correlations. In conclusion, the Reger direct numerical simulation-driven correlation for pressure drops should also be considered, given its best agreement with our calculations.

Mixed Convection↗
NASA NTRS2011

Photocatalytic Active Radiation Measurements and Use

Photocatalytic materials are being used to purify air, to kill microbes, and to keep surfaces clean. A wide variety of materials are being developed, many of which have different abilities to absorb various wavelengths of light. Material variability, combined with both spectral illumination intensity and spectral distribution variability, will produce a wide range of performance results. The proposed technology estimates photocatalytic active radiation (PcAR), a unit of radiation that normalizes the amount of light based on its spectral distribution and on the ability of the material to absorb that radiation. Photocatalytic reactions depend upon the number of electron-hole pairs generated at the photocatalytic surface. The number of electron-hole pairs produced depends on the number of photons per unit area per second striking the surface that can be absorbed and whose energy exceeds the bandgap of the photocatalytic material. A convenient parameter to describe the number of useful photons is the number of moles of photons striking the surface per unit area per second. The unit of micro-einsteins (or micromoles) of photons per m2 per sec is commonly used for photochemical and photoelectric-like phenomena. This type of parameter is used in photochemistry, such as in the conversion of light energy for photosynthesis. Photosynthetic response correlates with the number of photons rather than by energy because, in this photochemical process, each molecule is activated by the absorption of one photon. In photosynthesis, the number of photons absorbed in the 400 700 nm spectral range is estimated and is referred to as photosynthetic active radiation (PAR). PAR is defined in terms of the photosynthetic photon flux density measured in micro-einsteins of photons per m2 per sec. PcAR is an equivalent, similarly modeled parameter that has been defined for the photocatalytic processes. Two methods to measure the PcAR level are being proposed. In the first method, a calibrated spectrometer with a cosine receptor is used to measure the spectral irradiance. This measurement, in conjunction with the photocatalytic response as a function of wavelength, is used to estimate the PcAR. The photocatalytic response function is determined by measuring photocatalytic reactivity as a function of wavelength. In the second method, simple shaped photocatalytic response functions can be simulated with a broad-band detector with a cosine receptor appropriately filtered to represent the spectral response of the photocatalytic material. This second method can be less expensive than using a calibrated spectrometer.

Davis, Bruce A.↗
DOE OSTI2011

CAN OVERTURNING MOTIONS IN PENUMBRAL FILAMENTS BE DETECTED?

Numerical simulations indicate that the filamentation of sunspot penumbrae and the associated systematic outflow (the Evershed effect) are due to convectively driven fluid motions constrained by the inclined magnetic field. We investigate whether these motions, in particular the upflows in the bright filaments and the downflows at their edges, can be reliably observed with existing instrumentation. We use a snapshot from a sunspot simulation to calculate two-dimensional maps of synthetic line profiles for the spectral lines Fe I 7090.4 A and C I 5380.34 A. The maps are spatially and spectrally degraded according to typical instrument properties. Line-of-sight velocities are determined from line bisector shifts. We find that the detectability of the convective flows is strongly affected by spatial smearing, particularly so for the downflows. Furthermore, the line-of-sight velocities are dominated by the Evershed flow unless the observation is made very near the disk center. These problems may have compromised recent attempts to detect overturning penumbral convection. Lines with a low formation height are best suited for detecting the convective flows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

On- and off-axis spectral emission features from laser-produced gas breakdown plasmas

Laser-heated gas breakdown plasmas or sparks emit profoundly in the ultraviolet and visible region of the electromagnetic spectrum with contributions from ionic, atomic, and molecular species. Laser created kernels expand into a cold ambient with high velocities during its early lifetime followed by confinement of the plasma kernel and eventually collapse. However, the plasma kernels produced during laser breakdown of gases are also capable of exciting and ionizing the surrounding ambient medium. Two mechanisms can be responsible for excitation and ionization of surrounding ambient: viz. photoexcitation and ionization by intense ultraviolet emission from the sparks produced during the early times of its creation and/or heating by strong shocks generated by the kernel during its expansion into the ambient. In this study, an investigation is made on the spectral features of on- and off-axis emission features of laser-induced plasma breakdown kernels generated in atmospheric pressure conditions with an aim to elucidate the mechanisms leading to ambient excitation and emission. Pulses from an Nd:YAG laser emitting at 1064 nm with 6 ns pulse duration are used to generate plasma kernels. Laser sparks were generated in air, argon, and helium gases to provide different physical properties of expansion dynamics and plasma chemistry considering the differences in laser absorption properties, mass density and speciation. Point shadowgraphy and time-resolved imaging were used to evaluate the shock wave and spark self-emission morphology at early and late times while space and time resolved spectroscopy is used for evaluating the emission features as well as for inferring plasma fundaments at on- and off-axis. Structure and dynamics of the plasma kernel obtained using imaging techniques are also compared to numerical simulations using computational fluid dynamics code. The emission from the kernel showed that spectral features from ions, atoms and molecules are separated in time with an early time temperatures and densities in excess of 35000 K and 4×10 18 /cm 3 with an existence of thermal equilibrium. However, the emission from the off-kernel positions from the breakdown plasmas showed enhanced ultraviolet radiation with the presence of N 2 bands and represented by non-LTE conditions. Finally, our results also highlight that the ultraviolet radiation emitted during early time of spark evolution is the predominant source of the photo-excitation of the surrounding medium.

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

Load response of a floating wind turbine to turbulent atmospheric flow

The two turbulence-generation models [Kaimal Spectrum Exponential Coherence (KSEC) and Mann] specified in the international standard for wind turbine design assume neutral atmospheric conditions and are based on statistical and spectral methods. Mainly due to the lack of physics, the flow fields simulated with these models ultimately differ in their underlying structure, especially in terms of the spatial coherence of longitudinal velocity perturbations. While this may not be critical for smaller wind turbine rotors, it becomes important when rotor sizes increase. Furthermore, it might be especially important in the context of floating technologies as they are more sensitive to large turbulent coherent structures. Previous work found that these differences between KSEC and Mann can propagate to loads predictions and thereby affect the design space of the entire wind turbine system. It is therefore crucial to determine in which ways these two models are underperforming. Up until now, validation of these models had only been done in the vertical direction because it is extremely difficult to obtain atmospheric turbulence measurements separated laterally, and sampled at heights relevant to wind energy. In this work, we address the lack of measurements by using high-fidelity, high-resolution simulation data as a reference. We perform hour-long, large-eddy simulations of turbulent velocity fields that are stability-dependent and contain three-dimensional coherent structures. These flow fields are then used to investigate which stochastic model is a better predictor of loads on a realistic spar-system floating offshore wind turbine, and to quantify how the assumption of neutral stratification propagates to short-term load estimates. Both stochastic turbulence models are found to overpredict fatigue loading in high-wind scenarios (in some cases, by more than 25%) and underpredict it when the wind speed is low (by as much as 20%). The KSEC model matches the high-fidelity flow fields more closely than Mann at high wind speeds, and the opposite is true at low wind speeds. Finally, turbine loading is found to be sensitive to atmospheric stability even when the turbulence intensity remains fairly constant. This sensitivity is most pronounced at low wind speeds, when fatigue load estimates on the spar system can differ by 40%.

17 WIND ENERGY↗
DOE OSTI2016

Planck intermediate results: XXXIX. The Planck list of high-redshift source candidates

The Planck mission, thanks to its large frequency range and all-sky coverage, has a unique potential for systematically detecting the brightest, and rarest, submillimetre sources on the sky, including distant objects in the high-redshift Universe traced by their dust emission. In this paper, a novel method, based on a component-separation procedure using a combination of Planck and IRAS data, has been validated and characterized on numerous simulations, and applied to select the most luminous cold submillimetre sources with spectral energy distributions peaking between 353 and 857 GHz at 5' resolution. A total of 2151 Planck high-z source candidates (the PHZ) have been detected in the cleanest 26% of the sky, with flux density at 545 GHz above 500 mJy. Embedded in the cosmic infrared background close to the confusion limit, these high-z candidates exhibit colder colours than their surroundings, consistent with redshifts z > 2, assuming a dust temperature of T xgal = 35 K and a spectral index of β xgal = 1.5. Exhibiting extremely high luminosities, larger than 10 14 L ⊙ , the PHZ objects may be made of multiple galaxies or clumps at high redshift, as suggested by a first statistical analysis based on a comparison with number count models. Furthermore, first follow-up observations obtained from optical to submillimetre wavelengths, which can be found in companion papers, have confirmed that this list consists of two distinct populations. A small fraction (around 3%) of the sources have been identified as strongly gravitationally lensed star-forming galaxies at redshift 2 to 4, while the vast majority of the PHZ sources appear as overdensities of dusty star-forming galaxies, having colours consistent with being at z > 2, and may be considered as proto-cluster candidates. The PHZ provides an original sample, which is complementary to the Planck Sunyaev-Zeldovich Catalogue (PSZ2); by extending the population of virialized massive galaxy clusters detected below z < 1.5 through their SZ signal to a population of sources at z > 1.5, the PHZ may contain the progenitors of today’s clusters. Therefore the Planck list of high-redshift source candidates opens a new window on the study of the early stages of structure formation, particularly understanding the intensively star-forming phase at high-z.

79 ASTRONOMY AND ASTROPHYSICS↗