Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “density distribution”

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 379 records · Page 21

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗

Winds From Luminous Late-Type Stars: The Effects of Nonlinear Alfven Waves - 1

We present the results of magnetohydrodynamic (MHD) modeling of winds from luminous late-type stars using a 2.5-dimensional, nonlinear MHD computer code. We assume that the wind is generated within an initially hydrostatic atmosphere and is driven by torsional Alfven waves generated at the stellar surface. Two cases of atmospheric topology are considered: case I has longitudinally uniform density distribution and isotropic radial magnetic field over the stellar surface, and case II has an isotropic, radial magnetic field with a transverse density gradient, which we refer to as an "atmospheric hole." We use the same set of boundary conditions for both models. The calculations are designed to model a cool luminous star, for which we assume an initial hydrostatic pressure scale height of 0.072 Stellar Radius, an Alfven wave speed of 92 km/s at the surface, and a wave period of 76 days, which roughly corresponds with the convective turnover time. For case I the calculations produce a wind with terminal velocity of about 22 km/s and a mass loss rate comparable to the expected value of 10(exp -6) Solar Mass/yr. For case II we predict a two-component wind: a fast (25 km/s) and relatively dense wind outside of the atmospheric hole and a slow (1.5 km/s), rarefied wind inside of the hole.

Airapetian, V. S.↗

Rocket Measurement of a Daytime Electron Density Profile up to 620 Kilometers

On April 27, 1961 at 1502 EST a four-stage research rocket was fired from Wallops Island, Virginia, to measure the ionospheric electron density distribution by means of Seddon's CW propagation technique. This experimental technique is based upon the dispersive Doppler effect measured at two harmonically related frequencies, in this case f = 12.267 Mc and 6f = 73.6 Mc. The electron density profile measured above the peak of the F2 region is representative of a diffusive-equilibrium distribution in an isothermal ionosphere having a temperature of 1640 deg +/- 90 deg K. This result, when compared with satellite and other data, indicates that the upper ionosphere is in thermodynamic equilibrium.

Jackson, J. E.↗

Electron bunch dynamics and emission in particle-in-cell simulations of relativistic laser–solid interactions: On density artifacts, collisions, and numerical dispersion

Sub-optical-cycle dynamics of dense electron bunches in relativistic-intensity laser–solid interactions lead to the emission of high-order harmonics and attosecond light pulses. The capacity of particle-in-cell simulations to accurately model these dynamics is essential for the prediction of emission properties because the attosecond pulse intensity depends on the electron density distribution at the time of emission and on the temporal distribution of individual electron Lorentz-factors in an emitting electron bunch. Here, we show that in one-dimensional collisionless simulations, the peak density of the emitting electron bunch increases with the increase in the spatial resolution of the simulation grid. When collisions are added to the model, the peak electron density becomes independent of the spatial resolution. Collisions are shown to increase the spread of the peaks of Lorentz-factors of emitting electrons in time, especially in the regimes far from optimum generation conditions, thus leading to lower intensities of attosecond pulses as compared to those obtained in collisionless simulations.

Physics↗

Determination of plasma parameters from soft X-ray images for coronal holes /open magnetic field configurations/ and coronal large-scale structures /extended closed-field configurations/

In connection with high-quality solar soft X-ray images the 'quiet' features of the inner corona have been separated into two sharply different components, including the strongly reduced emission areas or coronal holes (CH) and the extended regions of looplike emission features or large-scale structures (LSS). Particular central meridian passage observations of the prominent CH1 on August 21, 1973, are selected for a quantitative study. Histogram photographic density distributions for full-disk images at other central meridian passages of CH 1 are also presented, and the techniques of converting low photographic density data to deposited energy are discussed, with particular emphasis on the problems associated with the CH data.

Maxson, C. W.↗

IMAGE EUV Observations and Modeling of the Plasmaspheric Density Trough Associated with the 24 May 2000 Geomagnetic Storm

