Search NASA⌕ Search

SEARCH · Search NASA

Results for “Activity coefficient”

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 955 records · Page 53

NORAD satellite tracking

NORAD routinely tracks about 6000 orbiting objects. During the last 30 days of orbital time, prior to reentry, special perturbations are used in the orbital update procedure. Besides routine orbit determination, NORAD does special tasks such as predicting satellite orbit conjunctions within 20 km, ephimerides of weather satellites, satellite decay predictions and other studies. Since their mission is operational, they do not store the data from their analyses. The ballistic coefficient is not known for most of the orbiting objects. If a ballistic coefficient were derived that was consistent with one density model, it might give erroneous results if used with a different density model. Given the ballistic coefficient, density values could, in principle, be obtained from their tracking data. The densities would represent an integrated mean over the orbital path near perigee. They would be model dependent and would not necessarily represent the real density. In summary, the primary need is for reliable forecasts of solar flux (F10.7) and geomagnetic activity (Ap) in the 1 to 4 week time scale. Forecasts over longer time spans would also be useful for special projects.

Liu, Joseph J. F.↗

Atmospheric drag model calibrations for spacecraft lifetime prediction

Although solar activity prediction uncertainty normally dominates decay prediction error budget for near-Earth spacecraft, the effect of drag force modeling errors for given levels of solar activity needs to be considered. Two atmospheric density models, the modified Harris-Priester model and the Jacchia-Roberts model, to reproduce the decay histories of the Solar Mesosphere Explorer (SME) and Solar Maximum Mission (SMM) spacecraft in the 490- to 540-kilometer altitude range were analyzed. Historical solar activity data were used in the input to the density computations. For each spacecraft and atmospheric model, a drag scaling adjustment factor was determined for a high-solar-activity year, such that the observed annual decay in the mean semimajor axis was reproduced by an averaged variation-of-parameters (VOP) orbit propagation. The SME (SMM) calibration was performed using calendar year 1983 (1982). The resulting calibration factors differ by 20 to 40 percent from the predictions of the prelaunch ballistic coefficients. The orbit propagations for each spacecraft were extended to the middle of 1988 using the calibrated drag models. For the Jaccia-Roberts density model, the observed decay in the mean semimajor axis of SME (SMM) over the 4.5-year (5.5-year) predictive period was reproduced to within 1.5 (4.4) percent. The corresponding figure for the Harris-Priester model was 8.6 (20.6) percent. Detailed results and conclusions regarding the importance of accurate drag force modeling for lifetime predictions are presented.

Binebrink, A. L.↗

Nonlinear ARMA models for the D(st) index and their physical interpretation

Time series models successfully reproduce or predict geomagnetic activity indices from solar wind parameters. A method is presented that converts a type of nonlinear filter, the nonlinear Autoregressive Moving Average (ARMA) model to the nonlinear damped oscillator physical model. The oscillator parameters, the growth and decay, the oscillation frequencies and the coupling strength to the input are derived from the filter coefficients. Mathematical methods are derived to obtain unique and consistent filter coefficients while keeping the prediction error low. These methods are applied to an oscillator model for the Dst geomagnetic index driven by the solar wind input. A data set is examined in two ways: the model parameters are calculated as averages over short time intervals, and a nonlinear ARMA model is calculated and the model parameters are derived as a function of the phase space.

Vassiliadis, D.↗

Evaluation of two fiber optic-based solar collection and distribution systems for advanced space life support

Growing plants in an enclosed controlled environment is crucial in developing bioregenerative life-support systems (BLSS) for space applications. The major challenge currently facing a BLSS is the extensive use of highly energy-intensive electric light sources, which leads to substantial energy wastes through heat dissipations by these lamps. An alternative lighting strategy is the use of a solar irradiance collection, transmission, and distribution system (SICTDS). Two types of fiber optic-based SICTDS, a Fresnel-lens Himawari and a parabolic-mirror optical waveguide (OW) lighting system, were evaluated. The overall efficiency for the OW SICTDS of 40.5% exceeded by 75% that for the Himawari of 23.2%. The spectral distributions of the light delivered by the Himawari and the OW SICTDS were almost identical and had practically no difference from that of terrestrial solar radiation. The ratios of photosynthetically active radiation (PAR) to total emitted radiation (k) of 0.39 +/- 0.02 for the Himawari and 0.41 +/- 0.04 for the OW SICTDS were statistically indistinguishable, were not significantly different from that of 0.042 +/- 0.01 for terrestrial solar radiation, and were comparable to that of 0.35 for a high-pressure sodium (HPS) lamp. The coefficients of variation (CV) of 0.34 and 0.39 for PPF distributions, both at 50 mm X 50 mm square grid arrays, corresponding to the Himawari and the OW SICTDS, respectively, were comparable with each other but were both significantly greater than the CV of 0.08 corresponding to the HPS lamp. The average fresh weight or dry weight of lettuce grown in the solar chamber with either the Himawari or the OW SICTDS showed no statistical difference from the average fresh weight or dry weight of lettuce grown in the reference chamber with the HPS lamp. The results of this study suggest that an SICTDS could help reduce the electric power demand in a BLSS.

