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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Computation of viscous flow in planar and axisymmetric ducts by an implicit marching procedure

A streamwise marching procedure, approximately 200 times faster than a full Navier-Stokes procedure with comparable accuracy, is presented for solving problems of compressible viscous subsonic flow. Results are presented and compared with experimental data for the cases of developing turbulent flow in a circular pipe; turbulent flow in a two-dimensional S-duct; and turbulent flow in a typical subsonic diffuser. Prior to each main marching step, a preliminary marching step is taken in which the integral continuity equation and an uncoupled form of the streamwise momentum equation are solved simultaneously to obtain the viscous pressure correction. During the main marching step the equations for continuity, streamwise momentum, cross-flow momentum, and energy are solved simultaneously as a coupled system using an implicit finite-difference method, with the viscous pressure correction treated as a source term. The analysis may be used for flows with both favorable and adverse pressure gradients and to predict the location of flow separation.

Towne, C. E.↗
NASA NTRS1984

A convective model for turbulent mixing in rotating convection zones

The effects of rotation are included in an analytical model for the convective motions in a plane-parallel layer of an ideal fluid. The turbulent stress tensor, formed by taking products and averages of the various velocity components, is calculated for an arbitrary eddy size and shape. Heuristic formulae presented for determining the size and shape of the dominant eddy then give a fully specified stress tensor. Applications for this stress tensor in problems of stellar internal dynamics, heat flow, scalar diffusion, and dynamo theory are suggested. The resultant stresses tend to produce differential rotation profiles with rapidly rotating equators and interiors. The dynamo activity associated with these convective motions tends to occur near the lower boundary of the convection zone.

Hathaway, D. H.↗
NASA NTRS1988

Analytical solutions of a class of multidimensional Fokker-Planck equations

Analytical representations of the solutions of Kolmogorov's forward (i.e., Fokker-Planck) and backward partial differential equations are obtained by employing Girsanov's (1960) transformation. The equations considered are restricted to those with the drift vector being the gradient of a function and the diffusion matrix being nonsingular and dependent upon time only.

Zhang, Weijian↗
NASA NTRS1989

Calculating Flows In Turbomachine Channels

Noniterative integral-entrainment method yields good approximations. Method of approximate calculation of flow in channel of turbomachine based on interaction of viscous flow in boundary layers with inviscid flow in core layer. Faster and more robust than other approximate methods of same type. Suitable for use in preliminary calculations for design and for off-design operation of turbomachinery. Flows in conical diffuser channels represented by two-dimensional boundary-layer and one-dimensional core flows described by equations of new method.

Schumann, Lawrence F.↗
NASA NTRS1990

Supersonic combustor modeling

The physical phenomena involved when a supersonic flow undergoes chemical reaction are discussed. Detailed physical models of convective and diffusive mixing, and finite rate chemical reaction in supersonic flow are presented. Numerical algorithms used to solve the equations governing these processes are introduced. Computer programs using these algorithms are used to analyze the structure of the reacting mixing layer. It is concluded that, as in subsonic flow, exothermic heat release in unconfined supersonic flows retards fuel/air mixing. Non mixing is shown to be a potential problem in reducing the efficiency of supersonic as well as subsonic combustion. Techniques for enhancing fuel/air mixing and combustion are described.

Riggins, David W.↗
NASA NTRS1990

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.↗
NASA NTRS1992

Collection Of Charge By Diffusion From Ion Track

Collection of electric charge carriers from track of energetic ion through silicon electronic device studied using diffusion model representing transport of charge from ion track to sink. Effects of ionizing radiation approximated by simple equations.

Edmonds, Larry D.↗
NASA NTRS1993

Renormalization group estimates of transport coefficients in the advection of a passive scalar by incompressible turbulence

The advection of a passive scalar by incompressible turbulence is considered using recursive renormalization group procedures in the differential sub grid shell thickness limit. It is shown explicitly that the higher order nonlinearities induced by the recursive renormalization group procedure preserve Galilean invariance. Differential equations, valid for the entire resolvable wave number k range, are determined for the eddy viscosity and eddy diffusivity coefficients, and it is shown that higher order nonlinearities do not contribute as k goes to 0, but have an essential role as k goes to k(sub c) the cutoff wave number separating the resolvable scales from the sub grid scales. The recursive renormalization transport coefficients and the associated eddy Prandtl number are in good agreement with the k-dependent transport coefficients derived from closure theories and experiments.

