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 415 records · Page 23

Large Eddy Simulation of Gravitational Effects on Transitional and Turbulent Gas-Jet Diffusion Flames

The basic objective of this work is to assess the influence of gravity on "the compositional and the spatial structures" of transitional and turbulent diffusion flames via large eddy simulation (LES), and direct numerical simulation (DNS). The DNS is conducted for appraisal of the various closures employed in LES, and to study the effect of buoyancy on the small scale flow features. The LES is based on our "filtered mass density function"' (FMDF) model. The novelty of the methodology is that it allows for reliable simulations with inclusion of "realistic physics." It also allows for detailed analysis of the unsteady large scale flow evolution and compositional flame structure which is not usually possible via Reynolds averaged simulations.

Givi, Peyman↗

Some Experiences with Nonoverlapping Schur Complement Parallel Preconditioning for CFD Calculations

In this work we consider solving matrices which arise from the discretization of advection-diffusion field equations on arbitrary triangulated domains using stabilized numerical methods. The talk will discuss several candidate matrix preconditioning algorithms based on the 2 x 2 block factorization induced by an apriori partitioning of the triangulated domain. Application of the 2 x 2 block preconditioner requires the formation and inversion of the Schur complement submatrix. We consider several strategies for simplifying this task: incomplete Schur complement factorizations, drop tolerance element filling, Schur complement probing, and localized Schur complement inversion. Numerical results will be shown comparing performance and efficiency of these approximations. The matrix preconditioner has also been embedded into a Newton algorithm for solving the nonlinear Euler and Navier-Stokes equations governing compressible flow. The remainder of the talk will show numerous examples in CFD to demonstrate the efficiency and robustness of the techniques.

Barth, Timothy J.↗

Low Stretch Solid-Fuel Flame Transient Response to a Step Change in Gravity

The effect of a step change in gravity level on the stability of low stretch diffusion flames over a solid fuel is studied both numerically and experimentally. Drop tower experiments have been conducted in NASA Glenn Research Center's 5.2 Zero Gravity Facility. In the experiments burning PMMA cylinders, a dynamic transition is observed when the steadily burning 1g flame is dropped and becomes a 0g flame. To understand the physics behind this dynamic transition, a transient stagnation point model has been developed which includes gas-phase radiation and solid phase coupling to describe this dynamic process. In this paper, the experimental results are compared with the model predictions. Both model and experiment show that the interior of the solid phase does not have time to change significantly in the few seconds of drop time, so the experimental results are pseudo-steady in the gas-phase, but the solid is inherently unsteady over long time scales. The model is also used to examine the importance of fractional heat losses on extinction, which clearly demonstrates that as the feedback from the flame decreases, the importance of the ongoing heat losses becomes greater, and extinction is observed when these losses represent 80% or more of the flame feedback.

Armstrong, J. B.↗

Multiple Scattering of Laser Pulses in Snow Over Ice: Modeling the Potential Bias in ICESat Altimetry

The primary goal of NASA's current ICESat and future ICESat2 missions is to map the altitude of the Earth's land ice with high accuracy using laser altimetry technology, and to measure sea ice freeboard. Ice however is a highly transparent optical medium with variable scattering and absorption properties. Moreover, it is often covered by a layer of snow with varying depth and optical properties largely dependent on its age. We describe a modeling framework for estimating the potential altimetry bias caused by multiple scattering in the layered medium. We use both a Monte Carlo technique and an analytical diffusion model valid for optically thick media. Our preliminary numerical results are consistent with estimates of the multiple scattering delay from laboratory measurements using snow harvested in Greenland, namely, a few cm. Planned refinements of the models are described.

Davis, A. B.↗

An improvement in the numerical integration procedure used in the NASA Marshall engineering thermosphere model

A proposed replacement scheme for the integration of the barometric and diffusion equations in the NASA Marshall Engineering Thermosphere (MET) model is presented. This proposed integration scheme is based on Gaussian Quadrature. Extensive numerical testing reveals it to be faster, more accurate and more reliable than the present integration scheme (a modified form of Simpson's Rule) used in the MET model. Numerous graphical examples are provided, along with a listing of a modified form of the MET model in which subroutine INTEGRATE (using Simpson's Rule) is replaced by subroutine GAUSS (which uses Gaussian Quadrature). It is recommended that the Gaussian Quadrature integration scheme, as used here, be used in the MET model.

Hickey, Michael Philip↗

Point source contributions to the extreme ultraviolet background

