Search NASA⌕ Search

SEARCH · Search NASA

Results for “batch simulation”

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 217 records · Page 12

Towards a Comprehensive Computational Simulation System for Turbomachinery

The objective of this work is to develop algorithms associated with a comprehensive computational simulation system for turbomachinery flow fields. This development is accomplished in a modular fashion. These modules includes grid generation, visualization, network, simulation, toolbox, and flow modules. An interactive grid generation module is customized to facilitate the grid generation process associated with complicated turbomachinery configurations. With its user-friendly graphical user interface, the user may interactively manipulate the default settings to obtain a quality grid within a fraction of time that is usually required for building a grid about the same geometry with a general-purpose grid generation code. Non-Uniform Rational B-Spline formulations are utilized in the algorithm to maintain geometry fidelity while redistributing grid points on the solid surfaces. Bezier curve formulation is used to allow interactive construction of inner boundaries. It is also utilized to allow interactive point distribution. Cascade surfaces are transformed from three-dimensional surfaces of revolution into two-dimensional parametric planes for easy manipulation. Such a transformation allows these manipulated plane grids to be mapped to surfaces of revolution by any generatrix definition. A sophisticated visualization module is developed to al-low visualization for both grid and flow solution, steady or unsteady. A network module is built to allow data transferring in the heterogeneous environment. A flow module is integrated into this system, using an existing turbomachinery flow code. A simulation module is developed to combine the network, flow, and visualization module to achieve near real-time flow simulation about turbomachinery geometries. A toolbox module is developed to support the overall task. A batch version of the grid generation module is developed to allow portability and has been extended to allow dynamic grid generation for pitch changing turbomachinery configurations. Various applications with different characteristics are presented to demonstrate the success of this system.

Shih, Ming-Hsin↗

Stochastic Approximation for Multi-period Simulation Optimization with Streaming Input Data

We consider a continuous-valued simulation optimization (SO) problem, where a simulator is built to optimize an expected performance measure of a real-world system while parameters of the simulator are estimated from streaming data collected periodically from the system. At each period, a new batch of data is combined with the cumulative data and the parameters are re-estimated with higher precision. The system requires the decision variable to be selected in all periods. Therefore, it is sensible for the decision-maker to update the decision variable at each period by solving a more precise SO problem with the updated parameter estimate to reduce the performance loss with respect to the target system. We define this decision-making process as the multi-period SO problem and introduce a multi-period stochastic approximation (SA) framework that generates a sequence of solutions. Two algorithms are proposed: Re-start SA (ReSA) reinitializes the stepsize sequence in each period, whereas Warm-start SA (WaSA) carefully tunes the stepsizes, taking both fewer and shorter gradient-descent steps in later periods as parameter estimates become increasingly more precise. We show that under suitable strong convexity and regularity conditions, ReSA and WaSA achieve the best possible convergence rate in expected sub-optimality either when an unbiased or a simultaneous perturbation gradient estimator is employed, while WaSA accrues significantly lower computational cost as the number of periods increases. In addition, we present the regularized ReSA, which obviates the need to know the strong convexity constant and achieves the same convergence rate at the expense of additional computation.

Computer Science↗

System time-domain simulation

Complex systems are simulated by engineers without extensive computer experience. Analyst uses free-form engineering-oriented language to input "black box" description. System Time Domain (SYSTID) Simulation Program generates appropriate algorithms and proceeds with simulation. Program is easily linked to postprocessing routines. SYSTID program is written in FORTRAN IV for batch execution and has been implemented on UNIVAC 1110 under control of EXEC 8, Level 31.

Dawson, C. T.↗

Evaluation of Heat Recuperation in a Concentric Hydrogen Reduction Reactor

