Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

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 739 records · Page 41

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗

Structural Fluctuations and Thermophysical Properties of Molten II-VI Compounds

The objectives of the project are to conduct ground-based experimental and theoretical research on the structural fluctuations and thermophysical properties of molten II-VI compounds to enhance the basic understanding of the existing flight experiments in microgravity materials science programs as well as to study the fundamental heterophase fluctuation phenomena in these melts by: 1) conducting neutron scattering analysis and measuring quantitatively the relevant thermophysical properties of the II-VI melts (such as viscosity, electrical conductivity, thermal diffusivity and density) as well as the relaxation characteristics of these properties to advance the understanding of the structural properties and the relaxation phenomena in these melts and 2) performing theoretical analyses on the melt systems to interpret the experimental results. All the facilities required for the experimental measurements have been procured, installed and tested. It has long been recognized that liquid Te presents a unique case having properties between those of metals and semiconductors. The electrical conductivity for Te melt increases rapidly at melting point, indicating a semiconductor-metal transition. Te melts comprise two features, which are usually considered to be incompatible with each other: covalently bound atoms and metallic-like behavior. Why do Te liquids show metallic behavior? is one of the long-standing issues in liquid metal physics. Since thermophysical properties are very sensitive to the structural variations of a melt, we have conducted extensive thermophysical measurements on Te melt.

Su, Ching-Hua↗

Long-time asymptotics of a system for plasma diffusion

A system of parabolic nonlinear equations that describe the diffusion of a fully collisional plasma across a strong magnetic field is discussed. It is demonstrated that the solution to this system tends to a time asymptotic state which is of space-time separable form, theta(t)f(x). Furthermore, f(x) is independent of the initial conditions and theta(t) depends slightly on the initial conditions. The rate of decay of the temporal part is governed by a nonlinear eigenvalue problem. Since the equations are considered in a bounded domain we are able to analyze the effect of boundary conditions on the evolution of the system. Additional effects on radiation, heating, and particle injection can also be accounted for. Essential differences between the behavior of a fully-coupled system and a scalar equation are observed.

Rosenau, P.↗

Time dispersion of energetic solar particles - Unexpected velocity and species dependence

The intensity-time behavior for protons and alpha-particles, as well as for carbon, oxygen, and iron ions, has been measured following the solar flare of September 19, 1974, for energies between 0.5 and about 5 MeV per nucleon. The profiles display a time dispersion which is inversely proportional to velocity for each individual species. In addition, at a given velocity the time dispersion also depends on the charge-to-mass ratio of the ion. Based on this latter dependence, it is concluded that while carbon and oxygen are essentially fully stripped, iron nuclei are not, having an effective charge of 10 plus or minus 5. The observed dispersion cannot be explained by purely rigidity-dependent diffusive propagation.

Ogallagher, J. J.↗

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗

Radiation Damage Workshop

The availability of data regarding the radiation behavior of GaAs and silicon solar cells is discussed as well as efforts to provide sufficient information. Other materials are considered too immature for reasonable radiation evaluation. The lack of concern over the possible catastrophic radiation degradation in cascade cells is a potentially serious problem. Lithium counterdoping shows potential for removing damage in irradiated P-type material, although initial efficiencies are not comparable to current state of the art. The possibility of refining the lithium doping method to maintain high initial efficiencies and combining it with radiation tolerant structures such as thin BSF cells or vertical junction cells could provide a substantial improvement in EOL efficiencies. Laser annealing of junctions, either those formed ion implantation or diffusion, may not only improve initial cell performance but might also reduce the radiation degradation rate.

Stella, P. M.↗

Solar flares and avalanches in driven dissipative systems

The contention of Lu and Hamilton (1991) that the energy release process in solar flares can be understood as avalanches of many small reconnection events is further developed. The dynamics of the complex magnetized plasma of solar active regions is modeled with a simple driven dissipative system, consisting of a vector field with local instabilities that cause rapid diffusion of the field. It is argued that the avalanches in this model are analogous to solar flares. The distributions of avalanches in this model are compared with the solar flare frequency distributions obtained from ISEE 3/ICE satellite observations. Quantitative agreement is found with the energy, peak luminosity, and duration distributions over four orders of magnitude in flare energy, from the largest flares down to the completeness limit of the observations. It is predicted that the power-law solar flare frequency distributions will be found to continue downward with the same logarithmic slopes to an energy of about 3 x 10 exp 25 ergs and duration of about 0.3 s, with deviations from power-law behavior below these values.

Lu, Edward T.↗

Numerical prediction of wakes in cascades and compressor rotors including the effects of mixing. II - Rotor passage flow and wakes including the effects of spanwise mixing

