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.

SourceResearch connectionFrom yearThrough year

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

At least 307 records · Page 17

All sources NASA NTRS DOE OSTI

NASA NTRS1980

Simultaneous wave and particle observations in the auroral kilometric radiation source region

Simultaneous observations from the topside sounder and the soft particle spectrometer onboard the ISIS 1 satellite reveal that very specific conditions on the local electron density and the energetic electron distributions must both occur in the auroral kilometric radiation (AKR) source region. Such regions are associated with inverted V electron precipitation and with depletions in the local electron density. The electron velocity distribution functions obtained near the inverted V peak were found to increase near several keV energy with increasing velocity as required for plasma instability. The electron density observed near the inverted V peak was too high to support AKR for three events investigated, however, and the AKR source was identified with the edge of the inverted V where the density was low (less than or equal to 30/cu cm) in each case. Whereas this density depletion can extend deep into the ionosphere (approximately 1500 km altitude), the severe depletion associated with the AKR density cavity is restricted to higher altitudes (greater than 2750 km for an event studied in detail).

Benson, R. F.↗

NASA NTRS1989

Direct simulation of compressible turbulence

Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Zang, T. A.↗

NASA NTRS2010

The Herschel ATLAS: Evolution of the 250 Micrometer Luminosity Function Out to z = 0.5

We have determined the luminosity function of 250 micrometer-selected galaxies detected in the approximately equal to 14 deg(sup 2) science demonstration region of the Herschel-ATLAS project out to a redshift of z = 0.5. Our findings very clearly show that the luminosity function evolves steadily out to this redshift. By selecting a sub-group of sources within a fixed luminosity interval where incompleteness effects are minimal, we have measured a smooth increase in the comoving 250 micrometer luminosity density out to z = 0.2 where it is 3.6(sup +1.4) (sub -0.9) times higher than the local value.

Dye, S.↗

DOE OSTI2025

A tale of two polarization paradoxes: The diamagnetic polarization paradox

An accurate calculation of the total polarization charge density in a plasma is essential for a self-consistent determination of the electric field. Yet, the “diamagnetic polarization paradox” refers to the fact that there is a paradoxical factor of 1/2 difference between the pressure-driven “diamagnetic polarization” density calculated using real space drift theory vs action-angle space guiding center and gyrokinetic theory that has not been explained before. Here, we show that both results can be made consistent with one another. Half of the diamagnetic polarization is due to the transformation from the guiding center density to the real space density. The other half is due to the fact that, within the drift kinetic ordering assumptions, the guiding center density should be expressed as the gyroaverage of the density in the limit of vanishing Larmor radius. A comprehensive review of polarization is presented to complete the derivation, and then we derive results that are required for the polarization calculations to agree. Expressions for the diamagnetic polarization density are given that are accurate to first order in amplitude and all orders in gyroradius within the gyrokinetic theory for a constant magnetic field. Applications to Maxwell–Boltzmann particle distribution functions (PDFs), including anisotropic temperature, are presented. Local invariants, like total energy and toroidal momentum, do not generate net polarization effects; the electric and thermodynamic polarizations must precisely cancel. In contrast, anisotropic dependence on the magnetic moment generates a net polarization proportional to the temperature anisotropy. Finally, when sources are present, the equilibrium PDF is approximately the ratio of two orbit averages.

adiabatic theory↗

NASA NTRS2005

The Far-Infrared Spectrum of Arp 220

