SEARCH · Search NASA

Results for “Kinetics Parameters”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

DOE OSTI2014

Modeling kinetics and equilibrium of membranes with fields: Milestoning analysis and implication to permeation

Coarse graining of membrane simulations by translating atomistic dynamics to densities and fields with Milestoning is discussed. The space of the membrane system is divided into cells and the different cells are characterized by order parameters presenting the number densities. The dynamics of the order parameters are probed with Milestoning. The methodology is illustrated here for a phospholipid membrane system (a hydrated bilayer of DOPC (1,2-dioleoyl-sn-glycero-3-phosphocholine) lipid molecules). Significant inhomogeneity in membrane internal number density leads to complex free energy landscape and local maps of transition times. Dynamics and distributions of cavities within the membrane assist the permeation of nonpolar solutes such as xenon atoms. It is illustrated that quantitative and detailed dynamics of water transport through DOPC membrane can be analyzed using Milestoning with fields. The reaction space for water transport includes at least two slow variables: the normal to the membrane plane, and the water density.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2013

Study on the kinetics and mechanism of Am(III) with TODGA

The removal of minor actinides (mainly Am, Cm and Np) from the high-level waste would remarkably reduce the necessary storage time. Recently developed, a tridentate ligand, N,N,N',N'-tetraoctyl diglycolamide (TODGA), has been identified as one of the most powerful extractants being considered for the partitioning of trivalent actinides and lanthanides. In the present paper, the extraction kinetics behavior of Am(III) in TODGA/n-dodecane-HNO{sub 3} system has been studied, dependent on various parameters such as the stirring speed, the interfacial area, the extractant concentration in n-dodecane, the extracted ions concentration, the acidity of aqueous phase and the temperature. The results show that: the extraction process is controlled by diffusion mode under 130 rpm of stirring speed and by chemical reaction mode above 150 rpm. The extraction rate equation is given in the paper and the apparent extraction rate constant of Am(III) by TODGA/n-dodecane in 170 rpm and at 25 Celsius degrees is k=(24.2±3.4)*10{sup -3} mol{sup -2.18}*L{sup 2.18}*min{sup -1}*cm. The initial extraction rate of Am(III) by TODGA/n-dodecane at 170 r/min and 25 Celsius degrees increases with increasing temperature, and the apparent active energy is calculated to be E{sub a}=(25.94±0.98) kJ/mol.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2022

Automatically generated model for light alkene combustion

Light alkenes are common combustion intermediates for a variety of fuels. Therefore, understanding their oxidation and pyrolysis chemistry is key to building detailed mechanisms for heavier fuels. This work was focused on the development and evaluation of a detailed kinetic mechanism suitable for the combustion of light alkenes up to C 4 without the use of tuned parameters, instead, the parameter values come from first principles or direct measurements. The generated mechanism accurately estimates the laminar burning velocity (S u ) and ignition delay time (IDT) of light alkenes available in the literature, which represent fundamental combustion properties at a wide range of conditions. Because each parameter is thought to have a physically realistic value, not tuned to these measurements, the new model could be used as a sub-mechanism in models for other applications. The reaction network was generated with the open-source Reaction Mechanism Generator (RMG) software. Sensitivity analyses were performed under wide ranges of temperatures and pressures, allowing for the identification of the most impactful species and reactions. Based on these, a comprehensive thermochemistry database, including calculations on 550 molecules performed in this work at the CBS-QB3 level of theory, and a kinetic library, including theoretically-derived reaction rates retrieved from the literature, were built and used in the mechanism generation. The developed mechanism was compared against several existing detailed kinetic mechanisms for ethene, propene, 1-butene, 2-butene, and isobutene. Here, the newly generated model is the most accurate among the ones analyzed, in terms of fractional bias and normalized mean square error. Hence, this new model was used to analyze the chemistry of alkene combustion. Key rate coefficients were compared, to identify the cause of deviations between the models and possible areas for further improvements.

09 BIOMASS FUELS↗

DOE OSTI2012

Predictive process simulation of cryogenic implants for leading edge transistor design

