Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Diffusion rate”

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 847 records · Page 47

Observations on the nature and mechanism of the bounded turbulent shear flow structure

A physical model of the bounded turbulent shear flow structure based on the viscous diffusion of vortex lines at the wall and their stretching by the strain rate field is presented. Evidence from the literature supporting the model is briefly reviewed, and new flow visualization evidence is presented. Progress toward quantifying the structure's characteristics with a multi-sensor hot-wire vorticity probe is described.

Wallace, J. M.↗

The evolution and growth of Al2O3 scales on beta-NiAl

The formation and growth of Al2O3 scales on (beta)-NiAl were studied using electron microscopy and other analytical techniques to gain an understanding of the oxidation properties of (beta)-NiAl and of alumina-forming alloys, in general. The transient and mature stages of oxidation were studied as well as the transformation stage during which the oxide scale transforms from metastable Al2O3 phases to the thermodynamically stable alpha-Al2O3 phase. The transient oxidation stages were studied at 800 deg C and for short times at 1100C. At 800C, the scales consist predominantly of delta-Al2O3 which forms by cation vacancy ordering in the defective spinel lattice of gamma-Al2O3. At 1100C, a fast-growth morphology of theta-Al2O3 forms as a surface layer over delta-Al2O3. For both oxidation temperatures, the scales are often epitaxially oriented with respect to the metal. The transient scales grow by outward cation diffusion as evidenced by surface growth morphologies. The transformation to alpha-Al2O3 occurs within 1 hour at 1100C by a nucleation and radial growth process. The large volume decrease associated with the transformation results in a highly strained alpha-Al2O3 microstructure. A change in scale growth mechanism from outward cation to inward anion diffusion allows transient surface morphologies to be smoothed by surface diffusion. The mature stage of oxidation involves the growth of an alpha-Al2O3 scale having the lacey morphology formed as a result of the gamma yields alpha transformation. Growth of the scale occurs by counterdiffusion along grain boundaries resulting in ridges formed by impingement of alpha-Al2O3 nuclei during the transformation stage. Also scale growth occurs by inward oxygen diffusion through healed cracks; the cracks result from transformation stresses. The measured growth rates of scales having the lacey morphology are an order of magnitude less than fine-grained alpha-Al2O3 scales. Metal orientations were found to have a large effect on oxide morpholoies during all stages of oxidation.

Doychak, J. K.↗

Grain boundary oxidation and its effects on high temperature fatigue life

Fatigue lives at elevated temperatures are often shortened by creep and/or oxidation. Creep causes grain boundary void nucleation and grain boundary cavitation. Grain boundary voids and cavities will accelerate fatigue crack nucleation and propagation, and thereby shorten fatigue life. The functional relationships between the damage rate of fatigue crack nucleation and propagation and the kinetic process of oxygen diffusion depend on the detailed physical processes. The kinetics of grain boundary oxidation penetration was investigated. The statistical distribution of grain boundary penetration depth was analyzed. Its effect on high temperature fatigue life are discussed. A model of intermittent micro-ruptures of grain boundary oxide was proposed for high temperature fatigue crack growth. The details of these studies are reported.

Liu, H. W.↗

Latitudinal variation of perturbation electric fields during magnetically disturbed periods - 1986 Sundial observations and model results

F-region incoherent scatter radar drift observations from Millstone Hill and Jicamarca, h-prime F observations from Huancayo, and high latitude ground-magnetometer measurements taken during the Sundial 1986 campaign are used to study the relationship between plasmaspheric electric field perturbations and high latitude currents during disturbed periods. The observations are in good agreement with numerical results from a Rice Covection Model run that involved a sharp increase in the polar cap potential drop followed by a subsequent decrease. The zonal disturbance electric field pattern is latitude independent, and the corresponding amplitudes change approximately as L exp n (where n is about 1.5). The meridional electric field patterns and amplitudes have larger latitudinal variations. The mid-, low, and equatorial electric fields from the Rice Convection Model are in good agreement with previous results from the semianalytic, Senior-Blanc (1987) model. Also discussed are three physical mechanisms (over-shielding, fossil winds, and magnetic reconfiguration) that contribute to the long lasting (1-2 h) equatorial zonal electric field perturbations associated with a sudden northward turning of the IMF. It is predicted that the penetration of high latitude electric fields to low latitudes should, in general, be closely related to the rate of motion of the shielding layer and the equatorward edge of the diffuse aurora.