Heat recuperation in an ISRU reactor system involves the recovery of heat from a reacted regolith batch by transferring this energy into a batch of fresh regolith. One concept for a hydrogen reduction reactor is a concentric chamber design where heat is transferred from the inner, reaction chamber into fresh regolith in the outer, recuperation chamber. This concept was tested and analyzed to define the overall benefit compared to a more traditional single chamber batch reactor. Data was gathered for heat-up and recuperation in the inner chamber alone, simulating a single chamber design, as well as recuperation into the outer chamber, simulating a dual chamber design. Experimental data was also used to improve two analytical models, with good agreement for temperature behavior during recuperation, calculated mass of the reactor concepts, and energy required during heat-up. The five tests, performed using JSC-1A regolith simulant, also explored the effectiveness of helium gas fluidization, hydrogen gas fluidization, and vibrational fluidization. Results indicate that higher hydrogen volumetric flow rates are required compared to helium for complete fluidization and mixing, and that vibrational fluidization may provide equivalent mixing while eliminating the need to flow large amounts of excess hydrogen. Analysis of the total energy required for heat-up and steady-state operations for a variety of conditions and assumptions shows that the dual-chamber concept requires the same or more energy than the single chamber concept. With no clear energy savings, the added mass and complexity of the dual-chamber makes it unlikely that this design concept will provide any added benefit to the overall ISRU oxygen production system.

Linne, Diane↗

Simulant Development of Potential 200 West Area Waste Feeds

Preliminary planning for retrieval, qualification, and pretreatment of waste in Hanford’s 200 West Area (200W) has begun as part of the West Area Risk Management project. Experimental studies to technically mature pretreatment process operations will likely be needed because of the uniqueness of 200W waste. Pacific Northwest National Laboratory formulated five simulants to represent 200W-qualified feed based on the preliminary flowsheet provided by Washington River Protection Solutions, LLC. The simulant recipes were devised using applicable historical information as a reference point to support the use of the flowsheet waste vectors, which were combined into five distinct groups. These five groups formed the basis for the liquid composition targets that were adapted into recipes using charged-balanced salt species. The liquid phase recipes were batched in 1-L quantities and analyzed at Pacific Northwest National Laboratory. Once confirmed to be stable, the liquid solutions were tested for compatibility with candidate solid components. Specific solid components were recommended based on cross-examining the proposed solid phases in the flowsheet with relevant data from the literature. Mixtures of solid components were added to aliquots of the liquid batches and sub-sampled to measure particle size distribution. The measured distribution was compared to independently created benchmark distributions appropriate for each simulant. This process was iterated until a solid phase composition that resulted in a representative particle size distribution was found. After the final compositions were confirmed, a suite of chemical and physical characterization data was collected. This report describes the simulant basis, formulation methodology, laboratory measurements, and data collected for the recipes recommended to represent 200W waste feeds.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Using Dynamic Interface Modeling and Simulation to Develop a Launch and Recovery Flight Simulation for a UH-60A Blackhawk

Joint Shipboard Helicopter Integration Process (JSHIP) is a Joint Test and Evaluation (JT&E) program sponsored by the Office of the Secretary of Defense (OSD). Under the JSHDP program is a simulation effort referred to as the Dynamic Interface Modeling and Simulation System (DIMSS). The purpose of DIMSS is to develop and test the processes and mechanisms that facilitate ship-helicopter interface testing via man-in-the-loop ground-based flight simulators. Specifically, the DIMSS charter is to develop an accredited process for using a flight simulator to determine the wind-over-the-deck (WOD) launch and recovery flight envelope for the UH-60A ship/helicopter combination. DIMSS is a collaborative effort between the NASA Ames Research Center and OSD. OSD determines the T&E and warfighter training requirements, provides the programmatics and dynamic interface T&E experience, and conducts ship/aircraft interface tests for validating the simulation. NASA provides the research and development element, simulation facility, and simulation technical experience. This paper will highlight the benefits of the NASA/JSHIP collaboration and detail achievements of the project in terms of modeling and simulation. The Vertical Motion Simulator (VMS) at NASA Ames Research Center offers the capability to simulate a wide range of simulation cueing configurations, which include visual, aural, and body-force cueing devices. The system flexibility enables switching configurations io allow back-to-back evaluation and comparison of different levels of cueing fidelity in determining minimum training requirements. The investigation required development and integration of several major simulation system at the VMS. A new UH-60A BlackHawk interchangeable cab that provides an out-the-window (OTW) field-of-view (FOV) of 220 degrees in azimuth and 70 degrees in elevation was built. Modeling efforts involved integrating Computational Fluid Dynamics (CFD) generated data of an LHA ship airwake and integrating a real-time ship motion model developed based on a batch model from Naval Surface Warfare Center. Engineering development and integration of a three degrees-of-freedom (DOF) dynamic seat to simulate high frequency rotor-dynamics dependent motion cues for use in conjunction with the large motion system was accomplished. The development of an LHA visual model in several different levels of resolution and an aural cueing system in which three separate fidelity levels could be selected were developed. VMS also integrated a PC-based E&S simFUSION system to investigate cost effective IG alternatives. The DIMSS project consists of three phases that follow an approved Validation, Verification and accreditation (VV&A) process. The first phase will support the accreditation of the individual subsystems and models. The second will follow the verification and validation of the integrated subsystems and models, and will address fidelity requirements of the integrated models and subsystems. The third and final phase will allow the verification and validation of the full system integration. This VV&A process will address the utility of the simulated WOD launch and recovery envelope. Simulations supporting the first two stages have been completed and the data is currently being reviewed and analyzed.

