Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation methods and programs”

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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

GeoMelt{sup R} In-Container Vitrification (ICV){sup TM} for Fukushima Daiichi Water Treatment Secondary Wastes - 20212

The Japanese government is supporting development work implemented jointly by Veolia subsidiaries Kurion Japan, K.K., Veolia Nuclear Solutions (VNS), Inc., and Veolia Nuclear Solutions Federal Services, LLC for treating radioactive waste generated from Fukushima Daiichi Nuclear Power Station (NPS) water treatment using the GeoMelt{sup TM} In-Container Vitrification (ICV){sup TM} technology. The initial work consisted of glass formulation and engineering-scale testing which was completed in 2018, in the frame of an IRID (International Research Institute for Nuclear Decommissioning) program as part of a project subsidized by Japan's Ministry of Economy, Trade and Industry (METI). The Fukushima Daiichi NPS Mid- and Long-Term Road-map requires investigation of methods to stabilize solid wastes (and to immobilize radioisotopes in the wastes) generated as a result of emergency response and decommissioning activities. Cooling water treatment has resulted in a significant amount of solid and slurry secondary wastes (mostly adsorbents and ion-exchange materials) which will require processing at some point. GeoMelt{sup R} ICV{sup TM} is a joule-heated melter technology which uses a refractory-lined single-use container combining the melter and disposal container. There is no pouring required nor concerns with refractory corrosion which allows the process to accommodate a wide range of waste chemistries and high waste loadings. The testing described here consisted of three engineering-scale melts, each processing between 212 kg and 240 kg of waste simulants, glass formers, and non-radioactive cesium (Cs) and strontium (Sr) tracers. Continuous isokinetic stack sampling of off-gas emissions was performed for each test in order to calculate Cs and Sr retention in the glass wasteform. Single-pass retention of Cs in the final glass wasteform ranged from 91.46 to 99.30%, and single-pass retention of Sr ranged from 99.76 to 100%. Planned particulate recycle will increase these retention levels. Melt 1 processed a mixture of KUR-EH (a zeolite-based ion-exchange material), simulated Advanced Liquid Processing System (ALPS) Carbonate and Iron Slurries, and glass additives. Melt 2 processed a mixture of KUR-EH, KUR-TSG (a titanate-based adsorbent), and glass additives. Melt 3 processed a mixture of KUR-EH, simulated barium sulfate/iron ferrocyanide sludge (AREVA sludge), and glass additives. Waste loadings for these melts ranged from 70 weight percent (wt%) to 82 wt%. Vitrification produces a waste form much denser than the stored water treatment secondary waste wastes, resulting in significant volume reduction. Volume reductions for the three tests ranged from 74 to 79 vol%. Vitrification produces a chemically durable wasteform. Pacific Northwest National Laboratory (PNNL) tested three glass samples from each engineering-scale melt) by the Materials Characterization Center 1 (MCC-1) test, an international standard leach test of the chemical durability of nuclear waste glasses. PNNL also obtained one U.S. reference glass (EA Glass) and two Japan reference glasses P0798) and tested these under the same MCC-1 conditions (90 deg. C, 10 m-1, DIW, and 7, 14, 28-day) as the GeoMelt{sup R} ICV{sup TM} glasses. The GeoMelt{sup R} ICV{sup TM} glasses exhibited lower total normalized releases and 14- to 28-day normalized release rates than the three reference glasses. These results suggest that the GeoMelt{sup R} ICV{sup TM} glasses have durabilities on par with high-level waste glasses under standard test conditions. Post-melt process sampling and analysis indicated no Cs migration into the melter refractory materials and very little deposition of Cs or Sr onto the melter hood or off-gas piping internals. The results of the testing indicated good Cs retention in the glass, high volume reduction and waste loadings, and excellent chemical durability. These factors are important to minimize treatment costs and to protect workers and the environment. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

LANL Critical Benchmark Comparison Study and Subsequent Revision

