Search NASA⌕ Search

SEARCH · Search NASA

Results for “Integral Benchmarks”

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,603 records · Page 89

Biological Research and Space Health Enabled by Machine Learning to Support Deep Space Missions

A key science goal of the NASA “Moon to Mars” campaign is to understand how biology responds to the Lunar, Martian, and deep space environments in order to advance fundamental knowledge, reduce risk, and support safe, productive human space missions. Through the powerful emerging computer science approaches of artificial intelligence (AI) and machine learning (ML), a paradigm shift has begun in biomedical science and engineered astronaut health systems, to enable Earth-independence and autonomy of mission operations. We present a decadal view of AI/ML architecture to support deep space mission goals, developed in concert with leaders in the field. We describe current AI/ML methods to support 1) fundamental biology, 2) in situ analytics, 3) high performance computing hardware, 4) automated science, 5) self-driving labs, 6) remote data management, 7) integrated real-time mission biomonitoring, and 8) a Precision Space Health system. Cutting-edge AI/ML approaches that can be integrated to support these domains include active learning, explainable AI, adaptive learning, causal inference, knowledge graphs, federated learning, transfer learning, and large language models. Finally, we present results from several current ML projects that are underway in the field to address key challenges of small sample n, high feature count, heterogeneity, and sparse data. These include 1) connecting omics data to phenotypic data using an ensemble model to infer causality of spaceflight rodent liver health disruption, 2) usage of explainable ML to interrogate the muscular underpinnings of spaceflight muscle atrophy, 3) ML models analyzing and determining directed acyclic graphs of human space health risk leveraging rodent bone datasets, 4) usage of large pre-trained models connecting biomedical knowledgebases with small spaceflight datasets to understand gene-to-gene interaction networks, and 5) a suite of benchmarked open science datasets (spaceflight mouse liver; radiation DNA damage) enabling programmers to identify the best ML algorithms to answer space biological science questions.

space biology↗

Biological Research and Space Health Enabled by Machine Learning to Support Deep Space Missions

A key science goal of the NASA “Moon to Mars” campaign is to understand how biology responds to the Lunar, Martian, and deep space environments in order to advance fundamental knowledge, reduce risk, and support safe, productive human space missions. Through the powerful emerging computer science approaches of artificial intelligence (AI) and machine learning (ML), a paradigm shift has begun in biomedical science and engineered astronaut health systems, to enable Earth-independence and autonomy of mission operations. We present a decadal view of AI/ML architecture to support deep space mission goals, developed in concert with leaders in the field. We describe current AI/ML methods to support 1) fundamental biology, 2) in situ analytics, 3) high performance computing hardware, 4) automated science, 5) self-driving labs, 6) remote data management, 7) integrated real-time mission biomonitoring, and 8) a Precision Space Health system. Cutting-edge AI/ML approaches that can be integrated to support these domains include active learning, explainable AI, adaptive learning, causal inference, knowledge graphs, federated learning, transfer learning, and large language models. Finally, we present results from several current ML projects that are underway in the field to address key challenges of small sample n, high feature count, heterogeneity, and sparse data. These include 1) connecting omics data to phenotypic data using an ensemble model to infer causality of spaceflight rodent liver health disruption, 2) usage of explainable ML to interrogate the muscular underpinnings of spaceflight muscle atrophy, 3) ML models analyzing and determining directed acyclic graphs of human space health risk leveraging rodent bone datasets, 4) usage of large pre-trained models connecting biomedical knowledgebases with small spaceflight datasets to understand gene-to-gene interaction networks, and 5) a suite of benchmarked open science datasets (spaceflight mouse liver; radiation DNA damage) enabling programmers to identify the best ML algorithms to answer space biological science questions.

space biology↗

Implementation and Validation of Reference Sensitivity Profile Calculations in TRIPOLI 4 - 14247

