Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion equation”

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 1,315 records · Page 73
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Energy Models for One-Carrier Transport in Semiconductor Devices

Moment models of carrier transport, derived from the Boltzmann equation, made possible the simulation of certain key effects through such realistic assumptions as energy dependent mobility functions. This type of global dependence permits the observation of velocity overshoot in the vicinity of device junctions, not discerned via classical drift-diffusion models, which are primarily local in nature. It was found that a critical role is played in the hydrodynamic model by the heat conduction term. When ignored, the overshoot is inappropriately damped. When the standard choice of the Wiedemann-Franz law is made for the conductivity, spurious overshoot is observed. Agreement with Monte-Carlo simulation in this regime required empirical modification of this law, or nonstandard choices. Simulations of the hydrodynamic model in one and two dimensions, as well as simulations of a newly developed energy model, the RT model, are presented. The RT model, intermediate between the hydrodynamic and drift-diffusion model, was developed to eliminate the parabolic energy band and Maxwellian distribution assumptions, and to reduce the spurious overshoot with physically consistent assumptions. The algorithms employed for both models are the essentially non-oscillatory shock capturing algorithms. Some mathematical results are presented and contrasted with the highly developed state of the drift-diffusion model.

Jerome, Joseph W.↗
NASA NTRS1999

Laser-Induced Fluorescence Measurements of Translational Temperature and Relative Cycle Number by use of Optically Pumped Trace-Sodium Vapor

Sodium fluorescence induced by a narrow-bandwidth tunable laser has been used to measure temperature, pressure, axial velocity, and species concentrations in wind tunnels, rocket engine exhausts, and the upper atmosphere. Optical pumping of the ground states of the sodium, however, can radically alter the shape of the laser-induced fluorescence excitation spectrum, complicating such measurements. Here a straightforward extension of rate equations originally proposed to account for the features of the pumped spectrum is used to make temperature measurements from spectra taken in pumped vapor. Also determined from the spectrum is the relative fluorescence cycle number, which has application to measurement of diffusion rate and transverse flow velocity, The accuracy of both the temperature and the cycle-number measurements is comparable with that of temperature measurements made in the absence of pumping.

Dobson, Chris C.↗
NASA NTRS2018

Effect of the Reconnection Electric Field on Electron Distribution Functions in the Diffusion Region of Magnetotail Reconnection

Electron distribution functions in the electron diffusion region during symmetric magnetic reconnection are investigated by means of theory and fully kinetic simulations. Crescent‐like striations are formed in distribution functions in the velocity plane perpendicular to the magnetic field. Using an analytical current sheet, we solve the equation of motion for electrons, and derive the shape of a crescent distribution, as a function of the distance from the neutral line, field gradients, and the reconnection electric field. Each crescent is tilted in the velocity plane because of the acceleration by the reconnection electric field, and multiple stripes appear due to multiple meandering bounces. Applying the theory to distribution functions observed in Earth's magnetotail, we deduce the amplitude of the reconnection electric field.

Bessho, N.↗
NASA NTRS1978

Whistler mode noise in Jupiter's inner magnetosphere

A study is made of the amplitude and spectral extent of whistler mode noise in the inner magnetosphere of Jupiter. It is found that the 'hat-shaped' pitch angle distributions of energetic electrons (21 and 31 MeV at L=3) are consistent with those predicted in the presence of a band-limited spectrum of whistler mode noise. The equatorial maximum linear growth rate of parallel propagating whistlers are consistent with those necessary to limit the energetic electron intensities by the whistler mode instability. It is noted that the wave phase speeds before wave reflection can occur at high latitudes and that wave growth is limited to a disk-like region centered around the magnetic equator. The frequency extent of the whistler mode noise spectrum may be estimated by the range of frequencies maximally unstable to equatorial linear growth. A value is found for the spectral density of the broadband whistler mode noise necessary to balance radial diffusion of energetic electrons above the critical range, and an expression is derived for the energetic electron system response to fluctuations about the limiting flux value.

Sentman, D. D.↗
NASA NTRS1979

Magnetic field effects on ionization inhomogeneities advected by motions of the neutral gas

