Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism allows the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In current research using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For 13 C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Sorption Modeling and Verification for Off-Gas Treatment

The project has made progress toward developing a comprehensive modeling capability for the capture of target species in off gas evolved during the reprocessing of nuclear fuel. The effort has integrated experimentation, model development, and computer code development for adsorption and absorption processes. For adsorption, a modeling library has been initiated to include (a) equilibrium models for uptake of off-gas components by adsorbents, (b) mass transfer models to describe mass transfer to a particle, diffusion through the pores of the particle and adsorption on the active sites of the particle, and (c) interconnection of these models to fixed bed adsorption modeling which includes advection through the bed. For single-component equilibria, a Generalized Statistical Thermodynamic Adsorption (GSTA) code was developed to represent experimental data from a broad range of isotherm types; this is equivalent to a Langmuir isotherm in the two-parameter case, and was demonstrated for Kr on INL-engineered sorbent HZ PAN, water sorption on molecular sieve A sorbent material (MS3A), and Kr and Xe capture on metal-organic framework (MOF) materials. The GSTA isotherm was extended to multicomponent systems through application of a modified spreading pressure surface activity model and generalized predictive adsorbed solution theory; the result is the capability to estimate multicomponent adsorption equilibria from single-component isotherms. This advance, which enhances the capability to simulate systems related to off-gas treatment, has been demonstrated for a range of real-gas systems in the literature and is ready for testing with data currently being collected for multicomponent systems of interest, including iodine and water on MS3A. A diffusion kinetic model for sorbent pellets involving pore and surface diffusion as well as external mass transfer has been established, and a methodology was developed for determining unknown diffusivity parameters from transient uptake data. Two parallel approaches have been explored for integrating the kernels described above into a mass-transport model for adsorption in fixed beds. In one, the GSTA isotherm kernel has been incorporated into the MOOSE framework; in the other approach, a focused finite-difference framework and PDE kernels have been developed. Issues, including oscillatory behavior in MOOSE solutions to advection-diffusion problems, and opportunities have been identified for each approach, and a path forward has been identified toward developing a stronger modeling platform. Experimental systems were established for collection of microscopic kinetics and equilibria data for single and multicomponent uptake of gaseous species on solid sorbents. The systems, which can operate at ambient temperature to 250°C and dew points from -69 to 17°C, are useful for collecting data needed for modeling performance of sorbents of interest. Experiments were conducted to determine applicable models and parameters for isotherms and mass transfer for water and/or iodine adsorption on MS3A. Validation experiments were also conducted for water adsorption on fixed beds of MS3A. For absorption, work involved modeling with supportive experimentation. A dynamic model was developed to simulate CO 2 absorption with chemical reaction using high alkaline content water solutions. A computer code was developed to implement the model based upon transient mass and energy balances. Experiments were conducted in a laboratory-scale column to determine model parameters. The influence of geometric parameters and operating variables on CO 2 absorption was studied over a wide range of conditions. This project has resulted in 7 publications, with 3 manuscripts in preparation. Also, 15 presentations were given at national meetings of ANS and AIChE and at Material Recovery and Waste Forms Campaign Working Group meetings.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Gamma Radiation-Induced Degradation of Acetohydroxamic Acid (AHA) in Aqueous Nitrate and Nitric Acid Solutions Evaluated by Multiscale Modelling

Acetohydroxamic acid (AHA) has been proposed for inclusion in advanced, single-cycle, used nuclear fuel reprocessing solvent systems for the reduction and complexation of plutonium and neptunium ions. For this application, a detailed description of the degradation of AHA in dilute aqueous nitric acid is required. To this end, we present a comprehensive, multiscale computer model for the coupled radiolytic and hydrolytic degradation of AHA in aqueous sodium nitrate and nitric acid solutions. Rate coefficients for the reactions of AHA and HA with the oxidizing nitrate radical were measured for the first time using pulsed electron radiolysis and used as inputs in the kinetic model. The computer model results are validated by comparison to experimental data from steady state gamma irradiations, for which the agreement is excellent. Here, the presented model accurately predicts the yields of the major degradation products of AHA: acetic acid, hydroxylamine, nitrous oxide, and molecular hydrogen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

Image-derived and arterial blood sampled input functions for quantitative PET imaging of the angiotensin II subtype 1 receptor in the kidney

