SEARCH · Search NASA

Results for “Density function theory”

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 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Ion mixing in the plasma sheet boundary layer by drift instabilities

The linear stability properties of collisionless drift instabilities are analyzed in a Harris equilibrium model of the plasma sheet boundary layer (PSBL). The strearmng ions with drift-type instabilities driven in the PSBL are considered. The fluid approximation leads to growth but predicts that the mode width approaches the gyroradius of the energetic ions. Thus an integral equation theory for the modes is developed taking into account that in the PSBL the curvature drift is weak compared with the grad-B drift. The exact wave particle resonance is kept in the nonlocal response functions. Plasma density, temperature, and magnetic gradient drift motions are taken into account. The drift modes produce an anomalous cross-field momentum transport mixing the PSBL ions on the time scale of tens of seconds. A nonlinear simulation is performed which shows the coalescence of the small scale, fast growing modes into large-scale vortices. The relation between these collective modes and plasma sheet transport phenomena is discussed including the comparison with the competing plasma mixing from single-particle stochasticity.

Horton, W.↗

NASA NTRS2005

An Ensemble Approach to Building Mercer Kernels with Prior Information

This paper presents a new methodology for automatic knowledge driven data mining based on the theory of Mercer Kernels, which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly dimensional feature space. we describe a new method called Mixture Density Mercer Kernels to learn kernel function directly from data, rather than using pre-defined kernels. These data adaptive kernels can encode prior knowledge in the kernel using a Bayesian formulation, thus allowing for physical information to be encoded in the model. Specifically, we demonstrate the use of the algorithm in situations with extremely small samples of data. We compare the results with existing algorithms on data from the Sloan Digital Sky Survey (SDSS) and demonstrate the method's superior performance against standard methods. The code for these experiments has been generated with the AUTOBAYES tool, which automatically generates efficient and documented C/C++ code from abstract statistical model specifications. The core of the system is a schema library which contains templates for learning and knowledge discovery algorithms like different versions of EM, or numeric optimization methods like conjugate gradient methods. The template instantiation is supported by symbolic-algebraic computations, which allows AUTOBAYES to find closed-form solutions and, where possible, to integrate them into the code.

Srivastava, Ashok N.↗

NASA NTRS1976

A quasi-Newton approach to optimization problems with probability density constraints

A quasi-Newton method is presented for minimizing a nonlinear function while constraining the variables to be nonnegative and sum to one. The nonnegativity constraints were eliminated by working with the squares of the variables and the resulting problem was solved using Tapia's general theory of quasi-Newton methods for constrained optimization. A user's guide for a computer program implementing this algorithm is provided.

Tapia, R. A.↗

NASA NTRS1986

Nonlocal screening in metal surfaces

Due to the effect of the nonuniform environment on the static screening of the Coulomb potential, the local-density approximation for the particle-hole interaction is found to be inadequate to determine the surface energy of simple metals. Use of the same set of single-particle states, and thus the same one-body density and the same work function, has eliminated the single-electron states in favor of the structure of the short-ranged correlations as the basis of this effect. A posteriori simplifications of the Fermi hypernetted-chain theory may be found to allow the same calculational accuracy with simpler computational tools.

Krotscheck, E.↗

NASA NTRS1989

A new adaptive control approach for aerospace vehicles with parameter uncertainties

A new stochastic adaptive control structure is developed for the problem of combined parameter estimation and control of aerospace vehicles with changing parameters. Parameter uncertainties are modeled as first-order Gauss-Markov processes, and are introduced to the system dynamics through a small parameter. It is assumed that an accurate inertial measurement unit gives perfect measurements of the state variables. Since the stochastic system is assumed to be Gauss-Markov, the density function of the parameters given these measurements is conditionally Gaussian. Based on this conditionally Gaussian density, the problem of minimizing a quadratic cost over an infinite time horizon can be set up within the framework of stochastic optimal control theory. The optimal feedback control law is derived from a straightforward expansion of the Hamilton-Jacobi-Bellman equation, based on the LQG solution. The resulting nonlinear controller is applied to the pitch axis control of a space platform with uncertain moments of inertia and is shown to produce marked improvement over a fixed controller.

