Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion coefficients”

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,243 records · Page 69

A Tomographic Method for the Reconstruction of Local Probability Density Functions

A method of obtaining the probability density function (PDF) of local properties from path integrated measurements is described. The approach uses a discrete probability function (DPF) method to infer the PDF of the local extinction coefficient from measurements of the PDFs of the path integrated transmittance. The local PDFs obtained using the method are compared with those obtained from direct intrusive measurements in propylene/air and ethylene/air diffusion flames. The results of this comparison are good.

Sivathanu, Y. R.↗

Finite Volume Methods: Foundation and Analysis

Finite volume methods are a class of discretization schemes that have proven highly successful in approximating the solution of a wide variety of conservation law systems. They are extensively used in fluid mechanics, porous media flow, meteorology, electromagnetics, models of biological processes, semi-conductor device simulation and many other engineering areas governed by conservative systems that can be written in integral control volume form. This article reviews elements of the foundation and analysis of modern finite volume methods. The primary advantages of these methods are numerical robustness through the obtention of discrete maximum (minimum) principles, applicability on very general unstructured meshes, and the intrinsic local conservation properties of the resulting schemes. Throughout this article, specific attention is given to scalar nonlinear hyperbolic conservation laws and the development of high order accurate schemes for discretizing them. A key tool in the design and analysis of finite volume schemes suitable for non-oscillatory discontinuity capturing is discrete maximum principle analysis. A number of building blocks used in the development of numerical schemes possessing local discrete maximum principles are reviewed in one and several space dimensions, e.g. monotone fluxes, E-fluxes, TVD discretization, non-oscillatory reconstruction, slope limiters, positive coefficient schemes, etc. When available, theoretical results concerning a priori and a posteriori error estimates are given. Further advanced topics are then considered such as high order time integration, discretization of diffusion terms and the extension to systems of nonlinear conservation laws.

Barth, Timothy↗

Scattering of light in a spherical multilayer atmosphere

The problem of radiation diffusion in the atmosphere of a planet illuminated by solar rays was examined. The basic equations of the problem were obtained and the solution was found for the case when the absorption coefficient in the atmosphere is constant. The assumption was introduced of exponential variation of absorptivity with altitude, and a first approximation solution was proposed.

Titarchuk, L. G.↗

Nonlinear diffusion in the solid phase

A Fourier series expansion method has been used to develop simple and relatively accurate solutions for nonlinear Fickian diffusion in the solid phase. For the fundamental harmonic, the first term in the Fourier expansion is qualitatively similar but quantitatively quite different from the conventional exponential function. In contrast, the Fourier coefficients of the higher order harmonics are combinations of fast-decay and slow-decay components. The simple analytical results agree well with numerical solutions for one sample calculation.

Tsang, Tung↗

Computational Fluid Dynamics (CFD) Analysis for the Reduction of Impeller Discharge Flow Distortion

The use of Computational Fluid Dynamics (CFD) in the design and analysis of high performance rocket engine pumps has increased in recent years. This increase has been aided by the activities of the Marshall Space Flight Center (MSFC) Pump Stage Technology Team (PSTT). The team's goals include assessing the accuracy and efficiency of several methodologies and then applying the appropriate methodology(s) to understand and improve the flow inside a pump. The PSTT's objectives, team membership, and past activities are discussed in Garcia1 and Garcia2. The PSTT is one of three teams that form the NASA/MSFC CFD Consortium for Applications in Propulsion Technology (McConnaughey3). The PSTT first applied CFD in the design of the baseline consortium impeller. This impeller was designed for the Space Transportation Main Engine's (STME) fuel turbopump. The STME fuel pump was designed with three impeller stages because a two-stage design was deemed to pose a high developmental risk. The PSTT used CFD to design an impeller whose performance allowed for a two-stage STME fuel pump design. The availability of this design would have lead to a reduction in parts, weight, and cost had the STME reached production. One sample of the baseline consortium impeller was manufactured and tested in a water rig. The test data showed that the impeller performance was as predicted and that a two-stage design for the STME fuel pump was possible with minimal risk. The test data also verified another CFD predicted characteristic of the design that was not desirable. The classical 'jet-wake' pattern at the impeller discharge was strengthened by two aspects of the design: by the high head coefficient necessary for the required pressure rise and by the relatively few impeller exit blades, 12, necessary to reduce manufacturing cost. This 'jet-wake pattern produces an unsteady loading on the diffuser vanes and has, in past rocket engine programs, lead to diffuser structural failure. In industrial applications, this problem is typically avoided by increasing the space between the impeller and the diffuser to allow the dissipation of this pattern and, hence, the reduction of diffuser vane unsteady loading. This approach leads to small performance losses and, more importantly in rocket engine applications, to significant increases in the pump's size and weight. This latter consideration typically makes this approach unacceptable in high performance rocket engines.