Sweeney, Christopher↗

Effect of cold cap coverage and emissivity on the plenum temperature in a pilot-scale waste vitrification melter

Integrated computational fluid dynamics (CFD) models are being developed to model the complex physics occurring within the high-level waste melter for vitrification of legacy tank waste at the Hanford site. Here, we present a validation of the integrated CFD model by using data from two experimental runs in a pilot-scale melter. While the model uses several simplifying assumptions (such as constant heat sinks from a cooling jacket and inleakage of ambient air, steady state feed-to-batch conversion heat, and a cold cap model with a simplified shape), it closely predicts the molten glass (1150°C and 1175°C) and plenum temperatures (550°C) obtained from thermocouples during two pilot-scale tests, with an average cold cap coverage of 80%. Additional simulations were performed to explore the sensitivity of the predicted plenum temperatures to variations in cold cap coverage (fraction of melt surface covered by the glass batch) and batch emissivity. The plenum temperature was found to be in the range of 606°C when cold cap coverage decreased from 95% to 70%. Cold cap emissivity had a smaller effect, increasing the plenum temperature by as much as 179°C when cold cap emissivity increased from 0.2 to 0.8. Maintaining a high cold cap coverage without overfeeding is important for a sustained melter operation with high glass throughput. This work provides a tool for achieving that goal in terms of correlating the plenum temperature with the cold cap coverage.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Solid phase characterization and transformation of illite mineral with gas-phase ammonia treatment

In situ remediation applications of ammonia (NH 3 ) gas have potential for sequestration of subsurface contamination. Ammonia gas injections initially increase the pore water pH leading to mineral dissolution followed by formation of secondary precipitates as the pH is neutralized. However, there is a lack of understanding of fundamental alteration processes due to NH 3 treatment. In these batch studies, phyllosilicate minerals (illite and montmorillonite) were exposed to NH 3 gas with subsequent aeration to simulate in situ remediation. Following treatments, solids were characterized using a variety of techniques, including X-ray diffraction, N 2 adsorption-desorption analysis for surface area, Fourier transform infrared (FTIR) spectroscopy, nuclear magnetic resonance (NMR), and microscopy methods to investigate physicochemical transformations. The results of this study indicate that, at high pH, the clays are altered as observed by differences in morphology and particle size via microscopy. However, the two clays interact differently with NH 3 . While montmorillonite interlayers collapsed due to intercalation, illite layers were unaffected as confirmed by FTIR analysis. Further, structural changes in silicate ([SiO 4 ] n- ) and aluminol (Al-OH) groups were identified by NMR and FTIR. This research showed that mineral alteration processes occur during and after NH 3 gas treatment which may be used to remove radionuclides from the aqueous phase through sorption, co-precipitation, and coating with secondary phyllosilicate alteration products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Adaptive Computing for Scale-Up Problems

Adaptive Computing is an application-agnostic outer loop framework to strategically deploy simulations and experiments to guide decision making for scale-up analysis. Resources are allocated over successive batches, which makes the allocation adaptive to some objective such as optimization or model training. The framework enables the characterization and management of uncertainties associated with predictive models of complex systems when scale-up questions lead to significant model extrapolation. A key advancement of this framework is its integration of multi-fidelity surrogate modeling, uncertainty management, and automated orchestration of various computing and experimentation resources into a single integrated software package. This enables efficient multi-fidelity modeling across multiple computing resources by incorporating real-world constraints such as relative queue times and throughput on individual machines into the multi-fidelity sampling decision. We discuss applications of this framework to problems in the renewable energy space, including biofuels production, material synthesis, perovskite crystal growth, and building electrical loads.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Dynamic simulation and stability analysis

