Search NASA⌕ Search

SEARCH · Search NASA

Results for “Diffusion problems”

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 433 records · Page 24

Verification of RadCAD: Specular Capabilities

As part of the RadCAD's development process, it is necessary to compare RadCAD's results with other radiation tools and exact solutions when and where possible. Form factor algorithms have been previously verified with exact solutions. This paper will consider RadCAD's specular capabilities. First, radiation exchange factors will be compared against exact solutions and results from TRASYS for various geometries. Critical dimensions and optical properties are changed for each geometry. Second, a specular adjunct plate system will be used to verify absorbed heat fluxes. This particular geometric problem has had some attention in the literature. Previous authors have used this problem to validate software results with exact analytical solution. This paper will compare absorbed heat rates against the exact solution and other published results from other thermal radiation tools. The agreement between RadCAD and the exact solutions is good. The maximum error for both specular and diffuse exchange factors for both geometries and all optical properties was 3%. The absorbed fluxes differed by a maximum of 4% for the adjunct plate problem.

Lucas, Sam↗

Hydrogen Effects on GaAs, Status and Progress

It has been observed that GaAs devices in hermetically sealed packages when exposed to temperatures as low as 125oC can exhibit unacceptable degradation in both RF and DC characteristics. The source of the problem was found to be hydrogen that has been absorbed in the package's metals (Kovar, plating, etc.) and converted into atomic hydrogen within the Pt metallization of the gate structure. Subsequently, atomic hydrogen diffuses into the channel region of the FET structure and neutralizes the Si donors, resulting in a degradation of the device characteristics. Most of the data indicates the onslaught of the problem to occur after 500 hours at 125oC. This onslaught has been observed to be dependent on the thickness and processing conditions of the passivation layer, the sealing environment, and the amount of Pt or Pd in the gate structure of the device.

Hydrogen GaAs↗

The Origin and Time Dependence of the Amount and Composition of Non-Constituent Gases Present in Crystal Growth Systems

Presence of different, non-constituent gases may be a critical factor in crystal growth systems. In Physical Vapor Transport processes the cras(es) can be used intentionally (to prevent excessively high, unstable growth conditions), or can evolve unintentionally during the course of the process (which may lead to undesired reduction in the -rowth rate). In melt growth, particularly under low gravity conditions (reduced hydrostatic pressure) the gas present in the system may contribute to formation of voids in the growing crystals and even to a separation of the crystal and the liquid phase [1]. On the other hand, some amount of gas may facilitate 'contactless' crystal growth particularly under reduced gravity conditions [2 - 6]. Different non-constituent gases may be present in growth ampoules, and their amount and composition may change during the crystallization process. Some gases can appear even in empty ampoules sealed originally under high vacuum: they may diffuse in from the outside, and/or desorb from the ampoule walls. Residual gases can also be generated by the source materials: even very high purity commercial elements and compounds may contain trace amounts of impurities, particularly oxides. The oxides may have low volatilities themselves but their reaction with other species, particularly carbon and hydrogen, may produce volatile compounds like water or carbon oxides. The non-constituent gases, either added initially to the system or evolved during the material processing, may diffuse out of the ampoule during the course of the experiment. Gases present outside (e.g. as a protective atmosphere or thermal conductor) may diffuse into the ampoule. In either case the growth conditions and the quality of the crystals may be affected. The problem is of a particular importance in sealed systems where the amount of the gases cannot be directly controlled. Therefore a reasonable knowledge and understanding of the origin, composition, magnitude, and change with time of gases present in sealed ampoules may be important for a meaningful control and interpretation of crystal growth processes. This problem is of a particular importance for processing of electronic materials in space because (i) safety considerations require using sealed systems only, and (ii) high cost of crystal growth experiments in microgravity calls for a throughout, accurate description of the processing conditions necessary for a meaningful, efficient, and conclusive interpretation of the space results. In this paper we present the results of our extensive studies on gases in closed crystal growth systems which include: (a) Degassing properties of fused silica; (b) Generation of inert gases by source materials (CdTe, ZnTe, CdZnTe, ZnSe, PbTe, PbSe, PbSeTe); (c) Diffusive cas losses from silica glass ampoules.

