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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

NASA NTRS1986

Diffusion flame extinction in slow convection flow under microgravity environment

A theoretical analysis is presented to study the extinction characteristics of a diffusion flame near the leading edge of a thin fuel plate in slow, forced convective flows in a microgravity environment. The mathematical model includes two-dimensional Navier-Stokes momentum, energy and species equations with one-step overall chemical reaction using second-order finite rate Arrhenius kinetics. Radiant heat loss on the fuel plate is applied in the model as it is the dominant mechanism for flame extinguishment in the small convective flow regime. A parametric study based on the variation of convective flow velocity, which varies the Damkchler number (Da), and the surface radiant heat loss parameter (S) simultaneously, is given. An extinction limit is found in the regime of slow convective flow when the rate of radiant heat loss from fuel surface outweighs the rate of heat generation due to combustion. The transition from existent envelope flame to extinguishment consists of gradual flame contraction in the opposed flow direction together with flame temperature reduction as the convective flow velocity decreases continuously until the extinction limit is reached. A case of flame structure subjected to surface radiant heat loss is also presented and discussed.

Chen, Chiun-Hsun↗

NASA NTRS1990

Filtering analysis of a direct numerical simulation of the turbulent Rayleigh-Benard problem

A filtering analysis of a turbulent flow was developed which provides details of the path of the kinetic energy of the flow from its creation via thermal production to its dissipation. A low-pass spatial filter is used to split the velocity and the temperature field into a filtered component (composed mainly of scales larger than a specific size, nominally the filter width) and a fluctuation component (scales smaller than a specific size). Variables derived from these fields can fall into one of the above two ranges or be composed of a mixture of scales dominated by scales near the specific size. The filter is used to split the kinetic energy equation into three equations corresponding to the three scale ranges described above. The data from a direct simulation of the Rayleigh-Benard problem for conditions where the flow is turbulent are used to calculate the individual terms in the three kinetic energy equations. This is done for a range of filter widths. These results are used to study the spatial location and the scale range of the thermal energy production, the cascading of kinetic energy, the diffusion of kinetic energy, and the energy dissipation. These results are used to evaluate two subgrid models typically used in large-eddy simulations of turbulence. Subgrid models attempt to model the energy below the filter width that is removed by a low-pass filter.

Eidson, T. M.↗

NASA NTRS1991

Improved millimeter-wave mixer performance analysis at cryogenic temperatures

The results of a 600-GHz mixer performance analysis using an improved model for computing the Schottky diode capacitance-voltage (C-V) relationship are presented. The computed C-V data for a realistic submillimeter-waver mixer diode are given as a function of physical temperature and compared to the standard analytic expression based on a solution of Poisson's equation. Both C-V relationships are used to predict the performance of an ideally terminated 600 GHz mixer operating at 300, 140, 80, and 30 K. It is shown that the drift-diffusion model more accurately describes the mixer performance when the physical temperature is reduced below roughly 100 K.

Siegel, Peter H.↗

NASA NTRS1939

The drag of airplane radiators with special reference to air heating : comparison of theory and experiment

This report contains a survey of past radiator research. This report also is intended as a systematic comparison of theoretical and experimental radiator drag, with the object of ascertaining the most important loss sources and their interaction in different cases of installation, and to separate the radiator systems which are amenable to calculation, both as regards axial flow and drag. The sources of loss due to the diffuser are to be looked into closely as in many cases they can be of preeminent magnitude and their customary appraisal, according to Fliegner's formula, does not meet actual conditions. Besides, generally applicable equations and charts are developed for the rapid determination of the heating effect of radiators as regards flow and drag, and then checked by routine tests on hot radiators.

Gothert, B↗

NASA NTRS1992

Calibration of a helium-cooled infrared spatial radiometer and grating spectrometer