In current neutron transport calculations, particularly on continuous-energy Monte-Carlo codes, the major source of uncertainty is coming from nuclear data. In order to improve nuclear data evaluations, integral experiments are used. Feedback is then precisely obtained by re-estimating cross-section levels and correlations according to their sensitivity. k{sub eff}-sensitivity coefficients are traditionally obtained using deterministic codes since hands-on methods were not yet implemented in continuous energy Monte Carlo codes. In this paper, an Iterated-Fission-Probability-based method newly developed in the continuous energy Monte Carlo code TRIPOLI-4 will be used to calculate sensitivity coefficient breakdown on isotopes and reactions and results will be compared to those of ERANOS deterministic code. Objective is twofold: to show the reliability of the new TRIPOLI-4 method and to identify where and why profiles are slightly different. Sensitivity profiles, forward and adjoint fluxes are thus meticulously compared using a fine 1968-group energy structure. An OECD international benchmark (Uncertainty Analysis for Criticality Safety Assessment - UACSA - Phase III) is taken as working example to enforce validation and to compare to independent results (e.g. MCNP-6). This benchmark is constituted of two bare spheres filled with homogeneous mixture of uranium and polyethylene (respectively 2% and 50% U-235 enrichments) and gives an overview of issues encountered in thermal reactor cores. Several energy-group structures (1968G, 172G, 33G) are also tested in ERANOS to identify multi-group and implicit effects. Lessons learnt are that continuous Monte Carlo methods seem to be mature and compare well with deterministic codes. Only damper the U-238 elastic - and on certain conditions the U-238 capture - sensitivity coefficient for which difference can reach 30% between TRIPOLI-4 and ERANOS. Indeed concerning U-238, implicit effects are emphasized and multi-group theory remains limited by hypothesis; even a 1968-group structure is not refined enough to properly calculate U-238 elastic scattering k{sub eff}-sensitivity. This reference 3D continuous-energy sensitivity capability in TRIPOLI4-IFP will be extensively used in the uncertainty automated tool RIB of the French Criticality-Safety package CRISTAL. (authors)

42 ENGINEERING↗

L3:RTM.MCH.P19.06: Neutronics UQ

This report summarizes some of the work performed in the last year on the nuclear data uncertainty tasks using the multipole representation. This approach was proposed to provide a simple and integrated way of evaluating the nuclear data uncertainty in a Monte Carlo simulation and evaluating their impact on multigroup cross sections. As previously described, the windowed multipole format allows for a resonance based representation of the cross sections with direct Doppler broadening. This representation also lends itself nicely to the evaluation of sensitivities of resonance parameters and the direct perturbation of resonance parameters in the simulation environment. Both of these approaches rely on the need of conversion of the covariance data (File 32 of the ENDF library) into a pole and residues form: i.e. into a multipole covariance. Section 5 details the current pathways that have been used so far to achieve this. In order to verify the conversion process of the covariance data, the sensitivies to the resonance parameters and the sampling from these parameters, an analytical benchmark of the slowing process was developed. The benchmark provides a closed-form solution for the flux in the presence of Hydrogen and multiple resonant nuclides. Section 2 details the analytical benchmark and provides the closed-form solutions for the flux and its sensitivity to resonance parameters. A verification of the proposed methodology is performed and compared with the analytic solution to demonstrate the accuracy of the benchmark and the applicability of our embedded resonance perturbation approach in a small Monte Carlo code. Papers are currently being drafted documenting the benchmark and the generation of sensitivities. As paths to best generate multipole covariance matrices are being compared, we also explored new methods to propagate nuclear data uncertainty across the Monte Carlo neutron transport calculation. Section 4 studies the use of multipole covariance matrices to perform sensitivity analysis. Section 3.2 documents the results of our investigations into a new way of propagating nuclear data uncertainty across Monte Carlo transport simulations by sampling new epistemic cross sections for each new neutron.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Autonomous Medical Care for Exploration

