Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 613 records · Page 34

Effect of molecular shape and flexibility on gamma-ray directional correlations.

Huntress' formalism (1968) is used to calculate the angular correlation behavior for an ensemble of radioactive nuclei, each of which is bound to an asymmetric-top molecule. The calculation is carried out both in an extreme narrowing limit (rapid motion) and in an adiabatic limit (slow motion). The two limiting cases span nearly the whole range of rates of rotational reorientation. It is shown that the problem of internal rotation is formally identical to that of rotation of a symmetric-top molecule. The experimentally observable time-integrated anisotropy is plotted for several rates of internal rotational diffusion.

Marshall, A. G.↗

Theoretical investigations of plasma processes

System analyses are presented for electrically sustained, collision dominated plasma centrifuges, in which the plasma rotates under the influence of the Lorentz forces resulting from the interaction of the current density fields with an external magnetic field. It is shown that gas discharge centrifuges are technically feasible in which the plasma rotates at speeds up to 1 million cm/sec. The associated centrifugal forces produce a significant spatial isotope separation, which is somewhat perturbed in the viscous boundary layers at the centrifuge walls. The isotope separation effect is the more pronounced. The induced magnetic fields have negligible influence on the plasma rotation if the Hall coefficient is small. In the technical realization of collision dominated plasma centrifuges, a trade-off has to be made between power density and speeds of rotation. The diffusion of sputtered atoms to system surfaces of ion propulsion systems and the deposition of the atoms are treated theoretically by means of a simple model which permits an analytical solution. The problem leads to an inhomogeneous integral equation.

Wilhelm, H. E.↗

Radio and EUV observations of a coronal hole

Current ideas and modeling assumptions concerning the structure of the transition region and corona are summarized. Observations of a specific coronal hole are reported which were made with the EUV spectroheliometer aboard Skylab and with high-resolution radio telescopes in Australia and West Germany. Theoretical EUV line intensities and radio brightness temperatures are calculated, and an attempt is made to analyze the observations on the basis of standard one-dimensional plane-parallel hydrostatic-equilibrium models for the transition region and corona. It is found that any given standard model predicts either higher radio brightnesses or lower EUV line intensities than are observed. Several potential problem areas are examined in order to resolve this conflict. It is concluded that a model involving departures from hydrostatic equilibrium, varying conductive flux, mass outflow, as well as effects of thermal diffusion and nonequilibrium ionization should be investigated in more detail.

Dulk, G. A.↗

Estimation of atmospheric species concentrations from remote sensing data

The atmospheric species concentration distribution-estimation problem that arises from the measurement configuration of remote sensing, which involves a field of view that moves across the region, is solved by means of the partial differential equation for atmospheric diffusion and the Wiener-Hopf theory. A numerical study is also undertaken for the estimation of concentration distribution downwind of a hypothetical, continuous ground-level source of pollutants.

Omatu, S.↗

Assessment of discretization schemes to reduce numerical diffusion in the calculation of complex flows

Comprehensive criteria for evaluating finite-difference schemes to reduce numerical diffusion in the solution of the Reynolds averaged Navier-Stokes equations have been established in this study. Several candidate schemes have been examined, based on these criteria. Taylor series analyses, model problem studies and examination of the difference equations have been used to assess the schemes. Two schemes, which showed the most promise, and satisfied a majority of the criteria, are recommended for incorporation into two dimensional and three dimensional computer codes for further evaluation.

Syed, S. A.↗

APART/PADE analytical evaluation of the diffuse infrared background experiment for NASA's cosmic background explorer

The stray light performance of NASA's Diffuse Infrared Background Experiment (DIRBE) has been calculated using the APART/PADE code. That code has been upgraded to handle off-axis optical systems such as DIRBE. Under observing conditions, sunlight is attenuated by 21 to 29 orders-of-magnitude; uniform diffuse illumination, by 6 to 9 orders of magnitude. The sunlight is attenuated by multiple diffraction at an external sun shield and by a forebaffle at the entrance to the DIRBE. Stray light is nearly a linear function of the Bidirectional Reflectance Distribution Function (BRDF) of the primary mirror at all wavelengths - or it can be made to be so by reducing the size of the field defining stop from 0.866 to 0.7 degrees square. For uniform diffuse illumination, half of the stray radiation comes from within 5 degrees of the center of the field-of-view. Particulate contamination of the primary mirror is expected to be a problem and special polishing and cleaning procedures are recommended - before mirror overcoating and again before flight.

