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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Inhomogeneous mean-field approach to collective excitations near the superfluid–Mott glass transition

Highlights: • Unconventional collective excitations across the superfluid–Mott glass transition. • Higgs mode localized in superfluid and insulating phases. • Broad uncritical amplitude response function at quantum phase transition. • Goldstone mode delocalizes in the superfluid phase. We develop an inhomogeneous quantum mean-field approach to the behavior of collective excitations across the superfluid–Mott glass quantum phase transition in two dimensions, complementing recent quantum Monte Carlo simulations (Puschmann et al. 2020). In quadratic (Gaussian) approximation, the Goldstone (phase) and Higgs (amplitude) modes completely decouple. Each is described by a disordered Bogoliubov Hamiltonian which can be solved by an inhomogeneous multi-mode Bogoliubov transformation. We find that the Higgs mode is spatially localized in both phases. The corresponding scalar spectral function shows a broad peak that is noncritical in the sense that its peak frequency does not soften but remains nonzero across the quantum phase transition. In contrast, the lowest-energy Goldstone mode delocalizes in the superfluid phase, leading to a zero-frequency spectral peak. We compare these findings to the results of the quantum Monte Carlo simulations. We also relate them to general results on the localization of bosonic excitations, and we discuss the limits and generality of our approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

CFD Evaluation of Pressure Change Along Coolant Passages in Sodium-Cooled Fast Reactor with Nek5000

To support the design efforts of advanced sodium-cooled fast reactors (SFRs), a series of computational fluid dynamics (CFD) simulations are performed to investigate the pressure change along various flow passages in the proposed SFR system. The simulations are carried out with the state-of-the-art spectral element flow solver, Nek5000. Two specific case studies are presented in this paper: the flow exiting the axial neutron reflector channels and the flow entering the fuel pin bundle. Due to the high Reynolds numbers expected, a Reynolds-averaged Navier-Stokes (RANS) approach is necessary to model the turbulence. A newly developed regularized k-omega RANS model is adopted in the related CFD calculations. The first case study explores the effect of Reynolds number on the pressure change when flow exits the reflector channels. The pressure change in this case has two major contributors: the change due to wall friction and the Bernoulli effect. It is noted that the nondimensional pressure loss follows a log-linear trend up to Re = 10 5 , and then the trend is flattened. In the second case study, the advanced NekNek coupling capability is tested where an integral domain can be divided into multiple subdomains with coupling interfaces, which would greatly ease the meshing process of complex engineering geometries and potentially save computational resources. The preliminary results obtained so far confirm the consistency between the NekNek results and those produced by regular Nek5000 simulation. The presented work demonstrates the readiness and flexibility of the related CFD techniques, which is part of the broader effort to leverage cutting-edge CFD to inform the advanced nuclear reactor designs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Capturing the interactions between ice sheets, sea level and the solid Earth on a range of timescales: a new “time window” algorithm

Retreat and advance of ice sheets perturb the gravitational field, solid surface and rotation of the Earth, leading to spatially variable sea-level changes over a range of timescales $\textit{O}$ (10 0–6 years), which in turn feed back onto ice-sheet dynamics. Coupled ice-sheet–sea-level models have been developed to capture the interactive processes between ice sheets, sea level and the solid Earth, but it is computationally challenging to capture short-term interactions $\textit{O}$ (10 0–2 years) precisely within longer $\textit{O}$ (10 3–6 years) simulations. The standard forward sea-level modelling algorithm assigns a uniform temporal resolution in the sea-level model, causing a quadratic increase in total CPU time with the total number of input ice history steps, which increases with either the length or temporal resolution of the simulation. In this study, we introduce a new “time window” algorithm for 1D pseudo-spectral sea-level models based on the normal mode method that enables users to define the temporal resolution at which the ice loading history is captured during different time intervals before the current simulation time. Utilizing the time window, we assign a fine temporal resolution $\textit{O}$ (10 0–2 years) for the period of ongoing and recent history of surface ice and ocean loading changes and a coarser temporal resolution $\textit{O}$ (10 3–6 years) for earlier periods in the simulation. This reduces the total CPU time and memory required per model time step while maintaining the precision of the model results. We explore the sensitivity of sea-level model results to the model temporal resolution and show how this sensitivity feeds back onto ice-sheet dynamics in coupled modelling. We apply the new algorithm to simulate sea-level changes in response to global ice-sheet evolution over two glacial cycles and the rapid collapse of marine sectors of the West Antarctic Ice Sheet in the coming centuries and provide appropriate time window profiles for each application. The time window algorithm reduces the total CPU time by ~ 50 % in each of these examples and changes the trend of the total CPU time increase from quadratic to linear. This improvement would increase with longer simulations than those considered here. Our algorithm also allows for coupling time intervals of annual temporal scale for coupled ice-sheet–sea-level modelling of regions such as West Antarctica that are characterized by rapid solid Earth response to ice changes due to the thin lithosphere and low mantle viscosities.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Self-radiography of imploded shells on OMEGA based on additive-free multi-monochromatic continuum spectral analysis

