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.

SourceResearch connectionFrom yearThrough year

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

At least 811 records · Page 45

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Dissipation of ionospheric irregularities by wave-particle and collisional interactions

The nonlinear dissipation of plasma irregularities aligned parallel to an ambient magnetic field is studied numerically using a model which employs both wave-particle and collisional diffusion. A wave-particle diffusion coefficient derived from a local theory of the universal drift instability is used. This coefficient is effective in regions of nonzero plasma gradients and produces triangular-shaped irregularities with spectra which vary as f to the -4th, where f is the spatial frequency. Collisional diffusion acts rapidly on the vertices of the irregularities to reduce their amplitude. The simultaneous action of the two dissipative processes is more efficient than collisions acting alone. In this model, wave-particle diffusion mimics the forward cascade process of wave-wave coupling.

Bernhardt, P. A.↗

NASA NTRS1982

Magnetograph response to canopy-type fields

The response of longitudinal-field magnetographs to magnetic fields which are semi-infinite or confined to a horizontal layer is discussed with respect to the interpretation of solar diffuse fields, observed towards the limb, in terms of magnetic canopy models. Numerical results are presented for several reference solar models and typical 'calibration' curves are shown for the CI 9111 A, Fe I 8688 A, and Ca II 8542 A lines in magnetostatic atmospheres derived from a mean model. A procedure is developed for determining the base heights of magnetic canopies from observations with an uncertainty not exceeding the order of a pressure scale height. Until definitive information regarding atmospheric structure inside flux tubes can be developed from theory or observation, reliable field strengths cannot be derived from the data.

Jones, H. P.↗

NASA NTRS1986

Applications of a new wall function to turbulent flow computations

A new wall function approach is developed based on a wall law suitable for incompressible turbulent boundary layers under strong adverse pressure gradients. This wall law was derived from a one-dimensional analysis of the turbulent kinetic energy equation with gradient diffusion concept employed in modeling the near-wall shear stress gradient. Numerical testing cases for the present wall functions include turbulent separating flows around an airfoil and turbulent recirculating flows in several confined regions. Improvements on the predictions using the present wall functions are illustrated. For cases of internal recirculating flows, one modification factor for improving the performance of the k-epsilon turbulence model in the flow recirculation regions is also included.

Chen, Y. S.↗

NASA NTRS1988

Effect of heat release on the spatial stability of a supersonic reacting mixing layer

A numerical study of the stability of compressible mixing layers in which a diffusion flame is embedded is described. The mean velocity profile has been approximated by a hyperbolic tangent profile and the limit of infinite activation energy taken, which reduces the diffusion flame to a flame sheet. The addition of combustion in the form of a flame sheet was found to have important, and complex, effects on the flow stability.

Jackson, T. L.↗

NASA NTRS1989

Ballistic transport in planetary ring systems due to particle erosion mechanisms. I - Theory, numerical methods, and illustrative examples

Ballistic transport, defined as the net radial transport of mass and angular momentum due to exchanges of meteoroid hypersonic-impact ejecta by neighboring planetary ring regions on time-scales orders-of-magnitude shorter than the age of the solar system, is presently considered as a problem in mathematical physics. The preliminary results of a numerical scheme for following the combined effects of ballistic transport and viscous diffusion demonstrate that ballistic transport generates structure near sharp edges already present in the ring-mass distribution; the entire ring system ultimately develops an undulatory structure whose length scale is typically of the order of the radial excursion of the impact ejecta.

Durisen, Richard H.↗

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 NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that removes the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional flow over a blunt biconic.

Turkel, E.↗

NASA NTRS1991

Multigrid for hypersonic viscous two- and three-dimensional flows

The use of a multigrid method with central differencing to solve the Navier-Stokes equations for hypersonic flows is considered. The time-dependent form of the equations is integrated with an explicit Runge-Kutta scheme accelerated by local time stepping and implicit residual smoothing. Variable coefficients are developed for the implicit process that remove the diffusion limit on the time step, producing significant improvement in convergence. A numerical dissipation formulation that provides good shock-capturing capability for hypersonic flows is presented. This formulation is shown to be a crucial aspect of the multigrid method. Solutions are given for two-dimensional viscous flow over a NACA 0012 airfoil and three-dimensional viscous flow over a blunt biconic.

Turkel, E.↗

NASA NTRS1990

Bridgman crystal growth in low gravity - A scaling analysis

