Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local density approximations”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Interstellar Extinction in the Vicinity of the Galactic Center

We present J (1.2 microns) H (1-6 microns) K' (2.11 microns) and Br(gamma) (2.166 microns) images from four large regions within the central 40 pc of the Galaxy. Localized variations in the extinction, as determined by observations of the stellar population, are examined using the median H-K' color as a function of position within each region. As the value of the derived extinction from the stars is dependent upon the intrinsic magnitude of the assumed stellar type, the J-H vs. H-K' diagrams are first used to investigate the distribution of stellar types in the four regions. We find that there is a distinct OB population, contrary to earlier assumptions, with the ratio of K and M giants and supergiants to OB stars approximately twice that of the solar neighborhood. Although the on the scale of approx. 1 min. fluctuations in the extinction are on the order of A(sub V) approx. greater than 2, throughout the entire region the extinction varies from A(sub V) approx. greater than 25 to A(sub V) approx. less than 40. We also examine whether there is any variation in the extinction and stellar number density relative to the unusual radio features in these regions and do not find a significant correlation. Spectral imaging in Br(gamma) 2.166 microns emission shows a strong morphological correspondence between the 6 cm radio images and the diffuse Br(gamma) emission. By comparing the theoretical Br(gamma) flux derived from the radio flux using recombination theory, with our measured Br(gamma) flux, we obtain a second, independent, estimate of the extinction. We compare the two data sets and find that the extinction as derived from the stars is consistently greater, sometimes by a factor of two, than the value of the extinction derived from the Br(gamma) images. The derived extinction in various regions, however, is insufficient for any of these regions to be located behind the molecular clouds as previously observed in the Galactic Center, consistent with the theory that the observed radio emission is produced on the foreground surface of these clouds.

Cotera, Angela S.↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2026

Model Parameter Development for Complex Materials: Species-Specific Diffusion Barriers in 316 Stainless Steel from Systematic DFT Calculations

Vacancy-mediated diffusion barriers in 316 stainless steel have been systematically calculated using density functional theory to provide essential parameters for mesoscale microstructure evolution models. A statistical sampling approach employing 210 nudged elastic band calculations across multiple special quasi-random structures captures the effects of local chemical environments in this concentrated alloy. The computational methodology addresses challenges specific to chemically disordered systems, including proper magnetic treatment throughout multi-step calculations and validation against experimental structural properties. The calculated activation barriers reveal clear species-dependent diffusion behavior with the hierarchy Ni >> Fe ˜ Cr >> Mo. Nickel exhibits the highest barriers (0.74–1.31 eV, mean 1.045 eV), confirming its role as the slowest-diffusing major component. Iron and chromium show similar moderate barriers averaging 0.587 eV and 0.522 eV, respectively. Remarkably, molybdenum demonstrates exceptionally low barriers (0.12–0.28 eV, mean 0.194 eV), suggesting much higher mobility than previously recognized and potentially significant implications for precipitation kinetics and microstructure evolution. The barrier ranges remain consistent across different 316 SS compositions, supporting parameter transferability for modeling applications. The overall mean barrier of 0.64 eV provides a practical approximation for phase field simulations, while species-specific values enable detailed treatments of diffusion-controlled processes. This systematic approach establishes a validated framework for generating diffusion parameters in other concentrated alloys where experimental data are limited, while providing the first systematic set of species-specific barriers for predictive modeling of 316 stainless steel microstructure evolution.

36 MATERIALS SCIENCE↗
DOE OSTI2024

The impact of urban configuration types on urban heat islands, air pollution, CO 2 emissions, and mortality in Europe: a data science approach