Radiographs of pure-DT cryogenic imploding shells provide critical validation of progress toward ignition-scalable performance of inertial confinement fusion implosions. Cryogenic implosions on the OMEGA Laser System can be self-radiographed by their own core spectral emission near ≈2 keV. Utilizing the distinct spectral dependences of continuum emissivity and opacity, the projected optical-thickness distribution of imploded shells, i.e., the shell radiograph, can be distinguished from the structure of the core emission distribution in images.Importantly, this can be done without relying on spectral additives (shell dopants), as in previous applications of implosion self-radiography. Furthermore, demonstrations with simulated data show that this technique is remarkably well-suited to cryogenic implosions and can also be applied to self-radiography of imploded room-temperature CH shells at higher spectral energy (hv ≈ 3–5 keV) based on the very similar continuum spectrum of carbon. Experimental demonstration of additive-free self-radiography with warm CH shell implosions on OMEGA will provide an important proof of principle for future applications to cryogenic DT implosions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Inferring the temperature profile of the radiative shock in the COAX experiment with shock radiography, Dante, and spectral temperature diagnostics

Predicting and modeling the behavior of experiments with radiation waves propagating through low-density foams require a detailed quantification of the numerous uncertainties present. In regimes where a prominent radiative shock is produced, key dynamical features include the shock position, temperature, and curvature and the spatial distribution and temperature of the corresponding supersonic radiation wave. The COAX experimental platform is designed to constrain numerical models of such a radiative shock propagating through a low-density foam by employing radiography for spatial and shock information, Dante for characterizing the x-ray flux from the indirectly driven target, and a novel spectral diagnostic designed to probe the temperature profile of the wave. In this work, we model COAX with parameterized 2D simulations and a Hohlraum-laser modeling package to study uncertainties in diagnosing the experiment. The inferred temperature profile of the COAX radiation transport experiments has been shown to differ from simulations more than expected from drive uncertainties that have been constrained by simultaneous soft x-ray flux and radiography measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Spectral Element Applications in Complex Nuclear Reactor Geometries: Tet-to-Hex Meshing

The spectral element code Nek5000 is an open-source, higher-order computational fluid dynamics code developed at Argonne National Laboratory. It is designed to solve incompressible Navier-Stokes equations, but it also has a low-Mach-number approximation feature available. Large eddy simulation is approached by explicit filtering of the velocity field (and other fields) to mimic the effect of dissipation due to the unresolved scale. The computational domain is decomposed into second-order hexahedral elements that conform to the boundaries. However, generating a high-quality pure-hexahedral mesh can be challenging for some problems. For simple geometries, traditional blocking methods can be used to decompose the domain into smaller blocks to generate a so-called structural mesh. A structural mesh can maintain good orthogonality but can have a highly skewed mesh to conform to the geometry, as well as unnecessary refinement in the far field. Moreover, for geometries with relative complexity, blocking the geometry becomes impossible. To address these issues, we adopted a let-to-hex strategy to generate a pure hexahedral mesh for Nek5000. First, we generate a pure tetrahedral mesh for the geometry; then we divide one tetrahedral element into four hexahedral elements. A pure tetrahedral mesh could be easily generated for complex geometries by using many current meshing codes. In this paper, we use the commercial codes ANSYS meshing and ANSYSICEM to generate the pure tetrahedral mesh and then convert it to a pure hexahedral mesh. Boundary layers are extruded in ANSYSICEM to maintain near-wall resolution.

ANSYS ICEM↗
DOE OSTI2020

Comparison of Eulerian Bin and Lagrangian Particle-Based Schemes in Simulations of Pi Chamber Dynamics and Microphysics

