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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Dynamic response of a pulsed Burke-Schumann diffusion flame

Turbulent flames are often envisioned as an ensemble of random vortices interacting with the combustion process. A better understanding of the vortex-flame interactions therefore would be useful in improving the modeling of turbulent diffusion flames. Substantial simplification may be made by investigating controlled interactions in a laminar flame, as opposed to random interactions in a turbulent flame. The general goals of the research project are to improve our understanding of (1) the influence of buoyancy on co-flow diffusion flames and (2) the effects of buoyancy on vortex-flame interactions in co-flow diffusion flames. As a first step toward objective (2), we conducted a joint experimental and numerical investigation of the vortex-flame interaction. Vortices were produced by mechanically pulsing the fuel flow at a low frequency, e.g., 10 Hz. Experiments were conducted using a nonflickering Burke-Schumann flame in both microgravity (mu-g) and normal gravity (1g) as a means of varying the buoyant force without modification of the pressure (i.e., density). The effects of buoyant convection may then be determined by a comparison of the mu-g and 1g results. The mu-g results may also reveal the important mechanisms which are masked or overwhelmed by buoyant convection in 1g. A numerical investigation was conducted using a validated, time-accurate numerical code to study the underlying physics during the flame interaction and to assist the interpretation of the experimental results.

Sheu, Jyh-Cherng↗
NASA NTRS1998

MAMA Detector

Work carried out under this grant led to fundamental discoveries and over one hundred publications in the scientific literature. Fundamental developments in instrumentation were made including all the instrumentation on the EUVE satellite, the invention of a whole new type of grazing instrument spectrometer and the development of fundamentally new photon counting detectors including the Wedge and Strip used on EUVE and many other missions and the Time Delay detector used on OREFUS and FUSE. The Wedge and Strip and Time Delay detectors were developed under this grant for less than two million dollars and have been used in numerous missions most recently for the FUSE mission. In addition, a fundamentally new type of diffuse spectrometer has been developed under this grant which has been used in instrumentation on the MMSAT spacecraft and the Lewis spacecraft. Plans are underway to use this instrumentation on several other missions as well.

Bowyer, Stuart↗
NASA NTRS1999

Pulse Phase Thermography for Defect Detection and Visualization

Pulse Phase Thermography (PPT) has been reported as a novel powerful technique of the thermal NDE. It employs application of the Discrete Fourier Transform (DFT) to thermal images obtained following flash heating of the front surface of a specimen. The computed phasegrams are excellent for defect visualization in a wide range of materials. This is in part due to their low sensitivity to uneven heating. This work is an attempt to analyze advantages and limitations of PPT. Results of application of the DFT to simulated temperature decays are presented. The temperature evolution on a surface has been simulated based on an analytical solution of the one-dimensional heat diffusion problem. A more sophisticated study has been done for different sizes of defects using numerical solution of the three-dimensional mathematical model. Capabilities of PPT for in-depth scanning and for monitoring of the material loss are discussed. The recommendations for the practical application of the PPT are presented. Experimental results obtained following these recommendations are reported.

Marinetti, Sergio↗
NASA NTRS2008

Vibroacoustic Response Data of Stiffened Panels and Cylinders

NASA has collected vibroacoustic response data on a variety of complex, aerospace structures to support research into numerical modeling of such structures. This data is being made available to the modeling community to promote the development and validation of analysis methods for these types of structures. Existing data from two structures is described, as well as plans for a data set from a third structure. The first structure is a 1.22 m by 1.22 m stiffened aluminum panel, typical of a commercial aircraft sidewall section. The second is an enclosed, stiffened aluminum cylinder, approximately 3.66 m long and 1.22 m in diameter, constructed to resemble a small aircraft fuselage with no windows and a periodic structure. The third structure is a filament-wound composite cylinder with composite stiffeners. Numerous combinations of excitation and response variables were measured on the structures, including: shaker excitation; diffuse acoustic field; velocity response from a laser vibrometer; intensity scans; and point acceleration.

Cabell, Randolph↗
NASA NTRS2012

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗
NASA NTRS1996

Kinetics of Diffusional Droplet Growth in a Liquid/Liquid Two-Phase System