Purpose: The radioligand {sup 11}C-KR31173 has been introduced for positron emission tomography (PET) imaging of the angiotensin II subtype 1 receptor in the kidney in vivo. To study the biokinetics of {sup 11}C-KR31173 with a compartmental model, the input function is needed. Collection and analysis of arterial blood samples are the established approach to obtain the input function but they are not feasible in patients with renal diseases. The goal of this study was to develop a quantitative technique that can provide an accurate image-derived input function (ID-IF) to replace the conventional invasive arterial sampling and test the method in pigs with the goal of translation into human studies. Methods: The experimental animals were injected with [{sup 11}C]KR31173 and scanned up to 90 min with dynamic PET. Arterial blood samples were collected for the artery derived input function (AD-IF) and used as a gold standard for ID-IF. Before PET, magnetic resonance angiography of the kidneys was obtained to provide the anatomical information required for derivation of the recovery coefficients in the abdominal aorta, a requirement for partial volume correction of the ID-IF. Different image reconstruction methods, filtered back projection (FBP) and ordered subset expectation maximization (OS-EM), were investigated for the best trade-off between bias and variance of the ID-IF. The effects of kidney uptakes on the quantitative accuracy of ID-IF were also studied. Biological variables such as red blood cell binding and radioligand metabolism were also taken into consideration. A single blood sample was used for calibration in the later phase of the input function. Results: In the first 2 min after injection, the OS-EM based ID-IF was found to be biased, and the bias was found to be induced by the kidney uptake. No such bias was found with the FBP based image reconstruction method. However, the OS-EM based image reconstruction was found to reduce variance in the subsequent phase of the ID-IF. The combined use of FBP and OS-EM resulted in reduced bias and noise. After performing all the necessary corrections, the areas under the curves (AUCs) of the AD-IF were close to that of the AD-IF (average AUC ratio =1 ± 0.08) during the early phase. When applied in a two-tissue-compartmental kinetic model, the average difference between the estimated model parameters from ID-IF and AD-IF was 10% which was within the error of the estimation method. Conclusions: The bias of radioligand concentration in the aorta from the OS-EM image reconstruction is significantly affected by radioligand uptake in the adjacent kidney and cannot be neglected for quantitative evaluation. With careful calibrations and corrections, the ID-IF derived from quantitative dynamic PET images can be used as the input function of the compartmental model to quantify the renal kinetics of {sup 11}C-KR31173 in experimental animals and the authors intend to evaluate this method in future human studies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Jacobian-free Newton–Krylov method for the simulation of non-thermal plasma discharges with high-order time integration and physics-based preconditioning

A preconditioning framework for the numerical simulation of non-thermal streamer discharges is developed using the Jacobian-free Newton-Krylov (JFNK) method. A reduced plasma fluid model is considered, consisting of electrons, one positive ion, one negative ion, and the electrostatic potential. Here, the plasma kinetics model includes ionization, electron-ion recombination, electron attachment, electron detachment, and ion-ion recombination. The governing equations are made dimensionless, discretized in space with finite differences, and integrated in time with a fully implicit method based on high-order backward differentiation formulas. The preconditioning framework is based on a linearized form of the governing equations and physics-based operator splitting. The efficiency of the preconditioning strategy is assessed through two test cases: streamer propagation between parallel plates and an axisymmetric pin-to-pin discharge. The fully implicit approach overcomes traditional restrictions in the time step size due to processes such as electron drift, electron diffusion, and dielectric relaxation. Excellent performance is observed through relevant statistics of the JFNK solver, although the number of linear iterations increases for the pin-to-pin discharge when nonlinear numerical boundary conditions are imposed at the electrodes. Performance studies show scalability with O(100-1000) processors for O(10M) unknowns with ample room for optimization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2018

Auto-ignition study of FACE gasoline and its surrogates at advanced IC engine conditions