The goal of Autonomous Medical Care (AMC) is to ensure a healthy, well-performing crew which is a primary need for exploration. The end result of this effort will be the requirements and design for medical systems for the CEV, lunar operations, and Martian operations as well as a ground-based crew health optimization plan. Without such systems, we increase the risk of medical events occurring during a mission and we risk being unable to deal with contingencies of illness and injury, potentially threatening mission success. AMC has two major components: 1) pre-flight crew health optimization and 2) in-flight medical care. The goal of pre-flight crew health optimization is to reduce the risk of illness occurring during a mission by primary prevention and prophylactic measures. In-flight autonomous medical care is the capability to provide medical care during a mission with little or no real-time support from Earth. Crew medical officers or other crew members provide routine medical care as well as medical care to ill or injured crew members using resources available in their location. Ground support becomes telemedical consultation on-board systems/people collect relevant data for ground support to review. The AMC system provides capabilities to incorporate new procedures and training and advice as required. The on-board resources in an autonomous system should be as intelligent and integrated as is feasible, but autonomous does not mean that no human will be involved. The medical field is changing rapidly, and so a challenge is to determine which items to pursue now, which to leverage other efforts (e.g. military), and which to wait for commercial forces to mature. Given that what is used for the CEV or the Moon will likely be updated before going to Mars, a critical piece of the system design will be an architecture that provides for easy incorporation of new technologies into the system. Another challenge is to determine the level of care to provide for each mission type. The level of care refers to the amount and type of care one will render based on perceived need and ability. This is in contrast to the standard of care which is the benchmark by which that care is provided. There are certainly some devices and procedures that have unique microgravity or partial gravity requirements such that terrestrial methods will not work. For example, performing CPR on Mars cannot be done in exactly the same way as on Earth because the reduced gravity causes too large a reduction in the forces available for effective compression of the chest. Likewise, fluid behavior in microgravity may require a specialized water filtration and mixing system for the creation of intravenous fluids. This paper will outline the drivers for the design of the medical care systems, prioritization and planning techniques, key system components, and long term goals.

Johnson-Throop, Kathy A.↗

Completing the Feedback Loop: The Impact of Chlorophyll Data Assimilation on the Ocean State

In anticipation of the integration of a full biochemical model into the next generation GMAO coupled system, an intermediate solution has been implemented to estimate the penetration depth (1Kd_PAR) of ocean radiation based on the chlorophyll concentration. The chlorophyll is modeled as a tracer with sources-sinks coming from the assimilation of MODIS chlorophyll data. Two experiments were conducted with the coupled ocean-atmosphere model. In the first, climatological values of Kpar were used. In the second, retrieved daily chlorophyll concentrations were assimilated and Kd_PAR was derived according to Morel et al (2007). No other data was assimilated to isolate the effects of the time-evolving chlorophyll field. The daily MODIS Kd_PAR product was used to validate the skill of the penetration depth estimation and the MERRA-OCEAN re-analysis was used as a benchmark to study the sensitivity of the upper ocean heat content and vertical temperature distribution to the chlorophyll input. In the experiment with daily chlorophyll data assimilation, the penetration depth was estimated more accurately, especially in the tropics. As a result, the temperature bias of the model was reduced. A notably robust albeit small (2-5 percent) improvement was found across the equatorial Pacific ocean, which is a critical region for seasonal to inter-annual prediction.

MERRA-OCEAN↗

Control rod treatment for FEM based radiation transport methods

This paper presents a novel control rod treatment for finite element method (FEM) based reactor physics simulations. The new method falls into the class of homogenization methods which treat material interfaces within mesh elements by defining a representative, homogenized nuclear cross section set. In contrast to existing homogenization based cusping treatments that attempt to approximate flux-volume weighting of the cross sections, we propose to exactly integrate terms appearing in the variational form of the transport or diffusion equation over mixed material elements. This goal is achieved by defining smoothly varying cross sections in elements that straddle material boundaries. It is shown that the smoothly varying cross sections can be obtained by projecting the piecewise set of cross sections onto the set of Legendre polynomials up to a finite order. Implementation of smoothly varying cross sections is supported within the Rattlesnake radiation multiphysics code that is currently under development at the Idaho National Laboratory. Homogenization methods traditionally suffer from cusping effects that originate from overestimation of absorption cross sections in strongly neutron-absorbing materials as found in control rods. We demonstrate that the presented method is (1) highly effective in removing cusping effects and (2) provides accurate results based on the C5G7-TD benchmark, a time-dependent version of the C5G7 benchmark.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Search for a non-standard-model Higgs boson decaying to a pair of new light bosons in four-muon final states

