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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

Thermal analysis of Bridgman-Stockbarger growth

A thermal analysis of a cylindrical HgCdTe sample in a Bridgman-Stockbarger crystal growth configuration was conducted with emphasis on the thermal profile, interface shape and position, and the thermal gradients at the liquid-solid interface. Alloys of HgTe and CdTe with compositions approximating 20 percent CdTe, 80 percent HgTe were used. This composition results in a bandgap suited for the detection of 10.6 micron CO2 radiation. The sensitivity of the sample thermal characteristics to important growth parameters, such as thermal diffusivities, thermal conductivities, furnace temperature profile, ampoule dimensions, and growth velocity was assessed. Numerical techniques and associated computational models necessary to analyze the heat transfer process within the sample and the Bridgman-Stockbarger boundary conditions were developed. This thermal analysis mode was programmed in FORTRAN V, and is currently operational on the MSFC Univac 1100 system.

Knopf, F. W.↗
NASA NTRS1983

Numerical solution of the problem of flame propagation by the use of the random element method

A numerical, grid-free algorithm is presented for one-dimensional reaction-diffusion model of laminar flame propagation in premixed gases. It is based on the random element method we developed for the analysis of diffusional processes. The effect of combustion is taken into account by applying the principle of fractional steps to separate the process of diffusion, modeled by the random walk of computational elements, from the exothermic effects of chemical reaction, monitoring their strength. The validity of the algorithm is demonstrated by application to flame propagation problems for which exact solutions exist. The flame speed evaluated by its use oscillates around the exact value at a relatively small amplitude, while the temperature and species concentration profiles are self-correcting in their convergence to the exact solution. A satisfactory resolution is obtained by the use of quite a small number of computational elements which automatically adjust their distribution of fit sharp gradients.

Ghoniem, A. F.↗
NASA NTRS1985

Effect of free stream turbulence on flow separation

An analysis is presented to study the influence of free stream turbulence on separation in the case of crossflow over a circular cylinder. An eddy diffusivity model has been formulated and the governing momentum equation has been solved numerically. A correlation parameter has been suggested to correlate the angle of separation. It has been found that the angle of separation increases with increasing free stream turbulence intensity.

Gorla, R. S. R.↗
NASA NTRS1990

SINDA '85/FLUINT, Version 2.2

SINDA, Systems Improved Numerical Differencing Analyzer, software system developed for solving physical problems governed by diffusion-type equations modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations, adapted to wide range of problems represented by such differential equations as Fourier, Poisson, or LaPlace. Software system updated. FLUINT, FLUid INTegrator, is advanced, one-dimensional fluid-analysis program solving equations of arbitrary fluid-flow networks. Written in FORTRAN and ASSEMBLER.

Cullimore, Brent↗
NASA NTRS1994

Effect of chemical heat release in a temporally evolving mixing layer

Two-dimensional numerical simulations of a temporally evolving mixing layer with an exothermic infinitely fast diffusion flame between two unmixed reactants have been carried out in the limit of zero Mach number to study the effect of the heat release on the early stages of the evolution of the flow. Attention has been directed to relatively large values of the oxidizer-to-fuel mass stoichiometric ratio typical of hydrocarbon flames, and initial vorticity distributions thicker than the temperature and species distributions have been chosen to mimic the situation at the outlet of a jet. The results show that, during the stages of the evolution covered by the present simulations, enhancement of combustion occurs by local stretching of the flame without much augmentation of its area. The rate of product generation depends strongly on the initial conditions, which suggests the possibility of controlling the combustion by acting on the flow. Rollup and vortex amalgamation still occur in these reacting flows but are very much affected by the production of new vorticity by baroclinic torques. These torques lead to counter rotating vortex pairs around the flame and, more importantly, in thin layers of light fluid that leave the vicinity of the flame when the Kelvin-Helmholtz instability begins to develop. Propelled by the vortex pairs, these layers wind around, split on reaching high pressure regions, and originate new vortex pairs in a process that ends up building large-scale vortices with a vorticity distribution more complex than for a constant density fluid.

Higuera, F. J.↗
NASA NTRS1997

Numerical Study of Stratified Charge Combustion in Wave Rotors

