Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,045 records · Page 58

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Shayesteh Zadeh, Armin (ORCID:0000000299896710)↗

Cell free biosynthesis of isoprenoids from isopentenol

Cell-free systems are growing in importance for the biosynthesis of complex molecules. These systems combine the precision of traditional chemistry with the versatility of biology in creating superior overall processes. Recently, a new synthetic pathway for the biosynthesis of isoprenoids using the substrate isopentenol, dubbed the isopentenol utilization pathway (IUP), was demonstrated to be a promising alternative to the native 2C-methyl-d-erythritol-4-phosphate (MEP) and mevalonate (MVA) pathways. This simplified pathway, which contains a minimum of four enzymes to produce basic monoterpenes and only depends on ATP and isopentenol as substrates, allows for a highly flexible approach to the commercial synthesis of isoprenoid products. In this work, we use metabolic reconstitution to characterize this new pathway in vitro and demonstrate its use for the cell-free synthesis of mono-, sesquit-, and diterpenoids. Kinetic modeling and sensitivity analysis were also used to identify the most significant parameters for taxadiene productivity, and metabolic control analysis was employed to elucidate protein-level interactions within this pathway, which demonstrated that the IUP enzymatic system is primarily controlled by the concentration and kinetics of choline kinase (CK) and not regulated by any pathway intermediates. We report this is a significant advantage over the natural MEP or MVA pathways as it greatly simplifies future metabolic engineering efforts, both in vitro and in vivo, aiming at improving the kinetics of CK. Finally, we used the insights gathered to demonstrate an in vitro IUP system that can produce 220 mg/L of the diterpene taxadiene, in 9 hr, almost 3-fold faster than any system reported thus far.

59 BASIC BIOLOGICAL SCIENCES↗

Constraining the Reaction Rate of Criegee Intermediates with Carboxylic Acids during the Multiphase Ozonolysis of Aerosolized Alkenes

Criegee intermediates (CIs) are central to the ozonolysis of unsaturated organic compounds, controlling the formation of small fragmentation reaction products and higher molecular weight oligomeric compounds through competing unimolecular and bimolecular pathways. In particular, the reaction of thermalized CI with carboxylic acids (RCOOH) to produce oligomeric α-acyloxyalkyl hydroperoxides (AAHPs) is of interest to the chemical transformation of organic compounds in the atmosphere. In the gas phase, the CI + RCOOH reaction proceeds at supercollisional rates, consistent with a barrierless reaction mechanism. However, the rate of this reaction in condensed phases remains uncertain. Here, we report the results of aerosol flow tube studies of the ozonolysis of submicron organic aerosol containing mixtures of an internal n-alkene (Z-9-tricosene, tri) and a saturated acid (2-hexyldecanoic acid, HDA). The decay of the reactants and formation of reaction products are monitored using atmospheric pressure chemical ionization mass spectrometry (APCI–MS). A 6-fold decrease in the reactive uptake of ozone is observed with increasing acid concentration, consistent with a previous study using alcohol additives. The decay of HDA with ozone exposure shows that HDA scavenges CIs at rates comparable to competing unimolecular pathways involving CIs. The decay kinetics of HDA are analyzed using a simple kinetic model in order to constrain the ratio of the CI + RCOOH and unimolecular CI rate constants. Empirical fitting to this model, combined with theoretical estimates of the unimolecular loss rate, yields a CI + RCOOH rate constant of 1.85 ± 0.27 × 10 –19 cm 3 molec –1 s –1 , which is six orders of magnitude slower than the expected diffusion limit in a liquid organic matrix. Stochastic kinetic simulations are conducted with different values of the CI + RCOOH rate constant to validate this result and show that the diffusion-limited rate is indeed much too rapid to reproduce the experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A quantitative and spatially resolved analysis of the performance-bottleneck in high efficiency, planar hybrid perovskite solar cells

Hybrid perovskites represent a potential paradigm shift for the creation of low-cost solar cells. Current power conversion efficiencies (PCEs) exceed 22%. However, despite this, record PCEs are still far from their theoretical Shockley–Queisser limit of 31%. To increase these PCE values, there is a pressing need to understand, quantify and microscopically model charge recombination processes in full working devices. Here, we present a complete microscopic account of charge recombination processes in high efficiency (18–19% PCE) hybrid perovskite (mixed cation and methylammonium lead iodide) solar cells. We employ diffraction-limited optical measurements along with relevant kinetic modeling to establish, for the first time, local photoluminescence quantum yields, trap densities, trapping efficiencies, charge extraction efficiencies, quasi-Fermi-level splitting, and effective PCE estimates. Correlations between these spatially resolved parameters, in turn, allow us to conclude that intrinsic electron traps in the perovskite active layers limit the performance of these state-of-the-art hybrid perovskite solar cells.