The world is becoming increasingly urbanized. As cities around the world continue to grow, it is important for urban planners and policymakers to understand how different urban configuration patterns affect the environment and human health. We aimed at identifying European urban configuration types, based on the Local Climate Zones categories and street design variables from Open Street Map, and evaluating their association with motorized traffic flows, Surface Urban Heat Island (SUHI) intensities, tropospheric nitrogen dioxide (NO 2 ), CO 2 per capita emissions and age-standardized mortality. We considered 946 European cities from 31 countries for the analysis defined in the 2018 Urban Audit database, of which 919 European cities were analysed. Data were collected at a 250 m × 250 m grid cell resolution. We divided all cities into five concentric rings based on the Burgess concentric urban planning model and calculated the mean values of all variables for each ring. First, to identify distinct urban configuration types, we applied the Uniform Manifold Approximation and Projection for Dimension Reduction method, followed by the k-means clustering algorithm. Next, statistical differences in exposures (including SUHI) and mortality between the resulting urban configuration types were evaluated using a Kruskal–Wallis test followed by a post-hoc Dunn's test. We identified four distinct urban configuration types characterising European cities: compact high density (n=246), open low-rise medium density (n=245), open low-rise low density (n=261), and green low density (n=167). Compact high density cities were a small size, had high population densities, and a low availability of natural areas. In contrast, green low-density cities were a large size, had low population densities, and a high availability of natural areas and cycleways. The open low-rise medium and low-density cities were a small to medium size with medium to low population densities and low to moderate availability of green areas. Motorised traffic flows and NO 2 exposure were significantly higher in compact high density and open low rise medium density cities when compared with green low density and open low-rise low density cities. Additionally, green low-density cities had a significantly lower SUHI effect compared with all other urban configuration types. Per person CO 2 emissions were significantly lower in compact high density cities compared with green low density cities. Lastly, green low density cities had significantly lower mortality rates when compared with all other urban configuration types. Our findings indicate that, although the compact city model is more sustainable, European compact cities still face challenges related to poor environmental quality and health. Our results have notable implications for urban and transport planning policies in Europe and contribute to the ongoing discussion on which city models can bring the greatest benefits for the environment, climate, and health.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Energy transfer between localized emitters in photonic cavities from first principles

Radiative and nonradiative resonant couplings between defects are ubiquitous phenomena in photonic devices used in classical and quantum information technology applications. In this work, we present a first-principles approach to enable quantitative predictions of the energy transfer between defects in photonic cavities, beyond the dipole-dipole approximation and including the many-body nature of the electronic states. As an example, we discuss the energy transfer from a dipolelike emitter to an 𝐹 center in MgO in a spherical cavity. We show that the cavity can be used to controllably enhance or suppress specific spin-flip and spin-conserving transitions. Specifically, we predict that an ∼10–100 enhancement in the resonant energy transfer rate can be gained in the case of the 𝐹 center in MgO at ∼10 nm distances from a dipolar source, using rather moderate cavity with quality factor 𝑄 ∼ 400. We also show that a similar suppression in the transfer rate can be achieved by off-tuning the cavity resonance relative to the emitter transition energy. The framework presented here is general and readily applicable to a wide range of devices where localized emitters are embedded in microspheres, core-shell nanoparticles, and dielectric Mie resonators. Hence, our approach paves the way to predict how to control energy transfer in quantum memories and in ultrahigh-density optical memories, and in a variety of quantum information platforms.

First-principles calculations↗
NASA NTRS2002

Glennan Microsystems Initiative

During the 2001-2002 award period, we performed research on Pt/Ti/bare 6H-SiC and bare 4H-SiC interfaces in order to identify their electronic properties as a function of surface preparation. The overall aim of this work is to optimize the electronic properties of metal contacts to SiC as well as the active SiC material itself as a function of surface preparation and subsequent processing. Initially, this work has involved identifying bare surface, subsurface, and metal induced gap states at the metal-SiC contact and correlating energies and densities of deep levels with Schottky barrier heights. We used low energy electron-excited nanoluminescence (LEEN) spectroscopy, X-ray photoemission spectroscopy (XPS), and Secondary Ion Mass Spectrometry (SIMS) in order to correlate electronic states and energy bands with chemical composition, bonding, and crystal structure. A major development has been the discovery of polytype transformations that occur in 4H-SiC under standard microelectronic process conditions used to fabricate SiC devices. Our results are consistent with the stacking fault generation, defect formation, and consequent degradation of SiC recently reported for state-of-the-art ABB commercial diodes under localized electrical stress. Our results highlight the importance of -optimizing process conditions and material properties - anneal times, temperatures and doping to control such structural changes within epitaxial SiC layers. Thus far, we have established threshold times and temperatures beyond which 4H-SiC exhibits 3C-SiC transformation bands for a subset of dopant concentrations and process conditions. On the basis of this temperature time behavior, we have been able to establish an activation energy of approximately 2.5 eV for polytype transformation and dislocation motion. Work continues to establish the fundamental mechanisms underlying the polytype changes and its dependence on material parameters.

