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 325 records · Page 18

`SkyPy`: A package for modelling the Universe

SkyPy is an open-source Python package for simulating the astrophysical sky. It comprises a library of physical and empirical models across a range of observables and a command-line script to run end-to-end simulations. The library provides functions that sample realisations of sources and their associated properties from probability distributions. Simulation pipelines are constructed from these models using a YAML-based configuration syntax, while task scheduling and data dependencies are handled internally and the modular design allows users to interface with external software. SkyPy is developed and maintained by a diverse community of domain experts with a focus on software sustainability and interoperability. By fostering development, it provides a framework for correlated simulations of a range of cosmological probes including galaxy populations, large scale structure, the cosmic microwave background, supernovae and gravitational waves. Version 0.4 implements functions that model various properties of galaxies including luminosity functions, redshift distributions and optical photometry from spectral energy distribution templates. Future releases will provide additional modules, for example, to simulate populations of dark matter halos and model the galaxy-halo connection, making use of existing software packages from the astrophysics community where appropriate.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

QTRAJ 1.0: A Lindblad equation solver for heavy-quarkonium dynamics

Here we introduce an open-source package called QTraj that solves the Lindblad equation for heavy-quarkonium dynamics using the quantum trajectories algorithm. The package allows users to simulate the suppression of heavy-quarkonium states using externally-supplied input from 3+1D hydrodynamics simulations. The code uses a split-step pseudo-spectral method for updating the wave-function between jumps, which is implemented using the open-source multi-threaded FFTW3 package. This allows one to have manifestly unitary evolution when using real-valued potentials. In this paper, we provide detailed documentation of QTraj 1.0, installation instructions, and present various tests and benchmarks of the code.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A holey cavity for single-transducer 3D ultrasound imaging with physical optimization

Within the compressive sensing (CS) framework, one effective way to increase the likelihood of successful signal reconstruction is to employ random processes in the construction of the sensing matrix. This work presents a 3D holey cavity, with diverse frequency modes, to spectrally code, that is, randomize, the ultrasound wave fields. The simulated results show that the use of such a cavity enables imaging simple or complex targets, such as spheres or the letter E, by only a single transceiver—something that is not possible without the use of a coding structure like the cavity. The effect of noise on imaging results and the size of the targets on the first-order Born approximation (BA) are also investigated. Moreover, this study attempts to optimize the cavity, based on a single numerical metric, such as the sum of singular values (SSV) or mutual coherence (MC). Yet, it will be shown that neither of these metrics can consistently compare the norm-one imaging performance between two cavities of different materials or hole sizes. This leaves finding a quantitative metric for these kinds of optimizations an open problem.

42 ENGINEERING↗

Investigating the Electronic Structure of Prospective Water-Splitting Oxide BaCe 0.25 Mn 0.75 O 3-δ before and after Thermal Reduction

BaCe 0.25 Mn 0.75 O 3-δ (BCM), a non-stoichiometric oxide with a layered perovskite-like crystal structure, has recently emerged as a prospective contender for application in renewable energy harvesting by solar thermochemical hydrogen generation. Using solar-thermal energy and a reducing environment, oxygen vacancies can be created in high-temperature BCM, and the reduced crystal so obtained can, in turn, produce H 2 by stripping oxygen from H 2 O. Therefore, a first step toward understanding the working mechanism and optimizing the performance of BCM is a thorough and comparative analysis of the electronic structure of the pristine and the reduced material. In this paper, we probe the electronic structure of BCM using the combined effort of first-principles calculations and experimental O K-edge X-ray absorption spectroscopy (XAS). The computed projected density of states (PDOS) and orbital plots are used to propose a simplified model for orbital mixing between the oxygen and metal atoms. With the help of state-of-the-art simulations, we are able to find the origins of the XAS peaks and categorize them on the basis of contribution from Ce and Mn. For the reduced crystal, the calculations show that the change in electron density resulting from the reduction is strongly localized around the oxygen vacancy. Experimental measurements reveal a marked lowering of the first O K-edge peak in the reduced crystal. Using theoretical analysis, this is shown to result from lifting of spin degeneracy in the absorption peaks as well as from a diminished O 2p contribution to the frontier unoccupied orbitals, in accordance with the tight binding scheme. The simulated results serve as a reference for the extent of spectral change as a function of the percentage of oxygen vacancies in the reduced crystal. Here, our study paves the way for the investigation of the working mechanism of BCM and for computational and experimental efforts aimed at design and discovery of efficient water-splitting oxides.

08 HYDROGEN↗

A 2D chiral microcavity based on apparent circular dichroism

