Search NASA⌕ Search

SEARCH · Search NASA

Results for “density evolution”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Comparison of Computed and Measured Performance of a Pulsed Inductive Thruster Operating on Argon Propellant

Pulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is charged to an initial voltage and then discharged through a coil as a high-current pulse that inductively couples energy into the propellant. The field produced by this pulse ionizes the propellant, producing a plasma near the face of the coil. Once a plasma is formed if can be accelerated and expelled at a high exhaust velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic field. A recent review of the developmental history of planar-geometry pulsed inductive thrusters, where the coil take the shape of a flat spiral, can be found in Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT)[2, 3] and the Faraday Accelerator with Radio-frequency Assisted Discharge (FARAD)[4]. There exists a 1-D pulsed inductive acceleration model that employs a set of circuit equations coupled to a one-dimensional momentum equation. The model was originally developed and used by Lovberg and Dailey[2, 3] and has since been nondimensionalized and used by Polzin et al.[5, 6] to define a set of scaling parameters and gain general insight into their effect on thruster performance. The circuit presented in Fig. 1 provides a description of the electrical coupling between the current flowing in the thruster I1 and the plasma current I2. Recently, the model was upgraded to include an equation governing the deposition of energy into various modes present in a pulsed inductive thruster system (acceleration, magnetic flux generation, resistive heating, etc.)[7]. An MHD description of the plasma energy density evolution was tailored to the thruster geometry by assuming only one-dimensional motion and averaging the plasma properties over the spatial dimensions of the current sheet to obtain an equation for the time-evolution of the total energy. The equation set governing the dynamics of the coupled electrodynamic-current sheet system is composed of first-order, coupled ordinary differential equations that can be easily solved numerically without having to resort to much more complex 2-D finite element plasma simulations.

Polzin, Kurt A.↗
NASA NTRS1994

Evolution of solitary density waves in stellar winds of early-type stars: A simple explanation of discrete absorption component behavior

We model the evolution of a density shell propagating through the stellar wind of an early-type star, in order to investigate the effects of such shells on UV P Cygni line profiles. Unlike previous treatments, we solve the mass, momentum, and energy conservation equations, using an explicit time-differencing scheme, and present a parametric study of the density, velocity, and temperature response. Under the assumed conditions, relatively large spatial scale, large-amplitude density shells propagate as stable waves through the supersonic portion of the wind. Their dynamical behavior appears to mimic propagating 'solitary waves,' and they are found to accelerate at the same rate as the underlying steady state stellar wind (i.e., the shell rides the wind). These hydrodynamically stable structures quantitatively reproduce the anomalous 'discrete absorption component' (DAC) behavior observed in the winds of luminous early-type stars, as illustrated by comparisons of model predictions to an extensive International Ultraviolet Explorer (IUE) time series of spectra of zeta Puppis (O4f). From these comparisons, we find no conclusive evidence indicative of DACs accelerating at a significantly slower rate than the underlying stellar wind, contrary to earlier reports. In addition, these density shells are found to be consistent within the constraints set by the IR observations. We conclude that the concept of propagating density shells should be seriously reconsidered as a possible explanation of the DAC phenomenon in early-type stars.

Waldron, Wayne L.↗
NASA NTRS1974

The evolution of electron density and temperature distributions in the topside ionosphere during magnetic storms

The latitudinal distributions of electron density and temperature during geomagnetic storms in the mid-latitude topside ionosphere are observed to change in a manner than can be related to the evolution of ring current particle populations. The region of auroral precipitation is characterized by correlated increases in electron temperature and density. Equatorwards of this region, there is a broad belt of elevated electron temperatures and depressed electron densities which is usually much broader than any stable auroral red arc distinguishable from the ground, but which is nevertheless the same basic physical phenomenon. The changes of position of this belt can be related to prior bursts of geomagnetic activity and injection of ring current particles into the magnetosphere.

Cole, K. D.↗
DOE OSTI2026

Evolution of charge density wave soft phonon modes in Pd 𝑥 ⁢ErTe 3

