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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

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↗
NASA NTRS1996

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↗
NASA NTRSDate not supplied

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↗
NASA NTRS2024

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↗
DOE OSTI2023

Pressure effects on reactivity and extinction of n -dodecane diffusion cool flame

Low-temperature combustion at high pressure has gained increasing interest owing to the growing demand for clean and efficient propulsion and energy systems. In particular, cool flame structure, extinction limits, and fuel reactivities at high pressure are critical in affecting the engine performance at near-limit conditions. Here, this paper investigates the effects of pressure on the cool flame extinction limit, structure, radical index, reactivity, and oxygen concentration dependence by experiments, analysis, and modeling. A large n-alkane, namely n-dodecane (n C 12 H 26 ), is selected to study its diffusion cool flame dynamics and reactivity in a high-pressure counterflow burner up to 10 atm. The experimental results show that higher pressure increases the cool flame extinction strain rates, and that the pressure-weighted extinction strain rate (aP) is proportional to the square of pressure, P 2 . A scaling analysis explains the relationship between the dependence of flame structure, heat release rate, and pressure-weighted cool flame strain rate on pressure. Furthermore, radical indexes at different pressures are measured by isolating the thermal and transport effects from the chemical contribution to diffusion cool flame extinction. The radical index clearly shows that the low-temperature reactivity increases with pressure. In addition, due to the critical role of multiple oxygen addition reactions in low-temperature chemistry, the relationship between the cool flame extinction limit and the oxygen concentration is explored. It is found that the cool flame extinction limits are proportional to the n th power of the oxygen concentration, [O 2 ] n , and increasing pressure leads to stronger extinction limit dependence (larger n) on the oxygen concentration. The present experiment and detailed kinetic analysis show clearly that increasing pressure promotes the low-temperature chemistry including the oxygen addition reactions, while the scaling analysis explains well the experimental results.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1983

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

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

Creep and stress rupture of a mechanically alloyed oxide dispersion and precipitation strengthened nickel-base superalloy

The creep and stress rupture behavior of a mechanically alloyed oxide dispersion strengthened (ODS) and gamma-prime precipitation strengthened nickel-base alloy (alloy MA 6000E) was studied at intermediate and elevated temperatures. At 760 C, MA 6000E exhibits the high creep strength characteristic of nickel-base superalloys and at 1093 C the creep strength is superior to other ODS nickel-base alloys. The stress dependence of the creep rate is very sharp at both test temperatures and the apparent creep activation energy measured around 760 C is high, much larger in magnitude than the self-diffusion energy. Stress rupture in this large grain size material is transgranular and crystallographic cracking is observed. The rupture ductility is dependent on creep strain rate, but usually is low. These and accompanying microstructural results are discussed with respect to other ODS alloys and superalloys and the creep behavior is rationalized by invoking a recently-developed resisting stress model of creep in materials strengthened by second phase particles.

Howson, T. E.↗
NASA NTRS1992

Co diffusion and growth of buried single-crystal CoSi2 in Si(111) by endotaxy

At a growth temperature of 800 C, Co deposited on a Si capping layer exhibits oriented growth on buried CoSi2 grains on Si(111), a process referred to as endotaxy. This occurs preferentially to surface nucleation of CoSi2, provided the thickness of the Si cap is less than a critical value of about 100 nm for deposition rates of 0.003-0.01 nm/s. The steady-state process is modeled using known values of the Co diffusion coefficient and solid solubility in Si, allowing some conclusions to be drawn regarding parameters relevant to CoSi2 epitaxy. Using this technique, single-crystal continuous layers of CoSi2 can be formed under a high-quality Si capping layer.

Fathauer, R. W.↗
NASA NTRS1986

Thermofluid analysis of the SSME preburner using a gas-gas diffusion model for oxygen and hydrogen combustion at supercritical pressures