The IMAGE EUV imager observed a plasmaspheric density trough in association with a geomagnetically active period on 24 May 2000. At EUV wavelengths, this density trough appeared as an Archimedes spiral extending from Earth's shadow to approximately 1800 MLT. We present an analysis of this density trough using simulated EUV images. Observational EUV images are subjected to edge analysis to establish the plasmapause L-shell and the location of the density trough in terms of L-shell, MLT extent, and radial width. The plasmaspheric density distribution is modeled using both static and dynamic models for the plasmasphere. The background plasmasphere is then numerically simulated using the 4-parameter plasmaspheric density model contained within the Global Core Plasma Model (GCPM) [Gallagher et al., 20001 and the Dynamic Global Core Plasma Model (DGCPM). Simulated EUV images of the model plasmasphere are produced once an artificial density depletion, matching the observed MLT extent and width, has been removed. Once the azimuthal extent and width of the trough have been simulated, the depth of the artificial density depletion is iteratively adjusted to produce simulated EUV images that approximate observation. The results of this analysis and discussion of possible origins for this density trough will be presented.

Adrian, M.L.↗

The Large-Scale Plasmaspheric Density Trough Associated With the 24 May 2000 Geomagnetic Storm: IMAGE EUV Observations and Global Core Plasma Modeling

The IMAGE EUV imager observed a plasmaspheric density, trough in association with a geomagnetically active period on 24 May 2000. In EUV, this density, trough appears as an Archimedes spiral extending from Earth's shadow to approximately 1900 MLT. We present an analysis of this density trough using simulated EUV images, Observational EUV images are subjected to edge analysis to establish the plasmapause L-shell and the location of the density trough in terms of L-shell, MLT extent, and radial width. The plasmaspheric density distribution is modeled using both static and dynamic models for the plasmasphere. The background plasmasphere is then numerically simulated using the 4-parameter plasmaspheric density model contained within the Global Core Plasma Model (GCPM) and the Dynamic Global Core Plasma Model (DGCPM). Simulated EUV images of the model plasmasphere are produced once an artificial density, depletion, matching the observed MLT extent and width, has been removed. Once the azimuthal extent and width of the trough have been simulated, the depth of the artificial density depletion is iteratively adjusted to produce simulated EUV images that approximate observation. The results of this analysis and discussion of possible origins for this density trough will be presented.

Adrian, M. L.↗

Density estimation via measure transport: Outlook for applications in the biological sciences

Abstract One among several advantages of measure transport methods is that they allow or a unified framework for processing and analysis of data distributed according to a wide class of probability measures. Within this context, we present results from computational studies aimed at assessing the potential of measure transport techniques, specifically, the use of triangular transport maps, as part of a workflow intended to support research in the biological sciences. Scenarios characterized by the availability of limited amount of sample data, which are common in domains such as radiation biology, are of particular interest. We find that when estimating a distribution density function given limited amount of sample data, adaptive transport maps are advantageous. In particular, statistics gathered from computing series of adaptive transport maps, trained on a series of randomly chosen subsets of the set of available data samples, leads to uncovering information hidden in the data. As a result, in the radiation biology application considered here, this approach provides a tool for generating hypotheses about gene relationships and their dynamics under radiation exposure.

gene expression data↗

Theoretical model of vertical distributions of CO and CH4 in the mesosphere and upper stratosphere.

Carbon compounds (CO and CO2) and methane are added to an existing computer program (Shimazaki and Laird, 1970, 1972), which calculates diurnal variations in density distributions of pure oxygen compounds, hydrogen compounds, and nitrogen compounds. For the present model calculations, the original program was modified to calculate the steady-state profiles for all constituents, and then these solutions were used as the initial condition for solving the time-dependent equations. The effects of chemical reactions and vertical eddy diffusion transport on the CH4 and CO concentrations above 40 km, where one-dimensional treatment may not be too inappropriate, are discussed. Some of the 15 constituents included in the computer program interact with CO and CH4 chemically, and CH4 is more sensitive to the eddy diffusion coefficient than any other constituent. Therefore it is easier to determine the most appropriate profile of the eddy diffusion coefficient by comparing the model profile with the observations for CH4.

Shimazaki, T.↗

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth-order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗

Electron Density Changes across the Pressure-Induced Iron Spin Transition

High-pressure single-crystal x-ray diffraction is used to experimentally map the electron-density distribution changes in (Fe,Mg)O as ferrous iron undergoes a pressure-induced transition from high- to low-spin states. As the bulk density and elasticity of magnesiowustite—one of the dominant mineral phases of Earth’s mantle—are affected by this electronic transition, our results have applications to geophysics as well as to validating first-principles calculations. Furthermore, the observed changes in diffraction intensities indicate a spin-transition-induced change in orbital occupancies of the Fe ion in general accord with crystal-field theory, illustrating the use of electron density measurements for characterizing high-pressure d-block chemistry and motivating further studies characterizing chemical bonding under pressure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Achieving Higher Order Accuracy in Space in Hydrodynamic Simulations of Self-Gravitating Gas

