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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Simulation of Ultra-Small MOSFETs Using a 2-D Quantum-Corrected Drift-Diffusion Model

We describe an electronic transport model and an implementation approach that respond to the challenges of device modeling for gigascale integration. We use the density-gradient (DG) transport model, which adds tunneling and quantum smoothing of carrier density profiles to the drift-diffusion model. We present the current implementation of the DG model in PROPHET, a partial differential equation solver developed by Lucent Technologies. This implementation approach permits rapid development and enhancement of models, as well as run-time modifications and model switching. We show that even in typical bulk transport devices such as P-N diodes and BJTs, DG quantum effects can significantly modify the I-V characteristics. Quantum effects are shown to be even more significant in small, surface transport devices, such as sub-0.1 micron MOSFETs. In thin-oxide MOS capacitors, we find that quantum effects may reduce gate capacitance by 25% or more. The inclusion of quantum effects in simulations dramatically improves the match between C-V simulations and measurements. Significant quantum corrections also occur in the I-V characteristics of short-channel MOSFETs due to the gate capacitance correction.

Biegel, Bryan A.↗
NASA NTRS1994

The Nonlinear Evolution of the Inviscid Secondary Instability of Streamwise Vortex Structures

The weakly nonlinear evolution of an inviscid marginally unstable wave growing on a boundary layer supporting a streamwise vortex structure is investigated. The nonlinear growth of the wave is found to be controlled by the diffusion layer located at the edge of the critical layer associated with the wave. The evolution equation is found to depend on the upstream history of the wave and the solution of the equation suggests that the wave either restructures the mean state so as to make it stable or develops a singularity at a finite distance downstream of the point of neutral stability.

Blackaby, Nicholas D.↗
NASA NTRS2005

Removal of organic contaminants by RO and NF membranes

Rejection characteristics of organic and inorganic compounds were examined for six reverse osmosis (RO) membranes and two nanofiltration (NF) membranes that are commercially available. A batch stirred-cell was employed to determine the membrane flux and the solute rejection for solutions at various concentrations and different pH conditions. The results show that for ionic solutes the degree of separation is influenced mainly by electrostatic exclusion, while for organic solutes the removal depends mainly upon the solute radius and molecular structure. In order to provide a better understanding of rejection mechanisms for the RO and NF membranes, the ratio of solute radius (r(i,s)) to effective membrane pore radius (r(p)) was employed to compare rejections. An empirical relation for the dependence of the rejection of organic compounds on the ratio r(i,s)/r(p) is presented. The rejection for organic compounds is over 75% when r(i,s)/r(p) is greater than 0.8. In addition, the rejection of organic compounds is examined using the extended Nernst-Planck equation coupled with a steric hindrance model. The transport of organic solutes is controlled mainly by diffusion for the compounds that have a high r(i,s)/r(p) ratio, while convection is dominant for compounds that have a small r(i,s)/r(p) ratio. c2005 Elsevier B.V. All rights reserved.

NASA Discipline Life Sciences Technologies↗
NASA NTRS2002

Computation of Turbulent Heat Transfer on the Walls of a 180 Degree Turn Channel With a Low Reynolds Number Reynolds Stress Model

The Low Reynolds number version of the Stress-omega model and the two equation k-omega model of Wilcox were used for the calculation of turbulent heat transfer in a 180 degree turn simulating an internal coolant passage. The Stress-omega model was chosen for its robustness. The turbulent thermal fluxes were calculated by modifying and using the Generalized Gradient Diffusion Hypothesis. The results showed that using this Reynolds Stress model allowed better prediction of heat transfer compared to the k-omega two equation model. This improvement however required a finer grid and commensurately more CPU time.

Ameri, A. A.↗
NASA NTRS2011

Analysis of the Radiative Transfer Equation with Highly Asymmetric Phase Function

This paper considers a scalar radiative transfer problem with high scattering anisotropy, Two computational methods are presented based on decomposition of the diffuse light field into a regular and anisotropic part. The first algorithm (DOMAS) singles out the anisotropic radiance in the forward scattering peak using the Small-Angle Modification of RTE. The second algorithm (DOM2+) separates the single scattering radiance as an anisotropic part, which largely defines the fine detail of the total radiance in the backscattering directions. In both cases, the anisotropic part is represented analytically. With anisotropy subtraction, the regular part of the signal. which requires a numerical solution, is essentially smoothed as a function of angles. Further, the transport equation is obtained for the regular part that contains an additional source function from the anisotropic part of the signal. This equation is solved with the discrete ordinates method. A conducted numerical analysis of this work showed that algorithm DOMAS has a strong advantage as compared to the standard discrete ordinates method for simulation of the radiance transmission, and DOM2 + is the best of the three for the reflection computations. Both algorithms offer at least a factor of three acceleration of convergence of the azimuthal series for highly anisotropic phase functions.