Fejer, B. G.↗

Dispersal of gases generated near a lunar outpost

The dispersal of gases evolved by prospective lunar base operations is presently modeled analytically using continuous (mining and habitat-venting) and impulsive (primarily rocket exhaust) injections of gases. In the case of impulsive injection, the neutral atmosphere and associated ionosphere both decay on time-scales of about 20 min; in that of continuous injection, the atmosphere near the base grows and achieves a steady state after about 20 min. Both direct and diffusive transport mechanisms are considered, and it is concluded that for the injection rates and assumptions presently employed the artificial lunar atmospheres produced are not significantly detrimental to astronomical observations and high-vacuum materials processing operations.

Burns, Jack O.↗

Modelling of the dynamics of HgCdTe growth by the vertical Bridgman method

The analysis of Kim and Brown (1989), which described the pseudo-steady state of the HgCdTe melt/crystal system grown by the vertical Bridgman method, is extended to the study of the transients during the crystal growth and, in particular, at the beginning of the growth. The time-dependent numerical simulations developed show the development of a diffusion layer in the melt and the direct relationship between this layer and the intensity and the pattern of convection in the melt. It is shown that the spread of the solute layer occurs on a diffusive time scale, with the initial transient to the steady-state growth rate lasting several hours for the slow growth rates that are typical of the HgCdTe growth. The simulations demonstrate all the features observed experimentally by Szofran et al. (1984).

Kim, Do H.↗

Effects of buoyancy-driven flow and thermal boundary conditions on physical vapor transport

A 2D numerical model was developed in order to ascertain if reduced gravity conditions are beneficial to physical vapor transport (PVT) and to determine its tolerance limits to residual accelerations. This was solved using the PHOENICS finite-volume code. Reduction of gravitational accelerations to less than 0.1 g0 was found to be sufficient to suppress buoyancy-driven convection to an extent that diffusion was the dominant transport mode, whence a greater uniformity in the growth rate could be obtained. It is shown that a uniform temperature gradient on the ampoule walls causes the vapor to be supersaturated throughout the ampoule, potentially resulting in undesirable nucleation at the walls. A 'hump' in the wall temperature profile can be used to avoid this. The prevailing transport conditions determine the size of the hump needed.

Nadarajah, Arunan↗

The application of preconditioning in viscous flows

The present time-derivative preconditioning algorithm is effective in flow conditions ranging from inviscid to very diffusive flows, as well as low subsonic to supersonic flow velocities. By means of a preconditioning matrix that (1) introduces well-conditioned eigenvalues and (2) avoids nonphysical time reversals for viscous flows, a mechanism is obtained which controls the inviscid and viscous time-step parameters at very diffusive flows. These capabilities are demonstrated for a variety of sample problems; convergence rates of solutions that are indistinguishable from those obtained without preconditioning are shown to be accelerated by as much as two orders of magnitude.

Choi, Y.-H.↗

Modeling of InP metalorganic chemical vapor deposition

The growth of InP by metalorganic chemical vapor deposition (MOCVD) in a horizontal reactor is being modeled with a commercially available computational fluid dynamics modeling code. The mathematical treatment of the MOCVD process has four primary areas of concern: 1) transport phenomena, 2) chemistry, 3) boundary conditions, and 4) numerical solution methods. The transport processes involved in CVD are described by conservation of total mass, momentum, energy, and atomic species. Momentum conservation is described by a generalized form of the Navier-Stokes equation for a Newtonian fluid and laminar flow. The effect of Soret diffusion on the transport of particular chemical species and on the predicted deposition rate is examined. Both gas-phase and surface chemical reactions are employed in the model. Boundary conditions are specified at the inlet and walls of the reactor for temperature, fluid flow and chemical species. The coupled set of equations described above is solved by a finite difference method over a nonuniform rectilinear grid in both two and three dimensions. The results of the 2-D computational model is presented for gravity levels of zero- and one-g. The predicted growth rates at one-g are compared to measured growth rates on fused silica substrates.

