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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

The Role of Cation Coordination in the Electrical and Optical Properties of Amorphous Transparent Conducting Oxides

Amorphous oxide semiconductor materials have demonstrated numerous advantages without compromise of electrical properties as compared to their crystalline counterparts, yet understanding of the fundamental principles allowing this has remained elusive. To study the origins of enhanced optoelectronic properties, we apply high-throughput, combinatorial sputtering, structural and spectral mapping, and computationally intensive ab initio molecular dynamics simulations with density functional theory to a ternary, post-transition metal oxide system, namely, zinc tin oxide. The deposited thin films exhibit a high figure of merit, achieving carrier densities in the range of 1019 to 1020 cm–3 and carrier mobilities up to 35 cm2/Vs. These results highlight the role of local distortions and cation coordination in determining the microscopic origins of carrier generation and transport. In particular, we identify the strong likelihood of Sn undercoordination in both Zn-poor and Zn-rich phases leading to the high carrier concentrations observed. This not only diverges from the still widespread historical indictment of oxygen vacancies controlling carrier population in crystalline oxides but also provides a comprehensive framework to describe the unique structure–property relationships using specific structural and electronic descriptors in disordered phase materials.

CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS,M↗
DOE OSTI2021

Separating Physics and Dynamics Grids for Improved Computational Efficiency in Spectral Element Earth System Models

Previous studies have shown that atmospheric models with a spectral element grid can benefit from putting physics calculations on a relatively coarse finite volume grid. Here we demonstrate an alternative high-order, element-based mapping approach used to implement a quasi-equal-area, finite volume physics grid in E3SM. Unlike similar methods, the new method in E3SM requires topology data purely local to each spectral element, which trivially allows for regional mesh refinement. Simulations with physics grids defined by 2 × 2, 3 × 3, and 4 × 4 divisions of each element are shown to verify that the alternative physics grid does not qualitatively alter the model solution. The model performance is substantially affected by the reduction of physics columns when using the 2 × 2 grid, which can increase the throughput of physics calculations by roughly 60%–120% depending on whether the computational resources are configured to maximize throughput or efficiency. A pair of regionally refined cases are also shown to highlight the refinement capability.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Statistical data analysis of x-ray spectroscopy data enabled by neural network accelerated Bayesian inference

Bayesian inference applied to x-ray spectroscopy data analysis enables uncertainty quantification necessary to rigorously test theoretical models. However, when comparing to data, detailed atomic physics and radiation transfer calculations of x-ray emission from non-uniform plasma conditions are typically too slow to be performed in line with statistical sampling methods, such as Markov Chain Monte Carlo sampling. Furthermore, differences in transition energies and x-ray opacities often make direct comparisons between simulated and measured spectra unreliable. Here, we present a spectral decomposition method that allows for corrections to line positions and bound–bound opacities to best fit experimental data, with the goal of providing quantitative feedback to improve the underlying theoretical models and guide future experiments. In this work, we use a neural network (NN) surrogate model to replace spectral calculations of isobaric hot-spots created in Kr-doped implosions at the National Ignition Facility. The NN was trained on calculations of x-ray spectra using an isobaric hot-spot model post-processed with Cretin, a multi-species atomic kinetics and radiation code. The speedup provided by the NN model to generate x-ray emission spectra enables statistical analysis of parameterized models with sufficient detail to accurately represent the physical system and extract the plasma parameters of interest.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

An enrichment wall modeling framework for spectral element methods

In the present work, a first-of-its-kind enrichment wall-model is developed within the spectral element method (SEM) framework for large-eddy simulations (LES) of wall-bounded turbulent flows. The method augments the polynomial solution in the wall-adjacent elements with an analytical law-of-the-wall enrichment function representing the mean velocity near the wall. In the solution representation, this enrichment function captures the large gradients in the boundary layer, which allows the polynomial modes to represent the turbulent fluctuations. The enriched solution is able to resolve the shear stress at the wall without any modification to the no-slip wall boundary conditions, which allows for greater accuracy in the near-wall region compared to traditional methods. The enrichment wall modeling approach is implemented in a high-order SEM computational fluid dynamics solver, Nek5000, and its performance is assessed in turbulent channel flow wall-modeled LES for a range of Reynolds numbers. It is demonstrated that the enrichment wall-model improves solution accuracy on under-resolved near-wall grids as compared to traditional shear stress wall-models.