Korkin, Sergey V.↗
NASA NTRS2005

Ring Current Dynamic in the Presence of EMIC Waves

The effect of EMIC waves, generated by a positive ion temperature anisotropy on Earth s RC ions is one of the best known examples of wave-particle interaction in the magnetosphere and the most controversial mechanism of RC losses. Under certain conditions, relativistic electrons with energy 21 MeV can be removed from the outer RB by EMIC wave scattering during a magnetic storm much faster than by any other loss mechanisms. That is why the calculation of EMIC waves is a very critical part of the NASA LWS program. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch- angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) LHW are included as an additional factor that contributes to saturation process of EMIC waves; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗
NASA NTRS1989

Assessment of numerical methods for the solution of fluid dynamics equations for nonlinear resonance systems

The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Przekwas, A. J.↗
NASA NTRS1977

SAS-2 observations of the diffuse gamma radiation in the galactic latitude interval from 10 to 90 deg in both hemispheres

An analysis of all the second Small Astronomy Satellite (SAS-2) gamma-ray data for galactic latitudes higher than 10 deg in both hemispheres has shown that the intensity varies with galactic latitude, being larger near 10 deg than 90 deg. For energies above 100 MeV the gamma-ray data are consistent with a latitude distribution of the form I(b) = C1 + C2/sin b, with the second term being dominant. This result suggests that the radiation above 100 MeV is coming largely from local regions of the galactic disk. Between 35 and 100 MeV, a similar equation is also a good representation of the data, but here the two terms are comparable. These results indicate that the diffuse radiation above 35 MeV consists of two parts, one with a relatively hard galactic component and the other an isotropic steep spectral component which extrapolates back well to the low-energy (less than 10 MeV) diffuse radiation. The steepness of the diffuse isotropic component places significant constraints on possible theoretical models of this radiation.

Fichtel, C. E.↗
NASA NTRS2005

Ring Current-Electromagnetic Ion Cyclotron Waves Coupling

The effect of Electromagnetic Ion Cyclotron (EMIC) waves, generated by ion temperature anisotropy in Earth s ring current (RC), is the best known example of wave- particle interaction in the magnetosphere. Also, there is much controversy over the importance of EMIC waves on RC depletion. Under certain conditions, relativistic electrons, with energies 21 MeV, can be removed from the outer radiation belt (RB) by EMIC wave scattering during a magnetic storm. That is why the calculation of EMIC waves must be a very critical part of the space weather studies. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch-angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) the EMIC waves saturation due to their modulation instability and LHW generation are included as an additional factor that contributes to this process; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗
NASA NTRS1982

Analytical modeling of operating characteristics of premixing-prevaporizing fuel-air mixing passages. Volume 1: Analysis and results

A model for predicting the distribution of liquid fuel droplets and fuel vapor in premixing-prevaporizing fuel-air mixing passages of the direct injection type is reported. This model consists of three computer programs; a calculation of the two dimensional or axisymmetric air flow field neglecting the effects of fuel; a calculation of the three dimensional fuel droplet trajectories and evaporation rates in a known, moving air flow; a calculation of fuel vapor diffusing into a moving three dimensional air flow with source terms dependent on the droplet evaporation rates. The fuel droplets are treated as individual particle classes each satisfying Newton's law, a heat transfer, and a mass transfer equation. This fuel droplet model treats multicomponent fuels and incorporates the physics required for the treatment of elastic droplet collisions, droplet shattering, droplet coalescence and droplet wall interactions. The vapor diffusion calculation treats three dimensional, gas phase, turbulent diffusion processes. The analysis includes a model for the autoignition of the fuel air mixture based upon the rate of formation of an important intermediate chemical species during the preignition period.

Anderson, O. L.↗
NASA NTRS1977

SAS-2 observations of the diffuse gamma radiation in the galactic latitude interval 10 deg absolute b or equal to 90 deg

An analysis of all of the second small astronomy satellite gamma-ray data for galactic latitudes with the absolute value of b 10 deg has shown that the intensity varies with galactic latitude, being larger near 10 deg than 90 deg. For energies above 100 MeV the gamma-ray data are consistent with a latitude distribution of the form I(b) = C sub 1 + C sub 2/sin b, with the second term being dominant. This result suggests that the radiation above 100 MeV is coming largely from local regions of the galactic disk. Between 35 and 100 MeV, a similar equation is also a good representation of the data, but here the two terms are comparable. These results indicate that the diffuse radiation above 35 MeV consists of two parts, one with a relatively hard galactic component and the other an isotropic, steep spectral component which extrapolates back well to the low energy diffuse radiation. The steepness of the diffuse isotropic component places significant constraints on possible theoretical models of this radiation.

Fichtel, C. E.↗
NASA NTRS1985

Higher order methods for convection-diffusion problems