We address the problem of diffusional interactions in a finite sized cluster of spherical particles for volume fractions, V(sub v) in the range 0-0.01. We determined the quasi-static monopole diffusion solution for n particles distributed at random in a continuous matrix. A global mass conservation condition is employed, obviating the need for any external boundary condition. The numerical results provide the instantaneous (snapshot) growth or shrinkage rate of each particle, precluding the need for extensive time-dependent computations. The close connection between these snapshot results and the coarsegrained kinetic constants are discussed. A square-root dependence of the deviations of the rate constants from their zero volume fraction value is found for the higher V(sub v) investigated. This behavior is consistent with predictions from diffusion Debye-Huckel screening theory. By contrast, a cube-root dependence, reported in earlier numerical studies, is found for the lower V(sub v) investigated. The roll-over region of the volume fraction where the two asymptotics merge depends on the number of particles, n, alone. A theoretical estimate for the roll-over point predicts that the corresponding V(sub v) varies as n(sup -2), in good agreement with the numerical results.

Glicksman, M. E.↗
NASA NTRS1979

High-altitude gas releases - Transition from collisionless flow to diffusive flow in a nonuniform atmosphere

The paper is motivated by a need for describing high-altitude gas releases from rockets or the Space Shuttle. In the tenuous upper atmosphere, the injected gases expand from a collisionless to a collision-dominated state. The analysis presented extends the efforts of Baum (1973-1974) to include the effects of a nonuniform background atmosphere and a velocity-dependent collision frequency. The unsteady expansion of gas releases is analyzed using gas kinetic theory. The Boltzmann equation with the Krook collision integral is solved numerically. At late times (after many collisions with the background atmosphere) the solution is identical to the one given by diffusion in a nonuniform atmosphere (Bernhardt, 1979). The theoretical model predicts elongation and heating of the vapor trails due to collisions.

Bernhardt, P. A.↗
NASA NTRS2014

High-Order Residual-Distribution Hyperbolic Advection-Diffusion Schemes: 3rd-, 4th-, and 6th-Order

In this paper, spatially high-order Residual-Distribution (RD) schemes using the first-order hyperbolic system method are proposed for general time-dependent advection-diffusion problems. The corresponding second-order time-dependent hyperbolic advection- diffusion scheme was first introduced in [NASA/TM-2014-218175, 2014], where rapid convergences over each physical time step, with typically less than five Newton iterations, were shown. In that method, the time-dependent hyperbolic advection-diffusion system (linear and nonlinear) was discretized by the second-order upwind RD scheme in a unified manner, and the system of implicit-residual-equations was solved efficiently by Newton's method over every physical time step. In this paper, two techniques for the source term discretization are proposed; 1) reformulation of the source terms with their divergence forms, and 2) correction to the trapezoidal rule for the source term discretization. Third-, fourth, and sixth-order RD schemes are then proposed with the above techniques that, relative to the second-order RD scheme, only cost the evaluation of either the first derivative or both the first and the second derivatives of the source terms. A special fourth-order RD scheme is also proposed that is even less computationally expensive than the third-order RD schemes. The second-order Jacobian formulation was used for all the proposed high-order schemes. The numerical results are then presented for both steady and time-dependent linear and nonlinear advection-diffusion problems. It is shown that these newly developed high-order RD schemes are remarkably efficient and capable of producing the solutions and the gradients to the same order of accuracy of the proposed RD schemes with rapid convergence over each physical time step, typically less than ten Newton iterations.

Mazaheri, Alireza R.↗
NASA NTRS1999

The Onset of Magnetic Reconnection in Tail-Like Equilibria

Magnetic reconnection is a fundamental mode of dynamics in the magnetotail, and is recognized as the basic mechanisms converting stored magnetic energy into kinetic energy of plasma particles. The effects of the reconnection process are well documented by spacecraft observations of plasmoids in the distant magnetotail, or bursty bulk flows, and magnetic field dipolarizations in the near Earth region. Theoretical and numerical analyses have, in recent years, shed new light on the way reconnection operates, and, in particular, which microscopic mechanism supports the dissipative electric field in the associated diffusion region. Despite this progress, however. the question of how magnetic reconnection initiates in a tail-like magnetic field with finite flux threading the current i.sheet remains unanswered. Instead, theoretical studies supported by numerical simulations support the point-of-view that such plasma and current sheets are stable with respect to collisionless tearing mode. In this paper, we will further investigate this conclusion, with emphasis on the question whether it remains valid in plasma sheets with embedded thin current sheets. For this purpose, we perform particle-in-cell simulations of the driven formation of thin current sheets, and their subsequent evolution either to equilibrium or to instability of a tearing-type mode. In the latter case we will pay particular attention to the nature of the electric field contribution which unmagnetizes the electrons.

