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,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

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.↗
DOE OSTI2015

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↗
DOE OSTI2025

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↗
DOE OSTI2018

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↗
DOE OSTI2019

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↗
DOE OSTI2014

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↗
DOE OSTI2017

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↗
DOE OSTI2017

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↗
DOE OSTI2020

Decoding the 2D IR spectrum of the aqueous proton with high-level VSCF/VCI calculations

The aqueous proton is a common and long-studied species in chemistry, yet there is currently intense interest devoted to understanding its hydration structure and transport dynamics. Typically described in terms of two limiting structures observed in gas-phase clusters, the Zundel H 5 O 2 + and Eigen H 9 O 4 + ions, the aqueous structure is less clear due to the heterogeneity of hydrogen bonding environments and room-temperature structural fluctuations in water. The linear infrared (IR) spectrum, which reports on structural configurations, is challenging to interpret because it appears as a continuum of absorption, and the underlying vibrational modes are strongly anharmonically coupled to each other. Recent two-dimensional IR (2D IR) experiments presented strong evidence for asymmetric Zundel-like motifs in solution, but true structure–spectrum correlations are missing and complicated by the anharmonicity of the system. In this study, we employ high-level vibrational self-consistent field/virtual state configuration interaction calculations to demonstrate that the 2D IR spectrum reports on a broad distribution of geometric configurations of the aqueous proton. We find that the diagonal 2D IR spectrum around 1200 cm –1 is dominated by the proton stretch vibrations of Zundel-like and intermediate geometries, broadened by the heterogeneity of aqueous configurations. There is a wide distribution of multidimensional potential shapes for the proton stretching vibration with varying degrees of potential asymmetry and confinement. In conclusion, we find specific cross peak patterns due to aqueous Zundel-like species. These studies provide clarity on highly debated spectral assignments and stringent spectroscopic benchmarks for future simulations.

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

Reconstruction of a broadband spectrum of Alfvénic fluctuations

Alfvénic fluctuations in the solar wind exhibit a high degree of velocities and magnetic field correlations consistent with Alfvén waves propagating away and toward the Sun. Two remarkable properties of these fluctuations are the tendencies to have either positive or negative magnetic helicity (–1 ≤ σ {sub m} ≤ +1) associated with either left- or right- topological handedness of the fluctuations and to have a constant magnetic field magnitude. This paper provides, for the first time, a theoretical framework for reconstructing both the magnetic and velocity field fluctuations with a divergence-free magnetic field, with any specified power spectral index and normalized magnetic- and cross-helicity spectrum field fluctuations for any plasma species. The spectrum is constructed in the Fourier domain by imposing two conditions—a divergence-free magnetic field and the preservation of the sense of magnetic helicity in both spaces—as well as using Parseval's theorem for the conservation of energy between configuration and Fourier spaces. Applications to the one-dimensional spatial Alfvénic propagation are presented. The theoretical construction is in agreement with typical time series and power spectra properties observed in the solar wind. The theoretical ideas presented in this spectral reconstruction provide a foundation for more realistic simulations of plasma waves, solar wind turbulence, and the propagation of energetic particles in such fluctuating fields.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Efficient Measurement-Driven Eigenenergy Estimation with Classical Shadows

Quantum algorithms exploiting real-time evolution under a target Hamiltonian have demonstrated remarkable efficiency in extracting key spectral information. However, the broader potential of these methods, particularly beyond ground-state calculations, is underexplored. In this work, we introduce the framework of multiobservable dynamic mode decomposition (MODMD), which combines the observable dynamic mode decomposition (DMD), a measurement-driven eigensolver tailored for near-term implementation, with classical shadow tomography. MODMD leverages random scrambling in the classical shadow technique to construct, with exponentially reduced resource requirements, a signal subspace that encodes rich spectral information. Notably, we replace typical Hadamard-test circuits with a protocol designed to predict low-rank observables, thereby broadening the use of classical shadow tomography for predicting many low-rank observables. We establish theoretical guarantees on the spectral approximation from MODMD, taking into account distinct sources of error. In the ideal case, we prove that the spectral error scales as exp (−Δ⁢𝐸⁢𝑡 max ), where Δ⁢𝐸 is the Hamiltonian spectral gap and 𝑡 max is the maximal simulation time. This analysis provides a rigorous justification of the rapid convergence observed across simulations. To demonstrate the utility of our framework, we consider its application to fundamental tasks, such as determining the low-lying, i.e., ground or excited, energies of representative many-body systems. Our work paves the path for efficient designs of measurement-driven algorithms on near-term and early fault-tolerant quantum devices.

quantum algorithms & computation↗
DOE OSTI2021

Random circuit block-encoded matrix and a proposal of quantum LINPACK benchmark

The LINPACK benchmark reports the performance of a computer for solving a system of linear equations with dense random matrices. Although this task was not designed with a real application directly in mind, the LINPACK benchmark has been used to define the list of TOP500 supercomputers since the debut of the list in 1993. We propose that a similar benchmark, called the quantum LINPACK benchmark, could be used to measure the whole machine performance of quantum computers. The success of the quantum LINPACK benchmark should be viewed as the minimal requirement for a quantum computer to perform a useful task of solving linear algebra problems, such as linear systems of equations. We propose an input model called the Random Circuit Block-Encoded Matrix (RACBEM), which is a proper generalization of a dense random matrix in the quantum setting. The RACBEM model is efficient to be implemented on a quantum computer and can be designed to optimally adapt to any given quantum architecture, with relying on a black-box quantum compiler. Besides solving linear systems, the RACBEM model can be used to perform a variety of linear algebra tasks relevant to many physical applications, such as computing spectral measures, time series generated by a Hamiltonian simulation, and thermal averages of the energy. We implement these linear algebra operations on IBM Q quantum devices as well as quantum virtual machines, and demonstrate their performance in solving scientific computing problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Isovector charges of the nucleon from 2+1-flavor QCD with clover fermions