A wave rotor may be used as a pressure-gain combustor effecting non-steady flow, and intermittent, confined combustion to enhance gas turbine engine performance. It will be more compact and probably lighter than an equivalent pressure-exchange wave rotor, yet will have similar thermodynamic and mechanical characteristics. Because the allowable turbine blade temperature limits overall fuel/air ratio to sub-flammable values, premixed stratification techniques are necessary to burn hydrocarbon fuels in small engines with compressor discharge temperature well below autoignition conditions. One-dimensional, unsteady numerical simulations of stratified-charge combustion are performed using an eddy-diffusivity turbulence model and a simple reaction model incorporating a flammability limit temperature. For good combustion efficiency, a stratification strategy is developed which concentrates fuel at the leading and trailing edges of the inlet port. Rotor and exhaust temperature profiles and performance predictions are presented at three representative operating conditions of the engine: full design load, 40% load, and idle. The results indicate that peak local gas temperatures will result in excessive temperatures within the rotor housing unless additional cooling methods are used. The rotor itself will have acceptable temperatures, but the pattern factor presented to the turbine may be of concern, depending on exhaust duct design and duct-rotor interaction.

Nalim, M. Razi↗
NASA NTRS2023

On the Onset Delays of Solar Energetic Electrons and Protons: Evidence for a Common Accelerator

The processes responsible for the acceleration of solar energetic particles (SEPs) are still not well understood, including whether SEP electrons and protons are accelerated by common or separate processes. Using a numerical particle transport model that includes both pitch-angle and perpendicular spatial diffusion, we simulate, among other quantities, the onset delay for MeV electrons and protons and compare the results to observations of SEPs from widely separated spacecraft. Such observations have previously been interpreted, in a simple scenario assuming no perpendicular diffusion, as evidence for different electron and proton sources. We show that, by assuming a common particle source together with perpendicular diffusion, we are able to simultaneously reproduce the onset delays for both electrons and protons. We argue that this points toward a common accelerator for these particles. Moreover, a relatively broad particle source is required in the model to correctly describe the observations. This is suggestive of diffusive shock acceleration occurring at large shock structures playing a significant role in the acceleration of these SEPs.

R D Strauss↗
DOE OSTI2018

Fluctuating hydrodynamics of reactive liquid mixtures

Fluctuating hydrodynamics (FHD) provides a framework for modeling microscopic fluctuations in a manner consistent with statistical mechanics and nonequilibrium thermodynamics. This paper presents an FHD formulation for isothermal reactive incompressible liquid mixtures with stochastic chemistry. Fluctuating multispecies mass diffusion is formulated using a Maxwell-Stefan description without assuming a dilute solution, and momentum dynamics is described by a stochastic Navier-Stokes equation for the fluid velocity. We consider a thermodynamically consistent generalization for the law of mass action for non-dilute mixtures and use it in the chemical master equation (CME) to model reactions as a Poisson process. The FHD approach provides remarkable computational efficiency over traditional reaction-diffusion master equation methods when the number of reactive molecules is large, while also retaining accuracy even when there are as few as ten reactive molecules per hydrodynamic cell. We present a numerical algorithm to solve the coupled FHD and CME equations and validate it on both equilibrium and nonequilibrium problems. We simulate a diffusively driven gravitational instability in the presence of an acid-base neutralization reaction, starting from a perfectly flat interface. We demonstrate that the coupling between velocity and concentration fluctuations dominates the initial growth of the instability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Mathematical model for the Bridgman-Stockbarger crystal growing system

In a major technical breakthrough, a computer model for Bridgman-Stockbarger crystal growth was developed. The model includes melt convection, solute effects, thermal conduction in the ampule, melt, and crystal, and the determination of the curved moving crystal-melt interface. The key to the numerical method is the use of a nonuniform computational mesh which moves with the interface, so that the interface is a mesh surface. In addition, implicit methods are used for advection and diffusion of heat, concentration, and vorticity, for interface movement, and for internal gracity waves. This allows large time-steps without loss of stability or accuracy. Numerical results are presented for the interface shape, temperature distribution, and concentration distribution, in steady-state crystl growth. Solutions are presented for two test cases using water, with two different salts in solution. The two diffusivities differ by a factor of ten, and the concentrations differ by a factor of twenty.

Roberts, G. O.↗
DOE OSTI2018

Compact Heat Exchanger Design and Testing for Advanced Reactors and Advanced Power Cycles