NASA Discipline Life Support Systems↗

Leveraging Pre-storm Soil Moisture Estimates for Enhanced Land Surface Model Calibration in Ungauged Hydrologic Basins

Despite long-standing efforts, hydrologists still lack robust tools for calibrating land surface model (LSM) streamflow estimates within ungauged basins. Using surface soil moisture estimates from the Soil Moisture Active Passive Level 4 Soil Moisture (L4_SM) product, precipitation observations, and streamflow gauge measurements for 617 medium-scale (200-10,000 km2) basins in the contiguous United States, we measure the temporal (Spearman) rank correlation between antecedent (i.e., pre-storm) surface soil moisture (ASM) and the storm-scale runoff coefficient (RC; the fraction of storm-scale precipitation accumulation converted into streamflow). In humid and semi-humid basins, this rank correlation is shown to be sufficiently strong to allow for the substitution of storm-scale RC observations (available only in basins that are both lightly regulated and gauged) with high-quality ASM values (available quasi-globally from L4_SM) in streamflow calibration procedures. Using this principle, we define a new, basin-wise LSM streamflow calibration approach based on L4_SM alone and successfully apply it to identify LSM configurations that produce a high rank correlation with observed RC. However, since the approach cannot detect RC bias, it is less successful in identifying LSM configurations with low mean-absolute error. Plain Text Summary Accurately forecasting the fraction of rainfall that runs off into streams, as opposed to infiltrates into the soil, is critical for flash-flood prediction, water-resource monitoring, and tracking the transport of nutrients from agricultural fields into local waterways. Such forecasting is typically performed by hydrologic models that attempt to represent the physical processes responsible for surface runoff generation. However, to provide accurate streamflow forecasts, these models typically need to be calibrated against actual streamflow observations. This is problematic given the relatively poor, and declining, global availability of stream gauges. This paper presents a novel model calibration strategy that uses soil moisture from remote sensing and numerical modeling in place of streamflow observations during calibration. This transition has significant practical advantages because, unlike streamflow observations, the soil moisture data are continuously available across space. Our results demonstrate that this new approach can significantly improve hydrologic models within humid and semi-humid basins lacking sufficient ground-based instrumentation for traditional streamflow calibration.

hydrologic models↗

Durable, Impermeable Brazes for Solid Oxide Fuel Cells (Final Report)