Zhou, YE↗
NASA NTRS1994

Approximate Solution Methods for Spectral Radiative Transfer in High Refractive Index Layers

Some ceramic materials for high temperature applications are partially transparent for radiative transfer. The refractive indices of these materials can be substantially greater than one which influences internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat conduction and heating at the boundaries by convection and radiation. Two-flux and diffusion methods are investigated here to obtain approximate solutions using a simpler formulation than required for exact numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent results for gray and two band spectral calculations. The diffusion method yields a good approximation for spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A hybrid spectral model is developed using the two-flux method in the optically thin bands, and radiative diffusion in bands that are optically thick.

Siegel, R.↗
NASA NTRS2004

Research on the Natural Variability of Climate

The scheme for ocean vertical mixing presented in Canuto, Howard, Cheng & Dubovikov and Canuto, Howard, Hogan, Cheng, Dubovikov & Montenegro [GISS mixing scheme] was extended to include the latitudinal dependence reported by Gregg, Sanford & Winkel of the input to the interior ocean turbulence field of energy from internal waves. The resulting latitude dependence was implemented in our module of code for the GISS vertical mixing scheme and tested in the global NCAR-CSM ocean model with 3 degree X 3 degree resolution and 25 levels in the vertical. With the latitude dependence, the background diffusivity decreases from approx. 0.1 sq cm/sec at midlatitudes to approx. less than 0.01 sq cm/sec at the equator. A significant improvement was seen in the tropical ocean, with tropical thermoclines being sharpened in agreement with the requirements of observational studies and the needs of ENS0 modeling. At the same time global measures nf performance such as the meridional overturning and northward heat transport were not adversely affected. Results were presented at the 2004 Layered Ocean Modeling meeting and a paper has been prepared and submitted to Geophysical Research Letters.

Source record↗
NASA NTRS2016

Budget of Turbulent Kinetic Energy in a Shock Wave Boundary-Layer Interaction

Implicit large-eddy simulation (ILES) of a shock wave/boundary-layer interaction (SBLI) was performed. Quantities present in the exact equation of the turbulent kinetic energy transport were accumulated and used to calculate terms like production, dissipation, molecular diffusion, and turbulent transport. The present results for a turbulent boundary layer were validated by comparison with direct numerical simulation data. It was found that a longer development domain was necessary for the boundary layer to reach an equilibrium state and a finer mesh resolution would improve the predictions. In spite of these findings, trends of the present budget match closely with that of the direct numerical simulation. Budgets for the SBLI region are presented at key axial stations. These budgets showed interesting dynamics as the incoming boundary layer transforms and the terms of the turbulent kinetic energy budget change behavior within the interaction region.

shock wave/boundary-layer interaction↗
NASA NTRS1994

A Note on the Wave Action Density of a Viscous Instability Mode on a Laminar Free-shear Flow

Using the assumptions of an incompressible and viscous flow at large Reynolds number, we derive the evolution equation for the wave action density of an instability wave traveling on top of a laminar free-shear flow. The instability is considered to be viscous; the purpose of the present work is to include the cumulative effect of the (locally) small viscous correction to the wave, over length and time scales on which the underlying base flow appears inhomogeneous owing to its viscous diffusion. As such, we generalize our previous work for inviscid waves. This generalization appears as an additional (but usually non-negligible) term in the equation for the wave action. The basic structure of the equation remains unaltered.

Balsa, Thomas F.↗
NASA NTRS1972

Purfication kinetics of beryllium during vacuum induction melting

The kinetics of evaporation in binary alloys were quantitatively treated. The formalism so developed works well for several systems studied. The kinetics of purification of beryllium was studied through evaporation data actually acquired during vacuum induction melting. Normal evaporation equations are shown to be generally valid and useful for understanding the kinetics of beryllium purification. The normal evaporation analysis has been extended to cover cases of limited liquid diffusion. It was shown that under steady-state evaporation, the solute concentration near the surface may be up to six orders of magnitude different from the bulk concentration. Corrections for limited liquid diffusion are definitely needed for the highly evaporative solute elements, such as Zn, Mg, and Na, for which the computed evaporation times are improved by five orders of magnitude. The commonly observed logarithmic relation between evaporation time and final concentration further supports the validity of the normal evaporation equations.

Mukherjee, J. L.↗
NASA NTRS1973

