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 379 records · Page 21

Reaction of atomic bromine with acetylene and loss rate of atmospheric acetylene due to reaction with OH, Cl, O, and Br

The first-order, diffusion, and bimolecular rate constants for the reaction Br + C2H2 yields C2H3Br are evaluated. The rate constants are measured at 210, 248, 298, and 393 K and at pressures between 15-100 torr Ar using flash photolysis combined with time-resolved detection of atomic bromine via Br resonance radiation. It is observed that the reaction is not affected by pressure or temperature and the bimolecular constant = (4.0 + or - 0.8) x 10 to the -15th cu cm/sec with an error of two standard deviations. The C2H2 + Br reaction rates are compared with reactions of C2H2 with Cl, OH, NH2, and H. The loss rates for atmospheric C2H2 for reactions with OH, Cl, O, and Br are calculated as a function of altitude.

Payne, W. A.↗

Silicon alloying enhances fast heating rate combustion of aluminum particles

Alloys and composites of aluminum (Al) have shown promise in regulating and enhancing particle combustion for energy generation applications. Recent work on aluminum-silicon (Al-Si) spherical alloy particles has shown improved combustion at low heating rates through enhanced diffusion accompanying a lower melting temperature. This study extends reactivity analysis to higher heating rates, comparing oxidation of Al-Si with Al. Flame speeds of Al-Si powder mixed with molybdenum trioxide (MoO 3 ) powder (Al-Si+MoO 3 ) exhibited a faster transition to steady propagation relative to Al+MoO 3 . Bomb calorimetry experiments revealed up to 5.8% greater early temperature rise for Al-Si powder. Rapid steady propagation for Al-Si mixture and faster temperature rise for Al-Si particles were attributed to accelerated kinetics evidenced in thermal equilibrium analysis of the mixtures using a differential scanning calorimeter (DSC). Larger DSC exotherms in the early stages of oxidation (i.e., 480–720 °C) correlated with early heat release and promoted steady flame propagation for Al-Si+MoO 3 compared with Al+MoO 3 . Furthermore, in multiple heating rate DSC studies, Al-Si+MoO 3 ignited at a lower heating rate (i.e., 15 ⁰C/min) than Al+MoO 3 (i.e., 20 ⁰C/min). Both mixtures ignited at temperatures less than 625 °C which is above the melting temperature of Al-Si (574 °C) and below Al (659 °C). Thus, alloy particle fuels with enhanced diffusion-controlled kinetics promote steady flame propagation and show promise for energy generation applications. This is especially promising for technologies driven by enhanced condensed phase combustion such as some propellant additives or applications dependent on reliable burn rates such as pyrotechnic time delay formulations and primers.

36 MATERIALS SCIENCE↗

Analysis of turbulent free jet hydrogen-air diffusion flames with finite chemical reaction rates

The nonequilibrium flow field resulting from the turbulent mixing and combustion of a supersonic axisymmetric hydrogen jet in a supersonic parallel coflowing air stream is analyzed. Effective turbulent transport properties are determined using the (K-epsilon) model. The finite-rate chemistry model considers eight reactions between six chemical species, H, O, H2O, OH, O2, and H2. The governing set of nonlinear partial differential equations is solved by an implicit finite-difference procedure. Radial distributions are obtained at two downstream locations of variables such as turbulent kinetic energy, turbulent dissipation rate, turbulent scale length, and viscosity. The results show that these variables attain peak values at the axis of symmetry. Computed distributions of velocity, temperature, and mass fraction are also given. A direct analytical approach to account for the effect of species concentration fluctuations on the mean production rate of species (the phenomenon of unmixedness) is also presented. However, the use of the method does not seem justified in view of the excessive computer time required to solve the resulting system of equations.

Sislian, J. P.↗

Analysis of turbulent free-jet hydrogen-air diffusion flames with finite chemical reaction rates

