Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic 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.

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

At least 73 records · Page 4

TRIPOLI-4 neutronics calculations for IAEA-CRP benchmark of CEFR start-up tests using new libraries JEFF-3.3 and ENDF/B-VIII

Sodium-cooled Fast Reactors (SFR) are one of the promising Generation IV fast reactors. The China Experimental Fast Reactor (CEFR) is a 65 MWth pool-type SFR with a high neutron leakage core using high enriched uranium oxide fuel. The CEFR start-up tests in 2010 consist of series of neutronics experiments. Essential experimental data are available from the 2018 IAEA-CEFR-CRP benchmark document and useful for the validation of neutron transport codes and nuclear data libraries. The TRIPOLI-4 Monte Carlo transport code is a general-purpose neutronics code using continuous-energy nuclear data libraries. It has a rich validation database covering different computational and experimental benchmark data sets to assure the accuracy and credibility of numerical studies. Previous TRIPOLI-4 SFR core physics calculations used mainly MOX fuels. SFR control rod worth studies with TRIPOLI-4 were mainly on big-size cores. To test the modeling capability of TRIPOLI-4 on fuel loading patterns and on variable control rods positions for SFRs and to check recent nuclear data libraries, it is interesting to investigate the CEFR start-up tests in this work, including core states from subcritical to supercritical, control rods and rod group worth, point kinetics parameters, and radial foil activation measurements. TRIPOLI-4 calculation results using new data libraries JEFF-3.3, ENDF/B-VIII, and those of ENDF/B-VII.1 for different CEFR core states are reported here. By means of different options of TRIPOLI-4, control-rod SAs reactivity worth, point kinetics parameters, and radial {sup 237}Np(n,f) fission rate distribution were successfully evaluated. Both ENDF/B-VIII and JEFF-3.3 nuclear data libraries provided reliable results with a difference of 170 +/- 11 pcm in k{sub eff}.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Verification of the PERSENT Software

Ongoing commercial design activities require a thorough verification of the Argonne Reactor Computation codes be performed. DIF3D is central to this system and substantial work has been done to verify its accuracy on several identified commercial needs. This manuscript details the verification work done on PERSENT which relies upon the DIF3D code for its forward and adjoint flux solution. Previous work identified the PERSENT features required to be verified to support commercial design activities, features of which are generally applicable to hexagonal-Z fast reactor designs. The scope of this verification effort includes verifying PERSENT’s ability to correctly calculate four key quantities: perturbation worth distributions, kinetics parameters, sensitivity coefficients, and cross section uncertainty quantification. This manuscript provides the verification tasks and their results with respect to these quantities needed for commercial design activities. For the perturbation worth distributions, hand calculations are deployed to verify the PERSENT calculated results. Similarly, hand calculation of the PERSENT computed kinetics parameters is also used to verify the PERSENT results. In both of these, the input to PERSENT is manipulated to ensure the hand calculation exactly matches the equations PERSENT is calculating. The sensitivity coefficients involve calculating the derivatives of a parameter (such as reactivity worth), with respect to the cross section data. Direct finite difference calculations with DIF3D are used to verify the PERSENT calculated results. For the uncertainty quantification, manufactured input to PERSENT is used to allow an exact hand calculation to reproduce the PERSENT calculated results. The work detailed in this report verified that significant issues were identified for earlier versions of PERSENT for sensitivity coefficients which were corrected in this work and thus version 12.1.0 of PERSENT must be used to reproduce all of the verified work in this report.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Morphology selection in dealloying: A phase field study of the coupling among kinetic mechanisms

A phase field model is used to investigate how the corrosion rate and morphology of a dealloying metal is controlled by the kinetic rates of bulk solid diffusion, interfacial diffusivity, and interface dissolution. A specific focus is investigating the manner by which variations in these kinetic parameters synergistically influence the onset of a morphological instability at the solid-liquid interface, which is known to precede rapid dealloying corrosion towards highly topologically complex structures such as bicontinuous porous solids. This work demonstrates that this instability can be suppressed by either an increase of interdiffusion in the solid alloy, or by a decreasing rate of dissolution at the solid-liquid interface. When the instability is suppressed, the result is a relatively planar, passivated dealloying front. Importantly, however, coupled changes in kinetic parameters can also promote the morphological instability and sustain fine-scale ligament and pore formation. Finally, the relevance of the solid-liquid interfacial diffusion for determining the rate of dealloying and ligament formation is highlighted. By elucidating the mechanistic influence of bulk and interfacial kinetic factors on the dealloyed morphology, this work helps to rationalize the contrasting morphologies observed in various dealloying systems.

