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 109 records · Page 6

Ionospheric density enhancement during relativistic electron precipitation

The temporal evolution of the ionospheric density enhancement produced by a widespread relativistic electron precipitation (REP) has been observed with the Chatanika Radar. The REP was associated with a substorm particle energization event, and both the ionospheric absorption and density perturbation exhibited an approximately 90 min periodicity associated with the particles' longitudinal drift. A 80-keV characteristic energy for the precipitating electrons is deduced from ground-based and satellite data. At the maximum of the event, electrons deposited approximately 50 ergs/sq cm per sec in the ionosphere, producing a peak density of 500,000/cu cm at 89 km altitude. At that time the radar observed densities greater than 100,000/cu cm between 70 km and 110 km altitude and riometer absorption at 30 MHz was approximately 12 db.

Foster, J. C.

Early time evolution of negative ion clouds and electron density depletions produced during electron attachment chemical release experiments

Two-dimensional electrostatic particle-in-cell simulations are used to study the early time evolution of electron depletions and negative ion clouds produced during electron attachment chemical releases in the ionosphere. The simulation model considers the evolution in the plane perpendicular to the magnetic field and a three-species plasma that contains electrons, positive ions, and also heavy negative ions that result as a by-product of the electron attachment reaction. The early time evolution (less than the negative ion cyclotron period) of the system shows that a negative charge surplus initially develops outside of the depletion boundary as the heavy negative ions move across the boundary. The electrons are initially restricted from moving into the depletion due to the magnetic field. An inhomogenous electric field develops across the boundary layer due to this charge separation. A highly sheared electron flow velocity develops in the depletion boundary due to E x B and Delta-N x B drifts that result from electron density gradients and this inhomogenous electric field. Structure eventually develops in the depletion boundary layer due to low-frequency electrostatic waves that have growth times shorter than the negative ion cyclotron period. It is proposed that these waves are most likely produced by the electron-ion hybrid instability that results from sufficiently large shears in the electron flow velocity.

Scales, W. A.

The evolution equation for the flame surface density in turbulent premixed combustion

The mean reaction rate in flamelet models for turbulent premixed combustion depends on two basic quantities: a mean chemical rate, called the flamelet speed, and the flame surface density. Our previous work had been primarily focused on the problem of the structure and topology of turbulent premixed flames, and it was then determined that the flamelet speed, when space-averaged, is only weakly sensitive to the turbulent flow field. Consequently, the flame surface density is the key quantity that conveys most of the effects of the turbulence on the rate of energy release. In flamelet models, this quantity is obtained via a modeled transport equation called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact but unclosed formulation for the turbulent Sigma-equation. In the exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the turbulent flame stretch as well as the flame surface density is not available from experiments, and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the alternative approach to get basic information on these fundamental quantities. In the present work, three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, our objective is to extract the turbulent flame stretch from the DNS data base and then perform extensive comparisons with flamelet models. Thanks to the detailed information produced by the DNS-based analysis, it is expected that this type of comparison will not only underscore the shortcomings of current models, but also suggest ways to improve them.

Trouve, Arnaud

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.

The post-pinatubo evolution of stratospheric aerosol surface area density as inferred from SAGE 2

Following the eruption of Mount Pinatubo in June of 1991, the aerosol mass loading of the stratosphere increased from -1 Mt to approximately 30 Mt. This change in aerosol loading was responsible for numerous radiative and chemical changes observed within the stratosphere. As a result, the ability to quantify aerosol properties on a global basis during this period is important. Aerosol surface area density is a critical parameter in governing the rates of heterogeneous reactions, such as ClONO2 plus H2O yields HNO3 plus HOCl, which influence the stratospheric abundance of ozone. Following the eruption of Mt. Pinatubo, measurements by the Stratospheric Aerosol and Gas Experiment (SAGE 2) indicated that the stratospheric aerosol surface area density increased by as much as a factor of 100. Using SAGE 2 multi-wavelength aerosol extinction data, aerosol surface area density as well as mass are derived for the period following the eruption of Mt. Pinatubo through the present.

Poole, L. R.

Spatially Resolved Spectroscopy of the SNR IC443