Results are reported from a search for non-standard-model Higgs boson decays to pairs of new light bosons, each of which decays into the μ + μ - final state. The new bosons may be produced either promptly or via a decay chain. The data set corresponds to an integrated luminosity of 5.3 fb -1 of proton–proton collisions at √s = 7 TeV, recorded by the CMS experiment at the LHC in 2011. Such Higgs boson decays are predicted in several scenarios of new physics, including supersymmetric models with extended Higgs sectors or hidden valleys. Thus, the results of the search are relevant for establishing whether the new particle observed in Higgs boson searches at the LHC has the properties expected for a standard model Higgs boson. No excess of events is observed with respect to the yields expected from standard model processes. A model-independent upper limit of 0.86±0.06 fb on the product of the cross section times branching fraction times acceptance is obtained. Finally, the results, which are applicable to a broad spectrum of new physics scenarios, are compared with the predictions of two benchmark models as functions of a Higgs boson mass larger than 86 GeV/c 2 and of a new light boson mass within the range 0.25–3.55 GeV/c 2 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electron temperatures in Orion A

We measured the hydrogen, helium, and carbon 106α radio recombination lines (RRLs), three different β and γ lines, and radio continuum at 4′ east, west, and south of the peak of the H ii region NGC 1976 (Messier 42, Orion A) with the NRAO Green Bank Telescope at 5 GHz. The α/β and α/γ integrated line intensities for the offset positions are in agreement with local thermodynamic equilibrium, and so the ratio of the integrated intensity from RRL to the thermal continuum allows for a simplified calculation of the electron temperature, T{sub e}. For the offset positions, the average value is T{sub e} = 7200 ± 100 K. The best value of T{sub e} for the peak continuum position, based on previous RRL data, is 8300 ± 200 K. Thus, there is a decrease in T{sub e} with angular offset. These data are consistent with the empirical model of Wilson et al., based on RRL data. There is good agreement between the T{sub e} values from RRLs and those from optical collisionally excited lines, [O iii], and Balmer decrement data, but the T{sub e} value from optical recombination lines is significantly lower. Analysis of this prototype H ii region provides additional inputs for the determination of element abundances and a benchmark for interpretations of more distant H ii regions.

79 ASTRONOMY AND ASTROPHYSICS↗

Auxiliary-state-facilitated phase synchronization phenomena in isolated spin systems

Extending classical synchronization to the quantum domain is of great interest both from the fundamental physics point of view and with a view toward quantum technology applications. This work characterizes phase synchronization of an effective spin-1 system, which is realized by coupling three quantum states with infinite lifetime to auxiliary excited states that have a finite lifetime. Integrating out the excited states, the effective spin-1 model features coherent and incoherent effective couplings. The following are our key findings. (i) Phase synchronization can be controlled by adjusting the phases of the couplings to the excited states. (ii) Unlike in the paradigmatic spin-1 system studied in the literature, where the dissipative couplings describe decay into the limit-cycle state, the effective spin-1 model investigated in this work is governed by a competition between dissipative decay into and out of the limit-cycle state, with the dissipative decay out of the limit-cycle state playing a critical role. (iii) We identify a parameter regime where phase synchronization of the effective spin-1 system is, in the absence of coherent effective couplings, governed entirely by the effective dissipators. The effective spin-1 model is benchmarked through comparisons with master-equation calculations for the full Hilbert space. Physical insights are gained through analytical perturbation theory calculations. In conclusion, our findings, which are expected to hold for a broad class of energy-level and coupling schemes, are examined using hyperfine states of 87 Rb as an example system.

Nonlinear optics↗

Uncertainty Analysis for a Criticality Benchmark [Slides]