This paper presents a comparison of simulations applying either a traditional Eulerian bin microphysics or a novel particle-based Lagrangian approach to represent CCN activation and cloud droplet growth. The Eulerian microphysics solve the evolution equation for the spectral density function, whereas the Lagrangian approach follows computational particles referred to as superdroplets. Each superdroplet represents a multiplicity of natural droplets that makes the Lagrangian approach computationally feasible. The two schemes apply identical representation of CCN activation and use the same droplet growth equation; these make direct comparison between the two schemes practical. The comparison, the first of its kind, applies an idealized simulation setup motivated by laboratory experiments with the Pi Chamber and previous model simulations of the Pi Chamber dynamics and microphysics. The Pi Chamber laboratory apparatus considers interactions between turbulence, CCN activation, and cloud droplet growth in moist Rayleigh–Bénard convection. Simulated steady-state droplet spectra averaged over the entire chamber are similar, with the mean droplet concentration, mean radius, and spectral width close in Eulerian and Lagrangian simulations. Small differences that do exist are explained by the inherent differences between the two schemes and their numerical implementation. The local droplet spectra differ substantially, again in agreement with the inherent limitations of the theoretical foundation behind each approach. There is a general agreement between simulations and Pi Chamber observations, with simplifications of the CCN activation and droplet growth equation used in the simulations likely explaining specific differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Dynamic control of the spatial frequency content of an intense laser via intra-beam energy transfer

We present theory and simulations demonstrating that bandwidth-inducing optical smoothing techniques such as smoothing by spectral dispersion (SSD) can initiate intra-beam energy transfer (IBET) within a single beam. IBET consists of energy transfer between the different frequency components of the beam and results in a simultaneous and correlated redistribution of the temporal and spatial frequencies of the laser's electric field. We identify a threshold in bandwidth for IBET to occur and validate the theory with simulations. The bandwidth and the color cycling of SSD can result in IBET in conditions relevant for high-energy density experiments at the National Ignition Facility and Omega Laser. We show that IBET could be used to manipulate a laser's spatial frequency content within a plasma. In particular, this technique could be used to increase a laser's effective f-number, with potential applications for self-guiding in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Machine learning aided line intensity ratio method for helium–hydrogen mixed recombining plasmas

The helium line intensity ratio (LIR) with the help of a collisional radiative (CR) model has long been used to measure the electron density, n e , and temperature, T e , and its potential and limitations for fusion applications have been discussed. However, it has been reported that the CR model approach leads to deviations in helium–hydrogen mixed plasmas and/or recombining plasmas. In this study, a machine learning (ML) aided LIR method is used to measure n e and T e from spectroscopic data of helium–hydrogen mixed recombining plasmas in the divertor simulator Magnum-PSI. To analyze mixed plasmas, which have more complex spectral shapes, the spectroscopy data were used directly for training instead of separating the intensities of each line. Finally, it is shown that the ML approach can provide a robust and simpler analysis method to deduce n e and T e from the visible emissions in helium–hydrogen mixed plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Influence of Electronic Polarization on the Spectral Density

Accurate spectral densities are necessary for computing realistic exciton dynamics and nonlinear optical spectra of chromophores in condensed-phase environments, including multichromophore pigment–protein systems. However, due to the significant computational cost of computing spectral densities from first principles, requiring many thousands of excited-state calculations, most simulations of realistic systems rely on treating the environment as fixed-point charges. In this work, using a number of representative systems ranging from solvated chromophores to the photoactive yellow protein (PYP), we demonstrate that the quantum mechanical (QM) electronic polarization of the environment is key to obtaining accurate spectral densities and line shapes within the cumulant framework. We show that the QM environment can enhance or depress the coupling of fast chromophore degrees of freedom to the energy gap, altering the electronic–vibrational coupling and the resulting vibronic progressions in the absorption spectrum. In analyzing the physical origin of peaks in the spectral density, we identify vibrational modes that couple the electron and the hole as being particularly sensitive to the QM screening of the environment. For PYP, we reveal the need for careful determination of the appropriate QM region to obtain reliable spectral densities. Our findings suggest that the QM polarization of the environment can be crucial not just for excitation energies but also for electronic–vibrational coupling in complex systems with implications for the correct modeling of linear and nonlinear optical spectroscopy in the condensed phase as well as energy transfer in pigment–protein complexes.

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

Performance assessment of near-fault buildings subjected to physics-based simulated earthquake ground motions with fling step