Black, Linda R.↗

Testing the Li-Strahler four-component canopy reflectance model in the HAPEX-Sahel shrub savanna sites using ground reflectance data

The rate at which a light field decays in water is characterized by the diffuse attenuation coefficient k. The Li-Strahler discrete-object canopy reflectance model was tested in two sites, a shrub grass savanna and a degraded shrub savanna on bare soil, in the proposed HAPEX (Hydrologic Atmospheric Pilot Experiment) II/Sahel study area in Niger, West Africa. Average site reflectance was predicted for each site from the reflectances and cover proportions of four components: shrub canopy, background (soil or grass and soil), shaded canopy, and shaded background. Component reflectances were sampled in the SPOT wavebands using a hand-held radiometer. Predicted reflectance was compared to average site reflectance measured using the same radiometer mounted on a backpack with measurements recorded every 5 m along two 1-km transects, also in the SPOT (Systeme Probatoire d'Observation de la Terre) bands. Measurements and predictions were made for each of the three days during the summer growing season, approximately two weeks apart. Red, near infrared reflectance, and the NDVI (normalized difference vegetation index) were all predicted with a high degree of accuracy for the shrub/grass site and with reasonable accuracy for the degraded shrub site.

Franklin, J.↗

Acoustic instability driven by cosmic-ray streaming

We study the linear stability of compressional waves in a medium through which cosmic rays stream at the Alfven speed due to strong coupling with Alfven waves. Acoustic waves can be driven unstable by the cosmic-ray drift, provided that the streaming speed is sufficiently large compared to the thermal sound speed. Two effects can cause instability: (1) the heating of the thermal gas due to the damping of Alfven waves driven unstable by cosmic-ray streaming; and (2) phase shifts in the cosmic-ray pressure perturbation caused by the combination of cosmic-ray streaming and diffusion. The instability does not depend on the magnitude of the background cosmic-ray pressure gradient, and occurs whether or not cosmic-ray diffusion is important relative to streaming. When the cosmic-ray pressure is small compared to the gas pressure, or cosmic-ray diffusion is strong, the instability manifests itself as a weak overstability of slow magnetosonic waves. Larger cosmic-ray pressure gives rise to new hybrid modes, which can be strongly unstable in the limits of both weak and strong cosmic-ray diffusion and in the presence of thermal conduction. Parts of our analysis parallel earlier work by McKenzie & Webb (which were brought to our attention after this paper was accepted for publication), but our treatment of diffusive effects, thermal conduction, and nonlinearities represent significant extensions. Although the linear growth rate of instability is independent of the background cosmic-ray pressure gradient, the onset of nonlinear eff ects does depend on absolute value of DEL (vector differential operator) P(sub c). At the onset of nonlinearity the fractional amplitude of cosmic-ray pressure perturbations is delta P(sub C)/P(sub C) approximately (kL) (exp -1) much less than 1, where k is the wavenumber and L is the pressure scale height of the unperturbed cosmic rays. We speculate that the instability may lead to a mode of cosmic-ray transport in which plateaus of uniform cosmic-ray pressure are separated by either laminar or turbulent jumps in which the thermal gas is subject to intense heating.

Begelman, Mitchell C.↗

Photochemical modeling of Titan's atmosphere

We have developed a new photochemical model of Titan's atmosphere which includes all the important compounds and reactions in spherical geometry from the surface to 1240 km. Compared to the previous model of Yung et al. (1984), the most significant recent change in the reactions used is the updated methane dissociation scheme (Mordaunt et al. 1993). Moreover, the transfer of the solar radiation in the atmosphere and the photolysis rates have been calculated by using a Monte Carlo code. Finally, the eddy diffusion coefficient profile is adjusted in order to fit the mean vertical distribution of HCN retrieved from millimeter groundbased observations of Tanguy et al. (1990); using new values for the boundary flux of atomic nitrogen (Strobel et al. 1992). We have run the model in both steady-state and diurnal modes, with 62 species involved in 249 reactions. There is little difference between diurnal and steady-state results. Overall our results are in a closer agreement with the abundances inferred from the Voyager infrared measurements at the equator than the Yung et al. results. We find that the catalytic scheme for H recombination invoked by Yung et al. only slightly improves the model results and we conclude that this scheme is not essential to fit observations.

Toublanc, D.↗

Annual Research Briefs, 1998

The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Spinks, Debra↗

Hydrocarbon photochemistry in the upper atmosphere of Jupiter

The hydrocarbon photochemistry in the upper atmosphere of Jupiter is investigated using a one-dimensional, photochemical-diffusive, and diurnally averaged model. The important chemical cycles and pathways among the major species are outlined and a standard model for the North Equatorial Belt region is examined in detail. It is found that several traditionally dominant chemical pathways among the C and C2 species are replaced in importance by cycles involving C-C4 species. The pressure and altitude profiles of mixing ratios for several observable hydrocarbon species are compared with available ultraviolet- and infrared-derived abundances. The results of sensitivity studies on the standard model with respect to variations in eddy diffusion profile, solar flux, atomic hydrogen influx, latitude, temperature, and important chemical reaction rates are presented. Measured and calculated airglow emissions of He at 584 angstroms and H at 1216 angstroms are also used to provide some constraints on the range of model parameters. The relevance of the model results to the upcoming Galileo mission is briefly discussed. The model is subject to considerable improvement; there is a great need for laboratory measurements of basic reaction rates and photodissociation quantum yields, even for such simple species as methylacetylene and allene. Until such laboratory measurements exist there will be considerable uncertainty in the understanding of the C3 and higher hydrocarbons in the atmospheres of the jovian planets.

NASA Discipline Number 52-20↗

Predicting Thermal Response in a Li-ion Cell on a UAV Fight Profile

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Also, a thermal runaway event can propagate from a single "failed" cell to the pack in a battery pack, leading to a more significant event. Mitigating a thermal runaway event is important in the commercial and automotive sectors. However, preventing such events in an electric aircraft (or air taxis) is paramount due to the lack of alternatives in the event of a failure. Battery prognostics algorithms allow the prediction of state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle) [1]. For this presentation, we will extend this two-level battery predictive algorithm to predict SOC, EOL, and estimated maximum temperature during a simulated flight. The model is extended by integrating a lumped physics-driven thermal model for high current densities [2]. The parameters used to control SOC and EOL are maximum storable charge, time constant for Li-ion diffusivity in the carbon particles, and internal cell resistance. Cycling leads to an increase in the heat generated by an aged Li-ion cell with a LiyCoO2 (LCO) cathode and a LixC6 (MCMB) anode. The aging of a cell leads to increase in SEI layer thickness, the diffusion time for the lithium ions inside the electrodes, and the local reaction rates, in addition to the thermodynamic abuse caused by fixed cycling voltages controlled by a Battery Management System. As the battery ages, the cell resistance increases, while the onset temperature of the thermal runaway decreases (depends on the cell chemistry and cell abuse history). Any large deviation of the cell temperature from the estimated (expected) value can identify a faulty cell. Since SEI decomposition has the lowest onset temperature in the series of reactions leading to thermal runaway, the model considers the self-heating rate of the SEI decomposition as onset temperature (similar to Ref. [3]). The parameters in the Arrhenius equation for the SEI heating rate depend on the number of cycles, the cell's operating temperature, and the cell's abuse history [4,5]. Coupling the electrochemical, thermal, and aging model allow the prognostic algorithm to estimate a typical cell voltage and temperature as a function of age (cycling and calendar), whose departure from measured values from the BMS is used to identify a safety event. In addition, we will present the results from two simulated flight scenarios for a UAV: typical and extreme, since the power requirements vary significantly during take-off, landing, and changing altitudes, while the power requirements remain low during the cruise. For this presentation, the power requirement for a battery pack in a UAV is scaled to a single cell. This cell is cycled through a simulated profile, and the data is collected and used to predict a safety event.