IC 443 is a supernova remnant of intermediate age, i.e. a few thousand years. It is especially interesting because part of its periphery is expanding into a molecular cloud while other sections are expanding into a typical interstellar medium of much lower density. Since the evolution of a supernova remnant through its various phases is affected by the density of the medium it expands into with the reasonable assumption that the supernova explosion was approximately symmetric we have an opportunity to observe a single object in two phases simultaneously. It was observed by ASCA in April, 1993 for a short period during the PV phase and more thoroughly in a 42 ksec exposure in March, 1994. The latter measurement provides most of the results that have been reported. Most of the analysis took place after the grant ended but is included here for completeness. The data was sent simultaneously to US and Japanese Pls. We worked independently. The software set of FTOOLs was used to construct images and spectra. They were judged to be rather unintuitive and not at all user friendly. I found I was using one FTOOL to read the header to obtain information that would only be provided to another FTOOL. The Japanese investigators were more successful. They analyzed the data and published results more rapidly. The scientific results summarized below are based primarily on their publications. Since IC 443 is an interesting example of a middle aged SNR in which a variety of processes are occurring it is one of a class. IC 443 exhibits shell-like emission in hard X-rays and extended soft X-rays with thin thermal spectra. It resembles SN 1006 in these respects. IC 443 contains hard X-rays in a semi-circular shell surrounding the thermal component. The total hard X-ray flux in the ASCA FOV is only a half of the Ginga hard component; which suggests that the hard X-rays are not confined only in the shell but some are extended larger than the ASCA FOV of eq 1 degree diameter. Japanese investigators examined the spatial structure of the thermal component and analyzed the GIS spectra with a non-equilibrium plasma model, and found no systematic variation of the interstellar absorption across the remnant. Evidence for shock acceleration of cosmic rays to high energies (10 TeV) was found by Keohane. X-ray imaging spectroscopy with ASCA reveals two regions of particularly hard emission: an unresolved source embedded in an extended emission region, and a ridge of emission coincident with the southeastern rim. Both features are located on part of the radio shell where the shock wave is interacting with molecular gas, and together they account for a majority of the emission at 7 keV. Though we would not have noticed it a priori, the unresolved feature is coincident with one resolved by the ROSAT HRI. The ASCA measurements were combined with higher energy data from the XTE and GRO missions and with radio and TeV gamma-ray data to produce a nonthermal multiwavelength spectrum for IC 443 which was fit with a cosmic ray interaction model. This model calculates the cynchrotron, bremsstrahlung, invers Compton, and neutral pion decay emission produced by locally accelerated cosmic ray interacting with ambient matter, soft photon fields, and magnetic fields.

Gorenstein, P.

NICER Observations Reveal that the X-ray Transient MAXI J1348-630 is a Black Hole X-ray Binary

We studied the outburst evolution and timing properties of the recently discovered X-ray transient MAXI J1348−630 as observed with NICER. We produced the fundamental diagrams commonly used to trace the spectral evolution, and power density spectra to study the fast X-ray variability. The main outburst evolution of MAXI J1348−630 is similar to that commonly observed in blackhole transients. The source evolved from the hard state (HS), through hard- and soft-intermediate states, into the soft state in the outburst rise, and back to the HS in reverse during the outburst decay. At the end of the outburst, MAXI J1348−630 underwent two reflares with peak fluxes approximately one and two orders of magnitude fainter than the main outburst, respectively. During the reflares, the source remained in the HS only, without undergoing any state transitions, which is similar to the so-called ‘failed outbursts’. Different types of quasi-periodic oscillations (QPOs) are observed at different phases of the outburst. Based on our spectral-timing results, we conclude that MAXI J1348−630 is a black hole candidate.

L Zhang

The spectrum of density perturbations in an expanding universe

The basic dynamic equations that govern the evolution of perturbations in a Friedmann-Lemaitre universe are derived. General solutions describing the evolution of adiabatic perturbations in the density of matter are obtained, and the choice of the appropriate initial conditions is examined. The various perturbation modes are compared, and the effects of decoupling on the perturbation spectrum are studied. The scheme used to follow the evolution of density perturbations through decoupling is based on an extension of the Eddington approximation to the radiative transfer equation, and is strictly valid in both optically thick and thin limits.

Silk, J.

Local primordial non-Gaussian bias from time evolution