ISO/LWS grating observations of the ultraluminous infrared galaxy Arp 220 shows absorption in molecular lines of OH, H(sub 2)O, CH, NH, and NH(sub 3), as well as in the [O I] 63 micron line and emission in the [C II] 158 micron line. We have modeled the continuum and the emission/absorption of all observed features by means of a non-local radiative transfer code. The continuum from 25 to 1300 microns is modeled as a warm (106 K) nuclear region that is optically thick in the far-infrared, attenuated by an extended region (size 2") that is heated mainly through absorption of nuclear infrared radiation. The molecular absorption in the nuclear region is characterized by high excitation due to the high infrared radiation density. The OH column densities are high toward the nucleus (2 - 6 x 10(exp 17) cm(exp -2)) and the extended region (approximately 2 x 10(exp 17) cm(exp -2)). The H(sub 2)O column density is also high toward the nucleus (2 - 10 x 10(exp 17) cm(exp -2)) and lower in the extended region. The column densities in a halo that accounts for the absorption by the lowest lying levels are similar to what are found in the diffuse clouds toward the star forming regions in the Sgr B2 molecular cloud complex near the Galactic Center. Most notable are the high column densities found for NH and NH(sub 3) toward the nucleus, with values of approximately 1.5 x 10(exp 16) cm(exp -2) and approximately 3 x 10(exp 16) cm(exp -2), respectively, whereas the NH(sub 2) column density is lower than approximately 2 x 10(exp 15) cm(exp -2). A combination of PDRs in the extended region and hot cores with enhanced H(sub 2)O photodissociation and a possible shock contribution in the nuclei may explain the relative column densities of OH and H(sub 2)O, whereas the nitrogen chemistry may be strongly affected by cosmic ray ionization. The [C II] 158 micron line is well reproduced by our models and its deficit relative to the CII/FIR ratio in normal and starburst galaxies is suggested to be mainly a consequence of the dominant non-PDR component of far-infrared radiation, although our models alone cannot rule out extinction effects in the nuclei.

Gonzalez-Alfonso, Eduardo↗

NASA NTRS1989

Probabilities for gravitational lensing by point masses in a locally inhomogeneous universe

Probability functions for gravitational lensing by point masses that incorporate Poisson statistics and flux conservation are formulated in the Dyer-Roeder construction. Optical depths to lensing for distant sources are calculated using both the method of Press and Gunn (1973) which counts lenses in an otherwise empty cone, and the method of Ehlers and Schneider (1986) which projects lensing cross sections onto the source sphere. These are then used as parameters of the probability density for lensing in the case of a critical (q0 = 1/2) Friedmann universe. A comparison of the probability functions indicates that the effects of angle-averaging can be well approximated by adjusting the average magnification along a random line of sight so as to conserve flux.

Isaacson, Jeffrey A.↗

NASA NTRS1985

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗

DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1988

The structure of shocks with thermal conduction and radiative cooling

A general analysis is presented of the structure of a steady state, plane-parallel shock wave in which both thermal conduction and radiative cooling are important. The fluid is assumed to have a perfect-gas equation of state, with radiative cooling a function only of its temperature and density. Conduction in both diffusive and saturated regimes is treated. For the case of a strong shock, with conductivity and cooling function varying as power laws in temperature, approximate analytic solutions describing the shock wave are derived. For a plasma of solar composition, conduction is found to have a significant effect on the shock temperature and overall thickness of the postshock layer only for shock velocities greater than about 30,000 km/s, corresponding to shock temperatures greater than about 10 to the 10th K, but it affects the local structure of parts of the shock wave at much lower velocities. The effects of conduction are greatly enhanced if the heavy-element abundance is increased.

Lacey, Cedric G.↗

NASA NTRS1993

On the synthesis of resonance lines in dynamical models of structured hot-star winds

We examine basic issues involved in synthesizing resonance-line profiles from 1-D, dynamical models of highly structured hot-star winds. Although these models exhibit extensive variations in density as well as velocity, the density scale length is still typically much greater than the Sobolev length. The line transfer is thus treated using a Sobolev approach, as generalized by Rybicki & Hummer (1978) to take proper account of the multiple Sobolev resonances arising from the nonmonotonic velocity field. The resulting reduced-lambda-matrix equation describing nonlocal coupling of the source function is solved by iteration, and line profiles are then derived from formal solution integration using this source function. Two more approximate methods that instead use either a stationary or a structured, local source function yield qualitatively similar line-profiles, but are found to violate photon conservation by 10% or more. The full results suggest that such models may indeed be able to reproduce naturally some of the qualitative properties long noted in observed UV line profiles, such as discrete absorption components in unsaturated lines, or the blue-edge variability in saturated lines. However, these particular models do not yet produce the black absorption troughs commonly observed in saturated lines, and it seems that this and other important discrepancies (e.g., in acceleration time scale of absorption components) may require development of more complete models that include rotation and other 2-D and/or 3-D effects.

Puls, J.↗

NASA NTRS2003

Toroidal wave frequency at L=6-10: AMPTE/CCE Observations and Comparison with Theoretical Model

Magnetospheric standing Alfven waves are guided along the ambient magnetic field and their frequency depends on the mass density of the plasma distributed along the field lines. These properties allow us to use Alfen waves to map time-dependent phenomena between space and ground and to estimate the mass density. In this paper we present a statistical study of the spatial variation of the fundamental frequency f(sub T1) of toroidal-mode standing Alfen waves in the L range from 6 to 10, where L indicates the maximum geocentric distance on the field line. The data used for this analysis are energetic particle flux anisotropy (proxy of transverse electric field) and magnetic field measurements from the Active Magnetospheric Particle Tracer Explorers Charge Composition Explorer (CCE) spacecraft. Using CCE data covering 4 years we obtained approximately 5000 20-min intervals containing a clear signature of toroidal waves. The median f(sub T1), is 6-10 mHz at L=7 and decreases to 4-8 mHz at L = . The frequency tends to be lower at noon than at midnight. The observed frequencies are compared with numerically derived frequencies using an empirical mass density model and a magnetic field model. We found a good agreement for 12-24 MLT but a large discrepancy near MLT= 3. This may indicate that the flux tubes at this local time are more heavily loaded than specified in the Gallagher et al. model.

Takahashi, Kazue↗

NASA NTRS2002

Heat Transfer in High-Temperature Fibrous Insulation

The combined radiation/conduction heat transfer in high-porosity, high-temperature fibrous insulations was investigated experimentally and numerically. The effective thermal conductivity of fibrous insulation samples was measured over the temperature range of 300-1300 K and environmental pressure range of 1.33 x 10(exp -5)-101.32 kPa. The fibrous insulation samples tested had nominal densities of 24, 48, and 72 kilograms per cubic meter and thicknesses of 13.3, 26.6 and 39.9 millimeters. Seven samples were tested such that the applied heat flux vector was aligned with local gravity vector to eliminate natural convection as a mode of heat transfer. Two samples were tested with reverse orientation to investigate natural convection effects. It was determined that for the fibrous insulation densities and thicknesses investigated no heat transfer takes place through natural convection. A finite volume numerical model was developed to solve the governing combined radiation and conduction heat transfer equations. Various methods of modeling the gas/solid conduction interaction in fibrous insulations were investigated. The radiation heat transfer was modeled using the modified two-flux approximation assuming anisotropic scattering and gray medium. A genetic-algorithm based parameter estimation technique was utilized with this model to determine the relevant radiative properties of the fibrous insulation over the temperature range of 300-1300 K. The parameter estimation was performed by least square minimization of the difference between measured and predicted values of effective thermal conductivity at a density of 24 kilograms per cubic meters and at nominal pressures of 1.33 x 10(exp -4) and 99.98 kPa. The numerical model was validated by comparison with steady-state effective thermal conductivity measurements at other densities and pressures. The numerical model was also validated by comparison with a transient thermal test simulating reentry aerodynamic heating conditions.

Daryabeigi, Kamran↗

NASA NTRS1994

Absorption in X-ray spectra of high-redshift quasars

We present evidence that X-ray absorption is common in high-redshift quasars. We have studied six high-redshift (z approximately 3) quasars with the ROSAT Position Sensitive Proportional Counter (PSPC) of which four are in directions of low Galactic N(sub H). Three out of these four show excess absorption, while only three in approximately 50 z approximately less than 0.4 quasars do, indicating that such absorption must be common, but not ubiquitous, at high redshifts, and that the absorbers must lie at z greater than 0.4. The six quasars were: S5 0014+81, Q0420-388, PKS 0438-436, S4 0636+680. PKS 2000-330, PKS 2126-158, which have redshifts between 2.85 and 3.78. PKS 0438-436 and PKS 2126-158 show evidence for absorption above the local Galactic value at better than 99.999% confidence level. If the absorber is at the redshift of the quasar, then values of N(sub H) = (0.86(+0.49, -0.28)) x 10(exp 22) atoms/sq cm for PKS 0438-436, and N(sub H) = (1.45(+1.20, -0.64)) x 10(exp 22) atoms/ sq cm for PKS 2126-158, are implied, assuming solar abundances. The spectrum of S4 0636+680 also suggests the presence of a similarly large absorption column density at the 98% confidence level. This absorption reverses the trend for the most luminous active galactic nuclei (AGN) to have the least X-ray absorption, so a new mechanism is likely to be responsible. Intervening absorption due to damped Lyman(alpha) systems is a plausible cause. We also suggest, as an intrinsic model, that intracluster material, e.g., a cooling flow, around the quasar could account for both the X-ray spectrum and other properties of these quasars. All the quasars are radio-loud and three are gigahertz peaked (two of the three showing absorption). No excess absorption above the Galactic value is seen toward Q0420-388. This quasar has two damped Lyman(alpha) systems at z = 3.08. The limit on the X-ray column density implies a low ionization fraction, N(H I)/N(H) approximately greater than 4 x 10(exp -3) (3 (sigma)), for solar abundances, for these systems, and can set a weak limit on the size of the absorber. In the emitted frame these PSPC spectra cover the band approximately 0.5-10 keV, which has been well observed for low-redshift quasars and AGN. Comparison of high and low-redshift spectra in this emitted band shows no change of mean spectral index greater than Delta alpha(sub E) greater than 0.3 (99% confidence) with either redshift or luminosity, for radio-loud quasars.

Elvis, Martin↗

DOE OSTI2025

Predicting interface structure using the minima hopping method

Here, we adapt the minima hopping method (MHM) to the problem of interfacial structure prediction and apply it to study a canonical problem, the tilt grain boundaries in SrTiO 3 . Our method employs a hybrid approach by first exploring the potential energy surface (PES) of different grain boundary samplings with an empirical force field, among which the fifteen candidates with lower energies are then refined using ab initio density functional theory (DFT) calculations. During the exploratory stage, we bias the search using a local order parameter to primarily sample various reconstructions in the vicinity of the interface, while preserving the crystallinity of the bulk regions. We further enhance the search by incorporating initial structures with rigid body displacements to account for translational variations between bulk phases, enabling the MHM to effectively generate both stoichiometric and nonstoichiometric SrTiO 3 Σ⁢3(111)[110] and Σ⁢3(112)[110] grain boundaries. From an algorithmic standpoint, MHM outperforms earlier studies based on genetic algorithms (GA) by identifying more stable interfacial structures of several SrTiO 3 grain boundaries. The performance of the present implementation of the MHM approach is primarily limited by exploring an approximate description of the PES with a rather simple Buckingham potential. This limitation leads to variations in performance when compared to approaches utilizing more advanced surrogate PES models, such as direct DFT-PES sampling or GA with the embedded atom method (EAM). Despite the present limitations, the MHM approach is able to yield interfacial structures with comparable or lower interfacial energies in specific cases, such as Σ⁢3(111)[110] Γ=1, ±0.5 and Σ⁢3(112)[110] Γ= ±1, −2, underscoring the robustness of the MHM approach even with a simple approximation of the DFT PES. The MHM interfacial structure prediction method thus offers an efficient approach to understanding the grain boundaries and heterointerfaces at the atomic scale, providing an important prerequisite for effective materials design.

density functional theory↗

NASA NTRS1978

Convective heat flux in a laser-heated thruster

An analysis is performed to estimate the convective heating to the wall in a laser-heated thruster on the basis of a solution of the laminar boundary-layer equations with variable transport properties. A local similiarity approximation is used, and it is assumed that the gas phase is in equilibrium. For the thruster described by Wu (1976), the temperature and pressure distributions along the nozzle are obtained from the core calculation. The similarity solutions and heat flux are obtained from the freestream conditions of the boundary layer, in order to determine if it is necessary to couple the boundary losses directly to the core calculation. In addition, the effects of mass injection on the convective heat transfer across the boundary layer with large density-viscosity product gradient are examined.

Wu, P. K. S.↗

NASA NTRS1997

Io's Sodium Corona and Spatially Extended Cloud: A Consistent Flux Speed Distribution

For Io neutral cloud calculations, an SO2 source strength of approximately 4x10(exp 27) molecules/sec was determined by successfully matching the SO2(+) density profile near the satellite deduced from magnetometer data acquired by the Galileo spacecraft during its close flyby on December 7, 1995. The incomplete collision source velocity distribution for SO2 is the same as recently determined for the trace species atomic sodium by Smyth and Combi (1997). Estimates for the total energy loss rate (i.e. power) of O and S atoms escaping Io were also determined and imply a significant pickup current and a significant reduction in the local planetary magnetic field near Io.

Smyth, William H.↗

NASA NTRS2023

Formation of the Lyman Continuum during Solar Flares

The Lyman continuum (LyC; <911.12 Å) forms at the top of the chromosphere in the quiet Sun, making LyC a powerful tool for probing the chromospheric plasma during solar flares. To understand the effects of nonthermal energy deposition in the chromosphere during flares, we analyzed LyC profiles from a grid of field-aligned radiative-hydrodynamic models generated using the RADYN code as part of the F-CHROMA project. The spectral response of LyC, the temporal evolution of the departure coefficient of hydrogen, b 1 , and the color temperature, T c , in response to a range of nonthermal electron distribution functions, were investigated. The LyC intensity was seen to increase by 4–5.5 orders of magnitude during solar flares, responding most strongly to the nonthermal electron flux of the beam. Generally, b 1 decreased from 10 2 –10 3 to closer to unity during solar flares, indicating a stronger coupling to local conditions, while T c increased from 8–9 to 10–16 kK. T c was found to be approximately equal to the electron temperature of the plasma when b1 was at a minimum. Both optically thick and optically thin components of LyC were found to be in agreement with the interpretation of recent observations. The optically thick layer forms deeper in the chromosphere during a flare compared to quiescent periods, whereas the optically thin layers form at higher altitudes due to chromospheric evaporation, in low-temperature, high-density regions propagating upward. We put these results in the context of current and future missions.

Solar flares↗

DOE OSTI2025

Simulating large one-dimensional neutral-atom quantum systems

While abstract models of quantum computation assume a closed system of two-level states, practical quantum devices inevitably couple to the environment in some way, creating sources of noise. Understanding the tolerance to noise of specific quantum algorithms run on specific devices is important for determining the feasibility of quantum computing in the current noisy intermediate-scale quantum era. Of particular interest is understanding the noise sensitivity of these devices as more qubits are added to the system. Classical simulations are a useful tool to understand the effects of this noise, but direct classical simulations of open quantum systems are burdened by an exponentially growing cost in the number of qubits and a large local Hilbert space dimension. For onedimensional, shallow circuits, using tensor networks can replace this exponential cost with a linear one and simulate far wider systems than what would normally be available. In this paper, we describe a tensor network simulation of a neutral atom quantum system under the presence of noise, while introducing a purity-preserving truncation technique that compromises between the simplicity of the matrix product state and the positivity of the matrix product density operator. We apply this simulation to a near-optimized iteration of the quantum approximate optimization algorithm on a transverse field Ising model in order to investigate the influence of large system sizes on the performance of the algorithm. We find that while circuits with a large number of qubits fail more often under noise that depletes the qubit population, their outputs on a successful measurement are just as robust under Rydberg atom dissipation or qubit dephasing as smaller systems. However, such circuits might not perform as well under coherent multiqubit errors such as Rydberg atom crosstalk. We also find that the optimized parameters are especially robust to noise, suggesting that a noisier quantum system can be used to find the optimal parameters before switching to a cleaner system for measurements of observables.

Allen, James↗