Two cryogenic implant TCAD-modules have been developed: (i) A continuum-based compact model targeted towards a TCAD production environment calibrated against an extensive data-set for all common dopants. Ion-specific calibration parameters related to damage generation and dynamic annealing were used and resulted in excellent fits to the calibration data-set. (ii) A Kinetic Monte Carlo (kMC) model including the full time dependence of ion-exposure that a particular spot on the wafer experiences, as well as the resulting temperature vs. time profile of this spot. It was calibrated by adjusting damage generation and dynamic annealing parameters. The kMC simulations clearly demonstrate the importance of the time-structure of the beam for the amorphization process: Assuming an average dose-rate does not capture all of the physics and may lead to incorrect conclusions. The model enables optimization of the amorphization process through tool parameters such as scan speed or beam height.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Producing ENDF/B-quality Evaluations of 239 Pu(n,f) and 235 U(n,f) Average Prompt Neutron Multiplicities using the CGMF Model

This report documents evaluations of the neutron-induced 239 Pu and 235 U average prompt neutron multiplicities, $\overline{v}_p$, using the CGMF model developed at LANL. These evaluations are not updates of existing ENDF/B-VIII.0 nuclear data, but were both re-done from "scratch". That is, all experimental data were extracted from EXFOR, re-normalized to the newest nuclear data representing monitor observables, and detailed uncertainties were estimated from information in EXFOR as well as from templates of expected measurement uncertainties. We also included experimental data that were not yet available for ENDF/B-VIII.0, for instance, the data of Marini et al. for 239 Pu(n,f) $\overline{v}_p$ or the data of Khoklov et al. for 235 U(n,f) $\overline{v}_p$. Finally, we include the CGMF model; it computes $\overline{v}_p$ based on model parameters that at the same time calculate fission fragments as a function of mass, Y (A), the average total kinetic energy as a function of incident-neutron energy, $\langle TKE \rangle$(E inc ), etc. Such a detailed fission model that ties together many fission quantities has not been used to date for any $\overline{v}_p$ evaluation in ENDF/B. Here, we show that we can get evaluated 235 U(n,f) and 239 Pu $\overline{v}_p$ that not only correspond well to experimental $\overline{v}_p$, but that the associated evaluated model parameters also yield parameterizations of Y (A), $\langle TKE \rangle$(E inc ), etc., that correspond well to their respective experimental data. The new evaluated 239 Pu(n,f) $\overline{v}_p$ was also combined into a 239 Pu test file with the newest nuclear data for the 239 Pu prompt-fission neutron spectrum and fission cross sections. This new 239 Pu file performed reasonably well in predicting PU-MET-FAST assemblies, reaction rates in Jezebel and Flattop and three LLNL pulsed spheres. Due to that, the new 239 Pu(n,f) $\overline{v}_p$ evaluation presented here is being considered for ENDF/B-VIII.1. Hence, we can show here that CGMF is able to produce ENDF/B quality $\overline{v}_p$ nuclear data.

239Pu, 235U, Average Prompt-fission Neutron Multip↗

DOE OSTI2024

In-Situ Species Concentration Measurements In Ammonia-Mix Flames Using Ftir Spectroscopy

Hydrogen and ammonia represent two carbon-free fuel sources that could be used in place of current fossil energy sources in combustion systems. To develop optimized ammonia combustion systems, validated modeling tools are needed. In the open literature, it has been shown that the complex chemistry associated with fuel-bound nitrogen contained in ammonia differs greatly from natural gas or hydrogen combustion. As a result, several new chemical kinetic mechanisms have been developed. Many of these mechanisms have been validated experimentally, however this has primarily focused on bulk parameters such as laminar flame speed and ignition delay time. Critically, high quality measurements of species concentrations are needed under controlled conditions which are easily represented by simple models. In this paper, direct, in-situ measurements of species concentrations and gas temperature are performed in a laminar flat-flame burner. This arrangement enables comparison with 1D model predictions, better isolating chemical kinetics from the fluid dynamics. Quantitative species concentrations are determined by absorption spectroscopy using an FTIR spectrometer. Fuel compositions representative of cracked ammonia (NH 3 /H 2 ) and ammonia-natural gas (NH 3 /CH 4 ) are considered for rich and lean equivalence ratios. A major focus of the paper is on the selection of spectral features for nitric oxide and ammonia and correcting for large amounts of baseline H 2 O absorption.

36 MATERIALS SCIENCE↗

DOE OSTI2024