Garcia, R.↗

A dynamically based transport parameterization for one-dimensional photochemical models of the stratosphere

The net vertical tracer flux in the stratosphere is due primarily to advection by the global-scale mean meridional circulation, not to diffusion by turbulent eddies. Using a simple model of this circulation, it is shown that the net flux can be approximated by a flux-gradient relationship in which the vertical 'eddy transport' coefficient is proportional to the square of the diabatic heating rate and inversely proportional to the rate at which meridional tracer gradients are destroyed by horizontal eddy mixing and chemical damping. Thus the transport coefficient will differ for tracers with differing chemical lifetimes. Profiles computed for various source gases (CH4 N2O, CF2Cl2, and CFCl3) show that this parameterization provides an improvement over conventional formulations that utilize a single transport coefficient for all species.

Holton, J. R.↗

Column Optical Depth (COD) Derived from CALIOP Ocean Surface Returns

The Cloud Aerosol Lidar with Orthogonal Polarization (CALIOP) on board the Cloud Aerosol Lidar Infrared Pathfinder Satellite Observations (CALIPSO) satellite has been making near-global measurements of clouds and aerosols since mid-June 2006. Among the properties reported in CALIPSO data products are estimates of total column optical depth (COD) obtained by integrating the retrieved extinction coefficients of detected particulates (i.e., aerosols and clouds) from 36 km to the surface. However, due to algorithm detection limits and instrument sensitivities, particulates can go undetected when the particulate loading is especially diffuse or in regions where the signal is attenuated by over-lying layers. Because extinction is not calculated where features are not detected, these undetected particulates can introduce low biases into the reported COD. Lidar ratio assumptions used in the extinction retrievals can also introduce errors. To minimize these biases, future releases of the CALIPSO lidar data products will implement an ocean-derived column optical depth (ODCOD) retrieval. The following paper briefly describes the algorithm, the current status of the algorithm development and verification efforts, and preliminary comparisons to col-located COD estimates derived using other retrieval schemes and obtained from other sensors.

Robert A Ryan↗

Influence of oscillation-induced diffusion on heat transfer in a uniformly heated channel

An analysis of the effect of flow oscillations on laminar flow heat transfer in a channel with uniform heat addition is presented. It is shown that the effect of flow oscillations will be to reduce the channel heat transfer coefficient. This effect is due to the fact that the heat addition along the channel wall produces an increasing fluid temperature along the channel length. The flow oscillations interacting with this positive temperature gradient will induce a heat flow back toward the channel inlet. This will tend to inhibit the heat transfer process and will raise the wall temperature required to transfer away a given amount of heat at the channel wall.

Siegel, R.↗

Use of a corrugated surface to enhance radiation tolerance in a GaAs solar cell

The use of a corrugated surface on a GaAs solar cell and its effects on radiation resistance were studied. A compute code was developed to determine the performance of the cell for various geometric parameters. The large optical absorption coefficient of GaAs allows grooves to be only 4-5 micrometers deep. Using accepted material parameters for GaAs solar cells the theoretical performances were compared for various corrugated cells before and after minority carrier diffusion length degradation. The total power output was maximized for both n(+)/p and p(+)/n cells. Optimum values of 1.0-1.5 and 5.0 micrometers for groove and ridge widths respectively were determined.

Leon, Rosa P.↗

Heterodyne efficiency for a coherent laser radar with diffuse or aerosol targets