Diffusion↗

Kinetic and X-ray Absorption Spectroscopic Analysis of Catalytic Redox Cycles over Highly Uniform Polymetal Oxo Clusters

Metal–organic framework materials (MOFs) offer an opportunity for investigating catalytic properties of polymetal oxo clusters that are highly well defined and uniform in nature, in contrast to other classes of catalysts that may exhibit a propensity toward active site heterogeneity. We report herein a kinetic and X-ray absorption spectroscopy (XAS) analysis of the two-electron oxidation of CO over divalent metal sites in MIL-100(M = Fe, Cr) (MIL = Materials of Institut Lavoisier) materials carrying μ 3 -oxo bridged trimers, and connect observations about the kinetic relevance of redox steps to density functional theory (DFT) predictions published previously. The high degree of uniformity evident from in situ titration measurements leads to a congruence in mechanistic inferences made from steady-state catalytic, transient stoichiometric, isotopic exchange, and isotopic tracer data that all point to a sequential mechanism comprised of separate oxidation and reduction half-cycle steps conjoined by an active oxygen intermediate. In situ XAS data reinforce mechanistic conclusions derived from kinetic analysis, and suggest that the active oxygen intermediate may be more appropriately characterized as an iron-oxyl (Fe 3+ –O – ) rather than an iron-oxo (Fe 4+ =O 2– ) species. The Cr analogue of MIL-100 exhibits contrasting rate features that can be rationalized using an identical sequence of steps as MIL-100(Fe), but with a highly dissimilar set of kinetic parameters that can also be validated using transient stoichiometric experiments. The larger coverages of active oxygen intermediates on MIL-100(Cr) are consistent with predictions from prior DFT studies that suggest more stable and less reactive active oxygen species for metals with lower d-electron counts, and point to metal identity as a lever for precisely controlling the kinetic relevance of oxidation and reduction half-cycles in catalytic redox sequences over polymetal oxo clusters. Importantly, the results presented point to the utility of developing broadly applicable structure–catalytic property relationships over MOF nodes specifically, and highly uniform catalysts more generally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Alternative bases in the RNA world: the prebiotic synthesis of urazole and its ribosides

Urazole is a five-membered heterocyclic compound which is isosteric with uracil's hydrogen-bonding segment. Urazole reacts spontaneoulsy with ribose (and other aldoses) to give a mixture of four ribosides: alpha and beta pyranosides and furanosides. This reaction occurs in aqueous solution at mild temperatures. Thermodynamic and kinetic parameters for the reaction of urazole with ribose were determined. In contrast, uracil is completely unreactive with ribose under these conditions. Urazole's unusual reactivity is ascribed to the hydrazine portion of the molecule. Urazole can be synthesized from biuret and hydrazine under prebiotic conditions. The prebiotic synthesis of guanazole, which is isosteric in part to diaminopyrimidine and cytosine, is accomplished from dicyandiamide and hydrazine. Kinetic parameters for both prebiotic reactions were measured. Urazole and guanazole are transparent in the UV, which would be a favorable property in the absence of an ozone layer on the early Earth. Urazole makes hydrogen bonds with adenine in DMSO similar to those of uracil, as established by H NMR. All of these properties make urazole an attractive potential precursor to uracil and guanazole a potential precursor to cytosine in the RNA or pre-RNA world.

Non-NASA Center↗

Reaction Mechanism Generator v3.0: Advances in Automatic Mechanism Generation

