SEARCH · Search NASA

Results for “Source term calculation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

NASA NTRS2023

Implementing a VLBI Time Delay Model for Earth-orbiting Satellites: Partial derivatives and Verification

This document describes the partial derivatives for commonly estimated parameters including antenna positions, satellite position, and satellite velocity from a near-field VLBI delay model for Earth-orbiting satellites. This model was presented in the Journal of Geodesy by Jaron & Nothnagel (2019). In this context, we present a streamlined version of the near-field VLBI delay. From this simplified model, we deduce a delay rate expression and calculate partial derivatives, maintaining only the terms with significant impact on the computed derivatives’ magnitude. To verify the accuracy of the simplified model and the partial derivatives computed from it, we have created a simple simulation in Matlab of an Earth-orbiting satellite at the altitude of a typical Global Navigation Satellite Systems (GNSS) satellite. From this simulation, we compare the simplified and original delays, and we verify the magnitude and direction of the partial derivatives against the numerically computed derivatives from the original delay model. The partial derivatives and simplified VLBI delay model detailed here are implemented in Fortran in the open-source library VTD.

Geodesy↗

DOE OSTI2024

Single-photon superradiance in individual caesium lead halide quantum dots

The brightness of an emitter is ultimately described by Fermi’s golden rule, with a radiative rate proportional to its oscillator strength times the local density of photonic states. As the oscillator strength is an intrinsic material property, the quest for ever brighter emission has relied on the local density of photonic states engineering, using dielectric or plasmonic resonators. By contrast, a much less explored avenue is to boost the oscillator strength, and hence the emission rate, using a collective behaviour termed superradiance. Recently, it was proposed that the latter can be realized using the giant oscillator-strength transitions of a weakly confined exciton in a quantum well when its coherent motion extends over many unit cells. Here we demonstrate single-photon superradiance in perovskite quantum dots with a sub-100 picosecond radiative decay time, almost as short as the reported exciton coherence time. The characteristic dependence of radiative rates on the size, composition and temperature of the quantum dot suggests the formation of giant transition dipoles, as confirmed by effective-mass calculations. The results aid in the development of ultrabright, coherent quantum light sources and attest that quantum effects, for example, single-photon emission, persist in nanoparticles ten times larger than the exciton Bohr radius.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2019

Spectral dynamics of shift current in ferroelectric semiconductor SbSI

Photoexcitation in solids brings about transitions of electrons/holes between different electronic bands. If the solid lacks an inversion symmetry, these electronic transitions support spontaneous photocurrent due to the geometric phase of the constituting electronic bands: the Berry connection. This photocurrent, termed shift current, is expected to emerge on the timescale of primary photoexcitation process. We observe ultrafast evolution of the shift current in a prototypical ferroelectric semiconductor antimony sulfur iodide (SbSI) by detecting emitted terahertz electromagnetic waves. By sweeping the excitation photon energy across the bandgap, ultrafast electron dynamics as a source of terahertz emission abruptly changes its nature, reflecting a contribution of Berry connection on interband optical transition. The shift excitation carries a net charge flow and is followed by a swing over of the electron cloud on a subpicosecond timescale. Understanding these substantive characters of the shift current with the help of first-principles calculation will pave the way for its application to ultrafast sensors and solar cells.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Automatically generated model for light alkene combustion

Light alkenes are common combustion intermediates for a variety of fuels. Therefore, understanding their oxidation and pyrolysis chemistry is key to building detailed mechanisms for heavier fuels. This work was focused on the development and evaluation of a detailed kinetic mechanism suitable for the combustion of light alkenes up to C 4 without the use of tuned parameters, instead, the parameter values come from first principles or direct measurements. The generated mechanism accurately estimates the laminar burning velocity (S u ) and ignition delay time (IDT) of light alkenes available in the literature, which represent fundamental combustion properties at a wide range of conditions. Because each parameter is thought to have a physically realistic value, not tuned to these measurements, the new model could be used as a sub-mechanism in models for other applications. The reaction network was generated with the open-source Reaction Mechanism Generator (RMG) software. Sensitivity analyses were performed under wide ranges of temperatures and pressures, allowing for the identification of the most impactful species and reactions. Based on these, a comprehensive thermochemistry database, including calculations on 550 molecules performed in this work at the CBS-QB3 level of theory, and a kinetic library, including theoretically-derived reaction rates retrieved from the literature, were built and used in the mechanism generation. The developed mechanism was compared against several existing detailed kinetic mechanisms for ethene, propene, 1-butene, 2-butene, and isobutene. Here, the newly generated model is the most accurate among the ones analyzed, in terms of fractional bias and normalized mean square error. Hence, this new model was used to analyze the chemistry of alkene combustion. Key rate coefficients were compared, to identify the cause of deviations between the models and possible areas for further improvements.