42 ENGINEERING↗
DOE OSTI2023

Spectral function of the $J_1$ – $J_2$ Heisenberg model on the triangular lattice

Spectral probes, such as neutron scattering, are crucial for characterizing excitations in quantum many-body systems and the properties of quantum materials. Among the most elusive phases of matter are quantum spin liquids, which have no long-range order even at zero temperature and host exotic fractionalized excitations with nontrivial statistics. These phases can occur in frustrated quantum magnets, such as the paradigmatic Heisenberg model with nearest- and next-nearest-neighbor exchange interactions on the triangular lattice, the so-called $J_1$ – $J_2$ model. In this work, we compute the spectral function using large-scale matrix product state simulations across the three different phases of this model's phase diagram, including a quantum spin-liquid phase at intermediate $J_2$/$J_1$. Despite a plethora of theoretical and experimental studies, the exact nature of this phase is still contested, with the dominant candidates being a gapped $\mathbb{Z}_2$, a gapless U(1) Dirac, and a spinon Fermi surface quantum spin-liquid state. We find a V-shaped spectrum near the center of the Brillouin zone ($Γ$ point), a key signature of a spinon Fermi surface, observed in prior neutron scattering experiments. However, we find a small gap near the $Γ$ point, ruling out such a phase. Furthermore, we find localized gapless excitations at the corner of the Brillouin zone boundary (K point) and the middle of the edge of the Brillouin zone boundary (M point), ruling out the gapped $\mathbb{Z}_2$ spin-liquid phase. Here, our results imply that the intermediate spin-liquid phase is a gapless $U(1)$ Dirac spin liquid, and provide clear signatures to detect this phase in future neutron scattering experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Startup regime of high-efficiency tapering-enhanced FEL oscillator

In this paper, we present a design of a high-efficiency high-gain free-electron laser oscillator based on the use of a strongly tapered undulator for extracting energy from high-brightness electron beams. We provide an analytical model of the setup followed by numerical simulations for lasing at the wavelength of 13.5 nm. We discuss the optimization of the system in steady state and the conditions necessary for the pass-per-pass buildup of the power from shot noise level. We propose the use of fast phase shifters as a way to accelerate the buildup. Finally, we present time-dependent simulations of the oscillator and discuss the role of spectral filtering. The optimized working point yields a total energy conversion efficiency from the electron beam to output radiation above 1% at the wavelength of 13.5 nm.

Beam dynamics↗
DOE OSTI2023

Real-time infrared spectroscopy coupled with blind source separation for nuclear waste process monitoring

On-line infrared absorbance spectroscopy enables rapid measurement of solution-phase molecular species. Many spectra-to-concentration models exist for spectral data, with some models able to handle overlapping spectral bands and nonlinearities. However, model accuracy is limited by the quality of training data used in model fitting. The process spectra of nuclear waste simulants at the Savannah River Site display incongruity between training and process spectra; the glycolate spectral signature in the training data does not match the glycolate signature in Savannah River National Laboratory process data. A novel blind source separation algorithm is proposed that preprocesses spectral data so that process spectra more closely resemble training spectra, thereby improving model quantification accuracy when unexpected sources of variation appear in process spectra. The novel blind source separation preprocessing algorithm is shown to improve nitrate quantification from an R 2 of 0.934 to 0.988 and from 0.267 to 0.978 in two instances analyzing nuclear waste simulants from the Slurry Receipt Adjustment Tank and Slurry Mix Evaporator cycle at the Savannah River Site.