In chemical kinetics research, kinetic models containing hundreds of species and tens of thousands of elementary reactions are commonly used to understand and predict the behavior of reactive chemical systems. Reaction Mechanism Generator (RMG) is a software suite developed to automatically generate such models by incorporating and extrapolating from a database of known thermochemical and kinetic parameters. Here, we present the recent version 3 release of RMG and highlight improvements since the previously published description of RMG v1.0. Most notably, RMG can now generate heterogeneous catalysis models in addition to the previously available gas- and liquid-phase capabilities. For model analysis, new methods for local and global uncertainty analysis have been implemented to supplement first-order sensitivity analysis. The RMG database of thermochemical and kinetic parameters has been significantly expanded to cover more types of chemistry. The present release includes parallelization for faster model generation and a new molecule isomorphism approach to improve computational performance. RMG has also been updated to use Python 3, ensuring compatibility with the latest cheminformatics and machine learning packages. Overall, RMG v3.0 includes many changes which improve the accuracy of the generated chemical mechanisms and allow for exploration of a wider range of chemical systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Rate Equations for Reversible Disproportionation Reactions and Fitting to Time-Course Data

Integrated rate equations are straightforward to fit to experimental data to verify a proposed mechanism and to extract kinetic parameters. Such equations are derived for reversible disproportionation/comproportionation reactions with any set of initial concentrations. Extraction of forward and reverse rate constants from experimental data by fitting the rate law to the data is demonstrated for the disproportionation of 2,2,6,6-tetramethyl-1-piperidinyl-N-oxyl (TEMPO) under acidic conditions where the approach to equilibrium is observed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Combined Experimental and Modeling Study of the Interactions of Acid Gas with Common Spacecraft Surfaces for Fire Safety Applications