This presentation discusses how KRUSTY was analyzed with LANL’s MCNP neutron transport code, Monte Carlo N-Particle. In the Monte Carlo method, neutrons are born in fission events and fly around the geometry having events (e.g., scatter, capture, fission, escape) according to the natural probabilities. “Natural probabilities” are given as neutron cross sections , which are measured and evaluated nuclear data. One of the goals of the benchmark project is to improve the nuclear cross sections (and the codes that use them). In this work, the response of interest is called k eff , which describes the neutron multiplication in a near critical system. Like most quantities, k eff is not measured directly; it is only inferred. In summation, MCNP results for k eff matched the measurements extremely well. The evaluated uncertainties are asymmetric and smaller than 0.00100 (e.g., + 0.00080 / – 0.00062). The largest uncertainties are due to alignment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Speeding Up Hartree–Fock in JuliaChem with Density Fitting

In this work, the density fitting (DF) approximation is added to the restricted Hartree–Fock (RHF) implementation in the JuliaChem computational chemistry code. Utilizing a DF algorithm that uses symmetry and integral screening, a significant reduction in time to compute the Fock matrix is achieved. The symmetry and screening DF-RHF techniques were adapted to be performed on graphics processing units (GPUs), which are well suited to perform the matrix multiplications that comprise the bulk of the Fock build time in DF-RHF. The JuliaChem DF-RHF GPU algorithm employs a novel approach that automatically switches between two DF-RHF algorithms depending on the number of basis functions in the calculation. The JuliaChem GPU DF-RHF implementation demonstrates up to 2× speedup for Fock build times compared to the existing best-in-class GPU DF-RHF implementation by operating directly on screened intermediate matrices. Due to the high portability of the Julia language code, the JuliaChem CPU and GPU DF-RHF implementations could be benchmarked on a variety of CPU and GPU architectures from multiple hardware vendors.

