SEARCH · Engineering Papers

Results for “complex fluids”

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.

SourceResearch connectionFrom yearThrough year

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

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Demonstration of Wake Steering Through Yaw Control in a Wind Plant Field Experiment: Cooperative Research and Development Final Report, CRADA Number CRD-16-00629

Over the last few decades, wind energy has evolved into a large international industry involving major players in the manufacturing, construction, and utility sectors. Coinciding with the industry’s growth, significant innovation in the technology has resulted in larger turbines with lower associated costs of energy and more complex designs in all subsystems. However, as the deployment of the technology has grown and its role within the electricity sector become more prominent, so have the expectations of the technology in terms of performance, reliability, and cost. The industry currently partitions its efforts into separate paths for turbine design, plant design and development, finance, grid interaction and operation, mitigation of adverse community and environmental impacts, and other areas. One prominent area where this partition is evident is in wind turbine control. Traditionally, each wind turbine in a wind plant has been controlled separately – via its own internal controller using only its own sensors. However, wind turbines in a plant interact with each other through the plant-level fluid dynamics. Wake losses (due to upstream turbines extracting energy from the winds and “waking” downstream turbines) can be up to 10% or even 20% of the gross energy production (if each turbine experienced the free stream wind inflow to the plant). A series of studies and experiments have demonstrated that there is potential for improving energy output at existing plants through plant control methods which seek to optimize total wind plant energy production over the current “greedy” approach where each turbine maximizes its own production. Wake steering induced by yaw offsets (turning the turbine to be out of the plane perpendicular to wind inflow) for upstream turbines has shown significant promise in simulations and wind tunnel experiments. In simulation studies, annual energy production has been shown to increase by 2% or more depending on the particular aspects of the wind plant (turbine spacing, meteorological conditions, etc). This project seeks to demonstrate the potential of plant-level controls via wake steering at a commercial wind plant. This is an important step towards commercialization and industry adoption of this plant-level modeling and analysis capability.

17 WIND ENERGY↗

DOE OSTI2022

Interface Capturing Simulations on Pool Boiling Performance With Multiple Nucleation Sites

High-fidelity pool boiling simulations can provide a numerical database for improving mechanistic boiling models by allowing for specific evaluation of interactions among bubbles. Previously published pool boiling simulations investigated two nucleation sites in which bubble growth at one site suppressed nucleation at the other site. Based on previous study results, more complicated interface-capturing simulations on pool boiling were conducted using PHASTA code with locally refined unstructured mesh.First, different boundary conditions (BCs) were assessed to support robustness and reproducibility of the boiling model. Then, a scale study was conducted at a larger domain with nine nucleation sites, four of which were activated. Involving more nucleation sites increased the complexity of bubble interactions from surrounding sites. Finally, bubble departure behavior influenced by wall heat flux was investigated. When heat flux was increased, the order of bubble departure changed, but diagonal bubbles always departed after one another. The departure time interval between the first and second bubble reduced as heat flux increased. The corresponding frequency was almost linearly proportional to the heat flux. In addition, bubble departure behavior was found to be greatly influenced by the nucleation site pattern. Multiple nucleation sites resulted in superimposed inhibitive effects from surrounding sites to each bubble, which extensively delayed the departure.The work presented here provides new insight on the fundamental understanding of boiling phenomena, contributes to the development of a 3D multiphase computational fluid dynamics (M-CFD) model, and provides a more comprehensive database for data-driven pool boiling studies.

Fan, Joy↗

DOE OSTI2025

Gaussian integral method for void fraction