Hesse, Michael↗
NASA NTRS1992

Numerical simulations of time-dependent cosmic ray mediated shocks

Recent research in unsteady shocks that are influenced by nonlinear aspects of diffusive particle acceleration and that involve multidimensional flows is outlined, and some of the difficulties in making further progress are discussed. It is suggested that in many situations it is probably not appropriate to assume steady states for astrophysical CR shocks. This results both from the need to be concerned about the finite rate at which CR pressure builds and modifies a shock and from the existence of instabilities. Even so, it seems likely that CR-modified shocks can develop to a state with substantial CR pressure. The overall acceleration efficiency and the CR particle momentum distribution produced in time-dependent flows typically depend on several details such as the ratio of various time scales and the initial proportion of relativistic to nonrelativistic particles in the CR populations.

Jones, T. W.↗
NASA NTRS2000

Laminar Diffusion Flame Studies (Ground- and Space-Based Studies)

Laminar diffusion flames are of interest because they provide model flame systems that are far more tractable for analysis and experiments than more practical turbulent diffusion flames. Certainly, understanding flame processes within laminar diffusion flames must precede understanding these processes in more complex turbulent diffusion flames. In addition, many properties of laminar diffusion flames are directly relevant to turbulent diffusion flames using laminar flamelet concepts. Laminar jet diffusion flame shapes (luminous flame boundaries) have been of particular interest since the classical study of Burke and Schumann because they are a simple nonintrusive measurement that is convenient for evaluating flame structure predictions. Thus, consideration of laminar flame shapes is undertaken in the following, emphasizing conditions where effects of gravity are small, due to the importance of such conditions to practical applications. Another class of interesting properties of laminar diffusion flames are their laminar soot and smoke point properties (i.e., the flame length, fuel flow rate, characteristic residence time, etc., at the onset of soot appearance in the flame (the soot point) and the onset of soot emissions from the flame (the smoke point)). These are useful observable soot properties of nonpremixed flames because they provide a convenient means to rate several aspects of flame sooting properties: the relative propensity of various fuels to produce soot in flames; the relative effects of fuel structure, fuel dilution, flame temperature and ambient pressure on the soot appearance and emission properties of flames; the relative levels of continuum radiation from soot in flames; and effects of the intrusion of gravity (or buoyant motion) on emissions of soot from flames. An important motivation to define conditions for soot emissions is that observations of laminar jet diffusion flames in critical environments, e.g., space shuttle and space station facilities, cannot involve soot emitting flames in order to ensure that test chamber windows used for experimental observations are not blocked by soot deposits, thereby compromising unusually valuable experimental results. Another important motivation to define conditions where soot is present in diffusion flames is that flame chemistry, transport and radiation properties are vastly simplified when soot is absent, making such flames far more tractable for detailed numerical simulations than corresponding soot-containing flames. Motivated by these observations, the objectives of this phase of the investigation were as follows: (1) Observe flame-sheet shapes (the location of the reaction zone near phi=1) of nonluminous (soot free) laminar jet diffusion flames in both still and coflowing air and use these results to develop simplified models of flame-sheet shapes for these conditions; (2) Observe luminous flame boundaries of luminous (soot-containing) laminar jet diffusion flames in both still and coflowing air and use these results to develop simplified models of luminous flame boundaries for these conditions. In order to fix ideas here, maximum luminous flame boundaries at the laminar smoke point conditions were sought, i.e., luminous flame boundaries at the laminar smoke point; (3) Observe effects of coflow on laminar soot- and smoke-point conditions because coflow has been proposed as a means to control soot emissions and minimize the presence of soot in diffusion flames.

Dai, Z.↗
NASA NTRS1986

Thermofluid analysis of the SSME preburner using a gas-gas diffusion model for oxygen and hydrogen combustion at supercritical pressures