A numerical analysis is presented of the nonequilibrium flow field resulting from the turbulent mixing and combustion of an axisymmetric hydrogen jet in a supersonic parallel ambient air stream. The effective turbulent transport properties are determined by means of a two-equation model of turbulence. The finite-rate chemistry model considers eight elementary reactions among six chemical species: H, O, H2O, OH, O2 and H2. The governing set of nonlinear partial differential equations was solved by using an implicit finite-difference procedure. Radial distributions were obtained at two downstream locations for some important variables affecting the flow development, such as the turbulent kinetic energy and its dissipation rate. The results show that these variables attain their peak values on the axis of symmetry. The computed distribution of velocity, temperature, and mass fractions of the chemical species gives a complete description of the flow field. The numerical predictions were compared with two sets of experimental data. Good qualitative agreement was obtained.

Sislian, J. P.↗

A kinetic model for the thermal nitridation of SiO2/Si

To explain the observed nitrogen distributions in thermally nitridated SiO2 films, a kinetic model is proposed in which the nitridation process is simulated, using the first-order chemical kinetics and Arrhenius dependence of the diffusion and reaction rates on temperature. The calculations show that initially, as the substrate reacts with diffusing nitrogen, a nitrogen-rich oxynitride forms at the SiO2-Si interface, while at nitridation temperatures above 1000 C, an oxygen-rich oxynitride subsequently forms at the interface, due to reaction of the substrate with an increasingly concentrated oxygen displaced by the slower nitridation of the SiO2. This sequence of events results in a nitrogen distribution in which the peak of the interfacial nitrogen concentration occurs away from the interface. The results are compared with the observed nitrogen distribution. The calculated results have correctly predicted the positions of the interfacial nitrogen peaks at the temperatures of 800, 1000, and 1150 C. To account for the observed width of the interfacial nitrogen distribution, it was found necessary to include in the simulations the effect of interfacial strain.

Vasquez, R. P.↗

Demonstrate Improved Ag Diffusion and Describe the Basis for Pd Penetration Modeling in SiC

In past work, an effective diffusivity coefficient was determined for Ag transport through the silicon carbide layer of a tristructural isotropic fuel particle. The effective diffusivity coefficient accounts for the microstructure of the silicon carbide and includes both bulk diffusion and grain boundary diffusion of Ag. In this report, the model has been improved by accounting for the enhanced concentration of vacancies in the bulk due to irradiation, which substantially influence bulk diffusivity at low temperatures. To improve the BISON model and make it fission rate dependent, effective diffusivity calculations have been performed that incorporate the radiation modified bulk diffusivity. The microstructure and irradiation-dependent effective diffusivity has also been implemented into BISON, and its predictions for Ag release from tristructural isotropic fuel have been successfully compared to AGR-1 post irradiation measurements. Moreover, a new feature has been developed in the Multiphysics Object-Oriented Simulation Environment (MOOSE) to account for different grain boundary types. The Ag diffusivity in 5 (210)/[001] grain boundaries has been computed and was found to be greater than in random high-angle grain boundaries. The presence of the fission product Pd can also have an important effect on the properties of the silicon carbide layer in tristructural isotropic particles. The penetration of Pd into the silicon carbide layer causes a corrosion reaction that can lead to the failure of the silicon carbide layer; however, this corrosion reaction is not well understood. To enable an improved understanding of the mechanism, ab-initio molecular dynamics simulations of Pd interaction with bulk silicon carbide have been performed. The improved understanding of the reaction will form a basis for future improvements to the BISON’s Pd penetration failure model

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

The use of the experimentally deduced Brunt-Vaisala frequency and turbulent velocity fluctuations to estimate the eddy diffusion coefficient

The determination of the turbulent energy dissipation rate or the eddy diffusion coefficient from radar observations can be done through the turbulence refractive index structure constant, deduced from calibrated echo power measurements, or through the turbulent velocity fluctuations, deduced from the echo spectrum width. Besides the radar parameters, power and spectrum width, the first approach needs knowledge of profiles of temperature and electron density in the mesosphere and the fraction of the radar volume filled with turbulence. The latter approach needs knowledge of the temperature profile, namely, the Brunt-Vaisala frequency. The use of this latter approach is demonstrated.

Rottger, J.↗