Hayes, John J. [Ames Laboratory, and Iowa State Un↗

Search for light Higgs bosons from supersymmetric cascade decays in $\text {pp}$ collisions at $\sqrt{s}=13\,\textrm{TeV}$

A search is reported for pairs of light Higgs bosons (H 1 ) produced in supersymmetric cascade decays in final states with small missing transverse momentum. A data set of LHC pp collisions collected with the CMS detector at $\sqrt{s}$ = 13 TeV and corresponding to an integrated luminosity of 138 fb -1 is used. The search targets events where both H 1 bosons decay into $\text{b}\bar{\text{b}}$ pairs that are reconstructed as large-radius jets using substructure techniques. No evidence is found for an excess of events beyond the background expectations of the standard model (SM). Results from the search are interpreted in the next-to-minimal supersymmetric extension of the SM, where a “singlino” of small mass leads to squark and gluino cascade decays that can predominantly end in a highly Lorentz-boosted singlet-like H 1 and a singlino-like neutralino of small transverse momentum. Upper limits are set on the product of the squark or gluino pair production cross section and the square of the branching fraction of the $\text{b}\bar{\text{b}}$ in a benchmark model containing almost mass-degenerate gluinos and light-flavour squarks. Under the assumption of an SM-like H 1 → $\text{b}\bar{\text{b}}$ branching fraction, H 1 bosons with masses in the range 40–120 GeV arising from the decays of squarks or gluinos with a mass of 1200–2500 GeV are excluded at 95% confidence level.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Observational Data for Next-Generation Climate Model Evaluation: Requirements, Considerations, and Best Practices

Climate model simulations are an important source of information about our planet’s climate system and also enable informed decision-making under different future scenarios. As a new archive of results from the next generation of climate models is anticipated to become available with the Coupled Model Intercomparison Project phase 7 (CMIP7), the need to develop efficient and robust methods to evaluate models is paramount. Observations are an integral part of model evaluation, providing a means to quantify and understand the degree to which climate models can faithfully reproduce Earth system processes. Such analysis is critical for constraining climate projections, identifying areas of focus for model development, and assisting analysts in deciphering the utility of models for specific applications. Observations of Earth system come from a diversity of sources, span different space–time domains, and are produced by different communities, and each dataset features different data structures and formats, metadata standards, and its own unique uncertainties. Uncertainties in an observational dataset may stem from gaps in temporal and spatial coverage, instrumentation errors, or assumptions in retrieval and processing methods. How then does one ensure that observational data are ready for use and utilized in the most appropriate way for robust, rapid, and routine climate model evaluation? The CMIP7 Model Benchmarking Task Team with input from the broader climate modeling, model evaluation, and observational data communities present a vision and considerations for best practices toward the optimal and appropriate use of observational data to support next-generation climate model evaluation.

Climate models↗

Development of a diffusive wave shallow water model with a novel stability condition and other new features

Traditionally, one of the approaches to flood modelling is numerical simulation of the diffusive wave approximation of the shallow water equations. Improving these models in various aspects is still an open area of research. In this study, a new diffusive wave model with explicit time integration was developed which includes some novel features: (1) time steps are determined using a novel stability criterion which resulted in more dynamic time steps (i.e., broader range) compared to the conventional Courant-Friedrichs-Lewy stability condition; (2) stability constraints are reduced, considering the flow processes within surface ponds; (3) besides Manning's formula, which is the common equation for computing velocities in diffusive wave models, the free fall velocity and a new equation for wave-front velocity are employed; and (4) the influence of upstream surface ponds on downstream flow is considered. This paper introduces the enhanced diffusive wave model, the socalled Overland Flow Simulator Cellular Automata (OFS-CA), and its results for five test cases. Available analytical solutions and an experimental study were used for verification. Two other shallow water models were used for comparison and benchmarking. Overall, good agreements were observed and OFS-CA was computationally less expensive compared to the other two shallow water models.

42 ENGINEERING↗

Development of Low-Cost, High-Performance, Easy-To-Apply, Non-Flammable, Inorganic Phase Change Material (PCM) Technology (Project Final Report)

This report describes a 45-months long research program focused on the development of novel, easy-to-apply, non-flammable, and high-performance inorganic phase change materials (PCMs) for building and industrial applications. The University of Massachusetts Lowell (UML) formed a world-class team consisting of researchers form InsolCorp (only N. American manufacturer of inorganic PCM systems for building applications), and a group of industrial advisors, to develop a universal/multipurpose, simple-to-manufacture and cost-effective PCM technology. The project team expects that the results of this work will spur in the future the adoption of thermal storage materials – a key building energy saving technology as identified by DOE BTO – for a variety of building envelope applications. The main goal of this project was to demonstrate a suite of low-cost, multipurpose, and durable inorganic PCM formulations with phase transition temperatures encompassing typical building applications (between +5 o C and +55 o C). The first objective was to design, fabricate, and experimentally validate a performance of inexpensive, durable, highly efficient, non-flammable, and easy to manufacture PCMs. To allow a variety of building applications, the project team focused on formulations that exhibit repeatable phase transitions between +5 o C and +55 o C. To follow the DOE BTO cost efficiency target without compromising thermal performance, our work was based on inorganic compounds (mostly salt hydrates) and their blends, which represent a fraction of the cost of most of organic PCMs with about twice as high density as well as significantly higher thermal conductivity and phase change enthalpy. The second objective was to develop easy-to-manufacture and -install packaging/encapsulation designs that are 1) a superior barrier to current state-of-the-art macro-packaging, which significantly reduces the risk of loss of hydration water and PCM leak, and 2) optimal in enhancing the heat exchange rates with the surroundings and within the PCM core to ensure complete charging/discharging of the entire PCM within the product. Finally, the project’s intend was to scale-up the fabrication process to demonstrate installation on system-scale applications, and to validate the performance under field conditions. This work aimed at developing low-cost, high-energy storage, and reliable latent heat storage technology for building applications. This development was realized by formulating and integrating the following two technology components: 1) inorganic salt hydrate based PCMs that have high latent enthalpies and are low-cost and durable, and 2) PCM encapsulation (packaging) technology that maximizes PCM concentration and enhances heat transport characteristics in the product and with the external environment/materials. High thermal storage capacity, low cost and fire resistance are key to the building market entry for PCM technology. Therefore, the project’s focus was on salt-hydrate-based formulations which satisfy all these criteria. Packaging and/or encapsulation of PCM is a key processing step. The project team recognized that a low-cost and simple-to-manufacture salt hydrate-based PCM technology holds the best chance to be successful in the building construction market, a market which is traditionally extremely sensitive to cost and where commodity thermal insulations are the benchmark for envelope-related energy saving measures. That is why, in this project, the main intention was to minimize the production cost and maximize the product energy storage density without sacrificing the PCM performance. It was achieved through: 1. Minimizing the non-PCM components (plastics, additives, packaging/encapsulation materials, etc.) because they are significantly more expensive than salt hydrates, 2. Using highly thermally conductive and lightweight PCM carrier (packaging material) to facilitate more complete phase cycling, and 3. Optimizing the thickness and minimizing air spaces in product design (such as in pouched PCM). For this purpose, our approach was to enable an easy system design, including selection of the PCM operating temperatures, optimizing the necessary heat storage capacity (by stacking together several layers of PCM products), and if needed, a synchronized usage of PCM products of different temperatures. A specially designed, robust, highly thermally conducting and highly impermeable packaging (to retain salt hydrate water during phase transition cycles) was designed and tested to increase the overall system thermal performance and durability. All PCM products developed during this project were tested in both lab scale and in full scale field conditions. It is expected that, after further developments and commercialization, the developed PCM technologies may be also applied in space conditioning, energy storage technologies, and heat transfer applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