The paper discusses the thermofluid analysis of the Space Shuttle Main Engine (SSME) fuelside preburner. The governing equations have been solved numerically to predict flow, heat transfer, mixing, and combustion. A two-fluid approach is adopted in which oxygen is regarded as one fluid and hydrogen is regarded as the other fluid. The chemical kinetics is assumed to be very fast so that combustion is primarily controlled by the rate of mixing between oxygen and hydrogen. The preburner pressure is much greater than the critical pressures of oxygen and hydrogen; hence, a gas-gas diffusion model (rather than an evaporation model) has been developed to compute the rate of interphase mixing. Empirical correlations have been incorporated to account for the effect of slip on the interphase exchange. A sensitivity study has been performed with various model parameters. It is observed that the model can predict possibility of incomplete combustion and local regions of high temperatures under steady operating conditions. Some of these anomalies have been observed in actual tests, and the numerical model is useful for understanding possible causes and remedies. At least some measurements are needed for quantitative verification of the model.

Prakash, C.↗
NASA NTRS1971

An analysis of the flow field near the fuel injection location in a gas core reactor.

An analytical study is presented which shows the effects of large energy release and the concurrent high acceleration of inner stream fluid on the coaxial flow field in a gas core reactor. The governing equations include the assumptions of only radial radiative transport of energy represented as an energy diffusion term in the Euler equations. The method of integral relations is used to obtain the numerical solution. Results show that the rapidly accelerating, heat generating inner stream actually shrinks in radius as it expands axially.

Weinstein, H.↗
NASA NTRS1974

Analysis of sound propagation in ducts using the wave envelope concept

A finite difference formulation is presented for sound propagation in a rectangular two-dimensional duct without steady flow for plane wave input. Before the difference equations are formulated, the governing Helmholtz equation is first transformed to a form whose solution does not oscillate along the length of the duct. This transformation reduces the required number of grid points by an order of magnitude, and the number of grid points becomes independent of the sound frequency. Physically, the transformed pressure represents the amplitude of the conventional sound wave. Example solutions are presented for sound propagation in a one-dimensional straight hard-wall duct and in a two-dimensional straight soft-wall duct without steady flow. The numerical solutions show evidence of the existence along the duct wall of a developing acoustic pressure diffusion boundary layer which is similar in nature to the conventional viscous flow boundary layer. In order to better illustrate this concept, the wave equation and boundary conditions are written such that the frequency no longer appears explicitly in them. The frequency effects in duct propagation can be visualized solely as an expansion and stretching of the suppressor duct.

Baumeister, K. J.↗
NASA NTRS1980

A simulation of the atmospheric cloud physics laboratory to aid in its design and the design of the experiments within the laboratory

Using the finite difference method with overrelaxation, numerical solutions of the steady-state vorticity transport equation were obtained for a continuous flow diffusion chamber of the Hudson-Squires type. The calculation neglected the effects due to temperature, gravity, and saturation. The size and shape of the manifold used to inject the aerosol laden flow were varied to obtain a design which would improve the performance of the chamber from strictly low Reynolds number (less than 20) fluid dynamical considerations.

Winchester, L. W., Jr.↗
NASA NTRS1984

Nonlinear evolution of the lower-hybrid drift instability

The results of simulations of the lower-hybrid drift instability in a neutral sheet configuration are described. The simulations use an implict formulation to relax the usual time step limitations and thus extend previous explicit calculations to weaker gradients, larger mass ratios, and long times compared with the linear growth time. The numerical results give the scaling of the saturation level, heating rates, resistivity, and cross-field diffusion and a demonstration by comparison with a fluid electron model that dissipation in the lower-hybrid drift instability is caused by electron kinetic effects.

Brackbill, J. U.↗
NASA NTRS1987

Magnetic braking in weakly ionized media

The combined magnetic braking-ambipolar diffusion problem in weakly ionized, rigidly rotating disks is studied. An analytical solution is presented for a disk whose angular velocity and magnetic yield vectors are aligned with the symmetry axis, illustrating the effects of the relative azimuthal drift of neutrals and ions. The effects of radial drift are added, commenting on the ratio of the characteristic ambipolar diffusion and magnetic braking time scales in high-mass and low-mass disks. A numerical calculation is used to show the combined action of these two processes.

Konigl, Arieh↗
NASA NTRS1989

Tau approximation techniques for identification of coefficients in parabolic PDE

A variant of the Tau method, called the weak Tau method, is developed on the basis of the weak form of the PDE for use in least-squares parameter estimation; also presented is a suitable abstract convergence framework. The emphasis is on the theoretical framework that allows treatment of the weak Tau method when it is applied to a wide class of inverse problems, including those for diffusion-advection equations, the Fokker-Planck model for population dynamics, and damped beam equations. Extensive numerical testing of the weak Tau method has demonstrated that it compares quite favorably with existing methods.

Banks, H. T.↗