Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “phase kinetics”

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 487 records · Page 27

Kinetics of diffusional droplet growth in a liquid/liquid two-phase system

A new powerful experimental technique based on holographic observations, developed at the NASA Marshall Space Flight Center, now permits observation of small liquid droplets coarsening. This technique was developed and used for mixed-dimensional coarsening studies. Experiments were conducted on an isopycnic two-phase alloy of succinonitrile and water, annealed isothermally over a four-month period. The succinonitrile-rich droplets precipitate from a water-rich liquid matrix having a density very close to that of the droplets. The matrix and droplets, however, have different optical indices. The results of these experiments, along with the results of computer simulation based on the quasi-static diffusion approximation developed at Rensselaer are reported. These results were published recently. Copies of these papers are attached to this report.

Glickman, M. E.↗

A three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawn and dusk sides of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H+ fluxes at tens of keV, which are always over-estimated. A newly-invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm-time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, j(sub o)(1+Ay(exp n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

Three-dimensional ring current decay model

This work is an extension of a previous ring current decay model. In the previous work, a two-dimensional kinetic model was constructed to study the temporal variations of the equatorially mirroring ring current ions, considering charge exchange and Coulomb drag losses along drift paths in a magnetic dipole field. In this work, particles with arbitrary pitch angle are considered. By bounce averaging the kinetic equation of the phase space density, information along magnetic field lines can be inferred from the equator. The three-dimensional model is used to simulate the recovery phase of a model great magnetic storm, similar to that which occurred in early February 1986. The initial distribution of ring current ions (at the minimum Dst) is extrapolated to all local times from AMPTE/CCE spacecraft observations on the dawnside and duskside of the inner magnetosphere spanning the L value range L = 2.25 to 6.75. Observations by AMPTE/CCE of ring current distributions over subsequent orbits during the storm recovery phase are compared to model outputs. In general, the calculated ion fluxes are consistent with observations, except for H(+) fluxes at tens of keV, which are always overestimated. A newly invented visualization idea, designated as a chromogram, is used to display the spatial and energy dependence of the ring current ion differential flux. Important features of storm time ring current, such as day-night asymmetry during injection and drift hole on the dayside at low energies (less than 10 keV), are manifested in the chromogram representation. The pitch angle distribution is well fit by the function, J(sub o)(1 + Ay(sup n)), where y is sine of the equatorial pitch angle. The evolution of the index n is a combined effect of charge exchange loss and particle drift. At low energies (less than 30 keV), both drift dispersion and charge exchange are important in determining n.

Fok, Mei-Ching↗

Development of an Advanced Computational Model for OMCVD of Indium Nitride

An advanced computational model is being developed to predict the formation of indium nitride (InN) film from the reaction of trimethylindium (In(CH3)3) with ammonia (NH3). The components are introduced into the reactor in the gas phase within a background of molecular nitrogen (N2). Organometallic chemical vapor deposition occurs on a heated sapphire surface. The model simulates heat and mass transport with gas and surface chemistry under steady state and pulsed conditions. The development and validation of an accurate model for the interactions between the diffusion of gas phase species and surface kinetics is essential to enable the regulation of the process in order to produce a low defect material. The validation of the model will be performed in concert with a NASA-North Carolina State University project.

Cardelino, Carlos A.↗

Upward Flame Spread Over Thin Solids in Partial Gravity

The effects of partial-gravity, reduced pressure, and sample width on upward flame spread over a thin cellulose fuel were studied experimentally and the results were compared to a numerical flame spread simulation. Fuel samples 1-cm, 2-cm, and 4-cm wide were burned in air at reduced pressures of 0.2 to 0.4 atmospheres in simulated gravity environments of 0.1-G, 0.16-G (Lunar), and 0.38-G (Martian) onboard the NASA KC-135 aircraft and in normal-gravity tests. Observed steady flame propagation speeds and pyrolysis lengths were approximately proportional to the gravity level. Flames spread more quickly and were longer with the wider samples and the variations with gravity and pressure increased with sample width. A numerical simulation of upward flame spread was developed including three-dimensional Navier-Stokes equations, one-step Arrhenius kinetics for the gas phase flame and for the solid surface decomposition, and a fuel-surface radiative loss. The model provides detailed structure of flame temperatures, the flow field interactions with the flame, and the solid fuel mass disappearance. The simulation agrees with experimental flame spread rates and their dependence on gravity level but predicts a wider flammable region than found by experiment. Some unique three-dimensional flame features are demonstrated in the model results.

Feier, I. I.↗

Finite Element Models for Electron Beam Freeform Fabrication Process

Electron beam freeform fabrication (EBF3) is a member of an emerging class of direct manufacturing processes known as solid freeform fabrication (SFF); another member of the class is the laser deposition process. Successful application of the EBF3 process requires precise control of a number of process parameters such as the EB power, speed, and metal feed rate in order to ensure thermal management; good fusion between the substrate and the first layer and between successive layers; minimize part distortion and residual stresses; and control the microstructure of the finished product. This is the only effort thus far that has addressed computer simulation of the EBF3 process. The models developed in this effort can assist in reducing the number of trials in the laboratory or on the shop floor while making high-quality parts. With some modifications, their use can be further extended to the simulation of laser, TIG (tungsten inert gas), and other deposition processes. A solid mechanics-based finite element code, ABAQUS, was chosen as the primary engine in developing these models whereas a computational fluid dynamics (CFD) code, Fluent, was used in a support role. Several innovative concepts were developed, some of which are highlighted below. These concepts were implemented in a number of new computer models either in the form of stand-alone programs or as user subroutines for ABAQUS and Fluent codes. A database of thermo-physical, mechanical, fluid, and metallurgical properties of stainless steel 304 was developed. Computing models for Gaussian and raster modes of the electron beam heat input were developed. Also, new schemes were devised to account for the heat sink effect during the deposition process. These innovations, and others, lead to improved models for thermal management and prediction of transient/residual stresses and distortions. Two approaches for the prediction of microstructure were pursued. The first was an empirical approach involving the computation of thermal gradient, solidification rate, and velocity (G,R,V) coupled with the use of a solidification map that should be known a priori. The second approach relies completely on computer simulation. For this purpose a criterion for the prediction of morphology was proposed, which was combined with three alternative models for the prediction of microstructure; one based on solidification kinetics, the second on phase diagram, and the third on differential scanning calorimetry data. The last was found to be the simplest and the most versatile; it can be used with multicomponent alloys and rapid solidification without any additional difficulty. For the purpose of (limited) experimental validation, finite element models developed in this effort were applied to three different shapes made of stainless steel 304 material, designed expressly for this effort with an increasing level of complexity. These finite element models require large computation time, especially when applied to deposits with multiple adjacent beads and layers. This problem can be overcome, to some extent, by the use of fast, multi-core computers. Also, due to their numerical nature coupled with the fact that solid mechanics- based models are being used to represent the material behavior in liquid and vapor phases as well, the models have some inherent approximations that become more pronounced when dealing with multi-bead and multi-layer deposits.

Chandra, Umesh↗

Sintering Characteristics of Multilayered Thermal Barrier Coatings Under Thermal Gradient and Isothermal High Temperature Annealing Conditions

Pyrochlore oxides have most of the relevant attributes for use as next generation thermal barrier coatings such as phase stability, low sintering kinetics and low thermal conductivity. One of the issues with the pyrochlore oxides is their lower toughness and therefore higher erosion rate compared to the current state-of-the-art TBC material, yttria (6 to 8 wt%) stabilized zirconia (YSZ). In this work, sintering characteristics were investigated for novel multilayered coating consisted of alternating layers of pyrochlore oxide viz Gd2Zr2O7 and t' low k (rare earth oxide doped YSZ). Thermal gradient and isothermal high temperature (1316 C) annealing conditions were used to investigate sintering and cracking in these coatings. The results are then compared with that of relevant monolayered coatings and a baseline YSZ coating.

Thermal Barrier Coatings↗

Design Study for Bio-Oil Production from Biomass Using a Dual Fluidized-Bed Reactor

To evaluate the re-design and reconfiguration of a dual-fluidized bed (DFB) gasification system into a recirculating pyrolysis reactor, Computation Fluid Dynamic (CFD) simulations of the system were conducted. The Barracuda Virtual Reactor® computational particle fluid dynamic code was used to perform simulations of the pyrolysis process. Modeling of the chemical reaction kinetics for both gas phase and solid particle phase were included. The recirculating pyrolysis reactor shown in Fig. 1a is based on a bubbling-bed biomass pyrolyzer and a riser combustor to convert the remaining char. The operational differences in the re-configuration of the DFB gasification system into a recirculating pyrolysis system for the production of bio-oil are (1) replacement of a low-surface area inert bed material with a high-surface-area bed material that has acidic properties to provide catalytic activity for the production of bio-oil with reduced oxygen content, (2) lower temperature and residence time for bio-oil production from pyrolysis, (3) replacement of the fluidization gas in the bubbling bed pyrolyzer from steam to nitrogen, and (4) the reduction of pyrolyzer freeboard volume. The bed material used for catalytic pyrolysis is Sasol 300 (300-micron dia., bulk density 0.94 kg/l, and surface area 130 m2/g) and is a theta-alumina with mild acidity. This is in comparison with previous standard bed material Carbo HSP (430-micron dia., bulk density 2.01 kg/l, and surface area 0.03 m2/g) used for gasification. For bio-oil production, pyrolysis in the bubbling bed requires temperatures in the range of 550 C in comparison with gasification temperatures near 850 C. To attain this lower temperature requires management of the energy mass balances, with control of the bed material recirculation rate between bubbling bed pyrolyzer and riser combustor, the introduction of a nitrogen purge in the pyrolyzer, and adjusting the pressure balance between the two vessels. To extract bio-oil from the pyrolysis reactor with a snorkel, two different freeboard configurations were evaluated. In Fig. 1 b the existing high freeboard configuration is shown and in Fig. 1 c the reduced freeboard design is presented. The introduction of a nitrogen purge for the high free board configuration provided the highest bio-oil production from the CFD simulations. To decrease the bed-material circulation rate, primary and secondary air on the combustor side were reduced, and the pressure on the combustor was slightly increased. A portion of the biomass and bio-oil was observed to be transported to the combustor, leading to a smaller pyrolysis yield. The control of the pyrolyzer temperature is performed by controlling the circulation rate. Good fluidization of the bubbling bed and cascade PID control are required to keep the temperature from oscillating due to large time delay experienced when changing primary and secondary air. (a) (b) (c) Figure 1. (a) Dual Fluidized Bed pyrolysis system configuration, (b) mole fraction of nonpolar biooil in high freeboard configuration, and (c) mole fraction of nonpolar bio-oil in low freeboard configuration with nitrogen purge introduced in both configurations.

BASIC BIOLOGICAL SCIENCES,BIOMASS FUELS↗

BETO 2021 Peer Review - Thermochemical Platform Analysis WBS: 2.1.0.302

The objective of the NREL Thermochemical Platform Analysis (WBS 2.1.0.302) project is to inform and guide R&D priorities for thermal and catalytic conversion processes by providing process design and techno-economic analysis (TEA). This is achieved through close collaboration with researchers and external experts, along with the use of both commercially available modeling tools and the development or use of collaboration-derived domain-specific tools and resources, such as refinery integration, kinetic and reactor models, phase equilibrium models, and pertinent bio-products market studies. This project is directly aligned with DOE-BETO goals, with the enabling of technology advancements and cost reduction for biomass derived biofuels being one of its primary objectives. TEA-guided research facilitated by this project has helped achieve significant modeled cost reductions for the ex situ catalytic fast pyrolysis (CFP) pathway and the indirect liquefaction (IDL) pathway for the conversion of syngas to high-octane gasoline (HOG). Cost reduction through refinery integration, development of valuable co-products, and other options are being identified for future research to help reduce the modeled MFSP to $2.50/GGE by 2030. Additional priorities anticipated in the future, such as the use of renewable electricity for liquid fuels and products, and emphasis on waste utilization are also being explored in conjunction with research on catalytic utilization of syngas and other gases (including CO2). Industry-relevant parameters are given deliberate attention as part of the work done under this project to help answer questions important for future commercialization and address associated risks.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Conversion Kinetics and Ionic Conductivity in Na-β”-Alumina + YSZ (Naβ”AY) Sodium Solid Electrolyte via Vapor Phase Conversion Process

Sodium ion batteries have been receiving increasing attention and may see potential revival in the near future, particularly in large-scale grid energy storage coupling with wind and solar power generation, due to the abundant sodium resources, low cost, and sufficiently high energy density. Among the known sodium ion conductors, the Na-β”-alumina electrolyte remains highly attractive because of its high ionic conductivity. This study focuses on the vapor phase synthesis of a Na-β”-Alumina + YSZ (Naβ”AY) composite sodium electrolyte, which has higher mechanical strength and stability than conventional single phase β”-Alumina. The objectives are the measurement of conversion kinetics through a newly developed weight-gain based model and the determination of sodium ionic conductivity in the composite electrolyte. Starting samples contained ~70 vol% α-Alumina and ~30 vol% YSZ (3 mol% Y 2 O 3 stabilized Zirconia) with and without a thin alumina surface layer made by sintering in air at 1600 °C. The sintered samples were placed in a powder of Na-β”-alumina and heat-treated at 1250 °C for various periods. Sample dimensions and weight were measured as a function of heat treatment time. The conversion of α-Alumina in the α-Alumina + YSZ composite into Naβ”AY occurred by coupled diffusion of sodium ions through Na-β”-alumina and of oxygen ions through YSZ, effectively diffusing Na2O. From the analysis of the time dependence of sample mass and dimensions, the effective diffusion coefficient of Na 2 O through the sample, D eff , was estimated to be 1.74 x 10 -7 cm2 s -1 , and the effective interface transfer parameter, k eff , was estimated as 2.33 x 10 -6 cm s -1 . By depositing a thin alumina coating layer on top of the bulk composite, the chemical diffusion coefficient of oxygen through single phase Na-β”-alumina was estimated as 4.35 x 10 -10 cm 2 s -1 . An AC impedance measurement was performed on a fully converted Naβ”AY composite, and the conductivity of the composite electrolyte was 1.3 x 10 -1 S cm -1 at 300 °C and 1.6 x 10 -3 S cm -1 at 25 °C, indicating promising applications in solid state or molten salt batteries at low to intermediate temperatures.

36 MATERIALS SCIENCE↗

Phase diagram of the three-dimensional subsystem toric code

Subsystem quantum error-correcting codes typically involve measuring a sequence of noncommuting parity check operators. They can sometimes exhibit greater fault tolerance than conventional codes, which use commuting checks. However, unlike subspace codes, it is unclear if subsystem codes—in particular their advantages—can be understood in terms of ground-state properties of a physical Hamiltonian. In this paper, we address this question for the three-dimensional subsystem toric code (3D STC), as recently constructed by Kubica and Vasmer [], which exhibits single-shot error correction. Motivated by a conjectured relation between single-shot properties and thermal stability, we study the zero- and finite-temperature phases of an associated noncommuting Hamiltonian. By mapping the Hamiltonian model to a pair of 3D Z 2 gauge theories coupled by a kinetic constraint, we find various phases at zero temperature, all separated by first-order transitions: There are 3D toric code-like phases with deconfined point-like excitations in the bulk, and there are phases with a confined bulk supporting a 2D toric code on the surface when appropriate boundary conditions are chosen. The latter is similar to the surface topological order present in 3D STC. However, the similarities between the single-shot correction in 3D STC and the confined phases are only partial: they share the same sets of degrees of freedom, but they are governed by different dynamical rules. Instead, we argue that the process of single-shot error correction can more suitably be associated with a path (rather than a point) in the zero-temperature phase diagram, a perspective, which inspires alternative measurement sequences enabling single-shot error correction. Moreover, since none of the above-mentioned phases survives at nonzero temperature, the single-shot error-correction property of the code does not imply thermal stability of the associated Hamiltonian phase. Published by the American Physical Society 2024

Li, Yaodong (ORCID:0000000337421944)↗

A Coupled Theoretical and Experimental Approach to Elucidating the Mechanisms of Methyl Esters

The goal of this project was to probe fundamental kinetics questions regarding the gas phase reactive behavior of oxygenates. In order to succeed, our program necessitated both the development of new experimental tools and the development of theory based kinetic models to explain the chemistry of oxygenate fuels. By using advanced hybrid additive manufacturing/ traditional manufacturing techniques to create specialized ceramic micro-reactors, the team was able to simulate extreme high-temperature environments with much greater precision and durability than was previously possible. These tiny, high-tech tubes allow researchers to capture and identify "fleeting" chemical species—molecules that exist for only a fraction of a millisecond—using sophisticated light sources and mass spectrometry. Additionally, the team supported the development of a tabletop VUV laser system for isomer detection without a synchrotron. Through this work, we have successfully mapped out the specific chemical pathways of various oxygen-rich fuels, solving long-standing mysteries about how these substances break down kinetically. Ultimately, this research provides the fundamental knowledge needed to design next-generation engines and fuels that are better for the environment.

09 BIOMASS FUELS↗

Photochemistry chemical kinetics in the interstellar medium.

Gas phase processes involving interstellar molecules are considered. The primary process is photodissociation in the ambient radiation field. In obscured regions, a combination of Johnson's (1965) and Stecher's (1969) Perseus extinction curves was used to obtain the radiation field as a function of the extinction at 5,500 A. On the basis of laboratory measurements, it is concluded that the molecules of water, ammonia, methane, and formaldehyde have lifetimes less than one hundred years.

Donn, B.↗

Calcium requirement of phytochrome-mediated fern-spore germination: no direct phytochrome-calcium interaction in the phytochrome-initiated transduction chain

Phytochrome-mediated germination of fern spores of Dryopteris paleacea Sw. was initiated by a saturating red-light (R) irradiation after 20 h of imbibition. For its realization external Ca2+ was required, with a threshold at a submicromolar concentration, and an optimum was reached around 10(-4) M. At concentrations > or = 10(-1) M only a reduced response was obtained, based probably on an unspecific osmotic or ionic effect. The germination response was inhibited by La3+, an antagonist of Ca2+. From these results it is concluded that Ca2+ influx from the medium into the spores may be an important event in phytochrome-mediated germination. In the absence of Ca2+ the R-stimulated system remained capable of responding to Ca2+, added as late as 40 h after R. Moreover, Ca2+ was effective even if added after the active form of phytochrome, Pfr, had been abolished by far-red (FR) 24 h after R. Thus, the primary effect of Pfr, that initiates the transduction chain, does not require calcium. "Coupling" of Pfr to subsequent dark reactions has been investigated by R-FR irradiations with various dark intervals. The resulting "escape kinetics" were characterized by a lag phase (6 h) and half-maximal escape from FR reversibility (19 h). These kinetics were not significantly changed by the presence or absence of calcium. Thus, direct interaction of Pfr and calcium is not a step in the transduction chain initiated by the active form of phytochrome.

Non-NASA Center↗

A critical analysis of U-Pu-Zr phase transitions using calorimetric, microstructural, and phase equilibria data

Metallic fuels consisting primarily of uranium, plutonium, and zirconium (U-Pu-Zr) are a leading material candidate for fast-spectrum nuclear reactors. Early demonstration programs proved the principle of safe and efficient fast reactor operation, however there is still considerable uncertainty regarding the phase equilibria and microstructural evolution across the ternary composition space. Quantitative phase formation and identification measurements are scarce and often incomplete, with studies reporting either phase transition temperatures or phase identification data, but not both from the same specimens. In this study, we critically compared experimental and calculated phase transition data and correlated with the microstructure and phase characterization data of as-cast and annealed U-Pu-Zr alloys. Differential scanning calorimetry (DSC) was used to measure phase transitions in the subsolidus regions (723−948 K) of three ternary U-Pu-Zr alloys with similar plutonium concentrations but various U/Zr ratios. Due to sluggish kinetics and narrow ranges of phase stability, complex peaks required the use of a Frazier-Suzuki peak fitting algorithm to deconvolute and calculate transition peak temperatures and enthalpies. We also identified trends of phase transition behavior by critically comparing our DSC data with previous phase transition measurements as well as historical and calculated phase equilibrium diagrams. In conclusion, this provides a critical approach for benchmarking and assessing the quality of new U-Pu-Zr phase equilibria data prior to its incorporation into nuclear material databases.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Hardenability and microstructural evolution of a precipitation strengthened Ni 50 Ti 21 Hf 25 Al 4 alloy

NiTi-based quaternary alloys are used in a variety of mechanical components, such as bearings, actuators, and dampers, owing to their good hardenability, wear resistance, and corrosion resistance. Additionally, one of the most notable characteristics of NiTi-based alloys is their shape memory effect and pseudoelastic properties. Connecting the macroscopic processing parameters employed in the design of new intermetallic alloys to the nanoscale structural characteristics dictating their behavior is crucial for improving their mechanical properties and expanding the spectrum of potential applications. Here, in this work, an arc melted Ni 50 Ti 21 Hf 25 Al 4 (at%) alloy was solution treated at 1050 °C followed by quenching and aging at 600 °C to investigate the effect of aging time on the microstructure and mechanical properties. Two types of nano-sized precipitates were observed and determined as face-centered orthorhombic H-phase (TiHf)Ni and L2 1 Heusler precipitates Ni 2 TiAl. The morphology and orientation of the H-phase were investigated using scanning and transmission electron microscopy (SEM and TEM), elucidating the coarsening kinetics and strengthening contribution of that phase to the intermetallic mechanical behavior. Following coarsening, the presence of Heusler nanoprecipitates was detected under overaged conditions through TEM imaging and nanobeam electron diffraction patterns. A peak hardness condition of 756 HV was achieved after 70 h of aging, indicating that the co-precipitation of H-phase and Heusler precipitates through a well-designed aging treatment can lead to optimal mechanical performance, thus elevating the alloy’s potential as a viable material for industrial applications.

36 MATERIALS SCIENCE↗

Metastable phase formation in undercooled liquid lead alloys

During non-equilibrium processes metastable phases are formed instead of stable phases due to the operation of various kinetic or structural constraints. By removing the most effective nucleation sites for the stable phase using emulsified droplet samples, stable phase formation can be prohibited in a broad range of undercooling and the phase space available in the metastable regime can be mapped out. With this method the thermodynamic properties of the undercooled liquid and the metastable phase boundaries corresponding to reversible metastable (solid + liquid) phase equilibria are examined from experimental data. The analysis reveals important implications for the nucleation kinetics and the choice of the kinetically most favored solidification path.

Fecht, Hans J.↗

The influence of cooling rate on condensation of iron, aluminum, and uranium oxide nanoparticles

Fundamental observations of particle size distributions are needed to develop models that predict the fate and transport of radioactive materials in the atmosphere following a nuclear incident. The extent of material transport is influenced by the time scales of particle formation processes (e.g., condensation, coagulation). In this study, we investigated the influence of cooling time scales on size distributions of uranium, aluminum, and iron oxide particles that are synthesized separately under identical run conditions inside the controlled environment of an argon plasma flow reactor. Two distinct temperature distributions are imposed along the flow reactor by varying the argon flow rate downstream of the plasma torch. The vaporized reactants of uranium, aluminum, and iron are cooled from about 5000K to 1000K before they are collected on silicon wafers for ex situ scanning electron microscope analysis. The microscope images show that the sizes of the largest aluminum and iron oxide particles heavily depend on the cooling time scales, whereas significant size variation with cooling rate is not observed for uranium oxide particles. In addition, the size distribution of aluminum oxide particles exhibits the broadest range among all three metal oxides studied. We performed simulations of particle size distributions using a kinetic model that couples gas phase oxidation chemistry with particle formation processes, including nucleation, condensation, and coagulation. The model results demonstrate the strong sensitivity of particle size distribution to different cooling histories (i.e., temperature vs residence time) along the flow reactor. In conclusion, the kinetic model also helps identify directions for future research to improve the predictions.

36 MATERIALS SCIENCE↗