A fire in a spacecraft poses detrimental consequences and risks mission success in addition to crew safety. This is compounded during long-duration missions when the crew has limited options to recover from a fire. A common spacecraft fire concern is the smoldering of wire insulation, typically made from Polyvinyl chloride (PVC) or Polytetrafluoroethylene (PTFE). This creates acid gases such as Hydrogen Chloride (HCl), Hydrogen Fluoride (HF) and Hydrogen Cyanide (HCN). These poisonous gases are hazardous to the crew. They also interact with common surfaces within the spacecraft more than dominant combustion products such as CO2 and H2O. This makes them more difficult to track for potential fire detection techniques, or for postfire clean-up. It is imperative to be able to understand and predict the fate of these poisonous species in a microgravity environment in order to design a safe vehicle. HCl interacts with a number of materials inside a spacecraft. Primary among these materials is aluminum, which is abundantly used due to its strong and light weight nature. Aluminum has a natural oxide layer that protects it from corrosion but is typically treated to enhance this oxide layer. Among these treatments is a chromate conversion coating (CCC), which provides a thin enough protective oxide layer to still conduct electricity, and a traditional anodized material that has a thicker oxide layer that does not conduct electricity. Nomex is another common material found inside a spacecraft. It is a flame-resistant woven polymer that is related to nylon. This commercially available material is used for cargo storage bags and as a fire barrier. Physics-based models were developed to predict the uptake of HCl by these materials. The ultimate objective of these models is to predict the fate of HCl within the spacecraft so that sensors can be placed in meaningful locations in future missions based on the model predictions. To support these modeling efforts, experiments were performed in a cast acrylic test cell that measured the difference between the inlet and outlet concentration of HCl after inserting a sample rod of the test material. Different uptake capacities were realized for each type of sample tested. A computational fluid dynamics model (CFD) model of the reactor was then constructed that used a one-step global reaction rate with calibratable reaction (or kinetic) constants. These constants were calibrated to match the HCl uptake on the CCC aluminum samples, and the same kinetic constants were then tested for the stock and anodized aluminum samples. Model predictions matched the experimental data for the stock aluminum, and to a much lesser extent, the anodized aluminum. The model was additionally validated at different flow rates, sample surface areas, and inlet concentrations, and showed good agreement for all stock and CCC samples. The model did not accurately predict the HCl uptake in the anodized samples compared to the other two types of aluminum. Adjusting the kinetic constants and transport properties did little to improve the prediction. X-Ray Photoelectron Spectroscopy (XPS) was used to determine that the oxide layer thickness of anodized aluminum is approximately 5,000 nm, compared to 250 nm for CCC and 50 nm for stock. XPS also revealed presence of chlorine further down in the aluminum oxide layer in anodized samples than CCC and stock samples after the samples were saturated with HCl, indicating that accounting for diffusion of HCl into the oxide layer is important for accurate prediction of HCl uptake onto anodized aluminum. Consequently, a multi-scale model was developed and tested. First, a single pore inside the anodized aluminum oxide layer was modeled and is referred to as the pore-scale model. In this model, HCl diffused through the pore and reacted with the aluminum oxide pore wall to create aluminum chloride. The sample was then saturated when the mass transfer resistance through the growing aluminum chloride layer became too large for the HCl to reach the aluminum oxide wall and continue the reaction. This pore-scale model was coupled to the reactor-scale model using a concentration-dependent diffusion coefficient, resulting in much more accurate predictions (approximately half the sum square error of the aforementioned reactor-scale model that produced good agreement for stock and CCC) for a variety of operating conditions. The amount of water vapor or relative humidity (RH) in the flow during a reactor experiment was determined to influence HCl uptake. Experiments were performed to understand the interaction of gaseous HCl with aluminum surfaces in the presence of water vapor. The results show that increasing levels of RH increased the capacity of aluminum to adsorb HCl but decreased the capacity of Nomex to uptake HCl. A series of tests were performed on individual aluminum samples after they had been saturated with a fixed concentration of HCl in dry air conditions with the goal of determining how their HCl uptake capacity changes after various treatments with water relative to the original saturation tests. HCl-saturated aluminum samples subjected to a second dry air flow at the same HCl concentration as the original test had an uptake of 23.5% of the original sample with no treatment in between. Saturated aluminum samples subjected to an air flow with a RH of 90% in between tests had an uptake of 35.6% of the original. Saturated aluminum samples submerged in distilled water for 12 hours in between tests had an uptake of 82.2% of the original sample. Previously saturated aluminum tested with HCl and a 50% RH air flow resulted in similar uptake characteristics in multiple repeated tests. The experiments show the profound effect water vapor has on HCl uptake onto aluminum surfaces. In the samples subjected to water vapor or liquid water, capillary condensation and capillary diffusion alters the transport of HCl significantly. A model was proposed that developed a relationship between RH and the coefficient of HCl diffusion in aluminum chloride. This produced an “S-shaped” curve with diffusion coefficient as a function of RH, with 45% RH represented as the point where the diffusion coefficient is halfway between no water saturation and 100% water saturation in the aluminum chloride product layer. No difference in uptake characteristics for the experiment or model were realized between 50% and 62% RH. The results from the large-scale microgravity experiment, Saffire, are discussed as they pertain to the fate of HCl throughout a spacecraft. HCl was released, both as a standalone event, and in concurrence with the burning of a structured cloth. These events only produced a small response in the far field HCl sensor, while a PMMA burn that did not produce HCl had a significantly greater response. A ground-based large-scale facility was constructed to flow acid gas at the scale and configuration realized in the Saffire experiments. A CFD model of this duct was constructed to test kinetic parameters developed in this work at a larger scale and different geometric configuration and to predict the results of the large-scale facility. The models developed in this work were used to interpret the results of the microgravity tests and lead the discussion on what further experiments and models are needed in order to predict the fate of acid gas in a spacecraft environment. To summarize, the major contributions of this work are as follows: the capacity to uptake HCl, with and without the presence of water vapor, was measured for a variety of real spacecraft surfaces. Several different models (single reactor-scale, multiscale, spacecraft-scale) were developed and with the aid of modeling, the rate of uptake for those surfaces was also predicted and validated. The kinetic parameters determined from the small-scale reactor experiments and models were used to predict large-scale and microgravity tests. Conclusions from this research will be used in the design of spacecraft vehicles and large-scale microgravity fire safety experiments. The models built by this work will aid designers in sensor placement and could be used to predict acid gas transport from fires in partial gravity, as would be seen in Lunar and Martian habitats.

fire safety↗

Electrochemical Study on the Kinetic Properties of Fe 2+ /Fe, Ni 2+ /Ni, Cr 2+ /Cr and Cr 3+ /Cr 2+ in Molten MgCl 2 -KCl-NaCl Salts