Robust surrogate formulation for gasoline fuels is challenging, especially in mimicking auto-ignition behavior observed under advanced combustion strategies including boosted spark-ignition and advanced compression ignition. This work experimentally quantifies the auto-ignition behavior of bi- and multi-component surrogates formulated to represent a mid-octane (Anti-Knock Index 91.5), full boiling-range, research grade gasoline (Fuels for Advanced Combustion Engines, FACE-F). A twin-piston rapid compression machine is used to achieve temperature and pressure conditions representative of in-cylinder engine operation. Changes in low- and intermediate-temperature behavior, including first-stage and main ignition times, are quantified for the surrogates and compared to the gasoline. This study identifies significant discrepancies in the first-stage ignition behavior, the influence of pressure for the bi- to ternary blends, and highlights that better agreement is achieved with multi-component surrogates, particularly at lower temperature regimes. A recently-updated detailed kinetic model for gasoline surrogates is also used to simulate the measurements. Finally, sensitivity analysis is employed to interpret the kinetic pathways responsible for reactivity trends in each gasoline surrogate.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Effect of Doppler-shifted photons on subnanosecond breakdown in high-voltage pulse discharge

The experiments in high-voltage open discharge in helium [1, 2] showed a controlled current growth rate of 500 A/(cm{sup 2}ns) for an applied voltage of 20 kV and gas pressure of 6 Torr. A kinetic model of the subnanosecond breakdown is developed to analyze the mechanism of current growth, which takes into account the kinetics of electrons, ions, fast atoms and photons with a Doppler shift (DS). DS photons appear in discharge due to collisions of heavy particles. Using particle in cell simulations, we show a critical role of DS photons in the electron emission from the cathode during the breakdown. Our experimental and calculation results show a decrease of the breakdown time with increasing gas pressure from 3 Torr to 16 Torr.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

SUMO-fusion, purification, and characterization of a (+)-zizaene synthase from Chrysopogon zizanioides

An uncharacterized plant cDNA coding for a polypeptide presumably having sesquiterpene synthase activity, was expressed in soluble and active form. Two expression strategies were evaluated in Escherichia coli. The enzyme was fused to a highly soluble SUMO domain, in addition to being produced in an unfused form by a cold-shock expression system. Yields up to ∼325 mg/L{sup −1} were achieved in batch cultivations. The 6x-His-tagged enzyme was purified employing an Ni{sup 2+}-IMAC-based procedure. Identity of the protein was established by Western Blot analysis as well as peptide mass fingerprinting. A molecular mass of 64 kDa and an isoelectric point of pI 4.95 were determined by 2D gel electrophoresis. Cleavage of the fusion domain was possible by digestion with specific SUMO protease. The synthase was active in Mg{sup 2+} containing buffer and catalyzed the production of (+)-zizaene (syn. khusimene), a precursor of khusimol, from farnesyl diphosphate. Product identity was confirmed by GC–MS and comparison of retention indices. Enzyme kinetics were determined by measuring initial reaction rates for the product, using varying substrate concentrations. By assuming a Michaelis–Menten model, kinetic parameters of K{sub M} = 1.111 μM (±0.113), v{sub max} = 0.3245 μM min{sup −1} (±0.0035), k{sub cat} = 2.95 min{sup −1}, as well as a catalytic efficiency k{sub cat}/K{sub M} = 4.43 × 10{sup 4} M{sup −1} s{sup −1} were calculated. Fusion to a SUMO moiety can substantially increase soluble expression levels of certain hard to express terpene synthases in E. coli. The kinetic data determined for the recombinant synthase are comparable to other described plant sesquiterpene synthases and in the typical range of enzymes belonging to the secondary metabolism. This leaves potential for optimizing catalytic parameters through methods like directed evolution. - Highlights: • Uncharacterized (+)-zizaene synthase from C. zizanoides was cloned and expressed. • Fusion to SUMO and cold-shock induction enhanced soluble yields in E. coli. • Ni{sup 2+}-IMAC purification of the SUMO-fused and unfused enzyme. • (+)-Zizaene identified as main cyclization product by GC–MS. • Enzyme kinetic parameters comparable to related sesquiterpene synthases.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Insights into the Adsorption of Carbon Dioxide in the Presence of Water Vapor Utilizing a Low Molecular Weight Polyethylenimine-Impregnated CARiACT Silica Sorbent