In-Situ Species Concentration Measurements in Ammonia-Mix Flames Using FTIR Spectroscopy

Hydrogen and ammonia represent two carbon-free fuel sources that could be used in place of current fossil energy sources in combustion systems. To develop optimized ammonia combustion systems, validated modeling tools are needed. In the open literature, it has been shown that the complex chemistry associated with fuel-bound nitrogen contained in ammonia differs greatly from natural gas or hydrogen combustion. As a result, several new chemical kinetic mechanisms have been developed. Many of these mechanisms have been validated experimentally, however this has primarily focused on bulk parameters such as laminar flame speed and ignition delay time. Critically, high quality measurements of species concentrations are needed under controlled conditions which are easily represented by simple models. In this paper, direct, in-situ measurements of species concentrations and gas temperature are performed in a laminar flat-flame burner. This arrangement enables comparison with 1D model predictions, better isolating chemical kinetics from the fluid dynamics. Quantitative species concentrations are determined by absorption spectroscopy using an FTIR spectrometer. Fuel compositions representative of cracked ammonia (NH3/H2) and ammonia-natural gas (NH3/CH4) are considered for rich and lean equivalence ratios. A major focus of the paper is on the selection of spectral features for nitric oxide and ammonia and correcting for large amounts of baseline H2O absorption.

Bedick, Clinton↗

DOE OSTI2020

Applying the Gyrokinetic Formulation of Magnetized Plasma Dynamics to Z-pinch plasmas

This project adapts and applies the gyrokinetic Eulerian code COGENT to an emerging area of highly-magnetized Z-pinch plasma systems. Simulation efforts for these plasmas currently use either magneto-hydrodynamic (MHD) fluid codes, which suffer from fidelity being restricted to narrow plasma conditions, or fully kinetic particle-in-cell (PIC) methods, which suffer from the run-time needed to resolve the very small spatial and temporal scales required for the plasma conditions. The progress with achieving detailed theoretical understanding of Z-pinch plasma systems can be greatly expedited using some combination of gyrokinetic and drift-ordered MHD models, developed in the tokamak community and implemented in the COGENT code. These formulations are applicable to a much larger parameter space than standard MHD models, while providing a substantial speed-up in run-time over fully kinetic codes. In particular, the gyrokinetic formalism is successfully applied to a flow-stabilized Z pinch (FSZP) system – a promising fusion concept where a sheared axial flow is used to provide stability against large-scale ideal-MHD modes. Results of gyrokinetic simulations with COGENT agree well with those from fully kinetic simulations, while requiring ~500x less run-time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2014

Hybrid Pressure Retarded Osmosis-Membrane Distillation System for Power Generation from Low-Grade Heat: Thermodynamic Analysis and Energy Efficiency

We present a novel hybrid membrane system that operates as a heat engine capable of utilizing low-grade thermal energy, which is not readily recoverable with existing technologies. The closed-loop system combines membrane distillation (MD), which generates concentrated and pure water streams by thermal separation, and pressure retarded osmosis (PRO), which converts the energy of mixing to electricity by a hydro-turbine. The PRO-MD system was modeled by coupling the mass and energy flows between the thermal separation (MD) and power generation (PRO) stages for heat source temperatures ranging from 40 to 80 degrees C and working concentrations of 1.0, 2.0, and 4.0 mol/kg NaCl. The factors controlling the energy efficiency of the heat engine were evaluated for both limited and unlimited mass and heat transfer kinetics in the thermal separation stage. In both cases, the relative flow rate between the MD permeate (distillate) and feed streams is identified as an important operation parameter. There is an optimal relative flow rate that maximizes the overall energy efficiency of the PRO-MD system for given working temperatures and concentration. In the case of unlimited mass and heat transfer kinetics, the energy efficiency of the system can be analytically determined based on thermodynamics. Our assessment indicates that the hybrid PRO-MD system can theoretically achieve an energy efficiency of 9.8% (81.6% of the Carnot efficiency) with hot and cold working temperatures of 60 and 20 degrees C, respectively, and a working solution of 1.0 M NaCl. When mass and heat transfer kinetics are limited, conditions that more closely represent actual operations, the practical energy efficiency will be lower than the theoretically achievable efficiency. In such practical operations, utilizing a higher working concentration will yield greater energy efficiency. Overall, our study demonstrates the theoretical viability of the PRO-MD system and identifies the key factors for performance optimization.