Here, we investigated the lattice dynamics of quasi-two-dimensional Pd-intercalated ErTe 3 in relation to its charge-density-wave (CDW) transitions by means of x-ray diffuse and meV-resolution inelastic x-ray scattering. In pristine ErTe 3 , CDW order develops at orthogonal in-plane wave vectors q$^𝑐_1$ =(0,0,0.29) (the 𝑐−CDW ) and q$^𝑎_2 $= (0.31,0,0) (the 𝑎−CDW ), with transition temperatures $𝑇^𝑐_1$ = 270K and $𝑇^𝑎_2$ = 160K, respectively. Remarkably, we observe diffuse x-ray scattering already near the higher transition temperature $𝑇^𝑐_1$ along 𝑎−CDW but at a slightly different wave vector q$^𝑎_1$ = (0.29,0,0). Inelastic x-ray scattering for Pd 0.01 ⁢ErTe 3 shows that a partial phonon softening at q$^𝑎_1$, underscoring the strong competition between ordering tendencies along the nearly equivalent in-plane axes of the orthorhombic lattice. For intercalation levels 𝑥 ≥ 0.02, the 𝑎−CDW state is suppressed. Nevertheless, a similar correlation between phonon softening and diffuse scattering persists along the [100] direction, again observed at q$^𝑎_1$ = (0.29,0,0) and $𝑇^𝑐_1$. These findings suggest that the 𝑎−CDW is fully suppressed for 𝑥 ≥ 0.02, and that the residual diffuse scattering at originates from the partial phonon softening associated with the 𝑐−CDW, reflected by the near equality of the absolute size of q$^𝑐_1$ and q$^𝑎_1$.