Methods used by the Space Dynamics Laboratory of Utah State University (SDL/USU) to calibrate infrared sensors are described, using the Infrared Background Signature Survey (IBSS) spatial radiometer and grating spectrometer as examples. A calibration equation and a radiometric model are given for each sensor to describe their responsivity in terms of individual radiometric parameters. The calibration equation terms include dark offset, linearity, absolute responsivity, and measurement uncertainty, and the radiometric model domains include spatial, spectral, and temporal domains. A portable calibration facility, designed and fabricated by SDL/USU, provided collimated, extended, diffuse scatter, and Jones sources in a single cryogenic dewar. This multi-function calibrator allowed calibration personnel to complete a full calibration of the IBSS infrared radiometer and spectrometer in two 15-day periods. A calibration data system was developed to control and monitor the calibration facility, and to record and analyze sensor data.

Jacobsen, Larry↗

NASA NTRS2006

Comparison of the Radiative Two-Flux and Diffusion Approximations

Approximate solutions are sometimes used to determine the heat transfer and temperatures in a semitransparent material in which conduction and thermal radiation are acting. A comparison of the Milne-Eddington two-flux approximation and the diffusion approximation for combined conduction and radiation heat transfer in a ceramic material was preformed to determine the accuracy of the diffusion solution. A plane gray semitransparent layer without a substrate and a non-gray semitransparent plane layer on an opaque substrate were considered. For the plane gray layer the material is semitransparent for all wavelengths and the scattering and absorption coefficients do not vary with wavelength. For the non-gray plane layer the material is semitransparent with constant absorption and scattering coefficients up to a specified wavelength. At higher wavelengths the non-gray plane layer is assumed to be opaque. The layers are heated on one side and cooled on the other by diffuse radiation and convection. The scattering and absorption coefficients were varied. The error in the diffusion approximation compared to the Milne-Eddington two flux approximation was obtained as a function of scattering coefficient and absorption coefficient. The percent difference in interface temperatures and heat flux through the layer obtained using the Milne-Eddington two-flux and diffusion approximations are presented as a function of scattering coefficient and absorption coefficient. The largest errors occur for high scattering and low absorption except for the back surface temperature of the plane gray layer where the error is also larger at low scattering and low absorption. It is shown that the accuracy of the diffusion approximation can be improved for some scattering and absorption conditions if a reflectance obtained from a Kubelka-Munk type two flux theory is used instead of a reflection obtained from the Fresnel equation. The Kubelka-Munk reflectance accounts for surface reflection and radiation scattered back by internal scattering sites while the Fresnel reflection only accounts for surface reflections.

Spuckler, Charles M.↗

NASA NTRS2011

Effect of Water on the Thermo-Mechanical Behavior of Carbon Cloth Phenolic

The results of thermo-mechanical experiments, which were conducted previously by one of the authors, are reviewed. The strain in the direction normal to the fabric plane was measured as a function of temperature for a variety of initial moisture contents and heating rates. In this paper, the general features of the thermo-mechanical response are discussed and the effect of heating rate and initial moisture content are highlighted. The mechanical interaction between the phenolic polymer and water trapped within its free volumes as the polymer is heated to high temperatures is discussed. An equation for the internal stresses which are generated within the polymer due to trapped water is obtained from the total stress expression for a binary mixture of polymer and water. Numerical solutions for moisture diffusion in the thermo-mechanical experiments were performed and the results of these solutions are presented. The results of the moisture diffusion solutions help to explain the effects of heating rate and moisture content on the strain behavior normal to the fabric plane.

Sullivan, Roy M.↗

NASA NTRS1979

Analysis of turbulent free-jet hydrogen-air diffusion flames with finite chemical reaction rates