As part of an international collaboration within the DOE Nuclear Criticality Safety Program (NCSP), LANL is involved in a comparison study to quantify differences in k-effective results from neutron transport simulations of critical benchmark experiments. The DOE NCSP Mission and Vision details the activity in which the French Institut De Radioprotection et De Sûreté Nucléaire (IRSN) leads the study with LANL and in conjunction with ORNL and LLNL to compare results of various neutron transport codes and nuclear data libraries to compute k-effective for ICSBEP benchmarks held in common by the entities. This report documents results obtained through partial completion of the overall effort with a focus on the changes made to LANL benchmarks modeled with MCNP6 using ENDF/B-VII.1 nuclear data that appeared to have discrepant results when compared with results of other codes. The feedback received through participation in the comparison collaboration has prompted an effort to review particular input files for benchmarks and revise when necessary. This report documents the results of review and revision of specific benchmarks highlighted as possibly discrepant in the comparison study. In addition, this effort prompted a new collaboration between LANL XCP and NCS Divisions in the development of a shared review/revision procedure and use of a new benchmark repository.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A symbolic framework to obtain mid-fidelity models of flexible multibody systems with application to horizontal-axis wind turbines

Abstract. The article presents a symbolic framework (also called computer algebra program) that is used to obtain, in symbolic mathematical form, the linear and nonlinear equations of motion of a mid-fidelity multibody system including rigid and flexible bodies. Our approach is based on Kane's method and a nonlinear shape function representation for flexible bodies. The shape function approach does not represent the state of the art for flexible multibody dynamics but is an effective trade-off to obtain mid-fidelity models with few degrees of freedom, taking advantage of the separation of space and time. The method yields compact symbolic equations of motion with implicit account of the constraints. The general and automatic framework facilitates the creation and manipulation of models with various levels of complexity by adding or removing degrees of freedom. The symbolic treatment allows for analytical gradients and linearized equations of motion. The linear and nonlinear equations can be exported to Python code or dedicated software. There are multiple applications, such as time domain simulation, stability analyses, frequency domain analyses, advanced controller design, state observers, and digital twins. In this article, we describe the method we used to systematically generate the equations of motion of multibody systems and present the implementation of the framework using the Python package SymPy. We apply the framework to generate illustrative land-based and offshore wind turbine models. We compare our results with OpenFAST simulations and discuss the advantages and limitations of the method. The Python implementation is provided as an open-source project.

Branlard, Emmanuel (ORCID:0000000277506128)↗
DOE OSTI2021

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

The Reactivity and Structural Dynamics of Supported Metal Nanoclusters Using Electron Microscopy, in-situ X-Ray Spectroscopy, Electronic Structure Theories, and Molecular Dynamics Simulations

This research will provide new methods to determine the structures and compositions of heterogeneous catalysts as they exist in reaction conditions. The new tools will make it possible to understand how extremely small material structures—ones in a size range spanning from single atoms to larger polyatomic clusters—make it possible to carry out highly selective chemical transformations required to synthesize useful materials as well as for efficient energy production and utilization. The tools will make it possible to see how the atomic structures of catalytic materials evolve over time under the conditions of their use, and in this way provide understandings to guide the design of new energy-efficient catalysts. The work will make it possible to answer long standing questions that include the following. How does catalyst structure help guide molecules along specific reaction pathways? What are the important intermediates involved in these transformations? How do the structures and reactivity of catalytic materials, as well as the chemistry of their interfaces, vary with reaction conditions? Is it possible to develop new catalytic materials using new knowledge of the active structures and reaction mechanisms? The work will provide heretofore unavailable experimental tools that can be applied in advanced programs of research addressing such questions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Optimal Renewable Resource Allocation and Load Scheduling of Resilient Communities

This paper presents a methodology for enhancing community resilience through optimal renewable resource allocation and load scheduling in order to minimize unserved load and thermal discomfort. The proposed control architecture distributes the computational effort and is easier to be scaled up than traditional centralized control. The decentralized control architecture consists of two layers: The community operator layer (COL) allocates the limited amount of renewable energy resource according to the power flexibility of each building. The building agent layer (BAL) addresses the optimal load scheduling problem for each building with the allowable load determined by the COL. Both layers are formulated as a model predictive control (MPC) based optimization. Simulation scenarios are designed to compare different combinations of building weighting methods and objective functions to provide guidance for real-world deployment by community and microgrid operators. The results indicate that the impact of power flexibility is more prominent than the weighting factor to the resource allocation process. Allocation based purely on occupancy status could lead to an increase of PV curtailment. Further, it is necessary for the building agent to have multi-objective optimization to minimize unserved load ratio and maximize comfort simultaneously.

