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,171 records · Page 65

Nanoflare Heating Frequency of an X-ray Bright Point Observed by MaGIXS

Nanoflares have been considered to be one of the most likely candidates for heating the solar corona to multi-million kelvin temperatures. Individual nanoflares are difficult to detect with today's instruments, but their presence may be established by comparing simulated nanoflare-heated plasma emissions to observed emissions. We present a simulation of emission from an X-ray Bright Point (XBP) detected by the MaGIXS, as well as simultaneous observations from SDO/AIA and Hinode/XRT. To simulate the XBP loops, we utilize the HYDRAD code. The length and magnetic field strength of these loops are determined using potential field extrapolation of SDO/HMI's observed photospheric magnetogram. Each loop is considered to be heated by random nanoflares, the amplitude and frequency of which are governed by the length of the loop and the strength of the magnetic field. The simulated outputs are used to estimate the intensity of spectrally pure maps of Fe-18, Fe-17, Ne-9, O-8, O-9, Ne-9, and so on, which is then compared to the intensity determined from MaGIXS observations. In addition, we derived the intensity maps obtained by AIA and XRT and compared them to the observed data. The composite distribution of the delay time of the nanoflares for which the simulated loops morphology and intensities match with observation shows a peak at 200s-500s, indicating that most of the nanoflares have a high/intermediate frequency.

coronal heating↗

THE INFRARED SPECTRUM OF PROTONATED OVALENE IN SOLID PARA-HYDROGEN AND ITS POSSIBLE CONTRIBUTION TO INTERSTELLAR UNIDENTIFIED INFRARED EMISSION

The mid-infrared emission from galactic objects, including reflection nebulae, planetary nebulae, proto-planetary nebulae, molecular clouds, etc, as well as external galaxies, is dominated by the unidentified infrared (UIR) emission bands. Large protonated polycyclic aromatic hydrocarbons (H{sup +}PAHs) were proposed as possible carriers, but no spectrum of an H{sup +}PAH has been shown to exactly match the UIR bands. Here, we report the IR spectrum of protonated ovalene (7-C{sub 32}H{sub 15} {sup +}) measured in a para -hydrogen ( p -H{sub 2}) matrix at 3.2 K, generated by bombarding a mixture of ovalene and p -H{sub 2} with electrons during matrix deposition. Spectral assignments were made based on the expected chemistry and on the spectra simulated with the wavenumbers and infrared intensities predicted with the B3PW91/6-311++G(2d,2p) method. The close resemblance of the observed spectral pattern to that of the UIR bands suggests that protonated ovalene may contribute to the UIR emission, particularly from objects that emit Class A spectra, such as the IRIS reflection nebula, NGC 7023.

79 ASTRONOMY AND ASTROPHYSICS↗

TPSAS-NF1676L-18357-DND

High-spectral resolution measurements (in the IR region) from space born infrared sensors can provide sufficient spectral information about cloud properties like cloud top pressure, cloud particle size and cloud optical thickness. Empirical orthogonal function (EOF) based regression scheme was used to retrieve cloud properties along with atmospheric state profiles. A channel based physical inversion scheme was also tried on cloud-contaminated radiances to improve retrieval accuracy achieved by regression scheme. Recently developed principle component (PC) based physical retrieval algorithm utilized the EOF transformation to compress both hyperspectral radiance spectra and atmospherical profiles into PC scores to greatly boost the computational speed while maintaining the information content. The accuracy of the (PC) based radiation transfer model (PCRTM) to simulate top-of the atmosphere(TOA) high-spectral resolution radiance spectra for both clear and cloudy sky conditions has been well tested by comparing with the other line-by-line radiation transfer model and using the observation data. The application of the PCRTM based physical inversion method on simulated IASI radiance under 'cloudy' conditions will be a good integrity test for the new algorithm.

X Liu↗

Aerosol Indirect Effects on Water Vapor in the UTLS of Typhoon Saomai (2006)

The evolution of landfalling Typhoon Saomai (2006) is simulated using the Weather Research and Forecasting model with spectral bin microphysics scheme, showing tropical cyclone (TC) could vertically transport water vapor from troposphere to upper troposphere and lower stratosphere (UTLS). Three aerosol sensitivity experiments are conducted to explore the indirect effects of aerosol served as cloud condensation nuclei (CCN) on overshooting convection (OTC) and water vapor content in UTLS, namely MAR (maritime), MIX (semi‐continental) and CON (continental). Results revealed that CON simulation with high initial aerosol concentration showing invigorated convection at the periphery of the TC, which interfered with the reformation of TC eyewall, leading to TC weakening finally. In addition, CON simulated the least water content in UTLS through suppressing the number of convection overshooting the tropopause level. Most overshooting events take place during TC developing stage in the eyewall region before landfall. Higher levels of CCN concentration could lead to decreased number of OTC in TC, consistently reduced the water vapor content in UTLS. CON and MIX experiments simulated respectively 21.9% (16.2%) and 9.96% (1.5%) less OTC (water vapor content in UTLS) than MAR simulation. Further, water vapor content in UTLS is affected by the interaction of OTC and ice‐phase microphysics. The interesting point is that the CCN sensitivity experiments show that water vapor content in UTLS of a TC is sensitive to the CCN concentration in the boundary layer.