Abstract Engineering asymmetric transmission between left-handed and right-handed circularly polarized light in planar Fabry–Pérot (FP) microcavities would enable a variety of chiral light-matter phenomena, with applications in spintronics, polaritonics, and chiral lasing. Such symmetry breaking, however, generally requires Faraday rotators or nanofabricated polarization-preserving mirrors. We present a simple solution requiring no nanofabrication to induce asymmetric transmission in FP microcavities, preserving low mode volumes by embedding organic thin films exhibiting apparent circular dichroism (ACD); an optical phenomenon based on 2D chirality. Importantly, ACD interactions are opposite for counter-propagating light. Consequently, we demonstrated asymmetric transmission of cavity modes over an order of magnitude larger than that of the isolated thin film. Through circular dichroism spectroscopy, Mueller matrix ellipsometry, and simulation using theoretical scattering matrix methods, we characterize the spatial, spectral, and angular chiroptical responses of this 2D chiral microcavity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Near-infrared C-term MCD spectroscopy of octahedral uranium (V) complexes

C-term magnetic circular dichroism (MCD) spectroscopy is a powerful method for probing d–d and f–f transitions in paramagnetic metal complexes. However, this technique remains underdeveloped both experimentally and theoretically for studies of U(V) complexes of Oh symmetry, which have been of longstanding interest for probing electronic structure, bonding, and covalency in 5f systems. In this study, C-term NIR MCD of the Laporte forbidden f–f transitions of [UCl 6 ] - and [UF 6 ] - are reported, demonstrating the significant fine structure resolution possible with this technique including for the low energy Γ 7 → Γ 8 transitions in [UF 6 ] - . The experimental NIR MCD studies were further extended to [U(OC 6 F 5 ) 6 ] - , [U(CH 2 SiMe 3 ) 6 ] - , and [U(NC(tBu)(Ph)) 6 ] - to evaluate the effects of ligand-type on the f–f MCD fine structure features. Theoretical calculations were conducted to determine the Laporte forbidden f–f transitions and their MCD intensity experimentally observed in the NIR spectra of the U(V) hexahalide complexes, via the inclusion of vibronic coupling, to better understand the underlying spectral fine structure features for these complexes. These spectra and simulations provide an important platform for the application of MCD spectroscopy to this widely studied class of U(V) complexes and identify areas for continued theoretical development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A novel approach for calculating galaxy rotation curves using spaxel cross-correlation and iterative smoothing

ABSTRACT Precise measurements of the internal dynamics of galaxies have proven of great importance for understanding the internal dark matter distribution of galaxies. We present a novel method for measuring the line-of-sight (LOS) velocities across the face of galaxies by cross-correlation of spectral pixels (spaxels) and an iterative method of smoothing. On simulated data the method can accurately recover the input LOS velocities for different types of spectra (absorption-line dominated, emission-line dominated, and differing shapes of the continuum), and can handle stellar population radial gradients. Most important of all, it continues to provide reliable measurements of LOS velocities with reasonable uncertainties even when the spectra are very low signal-to-noise ratio (approaching ∼1), which is a challenge for traditional template-fitting approaches. We apply our method to data from a real MaNGA galaxy as a demonstration and find promising results with good precision. This novel approach can be complementary to existing methods primarily based on template fitting.

79 ASTRONOMY AND ASTROPHYSICS↗

Shakeup and shakeoff spectra in the electron capture decay of atomic 7 Be

The most stringent laboratory-based experimental limits on the existence of submegaelectronvolt sterile neutrinos are currently set by decay spectroscopy of radioactive 7 Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shakeup and shakeoff spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wave functions ab initio and compute all single and double shake processes in the electron capture decay of atomic 7 Be . The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the 7 Be is embedded into. The new models also show that the 𝐿/𝐾 electron capture ratio of 7 Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Intergalactic magnetic field studies by means of γ -ray emission from GRB 190114C

The presence of delayed GeV emission after a strong transient, such as a gamma-ray burst (GRB), in the very-high energy (VHE) E > 100 GeV band can be the signature of a nonzero magnetic field in the intergalactic medium. Here we used a synchrotron self-Compton multiwavelength model to infer an analytical description of the intrinsic VHE spectrum (corrected for absorption by the extragalactic background light) of GRB 190114C to predict the light curves and spectral energy distributions of the delayed emission with Monte Carlo simulations for different intergalactic magnetic field (IGMF) configurations (strengths B = 8 × 10 -21 G, 10 -20 G, 3 × 10 -20 G and correlation length $\lambda$ > 1 Mpc), and compared them with the Fermi Large Area Telescope (Fermi LAT) limits computed for several exposure times. We found that the Fermi LAT is not sensitive enough to constrain any IGMF strengths using GRB 190114C.