Lin, SH↗

DOE OSTI2022

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. Here, the modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

In Vitro Sugar Interference Testing With Amperometric Glucose Oxidase Sensors

Background: Electrochemical enzymatic glucose sensors are intended to measure blood or interstitial fluid glucose concentrations. One class of these glucose sensors are continuous glucose monitors (CGMs), indicated for tracking and trending of glucose concentrations in interstitial fluid and as an adjunct to blood glucose testing. Currently approved CGMs employ a glucose oxidase (GOx) electrochemical detection scheme. Potential interfering agents can impact the accuracy of results obtained by glucose sensors, including CGMs. Methods: Seven sugars, seven sugar alcohols, and three artificial sweeteners were in vitro screened for interference with amperometric glucose oxidase (GOx) sensors at concentrations greater than physiologic concentrations. Galactose was investigated further at physiologically relevant concentrations using a custom amperometric system. Furthermore, glucose and galactose calibration experiments were conducted to facilitate multiple enzyme kinetic analysis approaches (Michaelis-Menten and Hill equation) to understand the potential source and mechanism of interference from galactose. Results: Under in vitro testing, except for galactose, xylose and mannose, all screened compounds exhibited interference bias, expressed in mean absolute relative difference (MARD), of ⩽ 20% even at concentrations significantly higher than normal physiologic concentrations. Galactose exhibited, CGM-dependent, MARD of 47-72% and was subjected to further testing. The highest recorded mean relative difference (MRD) was 6.9 ± 1.3% when testing physiologically relevant galactose concentrations (0.1-10 mg/dL). Enzyme kinetic analysis provided calculations of maximum reaction rates ([Formula: see text]), apparent Michaelis constants ([Formula: see text]), and Hill equation h parameters for glucose and galactose substrates for the enzymes in the CGMs. Conclusion: Under the conditions of in vitro screening, 14 of the 17 compounds did not exhibit measuarable interference. Galactose exhibited the highest interference during screening, but did not substantially interfere with CGMs under the conditions of in vitro testing at physiologically relevant concentrations. Enzyme kinetic analysis conducted with galactose supported the notion that (1) the reactivity of GOx enzyme toward nonglucose sugars and (2) the presence of enzymatic impurities (such as galactose oxidase) are two potential sources for sugar interference with GOx glucose sensors, and thus, should be considered during device development.

Boehm, Ryan↗

DOE OSTI2015

A New Approach to Reduce Uncertainties in Space Radiation Cancer Risk Predictions

The prediction of space radiation induced cancer risk carries large uncertainties with two of the largest uncertainties being radiation quality and dose-rate effects. In risk models the ratio of the quality factor (QF) to the dose and dose-rate reduction effectiveness factor (DDREF) parameter is used to scale organ doses for cosmic ray proton and high charge and energy (HZE) particles to a hazard rate for γ-rays derived from human epidemiology data. In previous work, particle track structure concepts were used to formulate a space radiation QF function that is dependent on particle charge number Z, and kinetic energy per atomic mass unit, E. QF uncertainties where represented by subjective probability distribution functions (PDF) for the three QF parameters that described its maximum value and shape parameters for Z and E dependences. Here I report on an analysis of a maximum QF parameter and its uncertainty using mouse tumor induction data. Because experimental data for risks at low doses of γ-rays are highly uncertain which impacts estimates of maximum values of relative biological effectiveness (RBE max ), I developed an alternate QF model, denoted QF γAcute where QFs are defined relative to higher acute γ-ray doses (0.5 to 3 Gy). The alternate model reduces the dependence of risk projections on the DDREF, however a DDREF is still needed for risk estimates for high-energy protons and other primary or secondary sparsely ionizing space radiation components. Risk projections (upper confidence levels (CL)) for space missions show a reduction of about 40% (CL~50%) using the QF γAcute model compared the QFs based on RBE max and about 25% (CL~35%) compared to previous estimates. In addition, I discuss how a possible qualitative difference leading to increased tumor lethality for HZE particles compared to low LET radiation and background tumors remains a large uncertainty in risk estimates.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2020

A pseudo-kinetic model to simulate phase changes in gas hydrate bearing sediments