Palosz, Witold↗

Vibration characteristics of the HPOTP (High-Pressure Oxygen Turbopump) of the SSME (Space Shuttle Main Engine)

Attention is given to rotor dynamic problems that have been encountered and eliminated in the course of Space Shuttle Main Engine (SSME) development, as well as continuing, subsynchronous problems which are being encountered in the development of a 109-percent power level engine. The basic model for the SSME's High Pressure Oxygen Turbopump (HPOTP) encompasses a structural dynamic model for the rotor and housing, and component models for the liquid and gas seals, turbine clearance excitation forces, and impeller diffuser forces. Linear model results are used to examine the synchronous response and stability characteristics of the HPOTP, with attention to bearing load and stability problems associated with the second critical speed. Differences between linear and nonlinear model results are discussed and explained in terms of simple models. Simulation results indicate that while synchronous bearing loads can be reduced, subsynchronous motion is not eliminated by seal modifications.

Childs, D. W.↗

Nonlinear theory of electron neutralization waves in ions beams with dissipation

An analytical theory of nonlinear neutralization waves generated by injection of electrons from a grid in the direction of a homogeneous ion beam of uniform velocity and infinite extension is presented. The electrons are assumed to interact with the ions through the self-consistent space charge field and by strong collective interactions, while diffusion in the pressure gradient is disregarded (zero-temperature approximation). The associated nonlinear boundary-value problem is solved in closed form by means of a von Mises transformation. It is shown that the electron gas moves into the ion space in the form of a discontinuous neutralization wave, which exhibits a periodic field structure (incomplete neutralization). This periodic wave structure is damped out by intercomponent momentum transfer - i.e., after a few relaxation lengths a quasi-neutral plasma results.

Wilhelm, H. E.↗

Cosmic-ray transport theory and out-of-the-ecliptic exploration

The reasons for studying cosmic-ray transport theory are summarized and the fundamentally three-dimensional nature of the process is pointed out. It is shown that observations in the solar ecliptic plane cannot unambiguously test transport theories since the solutions to the transport equations depend critically on boundary conditions and variation of parameters such as diffusion tensor out of the ecliptic. Sample calculations (Fokker-Planck coefficient) are shown which illustrate the problem. It is concluded that out-of-the-ecliptic observations are essential to further test transport theory.

Jokipii, J. R.↗

Evolution of weak disturbances in inert binary mixtures

The evolution of weak disturbances in inert binary mixtures is determined for the one-dimensional piston problem. The interaction of the dissipative and nonlinear mechanisms is described by Burgers' equation. The binary mixture diffusion mechanisms enter as an additive term in an effective diffusivity. Results for the impulsive motion of a piston moving into an ambient medium and the sinusoidally oscillating piston are used to illustrate the results and elucidate the incorrect behavior pertaining to the associated linear theory.

Rasmussen, M. L.↗

Development of low cost contacts to silicon solar cells

The results of the second phase of the program of developing low cost contacts to silicon solar cells using copper are presented. Phase 1 yielded the development of a plated Pd-Cr-Cu contact system. This process produced cells with shunting problems when they were heated to 400 C for 5 minutes. Means of stopping the identified copper diffusion which caused the shunting were investigated. A contact heat treatment study was conducted with Pd-Ag, Ci-Ag, Pd-Cu, Cu-Cr, and Ci-Ni-Cu. Nickel is shown to be an effective diffusion barrier to copper.

Tanner, D. P.↗

A theoretical study of the effects of interacting grain boundaries on electron-beam-induced currents