The effects of the co-seismic static offset (known as fling step) and associated velocity pulses on civil structures have been difficult to study because the static offset is typically removed during the processing of earthquake ground motion records. Simulated ground motions contain fling features and require no processing; therefore, they create new opportunities for representing fling features in seismic hazard analysis and assessing their influence on the seismic demands on near-fault structures. We use physics-based fault rupture simulations to study the characteristics of ground motions with fling step and the sensitivity of the near-fault structural demands to strong fling features. We uncover that simulated ground motions with a large fling step tend to have higher spectral intensity than those without a fling step at the same rupture distance, especially at periods longer than 2 s. As a result, the structural demands on flexible buildings tend to be the most sensitive to the fling features. Statistical analysis suggests that the ground motion spectral shape (represented by spectral accelerations at multiple periods) is—in most cases—a sufficient predictor of the structural demands on near-fault low-rise and mid-rise buildings at locations that are susceptible to strong fling effects. Finally, ground motion record selection experiments reveal that representing the spectral shape features at periods that are most relevant to a given structure may be an effective strategy to reduce the bias in the estimated demands on near-fault long-period structures when the available database of records is considered deficient in fling features.

Fling step↗
DOE OSTI2021

The Effect of Spectral Resolution of Light Sources on Photopic and a-opic Quantities

There is a growing interest in including light source spectrum in advanced lighting software tools and simulations. Given that traditional lighting software tools have been used primarily for calculating photometric quantities, their simplifying assumptions may not be suitable for calculating other a-opic quantities. Commonly used simulation tools use three values to represent the three primary colors (red, green, and blue), but others have expanded the number of bands (i.e., the spectral resolution) to 9, 27, or 81 in an attempt to more accurately represent variations in light source spectrum and downstream spectrally-derived quantities. It remains unclear, however, to what extent spectral resolution affects calculated quantities. To address this gap, a numerical analysis was completed using a large spectral power distribution database (n = 1,302 light sources). Calculated illuminance, a-opic irradiance, luminous efficacy of radiation (LER), melanopic efficacy of radiation (MER), and melanopic to photopic (M/P) ratio were compared for spectral resolutions of 3 and 9 bands compared to a baseline of 81 bands. Across all examined lighting quantities, considerable errors—a mean absolute percent error of 18.4%—were found when using 3 band calculations. These were reduced to 4% for 9 band calculations. The errors varied by metric and light source type. The results suggest that the use of 9 bands can more accurately characterize performance of different light source types across a range of metrics compared to using 3 bands.

Abboushi, Belal K.↗
DOE OSTI2025

Attosecond inner-shell lasing at ångström wavelengths

Since the invention of the laser, nonlinear effects such as filamentation, Rabi cycling and collective emission have been explored in the optical regime, leading to a wide range of scientific and industrial applications. X-ray free-electron lasers (XFELs) have extended many optical techniques to X-rays for their advantages of ångström-scale spatial resolution and elemental specificity. An example is XFEL-driven inner-shell Kα 1 (2p 3/2 → 1s 1/2 ) X-ray lasing in elements ranging from neon to copper, which has been used for nonlinear spectroscopy and development of new X-ray laser sources. Here, in this study, we show that strong lasing effects similar to those in the optical regime can occur at 1.5–2.1 Å wavelengths during high-intensity (>10 19 W cm −2 ) XFEL-driven Kα 1 lasing of copper and manganese. Depending on the temporal XFEL pump pulse substructure, the resulting X-ray pulses (about 10 6 −10 8 photons) can exhibit strong spatial inhomogeneities and spectral splitting, inhomogeneities and broadening. Three-dimensional Maxwell–Bloch calculations show that the observed spatial inhomogeneities result from X-ray filamentation and that the broad spectral features are driven by sub-femtosecond Rabi cycling. Our simulations indicate that these X-ray pulses can have pulse lengths of less than 100 attoseconds and coherence properties that provide opportunities for quantum X-ray optics applications.

47 OTHER INSTRUMENTATION↗
DOE OSTI2024

Simulations and models for the Richtmyer–Meshkov instability with broadband perturbations

The Richtmyer–Meshkov instability (RMI) is shock driven and affects many phenomena from inertial fusion to supernova explosions. The behavior of single-modes in the RMI has been studied extensively but less is known with the broadband perturbations that occur in applications. Here, we describe extensive numerical simulations and modeling of the RMI with broadband perturbations with an initial power spectrum of the form P ∝ km, where k is the wavenumber. The hydrodynamic condition is the same as in the θ-Group Collaboration [Thornber et al., “Late-time growth rate, mixing, and anisotropy in the multimode narrowband Richtmyer–Meshkov instability: The θ-group collaboration,” Phys. Fluids 29, 105107 (2017)] with a Mach ∼1.86 shock and gamma-law = 5/3 fluids with Atwood number A = −0.5. The bubble amplitude hB is found to grow in two stages. Initially, hB undergoes a phase-inversion and grows linearly in time (t) at a rate consistent with a linear theory. Asymptotically, hB grows as a power law ∼tθ when k|hB| > O(1), where k is the average wavenumber for the initial spectrum. The RMI behavior in simulations and modeling agree over a wide range of exponent m, spectral width, initial amplitude, and time. The agreement is quantified objectively using statistical analysis.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Population inversion X-ray laser oscillator