The performance of a Coherent Laser Radar is determined by the statistics of the coherent Doppler signal. The heterodyne efficiency is an excellent indication of performance because it is an absolute measure of beam alignment and is independent of the transmitter power, the target backscatter coefficient, the atmospheric attenuation, and the detector quantum efficiency and gain. The theoretical calculation of heterodyne efficiency for an optimal monostatic lidar with a circular aperture and Gaussian transmit laser is presented including beam misalignment in the far-field and near-field regimes. The statistical behavior of estimates of the heterodyne efficiency using a calibration hard target are considered. For space based applications, a biased estimate of heterodyne efficiency is proposed that removes the variability due to the random surface return but retains the sensitivity to misalignment. Physical insight is provided by simulation of the fields on the detector surface. The required detector calibration is also discussed.

Frehlich, R. G.↗

Internal temperature gradient of alloy semiconductor melts from interrupted growth experiments

Interrupted growth experiments on Hg(0.8)Cd(0.2)Se alloys were performed in an attempt to estimate the actual internal temperature gradient in the melt during directional solidification. The results have been analyzed in terms of a theoretical model which assumes an axial composition profile characteristic of diffusion controlled solute redistribution during growth. A comparison of the calculated and measured values suggests a reduction in the applied (empty furnace) temperature gradient by about a factor of three. Values of the interface segregation coefficient (k) determined from the phase diagram, predict values of solute concentration build up in the solid during growth interruptions which are inconsistent with measured results. It appears that the k values used tend to underestimate the actual k values for higher alloy compositions and overestimate k for the lower alloy compositions.

Andrews, R. N.↗

3-D Modeling of Directional Solidification of a Non-Dilute Alloy with Temperature and Concentration Fields Coupling via Materials Properties Dependence and via Double Diffusive Convection

Numerical simulation of the HgCdTe growth by the vertical Bridgman method was performed using FIDAP finite element code. Double-diffusive melt convection is analyzed, as the primary factor at controls inhomogeneity of the solidified material. Temperature and concentration fields in the model are also coupled via material properties, such as thermal and solutal expansion coefficients with the dependence on both temperature and concentration, and melting temperature evaluation from pseudobinary CdTe-HgTe phase diagram. Experimental measurements were used to obtain temperature boundary conditions. Parametric study of the melt convection dependence on the gravity conditions was undertaken. It was found, that the maximum convection velocity in the melt can be reduced under certain conditions. Optimal conditions to obtain a near flat solidified interface are discussed. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D calculations are compared with previous 2- D findings. A video film featuring 3-D melt convection will be presented.

Bune, Andris V.↗

Hall-Effect Thruster Simulations with 2-D Electron Transport and Hydrodynamic Ions

A computational approach that has been used extensively in the last two decades for Hall thruster simulations is to solve a diffusion equation and energy conservation law for the electrons in a direction that is perpendicular to the magnetic field, and use discrete-particle methods for the heavy species. This "hybrid" approach has allowed for the capture of bulk plasma phenomena inside these thrusters within reasonable computational times. Regions of the thruster with complex magnetic field arrangements (such as those near eroded walls and magnets) and/or reduced Hall parameter (such as those near the anode and the cathode plume) challenge the validity of the quasi-one-dimensional assumption for the electrons. This paper reports on the development of a computer code that solves numerically the 2-D axisymmetric vector form of Ohm's law, with no assumptions regarding the rate of electron transport in the parallel and perpendicular directions. The numerical challenges related to the large disparity of the transport coefficients in the two directions are met by solving the equations in a computational mesh that is aligned with the magnetic field. The fully-2D approach allows for a large physical domain that extends more than five times the thruster channel length in the axial direction, and encompasses the cathode boundary. Ions are treated as an isothermal, cold (relative to the electrons) fluid, accounting for charge-exchange and multiple-ionization collisions in the momentum equations. A first series of simulations of two Hall thrusters, namely the BPT-4000 and a 6-kW laboratory thruster, quantifies the significance of ion diffusion in the anode region and the importance of the extended physical domain on studies related to the impact of the transport coefficients on the electron flow field.

plasma↗

A frost formation model and its validation under various experimental conditions