The paper discusses the thermofluid analysis of the Space Shuttle Main Engine (SSME) fuelside preburner. The governing equations have been solved numerically to predict flow, heat transfer, mixing, and combustion. A two-fluid approach is adopted in which oxygen is regarded as one fluid and hydrogen is regarded as the other fluid. The chemical kinetics is assumed to be very fast so that combustion is primarily controlled by the rate of mixing between oxygen and hydrogen. The preburner pressure is much greater than the critical pressures of oxygen and hydrogen; hence, a gas-gas diffusion model (rather than an evaporation model) has been developed to compute the rate of interphase mixing. Empirical correlations have been incorporated to account for the effect of slip on the interphase exchange. A sensitivity study has been performed with various model parameters. It is observed that the model can predict possibility of incomplete combustion and local regions of high temperatures under steady operating conditions. Some of these anomalies have been observed in actual tests, and the numerical model is useful for understanding possible causes and remedies. At least some measurements are needed for quantitative verification of the model.

Prakash, C.↗
DOE OSTI2023

Charge decay in the spatial afterglow of plasmas and its impact on diffusion regimes

Abstract The spatial afterglow is a region at the boundary of a non-equilibrium plasma where charged species relax into ambient equilibrium. In many applications, the spatial afterglow is the part of the plasma that interacts with surfaces, such as suspended particles or a material substrate. However, compared to the bulk plasma, there has been little effort devoted to studying the properties of the spatial afterglow, and a fundamental analysis has not yet been developed. Here, we apply double Langmuir probe measurements and develop an advection-diffusion-recombination model to provide a detailed description of charged species in the spatial afterglow over a wide range of pressures, temperatures, plasma dimensions, and flow rates. We find that the density of charged species in the spatial afterglow decays by orders of magnitude, which leads to a transition from ambipolar to free diffusion. These insights can be used to explain or predict experimental observations of phenomena, such as the charging of dust grains and the dose of charged species to a biomaterial.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2000

Laminar Diffusion Flame Studies (Ground- and Space-Based Studies)

Laminar diffusion flames are of interest because they provide model flame systems that are far more tractable for analysis and experiments than more practical turbulent diffusion flames. Certainly, understanding flame processes within laminar diffusion flames must precede understanding these processes in more complex turbulent diffusion flames. In addition, many properties of laminar diffusion flames are directly relevant to turbulent diffusion flames using laminar flamelet concepts. Laminar jet diffusion flame shapes (luminous flame boundaries) have been of particular interest since the classical study of Burke and Schumann because they are a simple nonintrusive measurement that is convenient for evaluating flame structure predictions. Thus, consideration of laminar flame shapes is undertaken in the following, emphasizing conditions where effects of gravity are small, due to the importance of such conditions to practical applications. Another class of interesting properties of laminar diffusion flames are their laminar soot and smoke point properties (i.e., the flame length, fuel flow rate, characteristic residence time, etc., at the onset of soot appearance in the flame (the soot point) and the onset of soot emissions from the flame (the smoke point)). These are useful observable soot properties of nonpremixed flames because they provide a convenient means to rate several aspects of flame sooting properties: the relative propensity of various fuels to produce soot in flames; the relative effects of fuel structure, fuel dilution, flame temperature and ambient pressure on the soot appearance and emission properties of flames; the relative levels of continuum radiation from soot in flames; and effects of the intrusion of gravity (or buoyant motion) on emissions of soot from flames. An important motivation to define conditions for soot emissions is that observations of laminar jet diffusion flames in critical environments, e.g., space shuttle and space station facilities, cannot involve soot emitting flames in order to ensure that test chamber windows used for experimental observations are not blocked by soot deposits, thereby compromising unusually valuable experimental results. Another important motivation to define conditions where soot is present in diffusion flames is that flame chemistry, transport and radiation properties are vastly simplified when soot is absent, making such flames far more tractable for detailed numerical simulations than corresponding soot-containing flames. Motivated by these observations, the objectives of this phase of the investigation were as follows: (1) Observe flame-sheet shapes (the location of the reaction zone near phi=1) of nonluminous (soot free) laminar jet diffusion flames in both still and coflowing air and use these results to develop simplified models of flame-sheet shapes for these conditions; (2) Observe luminous flame boundaries of luminous (soot-containing) laminar jet diffusion flames in both still and coflowing air and use these results to develop simplified models of luminous flame boundaries for these conditions. In order to fix ideas here, maximum luminous flame boundaries at the laminar smoke point conditions were sought, i.e., luminous flame boundaries at the laminar smoke point; (3) Observe effects of coflow on laminar soot- and smoke-point conditions because coflow has been proposed as a means to control soot emissions and minimize the presence of soot in diffusion flames.