79 ASTRONOMY AND ASTROPHYSICS↗

Dynamic Phase Alignment in Navier-Stokes Turbulence

In Navier-Stokes turbulence, energy and helicity injected at large scales are subject to a joint direct cascade, with both quantities exhibiting a spectral scaling ∝k -5/3 . We demonstrate via direct numerical simulations that the two cascades are compatible due to the existence of a strong scale-dependent phase alignment between velocity and vorticity fluctuations, with the phase alignment angle scaling as cosα k ∝k -1 .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Spectral Rear Irradiance Testing and Modeling for Degradation and Performance of Solar Fields

This work investigates how the spectrum of irradiance incident on the rear of solar modules impacts the degradation and performance of backsheets. We model the spectral irradiance incident on the rear of modules through raytrace simulations and validate with measured field data collected from a 75kW single axis-tracked bifacial test site. A generic equation to estimate relative degradation is proposed, and we show that current acceleration factors for UV damage in chambers can be sub-estimate up to 4.5% absolute from the usually assumed 10% dosage on the rear surfaces.

bifacial performance↗

Rupture Jumping and Seismic Complexity in Models of Earthquake Cycles for Fault Stepovers with Off-Fault Plasticity

ABSTRACT Fault stepovers are prime examples of geometric complexity in natural fault zones that may affect seismic hazard by determining whether an earthquake rupture continues propagating or abruptly stops. However, the long-term pattern of seismicity near-fault stepovers and underlying mechanisms of rupture jumping in the context of earthquake cycles are rarely studied. Leveraging a hybrid numerical scheme combining the finite element and the spectral boundary integral methods, FEBE, we carry out fully dynamic simulations of sequences of earthquakes and aseismic slip for both compressive and tensile stepovers with off-fault plasticity. We consider a rate-and-state friction law for the fault friction and pressure-sensitive Drucker–Prager plasticity for the off-fault bulk response. We observe that the accumulation of plastic deformation, an indication of off-fault damage, is significantly different in the two cases, with more plastic deformation projected in the overlapping region for the tensile stepover. The seismic pattern for a tensile stepover is more complex than for a compressive stepover, and incorporating plasticity also increases complexity, relative to the elastic case. A tensile stepover with off-fault plasticity shows rupture segmentation, temporal clustering, and frequent rupture jumping from one fault to another. These results shed light on possible mechanisms of rupture jumping in fault stepovers as well as the long-term evolution of the fault zone.

Geochemistry & Geophysics↗

Assessing Impacts of Waves on Hub-Height Winds off the U.S. West Coast Using Lidar Buoys and Coupled Modeling Approaches

Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.

17 WIND ENERGY↗

X-ray Spectral-Timing Pipeline to Investigate an Electron-Scattering Time Delay in Black Hole Accretion Disks [Slides]

The soft lag in black hole X-ray binaries (BHXRBs) refers to the time delay for the soft, thermal disk to respond to hard, variable coronal irradiation. This time lag increases from less than 1 ms to ~10 ms during the hard-to-soft state transition. Interpretations of soft lag trends appeal to changing the light-travel path via an evolving coronal height and/or inner accretion disk radius. Both interpretations neglect a time delay contribution from the reprocessing of irradiation inside the disk, where electron-scattering opacity dominates. A new theory considering a thermalization (electron scattering) time delay in the disk can plausibly produce ~10 ms time delays in the intermediate state. To further investigate the impacts of adding a time delay component from electron scattering in the disk, we are creating a spectral-timing pipeline that can analyze NICER (Neutron Star Interior Composition Explorer) X-ray observations of BHXRBs in outburst. In the near future, we will apply this spectral timing pipeline to develop a reverberation lag model from simulations that include the thermalization time delay from electron scattering in the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

Shake-up and shake-off spectra in the electron capture decay of atomic $^7$Be

The most stringent laboratory-based experimental limits on the existence of sub-MeV sterile neutrinos are currently set by decay spectroscopy of radioactive $^7$Be embedded into superconducting sensors. The systematic uncertainties are dominated by the modeling of the electron shake-up and shake-off spectra that are not based on state-of-the-art atomic theory and do not include electron correlations or relativistic effects. We have used the multiconfiguration Dirac-Fock formalism to obtain correlated wavefunctions ab initio and compute all single and double shake processes in the electron capture decay of atomic $^7$Be. The simulations can explain some but not all of the observed spectral features, likely because the wave functions are modified by the Ta sensor material that the $^7$Be is embedded into. The new models also show that the L/K electron capture ratio of $^7$Be in Ta has previously been slightly underestimated revising the previous value of 0.070(7) to a new value of 0.0756(20).