09 BIOMASS FUELS↗

NASA NTRS1986

Calculations of cosmogenic nuclide production rates in the Earth's atmosphere and their inventories

The production rates of cosmogenic isotopes in the Earth's atmosphere and their resulting terrestrial abundances have been calculated, taking into account both geomagnetic and solar-modulatory effects. The local interstellar flux was assumed to be that of Garcia-Munoz, et al. Solar modulation was accounted for using the heliocentric potential model and expressed in terms of the Deep River neutron monitor count rates. The geomagnetic field was presented by vertical cutoffs calculated by Shea and Smart and the non-vertical cutoffs calculated using ANGRI. The local interstellar particle flux was first modulated using the heliocentric potential field. The modulated cosmic-ray fluxes reaching the earth's orbit then interacted with the geomagnetic field as though it were a high-pass filter. The interaction of the cosmic radiation with the Earth's atmosphere was calculated utilizing the Bolztmann transport equation. Spallation cross sections for isotope production were calculated using the formalism of Silberberg and Tsao and other cross sections were taken from standard sources. Inventories were calculated by accounting from the variation in solar modulation and geomagnetic field strength with time. Results for many isotope, including C-14, Be-7 and Be-10 are in generally good agreement with existing data. The C-14 inventory, for instance, amounts to 1.75/sq cm(e)/s, in excellent agreement with direct estimates.

Obrien, K.↗

NASA NTRS1998

NO(x) Concentrations in the Upper Troposphere as a Result of Lightning

Upper tropospheric NO(x) controls, in part, the distribution of ozone in this greenhouse-sensitive region of the atmosphere. Many factors control NO(x) in this region. As a result it is difficult to assess uncertainties in anthropogenic perturbations to NO from aircraft, for example, without understanding the role of the other major NO(x) sources in the upper troposphere. These include in situ sources (lightning, aircraft), convection from the surface (biomass burning, fossil fuels, soils), stratospheric intrusions, and photochemical recycling from HNO3. This work examines the separate contribution to upper tropospheric "primary" NO(x) from each source category and uses two different chemical transport models (CTMS) to represent a range of possible atmospheric transport. Because aircraft emissions are tied to particular pressure altitudes, it is important to understand whether those emissions are placed in the model stratosphere or troposphere and to assess whether the models can adequately differentiate stratospheric air from tropospheric air. We examine these issues by defining a point-by-point "tracer tropopause" in order to differentiate stratosphere from troposphere in terms of NO(x) perturbations. Both models predict similar zonal average peak enhancements of primary NO(x) due to aircraft (approx. = 10-20 parts per trillion by volume (pptv) in both January and July); however, the placement of this peak is primarily in a region of large stratospheric influence in one model and centered near the level evaluated as the tracer tropopause in the second. Below the tracer tropopause, both models show negligible NO(x) derived directly from the stratospheric source. Also, they predict a typically low background of 1 - 20 pptv NO(x) when tropospheric HNO3 is constrained to be 100 pptv of HNO3. The two models calculate large differences in the total background NO(x) (defined as the source of NO(x) from lightning + stratosphere + surface + HNO3) when using identical loss frequencies for NO(x). This difference is primarily due to differing treatments of vertical transport. An improved diagnosis of this transport that is relevant to NO(x) requires either measurements of a surface-based tracer with a substantially shorter lifetime than Rn-222 or diagnosis and mapping of tracer correlations with different source signatures. Because of differences in transport by the two models we cannot constrain the source of NO(x) from lightning through comparison of average model concentrations with observations of NO(x).