Solutions of the chemical kinetic equations for initially inhomogeneous mixtures.

Following the recent discussions by O'Brien (1971) and Donaldson and Hilst (1972) of the effects of inhomogeneous mixing and turbulent diffusion on simple chemical reaction rates, the present report provides a more extensive analysis of when inhomogeneous mixing has a significant effect on chemical reaction rates. The analysis is then extended to the development of an approximate chemical sub-model which provides much improved predictions of chemical reaction rates over a wide range of inhomogeneities and pathological distributions of the concentrations of the reacting chemical species. In particular, the development of an approximate representation of the third-order correlations of the joint concentration fluctuations permits closure of the chemical sub-model at the level of the second-order moments of these fluctuations and the mean concentrations.

Hilst, G. R.↗
NASA NTRS1975

A modified Henyey method for computing radiative transfer hydrodynamics

The implicit hydrodynamic code of Kutter and Sparks (1972), which is limited to optically thick regions and employs the diffusion approximation for radiative transfer, is modified to include radiative transfer effects in the optically thin regions of a model star. A modified Henyey method is used to include the solution of the radiative transfer equation in this implicit code, and the convergence properties of this method are proven. A comparison is made between two hydrodynamic models of a classical Cepheid with a 12-day period, one of which was computed with the diffusion approximation and the other with the modified Henyey method. It is found that the two models produce nearly identical light and velocity curves, but differ in the fact that the former never has temperature inversions in the atmosphere while the latter does when sufficiently strong shocks are present.

Karp, A. H.↗
NASA NTRS1988

Plasma instabilities stimulated by HF transmitters in space

Diffuse incoherent signal returns are often observed on Alouette and ISIS topside ionograms in addition to coherent echoes of electromagnetic and electrostatic waves. These diffuse signals, which at times can be the dominant features on topside ionograms, have been attributed to sounder-induced temperature anisotropies which drive the Harris instability. Previous theoretical investigations were based on the electrostatic approximation to the dispersion equation. The present paper will present calculations indicating that when the electromagnetic terms are retained in the dispersion equation and when the sounder-stimulated perpendicular electron temperature approaches 1 keV, then the whistler mode can have a temporal growth rate larger than the electrostatic electron cyclotron harmonic wave mode central to the diffuse resonance problem. Present sounders lack the power and antenna lengths to generate whistler mode waves in this manner. In addition, such waves would have large group velocities and would quickly leave the vicinity of the sounder. Experiments to investigate the wave growth, propagation, and damping of such stimulated waves are planned for the 1990s using a highly flexible sounder on the Space Shuttle and a receiver on a subsatellite.

Benson, Robert F.↗
NASA NTRS1977

Celestial diffuse gamma-ray emission observed by SAS-2 and its interpretation

A clearly established diffuse celestial gamma-ray component was seen by SAS-2 above 35 MeV, after examining several regions of the sky at different latitudes, including the north celestial pole. For energies above 100 MeV the gamma ray results are consistent with an equation of the form I(b)=C1+C2/sin b with the second term being dominant, suggesting that the radiation above 100 MeV comes largely from the local regions of the galactic disk. Between 35 and 100 MeV, a similar equation is also a reasonable representation of the data, but here the two terms are comparable, with the first, or isotropic term, actually being the larger one. In addition to indicating that the diffuse radiation is partially galactic, these results imply a steepness for the energy spectrum of the diffuse isotropic component which places significant constraints on possible theoretical models of this radiation.

Fichtel, C. E.↗
NASA NTRS1986

Thermophoretically augmented mass transfer rates to solid walls across laminar boundary layers

Predictions of mass transfer (heavy vapor and small particle deposition) rates to solid walls, including the effects of thermal (Soret) diffusion ('thermophoresis' for small particles), are made by numerically solving the two-dimensional self-similar forced convection laminar boundary-layer equations with variable properties, covering the particle size range from vapor molecules up to the size threshold for inertial (dynamical nonequilibrium) effects. The effect of thermophoresis is predicted to be particularly important for submicron particle deposition on highly cooled solid surfaces, with corresponding enhancement factors at atmospheric conditions being over a thousand-fold at T(w)/T(e) equal to about 0.6. As a consequence of this mass transfer mechanism, the particle size dependence of the mass transfer coefficient to a cooled wall will be much weaker than for the corresponding case of isothermal capture by Brownian-convective diffusion.

Gokoglu, S. A.↗