Atomic Physics (physics.atom-ph)↗

Improving the modeling of near-wall interphase heat transfer in porous media models of Pebble Bed Reactors

Here, this work aims to improve capabilities for modeling localized effects in porous media models of Pebble Bed Reactors. The wall-channeling effect is the primary local phenomenon of interest in a PBR, where the presence of the reflector wall disrupts the pebble packing, causing the pebbles near the wall to pack less efficiently and creating large void regions. Accurate modeling of the near-wall region is important as it will affect core bypass flow and temperature predictions. Porous media models are commonly used for design scoping and plant-level simulations of PBRs. Although these models have some capabilities to model the near-wall region, the correlations that are available in porous media codes are often inaccurate when a multi-region model is used to discretize the near-wall region. This work employs a high-to-low analysis to study the accuracy of available interphase heat transfer closures. NekRS, a spectral element computational fluid dynamics code, is used to perform Large Eddy Simulations. These LES simulation results are compared to porous media model results from the Pronghorn porous media code. The friction term of the KTA drag closure is first improved, reducing the error in the prediction of the near-wall velocity from over 50% to less than 5%. This is combined with improvements to the form term from previous works to produce a drag closure that is capable of accurately modeling the wall-channeling effect across a variety of flow conditions. The Nusselt number predictions of several heat transfer correlations are compared to the high-fidelity results where it is found that the KTA heat transfer correlation is capable of accurately predicting the local Nusselt numbers that were determined in the high-fidelity simulation. Comparison of the radial solid temperature profiles, however, reveal discrepancies between NekRS and Pronghorn. It is discovered that the implementation of the interphase heat transfer coefficient that exists in many current porous media codes is not valid when local porosities are modeled. Instead, it is suggested that the interphase heat transfer coefficient should be dependent on the local porosity, the Nusselt number, and the local solid surface-to-volume ratio. Implementation of this change produces improvement in the agreement between the results obtained by NekRS and Pronghorn while using the KTA heat transfer correlation.

interphase heat transfer↗

Error Bounds for Dynamical Spectral Estimation

Dynamical spectral estimation is a well-established numerical approach for estimating eigenvalues and eigenfunctions of the Markov transition operator from trajectory data. Although the approach has been widely applied in biomolecular simulations, its error properties remain poorly understood. Here we analyze the error of a dynamical spectral estimation method called “the variational approach to conformational dynamics" (VAC). We bound the approximation error and estimation error for VAC estimates. Our analysis establishes VAC's convergence properties and suggests new strategies for tuning VAC to improve accuracy.

97 MATHEMATICS AND COMPUTING↗

Consideration of temperature-dependent emissivity of selective emitters in thermophotovoltaic systems

Spectral emissivity control is paramount for designing a high-efficiency selective emitter surface required for thermophotovoltaic (TPV) applications. Owing to the temperature dependency of materials optical constants, the spectral properties of a selective emitter surface changes with the emitter temperature. This paper presents the fabrication of a multilayer metal-dielectric ( ${\rm Si}_{3}{\rm N}_{4}/\rm W/Si_{3}{\rm N_{4}}$ ) coated tungsten selective emitter aimed for GaSb-based TPV systems and studies the dependence of its surface spectral emissivity, $\varepsilon (\lambda)$ , upon a temperature ranging from 300 K to 1500 K. Both the simulation and experimental methods were used to characterize $\varepsilon (\lambda)$ as a function of temperature. For wavelengths less than 1.4 µm, $\varepsilon (\lambda)$ was found to have a minimal dependence on temperature. Beyond 1.4 µm, $\varepsilon (\lambda)$ increases with the temperature. At 1.55 µm, the simulation and experimental data estimated a ${\sim}{{4}}\%$ greater emissivity at 1500 K than at room temperature. At 1500 K, the increased $\varepsilon (\lambda)$ at longer wavelengths lowered the spectral conversion efficiency of the selective emitter from 58% to 47%. The output power density, sub-bandgap loss, and TPV conversion efficiency ( ${\eta _{\rm TPV}}$ ) for a GaSb cell illuminated by the selective thermal emitter at 1500 K were estimated. ${\eta _{\rm TPV}}$ drops from 13.7% to 11% due to the increased sub-bandgap emission at 1500 K. Essential approaches for mitigating the sub-bandgap losses to further improve ${\eta _{\rm TPV}}$ are also discussed.

42 ENGINEERING↗