This paper applies C1 cubic Hermite polynomials embedded in an orthogonal collocation scheme to the spatial discretization of the unsteady nonlinear Burgers equation as a model of the equations of fluid mechanics. The temporal discretization is carried out by means of either a noniterative finite difference or an iterative finite difference procedure. Results of this method are compared with those of a second-order finite difference scheme and a splined-cubic Taylor's series scheme. Stability limits are derived and the matrix structure of the several schemes are compared.

Murphy, J. D.↗
NASA NTRS1972

Diffusive and radiative effects on vaporization times of drops in film boiling

Diffusive and radiative effects are incorporated into an analysis for the vaporization time of drops in film boiling. The momentum, energy, and continuity equations are solved with some appropriate simplifications so as to obtain a simple closed form solution for the overall film boiling heat transfer coefficient. Next, a theoretical expression for the droplet vaporization time is developed and compared to the measured vaporization times of water droplets vaporizing into air, argon, nitrogen, and helium. The agreement between experiment and theory is good. Under the helium blanket, the diffusive evaporative component is significant in comparison to the film boiling component.

Baumeister, K. J.↗
NASA NTRS1973

An auroral breakup mechanism

A purely growing electrostatic drift instability driven by the electron temperature gradient at the inner edge of the plasma sheet can grow for large enough values of the temperature gradient. The parallel electric field associated with the instability is localized near the magnetic equator. The growth of the drift instability leads to enhanced whistler noise and increased electron pitch angle diffusion. If the current limit is exceeded in the ionosphere while the parallel electric field of the drift instability exists along the field line, rapid electron precipitation (the auroral breakup) can result.

Maggs, J. E.↗
NASA NTRS1974

U.S. National Congress of Applied Mechanics, 7th, University of Colorado, Boulder, Colo., June 3-7, 1974, Proceedings

A series of experimental, theoretical, and review papers dealing with various aspects of applied mechanics, including structural analysis, mathematical methods, and fluid dynamics. Some of the topics covered are: the direct formulation of two-dimensional mechanics problems, the mechanical aspects of athletics, Laplace's tidal equations, cochlear mechanics, dynamic constitutive equations, optical techniques for the determination of displacement and strain fields, chemical and biochemical problems involving nonlinear diffusion processes, thermal stresses, the rotating rod viscometer, finite-element methods in plasticity, ships of minimum resistance, and wide-band random vibration of structures. Individual items are announced in this issue.

Source record↗
NASA NTRS1975

Mathematical formulation of scatter-free propagation of solar cosmic rays

The observations of strong persistent velocity anisotropies in solar flare events demand a mathematical theory closer to the extreme of scatter-free (deterministic) propagation rather than diffusive (stochastic) transport, since the latter breaks down as inferred mean-free-paths exceed 0.1 AU. Equations are derived for the time-dependent phase-space density, and Laplace transform techniques are used to obtain solutions under rather general conditions. The case of an Archimedean spiral field has been solved numerically, and the results compared with observations from Mariner and Explorer spacecraft of nearly 0.4 MeV proton intensity and anisotropy histories. These can both be replicated if the inner boundary of the modulation region is placed beyond 2 AU.

Nolte, J. T.↗
NASA NTRS1976

High temperature reflectance of hyperpure slip cast silica

An experimental program to optically characterize hyperpure slip cast silica at elevated temperatures is described. This sintered ceramic material is being used in the development of a reflective heat shield for the Outer Planet Entry Probe. A unique high temperature integrating sphere reflectometer is described and reflectance data in the spectral range 0.23-2.3 micron are presented at temperatures up to 1427 C in air and helium environments. Kubelka-Munk (K-M) scattering and absorptance coefficients (K and S), using the Reichman extension of the K-M theory, are presented as is a unique method of obtaining K and S from the Reichman equations using only reflectance data. The hyperpure slip cast silica material was found to retain its high diffuse reflectance at temperatures up to 1204 C, with slight degradation occurring at higher temperatures.

Dillow, C. F.↗
NASA NTRS1978

The development of a three-dimensional partially elliptic flow computer program for combustor research

A three dimensional, partially elliptic, computer program was developed. Without requiring three dimensional computer storage locations for all flow variables, the partially elliptic program is capable of predicting three dimensional combustor flow fields with large downstream effects. The program requires only slight increase of computer storage over the parabolic flow program from which it was developed. A finite difference formulation for a three dimensional, fully elliptic, turbulent, reacting, flow field was derived. Because of the negligible diffusion effects in the main flow direction in a supersonic combustor, the set of finite-difference equations can be reduced to a partially elliptic form. Only the pressure field was governed by an elliptic equation and requires three dimensional storage; all other dependent variables are governed by parabolic equations. A numerical procedure which combines a marching integration scheme with an iterative scheme for solving the elliptic pressure was adopted.

Pan, Y. S.↗