The results of an order-of-magnitude or scaling analysis are compared with those of numerical simulations of the effects of steady low gravity on compositional nonuniformity in crystals grown by the Bridgman-Stockbarger technique. In particular, the results are examined of numerical simulations of the effect of steady residual acceleration on the transport of solute in a gallium-doped germanium melt during directional solidification under low-gravity conditions. The results are interpreted in terms of the relevant dimensionless groups associated with the process, and scaling techniques are evaluated by comparing their predictions with the numerical results. It is demonstrated that, when convective transport is comparable with diffusive transport, some specific knowledge of the behavior of the system is required before scaling arguments can be used to make reasonable predictions.

Alexander, J. I. D.↗

NASA NTRS1993

A flux-limited treatment for the conductive evaporation of spherical interstellar gas clouds

In this work, we present and analyze a new analytic solution for the saturated (flux-limited) thermal evaporation of a spherical cloud. This work is distinguished from earlier analytic studies by allowing the thermal conductivity to change continuously from a diffusive to a saturated form, in a manner usually employed only in numerical calculations. This closed form solution will be of interest as a computational benchmark. Using our calculated temperature profiles and mass-loss rates, we model the thermal evaporation of such a cloud under typical interstellar medium (ISM) conditions, with some restrictions. We examine the ionization structure of the cloud-ISM interface and evaluate column densities of carbon, nitrogen, oxygen, neon, and silicon ions toward the cloud. In accord with other investigations, we find that ionization equilibrium is far from satisfied under the assumed conditions. Since the inclusion of saturation effects in the heat flux narrows the thermal interface relative to its classical structure, we also find that saturation effects tend to lower predicted column densities.

Dalton, William W.↗

NASA NTRS1994

Optical observations of NGC 2915: A nearby blue compact dwarf galaxy

This paper presents B and R band Charge Coupled Device (CCD) images and medium resolution spectroscopy of NGC 2915, a relatively isolated BCD (blue compact dwarf) galaxy at a distance of approximately 5 Mpc. NGC 2915 contains two stellar populations: a high surface brightness blue core population and a red diffuse population. The core population contains all of the H II, and numerous embedded objects. It is the locus of current high mass star formation. The brightest embedded objects are likely to be young ionizing clusters, while many of the fainter objects are likely to be individual supergiant stars with masses up to approximately 25 solar mass, or blends of a few such stars. Curious aligned structures on the SE side of the galaxy are seen and their nature discussed. The spectrum of the core is dominated by bright narrow emission lines like that of a high excitation and low metallicity (less than half solar) H II region. The continuum is flat, with Balmer and Ca II features seen in absorption. The velocity of the Ca II features suggest contamination by galactic interstellar absorption. There is a significant velocity gradient in the spectra, probably indicative of rotation. Outside of its core, NGC 2915 resembles a dE (dwarf elliptical) galaxy, in that it has an exponential surface brightness profile, is red ((B-R)(sub 0) = 1.65), and has a low extrapolated central surface brightness (B(0)(sub c) = 22.44). NGC 2915's properties are compared with other BCDs, concentrating on two morphologically similar BCDs that are near enough to resolve into stars: NGC 1705 and NGC 5253. It is noted that the presence of winds in BCDs invalidates closed box chemical evolution models and the remaining constraints on star formation duration are relatively weak. Some BCDs, including NGC 2915, may be able to maintain their present star formation rate for Gyr time scales. This suggests that the overall evolution of these BCDs may be much slower than the approximately 10 Myr burst time scales commonly quoted. However, shortly after the formation of a massive (10(exp 6) solar mass) cluster a BCD will have all the properties of strong starburst galaxy).

Meurer, G. R.↗

NASA NTRSDate not supplied

Vapor Growth of ZnSe-based Compound Semiconductors in Low Gravity environment on International Space Station