Li-ion↗

SUPIN: A Tool for the Aerodynamic Design and Analysis of Supersonic Inlets

The Supersonic Inlet Design and Analysis Tool (SUPIN) is a computational program that performs geometric modeling and aerodynamic design and analysis of a set of supersonic inlets whose types include the axisymmetric pitot, three-dimensional pitot, axisymmetric, two-dimensional, two-dimensional bifurcated-duct, and streamline-traced inlets. The inlets are modeled by dividing the inlet flowpath into components that start at the freestream and ends at the engine face. The primary components of the inlet include the external supersonic diffuser, throat section, and subsonic diffuser. Each component is characterized by a set of inputs that include geometric and aerodynamic factors. Examples of the geometric factors are angles, lengths, area ratios, and coordinates that describe the geometry of the components. Examples of the aerodynamic factors are Mach numbers, flow ratios, and pressure ratios that set constraints or desired aerodynamic properties of the component. These input factors are specified within a text-based input data file that is read into SUPIN. SUPIN is a Fortran 95 program that performs the inlet design and analysis as a batch or non-interactive process. SUPIN writes output data files describing the inlet geometry and aerodynamic performance. SUPIN uses analytic, empirical, and numerical methods to perform the inlet design and analysis. The inlet geometry is constructed using planar and three-dimensional geometry constructs and is represented by coordinates, angles, areas, profiles, and surfaces. Surfaces are output in the form of a Plot3D surface grid file. SUPIN does have some capability to automatically generate planar and three-dimensional, multi-block, structured grids for computational fluid dynamics (CFD) analysis. The aerodynamic analysis performed by SUPIN computes quasi-one-dimensional properties, such as Mach numbers, pressures, and temperatures, through the inlet flowpath. Planar aerodynamic properties, such as those for shock and Mach waves, are computed for the supersonic compression of the external supersonic diffuser and throat section. The inlet performance is characterized by the inlet flow rates, total pressure recovery, and inlet drag. The primary intent of this document is to describe the usage of SUPIN. The methods used within SUPIN are documented and may serve as a useful reference for understanding the aerodynamic design and analysis of supersonic inlets.