A numerical analysis is presented of the nonequilibrium flow field resulting from the turbulent mixing and combustion of an axisymmetric hydrogen jet in a supersonic parallel ambient air stream. The effective turbulent transport properties are determined by means of a two-equation model of turbulence. The finite-rate chemistry model considers eight elementary reactions among six chemical species: H, O, H2O, OH, O2 and H2. The governing set of nonlinear partial differential equations was solved by using an implicit finite-difference procedure. Radial distributions were obtained at two downstream locations for some important variables affecting the flow development, such as the turbulent kinetic energy and its dissipation rate. The results show that these variables attain their peak values on the axis of symmetry. The computed distribution of velocity, temperature, and mass fractions of the chemical species gives a complete description of the flow field. The numerical predictions were compared with two sets of experimental data. Good qualitative agreement was obtained.

Sislian, J. P.↗

NASA NTRS1972

Surface composition of magnetic neutron stars.

The relative abundances of seven constituent nuclei, He4, C12, O16, Ne20, Mg24, Si28, and Fe56, are calculated as a function of time for neutron star atmospheres within which exist magnetic fields of the order of 10 to the 13th G. The opacity, equation of state of the electrons, and cooling rate of the magnetic star are discussed, and it is shown to be a reasonable approximation to assume an atmosphere to be isothermal. The effects of particle diffusion are included in the nuclear reaction network. Computations are performed both for a constant mass atmosphere and for an atmosphere in which mass is being ejected. It is found that the final abundances are model-independent, as long as the initial model contains predominantly He4. The relative abundances are compared to the cosmic ray spectrum. For both the constant-mass and mass-loss atmospheres, nucleosynthesis proceeds virtually completely to Fe56. However the outermost layers of the envelope, in which no mass is being ejected, are composed almost entirely of He4 with trace amounts of Fe56. After the loss of about 10 to the 21st g, only Fe56 is ejected from atmospheres expelling mass.

Rosen, L. C.↗

NASA NTRS1974

Development of a second order closure model for computation of turbulent diffusion flames

A typical eddy box model for the second-order closure of turbulent, multispecies, reacting flows developed. The model structure was quite general and was valid for an arbitrary number of species. For the case of a reaction involving three species, the nine model parameters were determined from equations for nine independent first- and second-order correlations. The model enabled calculation of any higher-order correlation involving mass fractions, temperatures, and reaction rates in terms of first- and second-order correlations. Model predictions for the reaction rate were in very good agreement with exact solutions of the reaction rate equations for a number of assumed flow distributions.

Varma, A. K.↗

NASA NTRS1984

Reducing numerical diffusion for incompressible flow calculations

A number of approaches for improving the accuracy of incompressible, steady-state flow calculations are examined. Two improved differencing schemes, Quadratic Upstream Interpolation for Convective Kinematics (QUICK) and Skew-Upwind Differencing (SUD), are applied to the convective terms in the Navier-Stokes equations and compared with results obtained using hybrid differencing. In a number of test calculations, it is illustrated that no single scheme exhibits superior performance for all flow situations. However, both SUD and QUICK are shown to be generally more accurate than hybrid differencing.

Claus, R. W.↗

NASA NTRS1986

Prediction of heat release effects on a mixing layer

A fully second-order closure model for turbulent reacting flows is suggested based on Favre statistics. For diffusion flames the local thermodynamic state is related to single conserved scalar. The properties of pressure fluctuations are analyzed for turbulent flows with fluctuating density. Closure models for pressure correlations are discussed and modeled transport equations for Reynolds stresses, turbulent kinetic energy dissipation, density-velocity correlations, scalar moments and dissipation are presented and solved, together with the mean equations for momentum and mixture fraction. Solutions of these equations are compared with the experimental data for high heat release free mixing layers of fluorine and hydrogen in a nitrogen diluent.

Farshchi, M.↗

NASA NTRS1988

Euler equation computations for the flow over a hovering helicopter rotor

A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.

Roberts, Thomas Wesley↗

NASA NTRS2004

Large-scale Atmospheric Transport Processes

