Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactivity Coefficients”

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 163 records · Page 9

On the evaluation of antenna quality factors

Antenna reactive energies and modal quality factors by compact expressions given in terms of polynomials with positive coefficients not involving spherical Bessel functions

Kalafus, R. M.↗

Developing Multilayer Thin Film Strain Sensors With High Thermal Stability

A multilayer thin film strain sensor for large temperature range use is under development using a reactively-sputtered process. The sensor is capable of being fabricated in fine line widths utilizing the sacrificial-layer lift-off process that is used for micro-fabricated noble-metal sensors. Tantalum nitride films were optimized using reactive sputtering with an unbalanced magnetron source. A first approximation model of multilayer resistance and temperature coefficient of resistance was used to set the film thicknesses in the multilayer film sensor. Two multifunctional sensors were fabricated using multilayered films of tantalum nitride and palladium chromium, and tested for low temperature resistivity, TCR and strain response. The low temperature coefficient of resistance of the films will result in improved stability in thin film sensors for low to high temperature use.

Wrbanek, John D.↗

Testing of an Integrated Reactor Core Simulator and Power Conversion System with Simulated Reactivity Feedback

A Direct Drive Gas-Cooled (DDG) reactor core simulator has been coupled to a Brayton Power Conversion Unit (BPCU) for integrated system testing at NASA Glenn Research Center (GRC) in Cleveland, Ohio. This is a closed-cycle system that incorporates an electrically heated reactor core module, turboalternator, recuperator, and gas cooler. Nuclear fuel elements in the gas-cooled reactor design are replaced with electric resistance heaters to simulate the heat from nuclear fuel in the corresponding fast spectrum nuclear reactor. The thermodynamic transient behavior of the integrated system was the focus of this test series. In order to better mimic the integrated response of the nuclear-fueled system, a simulated reactivity feedback control loop was implemented. Core power was controlled by a point kinetics model in which the reactivity feedback was based on core temperature measurements; the neutron generation time and the temperature feedback coefficient are provided as model inputs. These dynamic system response tests demonstrate the overall capability of a non-nuclear test facility in assessing system integration issues and characterizing integrated system response times and response characteristics.

Bragg-Sitton, Shannon M.↗

Testing of an Integrated Reactor Core Simulator and Power Conversion System with Simulated Reactivity Feedback

A Direct Drive Gas-Cooled (DDG) reactor core simulator has been coupled to a Brayton Power Conversion Unit (BPCU) for integrated system testing at NASA Glenn Research Center (GRC) in Cleveland, OH. This is a closed-cycle system that incorporates an electrically heated reactor core module, turbo alternator, recuperator, and gas cooler. Nuclear fuel elements in the gas-cooled reactor design are replaced with electric resistance heaters to simulate the heat from nuclear fuel in the corresponding fast spectrum nuclear reactor. The thermodynamic transient behavior of the integrated system was the focus of this test series. In order to better mimic the integrated response of the nuclear-fueled system, a simulated reactivity feedback control loop was implemented. Core power was controlled by a point kinetics model in which the reactivity feedback was based on core temperature measurements; the neutron generation time and the temperature feedback coefficient are provided as model inputs. These dynamic system response tests demonstrate the overall capability of a non-nuclear test facility in assessing system integration issues and characterizing integrated system response times and response characteristics.

Bragg-Sitton, Shannon M.↗

Far-Field Pattern Measurement and Simulation of VHF Antenna at 60 MHz for Europa Clipper Mission