Hahn, Yungsun↗

NASA NTRS2004

Putting Priors in Mixture Density Mercer Kernels

This paper presents a new methodology for automatic knowledge driven data mining based on the theory of Mercer Kernels, which are highly nonlinear symmetric positive definite mappings from the original image space to a very high, possibly infinite dimensional feature space. We describe a new method called Mixture Density Mercer Kernels to learn kernel function directly from data, rather than using predefined kernels. These data adaptive kernels can en- code prior knowledge in the kernel using a Bayesian formulation, thus allowing for physical information to be encoded in the model. We compare the results with existing algorithms on data from the Sloan Digital Sky Survey (SDSS). The code for these experiments has been generated with the AUTOBAYES tool, which automatically generates efficient and documented C/C++ code from abstract statistical model specifications. The core of the system is a schema library which contains template for learning and knowledge discovery algorithms like different versions of EM, or numeric optimization methods like conjugate gradient methods. The template instantiation is supported by symbolic- algebraic computations, which allows AUTOBAYES to find closed-form solutions and, where possible, to integrate them into the code. The results show that the Mixture Density Mercer-Kernel described here outperforms tree-based classification in distinguishing high-redshift galaxies from low- redshift galaxies by approximately 16% on test data, bagged trees by approximately 7%, and bagged trees built on a much larger sample of data by approximately 2%.

Srivastava, Ashok N.↗

NASA NTRS2018

Acoustic Radiation Pressure

The theoretical foundation of acoustic radiation pressure in plane wave beams is reexamined. It is shown from finite deformation theory and the Boltzmann-Ehrenfest Adiabatic Principle that the Brillouin stress tensor (BST) is the radiation stress in Lagrangian coordinates (not Eulerian coordinates) and that the terms in the BST are not the momentum flux density and mean excess Eulerian stress but are simply contributions to the variation in the wave oscillation period resulting from changes in path length and true wave velocity, respectively, from virtual variations in the strain. It is shown that the radiation stress in Eulerian coordinates is the mean Cauchy stress (not the momentum flux density, as commonly assumed) and that Langevin's second relation does not yield an assessment of the mean Eulerian pressure, since the enthalpy used in the traditional derivations is a function of the thermodynamic tensions - not the Eulerian pressure. It is shown that the transformation between Lagrangian and Eulerian quantities cannot be obtained from the commonly-used expansion of one of the quantities in terms of the particle displacement, since the expansion provides only the difference between the value of the quantity at two different points in Cartesian space separated by the displacement. The proper transformation is obtained only by employing the transformation coefficients of finite deformation theory, which are defined in terms of the displacement gradients. Finite deformation theory leads to the result that for laterally unconfined, plane waves the Lagrangian and Eulerian radiation pressures are equal with the value (1/4)(2K) along the direction of wave propagation, where (K) is the mean kinetic energy density, and zero in directions normal to the propagation direction. This is contrary to the Langevin result that the Lagrangian radiation pressure in the propagation direction is equal to (2K) and the BST result that the Eulerian radiation pressure in that direction is the momentum flux density.

Finite deformation theory↗

NASA NTRS1994

A time-dependent gyro-kinetic model of thermal ion upflows in the high-latitude F region

Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Loranc, M.↗

NASA NTRS1951

Boundary Layer Effect on the Surface Pressure of an Infinite Cone in Supersonic Flow

The theory of Taylor and Maccoll (Ref,1) gives the surface pressure on an infinite cone in supersonic flow as a function of the cone vertex angle and the free stream Mach number and static pressure for a gas of vanishing viscosity. When a slender conical probe is used together with an impact pressure probe to determine the static pressure and Mach number in a low density gas stream, it is desirable to have some theoretical estimate of the effect of viscous boundary layer on the probe readings. Theoretical and experimental results with respect to impact probes have been presented in Refs. 5 and 6. A simple approximation for a conical probe based on linearized supersonic flow and compressible boundary layer theory is presented here.

Lin, T. C.↗

NASA NTRS1985

Simulation of the electron acoustic instability in the polar cusp