model predictive control↗
DOE OSTI2021

Microstructure Investigation and Mechanical Properties of Coated Zircaloy Cladding

Accident Tolerant Fuel concepts have been developed and tested in diverse research programs presently around the world. Industry teams have developed proprietary coating compositions and geometries, but significant knowledge gaps exist in terms of the resulting properties of the coating and the cladding and how those components impact each other in reactor conditions. Furthermore, the national laboratory and university teams currently studying the industry concepts are limited in their ability to freely disseminate and discuss results obtained on commercial coatings. The Advanced Fuels Campaign is therefore working to develop an understanding of ‘generic’ coatings applied to commercial nuclear cladding alloys free of these constraints. This milestone report provides insights in the microstructure of a 7-micrometer Cr coating on Zircaloy-2 and a 5- micrometer coating on Zry-4 as reference materials. SEM and TEM methods were utilized to describe and identify the interface, while nanoindentation was used to determine hardness and E-Moduli of coating and cladding material. This milestone summarizes the acquired knowledge to build a necessary verified foundation for future simulation and qualification efforts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

CoSim: A Simulator for Co-Scheduling of Batch and On-Demand Jobs in HPC Datacenters

The increasing scale and complexity of scientific applications are rapidly transforming the ecosystem of tools, methods, and workflows adopted by the high-performance computing (HPC) community. Big data analytics and deep learning are gaining traction as essential components in this ecosystem in a variety of scenarios, such as, steering of experimental instruments , acceleration of high-fidelity simulations through surrogate computations, and guided ensemble searches. In this context, the batch job model traditionally adopted by the supercomputing in-frastructures needs to be complemented with support to schedule opportunistic on-demand analytics jobs, leading to the problem of efficient preemption of batch jobs with minimum loss of progress. In this paper, we design and implement a simulator, CoSim, that enables on-the-fly analysis of the trade-offs arising between delaying the start of opportunistic on-demand jobs, which leads to longer analytics latency, and loss of progress due to preemption of batch jobs, which is necessary to make room for on-demand jobs. To this end, we propose an algorithm based on dynamic programming with predictable performance and scalability that enables supercomputing infrastructure schedulers to analyze the aforementioned trade-off and take decisions in near real-time. Compared with other state-of-art approaches using traces of the Theta pre-Exascale machine, our approach is capable of finding the optimal solution, while achieving high performance and scalability.

batch job preemption↗
DOE OSTI2022

Smart Ventilation Controls for Occupancy and Auxiliary Fan Use Across U.S. Climates

