SEARCH · Engineering Papers

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 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

NASA NTRS1981

Theory of inhomogeneous fluids

A theory of nonuniform liquids is presented which is based on the Yvon-Born-Green equation for the one-particle density and the Ornstein-Zernike equation of an inhomogeneous system. The necessary closure is effected by exploiting the solution of a modified hypernetted-chain equation and making a local approximation on the (highly universal) bridge function. The theory is successfully applied to model problems that have also been studied by direct-simulation methods. A brief generalization to quantum liquids is also given.

Nieminen, R. M.↗

NASA NTRS1982

Low-energy ion distribution functions on a magnetically quiet day at geostationary altitude /L = 7/

Ion energy and pitch angle distribution functions are examined for a magnetically quiet day using averaged data from ATS 6. For both field-aligned and perpendicular fluxes, the populations have a mixture of characteristic energies, and the distribution functions can be fairly well approximated by Maxwellian distributions over three different energy bands in the range 3-600 eV. Pitch angle distributions varying with local time, and energy distributions are used to compute total ion density. Pitch angle scattering mechanisms responsible for the observed transformation of pitch angle distribution are examined, and it is found that a magnetic noise of a certain power spectral density belonging to the electromagnetic ion cyclotron mode near the ion cyclotron frequency can be effective in trapping the field aligned fluxes by pitch angle scattering.

Singh, N.↗

NASA NTRS1986

Plasma waves associated with the first AMPTE magnetotail barium release

Plasma waves observed during the March 21, 1985, AMPTE magnetotail barium release are described. Electron plasma oscillations provided local measurements of the plasma density during both the expansion and decay phases. Immediately after the explosion, the electron density reached a peak of about 400,000/cu cm, and then started decreasing approximately as t to the -2.4 as the cloud expanded. About 6 minutes after the explosion, the electron density suddenly began to increase, reached a secondary peak of about 240/cu cm, and then slowly decayed down to the preevent level over a period of about 15 minutes. The density increase is believed to be caused by the collapse of the ion cloud into the diamagnetic cavity created by the initial expansion. The plasma wave intensities observed during the entire event were quite low. In the diamagnetic cavity, electrostatic emissions were observed near the barium ion plasma frequency, and in another band at lower frequencies. A broadband burst of electrostatic noise was also observed at the boundary of the diamagnetic cavity. Except for electron plasma oscillations, no significant wave activity was observed outside of the diamagnetic cavity.

Gurnett, D. A.↗

NASA NTRS1995

Extreme ultraviolet spectroscopy and photometry of VV Puppis during a high accretion state

We determine the physical properties of the accretion region of the AM Her-type binary VV Puppis using extreme ultraviolet (EUV) medium-resolution spectroscopy and photometry obtained with the Extreme Ultraviolet Explorer (EUVE) observatory. The EUV continuum from VV Pup was detected in the wavelength range from 75 to 135 A and was simultaneously recorded with the Deep Survey/Spectrometer (DS/S) imaging telescope, allowing for the extraction of an accurate light curve. VV Pup appeared to have entered a high-accretion state just prior to the pointed EUVE observations. We use the EUV light curve to infer the diameter of the accretion region (d = 220 km) assuming a hemispherical geometry and a radius of 9000 km for the white dwarf. We perform a model atmosphere analysis and, based on the light curve properties and assuming a distance of 145 pc, we derive an effective temperature of the accretion region in the range 270,000 is less than T(sub eff) is less than 360,000 K and a neutral hydrogen column density in the local interstellar medium of n(sub H) = 1.9 - 3.7 x 10(exp 19)/sq cm. The total EUV/soft X-ray energy radiated by the accretion region is approximately 3.5 x 10(exp 32) ergs/s. Our results provide a first verification of past suggestions that deep heating of the white dwarf surface produces the soft X-ray flux from the polars. We present a possible detection of O VI absortion features, and we suggest that extensive EUVE observations targeting high-accretion events may result in oxygen and heavier element abundance determination in the accretion region.