Continuing earlier work, we continued an investigation of the seasonal behavior of the edges of the stratospheric surf zone. These edges form a barrier between the rapidly mixed surf zone and the relatively isolated tropics. In collaboration with Dr Lynn Sparling at GSFC, we used a statistical analysis of HALOE and CLAES trace gas data from UARS to identify and locate these edges during each UARS observing period. We found that the edges on both sides of the equator are present all year (a fact that is important for conceptual models of stratospheric transport), though that on the summer side of the equator is much less sharp than the winter edge. The edges migrate seasonally into the summer hemisphere. Their location also shows influence of the QBO, together with the SAO at higher altitudes. Comparisons with effective diffusivities, and the edge locations, suggest that the edge is sustained by surf zone entrainment during winter, but by the residual circulation during summer.

Plumb, R. Alan↗

NASA NTRS1993

A Realizable Reynolds Stress Algebraic Equation Model

The invariance theory in continuum mechanics is applied to analyze Reynolds stresses in high Reynolds number turbulent flows. The analysis leads to a turbulent constitutive relation that relates the Reynolds stresses to the mean velocity gradients in a more general form in which the classical isotropic eddy viscosity model is just the linear approximation of the general form. On the basis of realizability analysis, a set of model coefficients are obtained which are functions of the time scale ratios of the turbulence to the mean strain rate and the mean rotation rate. The coefficients will ensure the positivity of each component of the mean rotation rate. These coefficients will ensure the positivity of each component of the turbulent kinetic energy - realizability that most existing turbulence models fail to satisfy. Separated flows over backward-facing step configurations are taken as applications. The calculations are performed with a conservative finite-volume method. Grid-independent and numerical diffusion-free solutions are obtained by using differencing schemes of second-order accuracy on sufficiently fine grids. The calculated results are compared in detail with the experimental data for both mean and turbulent quantities. The comparison shows that the present proposal significantly improves the predictive capability of K-epsilon based two equation models. In addition, the proposed model is able to simulate rotational homogeneous shear flows with large rotation rates which all conventional eddy viscosity models fail to simulate.

Shih, Tsan-Hsing↗

NASA NTRS1972

Low-energy solar positrons.

Fluxes of solar particles, consistent with Kinsey's (1970) measurements and demodulated in a manner consistent with the cosmic-ray transport equations discussed by Gleeson and Axford (1967), are shown to be able to produce low-energy positrons in the solar atmosphere via unstable CNO beta-particle emitters. Computed positron fluxes are presented for various diffusion coefficient values and compared with measurements of Cline and Hones (1970).

Gregory, C. T.↗

NASA NTRS1975

Scatter-free collimated convection and cosmic-ray transport at 1 AU

Observations of solar cosmic rays between 0.1 MeV and 1 GeV are more consistent with scatter-free than diffusive propagation near 1 AU. Since these particles also are constrained to interplanetary magnetic field lines, the parallel and transverse flow are decoupled by collimated convection. Local propagation is therefore more appropriately described by Liouville's equation, and transport properties can be examined from the exact zeroeth and first moments of the equation averaged over all directions in momentum space. Observational results which may be explained directly without invoking local scattering include the decay of solar flare events as well as the radial gradient, diurnal and semidiurnal variation of galactic cosmic rays.

Roelof, E. C.↗

NASA NTRS1981

A high-resolution study of the gas flow in barred spirals

The considered investigation is concerned with conditions in that type of spiral galaxy which is known as a 'barred spiral'. The lack of star formation and the smoothness of the dark dust lanes may indicate that the bar region is pervaded by a more homogeneous interstellar medium with a lower viscosity than that characteristic of the outer regions. The effects of such a reduced viscosity are investigated on the basis of high-resolution time-dependent gasdynamical equations. Calculations modeling the gas flow in a given bar potential are presented. The calculations were performed on grids with various mesh sizes and consequently had different resolutions and numerical viscosities. A reduction in the numerical diffusion and viscosity was accompanied by changes in the flow pattern, including the appearance of a strong shock in the bar region and the occurrence of postshock outflow.

Van Albada, G. D.↗