Dai, Z.↗
NASA NTRS1980

The influence of a voltage ramp on the measurement of I-V characteristics of a solar cell

For efficiency and convenience the voltage applied to a Si solar cell is often fairly rapidly driven from zero to the open circuit value typically at a common rate of 1 V per millisecond. During this time the values of current are determined as a function of the instantaneous voltage thus producing an I-V characteristic. The present paper shows that the customary expressions for the current as a function of cell parameters still remain valid provided that the diffusion length in the expression for the dark current is changed from its steady state value L to the effective diffusion length L1 given by L1 = L(1 + qV/kT.tau) to the -1/2, where V is the ramp rate considered constant and tau is the lifetime of minority carriers. This result is true to a very good approximation provided that low level injection prevails.

Von Roos, O.↗
NASA NTRS1998

EUPDF: Eulerian Monte Carlo Probability Density Function Solver for Applications With Parallel Computing, Unstructured Grids, and Sprays

The success of any solution methodology used in the study of gas-turbine combustor flows depends a great deal on how well it can model the various complex and rate controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as convective and radiative heat transfer and other phenomena. The phenomena to be modeled, which are controlled by these processes, often strongly interact with each other at different times and locations. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and evaporation in many practical combustion devices. The influence of turbulence in a diffusion flame manifests itself in several forms, ranging from the so-called wrinkled, or stretched, flamelets regime to the distributed combustion regime, depending upon how turbulence interacts with various flame scales. Conventional turbulence models have difficulty treating highly nonlinear reaction rates. A solution procedure based on the composition joint probability density function (PDF) approach holds the promise of modeling various important combustion phenomena relevant to practical combustion devices (such as extinction, blowoff limits, and emissions predictions) because it can account for nonlinear chemical reaction rates without making approximations. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on the PDF method to unstructured grids, parallel computing, and sprays. EUPDF, which was developed by M.S. Raju of Nyma, Inc., was designed to be massively parallel and could easily be coupled with any existing gas-phase and/or spray solvers. EUPDF can use an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements. The application of the PDF method showed favorable results when applied to several supersonic-diffusion flames and spray flames. The EUPDF source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, M. S.↗
DOE OSTI2023

Solute effects upon dislocation motion and recovery in Mg alloys (Final Report)