Vennes, Stephane↗

DOE OSTI2021

Combined microwave and laser Rayleigh scattering diagnostics for pin-to-pin nanosecond discharges

In this work, the temporal decay of electrons produced by an atmospheric pin-to-pin nanosecond discharge operating in the spark regime was measured via a combination of microwave Rayleigh scattering (MRS) and laser Rayleigh scattering (LRS). Due to the initial energy deposition of the nanosecond pulse, a variation in the local gas density occurs on the timescale of electron decay. Thus, the assumption of a constant collisional frequency is no longer applicable when electron number data are extracted from MRS measurements. Furthermore, to recalibrate MRS measurements throughout the electron decay period, temporally resolved LRS measurements of the local gas density were performed over the event duration. The local gas density was calculated to be 30% of the ambient level during the later stages of electron decay, and it recovers at about 1 ms after discharge. A shock front traveling approximately 500 m/s was additionally observed. Coupled with plasma volume calibration via temporally resolved intensified charge-coupled device imaging, the corrected decay curves of the electron number and electron number density are presented with a measured peak electron number density of 4.5 × 1015 cm -3 and a decay rate of ~(0.1–0.35) × 10 7 s -1 . A hybrid MRS and LRS diagnostic technique can be applied for a broad spectrum of atmospheric-pressure microplasmas where a variation in the gas number density is expected due to energy deposition in the discharge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1972

Radiative transfer effects on reflected shock waves. II - Absorbing gas.

Radiative cooling effects behind a reflected shock wave are calculated for an absorbing-emitting gas by means of an expansion procedure in the small density ratio across the shock front. For a gray gas shock layer with an optical thickness of order unity or less the absorption integral is simplified by use of the local temperature approximation, whereas for larger optical thicknesses a Rosseland diffusion type of solution is matched with the local temperature approximation solution. The calculations show that the shock wave will attenuate at first and then accelerate to a constant velocity. Under appropriate conditions the gas enthalpy near the wall may increase at intermediate times before ultimately decreasing to zero. A two-band absorption model yields end-wall radiant-heat fluxes which agree well with available shock-tube measurements.

Su, F. Y.↗

NASA NTRS1994

Kinematics of the ion tail of comet P/Swift-Tuttle

We have obtained long-slit high resolution spectra of the H2O(+) 6199 A complex in the near tail of comet P/Swift-Tuttle. The observations were made using the Hamilton echelle spectrometer fed by the Lick Observatory 0.6 m coude auxiliary telescope. For most of our observations, the spectral slit was aligned along the Sun-tail axis and the cometary nucleus was placed at one end of the slit, giving us spectra having the spatial and spectral resolution needed to measure the radial velocity and velocity dispersion continuously down the cometary tail out to a distance of 4X10(exp 5) km. The radial velocities confirm the earlier more restricted observations by Rauer & Jockers (1993) and by Wyckoff & Lindholm (1994) showing that the tail motions are indeed bulk flows in the antisolar direction. Out to 3X10(exp 5) km in the tail typical bulk flows are at a speed of approximately 30 km/s. The velocity dispersion, (sigma(sub r)), of the H2O(+) lines follows a pattern that is quite systematic; sigma(sub r) is smallest near the cometary nucleus, and steadily increases down the tail. The highest velocity dispersions are found ahead of the nucleus and off the tail axis. These velocity dispersions are equivalent to ion temperatures ranging from 10(exp 5) to 10(exp 6) K. We note a clear anticorrelation between the H2O(+) line intensities (related to the ion density) and the bulk flow and dispersion velocities; direct mass loading of the solar wind by the observed water ions may be responsible. We discuss several approximate equipartition methods used to infer local magnetic fields induced by the interaction of the cometary ions with the solar wind particle/field stream. Typical fields derived are near 50 nT. The measured tailward accelerations are consistent with this order of magnitude B field.