The influence of several grain boundaries on the electron-beam-induced current in polycrystalline silicon is investigated. The simple geometry contemplated here, consisting of two plane parallel grain boundaries possessing different effective, albeit constant, surface recombination velocities, intersecting a planar junction at right angles allows for an exact solution of the corresponding boundary value problem. It is shown that if the distance between the grain boundaries w becomes less than the diffusion length for minority carriers L within the grain, the grain boundary with the smaller effective recombination velocity becomes less and less visible the smaller the ratio w/L.

Luke, K. L.↗

Solvability condition for needle crystals at large undercooling in a nonlocal model of solidification

It is explicitly shown that, in a realistic model of diffusion-controlled dendritic solidification, Ivantsov's continuous family of steady-state needle crystals is destroyed by the addition of surface tension. The starting point is in the exact integro-differential equation for the one-sided model, in two dimensions, in a moving frame of reference. In the limit of large undercooling, where the range of the diffusion field is much smaller than the radius of curvature of the tip of the needle, this problem is reduced to a linear, inhomogeneous differential equation of infinite order. A solvability condition for this equation is derived and it is shown that solutions cease to exist for arbitrarily small but finite isotropic surface tension.

Caroli, B.↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. The Roe approximate Riemann solution scheme or the computationally less expensive Advection Upstream Splitting Method (AUSM) flux-splitting scheme are used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passage and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

An Eulerian/Lagrangian coupling procedure for three-dimensional vortical flows

A coupled Eulerian/Lagrangian method is presented for the reduction of numerical diffusion observed in solutions of 3D vortical flows using standard Eulerian finite-volume time-marching procedures. A Lagrangian particle tracking method, added to the Eulerian time-marching procedure, provides a correction of the Eulerian solution. In turn, the Eulerian solution is used to integrate the Lagrangian state-vector along the particles trajectories. While the Eulerian solution ensures the conservation of mass and sets the pressure field, the particle markers describe accurately the convection properties and enhance the vorticity and entropy capturing capabilities of the Eulerian solver. The Eulerian/Lagrangian coupling strategies are discussed and the combined scheme is tested on a constant stagnation pressure flow in a 90 deg bend and on a swirling pipe flow. As the numerical diffusion is reduced when using the Lagrangian correction, a vorticity gradient augmentation is identified as a basic problem of this inviscid calculation.

Felici, Helene M.↗

Two-dimensional CFD modeling of wave rotor flow dynamics

A two-dimensional Navier-Stokes solver developed for detailed study of wave rotor flow dynamics is described. The CFD model is helping characterize important loss mechanisms within the wave rotor. The wave rotor stationary ports and the moving rotor passages are resolved on multiple computational grid blocks. The finite-volume form of the thin-layer Navier-Stokes equations with laminar viscosity are integrated in time using a four-stage Runge-Kutta scheme. Roe's approximate Riemann solution scheme or the computationally less expensive advection upstream splitting method (AUSM) flux-splitting scheme is used to effect upwind-differencing of the inviscid flux terms, using cell interface primitive variables set by MUSCL-type interpolation. The diffusion terms are central-differenced. The solver is validated using a steady shock/laminar boundary layer interaction problem and an unsteady, inviscid wave rotor passage gradual opening problem. A model inlet port/passage charging problem is simulated and key features of the unsteady wave rotor flow field are identified. Lastly, the medium pressure inlet port and high pressure outlet port portion of the NASA Lewis Research Center experimental divider cycle is simulated and computed results are compared with experimental measurements. The model accurately predicts the wave timing within the rotor passages and the distribution of flow variables in the stationary inlet port region.

Welch, Gerard E.↗

Obliquity variation in a Mars climate evolution model