Evans, D. C.↗

Coarsening in binary solid-liquid mixtures

A theory of Ostwald ripening has been developed for a solid-liquid mixture cosisting of a low volume fraction array of spherical solid particles in a liquid wherein the coarsening process proceeds via the transport of both heat and mass. It is found that the simultaneous transport of heat and mass during ripening does not alter the exponents of the temporal power laws governing the ripening process from their classical values but does alter the amplitudes of these power laws. The growth rate of the cube of the average particle radius, the rate constant, is found to depend both on the alloy solute concentration and the ratio of the thermal to solutal diffusivities. In most metallic systems, a large decrease in the rate constant can be expected with small additions of solute to a pure metal. Possible extensions of this theory to the analogous problem of ripening in isothermal ternary alloys are also discussed.

Voorhees, P. W.↗

Ignition and structure of a laminar diffusion flame in a compressible mixing layer with finite rate chemistry

The ignition and structure of a reacting compressible mixing layer is considered using finite rate chemistry lying between two streams of reactants with different freestream speeds and temperatures. Numerical integration of the governing equations show that the structure of the reacting flow can be quite complicated depending on the magnitude of the Zeldovich number. An analysis of both the ignition a diffusion flame regimes is presented using a combination of large Zeldovich number asymptotics and numerics. This allows to analyze the behavior of these regimes as a function of the parameters of the problem.

Grosch, C. E.↗

Ignition and structure of a laminar diffusion flame in a compressible mixing layer with finite rate chemistry

The ignition and structure of a reacting compressible mixing layer is considered using finite rate chemistry lying between two streams of reactants with different freestream speeds and temperatures. Numerical integration of the governing equations show that the structure of the reacting flow can be quite complicated depending on the magnitude of the Zeldovich number. An analysis of both the ignition and diffusion flame regimes is presented using a combination of large Zeldovich number asymptotics and numerics. This allows to analyze the behavior of these regimes as a function of the parameters of the problem.

Grosch, C. E.↗

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗

Internal low Reynolds number turbulent heat transfer

The results of a semi-analytical and experimental study of adiabatic tube flow and an analytical and experimental investigation of the thermal entrance region for gas flowing through electrically heated circular tubes are presented. Emphasis is placed on the low Reynolds number turbulent flow regime--defined as fully turbulent flow at bulk Reynolds numbers from 3,000 to about 15,000. Adiabatic air velocity and friction data and localized heat transfer measurements for air and helium, at low heating rates, are presented for this range. The adiabatic data were obtained in a 1.61 inch 1D tube for flow at bulk Reynolds numbers from 3,000 to 15,000. A continuous, Reynolds number-dependent, profile is developed from the data by using a modification of Reichardts' wall and middle law eddy diffusivity expressions. The velocity profile satisfies continuity. It is valid for all Reynolds numbers in excess of 3,000 for which the flow is fully turbulent and the Blasius friction factor expression is valid. The thermal entrance problem for a fully developed velocity profile is solved analytically by the method of Sparrow, Hallman, and Siegel. The solution is based on the profile developed from the velocity study. Tabular values of the eigenvalues and normalized Nusselt numbers for gases are presented for it range of Reynolds numbers from 3,000 to 50,000. The axial variation of Nusselt number is found to be correlated by [Nu/Nu_(∞)] = 1 + 0.8 (1+ 70,000 Re^(-3/2)) ((x/D)^(-1)) to within ± 5 per cent for x/D ≥2. The fully developed value agrees with the Dittus - Boelter correlation, Nu_(∞) = 0.021 Re^(0.8) Pr^(0.4) For the eigenvalues, λ^(2)_(n), and the associated constants, A_(n), correlations of the form λ^(2)_(n) =A_(1,n) Re^(-b_(1,n) + C_(1,n) Re^(-d_(1,n)) A_(n) = -A_(2,n) Re^(-b_(2,n) + C_(2,n) Re^(-d_(2,n)) are obtained. The coefficients and powers are presented in tabular form. Heat transfer data are presented, primarily for helium, for the low Reynolds number turbulent range. A one-quarter inch, resistively heated, vertical, circular tube was used for the study. The data cover an axial range from 1.2 to 96 diameters; wall-to-bulk temperature ratios vary from 1 to 1.4. In the low Reynolds number turbulent regime, these data clearly support the present analytical solution rather than the prediction obtained by applying the eddy diffusivity distribution used by Sparrow, Hallman, and Siegel.