Spinrad, Hyron↗

NASA NTRS2016

First-Principles Modeling Of Electromagnetic Scattering By Discrete and Discretely Heterogeneous Random Media

A discrete random medium is an object in the form of a finite volume of a vacuum or a homogeneous material medium filled with quasi-randomly and quasi-uniformly distributed discrete macroscopic impurities called small particles. Such objects are ubiquitous in natural and artificial environments. They are often characterized by analyzing theoretically the results of laboratory, in situ, or remote-sensing measurements of the scattering of light and other electromagnetic radiation. Electromagnetic scattering and absorption by particles can also affect the energy budget of a discrete random medium and hence various ambient physical and chemical processes. In either case electromagnetic scattering must be modeled in terms of appropriate optical observables, i.e., quadratic or bilinear forms in the field that quantify the reading of a relevant optical instrument or the electromagnetic energy budget. It is generally believed that time-harmonic Maxwell's equations can accurately describe elastic electromagnetic scattering by macroscopic particulate media that change in time much more slowly than the incident electromagnetic field. However, direct solutions of these equations for discrete random media had been impracticable until quite recently. This has led to a widespread use of various phenomenological approaches in situations when their very applicability can be questioned. Recently, however, a new branch of physical optics has emerged wherein electromagnetic scattering by discrete and discretely heterogeneous random media is modeled directly by using analytical or numerically exact computer solutions of the Maxwell equations. Therefore, the main objective of this Report is to formulate the general theoretical framework of electromagnetic scattering by discrete random media rooted in the Maxwell- Lorentz electromagnetics and discuss its immediate analytical and numerical consequences. Starting from the microscopic Maxwell-Lorentz equations, we trace the development of the first principles formalism enabling accurate calculations of monochromatic and quasi-monochromatic scattering by static and randomly varying multiparticle groups. We illustrate how this general framework can be coupled with state-of-the-art computer solvers of the Maxwell equations and applied to direct modeling of electromagnetic scattering by representative random multi-particle groups with arbitrary packing densities. This first-principles modeling yields general physical insights unavailable with phenomenological approaches. We discuss how the first-order-scattering approximation, the radiative transfer theory, and the theory of weak localization of electromagnetic waves can be derived as immediate corollaries of the Maxwell equations for very specific and well-defined kinds of particulate medium. These recent developments confirm the mesoscopic origin of the radiative transfer, weak localization, and effective-medium regimes and help evaluate the numerical accuracy of widely used approximate modeling methodologies.

Michael I Mishchenko↗

NASA NTRS1987

On the redistribution of OB star luminosity and the warming of nearby molecular clouds

Infrared Astronomy Satellite (IRAS) observations of the neighborhoods of six outer-Galaxy HII regions were combined with CO observations to show that most of the far infrared (FIR) luminosity from within approximately 25 to 75 pc of the ionizing stars is contributed by dust in molecular clouds, not by dust in the low-density ionized gas. Dust associated with the clouds is warmed by absorption of UV and visible light from the cluster of stars responsible for the ionization. Most of the OB cluster starlight is not absorbed locally. A fraction of the order of 10% of the OB cluster luminosity is absorbed by nearby molecular clouds and reradiated as FIR light. The luminosity per unit mass for the heated clouds is approximately 3 to 13 solar luminosity/solar mass, approximately one order of magnitude greater than the corresponding ratio for clouds found near clusters without O stars, and two orders of magnitude greater than the ratio for dark clouds heated primarily by the interstellar radiation field. If the observations of clouds near outer-Galaxy HII regions are used to characterize the molecular clouds heated by HII regions in the inner-Galaxy, then at most 30% of the Galaxy's molecular cloud mass is actively engaged in the formation of massive stars at the present time.

Leisawitz, D.↗

DOE OSTI2023