Here, the fundamental kinetic parameters of corrosion product ions (Fe, Ni and Cr ions) are of significant importance to understand the corrosion mechanisms of the structural materials in the molten salts. In the present study, the fundamental data of Fe 2+ , Ni 2+ , Cr 2+ and Cr 3+ (diffusion coefficient D, exchange current density i 0 , charge transfer coefficient α, limiting current density i L and standard rate constant k0) were measured at different concentrations (1.53 × 10 –6 –7.48 × 10 –4 mol cm –3 ) and temperatures (600 °C–800 °C). The values of D are independent of concentrations and follow Arrhenius law with temperature, which descend in the order of D Cr 2+ > D Fe 2+ > D Cr 3+ > D Ni 2+ with values ranging from 0.94 × 10 –5 to 3.31 × 10 –5 cm 2 s –1 . Both i 0 and k 0 depend on the temperatures, which also follow Arrhenius law.

14 SOLAR ENERGY↗

Assessment of Uncertainty in the Determination of Kinetic Reaction Parameters for Polymeric Materials

The use of thermogravimetric analysis to obtain data describing the thermal response of a polymeric resin with temperature. This data can then be used to obtain activation energy and pre-exponential factor used in an Arrhenius representation of material ablation. The methods which allow to assess the uncertainties associated with the experimental determination of activation energy and pre-exponential factor are employed.

Darby, Stephanie P.↗

Epidermal Homeostasis and Radiation Responses in a Multiscale Tissue Modeling Framework

The surface of skin is lined with several thin layers of epithelial cells that are maintained throughout life time by a small population of stem cells. High dose radiation exposures could injure and deplete the underlying proliferative cells and induce cutaneous radiation syndrome. In this work we propose a multiscale computational model for skin epidermal dynamics that links phenomena occurring at the subcellular, cellular, and tissue levels of organization, to simulate the experimental data of the radiation response of swine epidermis, which is closely similar to human epidermis. Incorporating experimentally measured histological and cell kinetic parameters, we obtain results of population kinetics and proliferation indexes comparable to observations in unirradiated and acutely irradiated swine experiments. At the sub-cellular level, several recently published Wnt signaling controlled cell-cycle models are applied and the roles of key components and parameters are analyzed. Based on our simulation results, we demonstrate that a moderate increase of proliferation rate for the survival proliferative cells is sufficient to fully repopulate the area denuded by high dose radiation, as long as the integrity of underlying basement membrane is maintained. Our work highlights the importance of considering proliferation kinetics as well as the spatial organization of tissues when conducting in vivo investigations of radiation responses. This integrated model allow us to test the validity of several basic biological rules at the cellular level and sub-cellular mechanisms by qualitatively comparing simulation results with published research, and enhance our understanding of the pathophysiological effects of ionizing radiation on skin.

Hu, Shaowen↗

Unexpected Thermomechanical Behavior of Off-Stoichiometry Epoxy/Amine Materials

Recent studies on off-stoichiometric thermosets reveal unique viscoelastic behavior derived from increased free volume and physical interactions between chain ends. To understand structural characteristics arising from cure and its effect on properties, we developed a Monte Carlo model based on step-growth polymerization. Our model accurately predicted structure–property trends for a two-component system of EPON 828 (EPON) and ethylenediamine. A second epoxy monomer, D.E.R. 732 (DER), was investigated to modulate T g . Binary mixtures of EPON and DER in off-stoichiometric, amine-rich formulations resulted in nonlinear evolution of thermomechanical properties with respect to initial formulation stoichiometry. Modifying our model with kinetic parameters allowing for differential epoxide/amine reaction kinetics only partially accounted for trends in T g , suggesting that spatiotemporal contributions-not captured by our model-were significant determinants of material properties compared to polymer architecture for three-component systems. Here these findings underpin the importance of spatial awareness in modeling to inform the development of dynamic thermosets.

36 MATERIALS SCIENCE↗

Microkinetic Modeling: A Tool for Rational Catalyst Design