54 ENVIRONMENTAL SCIENCES↗

Convective and Turbulent Motions in Nonprecipitating Cu. Part III: Characteristics of Turbulence Motions

Velocity field in a nonprecipitating Cu under BOMEX conditions, simulated by SAM with 10-m resolution and spectral bin microphysics is separated into the convective part and the turbulent part, using a wavelet filtering. In Part II of the study properties of convective motions of this Cu were investigated. Here in Part III of the study, the parameters of cloud turbulence are calculated in the cloud updraft zone at different stages of cloud development. The main points of this study are (i) application of a fine-scale LES model of a single convective cloud allowed a direct estimation of turbulence parameters using the resolved flow in the cloud and (ii) the separation of the resolved flow into the turbulence flow and the nonturbulence flow allowed us to estimate different turbulent parameters with sufficient statistical accuracy. We calculated height and time dependences of the main turbulent parameters such as turbulence kinetic energy (TKE), spectra of TKE, dissipation rate, and the turbulent coefficient. It was found that the main source of turbulence in the cloud is buoyancy whose contribution is described by the buoyancy production term (BPT). The shear production term (SPT) increases with height and reaches its maximum near cloud top, and so does BPT. In agreement with the behavior of BPT and SPT, turbulence in the lower cloud part (below the inversion level) is weak and hardly affects the processes of mixing and entrainment. The fact that BPT is larger than SPT determines many properties of cloud turbulence. For instance, the turbulence is nonisotropic, so the vertical component of TKE is substantially larger than the horizontal components. Another consequence of the fact that BPT is larger than STP manifests itself in the finding that the turbulence spectrum largely obeys the -11/5 Bolgiano–Obukhov scaling. The classical Kolmogorov -5/3 scaling dominates for the low part of a cloud largely at the dissolving stage of cloud evolution. Using the spectra obtained we evaluated an “effective” dissipation rate which increases with height from nearly zero at cloud base up to 20 cm 2 s -3 near cloud top. The coefficient of turbulent diffusion was found to increase with height and ranged from 5 m 2 s -1 near cloud base to 25 m 2 s -1 near cloud top. In conclusion, the possible role of turbulence in the process of lateral entrainment and mixing is discussed.

54 ENVIRONMENTAL SCIENCES↗

Advanced Simulation of ITER Core X-ray Crystal Spectroscopy

X-Ray Simulation Analysis (XRSA) is an analytical ray-tracing mixed code developed specifically for the ITER Core X-Ray Crystal Spectroscopy (XRCS-Core) diagnostic, which employs a dual-reflection configuration incorporating multiple pre-reflectors made of Highly Oriented Pyrolytic Graphite (HOPG) and spherically curved analyzing crystals. The ITER XRCS-Core is designed for high spectral resolution measurement in specific wavelength ranges, including narrow bands around 1.354 Å for W 64+ , 2.19 Å for Xe 51+ , and 2.555 Å for Xe 44+ and Xe 47+ , enabling diagnostic capability across a broad electron temperature range in the ITER plasma. XRSA facilitates efficient simulation of the spectral performance of this complex X-ray spectroscopic system. Recent updates to the XRSA code have incorporated two critical effects: auto-focusing, which specifically applies to HOPG, and polarization. These two effects are particularly important in the dual-reflection configuration used in the ITER XRCS-Core system to provide more accurate modeling results. Here, simulations conducted with the updated code demonstrate that polarization has a substantial impact on the performance of the dual-reflection system. Additionally, the combined influence of polarization and system layout introduces performance variations across channels through the same crystal.

Computer graphics↗

Numerical experiments in boundary-layer stability

Numerical solution of the three-dimensional incompressible Navier-Stokes equations is used to study the instability of a flat-plate boundary layer in a manner analogous to the vibrating-ribbon experiments. Flow-field structures are observed which are very similar to those found in the vibrating-ribbon experiment to which computational initial conditions have been matched. Streamwise periodicity is assumed in the simulation so that the evolution occurs in time, but the events which constitute the instability are so similar to the spatially occurring ones of the laboratory that it seems clear the physical processes involved are the same. A spectral and finite difference numerical algorithm is employed in the simulation.