Harold C. Reynolds, Jr.↗

Automatic phase control in solar power satellite systems

Various approaches to the problem of generating, maintaining and distributing a coherent, reference phase signal over a large area are suggested, mathematically modeled and analyzed with respect to their ability to minimize: phase build-up, beam diffusion and beam steering phase jitter, cable length, and maximize power transfer efficiency. In addition, phase control configurations are suggested which alleviate the need for layout symmetry.

Lindsey, W. C.↗

Mars - The role of the regolith in volatile storage and atmospheric evolution

Earth-based observations of the optical properties of Mars, laboratory alteration experiments, Mariner 9 data, and Viking lander observations are reviewed which indicate the existence of hydrated minerals on the Martian surface. The problem of estimating the total inventory of H2O, CO2, and other volatiles that may be stored in the Martian regolith is discussed, and morphological evidence is cited which supports the possible occurrence of hard-frozen permafrost in some areas on Mars. Physical adsorption of CO2 in the Martian regolith is considered, buffering of atmospheric pressure by regolith-atmosphere reequilibration is examined, and the problem of the exceedingly low inventory of Ar-36 in the Martian atmosphere is investigated. The possibility is considered that earth's degassing history may have had a larger component of long-diffusion-path degassing from a large hot body while the Martian surface inventory may have received a proportionately greater contribution from volatization of accreting particles.

Fanale, F. P.↗

Algorithms for the Fractional Calculus: A Selection of Numerical Methods

Many recently developed models in areas like viscoelasticity, electrochemistry, diffusion processes, etc. are formulated in terms of derivatives (and integrals) of fractional (non-integer) order. In this paper we present a collection of numerical algorithms for the solution of the various problems arising in this context. We believe that this will give the engineer the necessary tools required to work with fractional models in an efficient way.

Diethelm, K.↗

Microgravity Studies of Liquid-Liquid Phase Transitions in Alumina-Yttria Melts

The scientific objective of this research is to increase the fundamental knowledge base for liquid- phase processing of technologically important oxide materials. The experimental objective is to define conditions and hardware requirements for microgravity flight experiments to test and expand the experimental hypotheses that: 1. Liquid phase transitions can occur in undercooled melts by a diffusionless process. 2. Onset of the liquid phase transition is accompanied by a large change in the temperature dependence of melt viscosity. Experiments on undercooled YAG (Y3A15012)- and rare earth oxide aluminate composition liquids demonstrated a large departure from an Arrhenian temperature dependence of viscosity. Liquid YAG is nearly inviscid at its 2240 K melting point. Glass fibers were pulled from melts undercooled by ca. 600 K indicating that the viscosity is on the order of 100 Pans (1000 Poise) at 1600 K. This value of viscosity is 500 times greater than that obtained by extrapolation of data for temperatures above the melting point of YAG. These results show that the liquids are extremely fragile and that the onset of the highly non-Arrhenian viscosity-temperature relationship occurs at a temperature considerably below the equilibrium melting point of the solid phases. Further results on undercooled alumina-yttria melts containing 23-42 mole % yttrium oxide indicate that a congruent liquid-liquid phase transition occurs in the undercooled liquids. The rates of transition are inconsistent with a diffusion-limited process. This research is directed to investigation of the scientifically interesting phenomena of polyamorphism and fragility in undercooled rare earth oxide aluminum oxide liquids. The results bear on the technologically important problem of producing high value rare earth-based optical materials.