A numerical model that was used to calculate the frost properties for all regimes of frost growth is described. In the first regime of frost growth, the initial frost density and thickness was modeled from the theories of crystal growth. The 'frost point' temperature was modeled as a linear interpolation between the dew point temperature and the fog point temperature, based upon the nucleating capability of the particular condensing surfaces. For a second regime of frost growth, the diffusion model was adopted with the following enhancements: the generalized correlation of the water frost thermal conductivity was applied to practically all water frost layers being careful to ensure that the calculated heat and mass transfer coefficients agreed with experimental measurements of the same coefficients.

Dietenberger, M. A.↗

The O2 reduction at the IFC modified O2 fuel cell electrode

The International Fuel Corporation (IFC) state of the art (SOA) O2 electrode (Au-10 percent Pt electrocatalyst by weight) is currently being used in the alkaline H2-O2 fuel cell in the NASA Space Shuttle. Recently, IFC modified O2 electrode, as a possible replacement for the SOA electrode. In the present study, O2 reduction data were obtained for the modified electrode at temperatures between 23.3 and 91.7 C. BET measurements gave an electrode BET surface area of about 2070 sq. cm/sq. cm of geometric surface area. The Tafel data could be fitted to two straight line regions. The slope for the lower region, designated as the 0.04 V/decade region, was temperature dependent, and the transfer coefficient was about 1.5. The 'apparent' energy of activation for this region was about 19 kcal/mol. An O2 reduction mechanism for this 0.04 region is presented. In the upper region, designated as the 0.08 V/decade region, diffusion may be the controlling process. Tafel data are presented to illustrate the increase in performance with increasing temperature.

Fielder, William L.↗

Numerical simulation of life cycles of advection warm fog

The formation, development and dissipation of advection warm fog is investigated. The equations employed in the model include the equation of continuity, momentum and energy for the descriptions of density, wind component and potential temperature, respectively, together with two diffusion equations for the modification of water-vapor mixing ratio and liquid-water mixing ratios. A description of the vertical turbulent transfer of heat, moisture and momentum has been taken into consideration. The turbulent exchange coefficients adopted in the model are based on empirical flux-gradient relations.

Hung, R. J.↗

A model for microwave intensity propagation in an inhomogeneous medium

A combined analytic and phenomenological approach, utilizing Maxwell's equations in the Born approximation with radiative transfer theory, is used to describe the propagation of microwave intensity in a scattering medium characterized by three-dimensional random fluctuations in refractive index, as well as nonrandom variations in permittivity, temperature, and loss. This approach yields microwave intensities as a function of polarization, direction, and position. Numerical techniques are presented to solve the transport equations, which include cases of spatially varying coefficients, and highly peaked phase functions. Some computed results illustrating the behavior of microwave intensity in various media are presented. Included are the angular and frequency spectra of thermal emission from semi-infinite media, and the diffuse transmission and reflection response of a scattering layer. The effects of scatterer geometry and scale sizes; correlation function; and gradients in temperature, loss, and scattering parameters are also demonstrated. This model should be particularly useful in interpreting active and passive remote sensing data.

Fisher, A. D.↗

Theory of finite disturbances in a centrifugal compression system with a vaneless radial diffuser

A previous small perturbation analysis of circumferential waves in circumferential compression systems, assuming inviscid flow, is shown to be consistent with observations that narrow diffusers are more stable than wide ones, when boundary layer displacement effect is included. The Moore-Greitzer analysis for finite strength transients containing both surge and rotating stall in axial machines is adapted for a centrifugal compression system. Under certain assumptions, and except for a new second order swirl, the diffuser velocity field, including resonant singularities, can be carried over from the previous inviscid linear analysis. Nonlinear transient equations are derived and applied in a simple example to show that throttling through a resonant value of flow coefficient must occur in a sudden surge-like drop, accompanied by a transient rotating wave. This inner solution is superseded by an outer surge response on a longer time scale. Surge may occur purely as result of circumferential wave resonance. Numerical results are shown for various parametric choices relating to throttle schedule and the characteristic slope. A number of circumferential modes considered simultaneously is briefly discussed.

Moore, F. K.↗