Brillson, Leonard J.↗
DOE OSTI2026

Importance of finite-size corrections for accurate ab initio modeling of carrier capture at semiconductor defects: A case study of substitutional C N in GaN

In ab initio studies of carrier-capture processes in defective semiconductor materials, the single-effective-mode formalism and the static-coupling approximation have become the predominant theoretical approaches for determining carrier-capture coefficients. The single-mode formalism relies on accurate nonequilibrium defect energies obtained from density-functional theory (DFT), where required inputs are a series of configurationally displaced, defect-containing supercells obtained using an interpolative ansatz, and where the DFT outputs are corresponding total energies that have traditionally been postprocessed using a long-established ground-state formulation of finite-size corrections and defect-formation energies. This formulation remains commonly used even though the defects that form a configuration-coordinate (CC) diagram typically exist as structures that are displaced from the ground state. To remedy this inconsistency, Kumagai has recently proposed novel methods for implementing finite-size corrections specifically intended for DFT calculations of the defect energies used to construct CC diagrams and implement the single-mode formalism [Y. Kumagai, Phys. Rev. B 107, L220101 (2023)]. Kumagai's approach builds on the latest finite-size-correction methods introduced to describe vertical charge-state transitions for charge-localizing point defects in semiconductors and insulators [T. Gake et al., Phys. Rev. B 101, 020102 (2020); S. Falletta et al., Phys. Rev. B 102, 041115 (2020)]. The newly identified finite-size artifact treated in these studies is the polarization charge induced on a configurationally frozen defect and its subsequent interaction with a vertical transition in charge state. In this work, we evaluate Kumagai's proposed methodology by applying it in a high-precision DFT study of carrier capture by substitutional C N in GaN, a well-characterized and technologically relevant defect and material. We have rigorously calculated C N defect energies across various supercell sizes for each defect configuration and charge state on the hole-capture CC diagram of C N (𝑞=−1), enabling a direct comparison of the slopes of the defect energies versus inverse cell size with those predicted by Kumagai. The most consequential prediction of Kumagai's method is that these slopes distinctly vary as the square of the linear-interpolation parameter used to construct the nonequilibrium defect configurations. Our results quantitatively support this prediction. Moreover, with these new finite-size corrections and multiple-cell-size DFT calculations in place, we find that the classical energy barrier for hole capture by C N (𝑞=−1) in GaN decreases to 0.092–0.127 eV. This finding confirms the recent ≈ 0.1 eV prediction of Reshchikov based on the weak temperature dependence for hole capture observed in photoluminescence experiments [M. A. Reshchikov, J. Appl. Phys. 129, 121101 (2021)]. These results stand in stark contrast to previously calculated barriers of 0.486 and 0.73 eV, which also used the single-mode formalism but were obtained by instead using ground-state-based finite-size corrections. Our reduced classical barrier for capture increases the temperature-dependent hole-capture coefficient of a C N (𝑞=−1) defect by more than two to four orders of magnitude for temperatures of 100–600 K, compared to the previous 0.486 eV results. While other defects may not be as dramatically affected as here, we suggest that incorporating proper finite-size corrections for the vertical-transition-like states embedded within CC diagrams is an essential, yet previously unrecognized, component of accurate modeling of carrier-capture when using the single-effective-mode formalism.

dielectric properties↗
NASA NTRS2002

Image Dis-Integration for Improved Plasmasphere Visualization