Smart Ventilation has been developed as a way to reduce the energy associated with ventilation by changing when ventilation happens and how much ventilation occurs at any given time. In high performance buildings with low envelope and appliance related loads, ventilation is becoming a bigger part of the total building energy use and needs to be addressed if performance targets are to be met. This work explores the development and performance of smart ventilation controls based on occupancy and auxiliary fan operation that provide annual dwelling unit ventilation equivalence to ASHRAE Standard 62.2-2016. A prototype high performance home compliant with U.S. DOE Building America Zero Energy Ready program requirements was simulated using the REGCAP tool with and without smart controls across 15 U.S. DOE climate zones. Balanced and unbalanced IAQ fans were independently simulated, and all smart controlled fans had airflows double the reference 62.2-2016 airflow. Three idealized occupancy patterns were examined: 1st shift (a typical daily work/school absence), an extended 1st shift with more time spent out of the home evenings and weekends and 3rd shift (night work). The Occupancy controller turns the IAQ fan off during unoccupied periods, and it resumes ventilation upon their return. While unoccupied, contaminants are allowed to accumulate in the space, because occupants are not exposed to these contaminants. When occupants return home, they are exposed to these higher contaminant concentrations, and the controller increases the ventilation rate sufficiently to ensure equivalence with a continuous IAQ fan. Accounting for pollutant emissions that occur during unoccupied periods (as required by ASHRAE 62.2-2016), sharply distinguishes our occupancy controls from past Demand Controlled Ventilation systems. This new accounting method results in equivalent contaminant exposure, as well as lower reductions in average ventilation rates and lower energy savings. Smart controls were demonstrated that saved HVAC energy (i.e., avoided heating/cooling load, as well as fan energy)—averaging between 6 and 46% of ventilation-related energy use depending on the control strategy and occupancy pattern assumptions. The greatest savings were in the combined Auxiliary Fan + Occupancy control. Energy savings increased with climate heating demand and longer unoccupied time periods. The 3rd shift occupancy pattern had better performance, due to the thermal benefit of reducing the ventilation rate during the cold nighttime hours. This same effect provided a cooling energy benefit in hotter locations for the 1st shift. Overall, savings from occupancy-based smart controls were low, because of the recovery period required after occupants return home, during which the airflow is double the 62.2 reference. This recovery is required to maintain equivalence with the ASHRAE standard. Occupancy-based control performance was improved when combined with sensing auxiliary fans and when providing a pre-occupancy flush out of one- or two-hours. Performance was similarly improved if the ASHRAE Standard were to recognize that pollutant emissions are lower during unoccupied periods, iii allowing a lower target ventilation rate during unoccupied periods (not currently in the standard) (see Full vs. Half AEQ in this report). Finally, over-sized unbalanced fans that are cycled on-and-off by a smart controller (or timer) were found to substantially increase annual average air exchange and energy use relative to a continuous unbalanced fan due to the effects of superposition with natural infiltration.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

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 corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

kynema-fmb [SWR-23-07]

Kynema-FMB (FKA: Kynema) is an open-source performance portable flexible multibody (FMB) dynamics solver designed for time-domain simulations. While originally tailored for wind turbine structural dynamics, the formulation and implementation are those of a general flexible-multidbody dynamics solver that can readily be applied to a wide range of systems. Kynema was designed with a narrow focus, namely to provide a lightweight, fast, accurate FMD solver for coupling to computational-fluid-dynamics (CFD) codes, especially the CFD codes in the Kynema suite, for fluid-structure-interaction (FSI) simulations. Kynema-FMB is equipped to model systems that can be represented as a collection of beams and rigid bodies that are connected through constraints. Degrees of freedom are defined in the inertial/global frame of reference and include displacements and rotations (formally as rotation matrices, but stored as quaternions). The underlying formulation is built on a Lie-group time integrator designed for index-3 differential-algebraic equations, which is second-order accurate in time (Bruls et al., 2012). Beam models are based on geometrically exact beam theory and are discretized as high-order spectral finite elements similar to those in BeamDyn (Wang et al., 2017). The governing equations for a FMD system like a wind turbine constitute a highly nonlinear system of constrained partial-differential equations. Kynema-FMB uses analytical Jacobians in the nonlinear-system solves in each time step. Linear systems use sparse storage and several third-party sparse-linear-system solvers are enabled. Ill conditioning of linear systems is mitigated with preconditioning described in Bottasso et al, 2008. Kynema-FMB is integrated with a simple open-source controller (ROSCO). There is an application programming interface (API) for coupling to geometry-resolved CFD (like that in Sharma et al., 2023) and actuator-force CFD (like that in Kuhn et al., 2025). In the latter, for actuator-line models, Kynema-FMB includes an internal blade-element solver that depends on user-provided lookup tables for coefficients of lift and drag, i.e., aerodynamic polars. Kynema-FMB is written in C++ and leverages Kokkos and Kokkos-Kernels (KokkosEcosystem) as its performance portability layer enabling simulations on both CPU and GPU systems. The repository is equipped with extensive automated testing at the unit and regression/system levels. The following describes the high-level development objectives conceived for Kynema: *Kynema will follow modern software development best practices, including test-driven development (TDD), version control, hierarchical automated testing, and continuous integration (CI) for a robust development environment. *The core data structures are memory efficient and enable vectorization and parallelization at multiple levels. *Data structures are data-oriented to exploit methods for accelerated computing including high utilization of chip resources (e.g., single instruction multiple data (SIMD) instruction sets) and parallelization using GP-GPUs. *The computational algorithms incorporate robust open-source libraries for mathematical operations, resource allocation, and data management. *The API design considers multiple stakeholder needs and ensure integration with existing and future ecosystems for data science, machine learning, and AI. *Kynema-FMB is written in modern C++ and leverages Kokkos as its performance-portability library with inspiration from the kynema stack.