The design of heterogeneous catalysts relies on understanding the fundamental surface kinetics that controls catalyst performance, and microkinetic modeling is a tool that can help the researcher in streamlining the process of catalyst design. Microkinetic modeling is used to identify critical reaction intermediates and rate-determining elementary reactions, thereby providing vital information for designing an improved catalyst. In this review, we summarize general procedures for developing microkinetic models using reaction kinetics parameters obtained from experimental data, theoretical correlations, and quantum chemical calculations. We examine the methods required to ensure the thermodynamic consistency of the microkinetic model. We describe procedures required for parameter adjustments to account for the heterogeneity of the catalyst and the inherent errors in parameter estimation. We discuss the analysis of microkinetic models to determine the rate determining reactions using the degree of rate control and reversibility of each elementary reaction. We introduce incorporation of Brønsted–Evans–Polanyi relations and scaling relations in microkinetic models and the effects of these relations on catalytic performance and formation of volcano curves. We review the analysis of reaction schemes in terms of the maximum rate of elementary reactions, and we outline a procedure to identify kinetically significant transition states and adsorbed intermediates. Here, we explore the application of generalized rate expressions for the prediction of optimal binding energies of important surface intermediates and to estimate the extent of potential rate improvement. We also explore the application of microkinetic modeling in homogeneous catalysis, electro-catalysis, and transient reaction kinetics. We conclude by highlighting the challenges and opportunities in the application of microkinetic modeling for catalyst design.

42 ENGINEERING↗

Towards Self-Replicating Chemical Systems Based on Cytidylic and Guanylic Acids

This project is aimed towards a better understanding of template-directed reactions and, based on this, towards the development of efficient non-enzymatic RNA replicating systems. These systems could serve as models for the prebiotic synthesis of an RNA world. The major objectives of this project were: (a) To elucidate the mechanistic aspects of template-directed (TD) chemistry, (b) to identify the conditions, environmental and other, that favor "organized chemistry" and stereo selective polymerization of nucleotides and (c) to search and, hopefully, find catalysts that will improve the efficiency of these reactions. Enhanced efficiency is expected to facilitate the road towards a self-replicating chemical system based on all four nucleic acid bases. During the first nine months of the granting period from January 1997 to October 1997, we have made substantial progress towards the first two objectives. During this period our activities were directed towards (1) synthesizing activated nucleotides to be used as substrates, (2) using these substrates in order to determine the effect of the leaving group (imidazole (Im), 2-methylimidazole (2-MeIm), and 2,4-dimethylimidazole (2,4-diMeIm)) in the product distribution, (3) developing techniques for analysis of mixtures by LC/MS, (4) creating a protocol in order to obtain kinetic parameters of the dimerization reaction and (5) analyzing kinetic data obtained with the poly(C)/2-MeImpG system. With the exception of item (5), the experimental work for the projects (1) - (4) is still in progress. A list of publications and manuscripts resulted from this research is enclosed.

Kanavarioti, Anastassia↗

Comparative study of two Saccharomyces cerevisiae strains with kinetic models at genome-scale

The parameterization of kinetic models requires measurement of fluxes and/or metabolite levels for a base strain and a few genetic perturbations thereof. Unlike stoichiometric models that are mostly invariant to the specific strain, it remains unclear whether kinetic models constructed for different strains of the same species have similar or significantly different kinetic parameters. This important question underpins the applicability range and prediction limits of kinetic reconstructions. To this end, herein we parameterize two separate large-scale kinetic models using K-FIT with genome-wide coverage corresponding to two distinct strains of Saccharomyces cerevisiae: CEN.PK 113-7D strain (model k-sacce306-CENPK), and growth-deficient BY4741 (isogenic to S288c; model k-sacce306-BY4741). The metabolic network for each model contains 306 reactions, 230 metabolites, and 119 substrate-level regulatory interactions. The two models (for CEN.PK and BY4741) recapitulate, within one standard deviation, 77% and 75% of the fitted dataset fluxes, respectively, determined by 13 C metabolic flux analysis for wild-type and eight single-gene knockout mutants of each strain. Strain-specific kinetic parameterization results indicate that key enzymes in the TCA cycle, glycolysis, and arginine and proline metabolism drive the metabolic differences between these two strains of S. cerevisiae. Our results suggest that although kinetic models cannot be readily used across strains as stoichiometric models, they can capture species-specific information through the kinetic parameterization process.