Modern astrophysical simulation codes employ a variety of numerical algorithms capable of achieving higher-order accuracy in both space and time. Albeit they succeed in achieving an effective higher spatial resolution and in suppressing the numerical damping of waves, to our knowledge, all current astrophysical simulations invoking self-gravity are limited to second-order accuracy in space. If we can devise an algorithm to evaluate self-gravity with a higher-order spatial accuracy, we can better the evaluation of the gravitational acceleration and gravitational energy release which dictate the evolution of many astrophysical systems. Herein, we present a numerical algorithm for self-gravitating hydrodynamics capable of achieving fourth order accuracy for a given density distribution on a Cartesian uniform grid. First, we derive the cell-averaged gravitational potential at fourth-order accuracy from the cell-averaged density by solving the Poisson equation. Next, we obtain the cell average of the product of the density and gravitational acceleration, which differs from the cell-averaged density multiplied by the cell-averaged gravitational acceleration. We then show the verification of the algorithm by applying it to critical test problems: (1) maintaining equilibria of self-gravitating slabs, even upon advection, (2) evolving a polytropic sphere with a massive power-law envelope, and (3) conservation of specific entropy during the propagation of a sound wave.

79 ASTRONOMY AND ASTROPHYSICS↗

Exploring the Mass-loss Histories of the Red Supergiants

We report mid- to far-infrared imaging and photometry from 7 to 37 μm with SOFIA/FORCAST and 2 μm adaptive optics imaging with LBTI/LMIRCam of a large sample of red supergiants (RSGs) in four Galactic clusters: RSGC1, RSGC2 = Stephenson 2, RSGC3, and NGC 7419. The RSGs in these clusters cover their expected range in luminosity and initial mass from ≈9 to more than 25 M {sub ⊙}. The population includes examples of very late-type RSGs such as MY Cep, which may be near the end of the RSG stage, high-mass-losing maser sources, yellow hypergiants, and post-RSG candidates. Many of the stars and almost all of the most luminous have spectral energy distributions (SEDs) with extended infrared excess radiation at the longest wavelengths. To best model their SEDs, we use the DUSTY code with a variable radial density distribution function to estimate their mass-loss rates. Our M-dot –luminosity relation for 42 RSGs basically follows the classical de Jager curve, but at luminosities below 10{sup 5} L{sub ⊙}, we find a significant population of RSGs with M-dot below the de Jager relation. At luminosities above 10{sup 5} L {sub ⊙}, there is a rapid transition to higher mass-loss rates that approximates and overlaps the de Jager curve. We recommend that instead of using a linear relation or single curve, the empirical M-dot –luminosity relation is better represented by a broad band. Interestingly, the transition to much higher M-dot at about 10{sup 5} L {sub ⊙} corresponds approximately to an initial mass of 18–20 M {sub ⊙}, which is close to the upper limit for RSGs becoming Type II supernovae.

79 ASTRONOMY AND ASTROPHYSICS↗

The 2025 Evaluation of Experimental Thermonuclear Reaction Rates (ETR25)

This work describes the formalism for estimating thermonuclear reaction rates for astrophysical applications, emphasizing modern statistical approaches such as Monte Carlo sampling and Bayesian models. We discuss related topics including the calculation of resonance energies from nuclear Q values, indirect estimates of particle partial widths, and matching of reaction rates at elevated temperatures to statistical model results. We have evaluated available experimental data on cross sections, resonance energies and strengths, partial widths, lifetimes, spin-parities, and spectroscopic factors. Based on these results, we have estimated numerical values of 78 experimental charged-particle thermonuclear reaction rates for target nuclei in the A = 2–40 mass region, for temperatures ranging from 1 MK to 10 GK. For each reaction, three rate values are provided: low, median, and high, corresponding to the 16th, 50th, and 84th percentiles, respectively, of the cumulative reaction rate probability density distribution. Additionally, we present the factor uncertainty of each rate at each temperature grid point. These results enable users to sample the reaction rate probability density in nucleosynthesis calculations, facilitating uncertainty estimates of nuclidic abundances. The rates presented here refer to their laboratory values. For use in stellar model simulations, these values need to be corrected for the effects of thermal excitations of the interacting nuclei. For each reaction, we include graphs that illustrate the fractional contributions to the overall reaction rate along with the associated uncertainty. These visuals are designed to assist both stellar modelers and nuclear experimentalists by identifying the primary sources of rate uncertain=^texttx);ty at specific stellar temperatures. A graphical comparison with earlier Monte Carlo rates is also provided.