This paper presents measurements and simulations of a linearly polarized VHF folded dipole operating at 60 MHz for NASA’s upcoming Europa Clipper mission. The spacecraft contains a ground penetrating radar, consisting of an array of four VHF antennas, that will characterize the surface of Jupiter’s icy moon Europa. Measured far-field radiation patterns and reflection coefficients above a 50×80m perfect ground plane are compared to predicted results using commercial software. Normalized Far-Field radiation pattern cuts are measured utilizing a biconical antenna and receiver mounted on a flying drone. Simulations of the radiation pattern show that the number of lobes increases with the distance of the antenna from the ground due to the ground-bounce phenomenon. In addition, the measured and simulated reflection coefficients stay unchanged for height above 2 meters, which matches the predicted calculation of the reactive-near field region of this antenna. However, the antenna electrical parameters are optimized for use on the spacecraft and become detuned in the presence of a ground plane. Simulated and measured return loss discrepancy less than 0.5 dB is presented. Radiation pattern comparisons show similar agreement when taking into account the effect of far-field measurement uncertainties.

Decrossas, Emmanuel↗

Predicting CaO-(MgO)-Al2O3-SiO2 glass reactivity in alkaline environments from force field molecular dynamics simulations

In this investigation, force field-based molecular dynamics (MD) simulations have been employed to generate detailed structural representations for a range of amorphous quaternary CaO-MgO-Al{sub 2}O{sub 3}-SiO{sub 2} (CMAS) and ternary CaO-Al{sub 2}O{sub 3}-SiO{sub 2} (CAS) glasses. Comparison of the simulation results with select experimental X-ray and neutron total scattering and literature data reveals that the MD-generated structures have captured the key structural features of these CMAS and CAS glasses. Based on the MD-generated structural representations, we have developed two structural descriptors, specifically (i) average metal oxide dissociation energy (AMODE) and (ii) average self-diffusion coefficient (ASDC) of all the atoms at melting. Both structural descriptors are seen to more accurately predict the relative glass reactivity than the commonly used degree of depolymerization parameter, especially for the eight synthetic CAS glasses that span a wide compositional range. Hence these descriptors hold great promise for predicting CMAS and CAS glass reactivity in alkaline environments from compositional information.

36 MATERIALS SCIENCE↗

Energy, angle, and temperature dependencies of the sticking of D atoms on Li surfaces

Detailed experimental and computational information on the response of lithium surfaces to irradiation by slow hydrogenic particles (ions, atoms, molecules) is sparse and mainly speculative. In this work, we present a computational study of the reflection and retention of deuterium (D) atoms at crystalline and amorphous lithium surfaces at 300 and 500 K, where the D atoms have an impact energy in the range of 0.025–5 eV and incident angles of 0° (perpendicular incidence) or 85° (near-grazing incidence). Classical molecular dynamics simulations are performed with the reactive bond-order force field (ReaxFF) potentials. This study provides quantitative information on the deuterium sticking probability and recycling coefficient for lithium surfaces. Further, our results support the ongoing work at the Lithium Tokamak eXperiment-β fusion experiment as well as relevant experiments in the laboratory setting.

74 ATOMIC AND MOLECULAR PHYSICS↗

Diffusion-Limited Kinetics in Reactive Systems

A proper representation of chemical kinetics is vital to understanding, modeling, and optimizing many important chemical processes. In liquid and surface phases, where diffusion is slow, the rate at which the reactants diffuse together limits the overall rate of many elementary reactions. Commonly, the textbook Smoluchowski theory is utilized to estimate effective rate coefficients in the liquid phase. On surfaces, modelers commonly resort to much more complex and expensive Kinetic Monte Carlo (KMC) simulations. Here, in this study, we extend the Smoluchowski model to allow the diffusing species to undergo chemical reactions and derive analytical formulas for the diffusion-limited rate coefficients for 3D, 2D, and 2D/3D interface cases. With these equations, we are able to demonstrate that when species react faster than they diffuse they can react orders of magnitude faster than predicted by Smoluchowski theory, through what we term “the reactive transport effect”. We validate the derived steady-state equations against particle Monte Carlo (PMC) simulations, KMC simulations, and non-steady-state solutions. Furthermore, using PMC and KMC simulations, we propose corrections that agree with all limits and the computed data for the 2D and 2D/3D interface steady-state equations, accounting for unique limitations in the associated derived equations. Additionally, we derive equations to handle couplings between diffusion-limited rate coefficients in reaction networks. We believe these equations should make it possible to run much more accurate mean-field simulations of liquids, surfaces, and liquid–surface interfaces accounting for diffusion limitations and the reactive transport effect.