Wray, A.↗

Decoding the Atomic Structure of Ga 2 Te 5 Pulsed Laser Deposition Films for Memory Applications Using Diffraction and First-Principles Simulations

Neuromorphic computing, reconfigurable optical metamaterials that are operational over a wide spectral range, holographic and nonvolatile displays of extremely high resolution, integrated smart photonics, and many other applications need next-generation phase-change materials (PCMs) with better energy efficiency and wider temperature and spectral ranges to increase reliability compared to current flagship PCMs, such as Ge 2 Sb 2 Te 5 or doped Sb 2 Te. Gallium tellurides are favorable compounds to achieve the necessary requirements because of their higher melting and crystallization temperatures, combined with low switching power and fast switching rate. Ga 2 Te 3 and non-stoichiometric alloys appear to be atypical PCMs; they are characterized by regular tetrahedral structures and the absence of metavalent bonding. The sp 3 gallium hybridization in cubic and amorphous Ga 2 Te 3 is also different from conventional p-bonding in flagship PCMs, raising questions about its phase-change mechanism. Furthermore, gallium tellurides exhibit a number of unexpected and highly unusual phenomena, such as nanotectonic compression and viscosity anomalies just above their melting points. Using high-energy X-ray diffraction, supported by first-principles simulations, we will elucidate the atomic structure of amorphous Ga 2 Te 5 PLD films, compare it with the crystal structure of tetragonal gallium pentatelluride, and investigate the electrical, optical, and thermal properties of these two materials to assess their potential for memory applications, among others.

36 MATERIALS SCIENCE↗

NEOCAM: Near Earth Object Chemical Analysis Mission: Bridging the Gulf between Telescopic Observations and the Chemical and Mineralogical Compositions of Asteroids or Diogenes A: Diagnostic Observation of the Geology of Near Earth Spectrally-Classified Asteroids

Studies of meteorites have yielded a wealth of scientific information based on highly detailed chemical and isotopic studies possible only in sophisticated terrestrial laboratories. Telescopic studies have revealed an enormous (greater than 10(exp 5)) number of physical objects ranging in size from a few tens of meters to several hundred kilometers, orbiting not only in the traditional asteroid belt between Mars and Jupiter but also throughout the inner solar system. Many of the largest asteroids are classed into taxonomic groups based on their observed spectral properties and are designated as C, D. X, S or V types (as well as a wide range in sub-types). These objects are certainly the sources far the meteorites in our laboratories, but which asteroids are the sources for which meteorites? Spectral classes are nominally correlated to the chemical composition and physical characteristics of the asteroid itself based on studies of the spectral changes induced in meteorites due to exposure to a simulated space environment. While laboratory studies have produced some notable successes (e.g. the identification of the asteroid Vesta as the source of the H, E and D meteorite classes), it is unlikely that we have samples of each asteroidal spectral type in our meteorite collection. The correlation of spectral type and composition for many objects will therefore remain uncertain until we can return samples of specific asteroid types to Earth for analyses. The best candidates for sample return are asteroids that already come close to the Earth. Asteroids in orbit near 1 A.U. have been classified into three groups (Aten, Apollo & Amor) based on their orbital characteristics. These Near Earth Objects (NEOs) contain representatives of virtually all spectral types and sub-types of the asteroid population identified to date. Because of their close proximity to Earth, NEOs are prime targets for asteroid missions such as the NEAR-Shoemaker NASA Discovery Mission to Eros and the Japanese Hyabusa Mission to Itokawa. Also due to their close proximity to Earth, NEOs constitute the most likely set of celestial objects that will impact us in the relatively near future.

Nuth, Joseph A.↗

SOLAR FLARE CHROMOSPHERIC LINE EMISSION: COMPARISON BETWEEN IBIS HIGH-RESOLUTION OBSERVATIONS AND RADIATIVE HYDRODYNAMIC SIMULATIONS