The results of a numerical investigation to predict the flow field including wakes and mixing in axial-flow compressor rotors are presented. The wake behavior in a moderately loaded compressor rotor is studied numerically using a 3D incompressible Navier-Stokes solver with a high Reynolds number form of a turbulence model. The equations are solved using a time dependent implicit technique. The agreement between the measured data and the predictions is good; including the blade boundary-layer profiles, wake mean-velocity profiles, and decay. The ability of the pseudocompressibility scheme to predict the entire flow field including the near and far wake profiles and its decay characteristics, effect of loading, and the viscous losses of a 3D rotor flow field are demonstrated. The mixing in the downstream regions away from the hub and annulus walls is dominated by wake diffusion. In regions away from the walls the radial mixing is predominantly caused by the transport of mass, momentum, and energy by the radial component of velocity in the wake.

Lakshminarayana, B.↗

Extensional flow convecting a reactant undergoing a first order homogeneous reaction and diffusional mass transfer from a sphere at low to intermediate Peclet and Damkohler numbers

Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Shah, N. Y.↗

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

ionic liquids↗

Investigation of Structure and Transport in Li-Doped Ionic Liquid Electrolytes: [pyr14][TFSI], [pyr13][FSI], [EMIM][BF4]

Ionic liquid electrolytes have been proposed as a means of improving the safety and cycling behavior of advanced lithium batteries; however, the properties of these electrolytes under high lithium doping are poorly understood. Here, we employ both polarizable molecular dynamics simulation and experiment to investigate the structure, thermodynamics and transport of three potential electrolytes, N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N- methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-- methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li (-) salt concentration and temperature. Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi we find the presence of lithium aggregates. Furthermore, the computed density, diffusion, viscosity, and ionic conductivity show excellent agreement with experimental data. While the diffusion and viscosity exhibit a systematic decrease and increase, respectively, with increasing xLi, the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1-0.3 mScm. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions, which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, we comment on the relative kinetics of Li(+) transport in each liquid and we present strong evidence for transport through anion exchange (hopping) as opposed to the net motion of Li(+) with its solvation shell (vehicular).

molecular dynamics↗

Reliability Issues with PME and BME Ceramic Capacitors

Exposure of chip MnO2 tantalum capacitors to humid environments might result in increased ESR, leakage currents, and first turn-on failures. However, there is a lack of literature data on the effect of moisture on reverse bias behavior of the parts. The presence of moisture can also result in pop-corning when a high water vapor pressure develops when moisture absorbed in pores of tantalum slugs vaporizes instantly during soldering process resulting in damage to capacitors. A study of kinetics of moisture ingress to and release from active elements of capacitors would allow a better understanding of degradation mechanisms and is important for preventing failures. In this work, a technique for investigation of moisture sorption and desorption in solid chip tantalum capacitors that employs tantalum slugs as a humidity sensor have been developed and kinetics of the process analyzed for different types of capacitors at temperatures from room to 125 C. A model that relates diffusion characteristics of polymer cases and size of the slugs to characteristic times of moisture sorption has been developed. A strong effect of moisture on long-term degradation of reverse bias currents in MnO2 cathode capacitors has been demonstrated and physical mechanisms discussed.

degradation↗

Kinetics of Moisture Sorption and Reverse Bias Degradation in Chip Tantalum Capacitors

Exposure of chip MnO2 tantalum capacitors to humid environments might result in increased ESR, leakage currents, and first turn-on failures. However, there is a lack of literature data on the effect of moisture on reverse bias behavior of the parts. The presence of moisture can also result in pop-corning when a high water vapor pressure develops when moisture absorbed in pores of tantalum slugs vaporizes instantly during soldering process resulting in damage to capacitors. A study of kinetics of moisture ingress to and release from active elements of capacitors would allow a better understanding of degradation mechanisms and is important for preventing failures. In this work, a technique for investigation of moisture sorption and desorption in solid chip tantalum capacitors that employs tantalum slugs as a humidity sensor have been developed and kinetics of the process analyzed for different types of capacitors at temperatures from room to 125 ºC. A model that relates diffusion characteristics of polymer cases and size of the slugs to characteristic times of moisture sorption has been developed. A strong effect of moisture on long-term degradation of reverse bias currents in MnO2 cathode capacitors has been demonstrated and physical mechanisms discussed.

Teverovsky, Alexander↗

Kinetics of Moisture Sorption and Reverse Bias Degradation in Chip Tantalum Capacitors

Exposure of chip MnO2 tantalum capacitors to humid environments might result in increased ESR, leakage currents, and first turn-on failures. However, there is a lack of literature data on the effect of moisture on reverse bias behavior of the parts. The presence of moisture can also result in pop-corning when a high water vapor pressure develops when moisture absorbed in pores of tantalum slugs vaporizes instantly during soldering process resulting in damage to capacitors. A study of kinetics of moisture ingress to and release from active elements of capacitors would allow a better understanding of degradation mechanisms and is important for preventing failures. In this work, a technique for investigation of moisture sorption and desorption in solid chip tantalum capacitors that employs tantalum slugs as a humidity sensor have been developed and kinetics of the process analyzed for different types of capacitors at temperatures from room to 125 C. A model that relates diffusion characteristics of polymer cases and size of the slugs to characteristic times of moisture sorption has been developed. A strong effect of moisture on long-term degradation of reverse bias currents in MnO2 cathode capacitors has been demonstrated and physical mechanisms discussed.