Crouse, Steven H.↗
DOE OSTI2021

Thermal effects in molecular gas-filled hollow-core fibers

Few-cycle sources with high average powers are required for applications to attosecond science. Raman-enhanced spectral broadening of Yb-doped laser amplifiers in molecular gases can yield few-cycle pulses, but thermal excitation of vibrational and rotational degrees of freedom may preclude high-power operation. Here we investigate changes in the spectral broadening associated with repetitive laser interactions in an N 2 O -filled hollow-core fiber. By comparing experimental measurements of the spectrum associated with each laser pulse to simulations based on a density matrix model, we find that losses in a spectral bandwidth and transmission are largely dominated by thermal excitation of the gas.

Beetar, John E. (ORCID:000000032885852X)↗
DOE OSTI2025

Resolving turbulent magnetohydrodynamics: a hybrid operator-diffusion framework

We present a hybrid machine learning framework that combines physics-informed neural operators (PINOs) with score-based generative diffusion models to simulate the full spatio-temporal evolution of two-dimensional, incompressible, resistive magnetohydrodynamic turbulence across a broad range of Reynolds numbers (Re). The framework leverages the equation-constrained generalization capabilities of PINOs to predict coherent, low-frequency dynamics, while a conditional diffusion model stochastically corrects high-frequency residuals, enabling accurate modeling of fully developed turbulence. Trained on a comprehensive ensemble of high-fidelity simulations with Re ϵ {100, 250, 500, 750, 1000, 3000, 10000}, the approach achieves state-of-the-art accuracy in regimes previously inaccessible to deterministic surrogates. At Re = 1000 and 3000, the model faithfully reconstructs the full spectral energy distributions of both velocity and magnetic fields late into the simulation, capturing non-Gaussian statistics, intermittent structures, and cross-field correlations with high fidelity. At extreme turbulence levels (Re = 10 000), it remains the first surrogate capable of recovering the high-wavenumber evolution of the magnetic field, preserving large-scale morphology and enabling statistically meaningful predictions.

Diffusion-Integrated Neural Operators↗
DOE OSTI2020

Understanding cloud droplet spectral dispersion effect using empirical and semi-analytical parameterizations in NCAR CAM5.3

Five parameterizations of cloud droplet spectral shape are implemented in a global climate model to investigate the dispersion effect and aerosol indirect effect (AIE). We design a series of experiments by modifying the microphysical cloud scheme of NCAR CAM5.3 (Community Atmospheric Model Version 5.3). We employ four empirical (Martin94, RLiu03, PengL03, Liu08) and one semi-analytical (LiuLi15) expressions for cloud droplet spectral shape parameters. Analysis focuses on the instantaneous differences in the simulated cloud microphysical properties and the comparison between model output and satellite data. The results show that RLiu03, PengL03 and LiuLi15 produce wider droplet spectrum and faster autoconversion rate, but Liu08 has a narrower droplet spectrum and slower autoconversion rate than the default parameterization (Martin94) in CAM5.3. Global dispersion effects caused by the five parameterizations modify the aerosol indirect effect by -10% (counteract) to 13% (strengthen). The simulated AIEs and dispersion effects exhibit noticeably spatial inhomogeneity. In the sensitive regions of AIE (Southeast Asia, North Pacific and west coast of South America), we decompose the response of shortwave cloud forcing to the change in droplet number for analysis. The varying dispersion effects can be explained by different responses of cloud properties in different spectral parameterizations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Theoretical analysis of the role of complex transition dipole phase in XUV transient-absorption probing of charge migration