Solar flares involve impulsive energy release, which results in enhanced radiation over a broad spectral range and a wide range of heights. In particular, line emission from the chromosphere can provide critical diagnostics of plasma heating processes. Thus, a direct comparison between high-resolution spectroscopic observations and advanced numerical modeling results could be extremely valuable, but has not yet been attempted. In this paper, we present such a self-consistent investigation of an M3.0 flare observed by the Dunn Solar Telescope’s Interferometric Bi-dimensional Spectrometer (IBIS) on 2011 September 24 which we have modeled using the radiative hydrodynamic code RADYN. We obtained images and spectra of the flaring region with IBIS in Hα 6563 Å and Ca ii 8542 Å, and with RHESSI in X-rays. The latter observations were used to infer the non-thermal electron population, which was passed to RADYN to simulate the atmospheric response to electron collisional heating. We then synthesized spectral lines and compared their shapes and intensities to those observed by IBIS and found a general agreement. In particular, the synthetic Ca ii 8542 Å profile fits well to the observed profile, while the synthetic Hα profile is fainter in the core than for the observation. This indicates that Hα emission is more responsive to the non-thermal electron flux than the Ca ii 8542 Å emission. We suggest that it is necessary to refine the energy input and other processes to resolve this discrepancy.

79 ASTRONOMY AND ASTROPHYSICS↗

A high-order computational framework for particle-resolved simulations of disperse multiphase flows

This work presents a high-order numerical approach for particle-resolved simulations of disperse multiphase flows, where the Navier-Stokes equations for fluid flow are solved using a high-order spectral element method in the Eulerian framework, and the particle phase is directly simulated with a discrete element method. The coupling between particles and fluids is explicitly handled using an adapted direct-forcing immersed boundary method. Unlike the conventional schemes, a high-order barycentric Lagrange interpolation method and a Gaussian projection kernel are used to ensure accurate momentum exchange between local boundary points and surrounding fluid nodes in the framework of high-order fluid solver. Benchmark tests of increasing complexity are conducted to demonstrate the accuracy and efficiency of our method. Here, it is found that our approach exhibits an excellent convergence performance, as the fluid element/grid is refined and the number of boundary points increases. Compared to conventional low-order methods, the proposed high-order framework enables the use of substantially larger fluid elements while maintaining high accuracy in modeling fluid-particle interactions, owing to the enhanced resolution of high-order basis functions. Moreover, since the primary unknowns are stored at element or grid nodes, the high-order approach offers improved efficiency in both CPU memory usage and total computational cost.

42 ENGINEERING↗

Mixed convection heat transfer in a differentially heated cavity with two rotating cylinders

Mixed convection heat transfer in a differentially heated cavity with two rotating cylinders was studied numerically for various pertinent parameters such as Richardson number, Reynolds number, non-dimensional rotation speed of the cylinder, and the location of the cylinders. The results from these simulations were validated using an open-source spectral element code, Nek5000. The results of this study were presented in terms of streamlines, isotherms, local and average Nusselt number, and temperature profiles. The results illustrated that the average Nusselt number depends strongly on the rotation speed of the cylinder, Reynolds number, and Richardson number. However, the average Nusselt number was found to be independent of the rotation speed at high Richardson number (Ri = 10). In addition, the direction of the rotation speed of the cylinders has a profound effect on the average Nusselt. Finally, the study illustrated that the magnitude and direction of the rotation speed of the cylinders have a significant effect on the flow pattern and isotherms.

42 ENGINEERING↗

OpenACC acceleration for the P N – P N - 2 algorithm in Nek5000

Due to its high performance and throughput capabilities, GPU-accelerated computing is becoming a popular technology in scientific computing, in particular using programming models such as CUDA and OpenACC. The main advantage with OpenACC is that it enables to simply port codes in their ‘‘original’’ form to GPU systems through compiler directives, thus allowing an incremental approach. An OpenACC implementation is applied to the CFD code Nek5000 for simulation of incompressible flows, based on the spectral-element method. The work follows up previous implementations and focuses now on the P N - P N-2 method for the spatial discretization of the Navier–Stokes equations. Performance results of the ported code show a speed-up of up to 3.1 on multi-GPU for a polynomial order N > 11.

97 MATHEMATICS AND COMPUTING↗

Design of a correlation electron cyclotron emission diagnostic for Alcator C-Mod

A correlation electron cyclotron emission (CECE) diagnostic has been installed in Alcator C-Mod. In order to measure electron temperature fluctuations, this diagnostic uses a spectral decorrelation technique. Constraints obtained with nonlinear gyrokinetic simulations guided the design of the optical system and receiver. The CECE diagnostic is designed to measure temperature fluctuations which have k{sub {theta}}{<=} 4.8 cm{sup -1} (k{sub {theta}}{rho}{sub s} < 0.5) using a well-focused beam pattern. Because the CECE diagnostic is a dedicated turbulence diagnostic, the optical system is also flexible, which allows for various collimating lenses and antenna to be used. The system overview and the demonstration of its operability as designed are presented in this paper.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Impact of Xe partial pressure on the production of excimer vacuum ultraviolet emission for plasma display panels