Dynamic Interaction Simulation of Controls and Structure (DISCOS) program was developed for dynamic simulation and stability analysis of passive and actively controlled spacecraft. Program is written in FORTRAN IV for batch execution and requires access to finite-element structures program as NASTRAN for flexible-body input data.

Frisch, H. P.↗

Development and Validation of a Process Model and Open-Source Process Simulator for Microalgae-Based Tertiary Phosphorus Recovery

Microalgae-based tertiary wastewater treatment has the potential to meet stringent effluent phosphorus limits, with the added benefit of producing a marketable feedstock. However, the lack of validated mechanistic models and their implementation in process simulators have limited the adoption of this technology. In this study, an updated lumped pathway metabolic model (Phototrophic-Mixotrophic Process Model, PM 2 ), including both photoautotrophic and heterotrophic metabolisms of microalgae, was developed to predict effluent phosphorus concentration and biomass yield in response to dynamic influent and varying environmental conditions. The model was implemented in QSDsan – an open-source, Python-based design and simulation platform – for robust simulation under uncertainty. A global sensitivity analysis was performed to prioritize model parameters for calibration. The model was then calibrated and validated using batch experimental data and 45 days of continuous online monitoring data from a full-scale (568 m 3 ·d -1 ) microalgae-based tertiary wastewater treatment plant (EcoRecover process). In particular, along with dynamic influent composition, temperature and light intensity data with diel variation were provided as model inputs to reflect the microalgal behavior under day-night cycling. Overall, the QSDsan-based microalgae process simulator was able to predict effluent phosphorus within 0.02–0.04 mg-P·L -1 , while also capturing the general trends of state variables according to nutrient availability.

Lumped pathway metabolic model↗

Enabling machine learning-ready HPC ensembles with Merlin

With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computing (HPC) environment. Here, we present Merlin, a workflow framework to enable large ML-friendly ensembles of scientific HPC simulations. By augmenting traditional HPC with distributed compute technologies, Merlin aims to lower the barrier for scientific subject matter experts to incorporate ML into their analysis. As a producer–consumer workflow model, Merlin enables multi-machine, cross-batch job, dynamically allocated yet persistent workflows capable of utilizing surge-compute resources. Key features of Merlin are a flexible HPC-centric interface, low per-task overhead, multi-tiered fault recovery, and a hierarchical sampling algorithm that allows for $\mathscr{O}$(N) task execution and $\mathscr{O}$(N ln N) task queuing to ensembles of millions of tasks. In addition to Merlin’s design, we test the algorithm’s performance in an HPC center and demonstrate the ability to enqueue 40 million simulations in 100 s, with a 30 millisecond per-task overhead that is independent of ensemble size. Finally, we describe some example applications that Merlin has enabled on leadership-class HPC resources, such as the ML-augmented optimization of nuclear fusion experiments and the calibration of infectious disease models to study the progression of and possible mitigation strategies for COVID-19.

97 MATHEMATICS AND COMPUTING↗

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Effect of charge and anisotropic diffusivity on ion exchange kinetics in nuclear waste form materials

Motifs of radionuclide absorption particles (absorbers or fillers) for the nuclear waste treatment usually have crystal structures with nano-sized tunnels for fast ion diffusion. On the other hand, the accumulation of ions generates an electric field which affects the ion diffusion. Therefore, the strong anisotropic kinetic properties and the electric field affect ion exchange kinetics, especially in an assembly of these particles which may block the fast diffusion tunnels for each other. In this work, we developed a mesoscale model to describe the effect of anisotropic kinetic properties and electric field on ion exchange kinetics. With the model we simulated the effect of anisotropic diffusivity, electrochemical potential, particle morphology, size and aggregation on ion exchange kinetics during batch experiments. Here, tt is found that 1) ion exchange results in charge accumulation inside the particle because different ions have different diffusivities and different diffusion driving forces. The electric neutrality is not preserved; 2) the electric field speeds up the slower diffusion ions (either uptake or release) and slows down the faster diffusion ions to reduce charge accumulation and maintains charge neutrality; 3) the ion exchange becomes very slow at the later stages because a) diffusion driving force (electrochemical gradient, Cs + concentration in the solution and Na + inside the particle) continuously decreases, as observed in batch experiments; and b) the diffusivity decreases due to the Cs + concentration dependence of ion diffusivity. This can explain the observed experimental phenomena that the radionuclide capacity usually cannot be reached; 4) for needle-like particle and their clusters, both surface area and diffusion distance along the direction with the fastest diffusivities have significant impact on the ion exchange kinetics.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Batch Measurement Extremum Seeking Control of Distributed Energy Resources to Account for Communication Delays and Information Loss