Thermogravimetric analysis was employed to investigate the CO 2 and H 2 O adsorption rates and water vapor equilibrium on anhydrous and pre-hydrate linear polyethylenimine (LPEI) sorbent impregnated within a commercially functional CARiACT G10 (HPV) silica support. Water vapor experiments utilizing specific humidity of 2%, 8%, and 16% in contact with an anhydrous PEI sorbent resulted in proportional quantities of water vapor uptake. Subsequently, both anhydrous and pre-hydrated PEI-impregnated sorbents were made available to identical humidified gaseous streams containing a CO 2 concentration of 10% at 60oC. CO 2 capacity increased dramatically in the presence of different levels of humidity. Various kinetic models were systematically employed to interpret the experimental data including single and multiple-step models. The rate data was best represented by a reaction mechanism pathway involving the interplay of CO 2 with PEI-impregnated sorbents exhibited a quick adsorption phase followed by a slow approach to equilibrium. Moreover, a phenomenological rate model was developed to describe the dynamic H 2 O and CO 2 uptakes at specific humidity levels studied. The kinetic study showed good agreement with experimental data. Furthermore, the effects observed during the adsorption and hydration are shown to be complementary to known chemical and physical transformations within the polyethylenimine’s macromolecule.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2020

NUMERICAL INVESTIGATION OF FUEL PROPERTY EFFECTS ON MIXED-MODE COMBUSTION IN A SPARK-IGNITION ENGINE

In the present study, mixed-mode combustion of an E30 fuel in a direct-injection spark-ignition engine is numerically investigated at a fuel-lean operating condition using multidimensional computational fluid dynamics (CFD). A fuel surrogate matching Research Octane Number (RON) and Motor Octane Number (MON) of E30 is first developed using neural network based non-linear regression model. To enable efficient 3D engine simulations, a 164-species skeletal reaction mechanism incorporating NO, chemistry is reducedfrom a detailed chemical kinetic model. A hybrid approach that incorporates the G-equation model for tracking turbulentflame front, and the multi-zone well-stirred reactor model for predicting auto-ignition in the end gas, is employed to account for turbulent combustion interactions in the engine cylinder. Predicted in-cylinder pressure and heat release rate traces agree well with experimental measurements. The proposed modelling approach also captures moderated cyclic variability. Two different types of combustion cycles, corresponding to purely deflagrative and mixed-mode combustion, are observed. In contrast to the purely deflagrative cycles, mixed-mode combustion cycles feature early flame propagation followed by end-gas auto-ignition, leading to two distinctive peaks in heat release rate traces. The positive correlation between mixed-mode combustion cycles and early flame propagation is well captured by simulations. With the validated numerical setup, effects of NOx chemistry on mixed-mode combustion predictions are investigated. NOx chemistry is found to promote auto-ignition through residual gas recirculation, while the deflagrative flame propagation phase remains largely unaffected. Local sensitivity analysis is then performed to understand effects of physical and chemical properties of the fuel, i.e., heat of evaporation (HoV) and laminar flame speed (S-L). An increased HoV tends to suppress end-gas auto-ignition due to increased vaporization cooling, while the impact of HoV on flame propagation is insignificant. In contrast, an increased SL is found to significantly promote both flame propagation and auto-ignition. The promoting effect of SL on auto-ignition is not a direct chemical effect; it is rather caused by an advancement of the combustion phasing, which increases compression heating of the end gas.

detailed chemical kinetics↗
DOE OSTI2021

Pyrolysis of bio-derived dioxolane fuels: A ReaxFF molecular dynamics study

Alkyl-substituted 1,3-dioxolanes, including 4,5-dimethyl-2-pentan-3-yl-1,3-dioxolane (Fuel 1), 4,5-dimethyl-2-pentyl-1,3-dioxolane (Fuel 2), and 2-(heptan-3-yl)-4,5-dimethyl-1,3-dioxolane (Fuel 3), have been recently suggested as potential biodiesels. In this paper, we investigate the initial pyrolysis of the alkyl-substituted 1,3-dioxolanes at high temperatures using ReaxFF molecular dynamics (MD) simulations. We analyze the decomposition rate, reaction mechanism, and product distribution in the pyrolysis of the three alkyl 1,3-dioxlanes. The three fuels primarily decompose to 4,5-dimethyl-1,3-dioxolane radical and hydrocarbons derived from the alkyl side-chains. The further decomposition of 4,5-dimethyl-1,3-dioxolane radical primarily leads to 2-C 4 H 8 and CO 2 within a few decomposition steps. The hydrocarbon product distribution is significantly affected by the molecular structure of the alkyl side-chain, which would have a strong influence on the sooting tendency of these fuels. The ReaxFF simulations predict that the order of sooting tendency would be Fuel 3 > Fuel 1 > Fuel 2, which agrees with the measured sooting tendency trend. Based on the pyrolysis mechanism identified by ReaxFF, we propose a new alkyl dioxolane, 4-hexyl-5-methyl-1,3-dioxolane (Fuel 4), which might produce even less soot, by modifying the molecular structure of Fuel 2. Our ReaxFF simulation shows that Fuel 4 produce much less C 4 H 8 , an effective non-aromatic soot precursor, than Fuel 2. Moreover, more carbon atoms are bonded to each oxygen atom in Fuel 4 than Fuel 2, which would help reduce soot yield by removing more carbon atoms from the soot-producing pool of species. The major decomposition pathways identified in this work can be used to develop chemical kinetic models for 1,3-dioxolane based compounds, as biodiesel components, applicable to combustion engine simulations. We also demonstrate that the chemical kinetic insight offered by ReaxFF simulations can be used to design new fuel molecules with more desired properties.