A theory is developed for the ambipolar diffusion, drift, distortion, and production of ionization inhomogeneities within weakly ionized plasma in the presence of the geomagnetic field and motions of the neutral gas. The theory developed includes spatially varying velocity fields of the neutral gas, with a special emphasis on investigating the length scale restrictions necessary for the plasma to be treated as spatially homogeneous. Several effects of the magnetic field on ionization inhomogeneities are discussed. These include the anisotropy of ambipolar diffusion, the distortion and drift due to a gravitational field, the production of irregularities due to the turbulent acceleration field, the drift relative to the neutral gas, and the distortion of ionization inhomogeneities due to the effects of a uniform wind. Pertinent equations are derived for a multispecies plasma containing positive or negative ions, and for the case of electrons and one species of positive ions.

Hill, R. J.↗
NASA NTRS1986

Measurement and prediction of flow through a replica segment of a mildly atherosclerotic coronary artery of man

Pressure distributions were measured along a hollow vascular axisymmetric replica of a segment of the left circumflex coronary artery of man with mildly atherosclerotic diffuse disease. A large range of physiological Reynolds numbers from about 60 to 500, including hyperemic response, was spanned in the flows investigation using a fluid simulating blood kinematic viscosity. Predicted pressure distributions from the numerical solution of the Navier-Stokes equations were similar in trend and magnitude to the measurements. Large variations in the predicted velocity profiles occurred along the lumen. The influence of the smaller scale multiple flow obstacles along the wall (lesion variations) led to sharp spikes in the predicted wall shear stresses. Reynolds number similarity was discussed, and estimates of what time averaged in vivo pressure drop and shear stress might be were given for a vessel segment.

Back, L. H.↗
NASA NTRS1990

Radiometric tests on wet and dry antenna reflector surface panels

The results of X-band noise temperature tests on two types of antenna surface panels are presented. The first type tested was a solid antenna panel, while the second type was a perforated panel with 3/16-in.-diameter holes. Measurements were made at 8.45 GHz using an X-band radiometric system. Included in this article are measured noise temperature contributions from: (1) thermal diffusive white paint on solid and perforated panels, and (2) water sprayed on both painted and unpainted perforated panels. Experiments on perforated panels were restricted to the 3/16-in.-diameter hole panels formerly used on Deep Space Network 64-m antennas. Rigorous calibration equations, applicable to a variety of antenna panel and dichroic plate test configurations, are presented. It was demonstrated that an accurate, stable radiometric measurement system of the type used for the results of this research makes it possible to obtain information that would be much more difficult to obtain using other techniques.

Otoshi, T. Y.↗
NASA NTRS1990

A new class of solutions for interstellar magnetohydrodynamic shock waves

An analysis is presented of the equations of motion for steady MHD shock waves proopagating in interstellar clouds, for boundary conditions that preclude C shocks. In addition to J shocks, in which the neutral fluid component becomes subsonic at an adiabatic jump front, the equations admit a new class of solutions, called C-asterisk shocks, in which the transition to subsonic flow occurs continuously at a sonic point. Numerical methods are developed for computing the structure of J and C-asterisk shocks propagating in diffuse interstellar clouds. The effects of chemical, ionization, and recombination processes are included in this treatment. An alternative numerical method, which uses artificial viscosity to facilitate integration through sonic points, is analyzed and shown to be invalid. A set of exemplary solutions, computed for realistic shock parameters, shows that C-asterisk shocks occur for a broad range of conditions relevant to diffuse interstellar clouds.

Roberge, W. G.↗
NASA NTRS1998

Numerical Simulation of Combustion and Extinction of a Solid Cylinder in Low-Speed Cross Flow