Here, a novel method, the Gaussian Integral Method (GIM), is presented for calculating void fractions in Computational Fluid Dynamics–Discrete Element Method (CFD-DEM) simulations. GIM is versatile and applicable to various grid types, including structured and unstructured polyhedral meshes, without requiring special boundary treatments. An optimization technique is introduced to make GIM independent of grid resolution and type. The method is validated against experimental data from a fluidized bed, demonstrating that GIM produces realistic simulations closely resembling experimental observations. Additionally, unstructured polyhedral grids using GIM outperform structured grids of equivalent resolution, yielding results more aligned with experimental data. The gradient of the void fraction is computed in the CFD solver and utilized in the DEM solver for precise estimation at particle locations. Overall, GIM provides an effective solution for void fraction calculations in particulate media simulations with complex geometries, enhancing the accuracy and applicability of CFD-DEM simulations for industrial processes.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2023

Accelerating Scientific Simulations with Bi-Fidelity Weighted Transfer Learning

High-fidelity modeling is an essential design tool for many engineering applications. However, for complex systems, computational cost can be a limiting factor. Analyzing parameter sensitivity, uncertainty quantification, and design optimization require many model evaluations. Surrogate models are often used to develop the relationship between model parameters and quantities of interest. However, in the case of complex systems, surrogate models require several degrees of freedom and, thus, a large number of data points to determine the correct dependencies. For many applications, this may be prohibitively expensive. The reduction of computational requirements can be achieved by leveraging low-fidelity models. Low-fidelity models represent the system at a coarser resolution with the advantage of computational efficiency. Therefore, a bi-fidelity modeling paradigm, which augments the accuracy of a low-fidelity model in a computationally efficient manner by invoking limited runs of a high-fidelity model, can be leveraged to sufficiently balance the accuracy and computational requirements. In this work, a bi-fidelity weighted transfer learning method using neural networks was applied to a computational fluid dynamics heat transfer modeling problem. The transfer learning advantage was investigated as a function of hyperparameters. Our main finding is that the use of a bi-fidelity modeling paradigm achieves accuracy close to that of a high-fidelity Gaussian process model while significantly reducing computational cost. The bi-fidelity model achieves comparable performance with 90 high-fidelity samples-that is, 60% less than the samples needed to achieve similar accuracy without the use of bi-fidelity modeling,

Borowiec, Katarzyna↗

DOE OSTI2022

WORM (Write One, Read Many)

WORM (Write One, Run Many) is an easy to use, cross platform, embedded and extensible, functional programming language designed to facilitate the creation of input-decks for computer codes that use standard ASCII text files for input. WORM makes it easy to create generic (yet, complex and powerful) reusable models. Additionally its nature allows for complex calculations and routines to be coded once and easily reused, further simplifying the creation of input decks. WORM (Write One, Run Many) is a powerful and versatile tool designed to improve the efficiency of today’s criticality safety analyst by allowing: + input decks for parametric studies to be created quickly and easily, + calculations and variables to be imbedded into any input deck, thus allowing for meaningful parameter specifications, + problems to be specified using any combination of units, and + complex mathematically defined models to be created. A very simple syntax is employed, and therefore the WORM is easy to learn. A WORM model is essentially a standard input deck with some of its numerical values replaced by WORM code. WORM code may include and evaluate the following mathematical operators and functions: addition, subtraction, multiplication, division, exponentiation, modulus, sine, cosine, tangent, arcsine, arccosine, arctangent, the natural logarithm, logarithm base 10, integer truncation, absolute value, and random number. Several common constants, e.g., pi, e, and Avogadros’s Number (both as 6.022e23 and 0.6022), are predefined in WORM. Additionally, many unit conversion factors are also predefined: millimeters, meters, inches, feet, yards, and mils to centimeters; kilograms, pounds, and ounces to grams; liters, milliliters, gallons, and fluid ounces to cubic centimeters; and angular degrees to radians. For parametric studies, WORM supports various shorthand list specifications: the explicit step size, linear interpolation, and logarithmic interpolation. The list notation sequentially assigns multiple values to a name. WORM creates an input deck for each value of the name. If multiple lists are used, WORM steps through each list individually, i.e., WORM creates input decks corresponding to each and every permutation of the list values. Additionally, a library of standard material definitions and Perl subroutines are included. Any one of these files can be incorporated into the subject model with a simple WORM read command. WORM is completely written in Perl, the Practical Extraction and Reporting Language. Perl is one of the most portable programming languages available today. As such, the WORM works on practically any computer platform.