Sprague, MichaelA.↗
DOE OSTI2022

The outer stellar mass of massive galaxies: a simple tracer of halo mass with scatter comparable to richness and reduced projection effects

ABSTRACT Using the weak gravitational lensing data from the Hyper Suprime-Cam Subaru Strategic Program (HSC survey), we study the potential of different stellar mass estimates in tracing halo mass. We consider galaxies with log10(M⋆/M⊙) > 11.5 at 0.2 < z < 0.5 with carefully measured light profiles, and clusters from the redMaPPer and CAMIRA richness-based algorithms. We devise a method (the ‘Top-N test’) to evaluate the scatter in the halo mass–observable relation for different tracers, and to inter-compare halo mass proxies in four number density bins using stacked galaxy–galaxy lensing profiles. This test reveals three key findings. Stellar masses based on CModel photometry and aperture luminosity within R <30 kpc are poor proxies of halo mass. In contrast, the stellar mass of the outer envelope is an excellent halo mass proxy. The stellar mass within R = [50, 100] kpc, M⋆, [50, 100], has performance comparable to the state-of-the-art richness-based cluster finders at log10Mvir ≳ 14.0 and could be a better halo mass tracer at lower halo masses. Finally, using N-body simulations, we find that the lensing profiles of massive haloes selected by M⋆, [50, 100] are consistent with the expectation for a sample without projection or mis-centring effects. Richness-selected clusters, on the other hand, display an excess at R ∼ 1 Mpc in their lensing profiles, which may suggest a more significant impact from selection biases. These results suggest that M⋆-based tracers have distinct advantages in identifying massive haloes, which could open up new avenues for cluster cosmology. The codes and data used in this work can be found here:

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

DISCo-microbe: design of an identifiable synthetic community of microbes

Microbiomes are extremely important for their host organisms, providing many vital functions and extending their hosts’ phenotypes. Natural studies of host-associated microbiomes can be difficult to interpret due to the high complexity of microbial communities, which hinders our ability to track and identify individual members along with the many factors that structure or perturb those communities. For this reason, researchers have turned to synthetic or constructed communities in which the identities of all members are known. However, due to the lack of tracking methods and the difficulty of creating a more diverse and identifiable community that can be distinguished through next-generation sequencing, most such in vivo studies have used only a few strains. To address this issue, we developed DISCo-microbe, a program for the design of an identifiable synthetic community of microbes for use in in vivo experimentation. The program is composed of two modules; (1) create, which allows the user to generate a highly diverse community list from an input DNA sequence alignment using a custom nucleotide distance algorithm, and (2) subsample, which subsamples the community list to either represent a number of grouping variables, including taxonomic proportions, or to reach a user-specified maximum number of community members. As an example, we demonstrate the generation of a synthetic microbial community that can be distinguished through amplicon sequencing. The synthetic microbial community in this example consisted of 2,122 members from a starting DNA sequence alignment of 10,000 16S rRNA sequences from the Ribosomal Database Project. We generated simulated Illumina sequencing data from the constructed community and demonstrate that DISCo-microbe is capable of designing diverse communities with members distinguishable by amplicon sequencing. Using the simulated data we were able to recover sequences from between 97–100% of community members using two different post-processing workflows. Furthermore, 97–99% of sequences were assigned to a community member with zero sequences being misidentified. We then subsampled the community list using taxonomic proportions to mimic a natural plant host–associated microbiome, ultimately yielding a diverse community of 784 members. DISCo-microbe can create a highly diverse community list of microbes that can be distinguished through 16S rRNA gene sequencing, and has the ability to subsample (i.e., design) the community for the desired number of members and taxonomic proportions. Although developed for bacteria, the program allows for any alignment input from any taxonomic group, making it broadly applicable. The software and data are freely available from GitHub and Python Package Index (PYPI).

16S rRNA↗
DOE OSTI2024

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