Thermophoretically enhanced mass transport rates to solid and transpiration-cooled walls across turbulent (law-of-the-wall) boundary layers

Convective-diffusion mass transfer rate predictions are made for both solid wall and transpiration-cooled 'law-of-the-wall' nonisothermal turbulent boundary layers (TBLs), including the mechanism of thermophoresis, i.e., small particle mass transport 'down a temperature gradient'. The present calculations are confined to low mass-loading situations but span the entire particle size range from vapor molecules to particles near the onset of inertial ('eddy') impaction. It is shown that, when Sc is much greater than 1, thermophoresis greatly increases particle deposition rates to internally cooled solid walls, but only partially offsets the appreciable reduction in deposition rates associated with dust-free gas-transpiration-cooled surfaces. Thus, efficient particle sampling from hot dusty gases can be carried out using transpiration 'shielded' probe surfaces.

Gokoglu, Suleyman A.↗

Simulation of polymerization induced phase separation in model thermosets

Polymerization induced phase separation (PIPS) in a three component thermoset is studied using molecular dynamics simulations of a new coarse-grained thermoset model. The system includes two crosslinker molecules, which differ in their glass transition temperatures (T g ) and chain length and thus have the potential for phase separation. Here, one crosslinker has a high T g corresponding to a rubbery behavior, and simulations were performed for a short length (4 beads) and a long length (33 beads). The resin and other crosslinker have low T g . A coarse-grained model is developed with these features and with interaction parameters determined so that for either rubbery crosslinker length, the system is in the liquid state at the cure temperature. For sufficiently slow reaction rates, the long rubbery molecule exhibits PIPS into a bicontinuous array of nanoscale domains, but the short one does not, reproducing recent experimental results. The simulations demonstrate that the reaction rates must be slow enough to allow diffusion to yield phase separation. Particularly, the reaction rate corresponding to the secondary amine must be very slow, else the structure of crosslinked clusters and the substantially increased diffusion time will prevent PIPS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Response to the comments on creep-induced elemental redistribution at grain boundaries of 304L stainless steel – An experimental evidence for diffusional creep mechanisms

Sandström and He reanalyzed our 304L creep data (Acta Materialia 121137, 2025) at 10 MPa/750 °C and proposed that: (i) primary dislocation creep dominates based on the ϕ -model; (ii) measured creep rates deviate from published 304H data and Coble creep predictions using literature grain boundary (GB) diffusivity; and (iii) creep-induced elemental redistribution (CIER) could be explained by dislocation emission/absorption (DEAA) and GB creep (GBC), as opposed to diffusional creep. Here, in this response, we propose that: (i) the creep behavior at 15 MPa/750 °C—for which we demonstrated CIER—is only partially captured by the ϕ -model; (ii) direct comparison between 304L and 304H creep data is inappropriate due to carbides in 304H, and Coble creep predictions estimated from CALPHAD-based bulk diffusivity span experimental rates; and (iii) CIER exhibits distinct fingerprint for diffusional creep characterized by anisotropy and stress orientation dependence whereas neither DEAA nor GBC have been mechanistically substantiated.

36 - MATERIALS SCIENCE↗

Quantifying the hydrodynamic contribution to electrical transport in non-Brownian suspensions

Electrical transport in semiconducting and metallic particle suspensions is an enabling feature of emerging grid-scale battery technologies. Although the physics of the transport process plays a key role in these technologies, no universal framework has yet emerged. Here, we examine the important contribution of shear flow to the electrical transport of non-Brownian suspensions. We find that these suspensions exhibit a strong dependence of the transport rate on the particle volume fraction and applied shear rate, which enables the conductivity to be dynamically changed by over 10 7 decades based on the applied shear rate. We combine experiments and simulations to conclude that the transport process relies on a combination of charge and particle diffusion with a rate that can be predicted using a quantitative physical model that incorporates the self-diffusion of the particles.

42 ENGINEERING↗

Rate of energy gain and maximum energy in diffusive shock acceleration

