Engineering PapersSearch

SEARCH · Engineering Papers

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

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.

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

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.

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

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.

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

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.

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

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.

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.

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.

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.

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.

Measurement of interfacial thermal resistance in high-energy-density matter

Heat transport across interfaces is a ubiquitous phenomenon with many unresolved aspects. In particular, it is unknown if an interfacial thermal resistance (ITR) occurs in matter with high-energy-density where free electrons dominate the heat conduction. Here, we report on the first experimental evidence that a significant heat barrier is present between two different regions of high-energy-density matter: a strongly heated tungsten wire and a surrounding plastic layer that stays relatively cold. We use diffraction-enhanced imaging to track the time evolution of density discontinuities and reconstruct the temperature evolution in the quasi-stationary stage. The clear signatures of a temperature jump demonstrate the importance of the ITR for strongly heated systems with far-reaching implications for interpreting experiments and applications like inertial confinement fusion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

The impact of data density and data error on the evolution of mesoscale forecast error

The effects of data density and data errors on the accuracy of mesoscale weather forecasts were assessed by simulations of a period in February 1979 when a snowstorm occurred along the U.S. eastern seaboard. The simulations were initiated every 12 hr and the growth of the rms error was tracked as a function of time and varying data densities. A uniform grid of instrumentation separation and perfect boundary conditions were assumed for all the simulations. A bias error representation was defined to account for systematic measurement errors. A brief summary of features of the somewhat complex storm is provided, together with the parameters of the 11 different simulations with high, medium and low data densities. The accuracy of the forecasts were proportional to the data density, although the error decreased for all the simulations over a 24 hr period. The density of the vertical data had a significant impact on the accuracy of the forecast, whereas the horizontal data density did not. Finally, pathways by which errors were transferred among variables were identified.

Warner, T. T.

Temporal behavior of density perturbations in the polar wind

The time-dependent continuity and momentum equations for H(+) are solved in order to study the temporal evolution of density perturbations in the supersonic, collisionless polar wind. The results indicate that the sophisticated computer models of small-scale plasma expansions are indeed applicable to large-scale flow problems. The important prediction of the production of superthermal ions from small-scale simulations by the process of plasma expansion is borne out by the calculations in this paper. Energetic ions are produced for each type of density perturbation through the plasma expansion process. In the case of an extended density depletion, the forward-reverse shock pair forms when the depletion is very strong, but otherwise only the forward shock forms. In the case of localized density depletions, the density cavity evolves into two rarefaction waves joined by a forward-reverse shock pair. In the case of localized density enhancements, the shock pair evolves from compressive waves.

Singh, N.