MFANS 2024 - Formally Proving Characteristics of Cyber-Physical Systems

Cyber-physical systems (CPS) are engineered systems that rely on the smooth integration of computational algorithms and physical elements. This integration presents new challenges for verifying that systems will behave as expected. The goal of this presentation is to present current challenges and potential solutions for the formal verification of cyber-physical systems. For cyber systems, formal methods refer to systematically rigorous mathematical techniques employed in the specification, development, analysis, and verification of both software and hardware systems. Recent advancements in computer science have yielded sophisticated tools specifically designed to address challenges associated with formal methods in complex systems. These tools leverage various foundational concepts such as logic, formal languages, program semantics, type systems, type theory, and automata theory. A notable achievement in the application of formal methods is the seL4 microkernel, claimed to be the first general-purpose operating-system kernel to be verified. Its proof implies the absence of bugs and guarantees that the kernel meets specifications. For physical systems, dynamic and control theory has a history of using rigorous analytic techniques to prove functional correctness. Lyapunov, optimal, classical, modern, and robust control theories all provide rigorous mathematical methods both to analyze system performance and to design controller that can be guaranteed to meet certain objectives. Recent computational techniques like level set theory and reachability analysis provide assertions that a system's state will avoid unsafe regions. Even though success has been independently achieved for cyber systems and physical systems, the integration of such systems creates new challenges. In particular, there is an obvious discrepancy between finite-state machines and infinite-state systems, resulting in different approaches for modeling and analyzing these system. While it is possible to simulate hybrid systems, this provides only a demonstration of a performance and not proof. For hybrid systems, current formal methods and system analysis approaches typically require a workarounds to work on hybrid systems like CPS. This paper will outline the state of the art and limits of current practice for formally verifying CPS and will identify possible research directions that require attention.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Modeling Assessment of Tidal Energy Extraction in the Western Passage

Numerical models have been widely used for the resource characterization and assessment of tidal instream energy. The accurate assessment of tidal stream energy resources at a feasibility or project-design scale requires detailed hydrodynamic model simulations or high-quality field measurements. This study applied a three-dimensional finite-volume community ocean model (FVCOM) to simulate the tidal hydrodynamics in the Passamaquoddy–Cobscook Bay archipelago, with a focus on the Western Passage, to assist tidal energy resource assessment. IEC Technical specifications were considered in the model configurations and simulations. The model was calibrated and validated with field measurements. Energy fluxes and power densities along selected cross sections were calculated to evaluate the feasibility of the tidal energy development at several hotspots that feature strong currents. When taking both the high current speed and water depth into account, the model results showed that the Western Passage has great potential for the deployment of tidal energy farms. The maximum extractable power in the Western Passage was estimated using the Garrett and Cummins method. Different criteria and methods recommended by the IEC for resource characterization were evaluated and discussed using a sensitivity analysis of energy extraction for a hypothetical tidal turbine farm in the Western Passage.

16 TIDAL AND WAVE POWER↗
DOE OSTI2025

Developing a Nuclear Quality Assurance Compliant Design Methodology for Neutronic Analysis of Xe-100 Design

The primary objective of this work is to develop a design methodology compliant with nuclear quality assurance standards for the Xe-100 neutronic design verification studies. To achieve this, a Monte Carlo model of the Xe-100 reactor was constructed using the exclusion principle, transformation technique, and universe-based level specification following Idaho National Laboratory (INL) NQA level-1 compliant standards and an NQA-1 compliant version of MCNP6. The model encompasses the entire reactor core structures, including the upper plenum, core region, and lower plenum sections, along with all sub-components. The active core section was represented using the spectral regions, each comprising a particular fuel composition and temperature averaged over the considered zone, calculated by X-energy using Very Superior Old Programs (VSOP). Additionally, a component-wise temperature map was implemented into the model, not only for the core region but also for the structural components. Temperature-dependent cross-section libraries, along with thermal scattering law libraries, generated using INL NQA-1 compliant version of NJOY21, were utilized for each isotope in the burnt fuel and the structural materials. Furthermore, the volume of each modeled component was estimated using a stochastic approach with the ray tracing method in MCNP and criticality calculations were performed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