Nuclear astrophysics↗

The Aladdin experiment. I - Dynamics.

A coordinated rocket experiment was carried out on Nov. 20, 1970 at Eglin AFB, Fla. (30 deg N) to determine atmospheric layering and density distributions of ions and neutrals. The program measured vertical profiles between 70 and 150 km altitude of total neutrals, ions and electrons, individual species, temperature, density, diffusion coefficients, horizontal winds, and turbulent structure. A dynamic computer code model using these parameters was compared with the observed profiles. Wind fields were measured from seven chemical trails, separated up to 60 km horizontally and 600 sec in time. Kinetic temperatures were measured from fluorescence of chemically-released AlO, ranging from 440 to 680 K between 126 and 161 km. Densities derived from diffusion of the AlO ranged from 10 to 5 times 10 to the -12th power per cu cm between 126 and 139 km. A rocket-borne pitot probe showed wavelike density deviations from current standard atmosphere models, amounting to -25% at 80 and 110 km and to +25% at 100 and 125 km. The derived temperature profile from the pitot probe showed a double minimum at 74 and 102 km.

Rosenberg, N. W.↗

Dusty Plasmas in Planetary Magnetospheres Award

This is my final report for the grant Dusty Plasmas in Planetary Magnetospheres. The funding from this grant supported our research on dusty plasmas to study: a) dust plasma interactions in general plasma environments, and b) dusty plasma processes in planetary magnetospheres (Earth, Jupiter and Saturn). We have developed a general purpose transport code in order to follow the spatial and temporal evolution of dust density distributions in magnetized plasma environments. The code allows the central body to be represented by a multipole expansion of its gravitational and magnetic fields. The density and the temperature of the possibly many-component plasma environment can be pre-defined as a function of coordinates and, if necessary, the time as well. The code simultaneously integrates the equations of motion with the equations describing the charging processes. The charging currents are dependent not only on the instantaneous plasma parameters but on the velocity, as well as on the previous charging history of the dust grains.

Horanyi, Mihaly↗

Resolving Faint Structures in the Debris Disk Around TWA 7* Tentative Detections of an Outer Belt, a Spiral Arm, and a Dusty Cloud

Context. Debris disks are the intrinsic by-products of the star and planet formation processes. Most likely due to instrumental limitations and their natural faintness, little is known about debris disks around low mass stars, especially when it comes to spatially resolved observations. Aims. We present new VLT/SPHERE IRDIS dual-polarization imaging (DPI) observations in which we detect the dust ring around the M2 spectral type star TWA 7. Combined with additional angular differential imaging observations we aim at a fine characterization of the debris disk and setting constraints on the presence of low-mass planets. Methods. We modeled the SPHERE DPI observations and constrain the location of the small dust grains, as well as the spectral energy distribution of the debris disk, using the results inferred from the observations, and performed simple N-body simulations. Results. We find that the dust density distribution peaks at ~0.72′′ (25 au), with a very shallow outer power-law slope, and that the disk has an inclination of ~13° with a position angle of ~91° east of north. We also report low signal-to-noise ratio detections of an outer belt at a distance of ~1.5′′ (~52 au) from the star, of a spiral arm in the southern side of the star, and of a possible dusty clump at 0.11′′. These findings seem to persist over timescales of at least a year. Using the intensity images, we do not detect any planets in the close vicinity of the star, but the sensitivity reaches Jovian planet mass upper limits. We find that the SED is best reproduced with an inner disk at ~0.2′′ (~7 au) and another belt at 0.72′′ (25 au). Conclusions. We report the detections of several unexpected features in the disk around TWA 7. A yet undetected 100 Solar Mass planet with a semi-major axis at 20−30 au could possibly explain the outer belt as well as the spiral arm. We conclude that stellar winds are unlikely to be responsible for the spiral arm.

VLT/SPHERE IRDIS↗