To develop a new braze to replace the current state-of-the-art Ag-CuO reactive air braze (RAB), two approaches were proposed: 1) development of a new, self-passivating, silver-free braze and 2) development of a CuO-free, silver-based braze. For the first approach, an integrated computational-experimental method was used. Thermo-Calc® was utilized to quickly screen through hundreds of alloy systems to identify candidate braze compositions with the appropriate melting ranges. Material compositions, melting ranges, mechanical properties, oxidation resistance, and wetting characteristics of the candidate alloys were then analyzed experimentally. Substrate surface pre-treatments, active element additions, and novel brazing schemes were also investigated. Unfortunately, the current study failed to identify new silver-free braze compositions suitable for SOFC applications. To help measure the passivation characteristics of surface passivating braze oxide coatings, work was conducted to demonstrate, for the first time, that a variety of disparate and technologically-relevant thermal, mechanical, and electrochemical oxygen-exchange thin film material properties could all be obtained from in situ , current-collector-free wafer curvature measurements. Specifically, temperature or oxygen partial pressure induced changes in the curvature of 230 nm thick (100)-oriented Pr 0.1 Ce 0.9 O 1.95_x (10PCO) films atop 200 mm thick single crystal yttria stabilized zirconia or magnesium oxide substrates were used to measure the biaxial modulus, Young’s modulus, thermal expansion coefficient, thermochemical expansion coefficient, oxygen nonstoichiometry, chemical oxygen surface exchange coefficient, oxygen surface exchange resistance, thermal stress, chemical stress, thermal strain, and chemical strain of the model mixed ionic electronic conducting material 10PCO. The (100)-oriented thin film 10PCO thermal expansion coefficient, thermo-chemical expansion coefficient, oxygen nonstoichiometry, and Young’s modulus (which is essentially constant, at ~200 MPa, over the entire 280–700 °C temperature range in air) measured here were similar to those from other bulk and thin film 10PCO studies. In addition, the measured PCO10 oxygen surface coefficients were in agreement with those reported by other in situ , current collector- free techniques. However, the measured PCO10 oxygen surface exchange coefficients were significantly lower than those obtained from literature studies with large amounts of intentional or inadvertent precious metal surface coverage/contamination (suggesting that uncontaminated (100)-oriented 10PCO may not be a desirable SOFC cathode material). Taken together, this highlighted the advantages of using a sample’s mechanical response, instead of the more traditional electrical response, to probe the electrochemical properties of the ion-exchange materials used in solid oxide fuel cell, solid oxide electrolysis cell, gas-sensing, battery, emission control, water splitting, water purification, and other electrochemically-active devices. For the second approach, a novel silver-nickel brazing method was developed. It was demonstrated that transient porous nickel interlayers, instead of reactive element additions, could be used to promote Ag wetting on yttria stabilized zirconia (YSZ) and produce high-quality YSZ stainless steel (SS) braze joints. Mechanical tests on these reactive-element-free, silver-based SOFC braze joints, both before and after 500 hours of 750°C oxidation in air, showed that the braze and braze interface strength were higher than the underlying YSZ|NiO-YSZ substrate. The microstructural and compositional evolution of SS|Ag-Ni|YSZ braze joints exposed to dual atmosphere (air on one side and 4%H 2 -96%N 2 on the other side) for 300 hours of isothermal 750 °C aging and 300 25°C/min 35-830 °C rapid thermal cycles were compared to a Ag-3CuO braze joint. In contrast to conventional Ag-CuO RAB brazes, the Ag-Ni brazes remained pore free and well bonded, suggesting that Ag-Ni brazes may be more suitable for long term SOFC operation than Ag-CuO brazes. Other applications inspired by the silver-nickel wetting process on various ceramic substrates were also explored. In support of efforts to identify new CuO-free Ag-based brazes, density functional theory and ab initio molecular dynamics calculations were performed to 1) understand silver wetting on YSZ, and 2) identify new oxides to promote silver wetting on YSZ. These simulations found that while the formation of dissolved oxygen clusters within molten silver at the YSZ interface promoted wetting, much greater wetting angle reductions came from the formation of metal oxide interlayers at the Ag-YSZ interface. Through the development of a simple descriptor, many simple metal oxides (single cation) were examined. Unfortunately, their ability to promote Ag wetting were less than that of CuO. However, expanding the search to multi-cation oxides led to several promising candidates, such as CuAlO 2 , CuGaO 2 , and Cu 3 TiO 4 ; all of which are also stable in the reducing SOFC conditions. Depending upon their solubility in molten Ag, these newly-identified oxides could either be pre-applied as wetting promoting interlayers or directly incorporated into Ag to form new reactive air brazes.

36 MATERIALS SCIENCE↗

Hydrogen uptake in graphite matrix at high temperature