Here, we present high-statistics estimates of the isovector charges of the nucleon from four 2 + 1-flavor ensembles generated using Wilson-clover fermions with stout smearing and tree-level tadpole improved Symanzik gauge action at lattice spacings 𝑎 ≈ 0.127 and 0.09 fm and with 𝑀 𝜋 ≈ 280 and 170 MeV. The truncated solver method with bias correction and the coherent source sequential propagator construction are used to cost-effectively achieve 𝑂⁡(1⁢0 5 ) measurements on each ensemble. Using these data, the analysis of two-point correlation functions is extended to include four states in the fits, and of three-point functions to three states. Control over excited-state contamination in the calculation of the nucleon mass, the mass gaps between excited states, and in the matrix elements is demonstrated by the consistency of estimates using this multistate analysis of the spectral decomposition of the correlation functions and from simulations of the three-point functions at multiple values of the source-sink separation. The results for all three charges, 𝑔 𝐴 , 𝑔 𝑆 , and 𝑔 𝑇 , are in good agreement with calculations done using the clover-on-HISQ lattice formulation with similar values of the lattice parameters.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

High-Efficiency Thermoelectric Single-Photon Detector Based on Lanthanum and Cerium Hexaborides

A design for a high-efficiency single-photon detector based on lanthanum and cerium hexaborides, which operates from the infrared to ultraviolet spectral ranges, is suggested. The results of computer simulations of heat transfer in the sensitive element of the detector upon the absorption of photons with energies of 0.5–4.13 eV are presented. To attain a high efficiency of the system of photon detection in the wavelength range from near infrared to ultraviolet, lanthanum hexaboride is proposed as the absorber and heat-sink material in the sensitive element. It is shown that a sensitive element of both single- and three-layer design made entirely of hexaborides will possess a gigahertz counting rate and a detection efficiency exceeding 90%.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Characterizing the Turbulent Structure of the Convective Boundary Layer Using ARM/ASR Observations and LES Model Output

This project will ultimately combine observations from Raman lidar, Doppler lidar, AERI, eddy correlation measurements, and airborne and ground-based high spectral resolution lidar systems with many large eddy simulations (LES) of realistic boundary layer cases to study: The diurnal evolution of the boundary layer thermodynamic structure; The turbulent structure of the CBL and the entrainment zone; The ability of LES models to simulate the thermodynamic and turbulent structure of the CBL; How perturbations in the surface properties and CBL properties influence boundary layer clouds; The similarities and differences in the boundary layer and cloud interaction in the marine vs. continental boundary layer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Technical Report on Subsurface Monitoring of the Brady Hot Spring Geothermal Site, Nevada, based upon Full Waveform Inversion

Abilities to accurately characterize the subsurface in a geothermal setting is key to assess and support production. An important element of geothermal reservoir monitoring is also the ability to investigate fluid transport within fracture network. This report focuses on improving subsurface imaging and monitoring in geothermal settings using full waveform inversion based on the adjoint method and time-lapse imaging. To assess our method, we rely on a dense seismic dataset collected in 2016 at the Brady Hot Springs geothermal site in Nevada for the DOE-funded project Poroelastic Tomography by Adjoint Inverse Modeling of Data from Seismology, Geodesy, and Hydrology. This dataset captures subsurface changes across four stages of geothermal power plant operations, which involve varying rates of fluid injection and extraction. Two velocity models were previously derived from this dataset using different methods: one based on travel times and another on sweep interferometry. Our first step is to refine these models using adjoint tomography, which has been applied successfully at global and regional-scales but is less common at the reservoir-scale. Two approaches are then explored for time-lapse analysis: directly comparing refined tomographic models from different stages or backpropagating waveform differences relative to a baseline tomographic model. The main take away is that both approaches highlight similar reservoir behaviors, but the latter approach is more computationally effective in capturing small-scale changes in subsurface properties. For this work, we leverage the use of Salvus (www.mondaic.com), an end-to-end seismic imaging solution, relying on the spectral element method to compute forward and adjoint simulations, and developed by Mondaic Ltd. It includes integrated workflow management that handles waveform and metadata, launches simulations, computes waveform misfits and adjoint sources, and iterates for model updates by nonlinear optimization.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Temporal Properties of Compressible Magnetohydrodynamic Turbulence

Describing the temporal properties of compressible magnetohydrodynamic (MHD) turbulence is a fundamental problem that has important implications for particle acceleration and transport in astrophysical plasmas. Here, by carefully analyzing the spatial and temporal properties of compressible MHD turbulence, we derive a new spectral power density function that is supported by simulations. This new function reveals that the low-frequency fluctuations are dominated by modes with small parallel wavenumbers with respect to the mean background magnetic field. Furthermore, for fluctuations with dynamically significant parallel wavenumbers, broadening around their eigenfrequencies is described by this function, which is in close agreement with simulations. We use this formalism to present the scaling properties of individual MHD modes. Such broadening is a direct consequence of nonlinear processes and is different for the three fundamental MHD modes. Our results provide a new window to investigate the temporal properties of turbulence and will enable further studies on the interaction between compressible MHD turbulence and energetic plasmas.

79 ASTRONOMY AND ASTROPHYSICS↗