The problem of diffusive shock acceleration of fast charged particles is reexamined with emphasis on the rate of energy gain, and the maximum energy which can be attained in a given circumstance. The direction of the average magnetic field at the shock is shown to have a large effect. If the perpendicular diffusion coefficient is much smaller than the parallel coefficient, particles can gain much more energy if the shock is quasi-perpendicular than if it is quasi-parallel. The maximum energy attainable can be substantially higher (by a factor of 100 or more) than previous discussions would predict, in cases where the shock is quasi-perpendicular. The energy gain increases as kappa-perpendicular decreases. The principal limitation comes from the requirement that diffusion be a valid approximation to the particle motion, and that the particle be able to diffuse fast enough to encounter the shock many times.

Jokipii, J. R.↗

Mixing and the dynamics of the deep chlorophyll maximum in Lake Tahoe

Chlorophyll-temperature profiles were measured across Lake Tahoe about every 10 days from April through July 1980. Analysis of the 123 profiles and associated productivity and nutrient data identified three important processes in the formation and dynamics of the deep chlorophyll maximum (DCM): turbulent diffusion, nutrient supply rate, and light availability. Seasonal variation in these three processes resulted in three regimes: a diffusion-dominated regime with a weak DCM, a variable-mixing regime with a pronounced, nutrient supply-dominated DCM, and a stable regime with a deep, moderate light availability-dominated DCM. The transition between the first two regimes occurred in about 10 days, the transition between the last two more gradually over about 3 weeks. The degree of spatial variability of the DCM was highest in the second regime and lowest in the third. These data indicate that the DCM in Lake Tahoe is constant in neither time nor space.

Abbott, M. R.↗

Ablation, Thermal Response, and Chemistry Program for Analysis of Thermal Protection Systems

In previous work, the authors documented the Multicomponent Ablation Thermochemistry (MAT) and Fully Implicit Ablation and Thermal response (FIAT) programs. In this work, key features from MAT and FIAT were combined to create the new Fully Implicit Ablation, Thermal response, and Chemistry (FIATC) program. FIATC is fully compatible with FIAT (version 2.5) but has expanded capabilities to compute the multispecies surface chemistry and ablation rate as part of the surface energy balance. This new methodology eliminates B' tables, provides blown species fractions as a function of time, and enables calculations that would otherwise be impractical (e.g. 4+ dimensional tables) such as pyrolysis and ablation with kinetic rates or unequal diffusion coefficients. Equations and solution procedures are presented, then representative calculations of equilibrium and finite-rate ablation in flight and ground-test environments are discussed.

Milos, Frank S.↗

Water vapor diffusion membranes, 2

Transport mechanisms were investigated for the three different types of water vapor diffusion membranes. Membranes representing porous wetting and porous nonwetting structures as well as dense diffusive membrane structures were investigated for water permeation rate as a function of: (1) temperature, (2) solids composition in solution, and (3) such hydrodynamic parameters as sweep gas flow rate, solution flow rate and cell geometry. These properties were measured using nitrogen sweep gas to collect the effluent. In addition, the chemical stability to chromic acid-stabilized urine was measured for several of each type of membrane. A technology based on the mechanism of vapor transport was developed, whereby the vapor diffusion rates and relative susceptibility of membranes to fouling and failure could be projected for long-term vapor recovery trials using natural chromic acid-stabilized urine.

Holland, F. F.↗

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗

Cooling rates for lunar samples determined with a diffusion model for phosphide exsolution

A numerical model for diffusion-controlled phase growth has been applied to the exsolution of phosphide lamellae in lunar metal grains. Computer simulations reproduce observed composition profiles, and reveal the influence of cooling rate on dimensional and compositional parameters of phosphide and metal. At lower cooling rates, phosphide lamellae are larger and the concentration of P in the metal host close to the interface is lower. Cooling rates inferred for Apollo 16 samples, based on compositions and dimensions of the phosphide-metal grains, are mostly in the range 1-100 C/day. These rates correspond to burial depths of 5-0.5 m for melt rocks and 3-0.3 m for breccias. This is in good agreement with thicknesses of lunar cooling units determined by other techniques.

Hewins, R. H.↗