A computer simulation of the beam-driven electron acoustic instability in the polar cusp is used to investigate the temporal evolution of typical cusp electron distribution functions in self-consistently generated wave fields. The simulation results are compared with linear, second-order, and nonlinear theory. The behavior of the instability over a range of plasma parameters, such as beam speed and the cold electron density, is shown. The saturation mechanism of the instability is examined, and the relevance of the results for observations of cusp hiss are discussed.

Lin, C. S.↗

NASA NTRS1992

Electromagnetic ion/ion cyclotron instability - Theory and simulations

Linear theory and 1D and 2D hybrid simulations are employed to study electromagnetic ion/ion cyclotron (EMIIC) instability driven by the relative streaming of two field-aligned ion beams. The characteristics of the instability are studied as a function of beam density, propagation angle, electron-ion temperature ratios, and ion beta. When the propagation angle is near 90 deg the EMIIC instability has the characteristics of an electrostatic instability, while at smaller angles electromagnetic effects play a significant role as does strong beam coupling. The 2D simulations point to a narrowing of the wave spectrum and accompanying coherent effects during the linear growth stage of development. The EMIIC instability is an important effect where ion beta is low such as in the plasma-sheet boundary layer and upstream of slow shocks in the magnetotail.

Winske, D.↗

NASA NTRS1993

Small-scale density irregularities in the nightside Venus ionosphere - Comparison of theory and observations

The theory of the lower-hybrid-drift instability is compared with observations of plasma density irregularities. Marginal stability boundaries (gamma = 0) for the lower-hybrid-drift instability and the occurrence of small-scale density fluctuations are presented as a function of magnetic field B and density n. For plasma density gradient scale lengths in the range 2-10 km, 80-85 percent of the density fluctuations lie in the unstable B/n parameter regime (gamma greater than 0). Stability boundaries for the onset of instability at a wavelength for which the Pioneer Venus Orbiter could measure a Doppler-shifted frequency of 100 Hz are presented. The conditions on B and n for instability in this situation are more stringent than those for marginal stability, especially at low densities (n less than or equal to 5000/cu cm). In general, the instability is most likely to be active in regions of low beta: high magnetic field strength and low density, as found in ionospheric holes or troughs.

Huba, J. D.↗

NASA NTRS2009

The Principle of Energetic Consistency

A basic result in estimation theory is that the minimum variance estimate of the dynamical state, given the observations, is the conditional mean estimate. This result holds independently of the specifics of any dynamical or observation nonlinearity or stochasticity, requiring only that the probability density function of the state, conditioned on the observations, has two moments. For nonlinear dynamics that conserve a total energy, this general result implies the principle of energetic consistency: if the dynamical variables are taken to be the natural energy variables, then the sum of the total energy of the conditional mean and the trace of the conditional covariance matrix (the total variance) is constant between observations. Ensemble Kalman filtering methods are designed to approximate the evolution of the conditional mean and covariance matrix. For them the principle of energetic consistency holds independently of ensemble size, even with covariance localization. However, full Kalman filter experiments with advection dynamics have shown that a small amount of numerical dissipation can cause a large, state-dependent loss of total variance, to the detriment of filter performance. The principle of energetic consistency offers a simple way to test whether this spurious loss of variance limits ensemble filter performance in full-blown applications. The classical second-moment closure (third-moment discard) equations also satisfy the principle of energetic consistency, independently of the rank of the conditional covariance matrix. Low-rank approximation of these equations offers an energetically consistent, computationally viable alternative to ensemble filtering. Current formulations of long-window, weak-constraint, four-dimensional variational methods are designed to approximate the conditional mode rather than the conditional mean. Thus they neglect the nonlinear bias term in the second-moment closure equation for the conditional mean. The principle of energetic consistency implies that, to precisely the extent that growing modes are important in data assimilation, this term is also important.

Cohn, Stephen E.↗

NASA NTRS1999

A Plasma Instability Theory of Gamma-Ray Burst Emission