Robust scalable initialization for Bayesian variational inference with multi-modal Laplace approximations

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Local-Field Effects in Linear Response Properties within a Polarizable Frozen Density Embedding Method

In this work, we present a polarizable frozen density embedding (FDE) method for calculating polarizabilities of coupled subsystems. The method (FDE-pol) combines a FDE method with an explicit polarization model such that the expensive freeze/thaw cycles can be bypassed, and approximate nonadditive kinetic potentials are avoided by enforcing external orthogonality between the subsystems. To describe the polarization of the frozen environment, we introduce a Hirshfeld partition-based density-dependent method for calculating the atomic polarizabilities of atoms in molecules, which alleviates the need to fit the atomic parameters to a specific system of interest or to a larger general set of molecules. Further, we show that the Hirshfeld partition-based method predicts molecular polarizabilities close to the basis set limit, and thus, a single basis set-dependent scaling parameter can be introduced to improve the agreement against the reference polarizability data. To test the model, we characterized the uncoupled and coupled response of small interacting molecular complexes. Here, the coupled response properties include the perturbation of the frozen system due to the external perturbation which is ignored in the uncoupled response. We show that FDE-pol can accurately reproduce both the exact uncoupled polarizability and the coupled polarizabilities of the supermolecular systems. Using damped response theory, we also demonstrate that the coupled frequency-dependent polarizability can be described by including local field effects. The results emphasize the necessity of including local-field effects for describing the response properties of coupled subsystems, as well as the importance of accurate atomic polarizability models.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Tunable noninteracting free-energy density functionals for high-energy-density physics applications

In this work, we introduce the concept of a tunable noninteracting free-energy density functional and present two examples realized: (i) via a simple one-parameter convex combination of two existing functionals and (ii) via the construction of a generalized gradient approximation (GGA) enhancement factor that contains one free parameter and is designed to satisfy a set of incorporated constraints. Functional (i), constructed as a combination of the local Thomas–Fermi and a pseudopotential-adapted GGA for the noninteracting free-energy, has already demonstrated its practical usability for establishing the high temperature end of the equation of state of deuterium [Phys. Rev. B 104, 144104 (2021)] and CHON resin [Phys. Rev. E 106, 045207 (2022)] for inertial confinement fusion applications. Hugoniot calculations for liquid deuterium are given as another example of how the application of computationally efficient orbital-free density functional theory (OF-DFT) can be utilized with the employment of the developed functionals. Once the functionals have been tuned such that the OF-DFT Hugoniot calculation matches the Kohn–Sham solution at some low-temperature point, agreement with the reference Kohn–Sham results for the rest of the high temperature Hugoniot path is very good with relative errors for compression and pressure on the order of 2% or less.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1989

Multisatellite and ground-based observations of transient ULF waves

Transient ULF pulsations associated with variations in solar wind plasma density observed by the IMP 8 satellite are presently studied in light of observations obtained during a fortuitous alignment of the AMPTE and Viking satellites with respect to the EISCAT Magnetometer Cross. It is found that the isolated 10-min oscillation in solar wind plasma density produced magnetic field compression oscillations within the magnetosphere at the same frequency, thereby enhancing resonant oscillations at approximately twice the frequency which were already present. Support is seen for the periodic solar wind density variation's exciting of a tailward-traveling large-scale magnetosphere wave train which excites local field line resonant oscillations.

Potemra, T. A.↗

NASA NTRS1979

The growth of aspherical structure in the universe - Is the Local Supercluster an unusual system

