Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical diffusion”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

A semi-elliptic analysis for 2-D viscous flows through cascade configurations

A semielliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is negligible but which are significantly influenced by upstream interactions. A two-step alternating-direction-explicit numerical scheme is developed to solve the resulting governing equations efficiently. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Various simple channel and cascade flow configurations are examined for several values of Re, thickness ratio and Mach number, in order to establish the applicability of the IPNS model for these flows. With this model, solutions to compressible as well as nearly incompressible flows are obtained without any modification either in the analysis or in the solution procedure.

Ramamurti, R.↗
NASA NTRS1993

Diffusion of Na and K in the uppermost regolith of Mercury

We examine the conditions under which ions impacting the Hermean surface can act as a regional source of enhanced atmospheric column through ion implantation and subsequent release at the surface. At most latitudes, energetic (a few keV) Na or K ions which impact the nightside surface, are released quickly (well before noon) upon warming of the surface. We show that the relative sunrise/sunset difference produced by ion implantation is ((n)(sub SR) - (n)(sub SS)/(n)(sub Ave) = f, where f is the fraction of the photo-ions recycled, (n)(sub SR) is the average zenith column above the sunrise portion of the illuminated hemisphere, (n)(sub SS) is the average zenith column in the sunset portion, and (n)(sub Ave) is the average column over the sunlit disk. Thus, to produce a large sunrise/sunset difference via ion implantation and subsequent release requires efficient (close to total) recycling. We show that the most extensive set of available data reduced to Na column abundance does not show any sunrise enhancement. We argue that the K data do not permit an unambiguous interpretation in favor of sunrise/sunset differences. We further find that if an efficient surface loss process for the alkali is not operating after sunrise, the initial Na and K distribution will relax into the bulk of the solid. Preserved abundance gradients of Na and K in lunar glasses that the Arrhenius coefficients for impact glasses are likely more modest than those for laboratory glasses derived from rock by nonimpact processes, but the effect of the more modest diffusion rates is only to delay the efficient loss of Na by a few Earth days. We argue that implantation can lead to observable regional increases in the observed Na or K column densities only if it occurs at very high latitudes, where diffusion is slow. It is typically lower energy ions which impact at high latitude and these are both more numerous than the high energy ions and possess smaller average penetration depths; thus there are additional reasons to favor a high-latitude locus for any possible prompt return related increases in zenith column. We find that we cannot rule out sputtering as a source process, as the low-energy ions are efficient sputterers. Finally, we argue that the observed Na/K ratio in the atmopshere may be the expression of their different asymptotic rates of loss from the interiors of the regolith grains. If so, this allows us to fix the importance of sputtering and photon stimulated desorption relative to impact vaporization.

Killen, R. M.↗
NASA NTRS1999

A Textbook for a First Course in Computational Fluid Dynamics

This paper describes and discusses the textbook, Fundamentals of Computational Fluid Dynamics by Lomax, Pulliam, and Zingg, which is intended for a graduate level first course in computational fluid dynamics. This textbook emphasizes fundamental concepts in developing, analyzing, and understanding numerical methods for the partial differential equations governing the physics of fluid flow. Its underlying philosophy is that the theory of linear algebra and the attendant eigenanalysis of linear systems provides a mathematical framework to describe and unify most numerical methods in common use in the field of fluid dynamics. Two linear model equations, the linear convection and diffusion equations, are used to illustrate concepts throughout. Emphasis is on the semi-discrete approach, in which the governing partial differential equations (PDE's) are reduced to systems of ordinary differential equations (ODE's) through a discretization of the spatial derivatives. The ordinary differential equations are then reduced to ordinary difference equations (O(Delta)E's) using a time-marching method. This methodology, using the progression from PDE through ODE's to O(Delta)E's, together with the use of the eigensystems of tridiagonal matrices and the theory of O(Delta)E's, gives the book its distinctiveness and provides a sound basis for a deep understanding of fundamental concepts in computational fluid dynamics.

Zingg, D. W.↗
NASA NTRS1990

Thermal Analysis

Packaging engineers can avoid unnecessary and costly redesign of product through use of tools such as the SINDA '85/FLUINT software system for predicting thermal performance of electronic aviation equipment. SINDA '85/FLUINT is an acronym for Systems Improved Numerical Differencing Analyzer and Fluid Integrators. It is a program for solving the physical problems governed by flow and diffusion type equations and is most widely used as a general thermal analyzer.

Source record↗
NASA NTRS1995

Design, construction and testing of annular diffusers for high speed civil transportation combustor applications

A theoretical and an experimental design study of subsonic flow through curved-wall annular diffusers is being carried out in order to establish the most pertinent design parameters for such devices and the implications of their application in the design of engine components in the aerospace industries. This investigation consists of solving numerically the full Navier Stokes and Continuity equations for the time-mean flow. Various models of turbulence are being evaluated for adoption throughout the study and comparisons would be made with experimental data where they exist. Assessment of diffuser performance based on the dissipated mechanical energy would also be made. The experimental work involves the application of Computer Aided Design software tool to the development of a suitable annular diffuser geometry and the subsequent downloading of such data to a CNC machine at Central State University. The results of the investigations are expected to indicate that more cost effective component design of such devices as effective component design of such devices as diffusers which normally contain complex flows can still be achieved. In this regard a review paper was accepted and presented at the First International Conference on High Speed Civil Transportation Research held at North Carolina A&T in December of 1994.

Okhio, Cyril B.↗
NASA NTRS1990

Cellular structure of lean hydrogen flames in microgravity

Detailed, time-dependent, two-dimensional numerical simulations of premixed laminar flames have been used to study the initiation and subsequent development of cellular structures in lean hydrogen-air flames. The model includes detailed hydrogen-oxygen combustion with 24 elementary reactions of eight reactive species and a nitrogen diluent, molecular diffusion of all species, thermal conduction, viscosity, and convection. This model has been used to study the nonlinear evolution of cellular flame structure and shows that cell splitting, as observed in experiments, can be predicted numerically for sufficiently reactive mixtures. The structures that evolved also resembled the cellular structures observed in experiments. The present study shows that the 'cell-split limit' postulated from experimental observations is an intrinsic property of the mixture and that external factors such as heat losses are not necessary to cause this limit.

Patnaik, G.↗
NASA NTRS1999

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗
NASA NTRS1982

Effects of unsteady free-stream velocity and free-stream turbulence at a stagnation point

The combined effects of transient free stream velocity and turbulence at a stagnation point on a cylinder situated in a crossflow is investigated analytically, and a model is formulated for the eddy diffusivity induced by free-stream turbulence. The steepest descent method is used to integrate the governing momentum expression, and numerical solutions are given for the unsteady wall shear stress function for specific free-stream transients. It is found after correlation of the results by means of a new turbulence parameter that wall friction increases with increasing free-stream turbulence intensity, and that the friction factor increases with increasing reduced frequency of oscillation values.

Gorla, R. S. R.↗
NASA NTRS1982

Finite element analysis of the effect of a non-planar solid-liquid interface on the lateral solute segregation during unidirectional solidification

The effect of solid-liquid interface shape on lateral solute segregation during steady-state unidirectional solidification of a binary mixture is calculated under the assumption of no convection in the liquid. A finite element technique is employed to compute the concentration field in the liquid and the lateral segregation in the solid with a curved boundary between the liquid and solid phases. The computational model is constructed assuming knowledge of the solid-liquid interface shape; no attempt is made to relate this shape to the thermal field. The influence of interface curvature on the lateral compositional variation is investigated over a range of system parameters including diffusivity, growth speed, distribution coefficient, and geometric factors of the system. In the limiting case of a slightly nonplanar interface, numerical results from the finite element technique are in good agreement with the analytical solutions of Coriell and Sekerka obtained by using linear theory. For the general case of highly non-planar interface shapes, the linear theory fails and the concentration field in the liquid as well as the lateral solute segregation in the solid can be calculated by using the finite element method.

Carlson, F. M.↗
NASA NTRS1984

Structural Characterization of Solar Silicon

Compared to single crystal silicon, solar silicon generally contains large residual stresses and numerous structural defects. During subsequent processing, the defect structure can undergo further changes since, at the high temperatures required for diffusion, dislocations are sufficiently mobile to rearrange themselves in patterns which reduce long range residual stresses. This process, which is similar to the polygonization of strained metals, has no counterpart in single crystal silicon. Many of the solar silicon materials are grown from graphite dies and crucibles and therefore contain carbon concentrations in excess of 1E18. The interactions between carbon, crystal defects, intrinsic point defects, oxygen and diffusing dopants are discussed.

Ast, D. G.↗
NASA NTRS1986

Unimolecular diffusion-mediated reactions with a nonrandom time-modulated absorbing barrier

A diffusion-reaction model with time-dependent reactivity is formulated and applied to unimolecular reactions. The model is solved exactly numerically and approximately analytically for the unreacted fraction as a function of time. It is shown that the approximate analytical solution is valid even when the system is far from equilibrium, and when the reactivity probability is more complicated than a square-wave function of time. A discussion is also given of an approach to problems of this type using a stochastically fluctuating reactivity, and the first-passage time for a particular example is derived.

Bashford, D.↗
NASA NTRS1988

Simulation of 2-dimensional viscous flow through cascades using a semi-elliptic analysis and hybrid C-H grids

A semi-elliptic formulation, termed the interacting parabolized Navier-Stokes (IPNS) formulation, is developed for the analysis of a class of subsonic viscous flows for which streamwise diffusion is neglible but which are significantly influenced by upstream interactions. The IPNS equations are obtained from the Navier-Stokes equations by dropping the streamwise viscous-diffusion terms but retaining upstream influence via the streamwise pressure-gradient. A two-step alternating-direction-explicit numerical scheme is developed to solve these equations. The quasi-linearization and discretization of the equations are carefully examined so that no artificial viscosity is added externally to the scheme. Also, solutions to compressible as well as nearly compressible flows are obtained without any modification either in the analysis or in the solution process. The procedure is applied to constricted channels and cascade passages formed by airfoils of various shapes. These geometries are represented using numerically generated curilinear boundary-oriented coordinates forming an H-grid. A hybrid C-H grid, more appropriate for cascade of airfoils with rounded leading edges, was also developed. Satisfactory results are obtained for flows through cascades of Joukowski airfoils.

Ramamurti, R.↗
NASA NTRS1993

Analytical and Numerical Modeling of Flame-Balls in Hydrogen-Air Mixtures

Flame-balls are stationary spherical premixed flames observed in certain near-limit mixtures. It is believed that radiative heat losses are an important stabilizing influence. Numerical solutions of flame balls are constructed for hydrogen-air mixtures using an accurate description of the chemical kinetics, diffusive transport, and radiation losses. A lean limit equivalence ratio of 0.0866 is predicted and a rich limit of 2,828. For any equivalence ratio between the two limits there are two solutions. One is characterized by a small flame, incomplete reactant consumption, and negligible radiation losses. The other by a large flame, complete consumption of one of the reactants, and significant radiation losses. The maximum temperature varies between 1200 and 900 K as the two solution branches are traversed. Much of our discussion is a reprise and modification of previously published analytical results, for these provide physical insight into the nature of the solutions, and suggest that a portion of the large flame branch near the lean limit is stable and so corresponds to observable flames.

Buckmaster, J.↗
NASA NTRS1993

Meso-Beta scale numerical simulation studies of terrain-induced jet streak mass/momentum perturbations

The Control-B simulation experiment was designed to improve on the Control-A experiment performed with the GMASS model. This experiment addressed several inadequacies with the first smooth terrain numerical simulation by including: (1) increased nested-grid resolution to better define the simulated gravity waves, (2) increased horizontal diffusion to remove outflow boundary condition noise, and (3) the use of reanalyzed rawinsonde data and surface observations in the initial state to increase the definition of the observed jet streak as well as other low-level features. A smoothed-terrain dry simulation with the nested-grid GMASS model has revealed many important aspects of the processes which resulted in the generation of gravity waves in the region and time when and where they were observed. However, the vertical structure, number, and characteristics of the waves are still quite different from observed waves as diagnosed thus necessitating future improved simulations. However, this control simulation has produced substantial insight into processes which occur on many spatial scales over a 30 hour time period thus allowing one to draw promising inferences as to the mechanisms for the complex process which occurred in nature during the CCOPE case study. The theoretical aspects of the project have focus on understanding the nature of the ageostrophic circulations which are produced in idealized models of the atmosphere in which the troposphere is modeled in one of two ways. The first model assumes that the lower atmosphere can be represented as a single layer of homogeneous fluid whose upper surface is free to exhibit vertical displacement. Two-dimensional internal convergence (divergence) occurring during the adjustment to an asymptotic equilibrium state from an ageostrophic initial state whose momentum structure is representative of a midlatitude localized zonal wind anomaly will cause the free upper surface of the homogeneous atmosphere to rise (fall), and therefore the response can be viewed as being physically three-dimensional. The second model assumes that the troposphere can be represented by an unbounded continuously stratified Boussinesq fluid of constant Brunt-Vaisala frequency N, where the vertical gradient of the basic state potential temperature profile allows for the existence of vertically propagating internal inertia-gravity waves.

Lin, Yuh-Lang↗
NASA NTRS1993

Moisture Content and Migration Dynamics in Unsaturated Porous Media

Fundamental studies of fluid mechanics and transport in partially saturated soils are presented. Solution of transient diffusion problems in support of the development of probes for the in-situ measurement of moisture content is given. Numerical and analytical methods are used to study the fundamental problem of meniscus and saturation front propagation in geometric models of porous media.

Homsy, G. M.↗
NASA NTRS2005

Sol-Gel Glass Holographic Light-Shaping Diffusers

Holographic glass light-shaping diffusers (GLSDs) are optical components for use in special-purpose illumination systems (see figure). When properly positioned with respect to lamps and areas to be illuminated, holographic GLSDs efficiently channel light from the lamps onto specified areas with specified distributions of illumination for example, uniform or nearly uniform irradiance can be concentrated with intensity confined to a peak a few degrees wide about normal incidence, over a circular or elliptical area. Holographic light diffusers were developed during the 1990s. The development of the present holographic GLSDs extends the prior development to incorporate sol-gel optical glass. To fabricate a holographic GLSD, one records a hologram on a sol-gel silica film formulated specially for this purpose. The hologram is a quasi-random, micro-sculpted pattern of smoothly varying changes in the index of refraction of the glass. The structures in this pattern act as an array of numerous miniature lenses that refract light passing through the GLSD, such that the transmitted light beam exhibits a precisely tailored energy distribution. In comparison with other light diffusers, holographic GLSDs function with remarkably high efficiency: they typically transmit 90 percent or more of the incident lamp light onto the designated areas. In addition, they can withstand temperatures in excess of 1,000 C. These characteristics make holographic GLSDs attractive for use in diverse lighting applications that involve high temperatures and/or requirements for high transmission efficiency for ultraviolet, visible, and near-infrared light. Examples include projectors, automobile headlights, aircraft landing lights, high-power laser illuminators, and industrial and scientific illuminators.

Yu, Kevin↗
NASA NTRS1989

The numerical solution of the Navier-Stokes equations for 3-dimensional, unsteady, incompressible flows by compact schemes

The present numerical method for the solution of unsteady, incompressible three-dimensional flow Navier-Stokes equations using velocity-vorticity variables and irregular Cartesian grids proceeds by solving: (1) equations of Cauchy-Riemann type for the velocity; and (2) transport-diffusion equations for the vorticity, whose solenoidal vorticity components are generated by a Poisson equation for an appropriate scalar potential. Iterations are used to solve the finite difference equations, facilitating the use of vector and parallel-computing methods; numerical experiments have verified the method's second-order spatial and temporal accuracy.

Gatski, T. B.↗
NASA NTRS1972

Mass diffusion in a self-confined rotating flow

The effectiveness of fluid containment near an interior stagnation point and within a self-confined stagnation region is determined by numerically solving the species conservation equation for a bi-component mixture. The flow geometry is that of a swirling fluid stream containing a stationary eddy on the axis of rotation. The base flow is axisymmetric, and the Reynolds number is equal to 50. Schmidt numbers range from 0.1 to 10.

Torrance, K. E.↗