Six low-gravity experiments of crystal growth of ZnSe-based compound semiconductors by physical vapor transport (PVT) are scheduled to be processed inside Low Gradient Furnace (LGF) on Materials Science Research Rack (MSRR) in Material Science Laboratory (MSL) on International Space Station (ISS) starting from February 27th, 2023. The flight samples consist of three unseeded growths of ZnSe, Cr-doped ZnSe and ZnSeTe as well as growth of ZnSe on top of three single crystal ZnSe seeds with orientation of (100), (110) and (111). The main scientific objective of the flight experiments is to establish the effects of gravity-driven fluid flows on Earth to (i) the non-uniform incorporation of impurities and defects and (ii) the deviation from stoichiometry observed in the grown crystals caused by growth interface fluctuations and evolution as the results of buoyancy-driven convection and irregular fluid-flows. The investigation will also evaluate the effects of gravity on the PVT process by examining (i) the growth kinetics on various seed orientations, (ii) the dopant segregation and distribution in the Cr doped ZnSe, and (iii) the compositional segregation and distribution in the ternary ZnSeTe crystal grown by PVT. Additionally, from X-ray diffraction/topography and by measuring the density and distribution of various crystalline defects it will assess the self-induced strain effects developed during processing at elevated temperatures caused by the weight of the crystals. During the ground-based preparations, the following tasks have been accomplished: • Establish the thermodynamic properties of the ZnSe binary and related ternary systems by measuring the partial pressures of the vapor phase coexisting with the condensed phases and perform phase diagram calculations for comparison. • Evaluate the effects of gravity-driven convection in the growth process by performing experiments under various vapor transport orientations relative to gravity vector, i.e., horizontal, vertically stabilized and destabilized configurations. • Establish a quantitative correlation between growth parameters and the characteristics and properties of the grown crystals by characterizing the crystals using a variety of techniques. • Evaluate the fundamentals of the current theories on vapor transport and crystal growth kinetics by performing in-situ and real-time monitoring techniques during growth. • Develop a one-dimensional diffusion limited theory of mass transport and two- and three-dimensional numerical simulations of both mass transport and heat transfer and thus establish the microgravity requirements for the flight experiments and to predict the characteristics of the grown crystals such as the solid-vapor interface shapes… etc. and establish a fundamental understanding of the crystal growth process from the comparison with experimental results.

ZnSe↗

DOE OSTI2019

Practical Challenges and Failure Modes During Fabrication of Haynes 230 Micro-Pin Solar Receivers for High Temperature Supercritical Carbon Dioxide Operation

Printed circuit heat exchangers (PCHX) fabricated from high nickel alloys have shown promise as primary heat exchangers, recuperators, and solar receivers for high temperature and pressure supercritical carbon dioxide power cycles. There are numerous challenges in fabricating these devices including forming channel features, joining via diffusion or transient liquid phase bonding, and brazing/welding of headers. Commercial entities are understandably hesitant to share propriety best practices, and both commercial and non-commercial entities working on these components tend not to publically share failures and mistakes. However, sharing of this information could prevent similar costly errors and further the understanding of the coupled manufacturing, materials and mechanics issues in creating these components. Thus, in this paper, we document the challenges, failures, and mitigation methods uncovered in fabrication of prototype micro-pin based solar thermal receivers constructed of Haynes 230 and fabricated through a combination of photochemical machining (PCM), wire electrical discharge machining (EDM), transient liquid phase (TLP) bonding, vacuum brazing, and gas tungsten arc welding (GTAW). The receiver is designed to absorb concentrated solar fluxes greater than 140 W cm-2, while heating supercritical carbon dioxide from 550 °C to 720 °C at a pressure of 20 MPa to 25 MPa. The prototype receiver consists of a thin (~ 450 µm), Haynes 230 coversheet bonded to a 15 cm × 15 cm Haynes 230 micro-pin plate. The pieces are joined using a TLP bonding process with a nickel-phosphorus interlayer. Prior to bonding, micro-pins with height ~150 µm and diameter 300 µm are fabricated using PCM in the plate, and through slot features are made using wire EDM. Finally, flow headers are joined to the microchannel plate through a combination of vacuum brazing and GTAW. During hydrostatic proof testing, the prototype device failed when the coversheet delaminated from the pin array at a pressure of 290 bar. A failure analysis including scanning electron microscopy (SEM) to view failure sites and energy-dispersive X-ray spectroscopy (EDS) to evaluate elemental analysis of the failed areas was conducted. The failure modes can be broadly categorized as (1) failures potentially relating to reliquification of the transient liquid phase bonds between the micro-pins plate and coversheet during post-processing, (2) failures related to manufacturing defects, and (3) failures attributed to design.

14 SOLAR ENERGY↗

NASA NTRS1991

Examination of various turbulence models for application in liquid rocket thrust chambers

There is a large variety of turbulence models available. These models include direct numerical simulation, large eddy simulation, Reynolds stress/flux model, zero equation model, one equation model, two equation k-epsilon model, multiple-scale model, etc. Each turbulence model contains different physical assumptions and requirements. The natures of turbulence are randomness, irregularity, diffusivity and dissipation. The capabilities of the turbulence models, including physical strength, weakness, limitations, as well as numerical and computational considerations, are reviewed. Recommendations are made for the potential application of a turbulence model in thrust chamber and performance prediction programs. The full Reynolds stress model is recommended. In a workshop, specifically called for the assessment of turbulence models for applications in liquid rocket thrust chambers, most of the experts present were also in favor of the recommendation of the Reynolds stress model.