supersonic inlets↗

The origin of the Gum nebula

Obsrvations and theoretical investigations of the Gum nebula (GN) since about 1971 are reviewed. Direct observations of the GN, the Vela X supernova remnant (SNR), the Vela pulsar, and other stars in or near the GN are discussed with those of related phenomena such as the radio loops and known SNRs; the emphasis is on studies of the interstellar absorption lines, the evidence for hot gas in the GN, and the extended diffuse emission. The four basic models proposed for the GN are considered: a fossil Stromgren sphere, an old SNR, an H II region, or a superbubble. The GN physical parameters predicted by each model are listed in a table and compared. A minimum explanation which attributes the 36 x 36-deg filamentary structure and the 125-pc radius structure to the action of the stellar winds from Zeta Pup and Gamma-2 Vel (and perhaps the effect of a Vel X supernova explosion 20,000 years ago) is found most appropriate, at least until the questions of the net expansion rate of the GN (about 20 km/sec or about zero?) and the existence of the diffuse emission beyond the filamentary structure are resolved by observations.

Bruhweiler, F. C.↗

Computational and Experimental Study of Energetic Material in a Counterflow Microgravity Environment

Ground based (normal gravity) combustion studies can provide important information on the processes by which monopropellants and composite systems burn. The effects of gravitational forces, however, can often complicate the interpretation of the models and the implementation of experiments designed to help elucidate complex issues. We propose to utilize a combined computational/experimental approach in a microgravity environment to understand the interaction of oxidizer-binder diffusion flames in composite propellants. By operating under microgravity conditions we will be able to increase the length scales and suppress the gravitational forces on melting binders such that increased resolution of both major and minor species will be possible thus reducing the demands placed on both the computational and diagnostic tools. Results of a detailed transport/finite rate chemistry model will be compared with nonintrusive optical diagnostic measurements of the structure and extinction of diffusion flames in which oxidizers such as ammonium perchlorate (AP) and ammonium dinitramide (ADN) are counterflowed against realistic binders such as hydroxyl-terminated polybutadiene (HTPB) and 3,3-bis(azidomethyl)oxetane (BAMO). The work proposed herein represents a collaborative effort among the research groups at Yale University, Princeton University and the Combustion Diagnostics Laboratory at the Naval Air Warfare Center in China Lake, CA.

Smooke, Mitchell D.↗