The results of the Apollo-Soyuz extreme ultraviolet (EUV)/soft X-ray diffuse background survey are analyzed to assess the contribution of discrete-sources with blackbody or coronal plasma spectra. Blackbody emitters in the EUV must be less luminous and more numerous than hot white dwarfs in order to produce the bulk of the diffuse flux. Significant contributions to the EUV background appear likely only if sufficient numbers of white dwarfs possess low temperature coronae similar to the one suggested for Sirius B.

Stern, R.↗

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

Gamma-ray spectroscopy: The diffuse galactic glow

The goal of this project is the development of a numerical code that provides statistical models of the sky distribution of gamma-ray lines due to the production of radioactive isotopes by ongoing Galactic nucleosynthesis. We are particularly interested in quasi-steady emission from novae, supernovae, and stellar winds, but continuum radiation and transient sources must also be considered. We have made significant progress during the first half period of this project and expect the timely completion of a code that can be applied to Oriented Scintillation Spectrometer Experiment (OSSE) Galactic plane survey data.

Hartmann, Dieter H.↗

Analysis of redox additive-based overcharge protection for rechargeable lithium batteries

The overcharge condition in secondary lithium batteries employing redox additives for overcharge protection, has been theoretically analyzed in terms of a finite linear diffusion model. The analysis leads to expressions relating the steady-state overcharge current density and cell voltage to the concentration, diffusion coefficient, standard reduction potential of the redox couple, and interelectrode distance. The model permits the estimation of the maximum permissible overcharge rate for any chosen set of system conditions. Digital simulation of the overcharge experiment leads to numerical representation of the potential transients, and estimate of the influence of diffusion coefficient and interelectrode distance on the transient attainment of the steady state during overcharge. The model has been experimentally verified using 1,1-prime-dimethyl ferrocene as a redox additive. The analysis of the experimental results in terms of the theory allows the calculation of the diffusion coefficient and the formal potential of the redox couple. The model and the theoretical results may be exploited in the design and optimization of overcharge protection by the redox additive approach.

Narayanan, S. R.↗

Shapes of Buoyant and Nonbuoyant Methane Laminar Jet Diffusion Flames

Laminar gas jet diffusion flames represent a fundamental combustion configuration. Their study has contributed to numerous advances in combustion, including the development of analytical and computational combustion tools. Laminar jet flames are pertinent also to turbulent flames by use of the laminar flamelet concept. Investigations into the shapes of noncoflowing microgravity laminar jet diffusion flames have primarily been pursued in the NASA Lewis 2.2-second drop tower, by Cochran and coworkers and by Bahadori and coworkers. These studies were generally conducted at atmospheric pressure; they involved soot-containing flames and reported luminosity lengths and widths instead of the flame-sheet dimensions which are of Greater value to theory evaluation and development. The seminal model of laminar diffusion flames is that of Burke and Schumann, who solved the conservation of momentum equation for a jet flame in a coflowing ambient by assuming the velocity of fuel, oxidizer and products to be constant throughout. Roper and coworkers improved upon this model by allowing for axial variations of velocity and found flame shape to be independent of coflow velocity. Roper's suggestion that flame height should be independent of gravity level is not supported by past or present observations. Other models have been presented by Klajn and Oppenheim, Markstein and De Ris, Villermaux and Durox, and Li et al. The common result of all these models (except in the buoyant regime) is that flame height is proportional to fuel mass flowrate, with flame width proving much more difficult to predict. Most existing flame models have been compared with shapes of flames containing soot, which is known to obscure the weak blue emission of flame sheets. The present work involves measurements of laminar gas jet diffusion flame shapes. Flame images have been obtained for buoyant and nonbuoyant methane flames burning in quiescent air at various fuel flow-rates, burner diameters and ambient pressures. Soot concentrations were minimized by selecting conditions at low flowrates and low ambient pressures; this allows identification of actual flame sheets associated with blue emissions of CH and CO2. The present modeling effort follows that of Roper and is useful in explaining many of the trends observed.

Sunderland, Peter B.↗

Investigation of the flow in the diffuser section of the NASA Lewis icing research tunnel

The flow in the diffuser section of the Icing Research Wind Tunnel at NASA Lewis Research Center is investigated using both tunnel calibration measurements and numerical simulation techniques. Local pressure and temperature measurements are made to establish velocity and temperature profiles in the diffuser of the tunnel. These profiles are compared with similar measurements made prior to renovating the equipment which generates the tunnel's icing cloud. This comparison indicates the manner in which this change affected the flow. The measured data were also compared with a numerical simulation of the flow to help understand how such changes may favorably alter the tunnel flow.