Maity, Avishek [Technische Univ. München, Garching↗
NASA NTRS1995

Radial evolution of the energy density of solar wind fluctuations

On the basis of transport theories appropriate to a radially expanding solar wind, we describe new results for the radial evolution of the energy density in solar wind fluctuations at MHD scales. These models include the effects of 'mixing' and driving as well as the possibility of non-isotropic MHD turbulence. Implications of these results for solar wind heating, cosmic ray diffusion and interstellar pick-up ions will also be addressed.

Zank, G. P.↗
DOE OSTI2024

Determining the N -Representability of a Reduced Density Matrix via Unitary Evolution and Stochastic Sampling

The N-representability problem consists in determining whether, for a given p-body matrix, there exists at least one N-body density matrix from which the p-body matrix can be obtained by contraction, that is, if the given matrix is a p-body reduced density matrix (p-RDM). The knowledge of all necessary and sufficient conditions for a p-body matrix to be N-representable allows the constrained minimization of a many-body Hamiltonian expectation value with respect to the p-body density matrix and, thus, the determination of its exact ground state. However, the number of constraints that complete the N-representability conditions grows exponentially with system size, and hence, the procedure quickly becomes intractable for practical applications. This work introduces a hybrid quantum-stochastic algorithm to effectively replace the N-representability conditions. The algorithm consists of applying to an initial N-body density matrix a sequence of unitary evolution operators constructed from a stochastic process that successively approaches the reduced state of the density matrix on a p-body subsystem, represented by a p-RDM, to a target p-body matrix, potentially a p-RDM. The generators of the evolution operators follow the well-known adaptive derivative-assembled pseudo-Trotter method (ADAPT), while the stochastic component is implemented by using a simulated annealing process. The resulting algorithm is independent of any underlying Hamiltonian, and it can be used to decide whether a given p-body matrix is N-representable, establishing a criterion to determine its quality and correcting it. We apply the proposed hybrid ADAPT algorithm to alleged reduced density matrices from a quantum chemistry electronic Hamiltonian, from the reduced Bardeen–Cooper–Schrieffer model with constant pairing, and from the Heisenberg XXZ spin model. In all cases, the proposed method behaves as expected for 1-RDMs and 2-RDMs, evolving the initial matrices toward different targets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1978

Screening materials for smoke evolution

The density of smoke produced by burning or smoldering materials is a factor that affects the ability of occupants to escape from a burning structure and the effectiveness of firefighters. Accordingly, considerable effort is being spent on developing materials that minimize the amount of smoke without impairing their performance characteristics. In the present paper, laboratory tests on such materials are described and evaluated. It is seen that the National Bureau of Standards smoke test and the Ohio State University release-rate test are the most promising for screening materials for smoke evolution. It is shown that, for plastics, smoke obscuration is a more realistic concern than smoke toxicity.

Hilado, C. J.↗
NASA NTRS2010

Modeling of Alkane Oxidation Using Constituents and Species

It is currently not possible to perform simulations of turbulent reactive flows due in particular to complex chemistry, which may contain thousands of reactions and hundreds of species. This complex chemistry results in additional differential equations, making the numerical solution of the equation set computationally prohibitive. Reducing the chemical kinetics mathematical description is one of several important goals in turbulent reactive flow modeling. A chemical kinetics reduction model is proposed for alkane oxidation in air that is based on a parallel methodology to that used in turbulence modeling in the context of the Large Eddy Simulation. The objective of kinetic modeling is to predict the heat release and temperature evolution. This kinetic mechanism is valid over a pressure range from atmospheric to 60 bar, temperatures from 600 K to 2,500 K, and equivalence ratios from 0.125 to 8. This range encompasses diesel, HCCI, and gas-turbine engines, including cold ignition. A computationally efficient kinetic reduction has been proposed for alkanes that has been illustrated for n-heptane using the LLNL heptane mechanism. This model is consistent with turbulence modeling in that scales were first categorized into either those modeled or those computed as progress variables. Species were identified as being either light or heavy. The heavy species were decomposed into defined 13 constituents, and their total molar density was shown to evolve in a quasi-steady manner. The light species behave either in a quasi-steady or unsteady manner. The modeled scales are the total constituent molar density, Nc, and the molar density of the quasi-steady light species. The progress variables are the total constituent molar density rate evolution and the molar densities of the unsteady light species. The unsteady equations for the light species contain contributions of the type gain/loss rates from the heavy species that are modeled consistent with the developed mathematical forms for the total constituent molar density rate evolution; indeed, examination of these gain/loss rates shows that they also have a good quasi-steady behavior with a functional form resembling that of the constituent rate. This finding highlights the fact that the fitting technique provides a methodology that can be repeatedly used to obtain an accurate representation of full or skeletal kinetic models. Assuming success with the modified reduced model, the advantage of the modeling approach is clear. Because this model is based on the Nc rate rather than on that of individual heavy species, even if the number of species increases with increased carbon number in the alkane group, providing that the quasi-steady rate aspect persists, then extension of this model to higher alkanes should be conceptually straightforward, although it remains to be seen if the functional fits would remain valid or would require reconstruction.

Bellan, Jasette↗
NASA NTRS1986

Limits on coronal material in normal galaxies

Measurements of the X-ray surface brightness of a face on disk galaxy M101, have previously been used to place upper limits on the power radiated by a hot corona. Such analysis contrains the effective density of the disk; either it must be so low that the remnants drive a fast hot wind (low radiated power) or so high that the remnant temperature at overlap is low (low X-ray power). These X-ray measurements are here used to constrain the properties of the population of supernova remnants evolving in the disk. This adds a further constraint since young remnants evolving in higher density radiate more of their energy in X-rays, whether or not they eventually overlap to generate a hot corona. The strength of this second limit depends strongly on the density history of the remnants and on the assumed supernova rate. For evaporative evolution the analysis rules out McKee and Ostriker ISM model in particular and evaporative evolution in general unless the supernova rate is at least several times lower than current expectations. For standard Sedov evolutions, the density limit marginally admits evolution in 0.2 cu m, a popular alternative to the McKee and Ostriker model.

Mccammon, D.↗
DOE OSTI2025

Predicting the evolution of biomass bulk density through feedstock preprocessing: Discrete element modeling, regression analysis, and pilot-scale validation

Bulk density is an important material property of biomass feedstocks, influencing handling, storage, transport costs, and conversion efficiency. In this study, predictive regression models for loose and tapped bulk densities of Alamo and Cave-in-Rock switchgrass are developed using a comprehensive dataset generated via calibrated bonded-sphere discrete element method (DEM) simulations. Here, a key contribution of this study is the use of a DEM-based approach, which correlates density with moisture content and particle size distribution parameters and enables analysis across a continuous particle size range, overcoming limitations of purely experimental data. For comparison, regression models are also developed using only experimental data from pilot-scale runs at the Biomass Feedstock National User Facility at Idaho National Laboratory. Validation against pilot-scale data showed reasonable prediction accuracy for both model types, particularly for smaller particle sizes (post-secondary grinding). While the experimental model showed slightly better performance matching the validation data in some cases, the DEM-based model benefits from a much larger dataset, reduced predictor multicollinearity, and continuous parameter coverage, highlighting the utility of validated simulation models for developing robust predictive tools for biomass preprocessing applications.

09 - BIOMASS FUELS↗
NASA NTRS1986

Evolution of lightning flash density and reflectivity structure in a multicell thunderstorm

The radar reflectivity structure and the distribution of lightning in a storm cell was investigated using S-band and UHF-band radar data for six storm cells over Wallops Island. The S-band scans were vertical and continuous, while the UHF data were taken in steps of 2.5 deg elevation. The peak in lightning activity during the study corresponded to a merging of two storm cells. A minimum height of 7 km was found necessary for the appearance of a 40 dBZ core with lightning, which first appears in a multicell thunderstorm at the leading edge of the 50 dBZ core of the cell and between a cell and its decaying neighbor. The lightning moves further into the cell during cell decay and decreases in density. Finally, the lightning is offset horizontally from the precipitation core during cell growth but colocates with the precipitation core as the cell dissipates.

Mazur, V.↗
NASA NTRS1986

Limits on supernova rates and interstellar densities from X-ray observations of M101

Measurements of the X-ray surface brightness of a face-on disk galaxy M101 have previously been used to place upper limits on the power radiated by a hot corona. Such analysis constrains the effective density of the disk; either it must be so low that the remnants drive a fast hot wind (low radiated power), or so high that the remnant temperature at overlap is low (low X-ray power). These X-ray measurements are here used to constrain the properties of the population of supernova remnants evolving in the disk. This adds a further constraint since young remnants evolving in higher density radiate more of their energy in X-rays, whether or not they eventually overlap to generate a hot corona. The strength of this second limit depends strongly on the density history of the remnants (e.g., evaporative vs. nonevaporative evolution) and on the assumed supernova rate. For evaporative evolution the analysis rules out the McKee (1977) and Ostriker (1977) interstellar medium model in particular and evaporative evolution in general unless the supernova rate is at least several times lower than current expectations. For standard Sedov evolutions, the density limit marginally admits 0.2/cu cm, a popular alternative to the McKee and Ostriker model.

Cox, D. P.↗
NASA NTRS1975

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given; the equations are used to study the spacial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS1975

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitation and rotation terms in the equations are of first order in the space variables, the pressure gradient terms are of second order, and the turbulent viscosity term is of third order. The presence of a turbulent viscosity insures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the intial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial inflow plays an important role and allows collapse to occur even when the rotation is large. An approximate solution of the governing partial differential equations is also given, in order to study the spacial distributions of the density and velocity.

Deissler, R. G.↗
NASA NTRS1976

Gravitational collapse of a turbulent vortex with application to star formation

The gravitational collapse of a rotating cloud or vortex is analyzed by expanding the dependent variables in the equations of motion in two-dimensional Taylor series in the space variables. It is shown that the gravitational and rotational terms in the equations are of first order in the space variables, the pressure-gradient terms are of second order, and the turbulent-viscosity term is of third order. The presence of turbulent viscosity ensures that the initial rotation is solid-body-like near the origin. The effect of pressure on the collapse process is found to depend on the shape of the initial density disturbance at the origin. Dimensionless collapse times, as well as the evolution of density and velocity, are calculated by solving numerically the system of nonlinear ordinary differential equations resulting from the series expansions. The axial flow is always inward and allows collapse to occur (axially) even when the rotation is large. An approximate solution of the governing partial differential equations is also given in order to study the spatial distributions of the density and velocity.

Deissler, R. G.↗