Traditionally, study of the plasmasphere has involved terrestrial observation of local characteristics. Global modeling of the plasmasphere in such an observation regime made use of an ensemble of (sparse) local measurements. Recently, sensors aboard the IMAGE (Imager for Magnetopause-to-Aurora Global Exploration) satellite (in particular, the EUV (Extreme Ultra Violet) Imager) have created the potential for truly global study of the plasmasphere. IMAGE was launched in spring of 2000 in an orbit with apogee altitude 7.2 R E (Earth radii) and perigee altitude 1000 km. IMAGE's EUV sensor allows an external view of the distribution of cold plasma in the plasmasphere to be acquired. EUV is designed to image light emission at 30.4 nanometers, which is the emission wavelength of the He + ion in the presence of solar radiation. He + makes up approximately 15-20% of the plasma in the plasmasphere, thus imaging of He + enables determination of plasma distribution. The EUV instrument provides a 90° by 84° field of view which is imaged as an equally spaced 150x140 pixel array on a spherical imaging surface. The EUV produces an image approximately every 10 minutes when the sensor is operating. Since EUV images contain line-of-sight integrations of plasma distributions, they do not directly express equatorial plane density (which would enable comparison of observed plasma distributions with predictions from models). Furthermore, the plasma density at any point in three-space is not known. The goal of our work was development of a technique that can enable plasma density to be determined throughout three-space. Our approach to creation of a three-space representation of the plasma distribution involves disintegrating the EUV lines of sight to form a volumetric map of plasma densities.

Timothy S. Newman↗
NASA NTRS1994

The cosmological dependence of cluster density profiles

We use N-body simulations to study the shape of mean cluster density and velocity profiles in the nonlinear regime formed via gravitational instability. The dependence of the final structure on both cosmology and initial density field is examined, using a grid of cosmologies and scale-free initial power spectra P(k) varies as k(exp n). Einstein-de Sitter, open (Omega(sub 0) = 0.2 and 0.1) and flat, low density (Omega(sub 0) = 0.2 lambda(sub 0) = 0.8) models are examined, with initial spectral indices n = -2, -1 and 0. For each model, we stack clusters in an appropriately scaled manner to define an average density profile in the nonlinear regime. The profiles are well fit by a power law rho(r) varies as r(exp -alpha) for radii whereat the local density contrast is between 100 and 3000. This covers 99% of the cluster volume. We find a clear trend toward steeper slopes (larger alphas) with both increasing n and decreasing Omega(sub 0). The Omega(sub 0) dependence is partially masked by the n dependence; there is degeneracy in the values of alpha between the Einstein-de Sitter and flat, low-density cosmologies. However, the profile slopes in the open models are consistently higher than the Omega = 1 values for the range of n examined. Cluster density profiles are thus potentially useful cosmological diagnostics. We find no evidence for a constant density core in any of the models, although the density profiles do tend to flatten at small radii. Much of the flattening is due to the force softening required by the simulations. An attempt is made to recover the unsoftened profiles assuming angular momentum invariance. The recovered profiles in Einstein-de Sitter cosmologies are consistent with a pure power law up to the highest density contrasts (10(exp 6)) accessible with our resolution. The low-density models show significant deviation from a power law above density contrasts approximately 10(exp 5). We interpret this curvature as reflecting the non-scale-invariant nature of the background cosmology in these models. These results are at the limit of our resolution and so should be tested in the future using simulations with larger numbers of particles. Such simulations will also provide insight on the broader problem of understanding, in a statistical sense, the full phase space structure of collapsed, cosmological halos.

Crone, Mary M.↗
NASA NTRS2012

Joining and Assembly of Bulk Metallic Glass Composites Through Capacitive Discharge