The goal of the proposed research is to demonstrate the thermal hydraulic performance of innovative surface geometries in compact heat exchangers used as intermediate heat exchangers (IHXs) and recuperators for the supercritical carbon dioxide (s-CO 2 ) Brayton cycle. Printed-circuit heat exchangers (PCHEs) are the primary compact heat exchangers of interest. The overall objectives are: To develop optimized PCHE designs for different working fluid combinations including helium to s-CO 2 , liquid salt to s-CO 2 , sodium to s-CO 2 , and liquid salt to helium; To experimentally and numerically investigate thermal performance, thermal stress and failure mechanism of PCHEs under various transients; and To study diffusion bonding techniques for elevated-temperature alloys and examine post-test material integrity of the PCHEs. The project objectives were accomplished by defining and executing five different tasks corresponding to these specific objectives. The first task involved a thorough literature review and a selection of IHX candidates with different surface geometries as well as a summary of prototypic operational conditions. The second task involved optimization of PCHE design with numerical analyses of thermal-hydraulic performances and mechanical integrity. The subsequent task dealt with the development of testing facilities and engineering design of PCHE to be tested in s-CO 2 fluid conditions. The next task involved experimental investigation and validation of the thermal-hydraulic performances and thermal stress distribution of prototype PCHEs manufactured with particular surface geometries. The last task involved an investigation of diffusion bonding process and posttest destructive testing to validate mechanical design methods adopted in the design process. The experimental work utilized the two test facilities at The Ohio State University (OSU) including one existing High-Temperature Helium Test Facility (HTHF) and the newly developed s-CO 2 test loop (STL) facility and s-CO 2 test facility at University of Wisconsin – Madison (UW).

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1997

The Implicit and Explicit alpha-mu Schemes

Artificial numerical dissipation is an important issue in large Reynolds number computations. In such computations, the artificial dissipation inherent in traditional numerical schemes can overwhelm the physical dissipation and yield inaccurate results on meshes of practical size. In the present work, the space-time conservation element and solution element method is used to construct new and accurate numerical schemes such that artificial numerical dissipation will not overwhelm physical dissipation. Specifically, these schemes have the property that numerical dissipation vanishes when the physical viscosity goes to zero. These new schemes therefore accurately model the physical dissipation even when it is extremely small. The method of space-time conservation element and solution element, currently under development, is a nontraditional numerical method for solving conservation laws. The method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. Explicit solvers for model and fluid dynamic conservation laws have previously been investigated. In this paper, we introduce a new concept in the design of implicit schemes, and use it to construct two highly accurate solvers for a convection-diffusion equation. The two schemes become identical in the pure convection case, and in the pure diffusion case. The implicit schemes are applicable over the whole Reynolds number range, from purely diffusive equations to purely inviscid (convective) equations. The stability and consistency of the schemes are analyzed, and some numerical results are presented. It is shown that, in the inviscid case, the new schemes become explicit and their amplification factors are identical to those of the Leapfrog scheme. On the other hand, in the pure diffusion case, their principal amplification factor becomes the amplification factor of the Crank-Nicolson scheme. We also construct an explicit solver with the treatment of diffusion being based on that in the implicit solvers. The explicit solver has only a CFL stability limitation on the Courant number, yet it retains the second-order spatial accuracy of the implicit schemes.

Chang, Sin-Chung↗
DOE OSTI2018

The Effect of Flowing Velocity and Impact Angle on the Fluid-Accelerated Corrosion of X65 Pipeline Steel in a Wet Gas Environment Containing CO{sub 2}

This study investigated the effect of flowing velocity and impact angle on the fluid-accelerated corrosion (FAC) of X65 pipeline steel in a wet CO{sub 2}-containing environment by corrosion morphology analysis, electrochemical measurement, and hydromechanical numerical analysis. It has been found that the corrosion behavior was controlled by the ion diffusion rate and shear stress due to the flowing gas. With an increase in the flowing velocity, the pitting corrosion initially decreased and then increased when the flowing velocity exceeded 14 m/s; this increase is caused by the augmentation of shear stress. When the impact angle decreased, the position of serious corrosion gradually moved outwards.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Two-dimensional Radiative Magnetohydrodynamic Simulations of Partial Ionization in the Chromosphere. II. Dynamics and Energetics of the Low Solar Atmosphere