Tritium management is a critical challenge for the next generation of nuclear reactors, such as Fluoride Salt Cooled High Temperature Reactors (FHRs) and High Temperature Gas-cooled Reactors (HTGRs), due to the higher production rate (up to 10,000 times) than conventional Light Water Reactors (LWRs). Graphitic materials employed as moderator, reflector, and fuel pebbles offer a potential pathway for tritium recovery by serving as a sink for tritium. Prediction of uptake capacity under reactor relevant conditions remains a challenge due to a lack of low partial pressure data and significant inter-grade variability of graphite. This study addresses these gaps by providing a comprehensive characterization of hydrogen (as a tritium surrogate) uptake and release behavior in the A3-3 graphite matrix (GM) used in fuel pebbles. Uptake measurements are performed at reactor relevant temperatures of 600- 800 °C, 1-200 Torr hydrogen pressure, and 15-120 min equilibration time, followed by thermal desorption spectroscopy up to 1100 °C. Uptake experiments at different equilibration times demonstrate the role of kinetics in hydrogen uptake, which can be modeled as a diffusion-with-trapping process. In the thermodynamics limit, the Sips adsorption model is shown to capture the uptake in A3-3 GM well. Our campaign provides a set of new results for hydrogen uptake in A3-3, including limiting uptake capacity at 600 °C, apparent diffusion coefficient at 600 °C, and the first estimates of the FHR/HTGR relevant (600 °C, 20 Pa partial pressure) equilibrium uptake capacity and time to saturation. Desorption data highlights a new site for hydrogen uptake, not observed in nuclear graphite, which we attribute to the non-graphitized binder. Using the Kissinger method, we estimate activation energy for release from the desorption peaks, confirming the activation energy for release from the basal planes and providing the first estimate for the activation energy of release from the binder.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Vertical eddy diffusion coefficient from the LANDSAT imagery

Analysis of five stable cases of the smoke plumes that originated in eastern Cabo Frio (22 deg 59'S; 42 deg 02'W), Brazil using LANDSAT imagery is presented for different months and years. From these images the lateral standard deviation (sigma sub y) and the lateral eddy diffusion coefficient (K sub y) are obtained from the formula based on Taylor's theory of diffusion by continuous moment. The rate of kinetic energy dissipation (e) is evaluated from the diffusion parameters sigma sub y and K sub y. Then, the vertical diffusion coefficient (K sub z) is estimated using Weinstock's formulation. These results agree well with the previous experimental values obtained over water surfaces by various workers. Values of e and K sub z show the weaker mixing processes in the marine stable boundary layer. The data sample is apparently to small to include representative active turbulent regions because such regions are so intermittent in time and in space. These results form a data base for use in the development and validation of mesoscale atmospheric diffusion models.

Viswanadham, Y.↗

Application of Radon Transform to Multi-Angle Measurements Made by the Research Scanning Polarimeter: A New Approach to Cloud Tomography. Part I: Theory and Tests on Simulated Data.

The Research Scanning Polarimeter (RSP) is an airborne along-track scanner measuring the polarized and total reflectances in 9 spectral channels. The RSP was a prototype for the Aerosol Polarimetry Sensor (APS) launched on-board the NASA Glory satellite. Currently the retrieval algorithms developed for the RSP are being adopted for the measurements of the space-borne polarimeters on the upcoming NASA’s Plankton, Aerosol, Cloud Ocean Ecosystem (PACE)satellite mission. The RSP’s uniquely high angular resolution coupled with the high frequency of measurements allows for characterization of liquid water cloud droplet sizes using the polarized rainbow structure. It also provides geometric constraints on the cumulus cloud’s 2D cross section yielding the cloud’s geometric shape estimates. In this study we further build on the latter technique to develop a new tomographic approach to retrieval of cloud internal structure from remote sensing measurements. While tomography in the strict definition is a technique based on active measurements yielding a tomogram (directional optical thickness as a function of angle and offset of the view ray), we developed a “semi-tomographic” approach in which tomogram of the cloud is estimated from passive observations instead of being measured directly. This tomogram is then converted into 2D spatial distribution of the extinction coefficient using inverse Radon transform (filtered back projection) which is the standard tomographic procedure used e.g., in medical CT scans. This algorithm is computationally inexpensive compared to techniques relying on highly-multi-dimensional least-square fitting; it does not require iterative 3D RT simulations. The resulting extinction distribution is defined up to an unknown constant factor, so we discuss the ways to calibrate it using additional independent measurements. In the next step we use the profile of the droplet size distribution parameters from the cloud’s side (derived by fitting the polarized rainbows) to convert the 2D extinction distribution into that of the droplet number concentration. We illustrate and validate the proposed technique using 3D-RT-simulatedRSP observations of a LES-generated Cu cloud. Quantitative comparisons between the retrieved and the original optical and microphysical parameters are presented.

clouds↗

Effective Field Theory Analysis of CDMSlite Run 2 Data