The combustion and extinction behavior of a diffusion flame around a solid fuel cylinder (PMMA) in low-speed forced flow in zero gravity was studied numerically using a quasi-steady gas phase model. This model includes two-dimensional continuity, full Navier Stokes' momentum, energy, and species equations with a one-step overall chemical reaction and second-order finite-rate Arrhenius kinetics. Surface radiation and Arrhenius pyrolysis kinetics are included on the solid fuel surface description and a parameter Phi, representing the percentage of gas-phase conductive heat flux going into the solid, is introduced into the interfacial energy balance boundary condition to complete the description for the quasi-steady gas-phase system. The model was solved numerically using a body-fitted coordinate transformation and the SIMPLE algorithm. The effects of varying freestream velocity and Phi were studied. These parameters have a significant effect on the flame structure and extinction limits. Two flame modes were identified: envelope flame and wake flame. Two kinds of flammability limits were found: quenching at low-flow speeds due to radiative loss and blow-off at high flow speeds due to insufficient gas residence time. A flammability map was constructed showing the existence of maximum Phi above which the solid is not flammable at any freestream velocity.

Tien, J. S.↗
NASA NTRS1999

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗
NASA NTRS2011

Radiation Belt Environment Model: Application to Space Weather and Beyond

Understanding the dynamics and variability of the radiation belts are of great scientific and space weather significance. A physics-based Radiation Belt Environment (RBE) model has been developed to simulate and predict the radiation particle intensities. The RBE model considers the influences from the solar wind, ring current and plasmasphere. It takes into account the particle drift in realistic, time-varying magnetic and electric field, and includes diffusive effects of wave-particle interactions with various wave modes in the magnetosphere. The RBE model has been used to perform event studies and real-time prediction of energetic electron fluxes. In this talk, we will describe the RBE model equation, inputs and capabilities. Recent advancement in space weather application and artificial radiation belt study will be discussed as well.

Fok, Mei-Ching H.↗
NASA NTRS1976

On the structure of contact binaries. I - The contact discontinuity

The problem of the interior structure of contact binaries is reviewed, and a simple resolution of the difficulties which plague the theory is suggested. It is proposed that contact binaries contain a contact discontinuity between the lower surface of the common envelope and the Roche lobe of the cooler star. This discontinuity is maintained against thermal diffusion by fluid flow, and the transition layer is thin to the extent that the dynamical time scale is short in comparison with the thermal time scale. The idealization that the transition layer has infinitesimal thickness allows a simple formulation of the structure equations which are closed by appropriate jump conditions across the discontinuity. The further imposition of the standard boundary conditions suffices to define a unique model for the system once the chemical composition, the masses of the two stars, and the orbital separation are specified.

Shu, F. H.↗
NASA NTRS1978

Artificially created holes in the ionosphere

The artificial creation of ionospheric holes by the release of highly reactive molecules into the F region is investigated. Through ion-atom interchange or charge transfer reactions, H2 or H2O reacts with O(+) to form OH(+) or H2O(+), respectively, which subsequently dissociatively recombines with electrons at a very rapid rate. The diffusion of H2 is also modified by chemical loss to the ambient atomic oxygen atmosphere. The limited spatial and temporal extent of the hole-making process allows several approximations to be made which permit three-dimensional analytic solutions of the continuity equations for the released particles, the O(+) and e(-) densities, and the intermediary molecular ions. A versatile formalism is developed whereby the hole-making capability of virtually any spatial-temporal configuration of released particles can be determined by convolving a set of destruction operators which can be viewed as Green's functions for the problem. As a specific application of the techniques developed, the modification of a winter nighttime ionosphere is described by simulating the release of 1000 kg of water vapor near a height of 300 km.

Mendillo, M.↗
NASA NTRS1985

Application of closed-form solutions to a mesh point field in silicon solar cells

A computer simulation method is discussed that provides for equivalent simulation accuracy, but that exhibits significantly lower CPU running time per bias point compared to other techniques. This new method is applied to a mesh point field as is customary in numerical integration (NI) techniques. The assumption of a linear approximation for the dependent variable, which is typically used in the finite difference and finite element NI methods, is not required. Instead, the set of device transport equations is applied to, and the closed-form solutions obtained for, each mesh point. The mesh point field is generated so that the coefficients in the set of transport equations exhibit small changes between adjacent mesh points. Application of this method to high-efficiency silicon solar cells is described; and the method by which Auger recombination, ambipolar considerations, built-in and induced electric fields, bandgap narrowing, carrier confinement, and carrier diffusivities are treated. Bandgap narrowing has been investigated using Fermi-Dirac statistics, and these results show that bandgap narrowing is more pronounced and that it is temperature-dependent in contrast to the results based on Boltzmann statistics.