Herein, modeling of the phase transitions anticipated in gas hydrate bearing sediments (GHBS) is critical for a proper understanding of time-dependent changes in states and volumes (e.g. the production of methane from this type of soils). We propose a new pseudo-kinetic approach to simulate the typical phase changes anticipated in GHBS, using published experimental results involving gas hydrate dissociation that are the basis of a widely used kinetic model. The proposed pseudo-kinetic model is formulated in the pressure-temperature (P-T) plane and assumes a rate of gas hydrate dissociation (or formation) proportional to the distance between the current state and the phase boundary. The model consists of only one parameter and is simple to implement in numerical simulators. A similar concept is used to model ice formation/thawing phenomena, but based on the ice/liquid-water phase boundary. We implemented the pseudo-kinetic model in a fully coupled thermo-hydro-chemo-mechanical (THCM) finite element code and validated it against experimental results performed on the dissociation of synthetic gas hydrate. We also evaluated the pseudo-kinetic model using synthetic cases covering several scenarios associated with gas hydrate formation/dissociation and ice formation/thawing. The model successfully reproduced the gas production test from a natural GHBS core from Korea (scaled gas venting experiment over 14 h), and also the formation of gas hydrate and ice in permafrost in Alaska (over 2 × 10 6 years). -The analyses show the versatility of the proposed pseudo-kinetic approach by applying it to model the different types of phase transitions typically encounter in GHBS. The simple formulation, easy implementation in numerical simulator, and reduced number of parameters (only one per phase change) make this model an attractive option for simulating phase transformations in problems involving GHBS.

58 GEOSCIENCES↗

DOE OSTI2014

Reduced Activity of Double-Strand Break Repair Genes in Prostate Cancer Patients With Late Normal Tissue Radiation Toxicity

Purpose: To investigate clinical parameters and DNA damage response as possible risk factors for radiation toxicity in the setting of prostate cancer. Methods and Materials: Clinical parameters of 61 prostate cancer patients, 34 with (overresponding, OR) and 27 without (non-responding, NR) severe late radiation toxicity were assembled. In addition, for a matched subset the DNA damage repair kinetics (γ-H2AX assay) and expression profiles of DNA repair genes were determined in ex vivo irradiated lymphocytes. Results: Examination of clinical data indicated none of the considered clinical parameters to be correlated with the susceptibility of patients to develop late radiation toxicity. Although frequencies of γ-H2AX foci induced immediately after irradiation were similar (P=.32), significantly higher numbers of γ-H2AX foci were found 24 hours after irradiation in OR compared with NR patients (P=.03). Patient-specific γ-H2AX foci decay ratios were significantly higher in NR patients than in OR patients (P<.0001). Consequently, NR patients seem to repair DNA double-strand breaks (DSBs) more efficiently than OR patients. Moreover, gene expression analysis indicated several genes of the homologous recombination pathway to be stronger induced in NR compared with OR patients (P<.05). A similar trend was observed in genes of the nonhomologous end-joining repair pathway (P=.09). This is congruent with more proficient repair of DNA DSBs in patients without late radiation toxicity. Conclusions: Both gene expression profiling and DNA DSB repair kinetics data imply that less-efficient repair of radiation-induced DSBs may contribute to the development of late normal tissue damage. Induction levels of DSB repair genes (eg, RAD51) may potentially be used to assess the risk for late radiation toxicity.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2011

In Situ Hydrothermal Synthesis of LiFePO4 Studied by Synchrotron

The development of high capacity, safe lithium battery materials requires new tools to better understand how reaction conditions affect nucleation and crystallization, particle size, morphology, and defects. We present a general approach for studying the synthesis of Li battery electrode materials in real time. The formation of LiFePO{sub 4} was investigated by time-resolved in situ synchrotron X-ray diffraction under hydrothermal conditions, and the reaction kinetics were determined by changes of the Bragg reflections. We provide the first evidence in support of a dissolution-reprecipitation process for the formation of LiFePO{sub 4}, which occurs at temperatures as low as 105 C and appears to be a three-dimensional diffusion-controlled process. Lattice parameters and their evolution were monitored in situ, as well as the formation of antisite defects and their subsequent elimination under various synthesis conditions. The ability to characterize and tailor synthesis reactions in situ is essential for rapid optimization of the synthesis procedures and, ultimately, the development of new battery electrodes.