Bulk metallic glasses (BMGs), a class of amorphous metals defined as having a thickness greater than 1 mm, are being broadly investigated by NASA for use in spacecraft hardware. Their unique properties, attained from their non-crystalline structure, motivate several game-changing aerospace applications. BMGs have low melting temperatures so they can be cheaply and repeatedly cast into complex net shapes, such as mirrors or electronic casings. They are extremely strong and wear-resistant, which motivates their use in gears and bearings. Amorphous metal coatings are hard, corrosion-resistant, and have high reflectivity. BMG composites, reinforced with soft second phases, can be fabricated into energy-absorbing cellular panels for orbital debris shielding. One limitation of BMG materials is their inability to be welded, bonded, brazed, or fastened in a convenient method to form larger structures. Cellular structures (which can be classified as trusses, foams, honeycombs, egg boxes, etc.) are useful for many NASA, commercial, and military aerospace applications, including low-density paneling and shields. Although conventional cellular structures exhibit high specific strength, their porous structures make them challenging to fabricate. In particular, metal cellular structures are extremely difficult to fabricate due to their high processing temperatures. Aluminum honeycomb sandwich panels, for example, are used widely as spacecraft shields due to their low density and ease of fabrication, but suffer from low strength. A desirable metal cellular structure is one with high strength, combined with low density and simple fabrication. The thermoplastic joining process described here allows for the fabrication of monolithic BMG truss-like structures that are 90% porous and have no heat-affected zone, weld, bond, or braze. This is accomplished by welding the nodes of stacked BMG composite panels using a localized capacitor discharge, forming a single monolithic structure. This removes many complicated and costly fabrication steps. Moreover, the cellular structures detailed in this work are among the highest- strength and most energy-absorbent materials known. This implies that a fabricated structure made from these materials would have unequaled mechanical properties compared to other metal foams or trusses. The process works by taking advantage of the electrical properties of the matrix material in the metal-matrix composite, which in this case is a metallic glass. Due to the random nanoscale arrangement of atoms (without any grain boundaries), the matrix glass exhibits a near-constant electrical resistivity as a function of temperature. By placing the composite panels between two copper electrode plates and discharging a capacitor, the entire matrix of the panel can be heated to approximately 700 C in 10 milliseconds, which is above the alloy s solidus but below the liquidus. By designing the geometry of the panels into the shape of an egg box, the electrical discharge localizes only in the tips of each pyramidal cell. By applying a forging load during discharge, the nodes of the panels can be fused together into a single piece, which then dissipates heat through radiation back into a glassy state. This means that two panels can be metallurgically fused into one panel with no heat-affected zone, creating a seamless connection between panels. During the process, the soft metal particles (dendrites) that are uniformly distributed in the glassy matrix to increase the toughness are completely unaffected by the thermoplastic joining. The novelty is that a truss (or foam-like) structure can be formed with excellent energy- absorbing capabilities without the need for machining. The technique allows for large-scale fabrication of panels, well-suited for spacecraft shields or military vehicle door panels. Crystalline metal cellular structures cannot be fabricated using the thermoplastic joining technique described here. If metal panels were te assembled into a cellular structure, they would either have to be welded, brazed, bonded, or fastened together, creating a weak spot in the structure at each connection. Welded parts require a welding material to be added to the joint and exhibit a soft and weak heat-affected zone. Brazing and bonding do not form a metallurgical joint and thus exhibit low strengths, especially when the panels are pulled apart and fasteners require high-stress-concentration holes to be drilled. No equivalent rapid heating method exists for assembling metal panels together into cellular structures, and thus, those parts must be foamed, machined, or investment cast if they are to form a monolithic structure. If the crystalline panels were to be joined using capacitive discharge, as with a spot welder, their bond would be very weak, and the panels would have to be extremely thin. In contrast, the strength of joined BMG parts has been demonstrated to have strength comparable to the parent material. This technique opens up the possibility of using large-scale BMG hardware in spacecraft, military, or commercial applications.

Hofmann, Douglas C.↗
DOE OSTI2025

Molecular fluctuations inhibit intermittency in compressible turbulence

In the standard picture of fully developed turbulence, highly intermittent hydrodynamic fields are nonlinearly coupled across scales, where local energy cascades from large scales into dissipative vortices and large density gradients. Microscopically, however, constituent fluid molecules are in constant thermal (Brownian) motion, but the role of molecular fluctuations in large-scale turbulence is largely unknown, and with rare exceptions, it has historically been considered irrelevant at scales larger than the molecular mean free path. Recent theoretical and computational investigations have shown that molecular fluctuations can impact energy cascade at Kolmogorov length scales. Here, we show that molecular fluctuations not only modify energy spectrum at wavelengths larger than the Kolmogorov length in compressible turbulence, but also significantly inhibit spatio-temporal intermittency across the entire dissipation range. Using large-scale direct numerical simulations of computational fluctuating hydrodynamics, we demonstrate that the extreme intermittency characteristic of turbulence models is replaced by nearly Gaussian statistics in the dissipation range. These results demonstrate that the compressible Navier–Stokes equations should be augmented with molecular fluctuations to accurately predict turbulence statistics across the dissipation range. Our findings have significant consequences for turbulence modelling in applications such as astrophysics, reactive flows and hypersonic aerodynamics, where dissipation-range turbulence is approximated by closure models.