CDMSlite Run 2 was a search for weakly interacting massive particles (WIMPs) with a cryogenic 600 g Ge detector operated in a high-voltage mode to optimize sensitivity to WIMPs of relatively low mass from 2 - 20 GeV/$c^2$. In this article, we present an effective field theory (EFT) analysis of the CDMSlite Run 2 data using an extended energy range and a comprehensive treatment of the expected background. A binned likelihood Bayesian analysis was performed on the recoil energy data, taking into account the parameters of the EFT interactions and optimizing the data selection with respect to the dominant background components. Energy regions within 5$\sigma$ of known activation peaks were removed from the analysis. The Bayesian evidences resulting from the different operator hypotheses show that the CDMSlite Run 2 data are consistent with the background-only models and do not allow for a signal interpretation assuming any additional EFT interaction. Consequently, upper limits on the WIMP mass and coupling-coefficient amplitudes and phases are presented for each EFT operator. These limits improve previous CDMSlite Run 2 bounds for WIMP masses above 5 GeV/$c^2$.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Development of a water source heat pump hardware-in-the-loop (HIL) testing facility for smart building applications

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while maintaining the thermal comfort of building occupants. As the pivotalpoint between building electricity consumption and indoor thermal comfort, high-efficiency electrical heatpumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on the heatpump operation andperformance under different weather, building load, and grid requests calls for investigation and verification via experimental tests with actual heat pumps integrated with real-time building and grid responses. This study presents a Water-Source Heat Pump Hardware-in-The-Loop (HIL) Test Facility developed with the capability to perform such tests. The hardware configuration for this testfacility introduces a hydronic system that emulates the conditions for the heat pump water-side, and a duct system that emulates conditions for the heat pump airside. Both data acquisition and emulator control are implemented through the National Instruments (NI) LabVIEW software running on an NI PXIplatform. The HIL mechanism based on the hardware-software integration that allows the testbed to communicate with a generic simulation environment is also discussed. Currently, the test facility setup includes a single heatpump and virtual building model in EnergyPlus coupled with an occupant behavioral model in MATLAB. Preliminary test results of the current setup demonstrate the building load emulator's ability to track the simulated gone temperature with a Root Mean Square Deviation (RSME) below 0.12°C (0.216°F). An uncertainty analysis based on sensor accuracies shows that the heat pump coefficient of performance (COP) can be measured with a relative uncertainty of 10.4% in cooling and 3.7% in heating. Apartfrom the current testing on a single heat pump, the test facility also provides the flexibility to include additional heat pumps to form a heat pump cluster, as well as coupling the heat pump with active thermal storage to provide enhanced demandflexibility.

Calfa, Caleb↗

Is more always better? Tribofilm evolution and tribological behavior impacted by the concentration of ZDDP, ionic liquid, and ZDDP-Ionic liquid combination

Anti-wear additives (AWs), like zinc dialkyldithiophosphates (ZDDPs), have been used in lubricants for more than a half-century, and more effective candidate AWs, like oil-soluble ionic liquids (ILs), are emerging. However, there is a lack of fundamental understanding of how the AW concentration impacts the tribochemical activities at the contact interface. This study systematically investigated the friction and wear behavior in correlation to the tribofilm morphology and composition as a function of the concentration of a ZDDP, a phosphonium-alkylphosphate IL, and an IL+ZDDP combination. The ZDDP concentration rising from 0.4 to 6.4 wt% caused a proportionally increased friction coefficient, which was well correlated to a thicker, rougher, and more brittle tribofilm. This was further understood by the decreasing phosphates but increasing sulfur compounds in the ZDDP tribofilm. In contrast, the IL and IL+ZDDP maintained a low and stable friction coefficient at 0.5 wt% or above. There was no sulfur in the IL tribofilm, and the sulfur compounds were maintained at a low level in the tribofilm for IL+ZDDP even at high concentrations. Results suggested minimizing the sulfur content in the tribofilm would be an effective way to control the friction. An interesting V-shape relationship between the AW concentration and the wear loss was observed for both the IL-containing AWs, with an optimum at 2 wt% for the IL alone and 0.46–0.92 wt% for the IL+ZDDP, respectively. A unique heat map was created using the phosphorus binding energy distribution to reveal the phosphate polymerization at different depths of the tribofilms, which provided additional insight for the tribofilm evolution.

36 MATERIALS SCIENCE↗

In Situ High-Temperature Raman Spectroscopy of UCl 3 : A Combined Experimental and Theoretical Study