09 BIOMASS FUELS↗
DOE OSTI2019

Theoretical Study of the Reaction Mechanism and Kinetics of the Phenyl + Allyl and Related Benzyl + Vinyl Associations

Potential energy surfaces for the allyl + phenyl and benzyl + vinyl barrierless radical association reactions have been studied at the CCSD(T)-F12/cc-pVTZ-f12//B3LYP/6-311G** level of theory. Variable reaction coordinate transition state theory (VRC-TST) has been employed to evaluate high-pressure limit rate constants for the barrierless channels. Then, Rice–Ramsperger–Kassel–Marcus master equation (RRKM-ME) calculations have been performed to assess phenomenological rate constants and product branching ratios of various reaction channels at different temperatures and pressures. The initial step of both radical association reactions produces 3-phenylpropene which can further dissociate into a variety of bimolecular products including the indene precursor 1-phenylallyl + H. The results showed that at typical combustion conditions the collisional stabilization of 3-phenylpropene dominates both the phenyl + allyl and benzyl + vinyl reactions at temperatures below 1000 K and remains important at high pressures up to 2500 K. The main bimolecular products of the two reactions at high temperatures are predicted to be benzyl + vinyl and phenyl + allyl, respectively. The well-skipping mechanism to form 1-phenylallyl directly in the allyl + phenyl and benzyl + vinyl reactions appeared to be not significant, however, the reactions can provide some contributions into the formation of the indene precursor via the 3-phenylpropene stabilization/dissociation sequence and most of all, via the formation of 3-phenylpropene itself, which then can undergo H-abstraction by available radicals to produce 1-phenylallyl. The allyl + phenyl reaction can also contribute to the formation of two-ring PAH by producing benzyl radical at high temperatures, either by the well-skipping or stabilization/dissociation mechanisms; in turn, benzyl can readily react with acetylene or propargyl radical to form indene or naphthalene precursors, respectively. Finally, rate expressions for all important reaction channels in a broad range of temperatures and pressures have been generated for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Development of a Thermophysical Properties Model for Flowsheet Simulation of Biomass Pyrolysis Processes

A properties model was developed for use in commercial process simulators to model pyrolysis of lignocellulosic biomass. The component list was chosen to enable process simulations based on a recently published lumped pyrolysis kinetics model. Since many of the compounds involved in pyrolysis are not found in simulator databanks, estimation based on available literature data was used to establish missing parameters. Standard solid enthalpy of formation, solid heat capacity, and solid density estimates calculated from the limited experimental data available were prepared for six biomass constituents and nine intermediate and end-products of their pyrolysis. Ideal gas enthalpy of formation and heat capacity, critical property, and vapor pressure estimates were prepared for another four pyrolysis end-products and one biomass component. The estimates were all validated against the closest available experimental data in the literature. The addition of these new components and properties allows thermodynamically rigorous simulation of lumped biomass pyrolysis reactions with accurate energy balances. Furthermore enthalpies of reaction calculated from the properties model were compared with reported reaction enthalpies for the same lumped biomass pyrolysis reactions and found to be in general agreement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Accuracy of maximum likelihood estimates of a two-state model in single-molecule FRET