compressible turbulence↗
NASA NTRS1994

The Local Supercluster as a test of cosmological models

The Local Supercluster kinematic properties (the Local Group infall toward the Virgo Cluster and galaxy density distribution about the Virgo Cluster) in various cosmological models are examined utilizing large-scale N-body (PM) simulations 500(exp 3) cells, 250(exp 3) particles, and box size of 400 h(exp -1) Mpc) and are compared to observations. Five models are investigated: (1) the standard, Cosmic Background Explorer Satellite (COBE)-normalized cold dark matter (CDM) model with omega = 1, h = 0.5, and sigma(sub 8) = 1.05; (2) the standard Hot Dark Matter (HDM) model with omega = 1, h = 0.75, and sigma(sub 8) = 1; (3) the tilted CDM model with omega = 1, h = 0.5, n = 0.7, and sigma(sub 8) = 0.5; (4) a CDM + lambda model with omega = 0.3, lambda = 0.7, h = 2/3, and sigma(sub 8) = 2/3; (5) the PBI model with omega = 0.2, h = 0.8, x = 0.1, m = -0.5, and sigma(sub 8) = 0.9. Comparison of the five models with the presently available observational measurements v(sub LG) = 85 - 305 km/s (with mean of 250 km/s), delta(n(sub g))/(n(sub g)-bar) = 1.40 + or - 0.35) suggests that an open universe with omega approximately 0.5 (with or without lambda) and sigma(sub 8) approximately 0.8 is preferred, with omega = 0.3-1.0 (with or without lambda) and sigma(sub 8) = 0.7-1.0 being the acceptable range. At variance with some previous claims based on either direct N-body or spherical nonlinear approaches, we find that a flat model with sigma(sub 8) approximately 0.7-1.0 seems to be reasonably consistent with observations. However, if one favors the low limit of v(sub LG) = 85 km/s, then an omega approximately 0.2-0.3 universe seems to provide a better fit, and flat (omega = 1) models are ruled out at approximately 95% confidence level. On the other hand, if the high limit of v(sub LG) = 350 km/s is closer to the truth, then it appears that omega approximately 0.7-0.8 is more consistent. This test is insensitive to the shape of the power spectrum, but rather sensitive to the normalization of the perturbation amplitude on the relevant scale (e.g., sigma(sub 8)) and omega. We find that neither linear nor nonlinear relations (with spherical symmetry) are good approximations for the relation between radial peculiar velocity and density perturbation, i.e., nonspherical effects and gravitational tidal field are important. The derived omega using either of the two relations is underestimated. In some cases, this error is as large as a factor of 2-4.

Cen, Renyue↗
NASA NTRS2007

On the Nature of the Variability Power Decay towards Soft Spectral States in X-Ray Binaries. Case Study in Cyg X-1

A characteristic feature of the Fourier Power Density Spectrum (PDS) observed from black hole X-ray binaries in low/hard and intermediate spectral states is a broad band-limited noise, characterized by a constant below some frequency (a "break" frequency) and a power law above this frequency. It has been shown that the variability of this type can be produced by the inward diffusion of the local driving perturbations in a bounded configuration (accretion disk or corona). In the framework of this model, the perturbation diffusion time to is related to the phenomenological break frequency, while the PDS power-law slope above the "break" is determined by the viscosity distribution over the configuration. The perturbation diffusion scenario explains the decay of the power of X-ray variability observed in a number of compact sources (containing black hole and neutron star) during an evolution of theses sources from low/hard to high/soft states. We compare the model predictions with the subset of data from Cyg X-1 collected by the Rossi X-ray Time Explorer (RXTE). Our extensive analysis of the Cyg X-1 PDSs demonstrates that the observed integrated power P(sub x), decreases approximately as a square root of the characteristic frequency of the driving oscillations v(sub dr). The RXTE observations of Cyg X-1 allow us to infer P(sub dr), and t(sub o) as a function of v(sub dr). We also apply the basic parameters of observed PDSs, power-law index and low frequency quasiperiodic oscillations. to infer Reynolds (Re) number from the observations using the method developed in our previous paper. Our analysis shows that Re-number increases from values about 10 in low/hard state to that about 70 during the high/soft state. Subject headings: accretion, accretion disks-black hole physics-stars:individual (Cyg X-1) :radiation mechanisms: nonthermal-physical data and processes