Johnson, Matthew Sean↗

Modeling Stratospheric Constituents: Reactive Species That Regulate Ozone

Photochemical loss of stratospheric ozone occurs primarily by catalytic cycles whose rates are limited by the concentration of OH, HO2, NO2, ClO, and/or BrO as well as the concentration of either atomic oxygen or of ozone itself. Once the concentrations of these gases are established, the photochemical loss rate of O3 depends on the rate coefficient of only a handful of key reactions. We have developed a method for testing our understanding of stratospheric ozone photochemistry by comparing measured and modeled concentrations of reactive hydrogen, nitrogen, chlorine and bromine radicals using a photochemical steady state model constrained by observed concentrations of long-lived precursors (e.g., NO(y), Cl(y), Br(y), O3, H2O, CH4) and environmental parameters such as ozone column, reflectivity, and aerosol surface area. We will show based on analyses of observations obtained by aircraft, balloon, and satellite platforms during the POLARIS campaign that our overall understanding of the processes that regulate these radical species is very good. The most notable current discrepancies are the tendency to underestimate observed NO2 by 15 to 30% for air masses that experience near continuous solar illumination over a 24 hour period and the tendency to underestimate observed OH and H02 by about 10 to 20% during midday and by much larger amounts at high solar zenith angle (SZA > 85). Possible resolutions to these discrepancies will be discussed. This study was carried out in close collaboration with many members of the POLARIS science team.

Salawitch, Ross J.↗

Thermal-energy reactions of O2(2+) ions with O2, N2, CO2, NO, and Ne

The paper presents results of drift-tube mass-spectrometer studies of the reactivity of doubly charged molecular oxygen ions with several molecules and neon atoms. Thermal-energ rate coefficients for the reactions with the molecular reactants were found to be large, close to the limiting Langevin rates. Charge transfer with neon atoms was observed, but the measured rate coefficient was only a small fraction of the Langevin rate. It is concluded that the measured rate constants for the reactions considereed refer to vibrationally excited ions.

Chatterjee, B. K.↗

Temperature-Dependent Kinetics Studies of the Reactions Br((sup 2)P3/2) + H2S yields SH + HBr and Br((sup 2)P3/2) + CH3SH yields CH3S + HBr. Heats of Formation of SH and CH3S Radicals