A plasma instability theory is presented for the prompt radiation from gamma-ray bursts. In the theory, a highly relativistic shell interacts with the interstellar medium through the filamentation and the two-stream instabilities to convert bulk kinetic energy into electron thermal energy and magnetic field energy. The processes are not efficient enough to satisfy the Rankine-Hugoniot conditions, so a shock cannot form through this mechanism. Instead, the interstellar medium passes through the shell, with the electrons radiating during this passage. Gamma-rays are produced by synchrotron self-Compton emission. Prompt optical emission is also produced through this mechanism, while prompt radio emission is produced through synchrotron emission. The model timescales are consistent with the shortest burst timescales. To emit gamma-rays, the shell's bulk Lorentz factor must be $\simg 10(exp 3)$. For the radiative processes to be efficient, the interstellar medium density must satisfy a lower limit that is a function of the bulk Lorentz factor. Because the limits operate as selection effects, bursts that violate them constitute new classes. In particular, a class of optical and ultraviolet bursts with no gamma-ray emission should exist. The lower limit on the density of the interstellar medium is consistent with the requirements of the Compton attenuation theory, providing an explanation for why all burst spectra appear to be attenuated. Several tests of the theory are discussed, as are the next theoretical investigations that should be conducted.

Brainerd, Jerome J.↗

NASA NTRS1972

Viscosity and thermal conductivity of moderately dense gas mixtures.

Derivation of a simple, semitheoretical expression for the initial density dependence of the viscosity and thermal conductivity of gaseous mixtures in terms of the appropriate properties of the pure components and of their interaction quantities. The derivation is based on Enskog's theory of dense gases and yields an equation in which the composition dependence of the linear factor in the density expansion is explicit. The interaction quantities are directly related to those of the mixture extrapolated to zero density and to a universal function valid for all gases. The reliability of the formulation is assessed with respect to the viscosity of several binary mixtures. It is found that the calculated viscosities of binary mixtures agree with the experimental data with a precision which is comparable to that of the most precise measurements.

Wakeham, W. A.↗

NASA NTRS1974

Theory of waves incoherently scattered

Electromagnetic waves impinging upon a plasma at frequencies larger than the plasma frequency, suffer weak scattering. The scattering arises from the existence of electron density fluctuations. The received signal corresponds to a particular spatial Fourier component of the fluctuations, the wave vector of which is a function of the wavelength of the radiowave. Wavelengths short with respect to the Debye length of the medium relate to fluctuations due to non-interacting Maxwellian electrons, while larger wavelengths relate to fluctuations due to collective Coulomb interactions. In the latter case, the scattered signal exhibits a spectral distribution which is characteristic of the main properties of the electron and ion gases and, therefore, provides a powerful diagnosis of the state of the ionosphere.

Bauer, P.↗

NASA NTRS1975

An analytic method to account for drag in the Vinti satellite theory

A quadrature algorithm is presented which employs analytical expressions for the variations of satellite orbital elements caused by air drag. The Hamiltonian is formally preserved and the Jacobi constants of the motion are advanced with time through the variational equations. The atmospheric density profile is written as a fitted exponential function of the eccentric anomaly, which adheres to tabulated data at all altitudes and simultaneously reduces the variational equations to definite integrals with closed form evaluations, whose limits are in terms of the eccentric anomaly. Results are given for two intense air drag satellites and indicate that the satellite ephemerides produced by this method in conjunction with the Vinti program are of very high accuracy.

Watson, J. S.↗

NASA NTRS1976

A BBGKY framework for fluid turbulence

A framework is presented for a systematic kinetic theory of turbulence originating from the Liouville equation for the Fourier coefficients of fluid variables. The real and imaginary parts of these Fourier coefficients play the role that particle coordinates (positions and momenta) play in the BBGKY theory. The basic relations of the problem are the incompressible Navier-Stokes equations in two dimensions with zero viscosity, with the probability distributions of Fourier coefficients rather than moments being the basic variables of the theory. A kinetic equation is derived and shown to possess a number of requirements that any reasonable kinetic equation must have: conservation laws, positive-definite spectral densities, and an H-theorem. The major lack in the theory is any reliable information on the relaxation predicted by the complicated linear operator H. Closure of the hierarchy is achieved by the hypothesis that the five-coefficient correlation function is negligible. Problems associated with inclusion of viscosity and external driving forces are discussed.

Montgomery, D.↗