Photon sequences from single-molecule Förster resonance energy transfer (FRET) experiments can be analyzed using a maximum likelihood method. Parameters of the underlying kinetic model (FRET efficiencies of the states and transition rates between conformational states) are obtained by maximizing the appropriate likelihood function. In addition, the errors (uncertainties) of the extracted parameters can be obtained from the curvature of the likelihood function at the maximum. We study the standard deviations of the parameters of a two-state model obtained from photon sequences with recorded colors and arrival times. The standard deviations can be obtained analytically in a special case when the FRET efficiencies of the states are 0 and 1 and in the limiting cases of fast and slow conformational dynamics. These results are compared with the results of numerical simulations. The accuracy and, therefore, the ability to predict model parameters depend on how fast the transition rates are compared to the photon count rate. In the limit of slow transitions, the key parameters that determine the accuracy are the number of transitions between the states and the number of independent photon sequences. In the fast transition limit, the accuracy is determined by the small fraction of photons that are correlated with their neighbors. The relative standard deviation of the relaxation rate has a “chevron” shape as a function of the transition rate in the log-log scale. The location of the minimum of this function dramatically depends on how well the FRET efficiencies of the states are separated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Benchmarking and validation of global model code for negative hydrogen ion sources

Benchmarking and validation are prerequisites for using simulation codes as predictive tools. In this work, we have developed a Global Model for Negative Hydrogen Ion Source (GMNHIS) and performed benchmarking of the GMNHIS against another independently developed code, Global Enhanced Vibrational Kinetic Model (GEVKM). This is the first study to propose a quite comprehensive benchmarking test of this kind for models of negative hydrogen ion sources (NHIS), and excellent agreements have been achieved for collisional energy loss per electron-ion pair created, electron number density, electron temperature, densities of H$^{+}_{3}$ and H$^{+}_{2}$ ions, and densities of H(n = 1–3) atoms. Very small discrepancies in number densities of H - ions and H + ions, as well as the vibrational distribution function of hydrogen molecules, can be attributed to the differences in the chemical reactions datasets. The GEVKM includes additional chemical reactions that are more important at high pressures. Moreover, we validated the GMNHIS against experimental data obtained in an electron cyclotron resonance discharge used for H - production. The model qualitatively (and even quantitatively for certain conditions) reproduces the experimental H - number density. The H - number density as a function of pressure first increases at pressures below 1.6 Pa and then saturates for higher pressures. This dependence was analyzed by evaluating contributions from different reaction pathways to the creation and loss of the H - ions. The developed codes can be used for predicting the H - production, improving the performance of NHIS, and ultimately optimizing the parameters of negative ion beams for fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Kinetic Simulations of Cu Doping in Chlorinated CdSeTe PV Absorbers

This paper reports on 1D kinetic modeling of p-type doping formation in Cu-doped chlorinated CdSeTe photovoltaic (PV) absorbers with graded Se profiles. Following the recent progress in kinetic simulations of the defect chemistry in stoichiometric CdTe films, this work extends simulation capabilities to a domain of semiconductor alloy films with graded stoichiometries. The defect formation energies and ionization levels, as well as the diffusion and reaction barriers used in the reaction-diffusion equations are calculated from the first principles. A new formalism is employed to account for the uncertainty of the defect formation energies in semiconductor alloys caused by unknown local surrounding. Finally, the developed methodology is used to model a practical fabrication process of p-type doping formation in graded CdSeTe absorbers and to study the effect of alloy stoichiometry on the doping activation.

14 SOLAR ENERGY↗
DOE OSTI2014

Final Report: Development of a Chemical Model to Predict the Interactions between Supercritical CO2, Fluid and Rock in EGS Reservoirs

This report summarizes development of a coupled-process reservoir model for simulating enhanced geothermal systems (EGS) that utilize supercritical carbon dioxide as a working fluid. Specifically, the project team developed an advanced chemical kinetic model for evaluating important processes in EGS reservoirs, such as mineral precipitation and dissolution at elevated temperature and pressure, and for evaluating potential impacts on EGS surface facilities by related chemical processes. We assembled a new database for better-calibrated simulation of water/brine/ rock/CO2 interactions in EGS reservoirs. This database utilizes existing kinetic and other chemical data, and we updated those data to reflect corrections for elevated temperature and pressure conditions of EGS reservoirs.

15 GEOTHERMAL ENERGY↗