We theoretically investigate the role of complex dipole phase in the attosecond probing of charge migration. The iodobromoacetylene ion (ICCBr + ) is considered as an example, in which one can probe charge migration by accessing both the iodine and bromine ends of the molecule with different spectral windows of an extreme-ultraviolet (XUV) pulse. The analytical expression for transient absorption shows that the site-specific information of charge migration is encoded in the complex phase of cross dipole products for XUV transitions between the I-4 d and Br-3 d spectral windows. Ab-initio quantum chemistry calculations on ICCBr + reveal that there is a constant π phase difference between the I-4 d and Br-3 d transient-absorption spectral windows, irrespective of the fine-structure energy splittings. Transient absorption spectra are simulated with a multistate model including the complex dipole phase, and the results correctly reconstruct the charge-migration dynamics via the quantum beats in the two element spectral windows, exhibiting out-of-phase oscillations.

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

The Concept of a Quantum Edge Simulator: Edge Computing and Sensing in the Quantum Era

Sensors, enabling observations across vast spatial, spectral, and temporal scales, are major data generators for information technology (IT). Processing, storing, and communicating this ever-growing amount of data pose challenges for the current IT infrastructure. Edge computing—an emerging paradigm to overcome the shortcomings of cloud-based computing—could address these challenges. Furthermore, emerging technologies such as quantum computing, quantum sensing, and quantum communications have the potential to fill the performance gaps left by their classical counterparts. Here, we present the concept of an edge quantum computing (EQC) simulator—a platform for designing the next generation of edge computing applications. An EQC simulator is envisioned to integrate elements from both quantum technologies and edge computing to allow studies of quantum edge applications. The presented concept is motivated by the increasing demand for more sensitive and precise sensors that can operate faster at lower power consumption, generating both larger and denser datasets. These demands may be fulfilled with edge quantum sensor networks. Envisioning the EQC era, we present our view on how such a scenario may be amenable to quantification and design. Given the cost and complexity of quantum systems, constructing physical prototypes to explore design and optimization spaces is not sustainable, necessitating EQC infrastructure and component simulators to aid in co-design. We discuss what such a simulator may entail and possible use cases that invoke quantum computing at the edge integrated with new sensor infrastructures.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗
DOE OSTI2024

Simultaneous reproduction of experimental profiles, fluxes, transport coefficients, and turbulence characteristics via nonlinear gyrokinetic profile predictions in a DIII-D ITER similar shape plasma

Experimental conditions obtained on the DIII-D tokamak in the ITER Similar Shape (ISS) have been compared extensively with nonlinear gyrokinetic simulation using the CGYRO code [Candy et al., J. Comput. Phys. 324, 73–93 (2016)] with comparisons spanning ion and electron heat fluxes, electron and impurity particle transport, and turbulent fluctuation levels and characteristics. Bayesian optimization techniques [Rodriguez-Fernandez et al., Nucl. Fusion 62(7), 076036 (2022)], combined with nonlinear gyrokinetics, have been used to obtain simultaneously Qi, Qe, and Γe flux-matched profiles that are found to be in good agreement with experimental profile measurements. Synthetic diagnostics were used to compare measured beam emission spectroscopy and correlation electron cyclotron emission turbulent fluctuations with nonlinear simulation. Although some disagreements exist, nonlinear simulations are found to be in generally good agreement with measured fluctuation levels, spectral shapes, and measured radial trends in low-k δne/ne and δTe/Te. Low (Li and C) and mid-Z (Ca) impurity transport was also compared with these flux-matched simulations. Fully stripped, low-Z impurities are well reproduced by the gyrokinetic modeling while clear disagreement exists in comparisons with mid-Z impurities. Nonlinear gyrokinetic investigation into the Z dependence of impurity transport in the ISS conditions is also performed, demonstrating clear trends of impurity diffusion with impurity Z (both D∝Z and D∝1/Z) that vary with the radial location studied. These trends are shown to result from the local dominance of ion temperature gradient or ∇n driven trapped electron mode turbulence and may contribute to the disagreement between simulation and experiment in mid-Z impurity transport. The results of this work represent one of the most complete validation studies of the gyrokinetic model performed to date and provide an example of new capabilities for predicting performance in future fusion devices.

Physics↗
DOE OSTI2026

Improving Spectral Resolution from Real-time Evolution for Correlated Systems