25 ENERGY STORAGE↗

DOE OSTI2022

A numerical approach for detecting switch-like bistability in mass action chemical reaction networks with conservation laws

Theoretical analysis of signaling pathways can provide a substantial amount of insight into their function. One particular area of research considers signaling pathways capable of assuming two or more stable states given the same amount of signaling ligand. This phenomenon of bistability can give rise to switch-like behavior, a mechanism that governs cellular decision making. Investigation of whether or not a signaling pathway can confer bistability and switch-like behavior, without knowledge of specific kinetic rate constant values, is a mathematically challenging problem. Recently a technique based on optimization has been introduced, which is capable of finding example parameter values that confer switch-like behavior for a given pathway. Although this approach has made it possible to analyze moderately sized pathways, it is limited to reaction networks that presume a uniterminal structure. It is this limited structure we address by developing a general technique that applies to any mass action reaction network with conservation laws. In this paper we developed a generalized method for detecting switch-like bistable behavior in any mass action reaction network with conservation laws. The method involves (1) construction of a constrained optimization problem using the determinant of the Jacobian of the underlying rate equations, (2) minimization of the objective function to search for conditions resulting in a zero eigenvalue, (3) computation of a confidence level that describes if the global minimum has been found and (4) evaluation of optimization values, using either numerical continuation or directly simulating the ODE system, to verify that a bistability region exists. The generalized method has been tested on three motifs known to be capable of bistability. We have developed a variation of an optimization-based method for the discovery of bistability, which is not limited to uniterminal chemical reaction networks. Successful completion of the method provides an S-shaped bifurcation diagram, which indicates that the network acts as a bistable switch for the given optimization parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2014

Particle in cell simulations of Buneman instability of a current-driven plasma with q-nonextensive electron velocity distribution

The nonlinear evolution of low frequency Buneman instability in an unmagnetized current-driven plasma with q-nonextensive electron velocity distribution is investigated using particle in cell simulation. Simulation results show that the generation of electron phase space holes and the counter-streaming current induced in the plasma strongly depend on the q-parameter. It is found that by increasing the nonextensive parameter, the distribution of electron density becomes highly peaked. This density steepening or grating-like pattern occurs at the saturation time. In addition, a generalized dispersion relation is obtained using the kinetic theory. Analysis of the dispersion relation and the temporal evolution of the electric field energy density reveal that the growth rate of instability increases by increasing the q-parameter. Finally, the results of Maxwellian and q-nonextensive velocity distributions have been compared and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2011

The direct observation of secondary radical chain chemistry in the heterogeneous reaction of chlorine atoms with submicron squalane droplets

The reaction of Cl atoms, in the presence of Cl2 and O2, with sub-micron squalane particles is used as a model system to explore how surface hydrogen abstraction reactions initiate chain reactions that rapidly transform the chemical composition of an organic particle. The heterogeneous reaction is measured in a photochemical flow tube reactor in which chlorine atoms are produced by the photolysis of Cl2 at 365 nm. By monitoring the heterogeneous reaction, using a vacuum ultraviolet photoionization aerosol mass spectrometer, the effective reactive uptake coefficient and the distributions of both oxygenated and chlorinated reaction products are measured and found to depend sensitively upon O2, Cl2, and Cl concentrations in the flow reactor. In the absence of O2, the effective reactive uptake coefficient monotonically increases with Cl2 concentration to a value of ~;;3, clearly indicating the presence of secondary chain chemistry occurring in the condensed phase. The effective uptake coefficient decreases with increasing O2 approaching a diffusion corrected value of 0.65 ? 0.07, when 20percent of the total nitrogen flow rate in the reactor is replaced with O2. Using a kinetic model it is found that the amount of secondary chemistry and the product distributions in the aerosol phase are controlled by the competitive reaction rates of O2 and Cl2 with alkyl radicals. The role that a heterogeneous pathway might play in the reaction of alkyl radicals with O2 and Cl2 is investigated within a reasonable range of reaction parameters. These results show, more generally, that for heterogeneous reactions involving secondary chain chemistry, time and radical concentration are not interchangeable kinetic quantities, but rather the observed reaction rate and product formation chemistry depends sensitively upon the concentrations and time evolution of radical initiators and those species that propagate or terminate free radical chain reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