59 BASIC BIOLOGICAL SCIENCES↗

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Metastable Solution Thermodynamic Properties and Crystal Growth Kinetics

The crystal growth rates of NH4H2PO4, KH2PO4, (NH4)2SO4, KAl(SO4)2 central dot 12H2O, NaCl, and glycine and the nucleation rates of KBr, KCl, NaBr central dot 2H2O, (NH4)2Cl, and (NH4)2SO4 were expressed in terms of the fundamental driving force of crystallization calculated from the activity of supersaturated solutions. The kinetic parameters were compared with those from the commonly used kinetic expression based on the concentration difference. From the viewpoint of thermodynamics, rate expressions based on the chemical potential difference provide accurate kinetic representation over a broad range of supersaturation. The rates estimated using the expression based on the concentration difference coincide with the true rates of crystallization only in the concentration range of low supersaturation and deviate from the true kinetics as the supersaturation increases.

Kim, Soojin↗

Development of a Reference Dose for Perchlorate: Current Issues and Status

The perchlorate anion (ClO4) is typically manufactured as the ammonium salt. The most common use of ammonium perchlorate is in the aerospace program as a component of solid rocket fuel. The perchlorate anion is exceedingly stable under environmental conditions and has been found in ground and surface waters in CA, NV, UT, AZ, TX, AK, NY, MD, WV and FL. The National Center for Environmental Assessment (NCEA) of the U.S. Environmental Protection Agency (US EPA) is in the process of developing an oral reference dose (RfD) for perchlorate. An oral RfD is a body-weight-adjusted dose that can be consumed daily over an entire lifetime with the expectation of no adverse health effects. Once developed, the new RfD will be used by US EPA as the basis of a safe-drinking-water level (SDWL) guideline. US EPA and regional regulatory agencies will then jointly or separately propose clean-up action levels for ground and surface waters at contaminated sites. The toxicological database on CIO4- as of March 1997 was determined by an expert peer-review panel to be inadequate for the purpose of deriving an oral RfD. For example, little or no experimental data existed on the subchronic, reproductive, or developmental toxicity of perchlorate. To fill gaps in the toxicological database, eight animal studies were designed by a government-industry consortium that included US EPA and AFRL. These studies were performed in 1997-1998. It has been known for many years that in the thyroid, high doses of perchlorate block the function of iodide by competing for iodide binding sites. Perchlorate was used in the 1950s-60s as a treatment for Graves' disease (a hyperthyroid condition). Because of what was already known about the pharmacological mode of action of perchlorate, specific concerns addressed in the design of the recent animal studies included the potential for developmental toxicity, notably neurological development. Upon review of complete study reports from four of the studies and incomplete data from the other two, US EPA/NCEA issued a document in December 1998 stating that the critical effects were hormone and histology data from the neurodevelopmental/behavioral toxicity study that had been performed in rats; an RfD was proposed based upon the observation of thyroid-follicular-cell hypertrophy in the 5-day-old pups (PND) of dams given perchlorate in drinking water at 0.1 mg/kg-day (the lowest dose tested) in this study. US EPA/NCEA also focused attention on a nonsignificant increase in motor activity observed in 14-day-old male pups in the same study at the same dose. In February 1999, a public workshop was convened at which an external peer-review panel determined that there was still insufficient data to support the development of a RfD. Additional studies were recommended. Since that time a clinical study revealed significant depression of thyroidal uptake of iodide in volunteers receiving potassium perchlorate at 10 mg/day for 14 days. Follow-up studies currently underway include a pathology working-group review of histology slides from the all of the previous studies, a motor activity study similar to the first behavioral study, an additional immunotoxicity study, male and female rat kinetic studies, a hormone interlaboratory study, and an effects study for thyroid and brain in developing rats. These studies are expected to be complete by mid-2000. A clinical study is planned to determine inhibition of iodide uptake and kinetic parameters in humans following 14-day perchlorate ingestion. The kinetics studies will provide data for developing physiologically based pharmacokinetic (PBPK) models for rats and humans needed to establish the RfD.

Pleus, R. C.↗