14 SOLAR ENERGY↗

Molecular recognition of a host protein by NS1 of pandemic and seasonal influenza A viruses

The 1918 influenza A virus (IAV) caused the most severe flu pandemic in recorded human history. Nonstructural protein 1 (NS1) is an important virulence factor of the 1918 IAV. NS1 antagonizes host defense mechanisms through interactions with multiple host factors. One pathway by which NS1 increases virulence is through the activation of phosphoinositide 3-kinase (PI3K) by binding to its p85β subunit. Here we present the mechanism underlying the molecular recognition of the p85β subunit by 1918 NS1. Using X-ray crystallography, we determine the structure of 1918 NS1 complexed with p85β of human PI3K. We find that the 1918 NS1 effector domain (1918 NS1 ED ) undergoes a conformational change to bind p85β. Using NMR relaxation dispersion and molecular dynamics simulation, we identify that free 1918 NS1 ED exists in a dynamic equilibrium between p85β-binding–competent and –incompetent conformations in the submillisecond timescale. Moreover, we discover that NS1 ED proteins of 1918 (H1N1) and Udorn (H3N2) strains exhibit drastically different conformational dynamics and binding kinetics to p85β. These results provide evidence of strain-dependent conformational dynamics of NS1. Using kinetic modeling based on the experimental data, we demonstrate that 1918 NS1 ED can result in the faster hijacking of p85β compared to Ud NS1 ED , although the former has a lower affinity to p85β than the latter. Our results suggest that the difference in binding kinetics may impact the competition with cellular antiviral responses for the activation of PI3K. We anticipate that our findings will increase the understanding of the strain-dependent behaviors of influenza NS1 proteins.

59 BASIC BIOLOGICAL SCIENCES↗

Assessment the leaching characteristics and long-term leaching behavior of some radionuclides from synthesized zeolite cement matrix

Highlights: • Synthesis and characterization of zeolite cement synthesized from fly ash • Leaching tests were used to assess the release rate of {sup 85}Sr and {sup 60}Co radionuclides. • Simple mathematical models were used to detect leaching mechanisms. • Long-term leaching behavior of {sup 85}Sr and {sup 60}Co from cementitious matrix was assessed. The development of cementitious materials remains a vital goal to produce valuable products with good mechanical, physical and chemical properties suitable for the safe disposal of concentrated radionuclides resulted from the treatment of contaminated solutions. In this paper, zeolite cement with good properties was laboratory synthesized from industrial by-products fly ash and characterized using various characterization techniques. The leaching characteristics of strontium and cobalt radionuclides from the synthesized zeolite cement were scrutinized conforming with the standard leaching methodology of International Atomic Energy Authority (IAEA). A mechanical strength assessment was executed to characterize the extent of immobilization process of the solidified matrices. The cumulative leaching fraction (φ) of the two studied radionuclides was found to be less than 5% in all examined conditions, which implies the applicability of the IAEA-recommended methodology for estimating the diffusion coefficient. The experimental leaching data were regressed nonlinearly to various mathematical kinetic models to assess the controlling leaching mechanism and to determine the leaching parameters. The regression results indicated that strontium ({sup 85}Sr) and cobalt ({sup 60}Co) leaching resulted from two succeeding mechanisms: as first order kinetic reaction, and then diffusion. The calculated values of leachability indices signify that the performance of each of the studied matrices is within an acceptable range. A simplified mathematical model, rooted in the first order reaction and diffusion mechanisms, was simulated to predict the radionuclides leaching rates from zeolite cement matrices. By comparing the synthesized zeolite cement with other cementitious materials, it can be concluded that the synthesized material can be classified as an efficient material suitable to immobilize {sup 85}Sr and {sup 60}Co from radioactive wastes. The acquired findings demonstrated that the studied immobilized waste matrices have acceptable mechanical effectiveness.

36 MATERIALS SCIENCE↗

Gas Phase Chemical Evolution of Uranium, Aluminum, and Iron Oxides