Penner, Joyce E.↗

DOE OSTI2021

The effect of price-based demand response on carbon emissions in European electricity markets: The importance of adequate carbon prices

Price-based demand response (PBDR) has recently been attributed great economic but also environmental potential. However, the determination of its short-term effects on carbon emissions requires the knowledge of marginal emission factors (MEFs), which compared to grid mix emission factors (XEFs), are cumbersome to calculate due to the complex characteristics of national electricity markets. This study, therefore, proposes two merit order-based methods to approximate hourly MEFs and applies them to readily available datasets from 20 European countries for the years 2017–2019. Based on the calculated electricity prices, standardized daily load shifts were simulated which indicated that carbon emissions increased for 8 of the 20 countries and by 2.1% on average. Thus, under specific circumstances, PBDR leads to carbon emissions increases, mainly due to the economic advantage fuel sources such as lignite and coal have in the merit order. MEF-based load shifts reduced the mean resulting carbon emissions by 35%, albeit with 56% lower monetary cost savings compared to price-based load shifts. Finally, by repeating the load shift simulations for different carbon price levels, the impact of the carbon price on the resulting carbon emissions was analyzed. The Spearman correlation coefficient between carbon intensity and marginal cost along the German merit order substantially increased with increasing carbon price. The coefficients were -0.13 for the 2019 carbon price of 24.9 €/t, 0 for 42.6 €/t, and 0.4 for 100.0 €/t. Therefore, with adequate carbon prices, PBDR can be an effective tool for both economical and environmental improvement.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

NASA NTRS1990

The diffuse galactic gamma ray emission

The EGRET (Energetic Gamma-Ray Experiment Telescope) detector will provide a much more detailed view of the diffuse galactic gamma ray intensity in terms of higher resolution, greater statistical significance, and broader energy range than earlier missions. These observations will furnish insight into a number of very important questions related to the dynamics and structure of the Galaxy. A diffuse emission model is being developed that incorporates the latest information on matter distribution and source functions. In addition, it is tailored to the EGRET instrument response functions. The analysis code of the model maintains flexibility to accommodate the quality of the data that is anticipated. The discussion here focuses on the issues of the distributions of matter, cosmic rays, and radiation fields, and on the important source functions that enter into the model calculation of diffuse emission.

Bertsch, David L.↗

NASA NTRS1994

SeaWiFS technical report series. Volume 23: SeaWiFS prelaunch radiometric calibration and spectral characterization

Based on the operating characteristics of the Sea-viewing Wide Field-of-view Sensor (SeaWiFS), calibration equations have been developed that allow conversion of the counts from the radiometer into Earth-existing radiances. These radiances are the geophysical properties the instrument has been designed to measure. SeaWiFS uses bilinear gains to allow high sensitivity measurements of ocean-leaving radiances and low sensitivity measurements of radiances from clouds, which are much brighter than the ocean. The calculation of these bilinear gains is central to the calibration equations. Several other factors within these equations are also included. Among these are the spectral responses of the eight SeaWiFS bands. A band's spectral response includes the ability of the band to isolate a portion of the electromagnetic spectrum and the amount of light that lies outside of that region. The latter is termed out-of-band response. In the calibration procedure, some of the counts from the instrument are produced by radiance in the out-of-band region. The number of those counts for each band is a function of the spectral shape of the source. For the SeaWiFS calibration equations, the out-of-band responses are converted from those for the laboratory source into those for a source with the spectral shape of solar flux. The solar flux, unlike the laboratory calibration, approximates the spectral shape of the Earth-existing radiance from the oceans. This conversion modifies the results from the laboratory radiometric calibration by 1-4 percent, depending on the band. These and other factors in the SeaWiFS calibration equations are presented here, both for users of the SeaWiFS data set and for researchers making ground-based radiance measurements in support of Sea WiFS.

Barnes, Robert A.↗

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↗

NASA NTRS1981

Case studies of clear air turbulence using the diagnostic Richardson Number Tendency formulation