Titarchuk, Lev↗
NASA NTRS1994

The cosmic X-ray background-IRAS galaxy correlation and the local X-ray volume emissivity

We have cross-correlated the galaxies from the IRAS 2 Jy redshift survey sample and the 0.7 Jy projected sample with the all-sky cosmic X-ray background (CXB) map obtained from the High Energy Astronomy Observatory (HEAO) 1 A-2 experiment. We have detected a significant correlation signal between surface density of IRAS galaxies and the X-ray background intensity, with W(sub xg) = (mean value of ((delta I)(delta N)))/(mean value of I)(mean value of N)) of several times 10(exp -3). While this correlation signal has a significant implication for the contribution of the local universe to the hard (E greater than 2 keV) X-ray background, its interpretation is model-dependent. We have developed a formulation to model the cross-correlation between CXB surface brightness and galaxy counts. This includes the effects of source clustering and the X-ray-far-infrared luminosity correlation. Using an X-ray flux-limited sample of active galactic nuclei (AGNs), which has IRAS 60 micrometer measurements, we have estimated the contribution of the AGN component to the observed CXB-IRAS galaxy count correlations in order to see whether there is an excess component, i.e., contribution from low X-ray luminosity sources. We have applied both the analytical approach and Monte Carlo simulations for the estimations. Our estimate of the local X-ray volume emissivity in the 2-10 keV band is rho(sub x) approximately = (4.3 +/- 1.2) x 10(exp 38) h(sub 50) ergs/s/cu Mpc, consistent with the value expected from the luminosity function of AGNs alone. This sets a limit to the local volume emissivity from lower luminosity sources (e.g., star-forming galaxies, low-ionization nuclear emission-line regions (LINERs)) to rho(sub x) less than or approximately = 2 x 10(exp 38) h(sub 50) ergs/s/cu Mpc.

Miyaji, Takamitsu↗
NASA NTRS2020

Formation and Topology of Foreshock Bubbles

We use global and local hybrid (kinetic ions and fluid electrons) simulations to investigate the conditions under which foreshock bubbles (FBs) form and how their topology changes with solar wind conditions. FBs form as a result of the interaction between solar wind discontinuities and backstreaming ion beams in the foreshock. They consist of an outer shock and its associated sheath plasma and a low density high temperature core with low magnetic field strength. The structure of FBs is determined by the angle between the interplanetary magnetic field and the normal to the solar wind discontinuity. We show that interaction of rotational discontinuities with the foreshock during small angles between the interplanetary magnetic field and discontinuity normal results in the formation of a nearly spherical bubble with a radius that scales with the width of the foreshock. As this angle increases, FBs become more elongated and eventually become nearly planar structures with dimensions that scale with the length of the foreshock. Despite this transformation, the signatures of FBs in spacecraft time series data remain the same in agreement with the observations. Global simulation results show that FBs form when the solar wind flow speed corresponds to high or intermediate Alfvén Mach numbers (approximately >7 MA). In general, this is tied to the relative speed between the solar wind and ion beams and drop in density of the back streaming ions.

Nick Omidi↗
NASA NTRS1992

Lunar Observer Laser Altimeter observations for lunar base site selection

One of the critical datasets for optimal selection of future lunar landing sites is local- to regional-scale topography. Lunar base site selection will require such data for both engineering and scientific operations purposes. The Lunar Geoscience Orbiter or Lunar Observer is the ideal precursory science mission from which to obtain this required information. We suggest that a simple laser altimeter instrument could be employed to measure local-scale slopes, heights, and depths of lunar surface features important to lunar base planning and design. For this reason, we have designed and are currently constructing a breadboard of a Lunar Observer Laser Altimeter (LOLA) instrument capable of acquiring contiguous-footprint topographic profiles with both 30-m and 300-m along-track resolution. This instrument meets all the severe weight, power, size, and data rate limitations imposed by Observer-class spacecraft. In addition, LOLA would be capable of measuring the within-footprint vertical roughness of the lunar surface, and the 1.06-micron relative surface reflectivity at normal incidence. We have used airborne laser altimeter data for a few representative lunar analog landforms to simulate and analyze LOLA performance in a 100-km lunar orbit. We demonstrate that this system in its highest resolution mode (30-m diameter footprints) would quantify the topography of all but the very smallest lunar landforms. At its global mapping resolution (300-m diameter footprints), LOLA would establish the topographic context for lunar landing site selection by providing the basis for constructing a 1-2 km spatial resolution global, geodetic topographic grid that would contain a high density of observations (e.g., approximately 1000 observations per each 1 deg by 1 deg cell at the lunar equator). The high spatial and vertical resolution measurements made with a LOLA-class instrument on a precursory Lunar Observer would be highly synergistic with high-resolution imaging datasets, and will allow for direct quantification of critical slopes, heights, and depths of features visible in images of potential lunar base sites.