Abstract We use a recently developed plasma-flow reactor to experimentally investigate the formation of oxide nanoparticles from gas phase metal atoms during oxidation, homogeneous nucleation, condensation, and agglomeration processes. Gas phase uranium, aluminum, and iron atoms were cooled from 5000 K to 1000 K over short-time scales (∆t < 30 ms) at atmospheric pressures in the presence of excess oxygen. In-situ emission spectroscopy is used to measure the variation in monoxide/atomic emission intensity ratios as a function of temperature and oxygen fugacity. Condensed oxide nanoparticles are collected inside the reactor for ex-situ analyses using scanning and transmission electron microscopy (SEM, TEM) to determine their structural compositions and sizes. A chemical kinetics model is also developed to describe the gas phase reactions of iron and aluminum metals. The resulting sizes and forms of the crystalline nanoparticles (FeO-wustite, eta-Al 2 O 3 , UO 2 , and alpha-UO 3 ) depend on the thermodynamic properties, kinetically-limited gas phase chemical reactions, and local redox conditions. This work shows the nucleation and growth of metal oxide particles in rapidly-cooling gas is closely coupled to the kinetically-controlled chemical pathways for vapor-phase oxide formation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

QUANTITATIVELY ASSESSING THE ROLE OF CLOUDS IN THE TRANSMISSION SPECTRUM OF GJ 1214b

Recent observations of the super-Earth GJ 1214b show that it has a relatively featureless transmission spectrum. One suggestion is that these observations indicate that the planet's atmosphere is vertically compact, perhaps due to a water-rich composition that yields a large mean molecular weight. Another suggestion is that the atmosphere is hydrogen/helium-rich with clouds that obscure predicted absorption features. Previous models that incorporate clouds have included their effect without a strong physical motivation for their existence. Here, we present model atmospheres of GJ 1214b that include physically motivated clouds of two types. We model the clouds that are present in chemical equilibrium, as has been suggested to occur on brown dwarfs, which include KCl and ZnS for this planet. We also include clouds that form as a result of photochemistry, forming a hydrocarbon haze layer. We use a photochemical kinetics model to understand the vertical distribution and available mass of haze-forming molecules. We model both solar and enhanced-metallicity cloudy models and determine the cloud properties necessary to match observations. In enhanced-metallicity atmospheres, we find that the equilibrium clouds can match the observations of GJ 1214b if they are lofted high into the atmosphere and have a low sedimentation efficiency (f{sub sed} = 0.1). We find that models with a variety of hydrocarbon haze properties can match the observations. Particle sizes from 0.01 to 0.25 μm can match the transmission spectrum with haze-forming efficiencies as low as 1%-5%.

79 ASTRONOMY AND ASTROPHYSICS↗

Mechanism, thermochemistry, and kinetics of the reversible reactions: C 2 H 3 + H 2 ⇌ C 2 H 4 + H ⇌ C 2 H 5

High-level coupled cluster theory, in conjunction with Active Thermochemical Tables (ATcT) and E,J-resolved master equation calculations, was used in a study of the title reactions, which play an important role in the combustion of hydrocarbons. In the set of radical/radical reactions leading to soot formation in flames, the addition of H-atoms to alkenes is likely a common reaction, triggering the isomerization of complex hydrocarbons to aromatics. The heats of formation of C 2 H 3 , C 2 H 4 , and C 2 H 5 are established to be 301.26 ± 0.30 at 0 K (297.22 ± 0.30 at 298 K), 60.89 ± 0.11 (52.38 ± 0.11), and 131.38 ± 0.22 (120.63 ± 0.22) kJ mol -1 , respectively. The calculated rate constants from first principles agree well with experiments where they are available. Under conditions typical of high temperature combustion – where experimental work is very challenging with a consequent dearth of accurate data –here we provide high-level theoretical results for kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗

Numerical analysis of flows of rarefied gases in long channels with octagonal cross section shapes

Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modelling electron transfer in photosystem I: limits and perspectives

Uncovering the parameters underlying the electron transfer (ET) in photosynthetic reaction centres is of importance for understanding the molecular mechanisms underpinning their functionality. The reductive nature of most cofactors involved in photosynthetic ET makes the direct estimation of their properties difficult. Photosystem I (PSI) operates in a highly reducing regime, making the assessment of cofactor properties even more difficult. Kinetic modelling coupled to a non-adiabatic description of ET is a useful approach in overcoming this hindrance. Here we review the theory and modelling approaches that have been used in assessing parameters associated with ET reactions in PSI, with particular attention to ET reactions involving the phylloquinones and the iron–sulphur clusters. In most modelling studies, the goal is to estimate the driving force of ET, which is usually associated with the cofactor midpoint potentials. The driving force is sensitive to many factors, which define the ET rate, i.e. the reorganisation energy, the coupling with nuclear modes and the electronic matrix elements, which are explored and discussed here. The importance of an inclusive modelling of both forward and reverse ET processes is discussed and highlighted. It is shown that although estimates are indeed sensitive to the exact parameter sets employed in the modelling, a general consensus is still attained, pointing to a scenario where $ΔG^0_{A_{tA→Fx}}/ΔG^0_{A_{tB→Fx}}$ is weakly endergonic/exergonic, respectively. It is emphasised that to further refine those estimates, it will require a joint effort between computational modelling and more wide-ranging experimental studies.