Primordial non-Gaussianity (PNG) is a signature of fundamental physics in the early Universe that is probed by cosmological observations. Here, it is well known that the local type of PNG generates a strong signal in the two-point function of large-scale structure tracers, such as galaxies. This signal, often termed “scale-dependent bias” is a generic feature of modulation of gravitational structure formation by a large-scale mode. It is less well appreciated that the coefficient controlling this signal, b ϕ , is closely connected to the time evolution of the tracer number density. This correspondence between time evolution and local PNG can be simply explained for a universal tracer whose mass function only depends on peak height and, more generally, for nonuniversal tracers in the separate universe picture, which we validate in simulations. We also describe how to recover the bias of tracers subject to a survey selection function and perform a simple demonstration on simulated galaxies. Since the local PNG amplitude in n-point statistics ($f$ NL ) is largely degenerate with the coefficient b ϕ , this proof of concept study demonstrates that Galaxy survey data can allow for more optimal and robust extraction of local PNG information from upcoming surveys.

Sullivan, James M. [University of California, Berk

The Minimum-Mass Surface Density of the Solar Nebula using the Disk Evolution Equation

The Hayashi minimum-mass power law representation of the pre-solar nebula (Hayashi 1981, Prog. Theo. Phys.70,35) is revisited using analytic solutions of the disk evolution equation. A new cumulative-planetary-mass-model (an integrated form of the surface density) is shown to predict a smoother surface density compared with methods based on direct estimates of surface density from planetary data. First, a best-fit transcendental function is applied directly to the cumulative planetary mass data with the surface density obtained by direct differentiation. Next a solution to the time-dependent disk evolution equation is parametrically adapted to the planetary data. The latter model indicates a decay rate of r -1/2 in the inner disk followed by a rapid decay which results in a sharper outer boundary than predicted by the minimum mass model. The model is shown to be a good approximation to the finite-size early Solar Nebula and by extension to extra solar protoplanetary disks.

Davis, Sanford S.

Modeling the densification of metal matrix composite monotape

We present a first model that enables prediction of the density (and its time evolution) of a monotape lay-up subjected to a hot isostatic or vacuum hot pressing consolidation cycle. Our approach is to break down the complicated (and probabilistic) consolidation problem into simple, analyzable parts and to combine them in a way that correctly represents the statistical aspects of the problem, the change in the problem's interior geometry, and the evolving contributions of the different deformation mechanisms. The model gives two types of output. One is in the form of maps showing the relative density dependence upon pressure, temperature, and time for step function temperature and pressure cycles. They are useful for quickly determining the best place to begin developing an optimized process. The second gives the evolution of density over time for any (arbitrary) applied temperature and pressure cycle. This has promise for refining process cycles and possibly for process control. Examples of the models application are given for Ti3Al + Nb, gamma TiAl, Ti6Al4V, and pure aluminum.

Elzey, D. M.

Single-shot spatiotemporal plasma density measurements with a chirped probe pulse

In this work, we present the development and demonstration of a diagnostic for the measurement of the spatial and temporal evolution of plasma density in a single shot. Single-shot Advanced Plasma Probe HolographIc REconstruction (SAPPHIRE) utilizes a chirped probe pulse, a diffractive optical element, a self-referenced interferometer, and an interference bandpass filter to achieve high-fidelity electron density measurements suitable for underdense plasmas that exhibit cylindrical symmetry. The method overcomes limitations in conventional diagnostics, such as reliance on shot-to-shot reproducibility, while capturing plasma dynamics on picosecond timescales with micron-level spatial resolution. The capabilities of SAPPHIRE are demonstrated through measurements of laser-driven plasma channels in helium–nitrogen gas jets. SAPPHIRE demonstrates the formation and expansion of plasma channels in a single shot and the propagation of supersonic ionization fronts while revealing shot-to-shot variations in the plasma profiles. Experimental results are validated against theoretical models and scaling laws, underscoring the robustness and accuracy of this technique. By enabling ultrafast, high-resolution plasma diagnostics in a single exposure, SAPPHIRE represents a transformative advancement in plasma measurement technology.

Grace, Elizabeth S. [Lawrence Livermore National L

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS

Observation of room-temperature charge density wave correlations via coherent phonon spectroscopy in the Sn-doped kagome superconductor CsV 3 ⁢Sb 5

In this study, we perform ultrafast time-resolved reflectivity measurements to track the evolution of charge density wave (CDW) correlations in Sn-doped Kagome superconductor CsV 3⁢ Sb 5−𝑥 ⁢Sn 𝑥 . By extracting the coherent phonon spectrum, we evidence robust signatures of CDW correlations at temperature and doping ranges far beyond the phase boundary of long-range CDW order. Remarkably, we show that short-range CDW correlations survive up to room temperature in 𝑥 = 0.32 Sn-doped CsV 3 ⁢Sb 5 , supported by synchrotron x-ray diffraction measurements. We point out that the introduction of quenched disorder by Sn doping can pin the CDW and form static short-range CDW, which can explain the observed persistent CDW signatures. Our results thus corroborate the ubiquity and robustness of CDW correlations in Sn-doped CsV 3⁢ Sb 5 and provide new insights on the role of disorders on the CDW correlations in AV 3 ⁢Sb 5 family.

Charge density waves

A Model of Reduced Kinetics for Alkane Oxidation Using Constituents and Species for N-Heptane

The reduction of elementary or skeletal oxidation kinetics to a subgroup of tractable reactions for inclusion in turbulent combustion codes has been the subject of numerous studies. The skeletal mechanism is obtained from the elementary mechanism by removing from it reactions that are considered negligible for the intent of the specific study considered. As of now, there are many chemical reduction methodologies. A methodology for deriving a reduced kinetic mechanism for alkane oxidation is described and applied to n-heptane. The model is based on partitioning the species of the skeletal kinetic mechanism into lights, defined as those having a carbon number smaller than 3, and heavies, which are the complement of the species ensemble. For modeling purposes, the heavy species are mathematically decomposed into constituents, which are similar but not identical to groups in the group additivity theory. From analysis of the LLNL (Lawrence Livermore National Laboratory) skeletal mechanism in conjunction with CHEMKIN II, it is shown that a similarity variable can be formed such that the appropriately non-dimensionalized global constituent molar density exhibits a self-similar behavior over a very wide range of equivalence ratios, initial pressures and initial temperatures that is of interest for predicting n-heptane oxidation. Furthermore, the oxygen and water molar densities are shown to display a quasi-linear behavior with respect to the similarity variable. The light species ensemble is partitioned into quasi-steady and unsteady species. The reduced model is based on concepts consistent with those of Large Eddy Simulation (LES) in which functional forms are used to replace the small scales eliminated through filtering of the governing equations; in LES, these small scales are unimportant as far as the overwhelming part of dynamic energy is concerned. Here, the scales thought unimportant for recovering the thermodynamic energy are removed. The concept is tested by using tabular information from the LLNL skeletal mechanism in conjunction with CHEMKIN II utilized as surrogate ideal functions replacing the necessary functional forms. The test reveals that the similarity concept is indeed justified and that the combustion temperature is well predicted, but that the ignition time is over-predicted, a fact traced to neglecting a detailed description of the processes leading to the heavies chemical decomposition. To palliate this deficiency, functional modeling is incorporated into this conceptual reduction in addition to the modeling the evolution of the global constituent molar density, the enthalpy evolution of the heavies, the contribution to the reaction rate of the unsteady lights from other light species and from the heavies, the molar density evolution of oxygen and water, and the mole fractions of the quasisteady light species. The model is compact in that there are only nine species-related progress variables. Results are presented showing the performance of the model for predicting the temperature and species evolution. The model reproduces the ignition time over a wide range of equivalence ratios, initial pressure, and initial temperature.

Harstad, Kenneth G.

Control-Affine Schrödinger Bridge and Generalized Bohm Potential

From a stochastic control perspective, the Schrödinger bridge is a density-valued continuous curve parameterized by time that connects a given pair of initial and terminal probability densities via minimum effort controlled Brownian motion. The control-affine Schrödinger bridge extends this idea to a generic control-affine Itô diffusion, possibly with an additive state cost. Here, in this letter, we recast the necessary conditions of optimality for the control-affine Schrödinger bridge problem as a two point boundary value problem for a quantum mechanical Schrödinger PDE with complex potential. This complex-valued potential is a generalization of the real-valued Bohm potential in quantum mechanics. Our derived potential is akin to the optical potential in nuclear physics where the real part of the potential encodes elastic scattering (transmission of wave function), and the imaginary part encodes inelastic scattering (absorption of wave function). The key takeaway is that the process noise that drives the evolution of probability densities induces an absorbing medium in the evolution of wave function. These results make new connections between control theory and non-equilibrium statistical mechanics through the lens of quantum mechanics.

Markov processes