The results of four case studies of clear air turbulence (CAT) using the diagnostic Richardson number tendency (DRT) formulation are highlighted. The performance of this technique in resolving regions of documented CAT encounters is encouraging. Its operational adaptability appears particularly attractive in that input data can be supplied by the currently operational rawinsonde system. Two CAT indices are calculated deterministically, sensing synoptic scale changes in static stability and vertical wind shear conductive for supporting mesoscale CAT layers. These two indices reveal volumes of the troposphere which act as source regions for patches of CAT. The first, the time to reach the critical Richardson number necessary for the initiation of turbulent conditions, was first used by Oard (1974). A second index is devised which is an adaptation of Roach's (1970) work relating synoptic scale and mesoscale energetical coupling and gives more information on the relative intensity of these source regions. The output from the DRT computer module highlights specific regions of the atmosphere which can be interpreted operationally in terms of a CAT encounter probability.

Keller, J. L.↗

DOE OSTI2020

Tight-Binding-Based Method for Simulations of Complex Materials with Accompanying Database

The NRL-Tight-Binding (NRL-TB) method has been applied successfully to many monatomic, binary, and a few ternary systems. The method consists of a fit to Density Functional Theory (DFT) results and has transferability beyond the fitted DFT database including performing tight-binding molecular dynamics simulations. We are working on a new website with a detailed guide of how to use the method in performing the fit to DFT results (Fitting Code) and obtaining non-orthogonal TB parameters expressed in a polynomial form as environment-dependent Slater-Koster parameters (on-site terms), and neighbor distance-dependent Hamiltonian and Overlap terms. These parameters are then used in a Static Code to verify or make predictions for large systems often beyond the computational capabilities of standard DFT codes. Our starting point was a website developed at NRL some years ago which NRL no longer supports. We have now been augmenting that website with new examples of expanding the original capabilities of the method such as incorporating in the codes f-electron states. Our website (esd.cos.gmu.edu/tb) contains two separate branches: One named NRL-TB with the full documentation of the method including source codes and examples. The other branch, named ESD-dB, is a library of electronic structure calculations for all the monatomic elements (a subset of the elements also found in a recent book by the author of this proposal. See reference 22) and an approximately 700 hundred compounds in specific crystal structures. We discuss below the progress made so far in improving and expanding the capabilities of our TB method and the plans for further extensions of our approach.

36 MATERIALS SCIENCE↗

DOE OSTI2023

BUQEYE guide to projection-based emulators in nuclear physics

The BUQEYE collaboration (Bayesian Uncertainty Quantification: Errors in Your effective field theory) presents a pedagogical introduction to projection-based, reduced-order emulators for applications in low-energy nuclear physics. The term emulator refers here to a fast surrogate model capable of reliably approximating high-fidelity models. As the general tools employed by these emulators are not yet well-known in the nuclear physics community, we discuss variational and Galerkin projection methods, emphasize the benefits of offline-online decompositions, and explore how these concepts lead to emulators for bound and scattering systems that enable fast and accurate calculations using many different model parameter sets. We also point to future extensions and applications of these emulators for nuclear physics, guided by the mature field of model (order) reduction. All examples discussed here and more are available as interactive, open-source Python code so that practitioners can readily adapt projection-based emulators for their own work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1981

Observations and models of inertial waves in the deep ocean

A study of the structure of the inertial peak in deep ocean kinetic energy is presented, based on records taken from Polymode arrays deployed in the western North Atlantic Ocean. Results are interpreted in terms of both local sources and turning point effects on internal waves generated at lower latitudes, and it is found that three classes of environment and their corresponding spectra emerge from peak height variations: (1) the 1500-m level near the Mid-Atlantic Ridge, with the greatest peak height of 18 dB; (2) the upper and deep ocean over rough topography and the deep ocean underneath the Gulf Stream, with the intermediate peak height of 11.5 dB; and (3) the deep ocean over smooth topography, with the lowest peak height of 7.5 dB. Using the globally valid wave functions obtained by Munk and Phillips (1968), frequency spectra near f are calculated numerically. The model is latitudinally dependent, with the frequency shift and bandwidth of the inertial peak decreasing with latitude.

Fu, L.-L.↗

NASA NTRS1985

Infrared polar brightening on Jupiter. III - Spectrometry from the Voyager 1 IRIS experiment