We investigate the effects of interactions between ions and neutrals on the chromosphere and overlying corona using 2.5D radiative MHD simulations with the Bifrost code. We have extended the code capabilities implementing ion–neutral interaction effects using the generalized Ohm’s law, i.e., we include the Hall term and the ambipolar diffusion (Pedersen dissipation) in the induction equation. Our models span from the upper convection zone to the corona, with the photosphere, chromosphere, and transition region partially ionized. Our simulations reveal that the interactions between ionized particles and neutral particles have important consequences for the magnetothermodynamics of these modeled layers: (1) ambipolar diffusion increases the temperature in the chromosphere; (2) sporadically the horizontal magnetic field in the photosphere is diffused into the chromosphere, due to the large ambipolar diffusion; (3) ambipolar diffusion concentrates electrical currents, leading to more violent jets and reconnection processes, resulting in (3a) the formation of longer and faster spicules, (3b) heating of plasma during the spicule evolution, and (3c) decoupling of the plasma and magnetic field in spicules. Our results indicate that ambipolar diffusion is a critical ingredient for understanding the magnetothermodynamic properties in the chromosphere and transition region. The numerical simulations have been made publicly available, similar to previous Bifrost simulations. This will allow the community to study realistic numerical simulations with a wider range of magnetic field configurations and physics modules than previously possible.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2015

Numerical Study of the High-Speed Leg of a Wind Tunnel

The paper describes a numerical study of the high-speed leg of the NASA Langley 14 x 22-ft Low-Speed Wind Tunnel. The high-speed leg consists of the settling chamber, contraction, test section, and first diffuser. Results are shown comparing two different sources of surface geometry, and two different unstructured grid solvers for the flow characteristics. Numerical simulations of the flow on the tunnel centerline, boundary layer profiles on the floor, and wall static pressures have been compared with experiment. Flow angularities along the test section length have also been determined.

Nayani, Sudheer↗
NASA NTRS2001

Combustion of Metals in Carbon Dioxide and Reduced-Gravity Environments

Ongoing exploration and future mission2001110444 s to Mars have given impetus to research on the use of natural resources of the planet. Since carbon dioxide (CO2) constitutes approximately 95% of the Mars atmosphere and since it reacts directly and vigorously with several metals, this investigation focuses on metal-CO2 reactions as a possible combination for rocket-propellant production and energy generation. Magnesium (Mg) has been initially selected as the metal fuel owing to its low ignition temperature and high specific impulse and burning rate in CO2. Our studies in this field started with low gravity (g) combustion tests of Mg in O2, CO2, and CO. Reduced gravity provided a clear picture of the burning phenomena by eliminating the intrusive buoyant flows in high-temperature metal reactions and by removing the destructive effect of gravity on the shape of molten metal samples. Suspended cylindrical metal samples of 2, 3, and 4-mm in diameter and length were radiatively ignited in low-g to generate free-floating samples exhibiting a spherically symmetric flame with increasing metal-oxide accumulation in an outer shell. For the Mg-CO2 combination, burning times twice as long as in normal-g and five times longer than in Mg-O2 flames were observed, revealing a diffusion-controlled reaction. The burning time is proportional to the square of the sample diameter. In tests conducted with pure CO, combustion was not possible without constant heating of the sample due to the formation of a thick carbon-containing coating around the Mg sample generated by surface reactions. The following work presents two new studies that attempt to explain some of the low-g experimental observations. First, a simplified one-dimensional, quasi-steady numerical model is developed to obtain temperature, species concentrations, and burning rates of the spherically symmetric diffusion flame around the Mg sample burning in O2 and CO2. Second, a Planar Laser Induced Fluorescence (PLIF) technique is implemented to provide spatially resolved measurements of magnesium oxide (MgO) in the reaction zone of Mg samples burning in O2 and CO2. These experiments reveal fundamental differences between the two combustion systems.

Branch, M. C.↗
NASA NTRS1987

Development and evaluation of improved numerical schemes for recirculating flows

The paper examines the performance of the flux-spline scheme for convection-diffusion. Computations are presented for a number of test cases, both linear and nonlinear. It is shown that in all cases the flux-spline scheme yields results which are superior to those obtained with the lower-order formulations such as hybrid differencing. In order to improve the computational efficiency, the flux-spline scheme has been combined with a direct solution algorithm for the continuity and momentum equations. Such an approach eliminates the need for an equation for pressure or pressure correction and is found to be rapidly convergent.

Patankar, S. V.↗