The existence of layered terrain in both polar regions of Mars is strong evidence supporting a cyclic variation in climate. It has been suggested that periods of net deposition have alternated with periods of net erosion in creating the layered structure that is seen today. The cause for this cyclic climatic behavior is variation in the annually averaged latitudinal distribution of solar insolation in response to obliquity cycles. For Mars, obliquity variation leads to major climatological excursion due to the condensation and sublimation of the major atmospheric constituent, CO2. The atmosphere will collapse into the polar caps, or existing caps will rapidly sublimate into the atmosphere, dependent upon the polar surface heat balance and the direction of the change in obliquity. It has been argued that variations in the obliquity of Mars cause substantial departures from the current climatological values of the surface pressure and the amount of CO2 stored in both the planetary regolith and polar caps. In this new work we have modified the Haberle et al. model to incorporate variable obliquity by allowing the polar and equatorial insolation to become functions of obliquity, which we assume to vary sinusoidally in time. As obliquity varies in the model, there can be discontinuities in the time evolution of the model equilibrium values for surface pressure, regolith, and polar cap storage. The time constant, tau r, for the regolith to find equilibrium with the climate is estimated--depending on the depth, thermal conductivity, and porosity of the regolith--between 10(exp 4) and 10(exp 6) yr. Thus, using 2000-yr timesteps to move smoothly through the 0.1250 m.y. obliquity cycles, we have an atmosphere/regolith system that cannot be assumed in equilibrium. We have dealt with this problem by limiting the rate at which CO2, can move between the atmosphere and regolith, mimicking the diffusive nature and effects of the temperature and pressure waves, by setting the time rate of change of regolith storage proportional to the difference between equilibrium storage and current storage.

Tyler, D.↗

Filtration Combustion in Smoldering and SHS

Smolder waves and SHS (self-propagating high-temperature synthesis) waves are both examples of filtration combustion waves propagating in porous media. Smoldering combustion is important for the study of fire safety. Smoldering itself can cause damage, its products are toxic and it can also lead to the more dangerous gas phase combustion which corresponds to faster propagation at higher temperatures. In SHS , a porous solid sample, consisting of a finely ground powder mixture of reactants, is ignited at one end. A high temperature thermal wave, having a frontal structure, then propagates through the sample converting reactants to products. The SHS technology appears to enjoy a number of advantages over the conventional technology, in which the sample is placed in a furnace and "baked" until it is "well done". The advantages include shorter synthesis times, greater economy, in that the internal energy of the reactions is employed rather than the costly external energy of the furnace, purer products, simpler equipment and no intrinsic limitation on the size of the sample to be synthesized as exists in the conventional technology. When delivery of reactants through the pores to the reaction site is an important aspect of the combustion process, it is referred to as filtration combustion. The two types of filtration combustion have a similar mathematical formulation, describing the ignition, propagation and extinction of combustion waves in porous media. The goal in each case, however, is different. In smoldering the desired goal is to prevent propagation, whereas in SHS the goal is to ensure propagation of the combustion wave, leading to the synthesis of desired products. In addition, the scales in the two areas of application differ. Smoldering generally occurs at lower temperatures and propagation velocities than in SHS nevertheless, the two applications have much in common so that what is learned fit make application can be used to advantage in the other. In porous media, melting often occurs ahead of the propagating combustion wave. In certain cases there is so much melting that the porous solid structure is destroyed, e.g., by melting and a suspension arises, consisting of a liquid bath containing solid particles and/or gas bubbles. The resulting combustion wave is referred to as a liquid flame. We have considered a number of problems involving filtration combustion. Here, we describe four such studies: (A) rapid buoyant filtration combustion waves; (B) diffusion driven combustion waves; (C) rapidly propagating liquid flames in gravitational fields; and (D) gas-phase influence on liquid flames in gravitational fields.

Matkowsky, Bernard J.↗

L Band Brightness Temperature from Forest: Comparison of Approximate Techniques