The objective the research was to develop a firmer understanding of the interactions between substitutional solute atoms and dislocations within Mg alloys. These interactions govern the absolute and relative mobilities of various types of dislocations in Mg (e.g., basal < a>, non-basal < a>, and pyramidal < c+a>). Furthermore, they have an impact on dislocation recovery processes (e.g., cross-glide, climb, rearrangement, and annihilation). Ultimately, solute-dislocation interactions strongly impact a) strain hardening, b) strain rate sensitivity, c) plastic anisotropy of textured polycrystals, d) texture evolution, e) dislocation substructures evolution, and even f) recrystallization behavior. It is broadly known that Mg alloys exhibit poor low temperature formability. Since formability is largely governed by properties (a) – (c) in this list, it is critical to develop a better understanding of solute-dislocation interactions, if one hopes to improve the situation. While the theory of static solute strengthening is well developed, especially for alloys with face centered cubic (FCC) crystal structures, there are outstanding questions related to applications to hexagonal close packed (HCP) crystal structures and dynamic strain aging (DSA) in materials of various structures. DSA has far-reaching implications for metal formability and, in the case of Mg alloys, appears to correlate with the so-called rare-earth (RE) texture which has been shown to benefit formability. The Portevin-Le Chatelier (PLC) effect, and associated negative strain rate sensitivity, occur at higher temperatures (>100℃) in Mg alloys as compared with similar Al alloys, even though they have similar melting points and solute diffusivities. Our preliminary research has shown that modern, physics-based models of DSA can be tuned to describe the behavior of Mg alloys if a rather higher activation enthalpy is assumed for cross-core diffusivity. While this partially explains the delay in DSA to higher temperatures, it is also hypothesized that this delay is due in part to the intrinsically more thermally activated (rate sensitive) nature of non-basal < a> dislocation motion which is required for macroscopic flow of Mg alloys, whereas octahedral slip in many FCC metals like aluminum is essentially athermal at room temperature. It was originally proposed to employ a combination of in-situ diffraction-based experimental characterization to validate existing theory. It was envisioned to perform in-situ transmission electron microscopy (TEM) to assess individual dislocation behavior and in-situ high-energy X-ray diffraction (HEXRD) techniques which were showing great promise for elucidating collective dislocation behavior, 2 including recovery, especially if the contributions to various forms of diffraction peak broadening (𝜂𝜂,𝜔𝜔,and 2𝜃𝜃) can be effectively integrated. Finally, it was envisioned to perform discrete dislocation dynamics (DDD) modeling approaches to aide in the interpretation of both TEM and HEXRD experiments. In the end, mechanical tests were performed on more complex Mg alloys which exhibited evidence of dynamic strain aging, and this led to the establishment of another project. Mechanical test data obtained at McMaster University served as the basis of an assessment of the applicability of Bazinski’s “stress equivalence” theory of solute strengthening to polycrystalline alloys of Mg. Although we did not succeed in applying the approach to Mg alloys, we did develop expertise with the HEXRD approach using a BCC, β-Ti alloy and demonstrated numerous new capabilities that may be applied to any polycrystalline material in collaboration with researchers at CHESS and around the world: (1) assessment of details of the elastoplastic transition (yielding) using a combination of HEXRD and full-field polycrystal plasticity modeling, (2) the first-ever experimental observation of strong stress rotation within the individual grains of a polycrystalline material, and (3) a comprehensive analysis of the grain-level dislocation density evolution based upon diffraction peak broadening along 𝜂𝜂,𝜔𝜔,and 2𝜃𝜃 directions. This final aspect allowed us to confirm that dislocations were gliding on multiple plane types and not restricted to {110} type planes, and it also provided clues as to why some grains were unloading during straining, with surprising implications for our understanding of the effects of geometrically necessary dislocations (GNDs). Finally, graduate student, Mohammed Shabana, developed a MATLAB code which confirmed the conclusions of Prof. Catalin Picu (Rensallear Polytechnic Institute, RPI) regarding the effect of solute-trapped, forest dislocations on the breaking stress of Lomer lock junctions in FCC metal alloys. He applied the same anisotropic line-tension model to a variety of dislocation junction configurations and found an inconsistency in the widely cited results of Dupuy and Fivel regarding the Hirth Lock, and he outlined an approach to extend these finding to HCP Mg alloys that we are still pursuing with discretionary fundings at UVA.

36 MATERIALS SCIENCE↗
NASA NTRS1990

Evolution of magnetized, rotating, isothermal clouds

Molecular cloud cores, in which star formation process now proceeds, are often found with elongated shape. This suggests that the core collapsed preferentially along the direction parallel to the global magnetic field and/or parallel to the cloud's initial angular momentum. Actually the magnetic field strength in the cloud has been measured recently with the Zeeman splitting (Goodman et al. 1989). The authors indicate the magnetic field of 10 to 30 mu G exists in the cloud. Further, the observation of polarization in the near IR from background stars shows that the magnetic field runs perpendicularly to the major axis of the cloud (Tamura et al. 1987). As for the angular momentum, the rotation rate of 0.2 to 6 kms s(exp -1) pc(exp -1) is reported (Goldsmith and Arquilla 1984) in 16 dark cloud regions. If the cloud collapses from the diffuse cloud with density n approx. 1 cm (exp -3) with strictly conserving the angular momentum which was shared from the galactic rotation, the rotation rate of the cloud will be x mega sub j = const approx. 3(n/1000 cm exp -3) exp 2/3 km s(exp -1) pc (exp -1) (Mouschovias 1987). The evolution of the rotating magnetized cloud is discussed here.