The growth and subsequent collapse of homogeneous ellipsoidal perturbations in a uniform expanding background is considered as a simple model for the formation of large-scale aspherical structures in the observed universe. Numerical calculations of the evolution of such perturbations turn out to be well described by an approximate analytic solution of the equations of motion, and simple relationships are found between the initial shape of a perturbation and its shape and kinematic properties at the time of collapse. Perturbations do not change their shape significantly until they reach a density contrast of order unity. As a result, structures with the kinematic properties of the Local Supercluster should form much more commonly in a low-density universe than in a flat universe. The homogeneity of the local Hubble flow, the motion of the Milky Way with respect to the microwave background, and the flattening of the Local Supercluster can be successfully accounted for by these models, provided that the initial perturbation is sufficiently flattened. Viable models are obtained only if the ratio of the lengths of the two smaller axes of the initial perturbation is at least 3:1 in an Einstein-de Sitter universe or at least 1.8:1 in a universe for which the density parameter (Omega) is of order 0.1, when the protocluster pancakes.

White, S. D. M.↗

NASA NTRS1976

The lower ionosphere of Titan

Ionization of the atmosphere of Titan by galactic cosmic rays is a very significant process throughout the altitude range from 100 to 400 km. An approximate form of the Boltzmann equation for cosmic-ray transport has been used to obtain local ionization rates. Models of both ion and neutral chemistry have been employed to compute electron and ion density profiles for three different values of the H2/CH4 abundance ratio. The peak electron density is of the order of 1000 per cu cm. The most abundant positive ions are C2/H9(+) and C3H9(+), while the predicted densities of the negative ions H(-) and CH3(-) are very small (less than one ten-thousandth that of the positive ions). It is suggested that inclusion of the ion chemistry is important in the computation of the H and CH3 density profiles in the lower ionosphere.

Capone, L. A.↗

NASA NTRS1986

Wave speeds in the corona and the dynamics of mass ejections

A disturbance or coronal mass ejection being advected by the solar wind will expand at the fastest local characteristic speed - typically approximately the fast-mode speed. To estimate this characteristic wave speed and the velocity field in the ambient corona, it is necessary to know the magnetic field, temperature, and density. Only the density is known from coronal observations. The temperature, magnetic field, and velocity are not yet directly measured in the outer corona and must be estimated from a model. In this study, it is estimated that the magnetic field, solar wind velocity, and characteristic speeds use the MHD model of coronal expansion between 1 and 5 solar radii (R solar radii) with a dipole magnetic field at the base. This model, for a field strength of about 2 gauss at the base, gives flow speeds at low latitudes (near the heliospheric current sheet) of 250 km/s at 5 R solar radii and, 50 km/s at 2 solar radii, and fast-mode speeds to 400 to 500 km/s everywhere between 2 and 5 solar radii. This suggests that the outer edge of a velocity of mass ejection reported by MacQueen and Fisher (1983) and implies that the acceleration mechanism for coronal mass ejections is other than simple entrainment in the solar wind.

Suess, S. T.↗

NASA NTRS1947

On the one-dimensional theory of steady compressible fluid flow in ducts with friction and heat addition

Steady, diabatic (nonadiabatic), frictional, variable-area flow of a compressible fluid is treated in differential form on the basis of the one-dimensional approximation. The basic equations are first stated in terms of pressure, temperature, density, and velocity of the fluid. Considerable simplification and unification of the equations are then achieved by choosing the square of the local Mach number as one of the variables to describe the flow. The transformed system of equations thus obtained is first examined with regard to the existence of a solution. It is shown that, in general, a solution exists whose calculation requires knowledge only of the variation with position of any three of the dependent variables of the system. The direction of change of the flow variables can be obtained directly from the transformed equations without integration. As examples of this application of the equations, the direction of change of the flow variables is determined for two special flows. In the particular case when the local Mach number m = 1, a special condition must be satisfied by the flow if a solution is to exist. This condition restricts the joint rate of variation of heating, friction, and area at m = 1. Further analysis indicates that when a solution exists at this point it is not necessarily unique. Finally it is shown that the physical phenomenon of choking, which is known to occur in certain simple flow situations, is related to restrictions imposed on the variables by the form of the transformed equations. The phenomenon of choking is thus given a more general significance in that the transformed equations apply to a more general type of flow than has hitherto been treated. (author)

HIGH-SPEED↗