In this work, the effect of the Xe partial pressure on the excimer vacuum ultraviolet (VUV) emission intensity of the plasma display panels is investigated, both by measuring the spectral emission directly and by two-dimensional simulations. Experimentally, we find that at the high Xe partial pressure levels, there is an supra-linear increase of excimer VUV radiation and that determines the strong increase of luminance at the high pressures and high voltage. Due to the increase of the luminance and the almost unchanged discharge current, the luminous efficacy strongly increases with the Xe partial pressure. In addition, we also investigated the dynamics of the VUV generation, by measuring the decay time of the excimer VUV light as a function of the gas pressure. It is found that the decay time decreases with the increase of gas pressure. The spatial characteristics of the excimer VUV emission are also discussed. Different from the Ne and near-infrared emission, the excimer VUV emission is generated near the surface of the electrodes and increases uniformly on both sides of the anode and cathode (i.e., the bulk plasma region). Most importantly, it is found that the VUV production occurs during the afterglow period, while it is almost zero at the moment of the discharge itself. From the simulations, it can be seen that the Xe{sub 2}*({sup 3}{Sigma}{sub u}{sup +}) excimer species, which are generated from Xe*(1s{sub 5}), play a dominant role in the excimer VUV emission output at the high Xe partial pressure. The two-dimensional simulations also show that the strong increase of Xe excimer excitation states in the case of high pressure is mainly the result of the high conversion efficiency of the Xe excimer states, especially in the afterglow period. Due to the high conversion efficiency of Xe excitation species to Xe excimer species by the high collision rate in the case of high pressure, there is a strong increase of excimer VUV production, especially from the cathode.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of the LH wave for demo in pulsed and steady state scenario

The Lower Hybrid Current Drive (LHCD) has been analysed in DEMO tokamak plasma in the 'pulsed and steady state regime' considering two plasma scenarios characterized, respectively, by flat density profile and peaked density profiles. We have obtained LH deposition profiles in cases of neglecting the effect of spectral broadening produced by PI at the edge. By comparing the Power Deposition Profiles for both DEMO scenarios ('flat' and 'peaked'), the SOL of DEMO does not play any role in the absorption of the LH wave. In all cases the deposition is localized inside the separatrix layer r/a≤1. By lowering the parallel wave-number peak of the power spectrum from 1.8 to 1.5, the accessibility condition in both case prevents the power from reaching the deposition layer apart from a small fraction which pertains to the higher n∥ of the power spectrum. The spectrum centred at 1.8 is suggested to be useful in DEMO. More realistically, as supported by available data of LHCD in a wide range of operating densities, the effect of parametric decay instability (PDI) can produce a spectral broadening which should be included in the simulations. Further studies would be necessary for assessing the temperature profiles in the SOL at reactor-graded conditions. This is because, if the SOL temperature is at least of the order of 50 to 100 eV, the effect of PDI broads the spectrum up to n∥≤10, and the deposition profile is slightly wider but not much shifted outwards.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact of temperature-velocity distribution on fusion neutron peak shape

Doppler broadening of the 14 MeV DT and 2.45 MeV DD fusion neutron lines has long been our best measure of temperature in a burning plasma. At the National Ignition Facility (NIF), yields are high enough and our neutron spectrometers accurate enough that we see finer details of the peak shape. For example, we can measure the shift of the peak due to the bulk motion of the plasma, and we see indications of non-thermal broadening, skew, and kurtosis of the peak caused by the variations of temperature and fluid velocity during burn. We can also distinguish spectral differences among several lines of sight. Finally, this paper will review the theory of fusion neutron line shape, show examples of non-Gaussian line shapes and directional variations in NIF data, and describe detailed spectral shapes we see in radiation-hydrodynamics simulations of implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The two-dimensional kinetic ballooning theory for ion temperature gradient mode in tokamak

The two-dimensional (2D) kinetic ballooning theory is developed for the ion temperature gradient mode in an up-down symmetric equilibrium (illustrated via concentric circular magnetic surfaces). The ballooning transform converts the basic 2D linear gyro-kinetic equation into two equations: (1) the lowest order equation (ballooning equation) is an integral equation essentially the same as that reported by Dong et al., [Phys. Fluids B 4, 1867 (1992)] but has an undetermined Floquet phase variable, (2) the higher order equation for the rapid phase envelope is an ordinary differential equation in the same form as the 2D ballooning theory in a fluid model [Xie et al., Phys. Plasmas 23, 042514 (2016)]. The system is numerically solved by an iterative approach to obtain the (phase independent) eigen-value. The new results are compared here to the two earlier theories. We find a strongly modified up-down asymmetric mode structure, and non-trivial modifications to the eigen-value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