Tomisaka, Kohji↗
DOE OSTI2020

Polysulfide entrapment and retardation in gel electrolyte Li–S batteries: experiments and modeling

Gelling the electrolyte in lithium–sulfur (Li–S) batteries provides advantages in terms of safety and mitigating polysulfide shuttling leading to reduced side reactions and dendrite growth on the anode side. In the current study, organic and ceramic crosslinkers have been used to improve the entrapment of polysulfides and prevent them from shuttling to the anode side. However, the performance of gel electrolytes is still far from being on par with liquid electrolytes in terms of raw capacity. The discharge curves exhibit an extended first plateau for gel electrolyte batteries compared to liquid electrolytes, showing that polysulfides are better trapped in the gel system. Yet, in the conversion of soluble polysulfides into insoluble polysulfides, the second plateau was delayed. Qualitative and quantitative analysis of the polysulfide diffusion in the liquid electrolyte and gel electrolyte systems were also conducted to predict the relative diffusivity ratio between the liquid electrolyte and gel electrolyte for simulation. To probe what happens at the end of the discharge, where the mass transport in the gel electrolyte appears to be slow, a modified discharge procedure was used where the cell is first paused after discharge and then discharged a second time before recharging normally. Coin cell tests for two gel electrolytes in addition to liquid electrolyte show up to 7% capacity recovery after the pause, which indicates that the diffusion of species during the pause is responsible for the failure recovery. Additionally, the recovery is higher for gel electrolytes compared to liquid, for higher C-rates, and for longer pause times. To better understand this behavior, a Li–S numerical model with mass transport limited reactions was used to examine the polysulfide diffusion coefficients expected in different electrolyte systems. The model is able to reproduce the trends seen in experiments and yields a higher recovery for smaller Li 2 S 4 diffusion coefficients, suggesting that insufficient Li 2 S 4 mass transport is responsible for failure at the end of discharge. By identifying the effect that gel electrolytes have on the state of the battery during discharge using both experimental and simulation methods, this work is a step towards leveraging the advantages of gel electrolytes and mitigating their weaknesses.

25 ENERGY STORAGE↗
NASA NTRS1994

A Coupled Soil-Atmosphere Model of H2O2 on Mars

The Viking Gas Chromatograph Mass Spectrometer failed to detect organic compounds on Mars, and both the Viking Labeled Release and the Viking Gas Exchange experiments indicated a reactive soil surface. These results have led to the widespread belief that there are oxidants in the martian soil. Since H2O2 is produced by photochemical processes in the atmosphere of Mars, and has been shown in the laboratory to reproduce closely the Viking LR results, it is a likely candidate for a martian soil oxidant. Here, we report on the results of a coupled soil/atmosphere transport model for H202 on Mars. Upon diffusing into the soil, its concentration is determined by the extent to which it is adsorbed and by the rate at which it is catalytically destroyed. An analytical model for calculating the distribution of H202 in the martian atmosphere and soil is developed. The concentration of H202 in the soil is shown to go to zero at a finite depth, a consequence of the nonlinear soil diffusion equation. The model is parameterized in terms of an unknown quantity, the lifetime of H202 against heterogeneous catalytic destruction in the soil. Calculated concentrations are compared with a H202 concentration of 30 nmoles/cu cm, inferred from the Viking Labeled Release experiment. A significant result of this model is that for a wide range of H202 lifetimes (up to 105 years), the extinction depth was found to be less than 3 m. The maximum possible concentration in the top 4 cm is calculated to be approx. 240 nmoles/cu cm, achieved with lifetimes of greater than 1000 years. Concentrations higher than 30 nmoles/cu cm require lifetimes of greater than 4.3 terrestrial years. For a wide range of H202 lifetimes, it was found that the atmospheric concentration is only weakly coupled with soil loss processes. Losses to the soil become significant only when lifetimes are less than a few hours. If there are depths below which H202 is not transported, it is plausible that organic compounds, protected from an oxidizing environment, may still exist. They would have been deposited by meteors, or be the organic remains of past life.

Bullock, Mark A.↗