Time resolved resonance fluorescence detection of Br(sup 2)P3/2) atom disappearance or appearance following 266-nm laser flash photolysis of CF2Br2/H2S/H2/N2, CF2Br2/CH3SH/H2/N2, Cl2CO/H2S/HBr/N2, and CH3SSCH3/HBr/H2/N2 mixtures has been employed to study the kinetics of the reactions Br((sup 2)P3/2) + H2S = SH + HBr (1,-1) and Br((sup2)P3/2) + CH3SH = CH3S + HBr (2, -2) as a function of temperature over the range 273-431K. Arrhenius expressions in units of 10(exp -12) cu cm/molecule/s which describe the results are k1 = (14.2 +/- 3.4) exp[(-2752 +/- 90)/T],(k-1) = (4.40 +/- 0.92) exp[(-971 +/- 73)/T],k(2) = (9.24 +/- 1.15) exp[(-386 +/- 41)/T], and k(-2) = (1.46 +/-0.21) exp[(-399 +/-41)/T; errors are 2 sigma and represent precision only. By examining Br((sup 2)P3/2) equilibrium kinetics following 355nm laser flash photolysis of Br2/CH3SH/H2/N2 mixtures, a 298 K rate coefficient of (1.7 +/- 0.5) x 10(exp -10) cu cm/molecule/s has been obtained for the reaction CH3S + Br2 yields CH3SBr + Br. To our knowledge, these are the first kinetic data reported for each of the reactions studied. Measured rate coefficients, along with known rate coefficients for similar radical + H2S, CH3SH, HBr,Br2 reactions are considered in terms of possible correlations of reactivity with reaction thermochemistry and with IP - EA, the difference between the ionization potential of the electron donor and the electron affinity of the electron acceptor. Both thermochemical and charge-transfer effects appear to be important in controlling observed reactivities. Second and third law analyses of the equilibrium data for reactions 1 and 2 have been employed to obtain the following enthalpies of reaction in units of kcal/mol: for reaction 1, Delta-H(298) = 3.64 +/- 0.43 and Delta-H(0) = 3.26 +/-0.45; for reaction 2, Delta-H(298) = -0.14 +/- 0.28 and Delta-H(0) = -0.65 +/- 0.36. Combining the above enthalpies of reaction with the well-known heats of formation of Br, HBr, H2S, and CH3SH gives the following heats of formation for the RS radicals in units of kcal/mol: Delta-H(sub f)(sub 0)(SH) = 34.07 +/- 0.72, Delta-H(sub f)(sub 298)(SH) = 34.18 +/- 0.68, Delta-H(sub f)(sub 0)(CH3S) = 31.44 +/- 0.54, Delta-H(sub f)(sub 298)(CH3S) = 29.78 +/- 0.44; errors are 2 sigma and represent estimates of absolute accuracy. The SH heat of formation determined from our data agrees well with literature values but has reduced error limits compared to other available values. The CH3S heat of formation determined from our date is near the low end of the range of previous estimates and is 3-4 kcal/mol lower than values derived from recent molecular beam photofragmentation studies.

Nicovich, J. M.↗

The nature and role of quantized transition states in the accurate quantum dynamics of the reaction O + H2 yields OH + H

Accurate quantum mechanical dynamics calculations are reported for the reaction probabilities of O(3P) + H2 yields OH + H with zero total angular momentum on a single potential energy surface. The results show that the reactive flux is gated by quantized transition states up to the highest energy studied, which corresponds to a total energy of 1.90 eV. The quantized transition states are assigned and compared to vibrationally adiabatic barrier maxima; their widths and transmission coefficients are determined; and they are classified as variational, supernumerary of the first kind, and supernumerary of the second kind. Their effects on state-selected and state-to-state reactivity are discussed in detail.

Chatfield, David C.↗

Inelastic, exchange, and reactive processes in rovibrationally excited collisions of HD with H

ABSTRACT The HD molecule is an important coolant in early universe chemistry models and a tracer of H2 in star-forming regions. Rate coefficients for collisional excitation and de-excitation of HD rotational and vibrational levels form important ingredients in astrophysical models. While collisions with He, H2, and H are the most important, available data for H + HD collisions are largely limited to temperatures less than 1000 K for the vibrational ground state, low-lying rotational levels of the v = 1 HD vibrational level, or computed without reactive contributions. Here, through explicit quantum scattering calculations, we report extensive data for rovibrational transitions in HD induced by H atoms for a range of rotational levels in v = 1 and some v = 0 levels for temperatures up to 1000 K. The significance of the computed results for astrophysical modeling is discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion of Atoms in Glassy Mixtures of Deuterium and Lithium

The diffusion coefficients of D and Li were calculated in amorphous, glassy Li+D mixtures for various concentrations of D in lithium over the temperature range 100–1000 K. The densities of the mixtures as a function of temperature were also determined. The diffusion coefficients were obtained by the analysis of the mean-squared displacement using molecular dynamics with ReaxFF, a reactive force field. In conclusion, the diffusion and density data obtained for Li+D mixtures were compared with the available experimental and calculated data reported in the literature.

Diffusion↗

Structural isomers of C2N(+) - A selected-ion flow tube study

Reactivities of the structural isomers CCN(+) and CNC(+) were examined in a selected-ion flow tube at 300 + or - 5 K. The less reactive CNC(+) isomer was identified as the product of the reactions of C(+) + HCN and C(+) + C2N2; in these reactions only CNC(+) can be produced because of energy constraints. Rate coefficients and branching ratios are reported for the reactions of each isomer with H2, CH4, NH3, H2O, C2H2, HCN, N2, O2, N2O, and CO2. Ab initio calculations are presented for CCN(+) and CNC(+); a saddle point for the reaction CCN(+) yielding CNC(+) is calculated to be 195 kJ/mol above CNC(+). The results provide evidence that the more reactive CCN(+) isomer is unlikely to be present in measurable densities in interstellar clouds.

Knight, J. S.↗

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Xenon compatibility in magmatic processes: Hadean to current contexts

Xenon (Xe) behaviour in petrological processes, albeit essential to constrain mantle ingassing and degassing models, is elusive due to its volatile nature, and lack of direct investigation at the pressures (P) and temperatures (T) relevant to magma formation and crystallisation at depth. Xenon stands out amongst noble gases due to its unique reactivity with silicates of the lower crust and upper mantle, which could at least partially explain that published mineral/melt partitioning coefficients span up to six orders of magnitude. We report partition coefficients of Xe using in situ X-ray fluorescence at high P and T, and mass spectrometry analyses. Xenon is found to be moderately incompatible in anorthite-clinopyroxene mix in equilibrium with basalt (partition coefficient value of 0.16 ± 0.06), and compatible in olivine in equilibrium with basalt (partition coefficients in the range 88 ± 22 to 302 ± 46). While Xe is, thus, concentrated in basaltic melts coexisting with crystallising pyroxenes and feldspars, it is strongly retained in olivine at depth. Consequently, Xe originally contained in solid Earth has been preferentially retained at depth throughout Earth’s history, from the magma ocean stages to present day partial mantle melting processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermophysical Property Measurements of Molten Semiconductors in 1-g and Reduced-g Condition

Understanding and controlling the formation kinetics of varieties of crystal imperfections such as point defects, non uniform distribution of doping atoms, and impurity atoms in growing crystals are very important. Theoretical (numerical) modeling of the crystal growth process is an essential step to achieving these objectives. In order to obtain reliable modeling results, input parameters, i.e. various thermophysical parameters, must be accurate. The importance of accurate thermophysical properties of semiconductors in crystal growth cannot be overly emphasized. The total hemispherical emissivity, for instance, has a dramatic impact on the thermal environment. It determines the radiative emission from the surface of the melt which determines to a large extent the profile of the solidified crystal. In order to understand the convection and the turbulence in a melt, viscosity becomes an important parameter. The liquid surface tension determines the shape of the liquid-atmosphere interface near the solid-liquid-atmosphere triple point. Currently used values for these parameters are rather inaccurate, and this program intends to provide more reliable measurements of these thermophysical properties. Thus, the objective of this program is in the accurate measurements of various thermophysical properties which can be reliably used in the modeling of various crystal growth processes. In this program, thermophysical properties of molten semiconductors, such as Si, Ge, Si-Ge, and InSb will be measured as a function of temperature using the High Temperature Electrostatic Levitator at JPL. Each material will be doped by different kinds of impurities at various doping levels. Thermophysical properties which will be measured include: density, thermal expansion coefficient, surface tension, viscosity, specific heat, hemispherical total emissivity, and perhaps electrical and thermal conductivities. Many molten semiconductors are chemically reactive with crucibles. As a result, these dispersed impurities in the melts tend to substantially modify the properties of pure semiconductors. Sample levitation done in a vacuum clearly helps maintain the sample purity. However, in the 1-g environment, all gravity caused effects such as convection, sedimentation and buoyancy are still present in the sample. In addition, large forces needed to levitate a sample in the presence of the gravity can cause additional flows in the melt. The use of the High Temperature Electrostatic Levitator (HTESL) for the present research is a recent development and little is known about the flows induced by the electrostatic forces. In this ground base program, we will define the limits of HTESL technology as various thermophysical properties of molten semiconductors are measured.

Rhim, Won-Kyu↗