Distributed Energy Resources (DER) have great potential to enhance the operation of electric power distribution systems. Previously, we explored the use of 2 Dimensional Extremum Seeking (2D-ES) control algorithms to enable model-free optimal control of DER to provide grid services to both the distribution and transmissions systems. Motivated by preliminary deployments of DER managed by 2D-ES algorithms in hardware-in-the-loop tests and in operational distribution grids, in this work, we extend the control scheme to accommodate communication delays and information loss. We propose a modification to the 2D-ES scheme to allow for the processing of batches of possibly noncontiguous objective function measurements at unknown and possibly uneven intervals. We provide a proof of the convergence of the batch 2D-ES (2D-BES) scheme when optimizing a generic convex objective function, as well as simulation results that demonstrate the suitability of the approach for substation active and reactive power target tracking.

Sankur, Michael↗

Preliminary Investigation of Target keff Search for Fuel Feed in MSRs

Fluid fuel molten salt reactor (MSR) concepts may achieve continuous-at-power operation with continuous or batch fresh fuel feed to counter reactivity loss from fuel burnup and fission product buildup. Accurate modeling and simulation of MSRs requires consideration of online fuel feed during depletion as well as a capability to determine the required feed amount to maintain a target multiplication factor (keff). As a first step supporting a full implementation of a new search capability into the SCALE code system’s TRITON reactor physics and depletion sequence, a preliminary investigation was performed to identify the best approach for a critical feed search capability considering computation time and accuracy. Using a simple iterative approach in which initial feed amounts were guessed based on previous depletion steps, the main drivers of computation time were found as the Monte Carlo neutron transport settings and the applied tolerance for the target keff. In the applied model, with tight convergence criteria (keff convergence to 8 pcm and a tolerance of ±15 pcm), the critical search calculation required twice the number of neutron transport calculations compared to a depletion calculation without critical search. In contrast, with a relaxed tolerance of ±50 pcm, the overall computation time for the critical search for merely increased by 30-40%. Further studies will be performed before implementing the critical search capability into TRITON.

Faulkner, Jonathon [ORNL]↗

Combined Technologies for In Situ Remediation of Tc-99 and U in Subsurface Sediments

In this study, combinations of chemical remedies were tested in bench-scale batch experiments to evaluate a two-step reduction-sequestration approach to effectively stabilize high concentrations of inorganic contaminant mixtures. Bench tests simulated contaminant and geochemical conditions of a perched aquifer located within the Central Plateau at the Hanford Site, located in southeastern Washington State (USA). Pairwise combinations of a reductant [e.g., zero valent iron, sulfur modified iron (SMI), or calcium polysulfide] and a sequestering agent [e.g., calcite, apatite, or dilute alkaline solution (e.g., NaOH)] were evaluated for immobilization and stabilization of technetium (Tc) (50,000 pCi/L), uranium (U) (150 mg/L), and nitrate (NO 3 ) (200 mg/L) in high ionic strength groundwater. The results of these batch studies demonstrated that reduction by SMI and sequestration in apatite or calcite are the most effective combination for these contaminant mixtures and conditions. Aqueous concentrations of Tc and U decreased by 95.6% ± 2.5% and 101.1% ± 5.2%, respectively, with SMI-apatite and 98.3% ± 0.0% and 101.2% ± 5.2%, respectively, with SMI-calcite. Sequential extractions showed that sequestered contaminants had limited capacity for re-oxidation; in fact, less than 10% of immobilized Tc and U was recovered by selective extraction of mineral phases most susceptible to oxidation. In addition, X-ray absorption near edge structure analysis of the sediment samples treated with SMI-calcite showed the presence of only U(IV), while both U(IV) and U(VI) were present in the SMI apatite combination [ratio of 0.43 U(IV):0.59 U(VI)].This study describes preliminary results that a two-step approach for stabilizing contaminant mixtures of long-lived radionuclides can be effective at reducing contaminant fluxes to groundwater from vadose and perched water zones.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