Sartor, Raymond↗

DOE OSTI2025

Role of Diurnal Cycle of Insolation on the MJO Propagation in the Maritime Continent

The diurnal cycle of convection in the Maritime Continent (MC) has been hypothesized to act as a barrier to the eastward propagation of the Madden‐Julian oscillation (MJO). To test this hypothesis, we use a regional model with realistic MJO to simulate an event from the boreal spring of 2013 that weakened and stalled over the MC. Two simulations are conducted: one that includes the diurnal cycle of insolation (CTL), and another without it (NO_DC). The MJO in the simulations was identified and tracked using a large‐scale precipitation tracking method that distinguishes propagation and non‐propagation unlike the usual Real‐time Multivariate MJO method. In the NO_DC simulation, the absence of diurnal heating reduces land precipitation, allowing more continuous eastward MJO propagation. An analysis of moist static energy budget reveals that MJO maintenance in NO_DC is due to increased longwave heating and reduced advection, whereas the persistent MJO propagation in NO_DC is due to increased advection and reduced longwave heating and surface latent heat flux. These processes, however, may vary across different parts of the MC, emphasizing the complexity of MJO propagation across the MC.

Zhou, Xin [National Center for Atmospheric Researc↗

DOE OSTI2021

Volatile Element Evidence of Local MORB Mantle Heterogeneity Beneath the Southwest Indian Ridge, 48°–51°E

Abstract The mantle source beneath the Southwest Indian Ridge (SWIR) reflects a complex history of contamination. Magmatic volatile contents are vital tracers of these kinds of heterogeneities, which may fractionate otherwise constant volatile/non‐volatile elemental ratios, such as the H 2 O/Ce ratio. Although several studies have recently used trace element and isotopic data to address mantle source heterogeneity and magmatic processes at SWIR 48°–51°E region, volatile element constraints provide a valuable test of models for the origins of mantle heterogeneities in this region. Here, we present new data for nine rare basaltic glass samples from the 48°–51°E region, which enable careful assessment of the effects of primary versus secondary processes on the glass volatile contents. These samples are strongly affected by variable extents of carbon degassing, and shallow assimilation of Cl‐rich fluid, but also reveal consistently high H 2 O/Ce ratios (458.8 ± 14.9), among the highest in MORBs, that cannot be explained by late‐stage secondary processes, crustal assimilation, or simple melting of peridotite mantle at variable depths. Instead, the high H 2 O/Ce ratios are features of the mantle source composition. The 48°–51°E region is notably more depleted in highly incompatible trace elements relative to other regions of the SWIR, although this depletion is not apparent in H 2 O, which is similarly abundant throughout the SWIR. We link the high H 2 O/Ce ratios in these glasses with other trace element characteristics diagnostic of subduction and fluid addition, suggesting that the mantle source reflects signatures of a refractory mantle residue that previously melted within a subduction zone.

Wang, Wei↗

DOE OSTI2020

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Insights into scale translation of methane transport in nanopores

Accurate prediction of flow behavior in shale matrix is critical for efficient development of shale gas reservoirs. In these systems, the majority of pores are in the nano-size range. As a result, continuum-based approaches may not be appropriate to simulate flow in such systems. Molecular dynamics (MD) simulations are capable of capturing the relevant microscale physics. Their relatively high computational expense, however, restricts MD simulations to rather small systems and domains. This limitation creates a gap between computational need of macroscale systems and capabilities of MD simulations. The lattice Boltzmann method (LBM) is a suitable candidate to bridge this gap. In this work, the multiple-relaxation-time (MRT)-LBM is used to study methane transport in nano-size pores. Adsorption effects near solid boundaries, as well as non-ideal behavior of fluids, are accounted for via incorporating appropriate force terms in LBM. In this work, parameters associated with the force terms in the equation of state are studied in detail, and a workflow is proposed to determine optimal values of these parameters for gas flow in slit pores. Specifically, we establish these parameters such that the range of density values that the model is able to simulate is maximized. We demonstrate this workflow by simulating gas flow where velocity and density profiles from MD simulations are used as reference data. Results from LBM simulations are in good agreement with MD reference data for pores that are 4 nm in width or larger. Moreover, we propose a preconditioning scheme to improve the stability of LBM in dealing with complex geometries. The robustness of this scheme is demonstrated by simulating several roughness geometries. This work motivates the use of LBM in scale translation of the physics of mass transport in more complex permeable media.

03 NATURAL GAS↗

DOE OSTI2020

Dimits shift, avalanche-like bursts, and solitary propagating structures in the two-field flux-balanced Hasegawa–Wakatani model for plasma edge turbulence

We show that the recently introduced two-field flux-balanced Hasegawa–Wakatani (BHW) model captures the key features of drift-wave turbulent transport mediated by zonal flows observed in more complete and accurate gyrokinetic simulations, such as the existence of a nonlinear upshift of the threshold for drift wave turbulence driven transport, often called the Dimits shift, as well as non-local transport with avalanche bursts and solitary propagating structures. Because of the approximations made in the BHW model, these observations are made for the particle flux instead of the heat flux more commonly studied in ion temperature gradient (ITG) driven turbulence in fluid or gyrokinetic codes. Many of these features are not seen in other Hasegawa–Wakatani models, which confirm the critical role of the electron dynamics parallel to the magnetic field lines. To address questions regarding the role of boundary conditions on the drift-wave zonal flow dynamics, we apply our model to both a channel domain geometry and the more typical doubly periodic geometry. We only observe strong soliton-like solutions in the particle flux for the channel geometry, in the vicinity of the boundaries, where strong velocity shear and density gradients are generated, which are absent in the doubly periodic simulations. Changing the aspect ratio of the simulation domain also has a significant effect. In domains which are elongated in the radial direction, more complex multi-scale dynamics takes place, with multiple zonal jets interacting with each other, and large scale avalanches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Development of a machine learning model for polyethylene pyrolysis using a detailed reaction mechanism

Waste plastics have recently received significant attention as the issue of waste generation continues to increase. Thermal conversion processes, such as pyrolysis and gasification, are attractive potential technologies for utilizing waste plastics and reducing overall waste generation. Efficient utilization of plastics requires a detailed understanding of the conversion process such as pyrolysis and gasification. However, a mechanistic understanding of these processes lead to large and complex kinetic schemes that are not suited for large-scale and long-time simulation methods. Currently, most modeling approaches for pyrolysis and gasification rely on globally lumped, simplified kinetic schemes that provide results that are classified by their product type and not individual species, which limit the level of fidelity achieved via modeling. A machine learning (ML) model has been developed for the primary reactions of high-density polyethylene (HDPE) in an attempt to increase computational efficiency while still maintaining a high level of detail and accuracy. The ML model is trained on a detailed reaction mechanism containing 42 total species and 737 chemical reactions. A DeepONet branch and trunk architecture was adopted to train the model using time-steps relevant to computational fluid dynamics simulations. The ML used physics-informed loss functions to ensure mass conservation. The surrogate model has been deployed in simple MFiX CFD simulations, single particle and an experimental drop tube reactor, and has shown promising performance compared to the original scheme.

Houston, Ross↗

DOE OSTI2021

Analysis of Completion Design Impact on Cluster Efficiency and Pressure-Based Well Communication in HFTS-2 Delaware Basin

The Hydraulic Fracturing Test Site 2 (HFTS-2) is a joint industry project in the Delaware basin to advance hydraulic fracturing understanding and improve productivity in shale reservoirs. The project integrates multi-disciplinary approaches to evaluate different completion designs, well spacing, inter-well communication, and stimulated rock volume, among other factors. This paper focused on two major areas related to hydraulic fracture performance. First, an analysis of different completion designs on cluster efficiency based on near well Distributed Acoustic Sensing (DAS). Second, an evaluation of well and completion designs on inter-well fracture driven interactions (FDIs) based on downhole pressure monitoring. Fluid/sand distribution and cluster efficiency analyses were based on near wellbore DAS data collected from two adjacent horizontal wells completed in two different landing zones in the Wolfcamp formation. These wells had different completion designs aiming to evaluate the effect of normal vs. extended stage lengths, perforation hole tapering and limited entry. Standard deviation from ideal fluid/sand distribution and waterfall plots were used to evaluate cluster efficiency for each design and stage. Inter-well FDIs analysis was conducted among the horizontal wells and a vertical monitor well. One horizontal well served as the monitor well while the other horizontal well was being treated. The vertical well was instrumented with downhole pressure and temperature gauges to aid monitoring fracture height growth. The pressure response during and after fracturing was characterized based on maximum pressure increase value and slope. Pressure response vs. FDIs trigger factors such distance, cluster efficiency and stage fluid volume were also analyzed. Based on the different completion and perforation designs tested, DAS analysis suggests that limited entry design worked best. Extreme limited entry showed the potential of high perforation erosion and reduced cluster efficiency. Here, the limited entry and tapered perforation design demonstrated potential to improve the cluster efficiency for extended stage lengths. Pressure monitoring across formation units proved to be critical to understand fracture interactions and fracture vertical growth. Pressure communication across different formation units during hydraulic fracturing operation indicate fractures grew upwards during Wolfcamp wells fracturing. However, this pressure communication dissipated over time. High intensity FDIs were recorded when the frac stages Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S018R001/2477607/urtec-2021-5289-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5289 were closer to the pressure gauge location in the monitor wells. Some of these stages that produced high intensity FDIs also had high fluid volume per cluster and low cluster efficiency. The multi-disciplinary and high-quality data collected from HFTS-2 helped to further understand why completion approaches such as limited entry and tapered perforation design are successful in improving cluster efficiency. The DAS data combined with downhole high-resolution pressure measurements also helped to quantify the effect of lower cluster efficiency data on the incidence and intensity of FDIs.

58 GEOSCIENCES↗

DOE OSTI2020

An Adaptive-Mesh-Refinement Based Computational Tool for Simulating Catalysis at Mesoscale

In this work, we present a computational tool for mesoscale applications using open-source exascale- computing compatible adaptive-mesh-refinement (AMR) library, AMReX [2]. AMReX is software library that enables development of application solvers with block-structured Cartesian AMR. Our tool has capabilities to include realistic geometry representation, chemical species transport, reactions and thermodynamics that are critical for capturing mesoscale physics. A significant achievement is the ability of our solver to automatically import electron microscopy data in the form of a stereolithography (STL) or pixelated file format (mrc, tiff) without undergoing the tedious task of unstructured mesh generation. This feature allows for rapid simulation of catalyst particles with complex morphologies using an immersed-boundary formulation. The use of AMR allows for higher resolutions at catalyst surface interfaces, which in turn provides an accurate description of surface reactions and transport. Our solver uses a hybrid distributed and shared memory parallelism (OpenMP/GPU-based) with which strong scaling up to 10,000 processors for realistic catalyst particle simulations have been demonstrated.

BIOMASS FUELS,MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Fracture Characterization Via AI‐Assisted Analysis of Temperature Logs

Abstract Fractures control fluid flow, mass transport, and heat transfer in a geothermal reservoir. This makes accurate characterization of fracture networks a prerequisite for optimal design and control of a reservoir's exploitation. We develop a deep‐learning procedure to identify fracture locations via interpretation of temporally and spatially continuous downhole temperature measurements. A long short‐term memory fully convolutional network (LSTM‐FCN) is used both to capture long‐term dependencies in sequential temperature data and to distill local features around fractures. A wellbore and fractured‐reservoir thermal model is established to generate temperature data for network training. The trained LSTM‐FCN exhibits a unique ability to detect multiple fractures intersecting a borehole. We use the LSTM‐FCN algorithm to evaluate the effectiveness of different‐stage wellbore temperature measurements on fracture detection in a complex fractured system. Our experiments reveal that the use of various‐stage temperature information as an input feature set improves the robustness of fracture detection to noise interference. This study indicates the practical feasibility of obtaining accurate fracture‐network reconstructions from temperature signals, at reasonable computational cost.

Yang, Xiaoyu↗

DOE OSTI2023

Extremely rapid, Lagrangian modeling of 2D flooding: A rivulet-based approach

Estimates of potential flood inundation areas and depths are critical to informing the preparedness, response, and investment decisions of many government agencies and private sector organizations, especially under a changing climate. The standard modeling approaches, however, are often either computationally intensive or constrained in their accuracy or applicability. A novel, rivulet-based, 2D model of pluvial flooding is described in this article that is 10,000 to 10 million times less computationally complex than the full solution of the shallow water equations, yet achieves inundation area hit rates of between 0.8 and 0.9 and relative absolute mean errors of 10%-20% across a wide range of flow depths. This combination of accuracy and efficiency will enable real-time depth estimates during flood events, detailed sensitivity analyses, and the generation of large ensembles to support broad uncertainty analyses.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Dominant balance-based adaptive mesh refinement for incompressible fluid flows

This work introduces a novel adaptive mesh refinement (AMR) method that utilizes dominant balance analysis (DBA) for efficient and accurate grid adaptation in computational fluid dynamics (CFD) simulations. The proposed method leverages a Gaussian mixture model (GMM) to classify grid cells into active and passive regions based on the dominant physical interactions within the equation space. By modeling truncation error probabilistically from discretized terms, the method identifies regions of high interaction where numerical accuracy is most sensitive to resolution. Unlike traditional AMR strategies, this approach does not rely on heuristic-based sensors or user-defined thresholds, providing a fully automated and problem-independent framework for AMR. Applied to the incompressible Navier-Stokes equations for steady and unsteady flow past a cylinder, the DBA-based AMR method achieves comparable accuracy to high-resolution grids while reducing computational costs by up to 70 %. The validation highlights the method’s effectiveness in capturing complex flow features while minimizing grid cells, directing computational resources toward regions with the most critical dynamics. This modular and scalable strategy is adaptable to a wide range of applications, presenting a promising tool for efficient high-fidelity simulations in CFD and other multiphysics domains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Direct Measurement of Cl − Activity in Metal Chloride Molten Salts using a Cl 2 /Cl − Electrode

Molten metal chloride salts are promising candidates for advanced heat transfer fluids in generation IV nuclear reactors and beyond. The activity of the chloride ion in the salt has a large influence on the redox characteristics of the corresponding molten salt. There exists a knowledge gap for the direct measurement of chlorobasicity (Lewis basicity) in molten-salt mixtures. Here, the focus of this work was to develop a Cl 2 /Cl − electrode for the direct measurement of chloride activity in molten metal-chloride salt mixtures. The LiCl–KCl eutectic system was utilized as the reference melt to determine the change in chloride ion activity as increasing concentrations of MgCl 2 are added to the LiCl–KCl eutectic working electrode. The results of the measurements indicated a reduction in chloride ion activity as the concentration of magnesium chloride increased, consistent with the complexation of free chlorides by Lewis acidic magnesium cationic species. The Temkin model was used to estimate the thermodynamic properties of MgCl 4 2− complex.

Chloride↗

DOE OSTI2021

Molten Salt Reactor Signatures and Modeling Study

Molten salt reactor (MSR) technologies, either liquid fueled and cooled or only liquid cooled, pose specific, unique challenges for safeguards of the special nuclear material during the operation, fueling, and maintenance of the reactor. MSRs are one type of Generation IV technologies being invested in and considered for U.S. domestic fabrication primarily for electricity and process heat production. These designs have generated growing commercial interest for several reasons, including high (≈40%) thermal efficiency, ease of fueling, improved use of uranium fuel, potential utilization of thorium fuel, and proposed inherent safety features. U.S. companies have several planned designs that differ in fuel, cooling, and neutron energy spectrum. Driven by commercial interest and the intent of licensing MSRs, the U.S. Nuclear Regulatory Commission (NRC) has developed a vision and strategy to accommodate non-light water reactors (LWR), which include MSRs (U.S. Nuclear Regulatory Commission 2019). In addition, the Department of Energy (DOE) through the Office of Nuclear Energy (NE) Office of Advanced Reactor Technologies (ART), …sponsors research, development and deployment (RD&D) activities through its Next Generation Nuclear Plant (NGNP), Advanced Reactor Concepts (ARC), and Advanced Small Modular Reactor (aSMR) programs to promote safety, technical, economical, and environmental advancements of innovative Generation IV nuclear energy technologies. Reactor types considering the use of salts, liquid metals, or gases for coolant fall under both ARC and aSMR. Therefore, Research Design & Development is being pursued by DOE-NE through national laboratories, universities, and international and industrial collaborations. Additionally, the U.S. is a member of the Gen IV International Forum (GIF). The GIF is a cooperative, multinational organization to guide and carry out research and development needed for the GEN IV reactor systems (Forum 2018). GIF evaluated numerous reactor concepts and down-selected to the six most feasible advanced reactor technologies: gas-cooled fast reactor (GFR), lead-cooled fast reactor (LFR), MSR, supercritical watercooled reactor (SCWR), sodium-cooled fast reactor (SFR), and very high temperature reactor (VHTR). In support of the growing interest domestically and internationally, the Materials Protection, Accounting, and Control Technologies (MPACT) campaign, under the DOE-NE Fuel Cycle Technologies (FCT) program, engages in R&D activities by developing advanced instrumentation and analysis for safeguards and security of modern, advanced nuclear fuel cycle (non-LWR) facilities. Because of the historic experience in the operation of the Aircraft Reactor Experiment (ARE) and the Molten Salt Reactor Experiment (MSRE) (Robertson, MSRE Design and Operations Report Part I 1965), Oak Ridge National Laboratory (ORNL) is heavily engaged in the various R&D activities through the DOE complex related to MSRs including national technical leadership of the DOE-NE MSR campaign. This report discusses and presents the outcomes of the FY19 MPACT MSR Safeguards task. The challenges presented by MSRs for nuclear material accountancy and control (NMAC) and associated safeguards will be investigated. The objective of this research is to explore and compile the safeguards requirements and identify measurement signatures through an initial high-level MSR design and develop complementary advanced simulation and modeling capabilities. A high-level ORNL-developed MSR design called the Molten Salt Demonstration Reactor (MSDR) (Bettis, Alexander and Watts 1972) was used as the target reactor design for this research. The MSDR model incorporates technology from the MSRE and the Molten Salt Breeder Reactor (Robertson, Conceptual Design of a Single-Fluid Molten-Salt Breeder Reactor 1971). But the MSDR is a 750 MWth graphite moderated liquid fueled (low-enriched uranium) MSR compared to the MSRE’s of 7.5 MWth. The focus of this report is to discuss the evaluation of novel signatures, correlations, and indicators to understand the applicability of current safeguards instrumentation to MSRs using the modeling results from the MSDR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