Addy, Harold E., Jr.↗

Investigation of the flow in the diffuser section of the NASA Lewis Icing Research Tunnel

The flow in the diffuser section of the Icing Research Wind Tunnel at NASA Lewis Research Center is investigated using both tunnel calibration measurements and numerical simulation techniques. Local pressure and temperature measurements are made to establish velocity and temperature profiles in the diffuser of the tunnel. These profiles are compared with similar measurements made prior to renovating the equipment which generates the tunnel's icing cloud. This comparison indicates the manner in which this change affected the flow. The measured data were also compared with a numerical simulation of the flow to help understand how such changes may favorably alter the tunnel flow.

Addy, Harold E., Jr.↗

Inlet Flow Control and Prediction Technologies for Embedded Propulsion Systems

Fail-safe, hybrid, flow control (HFC) is a promising technology for meeting high-speed cruise efficiency, low-noise signature, and reduced fuel-burn goals for future, Hybrid-Wing-Body (HWB) aircraft with embedded engines. This report details the development of HFC technology that enables improved inlet performance in HWB vehicles with highly integrated inlets and embedded engines without adversely affecting vehicle performance. In addition, new test techniques for evaluating Boundary-Layer-Ingesting (BLI)-inlet flow-control technologies developed and demonstrated through this program are documented, including the ability to generate a BLI-like inlet-entrance flow in a direct-connect, wind-tunnel facility, as well as, the use of D-optimal, statistically designed experiments to optimize test efficiency and enable interpretation of results. Validated improvements in numerical analysis tools and methods accomplished through this program are also documented, including Reynolds-Averaged Navier-Stokes CFD simulations of steady-state flow physics for baseline, BLI-inlet diffuser flow, as well as, that created by flow-control devices. Finally, numerical methods were employed in a ground-breaking attempt to directly simulate dynamic distortion. The advances in inlet technologies and prediction tools will help to meet and exceed "N+2" project goals for future HWB aircraft.

McMillan, Michelle L.↗

Fire Suppression in Low Gravity Using a Cup Burner

Longer duration missions to the moon, to Mars, and on the International Space Station increase the likelihood of accidental fires. The goal of the present investigation is to: (1) understand the physical and chemical processes of fire suppression in various gravity and O2 levels simulating spacecraft, Mars, and moon missions; (2) provide rigorous testing of numerical models, which include detailed combustion-suppression chemistry and radiation sub-models; and (3) provide basic research results useful for advances in space fire safety technology, including new fire-extinguishing agents and approaches.The structure and extinguishment of enclosed, laminar, methane-air co-flow diffusion flames formed on a cup burner have been studied experimentally and numerically using various fire-extinguishing agents (CO2, N2, He, Ar, CF3H, and Fe(CO)5). The experiments involve both 1g laboratory testing and low-g testing (in drop towers and the KC-135 aircraft). The computation uses a direct numerical simulation with detailed chemistry and radiative heat-loss models. An agent was introduced into a low-speed coflowing oxidizing stream until extinguishment occurred under a fixed minimal fuel velocity, and thus, the extinguishing agent concentrations were determined. The extinguishment of cup-burner flames, which resemble real fires, occurred via a blowoff process (in which the flame base drifted downstream) rather than the global extinction phenomenon typical of counterflow diffusion flames. The computation revealed that the peak reactivity spot (the reaction kernel) formed in the flame base was responsible for attachment and blowoff of the trailing diffusion flame. Furthermore, the buoyancy-induced flame flickering in 1g and thermal and transport properties of the agents affected the flame extinguishment limits.

Takahashi, Fumiaki↗

Fire Suppression in Low Gravity Using a Cup Burner

Longer duration missions to the moon, to Mars, and on the International Space Station increase the likelihood of accidental fires. The goal of the present investigation is to: (1) understand the physical and chemical processes of fire suppression in various gravity and O2 levels simulating spacecraft, Mars, and moon missions; (2) provide rigorous testing of numerical models, which include detailed combustion suppression chemistry and radiation sub-models; and (3) provide basic research results useful for advances in space fire safety technology, including new fire-extinguishing agents and approaches. The structure and extinguishment of enclosed, laminar, methane-air co-flow diffusion flames formed on a cup burner have been studied experimentally and numerically using various fire-extinguishing agents (CO2, N2, He, Ar, CF3H, and Fe(CO)5). The experiments involve both 1g laboratory testing and low-g testing (in drop towers and the KC-135 aircraft). The computation uses a direct numerical simulation with detailed chemistry and radiative heat-loss models. An agent was introduced into a low-speed coflowing oxidizing stream until extinguishment occurred under a fixed minimal fuel velocity, and thus, the extinguishing agent concentrations were determined. The extinguishment of cup-burner flames, which resemble real fires, occurred via a blowoff process (in which the flame base drifted downstream) rather than the global extinction phenomenon typical of counterflow diffusion flames. The computation revealed that the peak reactivity spot (the reaction kernel) formed in the flame base was responsible for attachment and blowoff of the trailing diffusion flame. Furthermore, the buoyancy-induced flame flickering in 1g and thermal and transport properties of the agents affected the flame extinguishment limits.