Oscillators are at the heart of optical lasers, providing stable, transform-limited pulses. Until now, laser oscillators have been available only in the infrared to visible and near-ultraviolet (UV) spectral region. In this paper, we present a study of an oscillator operating in the 5- to 12-keV photon-energy range. We show that, using the Kα 1 line of transition metal compounds as the gain medium, an X-ray free-electron laser as a periodic pump, and a Bragg crystal optical cavity, it is possible to build X-ray oscillators producing intense, fully coherent, transform-limited pulses. As an example, we consider the case of a copper nitrate gain medium generating ~ 5 × 10 10 photons per pulse with 37-fs pulse length and 48-meV spectral resolution at 8-keV photon energy. Our theoretical study and simulation of this system show that, because of the very large gain per pass, the oscillator saturates and reaches full coherence in four to six optical-cavity transits. We discuss the feasibility and design of the X-ray optical cavity and other parts of the oscillator needed for its realization, opening the way to extend X-ray–based research beyond current capabilities.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Delensing the CMB with the cosmic infrared background: the impact of foregrounds

ABSTRACT The most promising avenue for detecting primordial gravitational waves from cosmic inflation is through measurements of degree-scale cosmic microwave background (CMB) B-mode polarization. This approach must face the challenge posed by gravitational lensing of the CMB, which obscures the signal of interest. Fortunately, the lensing effects can be partially removed by combining high-resolution E-mode measurements with an estimate of the projected matter distribution. For near-future experiments, the best estimate of the latter will arise from co-adding internal reconstructions (derived from the CMB itself) with external tracers such as the cosmic infrared background (CIB). In this work, we characterize how foregrounds impact the delensing procedure when CIB intensity, I, is used as the matter tracer. We find that higher point functions of the CIB and Galactic dust such as 〈BEI〉c and 〈EIEI〉c can, in principle, bias the power spectrum of delensed B-modes. To quantify these, we first estimate the dust residuals in currently available CIB maps and upcoming, foreground-cleaned Simons Observatory CMB data. Then, using non-Gaussian simulations of Galactic dust – extrapolated to the relevant frequencies, assuming the spectral index of polarized dust emission to be fixed at the value determined by Planck – we show that the bias to any primordial signal is small compared to statistical errors for ground-based experiments, but might be significant for space-based experiments probing very large angular scales. However, mitigation techniques based on multifrequency cleaning appear to be very effective. We also show, by means of an analytical model, that the bias arising from the higher point functions of the CIB itself ought to be negligible.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Single Shot Characterization of High Transformer Ratio Wakefields in Nonlinear Plasma Acceleration

We report that plasma wakefields can enable very high accelerating gradients for frontier high energy particle accelerators, in excess of 10 GeV/m. To overcome limits on single stage acceleration, specially shaped drive beams can be used in both linear and nonlinear plasma wakefield accelerators (PWFA), to increase the transformer ratio, implying that the drive beam deceleration is minimized relative to acceleration obtained in the wake. In this Letter, we report the results of a nonlinear PWFA, high transformer ratio experiment using high-charge, longitudinally asymmetric drive beams in a plasma cell. An emittance exchange process is used to generate variable drive current profiles, in conjunction with a long (multiple plasma wavelength) witness beam. The witness beam is energy modulated by the wakefield, yielding a response that contains detailed spectral information in a single-shot measurement. Using these methods, we generate a variety of beam profiles and characterize the wakefields, directly observing transformer ratios up to R=7.8. Lastly, a spectrally based reconstruction technique, validated by 3D particle-in-cell simulations, is introduced to obtain the drive beam current profile from the decelerating wake data.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Relativistic Asymmetric Magnetic Reconnection

We derive basic scaling equations for relativistic magnetic reconnection in the general case of asymmetric inflow conditions and obtain predictions for the outflow Lorentz factor and the reconnection rate. Kinetic particle-in-cell simulations show that the outflow speeds as well as the nonthermal spectral index are constrained by the inflowing plasma with the weaker magnetic energy per particle, in agreement with the scaling predictions. These results are significant for understanding nonthermal emission from reconnection in magnetically dominated astrophysical systems, many of which may be asymmetric in nature. The results provide a quantitative approach for including asymmetry on reconnection in the relativistic regime.

79 ASTRONOMY AND ASTROPHYSICS↗