59 BASIC BIOLOGICAL SCIENCES↗

A deep learning approach to programmable RNA switches

Engineered RNA elements are programmable tools capable of detecting small molecules, proteins, and nucleic acids. Predicting the behavior of these synthetic biology components remains a challenge, a situation that could be addressed through enhanced pattern recognition from deep learning. Here, we investigate Deep Neural Networks (DNN) to predict toehold switch function as a canonical riboswitch model in synthetic biology. To facilitate DNN training, we synthesize and characterize in vivo a dataset of 91,534 toehold switches spanning 23 viral genomes and 906 human transcription factors. DNNs trained on nucleotide sequences outperform (R 2 = 0.43–0.70) previous state-of-the-art thermodynamic and kinetic models (R 2 = 0.04–0.15) and allow for human-understandable attention-visualizations (VIS4Map) to identify success and failure modes. This work shows that deep learning approaches can be used for functionality predictions and insight generation in RNA synthetic biology.

59 BASIC BIOLOGICAL SCIENCES↗

Formation of solitary zonal structures via the modulational instability of drift waves

{\rtf1\ansi\ansicpg1252\cocoartf1561\cocoasubrtf600{\fonttbl\f0\fswiss\fcharset0 Helvetica;}{\colortbl;\red255\green255\blue255;\red0\green0\blue0;}{\*\expandedcolortbl;;\cssrgb\c0\c0\c0;}\margl1440\margr1440\vieww10800\viewh8400\viewkind0\pard\tx887\tx1775\tx2662\tx3550\tx4438\tx5325\tx6213\tx7101\tx7988\tx8876\tx9764\tx10651\tx11539\tx12427\tx13314\tx14202\tx15090\tx15977\tx16865\tx17753\tx18640\tx19528\tx20416\tx21303\tx22191\tx23079\tx23966\tx24854\tx25742\tx26629\tx27517\tx28405\tx29292\tx30180\tx31067\tx31955\tx32843\tx33730\tx34618\tx35506\tx36393\tx37281\tx38169\tx39056\tx39944\tx40832\tx41719\tx42607\tx43495\tx44382\tx45270\tx46158\tx47045\tx47933\tx48821\tx49708\tx50596\tx51484\tx52371\tx53259\tx54147\tx55034\tx55922\tx56810\tx57697\tx58585\tx59472\tx60360\tx61248\tx62135\tx63023\tx63911\tx64798\tx65686\tx66574\tx67461\tx68349\tx69237\tx70124\tx71012\tx71900\tx72787\tx73675\tx74563\tx75450\tx76338\tx77226\tx78113\tx79001\tx79889\tx80776\tx81664\tx82552\tx83439\tx84327\tx85215\tx86102\tx86990\tx87877\tx88765\slleading20\pardirnatural\partightenfactor0\f0\fs38 \cf2 The dynamics of the radial envelope of a weak coherent drift wave is approximately governed by a nonlinear Schr\'f6dinger equation, which emerges as a limit of the modified Hasegawa\'97Mima equation. The nonlinear Schr\'f6dinger equation has well-known soliton solutions, and its modulational instability can naturally generate solitary structures. In this paper, we demonstrate that this simple model can adequately describe the formation of solitary zonal structures in the modified Hasegawa\'97Mima equation, but only when the amplitude of the coherent drift wave is relatively small. At larger amplitudes, the modulational instability produces stationary zonal structures instead. Furthermore, we find that incoherent drift waves with beam-like spectra can also be modulationally unstable to the formation of solitary or stationary zonal structures, depending on the beam intensity. Notably, we show that these drift waves can be modeled as quantumlike particles (\'93driftons\'94) within a recently developed phase-space (Wigner\'97Moyal) formulation, which intuitively depicts the solitary zonal structures as quasi-monochromatic drifton condensates. Quantumlike effects, such as diffraction, are essential to these condensates; hence, the latter cannot be described by wave-kinetic models that are based on the ray approximation.\}

Solitons↗

BETO 2021 Peer Review - 3.4.2.302 - Process Scale Up for Production Environments

This project supports the Department of Energy's Bioenergy Technology Office mission of transitioning bioenergy technologies to market by de-risking integration, developing scaling relations and modeling that de-risks scale up, and providing a feedback loop with lower TRL projects to understand process and technology fundamentals. Two objectives were pursued during this merit cycle. First, the Thermal and Catalytic Process Development Unit (TCPDU), a 0.5 ton/day pilot plant, was used to validate kinetic models developed by the Consortium for Computational Physics and Chemistry (CCPC). Validating the models accurately predict product yields and composition is a first step toward de-risking a common industry failure - taking too large of steps between scales. In addition, these models help with design and troubleshooting of new unit operations for the TCPDU. The second objective was to conduct the FY22 verification campaign around ex situ CFP. The primary challenge of this effort was the coordination of multiple national labs and projects. This was overcome by a dedicated leadership role and effective communication strategy between the projects. Currently the project is undergoing a pivot based on a stagegate decision to not conduct the verification campaign in the TCPDU and it is focused on closing out the CFP technology at a smaller scale by the end of the current fiscal year.

BIOMASS FUELS↗

Developing the test conditions for the first fueled transients in the TREAT reactor

The restart of the Transient Reactor Test (TREAT) facility at Idaho National Laboratory, re-establishes a national in-pile testing capability for nuclear fuels in accident scenarios. The purpose of the first tests in TREAT are to study the response of nuclear fuel durring a reactivity initiated accident (RIA). The initial objective is to demonstrate the ability to create a desired enthalpy increase in a representative fuel test specimen over a desired time step and to better understand the uncertainty. The five test rods initially available for these irradiations are UO2-Zircaloy with standard Pressurized Water Reactor (PWR) geometry. The rods are 15 cm long with 9 cm of active fuel. The initial testing environment is dry helium. BISON fuel performance calculations were run to model the heat transfer during postulated transients. Based on those results energy deposition targets were selected around the range of observable phenomena in the testing environment. Reactor / experiment power coupling factors, step insertion size and transient clipping response times were then calculated using results from a RELAP point kinetics model. Calculated energy deposition targets are 150, 350, 650, 960, and 1300 J/gUO2. The energy coupling factor for the experiments was calculated to be ~1.5 J/gUO2-MJ. Two step insertion sizes were selected: a 0.6% ?k/k step for the 150 J/g transient and a 1.8% ?k/k step for the other transients. Clipping response times ranged from 1.3 seconds to unclipped transients. The experiments will be irradiated in the TREAT facility in the fall of 2018. The first transient will be the 150 J/g pulse with a 0.6% ?k/k step that is unclipped. Follow on tests may modify the reactor operating parameters based on results from previous transients. Any preliminary data available on the results of the transients at the time of the NUMAT Conference will be presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Point kinetics calculations with fully coupled thermal fluids reactivity feedback

The point kinetics model has been widely used in the analysis of the transient behavior of a nuclear reactor. In the traditional nuclear reactor system safety analysis codes such as RELAP5, the reactivity feedback effects are calculated in a loosely coupled fashion through operator splitting approach. This paper discusses the point kinetics calculations with the fully coupled thermal fluids and fuel temperature feedback implemented into the RELAP-7 code currently being developed with the MOOSE framework. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Detailed Reaction Mechanism for Thiosulfate Oxidation by Ozone in Aqueous Environments

The ozone oxidation, or ozonation, of thiosulfate is an important reaction for wastewater processing, where it is used for remediation of mining effluents, and for studying aerosol chemistry, where its fast reaction rate makes it an excellent model reaction. Although thiosulfate ozonation has been studied since the 1950s, challenges remain in developing a realistic reaction mechanism that can satisfactorily account for all observed products with a sequence of elementary reaction steps. Here, we present novel measurements using trapped microdroplets to study the pH-dependent thiosulfate ozonation kinetics. We detect known products and intermediates, including SO32-, SO42-, S3O62-, and S4O62-, establishing agreement with the literature. However, we identify S2O42- as a new reaction intermediate and find that the currently accepted mechanism does not directly explain observed pH effects. Thus, we develop a new mechanism, which incorporates S2O42- as an intermediate and uses elementary steps to explain the pH dependence of thiosulfate ozonation. The proposed mechanism is tested using a kinetic model benchmarked to the experiments presented here, then compared to literature data. We demonstrate good agreement between the proposed thiosulfate ozonation mechanism and experiments, suggesting that the insights in this paper can be leveraged in wastewater treatment and in understanding potential climate impacts.

Deal, Alexandra M↗