Guynes, Buddy↗

Emittance Measurements for a Thin Liquid Sheet Flow

The Liquid Sheet Radiator (LSR) is an external flow radiator that uses a triangular-shaped flowing liquid sheet as the radiating surface. It has potentially much lower mass than solid wall radiators such as pumped loop and heat pipe radiators, along with being nearly immune to micrometeoroid penetration. The LSR has an added advantage of simplicity. Surface tension causes a thin (100-300 microns) liquid sheet to coalesce to a point, causing the sheet flow to have a triangular shape. Such a triangular sheet is desirable since it allows for simple collection of the flow at a single point. A major problem for all external flow radiators is the requirement that the working fluid be of very low (approx. 10(sup -8) torr) vapor pressure to keep evaporative losses low. As a result, working fluids are limited to certain oils (such as used in diffusion pumps) for low temperatures (300-400 K) and liquid metals for higher temperatures. Previous research on the LSR has been directed at understanding the fluid mechanics of thin sheet flows and assessing the stability of such flows, especially with regard to the formation of holes in the sheet. Taylor studied extensively the stability of thin liquid sheets both theoretically and experimentally. He showed that thin sheets in a vacuum are stable. The latest research has been directed at determining the emittance of thin sheet flows. The emittance was calculated from spectral transmittance data for the Dow Corning 705 silicone oil. By experimentally setting up a sheet flow, the emittance was also determined as a function of measurable quantities, most importantly, the temperature drop between the top of the sheet and the temperature at the coalescence point of the sheet. Temperature fluctuations upstream of the liquid sheet were a potential problem in the analysis and were investigated.

Englehart, Amy N.↗

Numerical modeling of physical vapor transport under microgravity conditions: Effect of thermal creep and stress

One of the most promising applications of microgravity (micro-g) environments is the manufacture of exotic and high-quality crystals in closed cylindrical ampoules using physical vapor transport (PVT) processes. The quality enhancements are believed to be due to the absence of buoyant convection in the weightless environment - resulting in diffusion-limited transport of the vapor. In a typical experiment, solid-phase sample material is initially contained at one end of the ampoule. The sample is made to sublime into the vapor phase and deposit onto the opposite end by maintaining the source at an elevated temperature with respect to the deposit. Identification of the physical factors governing both the rates and uniformity of crystal growth, and the optimization of the micro-g technology, will require an accurate modeling of the vapor transport within the ampoule. Previous micro-g modeling efforts have approached the problem from a 'classical' convective/diffusion formulation, in which convection is driven by the action of buoyancy on thermal and solutal density differences. The general conclusion of these works have been that in low gravity environments the effect of buoyancy on vapor transport is negligible, and vapor transport occurs in a diffusion-limited mode. However, it has been recently recognized than in the non-isothermal (and often low total pressure) conditions encountered in ampoules, the commonly-assumed no-slip boundary condition to the differential equations governing fluid motion can be grossly unrepresentative of the actual situation. Specifically, the temperature gradients can give rise to thermal creep flows at the ampoule side walls. In addition, temperature gradients in the vapor itself can, through the action of thermal stress, lead to bulk fluid convection.

Mackowski, Daniel W.↗

An inverse problem design method for branched and unbranched axially symmetrical ducts

This paper concerns the potential flow design of axially symmetrical ducts of both circular and annular cross section with or without wall suction or blowing slots. The objective of the work was to develop a method by which such ducts could be designed with directly prescribed wall pressure variation. Previous axially symmetrical design methods applied only to circular cross sectional ducts and required that the pressure distribution be prescribed along the duct centerline and not along the duct wall. The present method uses an inverse problem approach which extends the method of Stanitz to the axially symmetrical case, and an approximation is used to account for the stagnation point in branched duct designs. Two examples of successful designs of diffusers with suction slots are presented.

Nelson, C. D.↗