Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reaction Rates”

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 667 records · Page 37

Gamma-ray measurements for inertial confinement fusion applications

High-energy γ rays generated from inertial confinement fusion (ICF) experiments have become an important signature for studying the dynamics of implosion processes. Due to their high-energy and penetrating nature, γ rays are the most unperturbed fusion products, which can preserve the original birth information of the fusion process. Fusion γ rays provide a direct measure of nuclear reaction rates (unlike x rays) without being compromised by Doppler spreading (unlike neutrons). However, unambiguous γ-ray measurements for ICF study further required a decade-long period of technological development, which included a deepening understanding of fusion γ-ray characteristics and innovations in instrument performance. This review article introduces the production mechanism of the prompt and secondary γ rays and various ICF performance parameters (e.g., bang time and burn width), which can be derived from γ-ray measurement. A technical overview will be followed by summarizing γ-ray detectors fielded or proposed, especially for high-yield ICF experiments at the Omega Laser Facility and National Ignition Facility. Over the past few years, γ-ray diagnostic technologies have been extended beyond ICF research. A few examples of non-ICF applications of γ-ray detectors are introduced at the end of this article.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Temperature-based reactive flow model for triaminotrinitrobenzene (TATB) plastic bonded explosives

A new reactive flow model is presented for triaminotrinitrobenzene (TATB)-based plastic bonded explosives, applicable to shock initiation and steady detonation problems of differing initial temperature. Temperature disequilibrium is assumed between unreacted explosive, material in the vicinity of compressed defects (called hot spots), and reaction products. The model incorporates temperature-dependent decomposition reaction rates. Particularly, Arrhenius model parameters were derived from quantum-based molecular dynamics simulations of TATB decomposition. Further, a model of detonation carbon aggregation is incorporated, describing the slow release of energy inherent to detonation in TATB-based materials. Model parameters were calibrated against gas gun shock initiation experiments and steady detonation rate stick tests. The predictive ability of the model in the shock initiation regime is tested against recent thin pulse experiments. The model is found to perform equally well in predicting the size-effect curve of ambient, cold, and hot rate sticks. The present work demonstrates the viability of incorporating results from subscale simulations into a continuum-scale reactive flow model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structure–property relations of binary ferrite melts

Molten ferrite systems are used in the smelting and refining processes in steelmaking, to reduce the loss of metals in slags and to accelerate reaction rates. Here, high-energy x-ray diffraction experiments have been performed on aerodynamically levitated molten spheres of 43BaO–57FeO X and 43SrO–57FeO X at 1873 K using laser beam heating. The composition was varied within the range of x = 1–1.5 by changing the oxygen partial pressure of the levitation gas. The corresponding x-ray pair distribution functions have been interpreted using empirical potential structure refinement (EPSR) modeling. In oxygen-rich melts (x = 1.5), our EPSR models indicate very similar structures for the different alkaline-earth liquids, with both the Ba–O and Sr–O coordination numbers to be ∼8.4 and the total Fe–O coordination numbers ∼5.7. However, our models show that in reducing environments, the Fe 3+ and Fe 2+ ions exhibit very different behaviors in the Ba- and Sr-ferrite liquids. In the Ba-ferrite melt, the Fe 3+ –O coordination number decreases from 5.7 (at x = 1.5) to 5.2 (at x = 1.07), whereas Fe 2+ –O remains constant at ∼5.0 across the same compositional range. In the Sr melts, both the Fe 2+ –O and Fe 3+ –O coordination numbers rise from ∼5.7 (at x = 1.5) to 6.3 (at x = 1.07). All models show the structures to be heterogeneous with intertwined nanometer sized clusters or channels of Ba/Sr–O and Fe–O polyhedra that grow as oxygen content is reduced. Changes in the viscosity and electrical properties are interpreted in terms of the number of bridging and non-bridging oxygens associated with FeO 4 tetrahedra and concentration of charge carriers, respectively.