In this paper, three approximate physical microwave radiometry models have been used to calculate brightness temperatures from a forest canopy at L-band. These models are (1) tau-omega model (zero order scattering approximation to radiative transfer equations), (2) successive order of scattering model up to first order (first order scattering approximation to the radiative transfer equations), and (3) Peak technique utilizing the active solution obtained from the Distorted Born Approximation (DBA). These models are physically-based and treat vegetation as a layer of discrete scatterers over a rough surface. Vegetation components within the canopy are represented by canonical shapes such as dielectric discs and cylinders. The tau-omega model is based on a zero-order solution to the radiative transfer (RT) equations. The model ignores scattering except for the effect of the scatterers in the attenuation of the emission through the vegetation. Application of the tau-omega model to data acquired during airborne and ground-based campaigns over the years has solidified scientific understanding of microwave interactions with different landscapes. In particular, shrubland, grasslands, agricultural crops, and light to moderate vegetation have been investigated. Its applicability to areas with a significant tree fraction is unknown. The first order scattering model is based on an iterative solution of the RT equation up to the first order. The first order solution is obtained by substituting the zeroth-order solution into the scattering source term and then solving the resulting radiative transfer equations. This formulation adds a new scattering term to the tau-omega model. It represents emission by particles in the layer and emission by the ground that is scattered once by particles in the layer. The resulting model represents an improvement over the standard zero-order solution (the tau-omega model) since it accounts for the scattered vegetation and ground radiation that can have a pronounced effect on the observed brightness temperature. The third model is based on the Peake formulation in conjunction with the DBA. The procedure for calculation of forest emission is accomplished by first calculating the bistatic scattering cross section for each type of scatterer, then by using the DBA to calculate specular albedo of the ground and the diffused albedo of the layer. Once the albedos are determined, Peake s principle relating active and passive problems can be used to determine the effective emissivity of the forest layer.

Kurum, Mehmet↗

VIIRS On-Orbit Optical Anomaly - Investigation, Analysis, Root Cause Determination and Lessons Learned

A gradual, but persistent, decrease in the optical throughput was detected during the early commissioning phase for the Suomi National Polar-Orbiting Partnership (SNPP) Visible Infrared Imager Radiometer Suite (VIIRS) Near Infrared (NIR) bands. Its initial rate and unknown cause were coincidently coupled with a decrease in sensitivity in the same spectral wavelength of the Solar Diffuser Stability Monitor (SDSM) raising concerns about contamination or the possibility of a system-level satellite problem. An anomaly team was formed to investigate and provide recommendations before commissioning could resume. With few hard facts in hand, there was much speculation about possible causes and consequences of the degradation. Two different causes were determined as will be explained in this paper. This paper will describe the build and test history of VIIRS, why there were no indicators, even with hindsight, of an on-orbit problem, the appearance of the on-orbit anomaly, the initial work attempting to understand and determine the cause, the discovery of the root cause and what Test-As-You-Fly (TAYF) activities, can be done in the future to greatly reduce the likelihood of similar optical anomalies. These TAYF activities are captured in the lessons learned section of this paper.

Iona, Glenn↗

Survey of the Relationship Between Theory and Experiment for Convective Heat Transfer From Rocket Combustion Gases

The problem of predicting the distribution of heat flux in a rocket engine is complicated by a flow field and thermodynamics so perturbed by as yet indescribable combustion effects that it continues to be necessary to utilize a simplified theoretical model against which to compare experimental results in order to gain insight. The most successful model applied to this problem is that of a compressible, turbulent, accelerating boundary layer flow with an arbitrary point of origin, and chemical reaction restricted to equilibrium recombination or diffusion-limited recombination on a catalytic wall. Analysis of this model was not possible without making several assumptions about the nature of the skin friction coefficient relationship and Reynold's analogy in an accelerating boundary layer flow. The paper summarizes the analysis of this model (together with detailed treatment of the special problems of variable properties, chemical reaction, and sources of transport properties), and describes several closed-form equations that have also been developed and applied to the problem. Described are results of several experiments with hot air flow made for the purpose of checking the validity of the assumptions in the boundary layer analysis. Most significant, however, are presentations of sample distributions of local heat flux measured in rocket engines operating with a wide range of conditions, propellants, injectors and sizes. When compared with predictions, these data show the significant range of deviation in heat flux that results from the combustion effects and which precludes a generally accurate prediction theory or correlation equation until the combustion process itself is correlated.

D R Bartz↗