Hung, R. J.↗

DOE OSTI2024

Thermal diffusion, exhaust gas recirculation and blending effects on lean premixed hydrogen flames

Thermodiffusively-unstable lean premixed hydrogen flames are investigated using two-dimensional direct numerical simulation employing finite-rate chemical kinetics. Three databases are generated focussing on the inclusion of the Soret effect, the recirculation of exhaust gas, and blending with methane. A simple rescaling of a classic thermal diffusion model is presented and shown to mimic multicomponent diffusion with very low computational cost and little-to-no loss in accuracy. It is also shown that a previously developed model for mean local flame speeds in lean premixed hydrogen flames can still be used provided Soret effects are taken into account in one-dimensional calculations. The addition of exhaust gas to the unburned mixture is found to enhance thermodiffusive instability; the primary mechanism for this was shown to be the highly-efficient third-body nature of water, with the reduction of adiabatic flame temperature a second-order effect. Again, the existing mean local flame speed model proved sufficient. Finally, blending with methane was found to reduce the thermodiffusive response of the flame, more so than the existing model suggests, despite adjustment of the fuel Lewis number; an adapted model is presented to account for this.

08 HYDROGEN↗

DOE OSTI2020

Lyapunov and Clearing Timescales in Planetary Chaotic Zones

Debris disks in planetary systems are known to possess gaps “opened” by planets in their close-to-coorbital neighborhoods; the gaps are free from low-mass material such as planetesimals. We consider the escape process of the initially zero-eccentric particles from the gap zone in a scenario consisting of two basic consequent stages: (1) the particle orbit's eccentricity is inflated, while its specific energy is relatively constant, and (2) the particle orbit's semimajor axis is inflated, while its specific angular momentum is relatively constant. During the first stage, the chaotic transport (diffusion) proceeds along the “staircase” of (p + 1):p (with integer p ≫ 1) overlapping particle–planet resonances, and during the second stage, it proceeds along the p:1 particle–planet resonance staircase. A specialized two-dimensional area-preserving map (a “tokamap” version) is derived to describe the first stage, while a multiharmonic Kepler map theory is used to describe the second stage. We derive the μ (mass parameter) scaling for the clearing timescale in the planetary chaotic zone. At small μ, it turns out to be inverse-quadratic in μ. The μ scaling for the Lyapunov timescale inside the gap is also derived. The derived scalings are in qualitative agreement with available numerical-experimental data.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1974

Effect of gravity on methane-air combustion

Analytical and numerical techniques dealing with the theoretical description of the influence of zero and reduced gravitational acceleration on diffusion flames, with a view to improving understanding of fires in space vehicles, were developed in support of experimental work performed in this area. This was done in order to confirm qualitative understanding of the process, to determine the quantitative accuracy of numerical predictions, and to establish a mathematical model of the process for subsequent use as a predictive and exploratory tool. The following results were accomplished: (1) derivation of differential equations and boundary conditions describing the system, (2) details of the computations, using a FORTRAN computer program, for calculating the flow and heat and mass transfer in two dimensions (both steady and unsteady). It was shown that the experimental behavior can be reproduced with fair accuracy, provided that the time step is sufficiently short.

Elghobashi, S.↗

NASA NTRS1998

3-D Modeling of Double-Diffusive Convection During Directional Solidification of a Non-Dilute Alloy with Application to the HgCdTe Growth Under Microgravity Conditions

A numerical calculation for a non-dilute alloy solidification was performed using the FIDAP finite element code. For low growth velocities plane front solidification occurs. The location and the shape of the interface was determined using melting temperatures from the HgCdTe liquidus curve. The low thermal conductivity of the solid HgCdTe causes thermal short circuit through the ampoule walls, resulting in curved isotherms in the vicinity of the interface. Double-diffusive convection in the melt is caused by radial temperature gradients and by material density inversion with temperature. Cooling from below and the rejection at the solid-melt interface of the heavier HgTe-rich solute each tend to reduce convection. Because of these complicating factors dimensional rather then non-dimensional modeling was performed. Estimates of convection contributions for various gravity conditions was performed parametrically. For gravity levels higher then 1 0 -7 of earth's gravity it was found that the maximum convection velocity is extremely sensitive to gravity vector orientation and can be reduced at least by factor of 50% for precise orientation of the ampoule in the microgravity environment. The predicted interface shape is in agreement with one obtained experimentally by quenching. The results of 3-D modeling are compared with previous 2-D finding. A video film featuring melt convection will be presented.

Bune, Andris V.↗