Abstract The quality of numerically simulated spectra using real-time evolution methods for strongly correlated systems is affected by both the length of simulation time and the system size, limiting resolution in both frequency and momentum. In this work, we propose a computationally cheap, linear autoregressive machine learning-based framework to extend short-time and short-distance results over a wider range. We use the proposed method to extend the lesser Green’s function for both the Hubbard model and the much more computationally challenging Hubbard-extended Holstein model. This technique significantly improves both the frequency and momentum resolution of the single-particle removal spectrum $${\mathcal{A}}(k,\omega )$$ A ( k , ω ) , allowing the observation of otherwise obscured spectral features due to electron-phonon coupling.

Tang, Ta↗
DOE OSTI2025

A discontinuous Galerkin spectral element method for compressible reacting flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large-eddy simulations because of their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reacting Navier-Stokes equations. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of the DG approach. The framework, implemented in the spectral element code Nek5000, is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. An entropy-residual based artificial viscosity is added to smooth shocked regions of flow, and a positivity-preserving limiter is implemented to suppress non-physical oscillations. These enhancements support the numerical stability of the hydrodynamic sub-step, which is decoupled from the chemistry integration through a second-order operator splitting method. Here, a series of smooth and discontinuous validation cases are presented in increasing physical and computational complexity for both inviscid and viscous flows. In particular, simulations of canonical one-dimensional and two-dimensional detonations are performed, and the high-order numerical results are validated against available literature data. Additional validation studies are carried out for classical three-dimensional numerical simulations of incompressible and compressible turbulent flows.

Compressible reacting flows↗
DOE OSTI2026

Stochastic equilibrium Raman spectroscopy (STERS)

In this manuscript, we propose a new method for cavity- and surface-enhanced Raman spectroscopy (SERS) with improved temporal resolution in the measurement of stochastic Raman spectral fluctuations. Our approach combines Fourier spectroscopy and photon correlation to decouple the integration time from the temporal resolution. Using statistical optics Monte Carlo simulations, we establish the relationship between time resolution and Raman signal strength, revealing that typical Raman spectral fluctuations, commensurate with molecular conformational dynamics, can theoretically be resolved on micro- to millisecond timescales. The method can further extract average single-molecule dynamics from small sub-ensembles, thereby potentially mitigating challenges in achieving strictly single-molecule isolation on SERS substrates.

Cobb-Bruno, Colburn [University of California, Ber↗
DOE OSTI2025

Angular-spatial hp -adaptivity for radiative transfer with discontinuous Galerkin spectral element methods

Radiative transfer is important for many science and engineering applications, and numerical simulations of radiative transfer can be challenging. For instance, the radiation field is seven-dimensional – three spatial, two angular, one wavelength, and one temporal – and often features steep gradients. Therefore, memory usage is a key issue. To reduce memory, some past work has investigated the use of adaptive mesh refinement (AMR), typically for either the spatial or angular coordinate, and typically for only h -adaptivity. Here, we propose the use of AMR for the spatial and angular coordinates together, and the use of h - and p -adaptivity together as hp -AMR for the potential for further memory savings. We implemented the proposed method for several test cases in two spatial and one angular dimension, with the discontinuous Galerkin spectral element method. These test cases featured highly anisotropic angular radiation, with or without steep spatial gradients. Our primary findings from these test cases were: (1) Angular hp -adaptivity can deliver the radiation solution with the same accuracy as, and with much less computational memory than, uniform angular h - or p -refinements, or angular h -adaptivity alone. This is most obvious when the incoming radiation is highly anisotropic, in which case the savings can be orders of magnitude. (2) Full spatial-angular hp -adaptivity is more efficient in solution representation, compared to solely spatial or solely angular -adaptivity. This is most evident when steep gradients are present in both the spatial and angular distribution. These results suggest that adaptive spatial- hp angular-refinement may perform well in large-scale seven-dimensional applications.

Adaptive refinement↗