Spectra from the Voyager 1 IRIS experiment confirm the existence of enhanced infrared emission near Jupiter's north magnetic pole in March 1979. The spectral characteristics of the enhanced emission are consistent with a Planck source function. A temperature-pressure profile is derived for the region near the north magnetic pole, from which quantitative abundance estimates of minor species are made. Some species previously detected on Jupiter, including CH3D, C2H2, and C2H6, have been observed again near the pole. Newly discovered species, not previously observed on Jupiter, include C2H4, C3H4, and C6H6. All of these species except CH3D appear to have enhanced abundances at the north polar region with respect to midlatitudes. Upper limits are determined for C4H2 and C3H8. The quantitative results are compared with model calculations based on ultraviolet results from the IUE satellite. The plausibility of the C6H6 identification is discussed in terms of the literature on C2H2 polymerization. The relation of C6H6 to cuprene is also discussed.

Kim, S. J.↗

NASA NTRS1993

Gas distribution and starbursts in shell galaxies

Detailed maps of most elliptical galaxies reveal that, whereas the greatest part of their luminous mass originates from a smooth distribution with a surface brightness approximated by a de Vaucouleurs law, a small percentage of their light is contributed by low surface brightness distortions termed 'fine structures'. The sharp-edged features called 'shells' are successfully reproduced by merger and infall models involving accretion from less massive companions. In this context, dwarf spheroidal and compact disk galaxies are likely progenitors of these stellar phenomena. However, it is probable that the sources of shell-forming material also contain significant amounts of gas. This component may play an important role in constraining the formation and evolution of shell galaxies. To investigate the effects of the gaseous component, numerical simulations were performed to study the tidal disruption of dwarf galaxies containing both gas and stars by more massive primaries, and the evolution of the ensuing debris. The calculations were performed with a hybrid N-body/hydrodynamics code. Collisionless matter is evolved using a conventional N-body technique and gas is treated using smoothed particle hydrodynamics in which self-gravitating fluid elements are represented as particles evolving according to Lagrangian hydrodynamic equations. An isothermal equation of state is employed so the gas remains at a temperature 104 K. Owing to the large mass ratio between the primary and companion, the primary is modeled as a rigid potential and the self-gravity of both galaxies is neglected.

Weil, Melinda L.↗

DOE OSTI2021

Saturn Radiation Dose Environment Characterization

To understand the environment where a time-resolved hard x-ray spectrometer (AXIOM) might be fielded, experiments and simulations were performed to analyze the radiation dose environment underneath the Saturn vacuum dome. Knowledge of this environment is critical to the design and placement of the spectrometer. Experiments demonstrated that the machine performance, at least in terms of on-axis dose, has not significantly changed over the decades. Simulations of the off-axis dose were performed to identify possible spectrometer locations of interest. The effects from the source and dome hardware as well as source distributions and angles of incidence on the radiation environment were also investigated. Finally, a unified radiation transport model was developed for two widely used radiation transport codes to investigate the off-axis dose profiles and the time-dependent x-ray energy spectrum. The demonstrated equivalence of the unified radiation transport model between the radiation transport codes allows the team to tie future time-dependent x-ray environment calculations to previous integral simulations for the Saturn facility.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Recent Improvements to the LAURA and HARA Codes

This paper describes recent improvements to the LAURA and HARA codes. LAURA is a CFD code for aerothermodynamics, and HARA evaluates the shock-layer radiation that provides the radiative source term for the flowfield energy equations and radiative heating to a surface. The next release of LAURA and HARA includes a variety of new capabilities. These new capabilities include an automated uncertainty quantification workflow for radiative heat transfer, options for specifying surface roughness and turbulent transition location in the algebraic turbulence models, and improved grid and solution interpolation techniques. Additionally, the computational efficiency of both LAURA and HARA have been improved. Optimization of the MPI communication routines in LAURA are shown to improve the parallel efficiency of the primary flow when running with multiple processes per block, and recent optimization of HARA leverage graphics processing unit (GPU) acceleration in the radiation calculations. Using GPU acceleration of HARA is shown to decrease the cost of the radiation line-of-sight calculation by approximately one order of magnitude for a 10.5 km/s Earth entry simulation.

LAURA HARA CFD 5.6↗