Lamorte, M. F.↗
NASA NTRS1993

Time-dependent modulation of galactic cosmic rays by merged interaction regions

Models that solve the one-dimensional, solar modulation equation have reproduced the 11-year galactic cosmic ray using functional representations of global merged interaction regions (MIRs). This study extends those results to the solution of the modulation equation with explicit time dependence. The magnetometers on Voyagers 1 and 2 provide local magnetic field intensities at regular intervals, from which one calculates the ratio of the field intensity to the average local field. These ratios in turn are inverted to form diffusion coefficients. Strung together in radius and time, these coefficents then fall and rise with the strength of the interplanetary magnetic field, becoming representations of MIRs. These diffusion coefficients, calculated locally, propagate unchanged from approx. 10 AU to the outer boundary (120 AU). Inside 10 AU, all parameters, including the diffusion coefficient are assumed constant in time and space. The model reproduces the time-intensity profiles of Voyager 2 and Pioneer 10. Radial gradient data from 1982-1990 between Pioneer 10 and Voyager 2 are about the same magnitude as those calculated in the model. It is also shows agreement in rough magnitude with the radial gradient between Pioneer 10 and 1 AU. When coupled with enhanced, time-dependent solar wind speed at the probe's high latitude, as measured by independent observers, the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source the model also follows Voyager 1's time-intensity profile reasonably well, providing a natural source for the observed negative latitudinal gradients. The model exhibits the 11-year cyclical cosmic ray intensity behavior at all radii, including 1 AU, not just at the location of the spacecraft where the magnetic fields are measured. In addition, the model's point of cosmic ray maximum correctly travels at the solar wind speed, illustrating the well-known propagation of modulation. Finally, at least in the inner heliosphere this model accounts for the delay experienced by lower-rigidity protons in reaching their time-intensity peak. The actual delays in this model, however, are somewhat smaller than the data. In the outer heliosphere the models sees no delays, and the data are ambiguous as to their existence. It appears that strong magnetic field compression regions (merged interaction regions) that are 3-4 times the average field strength can, at least in a helioequatorial band, disrupt effects, such as drifts, that could dominate in quieter magnetic fields. The question remains: Is the heliosphere ever quiet enough to allow such effects to be unambiguously measured, at least in the midlatitudes?

Perko, J. S.↗
NASA NTRS1999

Hierarchical Approach to 'Atomistic' 3-D MOSFET Simulation

We present a hierarchical approach to the 'atomistic' simulation of aggressively scaled sub-0.1 micron MOSFET's. These devices are so small that their characteristics depend on the precise location of dopant atoms within them, not just on their average density. A full-scale three-dimensional drift-diffusion atomistic simulation approach is first described and used to verify more economical, but restricted, options. To reduce processor time and memory requirements at high drain voltage, we have developed a self-consistent option based on a solution of the current continuity equation restricted to a thin slab of the channel. This is coupled to the solution of the Poisson equation in the whole simulation domain in the Gummel iteration cycles. The accuracy of this approach is investigated in comparison to the full self-consistent solution. At low drain voltage, a single solution of the nonlinear Poisson equation is sufficient to extract the current with satisfactory accuracy. In this case, the current is calculated by solving the current continuity equation in a drift approximation only, also in a thin slab containing the MOSFET channel. The regions of applicability for the different components of this hierarchical approach are illustrated in example simulations covering the random dopant-induced threshold voltage fluctuations, threshold voltage lowering, threshold voltage asymmetry, and drain current fluctuations.

Asenov, Asen↗
NASA NTRS1993

NASA/DOD Aerospace Knowledge Diffusion Research Project. XXIV - A general approach to measuring the value of aerospace information products and services

This paper discusses the various approaches to measuring the value of information, first defining the meanings of information, economics of information, and value. It concludes that no general model of measuring the value of information is possible and that the usual approaches, such as cost/benefit equations, have very limited applications. It also concludes that in specific contexts with given goals for newly developed products and services or newly acquired information, there is a basis for its objective valuation. The axioms and inputs for such a model are described and directions for further verification and analysis are proposed.

Brinberg, Herbert R.↗