MOSCATO Development and Integration in Fiscal Year 2024

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy evolution for molten salt equipment. In FY24, significant developments were made to the code package, enhancing its capabilities in many aspects. The improvements and advancements can be summarized as follows: 1. Implementation of tritium transport capabilities and validation with experimental data: To enable modeling of tritium and other fission gases within MSRs, we implemented gas transport within MOSCATO via inclusion of couple mass transport equations within the salt and structural alloys. Comparisons to experimental data from literature showed good agreement with respect to tritium release rates. 2. Preliminary implementation of two-phase flow models in MOSCATO: To model tritium and other gases above their solubility limits, we implemented preliminary two-phase flow models within MOSCATO to account for bubble transport. The first model adopted was the Level-Set approach, which can handle the high void fraction regime, but with a requirement for high mesh resolution thus high computational expense. In this report, we present a verification of the Level-Set method using a simple benchmark case. We also performed a demonstration of the code as applied to an experimental case involving cover gas flow through salt in an experimental vessel. The second model adopted was the Eulerian-Eulerian dispersed flow model, which is computationally cheaper but limited to low void fraction regimes, such as bubbly flow. Validation and verification have not yet been performed for the Eulerian-Eulerian approach, but a preliminary implementation was completed. 3. Validation with static corrosion experiments: Static corrosion experimental data for stainless steel coupons within molten salts was used to further validate the corrosion model in MOSCATO. To do so, we leveraged the existing models in MOSCATO and simulated the sample mass loss and mass gain phenomena. Several ion species, including Cr 2+ , Fe 2+ and H + , were simulated in salt using the PNP solver, while Cr 0 and Fe 0 were simulated with a diffusion solver in stainless steel. The mass loss of the samples was compared with experimental data, and good agreement was achieved. These combined activities served to further expand the capabilities of MOSCATO and make it more generally applicable to the full range of phenomena that can control chemistry and corrosion in molten salt reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

First-principles prediction of the softening of the silicon shock Hugoniot curve

Here, whock compression of silicon (Si) under extremely high pressures (>100 Mbar) was investigated by using two first-principles methods of orbital-free molecular dynamics (OFMD) and path integral Monte Carlo (PIMC). While pressures from the two methods agree very well, PIMC predicts a second compression maximum because of 1s electron ionization that is absent in OFMD calculations since Thomas–Fermi-based theories lack inner shell structure. The Kohn–Sham density functional theory is used to calculate the equation of state (EOS) of warm dense silicon for low-pressure loadings (P < 100 Mbar). Combining these first-principles EOS results, the principal Hugoniot curve of silicon for pressures varying from 0.80 Mbar to above ~10 Gbar was derived. We find that silicon is ~20% or more softer than what was predicted by EOS models based on the chemical picture of matter. Existing experimental data (P ≈ 1–2 Mbar) seem to indicate this softening behavior of Si, which calls for future strong-shock experiments (P > 10 Mbar) to benchmark our results.

36 MATERIALS SCIENCE↗