Teverovsky, Alexander↗

Flight Experiment to Study Double-Diffusive Instabilities in Silver-Doped Lead Bromide Crystals

A detailed study on the effect of convection on crystal quality was carried out by growing lead bromide crystals in transparent Bridgman furnace. Direct observations were made on the solid-liquid interface and a new kind of instability was observed. This could be explained on the basis of toroidal flow in the AgBr-doped lead bromide sample. With the increasing translation velocity, the interface changed from flat to depressed, and then formed a cavity in the center of the growth tube. The crystal grown at the lowest thermal Rayleigh number showed the highest quality and crystal grown at the largest thermal Rayleigh number showed the worst quality. Numerical studies were carried out to provide a framework for interpreting the observed convective and morphological instabilities, and to determine the critical (limiting) concentration of dopant for a particular growth velocity and gravity level. Theoretical instability diagrams were compared with data obtained from the experimental studies. These studies provided basic data on convective behavior in doped lead bromide crystals grown by the commercially important Bridgman process.

Singh, N. B.↗

Calculations of turbulent mass transport in a bluff-body diffusion-flame combustor

Experimental measurements of turbulent mass transit in a bluff-body diffusion-flame combustor (using CO2 instead of fuel) are analyzed to evaluate the accuracy of physical models used in flow computations. The data of Lightman and Magill (1981), Lightman et al. (1983), and Roquemore et al. (1983) are used to calculate apparent turbulent Schmidt numbers (TSN) for a series of flow conditions by a modified TEACH viscous-flow code. The modeling principles, calculation grid, and boundary conditions are discussed, and the results are presented in graphs comparing calculated and measured values. It is shown that models such as the two-equation (K - epsilon) model which use a single value of the TSN are inappropriate for conditions typical of gas-turbine combustors: local TSN variations due to turbulence behavior must be taken into account in improved models.

Sturgess, G. J.↗

Creep Resistance of ZrO2 Ceramic Improved by the Addition of a Small Amount of Er2O3

Zirconia (ZrO2) has great technological importance in structural, electrical, and chemical applications. It is the crucial component for state-of-the art thermal barrier coatings and an enabling component as a solid electrolyte for solid-oxide fuel cell systems. Pure ZrO2 is of limited use for industrial applications because of the phase transformations that occur. Upon the addition of stabilizers, cubic (c-ZrO2) and tetragonal (t-ZrO2) forms can be preserved. It is the stabilized and partially stabilized forms of zirconia that function as thermal barrier coatings, solid electrolytes, and oxygen sensors and that have numerous applications in the electrochemical industry. The cubic form of ZrO2 is typically stabilized through Y2O3 additions. However, Y2O3-stabilized zirconia is susceptible to deformation at high temperatures (greater than 900 C) because of the large number of slip systems and the high oxygen diffusion rates, which result in high creep rates at high temperatures. Successful use of ZrO2 at high temperatures requires that new dopant additives be found that will retain or enhance the desirable properties of cubic ZrO2 and yet produce a material with lower creep rates. At the NASA Glenn Research Center, erbium oxide (Er2O3) was identified as a promising dopant for improving the creep resistance of. ZrO2. The selection of Er2O3 was based on the strong interactions of point defects and dislocations. Single crystals of 5 mol% Er2O3- doped ZrO2 rods (4 mm in diameter) and monofilaments (200 to 300 mm in diameter and 30 cm long) were grown using the laser-heated float zone technique, and their creep behavior was measured as a function of temperature. The addition of 5 mol% Er2O3 to single-crystal ZrO2 improved its creep resistance at high temperatures by 2 to 3 orders of magnitude over state-of-the-art Y2O3-doped crystals. Detailed microstructural characterization of ZrO2-Er2O3 single crystals has identified new mechanisms for improving the creep resistance of this class of materials. Adding Er2O3 to ZrO2 results in microstructure of stable and metastable tetragonal precipitates that with thermal treatment evolve to a tweed structure of nanosize tetragonal lamellae. The superior high-temperature creep resistance of Er2O3-doped ZrO2 is attributed to nanoscale precipitation hardening. Doping with Er2O3 will significantly increase the upper-use temperature limit of ZrO2. Potential applications include using Er2O3-doped ZrO2 as a high-temperature fiber for structural applications and adding Er2O3 to reduce the sintering rates of ZrO2 thermal barrier coatings. This work was conducted at Dpto. de F sica de la Materia Condensada, Universidad de Sevilla, Spain, and at NASA Glenn.

Martinez-Fernandez, Julian↗