Takahashi, Fumiaki↗

Variations of the cosmic-ray flux with time

The physical foundations of the modulation of the galactic cosmic ray flux by the sun are reviewed and related to heliospheric structure and dynamics. The basic physical effects - diffusion, convection, adiabatic cooling and drifts - are evaluated and shown to be all important. The results of numerical models are briefly presented and compared with observations. Present-day modulation is shown to reflect the combination of diffusion-convection effects and drift effects. The Maunder minimum is conjectured to be a period when a quiet sun and smooth solar wind resulted in the dominance of drift effects.

Jokipii, J. R.↗

Binary Oscillatory Crossflow Electrophoresis

Electrophoresis has long been recognized as an effective analytic technique for the separation of proteins and other charged species, however attempts at scaling up to accommodate commercial volumes have met with limited success. In this report we describe a novel electrophoretic separation technique - Binary Oscillatory Crossflow Electrophoresis (BOCE). Numerical simulations indicate that the technique has the potential for preparative scale throughputs with high resolution, while simultaneously avoiding many problems common to conventional electrophoresis. The technique utilizes the interaction of an oscillatory electric field and a transverse oscillatory shear flow to create an active binary filter for the separation of charged protein species. An oscillatory electric field is applied across the narrow gap of a rectangular channel inducing a periodic motion of charged protein species. The amplitude of this motion depends on the dimensionless electrophoretic mobility, alpha = E(sub o)mu/(omega)d, where E(sub o) is the amplitude of the electric field oscillations, mu is the dimensional mobility, omega is the angular frequency of oscillation and d is the channel gap width. An oscillatory shear flow is induced along the length of the channel resulting in the separation of species with different mobilities. We present a model that predicts the oscillatory behavior of charged species and allows estimation of both the magnitude of the induced convective velocity and the effective diffusivity as a function of a in infinitely long channels. Numerical results indicate that in addition to the mobility dependence, the steady state behavior of solute species may be strongly affected by oscillating fluid into and out of the active electric field region at the ends of the cell. The effect is most pronounced using time dependent shear flows of the same frequency (cos((omega)t)) flow mode) as the electric field oscillations. Under such conditions, experiments indicate that solute is drawn into the cell from reservoirs at both ends of the cell leading to a large mass build up. As a consequence, any initially induced mass flux will vanish after short times. This effect was not captured by the infinite channel model and hence numerical and experimental results deviated significantly. The revised model including finite cell lengths and reservoir volumes allowed quantitative predictions of the time history of the concentration profile throughout the system. This latter model accurately describes the fluxes observed for both oscillatory flow modes in experiments using single protein species. Based on the results obtained from research funded under NASA grant NAG-8-1080.S, we conclude that binary separations are not possible using purely oscillatory flow modes because of end effects associated with the cos((omega)t) mode. Our research shows, however, that a combination of cos(2(omega)t) and steady flow should lead to efficient separation free of end effects. This possibility is currently under investigation.

Molloy, Richard F.↗

Short efficient ejector systems

A research program was conducted to improve the performance of low pressure ratio ejector systems. The results show that short, efficient ejectors operating at nearly ideal performance are possible through the use of forced mixer lobes. Forced mixer lobes generate large scale axial vorticity which results in rapid mixing and improved diffuser performance. Ejector testing was conducted using both an ejector wind tunnel. Numerous mixer lobe variations were tested to develop lobe design guidelines. The improved performance of mixer-ejectors is presented over a range of operating conditions. Results of mixer lobe angle, penetration, and alignment are presented. Lobe angles of up to 25 deg, coupled with diffuser wall angles over 20 deg, operate without separation, allowing efficient, short ejector systems. Both warm and cold flow test results are presented. Temperature similarity expressions previously developed are further formulated and verified.

Presz, Walter M., Jr.↗