Benmore, Chris J. [Argonne National Laboratory (AN↗

Chemical pumps and flexible sheets spontaneously form self-regulating oscillators in solution

The synchronization of self-oscillating systems is vital to various biological functions, from the coordinated contraction of heart muscle to the self-organization of slime molds. Through modeling, we design bioinspired materials systems that spontaneously form shape-changing self-oscillators, which communicate to synchronize both their temporal and spatial behavior. Here, catalytic reactions at the bottom of a fluid-filled chamber and on mobile, flexible sheets generate the energy to “pump” the surrounding fluid, which also transports the immersed sheets. The sheets exert a force on the fluid that modifies the flow, which in turn affects the shape and movement of the flexible sheets. This feedback enables a single coated (active) and even an uncoated (passive) sheet to undergo self-oscillation, displaying different oscillatory modes with increases in the catalytic reaction rate. Two sheets (active or passive) introduce excluded volume, steric interactions. This distinctive combination of the hydrodynamic, fluid–structure, and steric interactions causes the sheets to form coupled oscillators, whose motion is synchronized in time and space. We develop a heuristic model that rationalizes this behavior. These coupled self-oscillators exhibit rich and tunable phase dynamics, which depends on the sheets’ initial placement, coverage by catalyst and relative size. Moreover, through variations in the reactant concentration, the system can switch between the different oscillatory modes. This breadth of dynamic behavior expands the functionality of the coupled oscillators, enabling soft robots to display a variety of self-sustained, self-regulating moves.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the mechanisms of ion adsorption to aqueous interfaces: air-water vs. oil-water

The adsorption of ions to water-hydrophobe interfaces influences a wide range of phenomena, including chemical reaction rates, ion transport across biological membranes, and electrochemical and many catalytic processes; hence, developing a detailed understanding of the behavior of ions at water-hydrophobe interfaces is of central interest. Here, we characterize the adsorption of the chaotropic thiocyanate anion (SCN − ) to two prototypical liquid hydrophobic surfaces, water-toluene and water-decane, by surface-sensitive nonlinear spectroscopy and compare the results against our previous studies of SCN − adsorption to the air-water interface. For these systems, we observe no spectral shift in the charge transfer to solvent spectrum of SCN − , and the Gibb’s free energies of adsorption for these three different interfaces all agree within error. We employed molecular dynamics simulations to develop a molecular-level understanding of the adsorption mechanism and found that the adsorption for SCN − to both water-toluene and water-decane interfaces is driven by an increase in entropy, with very little enthalpic contribution. This is a qualitatively different mechanism than reported for SCN − adsorption to the air-water and graphene-water interfaces, wherein a favorable enthalpy change was the main driving force, against an unfavorable entropy change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A theory of phonon-induced friction on molecular adsorbates

In this manuscript, we provide a general theory for how surface phonons couple to molecular adsorbates. Our theory maps the extended dynamics of a surface’s atomic vibrational motions to a generalized Langevin equation, and by doing so captures these dynamics in a single quantity: the non-Markovian friction. The different frequency components of this friction are the phonon modes of the surface slab weighted by their coupling to the adsorbate degrees of freedom. Using this formalism, we demonstrate that physisorbed species couple primarily to acoustic phonons while chemisorbed species couple to dispersionless local vibrations. We subsequently derive equations for phonon-adjusted reaction rates using transition state theory and demonstrate that these corrections improve agreement with experimental results for CO desorption rates from Pt(111).

Science & Technology - Other Topics↗

Cross-Code Verification of Neutronics Analysis Tools at INL Applied for 238 Pu Production in the Advanced Test Reactor

Here, analyses are completed for experiments prior to experiment irradiation in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL). Various codes are used to qualify all experiments planned for insertion in the reactor, thereby ensuring that all safety and programmatic requirements are satisfied preirradiation. Among the common experiment analysis tools at INL are MCNP5 coupled to ORIGEN2 (MOPY) and MC21. MOPY uses MCNP5 for transport calculations along with calculations for fluxes and select reaction rates, and then ORIGEN2 handles the step-by-step and postirradiation depletion. MC21 handles all in-reactor transport and step-by-step, during-irradiation, depletion calculations, and then ORIGEN (SCALE 6.2.3) is used for decay and dose calculations postirradiation. The MOPY results, along with those obtained via two variations of the MC21 model, were compared in terms of 238 Pu production in the ATR’s H10 position. For the MOPY model, the MC21 model utilizing the HELIOS-based fission product (FP) library, and the MC21 model utilizing the expanded 1300 FP library, the during-cycle irradiation in-core heating results were sufficiently equivalent; however, the MOPY model and the MC21 model with the HELIOS library showed some differences relating to the respective FP libraries. Ultimately, the MC21 model with a 1300 FP library produced the most consistent results throughout the cycle, whereas the MC21 model that utilized the (smaller) HELIOS library was able to handle during-irradiation analysis but lacked certain short-lived FPs that significantly contributed to the total decay heat at shutdown. MOPY, on the other hand, was found to overpredict fission gas production, as a result of limitations in the ORIGEN2 code.

ATR↗

A Monte Carlo Laplace Transform Estimator for Radiation Transport

This work formulates and implements a Laplace transform estimator in a simple Monte Carlo radiation transport code. The estimator maps flux-based quantities of interest, like reaction rates, from a desired phase-space dimension to the complex Laplace domain. This on-the-fly Monte Carlo integration technique enables the spectral analysis of arbitrary nuclear systems via the Laplace transform. A simple code tests the estimator in neutron slowing-down problems across various infinite media, and the results compare well with Ganapol’s uninverted analytical solution of the neutron slowing-down equation.

97 MATHEMATICS AND COMPUTING↗

Vented and sealed cookoff of powdered and pressed ε-CL-20

In this work, we have completed a series of vented and sealed cookoff experiments of the ε-polymorph of CL-20 in our Sandia Instrumented Thermal Ignition (SITI) apparatus using both powder and pressed pellets at nominal densities of 313 ± 8 kg/m 3 and 1030 ± 4 kg/m 3 , respectively. The boundary temperature of our aluminum confinement cylinder was ramped in 10 minutes from room temperature to a prescribed set-point temperature ranging between 448 nd 468 K and held at the set-point temperature until ignition. A universal cookoff model (UCM) has been calibrated using the ε-CL-20 SITI data to predict pressurization and thermal ignition of ε-CL-20. The ignition model was validated by using one-dimensional time-to-explosion (ODTX) ignition data from a different laboratory. We found that a thirtyfold increase in the reaction rates due to liquefaction at 520 K could explain the high temperature ODTX cookoff data. The model gives a plausible explanation of why melting is important in fast cookoff events involving CL-20. Our model also gives support to 520 K as the liquefaction point of CL-20, which has different values in the literature.

36 MATERIALS SCIENCE↗

Importance of gas heating in capacitively coupled radiofrequency plasma-assisted synthesis of carbon nanomaterials

In pursuit of diamond nanoparticles, a capacitively-coupled radio frequency flow-through plasma reactor was operated with methane-argon gas mixtures. Signatures of the final product obtained microscopically and spectroscopically indicated that the product was an amorphous form of graphite. This result was consistent irrespective of combinations of the macroscopic reactor settings. To explain the observed synthesis output, measurements of C 2 and gas properties were carried out by laser-induced fluorescence and optical emission spectroscopy. Strikingly, the results indicated a strong gas temperature gradient of 100 K per mm from the center of the reactor to the wall. Based on additional plasma imaging, a model of hot constricted region (filamentation region) was then formulated. It illustrated that, while the hot constricted region was present, the bulk of the gas was not hot enough to facilitate diamond sp 3 formation: characterized by much lower reaction rates, when compared to sp 2 , sp 3 formation kinetics are expected to become exponentially slow. This result was further confirmed by experiments under identical conditions but with a H 2 /CH 4 mixture, where no output material was detected: if graphitic sp 2 formation was expected as the main output material from the methane feedstock, atomic hydrogen would then be expected to etch it away in situ, such that the net production of that sp 2 -hybridized solid material is nearly a zero. Finally, the crucial importance of gas heating was corroborated by replacing RF with microwave source whereby facile sp 3 production was attained with H 2 /CH 4 gas mixture.

36 MATERIALS SCIENCE↗

16 O + 16 O fusion at sub-coulomb energies

A simple optical model (OPM) method using non-monotonic (NM) potentials characterized by a repulsive core from the microscopic theory of the Pauli-led energy-density functional (EDF) has been developed to investigate the 16 O + 16 O fusion at sub-Coulomb energies relevant to the oxygen burning. Here, the study involves the analysis of the experimental fusion cross-section (FCS) data in energy range 6.92 ≤ E cm ≤ 13.83 MeV, which includes the Coulomb barrier region E cm = 10.0–11.92 MeV. Apart from an excellent description of the existing FCS data in the studied energy range, the associated hindrance, characterized by the S-factor reaching a maximum and then gradually falling off at lower energies, so far suggested empirically for the system, is reproduced down to 4 MeV for the first time in the simple OPM. An S-factor maximum of S 0 = 3.15 × 10 25 MeV.b at E 0 = 7.14 MeV is observed with T ≃ 2.6 GK, which conforms to the values reported for quiescent and explosive burning. Our reaction rate, deduced from the NM potential, compares well with the Caughlan and Fowler data. Dominant partial waves implicit in the observed maximum S-factor in the studied Gamow range are also explored. Our present findings, with the success of NM potentials, suggest that the nucleus–nucleus potential is non-monotonic.

16O+16O↗

Fabrication and field emission properties of vertical, tapered GaN nanowires etched via phosphoric acid

The controlled fabrication of vertical, tapered, and high-aspect ratio GaN nanowires via a two-step top-down process consisting of an inductively coupled plasma reactive ion etch followed by a hot, 85% H 3 PO 4 crystallographic wet etch is explored. The vertical nanowires are oriented in the $[0001]$ direction and are bound by sidewalls comprising of $\{33\bar{6}2\}$ semipolar planes which are at a 12° angle from the [0001] axis. High temperature H3PO4 etching between 60 °C and 95 °C result in smooth semipolar faceting with no visible micro-faceting, whereas a 50 °C etch reveals a micro-faceted etch evolution. High-angle annular dark-field scanning transmission electron microscopy imaging confirms nanowire tip dimensions down to 8–12 nanometers. The activation energy associated with the etch process is 0.90 ± 0.09 eV, which is consistent with a reaction-rate limited dissolution process. The exposure of the $\{33\bar{6}2\}$ type planes is consistent with etching barrier index calculations. The field emission properties of the nanowires were investigated via a nanoprobe in a scanning electron microscope as well as by a vacuum field emission electron microscope. Finally, the measurements show a gap size dependent turn-on voltage, with a maximum current of 33 nA and turn-on field of 1.92 V nm -1 for a 50 nm gap, and uniform emission across the array.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Present understanding of ignition and gain using indirect-drive inertial confinement fusion target designs on the U.S. National Ignition Facility

For many decades, the running joke in fusion research has been that 'fusion' is thirty years away and always will be. Yet, these past few years we find ourselves in a position where we can now talk about the milestones of burning plasmas, fusion ignition, and target energy gain greater than unity (scientific breakeven) in the past tense. Fusion is no longer a joke! Yet getting to fusion ignition, the tipping-point of thermonuclear instability resulting in an explosive increase in ion thermal temperature and fusion reaction-rate, and scientific breakeven (target gain, $G_{target} =$ fusion yield/deposited laser energy >1, in the laser-driven inertial confinement fusion context) has not been easy. Here, in this publication, we discuss our present understanding of the physics and technological challenges surrounding ignition and Gain as well as highlight some outstanding problems that still need resolution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Characterization of uncertainties in electron-argon collision cross sections

Abstract The predictive capability of a plasma discharge model depends on accurate representations of electron-impact collision cross sections, which determine the corresponding reaction rates and electron transport properties. The values of cross sections can be known only approximately either through experiments or simulations and are thus subject to uncertainties. Quantifying the uncertainties in plasma simulations allows us to assess the reliability of simulations and to provide a basis for interpreting discrepancies between simulations and experiments. For such uncertainty quantification of plasma simulations, it is essential to quantify the uncertainties of the underlying cross sections. Although much effort has been committed to calibrate the cross section values, their uncertainties are not well investigated. We characterize uncertainties in electron-argon atom collision cross sections using a Bayesian framework. Six collision processes—elastic momentum transfer, ionization, and four excitations—are characterized with semi-empirical models, which effectively capture the features important to the macroscopic properties of the plasma. A probability model for the uncertain parameters of these semi-empirical models is developed. Specifically, a Gaussian-process likelihood model is proposed to capture discrepancies among data sets, as well as the model-form inadequacies of the semi-empirical models. Two other likelihood models are compared with the proposed Gaussian-process model, to illustrate the importance of the choice of the likelihood model. The cross section models are calibrated using the electron-beam experiments and ab-inito quantum simulations. The resulting calibrated uncertainties capture well the scattering among the data sets. The calibrated cross section models are further validated against swarm-parameter experiments and zero-dimensional Boltzmann equation simulations of widely used cross section datasets.

Chung, Seung Whan (ORCID:0000000302501549)↗

Successful νp -process in neutrino-driven outflows in core-collapse supernovae

The origin of the solar system abundances of several proton-rich isotopes, especially 92,94 Mo and 96,98 Ru, has been an enduring mystery in nuclear astrophysics. An attractive proposal to solve this problem is the νp-process, which can operate in neutrino-driven outflows in a core-collapse supernova after the shock is launched. Years of detailed studies, however, have cast doubt over the ability of this process to generate sufficiently high absolute and relative amounts of various p-nuclei. The νp-process is also thought to be excluded by arguments based on the long-lived radionuclide 92 Nb. Here, we present explicit calculations, in which both the abundance ratios and the absolute yields of the p-nuclei up to A ≲ 105 are successfully reproduced, even when using the modern (medium enhanced) triple-α reaction rates. The process is also shown to produce the necessary amounts of 92 Nb. The models are characterized by subsonic outflows and by the protoneutron star masses in the ≳ 1.7 M ⊙ range. In conclusion, this suggests that the Mo and Ru p-nuclides observed in the Solar System were made in CCSN explosions characterized by an extended accretion stage.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of the inter-species power balance in JET plasmas

Most auxiliary heating methods provide heating to more than one particle species (electrons, ions, impurities) in a fusion plasma. This can lead to substantial temperature differences between species, depending on conditions such as heating power to the different species and collisionality, with temperature differences between species limited by inter-species thermal equipartition and transport. The analysis of the steady-state electron-ion and impurity-ion power balances presented in this paper are used for consistency-checking experimental ion and electron temperature measurements and for inferring the main ion temperature from measured impurity temperatures. As ion temperature measurements by charge exchange spectroscopy (CXS) based on impurity ions have become more difficult and time-consuming since the installation of the ITER-like wall (ILW) with Be and W PFC's, knowing the maximum sustainable temperature difference between ions and electrons, |T i -T e | allows rejecting erroneous measurements. It also obviates the need for an ion temperature measurement, if an electron temperature measurement is available and |T i - T e | cannot be larger than the combined errors of the underlying measurements. A power balance analysis is also required for estimating the errors of the ion and electron heat fluxes prior to any species-resolved transport analysis. The ion-impurity temperature differences are usually found to be small due to strong thermal equipartition between ion species. However, they can approach 10% in JET-ILW low density, high power discharges, such the ones under development for a future JET deuterium–tritium campaign (Joffrin et al 2019 Nucl. Fusion). This has a generally small, but not always negligible effect on the calculation of fusion reaction rates, which depend on main ion temperatures. An important outcome of this analysis is that temperature differences between impurity species are always much smaller than between the impurities and hydrogenic species and can usually be neglected. The article presents two methods for calculating the impurity-to-main ion temperature ratio. Finally, this analysis leads to a method for the reconstruction ion temperature profiles from ion temperature data available at only one or a small number of spatial locations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear burn control in ITER using adaptive allocation of actuators with uncertain dynamics

Abstract ITER will be the first tokamak to sustain a fusion-producing, or burning, plasma. If the plasma temperature were to inadvertently rise in this burning regime, the positive correlation between temperature and the fusion reaction rate would establish a destabilizing positive feedback loop. Careful regulation of the plasma’s temperature and density, or burn control, is required to prevent these potentially reactor-damaging thermal excursions, neutralize disturbances and improve performance. In this work, a Lyapunov-based burn controller is designed using a full zero-dimensional nonlinear model. An adaptive estimator manages destabilizing uncertainties in the plasma confinement properties and the particle recycling conditions (caused by plasma–wall interactions). The controller regulates the plasma density with requests for deuterium and tritium particle injections. In ITER-like plasmas, the fusion-born alpha particles will primarily heat the plasma electrons, resulting in different electron and ion temperatures in the core. By considering separate response models for the electron and ion energies, the proposed controller can independently regulate the electron and ion temperatures by requesting that different amounts of auxiliary power be delivered to the electrons and ions. These two commands for a specific control effort (electron and ion heating) are sent to an actuator allocation module that optimally maps them to the heating actuators available to ITER: an electron cyclotron heating system (20 MW), an ion cyclotron heating system (20 MW), and two neutral beam injectors (16.5 MW each). Two different actuator allocators are presented in this work. The first actuator allocator finds the optimal mapping by solving a convex quadratic program that includes actuator saturation and rate limits. It is nonadaptive and assumes that the mapping between the commanded control efforts and the allocated actuators (i.e. the effector model) contains no uncertainties. The second actuator allocation module has an adaptive estimator to handle uncertainties in the effector model. This uncertainty includes actuator efficiencies, the fractions of neutral beam heating that are deposited into the plasma electrons and ions, and the tritium concentration of the fueling pellets. Furthermore, the adaptive allocator considers actuator dynamics (actuation lag) that contain uncertainty. This adaptive allocation algorithm is more computationally efficient than the aforementioned nonadaptive allocator because it is computed using dynamic update laws so that finding the solution to a static optimization problem is not required at every time step. A simulation study assesses the performance of the proposed adaptive burn controller augmented with each of the actuator allocation modules.

Physics↗