Garvin, James B.↗
DOE OSTI2025

Assessment of integral models for non-Boussinesq lazy plumes using numerical simulations

Integral modelling of turbulent buoyant plumes is crucial for rapid predictions of plume characteristics. While the governing equations are typically derived using self-similarity and a Boussinesq approximation, these assumptions may not hold for plumes originating from finite-area sources with large density ratios. Here, this work evaluates the accuracy of integral-scale models for non-Boussinesq lazy plumes using high-fidelity numerical simulations of turbulent helium plumes. We analyse the plume kinematics by computing vertical fluxes, plume radius and radial profiles, establishing some disparities between common practice and physical accuracy. We identify how the definition of the plume radius changes the perception of the plume structure when the flow is not self-similar and derive a relationship between the flux-based and threshold-based definitions without requiring self-similarity. We then examine the plume dynamics by evaluating the source terms from the governing plume equations. Our results support neglecting diffusive and viscous effects but emphasise the importance of the mean pressure gradient, even in the self-similar regime. Two coefficients need to be modelled: the well-known entrainment coefficient and the lesser-known momentum correction coefficient, which is a correction required for the momentum equation to account for self-similar and slender approximations. The momentum correction coefficient is found to be approximately constant and slightly greater than the assumed value of 1. The standard entrainment coefficient models perform well up to a local Richardson number three times the asymptotic value but overpredict entrainment for larger Richardson numbers. We propose a correction using the known finite limit of entrainment at infinite Richardson number.

Meehan, Michael Alexander [Sandia National Laborat↗
NASA NTRS2010

EPR Imaging at a Few Megahertz Using SQUID Detectors

An apparatus being developed for electron paramagnetic resonance (EPR) imaging operates in the resonance-frequency range of about 1 to 2 MHz well below the microwave frequencies used in conventional EPR. Until now, in order to obtain sufficient signal-to-noise radios (SNRs) in conventional EPR, it has been necessary to place both detectors and objects to be imaged inside resonant microwave cavities. EPR imaging has much in common with magnetic resonance imaging (MRI), which is described briefly in the immediately preceding article. In EPR imaging as in MRI, one applies a magnetic pulse to make magnetic moments (in this case, of electrons) precess in an applied magnetic field having a known gradient. The magnetic moments precess at a resonance frequency proportional to the strength of the local magnetic field. One detects the decaying resonance-frequency magnetic- field component associated with the precession. Position is encoded by use of the known relationship between the resonance frequency and the position dependence of the magnetic field. EPR imaging has recently been recognized as an important tool for non-invasive, in vivo imaging of free radicals and reduction/oxidization metabolism. However, for in vivo EPR imaging of humans and large animals, the conventional approach is not suitable because (1) it is difficult to design and construct resonant cavities large enough and having the required shapes; (2) motion, including respiration and heartbeat, can alter the resonance frequency; and (3) most microwave energy is absorbed in the first few centimeters of tissue depth, thereby potentially endangering the subject and making it impossible to obtain adequate signal strength for imaging at greater depth. To obtain greater penetration depth, prevent injury to the subject, and avoid the difficulties associated with resonant cavities, it is necessary to use lower resonance frequencies. An additional advantage of using lower resonance frequencies is that one can use weaker applied magnetic fields: For example, for a resonance frequency of 1.4 MHz, one needs a magnetic flux density of 0.5 Gauss approximately the flux density of the natural magnetic field of the Earth.

Hahn, Inseob↗