Uranium trichloride (UCl 3 ) has received growing interest for its use in uranium-fueled molten salt reactors and in the pyrochemical processing of used fuel. In this paper, we report for the first time the experimentally determined Raman spectra of UCl 3 , at both ambient condition and in situ high temperatures up to 871 K. The frequencies of five of the Raman-active vibrational modes (v i ) of UCl 3 exhibit a negative temperature derivative ((∂ν i /∂T) P ) with increasing temperature. Here, this red-shift behavior is likely due to the elongation of U–Cl bonds. The average isobaric mode Grüneisen parameter (γ iP = 0.91 ± 0.02) of UCl 3 was determined through use of the coefficient of thermal expansion published in Vogel et al. (2021) and the (∂ νi /∂T) P values determined in this study. These results are in general agreement with those calculated here by density functional theory (DFT+U). Finally, a comparison of the ambient band positions of UCl 3 to those of isostructural lanthanide (La–Eu) and actinide chlorides (Am–Cf) has been made.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Two-dimensional ring polymer molecular dynamics determination of the MnO + + H 2 /D 2 reaction rates on a Δ-machine learned potential energy surface

In this work, we investigate the impact of nuclear quantum effects in the kinetics of the MnO + + H 2 reaction, a prototypical system for gas-phase H 2 activation by transition metal oxide ions. The DFT based potential energy surfaces (PESs) for the lowest-lying quintet and septet spin states reported in our previous work are improved by 2953 newly calculated CCSD(T)/AVDZ points using a delta-machine learning (Δ-ML) method. To examine nuclear quantum effects, the rate coefficients are computed using ring-polymer molecular dynamics (RPMD). Due to the mechanistic complexity of the reaction, two reaction coordinates are necessary to map out the free-energy surface and an extended RPMD rate theory is developed. The calculated RPMD rate coefficients on the new PES are in better agreement with experimental data. The calculated kinetic isotope effects (KIEs) range from 1.6 to 1.8, also in good agreement with the experimental values. Additionally, the converged RPMD rate coefficients are 2.3 to 3.1 times higher than their classical counterparts across the studied temperature range, suggesting the presence of moderate nuclear quantum effects in this reaction.

Liu, Yang [University of New Mexico, Albuquerque, ↗

Improved Modeling of Waste Isolation Pilot Plant Disposal Room Porosity

Disposal rooms at the Waste Isolation Pilot Plant (WIPP) contain waste and gas, and their porosity evolves over time. This report presents several improvements to the disposal room porosity model and presents new porosity predictions for use in future WIPP Performance Assessment activities. The improvements pertain to three sub-models: the geomechanical model, the waste compaction model, and the gas pressurization model. The impacts of each major improvement were quantified and the new porosity predictions were shown to be both mesh and domain size converged. Also, sensitivity studies on the disposal room horizon, clay seam friction coefficients, and homogenized waste representation were performed to support assumptions in the disposal room porosity model. To compare the legacy and new porosity predictions, the simulation results were plotted as a response surface, where gas pressure and time are inputs and porosity is the output. The new porosity response surface is insensitive to pressures beneath lithostatic pressure and highly sensitive to pressures above lithostatic pressure. The legacy porosity response surface, on the other hand, has moderate porosity gradients over all pressures. The new porosity response surface has a stronger scientific foundation than the legacy surface and may now be used for Compliance Decision Analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Towards the Development of Molten Salt Reactor Multiphysics Modeling Capabilities: Depletion, Thermochemistry, Reactor Dynamics, and Species Tracking

This report documents the simulation and modeling work performed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing the tools and codes developed under the program. The focus in this Fiscal Year (FY) is on the depletion, thermochemistry, reactor dynamics, and species transport in MSRs. The fission products travelling in the entire primary loop is unique in MSRs. There are many physical/chemical mechanisms that can significantly affect the fission products concentration in the fuel salt, e.g., diffusion into gas bubbles, diffusion in the core moderators, deposition on the surface of structures/pipes, etc. The fission product concentration in the fuel salt can potentially affect the neutronics, system dynamics, chemical potentials, etc., so accurately predicting the fission products distribution in the primary loop is important. Under the main objective, different activities are pursued to assess the NEAMS tools applicability for a variety of areas on species transport, including the depletion analyses using Griffin (fission product source), the redox potential assessment at different depletion points, the assessment of xenon transport effect on the neutronics and reactor dynamics, a survey of correlations for calculating the mass transfer coefficient of different fission products, etc.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