SEARCH · Engineering Papers

Results for “flow analysis”

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 613 records · Page 34

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗

DOE OSTI2024

Higher-order tails and RG flows due to scattering of gravitational radiation from binary inspirals

Abstract We establish and develop a novel methodology to treat higher-order non-linear effects of gravitational radiation that is scattered from binary inspirals, which employs modern scattering-amplitudes methods on the effective picture of the binary as a composite particle. We spell out our procedure to study such effects: assembling tree amplitudes via generalized-unitarity methods and employing the closed-time-path formalism to derive the causal effective actions, which encompass the full conservative and dissipative dynamics. We push through to a new state of the art for these higher-order effects, up to the third subleading tail effect, at order$$ {G}_N^5 $$ G N 5 and the 5-loop level, which corresponds to the 8.5PN order. We formulate the consequent dissipated energy for these higher-order corrections, and carry out a renormalization analysis, where we uncover new subleading RG flow of the quadrupole coupling. For all higher-order tail effects we find perfect agreement with partial observable results in PN and self-force theories, where available.

Physics↗

DOE OSTI2025

Observation of Nonaxisymmetric Standard Magnetorotational Instability Induced by a Free-Shear Layer

The standard magnetorotational instability (SMRI) with a magnetic field component parallel to the rotation axis is widely believed to be responsible for the fast accretion in astronomical disks. In conventional base flows with a Keplerian profile or an ideal Couette profile, most studies focus on axisymmetric SMRI, since excitation of nonaxisymmetric SMRI in such flows requires a magnetic Reynolds number (Rm) more than an order of magnitude larger. Here, we report that, in a magnetized Taylor-Couette flow, nonaxisymmetric SMRI with an azimuthal mode number m = 1 can be triggered by a free-shear layer in the base flow at Rm ≳1, the same threshold as for axisymmetric SMRI. Global linear analysis reveals that the free-shear layer reduces the required Rm, possibly by introducing an extremum in the vorticity of the base flow. Nonlinear simulations validate the results from linear analysis and confirm that a novel instability recently discovered experimentally [Wang et al., Nat. Commun. 13, 4679 (2022)] is the nonaxisymmetric m =1 SMRI. Further, our finding has astronomical implications as free-shear layers are ubiquitous in celestial systems, such as the disk-star boundary layer, the solar tachocline, and the edge of planet-opened gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Geomagnetically Induced Current Field Test on Large Grid-Connected Power Transformers: Analysis, Model Development, and Simulations

Geomagnetic-induced current (GIC) flow in power grids can cause undesirable effects such as transformer overheating, harmonics, higher reactive power demand, etc. Many simulation models have been developed to study these effects, but real-world verification on modern transformer designs is rare. Here, this paper presents the first long-duration GIC field test in the U.S. performed on high-voltage, grid-connected transformers featuring winding clamps and tie rods instead of conventional tie bars. Field measurements were taken to evaluate GIC effects. These measurements also aided in developing and validating thermal and electromagnetic transient (EMT) models of the transformers. During the test, significant current and voltage distortions were observed along with considerable transformer reactive power losses. Analysis of the field measurements showed that the transformers’ hottest spot was at the inner windings, and their k-factors were close to factory test and software default values. Thermal simulations indicated that the transformers would not violate their thermal limits even for a GIC waveform that peaks at about 200 A/phase. EMT simulations revealed that increased transformer loading may reduce GIC-induced reactive power demand and harmonics in certain scenarios. The study also highlighted potential inaccuracies in using the k-factor method to calculate transformer reactive power losses.

EMTDC↗

DOE OSTI2021

Experimental Measurements of the Wake of a Sphere at Subcritical Reynolds Numbers

This work experimentally investigated the flow phenomena and vortex structures in the wake of a sphere located in a water loop at Reynolds numbers of Re = 850, 1,250, and 1,700. Velocity fields in the wake region were obtained by applying the time-resolved stereoscopic particle image velocimetry (TR-SPIV) technique. From the acquired TR-SPIV velocity vector fields, the statistical values of mean and fluctuating velocities were computed. Spectral analysis, two-point velocity–velocity cross-correlation, proper orthogonal decomposition (POD) and vortex identification analyses were also performed. The velocity fields show a recirculation region that decreases in length with an increase of Reynolds numbers. The power spectra from the spectral analysis had peaks corresponding to a Strouhal number of St = 0.2, which is a value commonly found in the literature studies of flow over a sphere. The two-point cross-correlation analysis revealed elliptical structures in the wake, with estimated integral length scales ranging between 12% and 63% of the sphere diameter. The POD analysis revealed the statistically dominant flow structures that captured the most flow kinetic energy. It is seen that the flow kinetic energy captured in the smaller scale flow structures increased as Reynolds number increased. The POD modes contained smaller structure as the Reynolds number increased and as mode order increased. In addition, spectral analysis performed on the POD temporal coefficients revealed peaks corresponding to St = 0.2, similar to the spectral analysis on the fluctuating velocity. The ability of POD to produce low-order reconstructions of the flow was also utilized to facilitate vortex identification analysis, which identified average vortex sizes of 0.41D for Re1, 0.33D for Re2, and 0.32D for Re3.

Engineering↗

DOE OSTI2022

Determination of fundamental mechanical properties of biomass using the cubical triaxial tester to model biomass flow

The flowability of biomass is an indicator of how amenable biomass is to handling. The flowability measurements are often determined with empirical or tertiary experiments. This status quo impedes employing engineering principles to advance the design and operation of biomass handling systems. It is imperative to establish an experimental protocol minimizing empirical aspects of flowability characterization to account for the variability of biomass. This study demonstrates the operational principles of the large chamber cubical triaxial tester and the procedure of triaxial tests to quantify bulk flow behaviors of two milled biomass feedstocks, namely corn stover 2 mm (CS) and Douglas fir 1 mm (DF). The Mohr-Coulomb (MC), Drucker-Prager (DP) and modified Cam Clay (mCC) models are calibrated for both biomasses. Analysis indicates that CS will exhibit a cohesive flow with a larger cohesion coefficient of MC (3.8 ± 3.5 kPa) than DF (0.42 ± 0.9 kPa) and a larger d value of DP (6.9 kPa) than DF (0.0 kPa), respectively. CS exhibits a higher spring-back index of mCC (0.39 ± 0.05) than DF (0.27 ± 0.05) also suggesting handling issues, which agrees with the experiences in the industry. This study demonstrates the capability of a CTT in a quantitative investigation of biomass handling characteristics.

09 BIOMASS FUELS↗

DOE OSTI2022

SAS4A/SASSYS-1 Validation with EBR-II Tests Performed During the SHRT Testing Program

The EBR-II sodium-cooled fast reactor supported a wide range of fast reactor research activities. The final fifteen years of operations at EBR-II were used for experiments and tests to demonstrate the importance of passive safety in liquid metal reactors and the capability of such systems to provide strong passive responses during off-normal conditions. The Shutdown Heat Removal Test program at EBR-II provided test data supporting the validation of computer codes for design, licensing, and operation of LMRs, among other objectives. The program included nearly sixty tests, most of which were protected and unprotected loss of flow and loss of heat sink tests run at various initial powers and flow rates. This report documents the results of modeling and simulation of tests from the Shutdown Heat Removal Test program to support validation of Argonne’s SAS4A/SASSYS-1 fast reactor safety analysis code. SHRT-17 was a protected loss of flow test where a loss of electrical power to all sodium coolant pumps was simulated to demonstrate the effectiveness of natural circulation cooling characteristics. SHRT-45R was an unprotected loss of flow test where the control rod scram function of the plant protection system was disabled to demonstrate the effectiveness of EBR-II’s passive reactivity feedbacks. BOP-301 and BOP-302R were unprotected loss of heat sink tests that demonstrated how reactivity feedback effects driven by changes to the core inlet temperature can shut down the fission process. This validation activity leverages modeling and simulation efforts that originated under an International Atomic Energy Agency Coordinated Research Project for benchmark analysis of the EBR-II Shutdown Heat Removal Tests. The SAS4A/SASSYS-1 models originally developed for the CRP have been modernized to utilize new code features and modeling practices. Differences between the predicted and measured test data were explored and the causes of these differences were identified as being due to modeling approximations, uncertainty in transient component behavior, and instrumentation uncertainties. Overall, it was concluded that for SHRT-17, SHRT-45R, BOP-301, and BOP-302R there was good agreement between the flow, power, and temperature predictions and the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2021

Assimilation of a Coordinated Fleet of Uncrewed Aircraft System Observations in Complex Terrain: EnKF System Design and Preliminary Assessment

Abstract Uncrewed aircraft system (UAS) observations collected during the 2018 Lower Atmospheric Process Studies at Elevation—a Remotely Piloted Aircraft Team Experiment (LAPSE-RATE) field campaign were assimilated into a high-resolution configuration of the Weather Research and Forecasting Model using an ensemble Kalman filter. The benefit of UAS observations was assessed for a terrain-driven (drainage and upvalley) flow event that occurred within Colorado’s San Luis Valley (SLV) using independent observations. The analysis and prediction of the strength, depth, and horizontal extent of drainage flow from the Saguache Canyon and the subsequent transition to upvalley and up-canyon flow were improved relative to that obtained both without data assimilation (benchmark) and when only surface observations were assimilated. Assimilation of UAS observations greatly improved the analyses of vertical variations in temperature, relative humidity, and winds at multiple locations in the northern portion of the SLV, with reductions in both bias and the root-mean-square error of roughly 40% for each variable relative to the benchmark run. Despite these noted improvements, some biases remain that were tied to measurement error and/or the impact of the boundary layer parameterization on vertically spreading the observations, both of which require further exploration. The results presented here highlight how observations obtained with a fleet of profiling UAS improve limited-area, high-resolution analyses and short-term forecasts in complex terrain.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Physical Controls on Irrigation Return Flow Contributions to Stream Flow in Irrigated Alluvial Valleys

Irrigation can be a significant source of groundwater recharge in many agricultural regions, particularly in arid and semi-arid climates. Once infiltrated, irrigation recharge can travel via subsurface flowpaths that return to the river system in a lagged manner, supplementing natural streamflow weeks, months, or even years from when the irrigation was applied. In regions that experience low flows during summer and early fall, return flows can be a significant source of supplementary streamflow. Many water planning and operations models either ignore return flows or roughly approximate them with analytical solutions. Thus, return flows represent an important but often overlooked component of the hydrological exchange and overall water balance in agricultural regions. This study uses groundwater models to explore a wide range of factors that control irrigation return flow timing in irrigated alluvial valleys. A sensitivity analysis approach is used to assess how factors such as the extent of irrigated land adjacent to a stream, irrigation recharge rate, aquifer hydraulic conductivity, aquifer thickness, water table configuration, and seasonal fluctuations in stream stage control the timing of subsurface return flows. Modeling is conducted using MODFLOW models representing an irrigated alluvial valley adjacent to a stream. While a simplification of the full complexity in real systems, the models are a significant advancement from the analytical solution and provide new insight into the timescales of return flows over a broad range of possible conditions. To contextualize our modeling results, they are compared to an analytical solution commonly used for approximating return flows to evaluate its performance. Our findings show what factors and conditions influence return flow timing and control whether they contribute to streamflow over short term (months) or longer term (seasonal) time scales.

Ferencz, Stephen B.↗

DOE OSTI2023

Dynamic systems modeling of the spallation neutron source cryogenic moderator system to optimize transient control and prepare for power upgrades

Through support of the US Department of Energy's Office of Basic Energy Sciences, Oak Ridge National Laboratory has begun applying machine learning methods to improve accelerator and target performance of the Spallation Neutron Source (SNS). These methods are being applied to the control optimization and power upgrade of the SNS Cryogenic Moderator System (CMS). A numerical model of the CMS has been developed to study these optimizations and system modifications using EcosimPro. This paper compares steady-state and transient numerical results with experimental data. Control optimization studies focused on dampening mass flow, temperature, and pressure fluctuations during sudden losses of accelerator beam power. This analysis was conducted by adjusting five proportional-integral-derivative controllers connected to four flow control valves and one heater. Future efforts include power uprate studies focused on increasing the CMS cooling capacity. The current accelerator power is 1.4 MW; the first target station is being upgraded to 2.0 MW as part of the Proton Power Upgrade effort. The CMS cooling capacity is sufficient for 2.0 MW operation.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Microextraction-Single Particle-Inductively Coupled Plasma-Mass Spectrometry for the Direct Analysis of Nanoparticles on Surfaces

A novel employment of single particle-inductively coupled plasma-mass spectrometry (SP-ICP-MS) was developed, where a microextraction (ME) probe is used to sample nanoparticles from a surface and analyze them in a single analytical step. The effects of several parameters on the performance of ME-SP-ICP-MS were investigated, including the flow rate, choice of carrier solution, particle size, and the design of the microextraction probe head itself. The optimized ME-SP-ICP-MS technique was used to compare the extraction efficiency (EE, defined as the ratio of particles measured to particles deposited on the surface) of the commercial probe head to a newly designed SP polyether ether ketone (PEEK) probe head. The SP PEEK probe head was found to have increased EE compared to the commercial probe head (8.5 ± 3% vs 3.9 ± 3%, respectively). Increasing the carrier solution flow rate was found to decrease the total analysis time at the cost of decreasing EE. Extraction efficiencies for ME-SP-ICP-MS were typically 4–10%, which is similar to transport efficiencies (1–10%) for conventional SP-ICP-MS. Lastly, ME-SP-ICP-MS was employed for the analysis of nano- and microparticles. The sizes of gold nanoparticles, 30 ± 3 and 51 ± 1.9 nm (certified sizes), and iron-based microparticles, 1000 ± 50 nm (certified size), were accurately determined to be 32.2 ± 2.5, 50.8 ± 3.4, and 1030 ± 57 nm, respectively, by ME-SP-ICP-MS. Further, this work demonstrates the potential of ME-SP-ICP-MS for the direct analysis of particles on common collection surfaces (GSR tabs, carbon planchettes, etc.) while retaining spatial information on particle distribution across the surface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Differential Impact of Climate Change on the Hydropower Economies of Karakoram and Himalayan River Basins

Water stored in the form of snow and glaciers in the High Mountain Asia (HMA) region regulates the water supply, and resultant water-based economies, that support the livelihoods of millions of people. Trends in the seasonal and long-term melting of snow and glaciers, governed by initial ice reserves, meteorological factors and geographic features, vary across sub-basins in the HMA region. We examined the economic impacts of climate-led changes in river flow in two drainage basins, one each from the Karakoram and Central Himalaya region. We used an integrated assessment framework to estimate the changes in economic value of the hydropower generation from hydropower plants on rivers fed by snow and glacier melt in the two sub-basins. The framework, developed under a NASA High Mountain Asia project, coupled biophysical models (a suite of climate models, snow/glacier-hydrology, and hydropower model) with economic analysis. We compared the differences in estimated river flow over historic and future time using the water balance model in sixteen scenarios (eight climate models and two emissions scenarios) for rivers upstream of hydropower plants in each sub-basin. Using the hydropower model we developed, we estimated the changes in hydropower generation at the Naltar IV hydropower plant, with an 18 MW capacity, located in Hunza, Karakoram, and the Trishuli hydropower plant, with a 19.6 MW capacity, in Trishuli, Central Himalaya. When compared to their baselines, the estimated impact of climate change and temporal variability were higher for the Naltar plant than for the Trishuli plant. Our sensitivity analysis shows that hydropower plants with water storage facilities help reduce the impact of changes, but the estimated impacts are higher for the higher capacity plants. This study provides an example of the differential impacts of climate change on hydropower plants located in rivers fed by varying amounts of snow and glacier melt at different decades in this century. This type of integrated assessment of climate change impact will support the scientific understanding of hydrologic flow and its impacts on a hydropower economy under various climate scenarios, as well as generate information about water resource management in a changing climate.

climate change↗

DOE OSTI2025

Thermal oxidation of nuclear graphite and pyrolytic carbon coatings

The oxidation of pyrolytic carbon (PyC) deposited via fluidized bed chemical vapor deposition was characterized and compared with that of standard nuclear-grade graphite. The materials were heated at 700 to 1000 °C in a thermogravimetric analysis system under 20% v/v O 2 flow, allowing for direct comparison of dynamic oxidative mass change in each material. Further, three different PyC samples fabricated under different conditions exhibited variation in total mass loss and mass loss rate, varying by as much as 709 mg/cm 2 in total mass loss and 14.2 (mg/cm 2 )/min in mass loss rate at a single temperature. These variations highlight the correlation between PyC microstructure/defect density and oxidation susceptibility. Additionally, changes in the microstructure and composition between PyC and graphite were characterized via scanning electron microscopy and correlated to the mass loss results. The results of this work have implications toward the safety of tristructural isotropic (TRISO) and other coated particle fuels, especially under off-normal conditions, given the limited information that exists about the oxidation behavior of PyC.

36 MATERIALS SCIENCE↗

DOE OSTI2022

ASDFL: An adaptive super‐pixel discriminative feature‐selective learning for vehicle matching

Abstract There are a large number of cameras in modern transportation system that capture numerous vehicle images continuously. Therefore, automatic analysis of these vehicle images is helpful for traffic flow management, criminal investigations and vehicle inspections. Vehicle matching, which aims to determine whether two input images depict an identical vehicle, is one of the core tasks in vehicle analysis. Recent relevant studies have focused on local feature extraction instead of global extraction, since local details can provide crucial cues to distinguish between cars. However, these methods do not select local features; that is, they do not assign weights to local features. Therefore, in this research, we systematically study the vehicle matching task, and present a novel annotation‐free local‐based deep learning method called Adaptive super‐pixel discriminative feature‐selective learning (ASDFL) to address this issue. In ASDFL, vehicle images are segmented into clusters of super‐pixels of similar size by considering the location and colour similarities of pixels without using any component‐level annotation. These super‐pixels are deemed to be the virtual components of vehicles. Moreover, a convolutional neural network is used to extract the deep features of these virtual components. Thereafter, an instance‐specific mask generation module driven by the extracted global features is enhanced to produce a mask to select the most distinctive virtual components of each vehicle image pair in the feature space. Finally, the vehicle matching task is accomplished by classifying the selected virtual component features of each imaged vehicle pair. Extensive experiments on two popular vehicle identification benchmarks demonstrate that our method is 1.57% and 0.8% more accurate than the previous baselines in a vehicle matching task on the VeRi and VehicleID datasets, respectively, which demonstrates the effectiveness of our method.

Qin, Rong↗

DOE OSTI2023

Turbulence Modeling with Nek5000/RS, SOD2D and Alya

We present Validation and Verification (V & V) study of two high-order spectral element method (SEM) based Computational Fluid Dyanimcs (CFD) codes that we will utilize for work on turbulence modeling: Nek5000 and SOD2D. While the former solves the incompressible form of Navier-Stokes equation, the latter works with compressible set of equations and uses an entropy-viscosity formulation to account for the discontinuities for high Mach number flows. We demonstrate the accuracy of these codes for two benchmark sub-sonic turbulent flows: periodic channel and pipe flow, by carrying out first and second order statistical analysis, grid convergence and turbulent structure analysis using wall-resolved large eddy simulations (WRLES). Later, we report on the implementation and testing of various wall-modeling strategies for large eddy simulation of turbulent flows in Nek5000. These include both classical log-law based and decision-tree based machine-learning models. Accuracy of these closure strategies are analyzed and necessary future work is outlined.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Epsilon Metal: A Waste Form for Noble Metals from Used Nuclear Fuel

Epsilon metal (e-metal) is the metallic phase that forms as inclusions at the grain boundaries in the UO2 fuel during reactor operation. This metal is composed of Pd, Mo, Rh, Ru, and Tc. These metallic inclusions are insoluble in strong acid and remnants of these metallic inclusions have been found in the UO2 matrix that remains from the natural reactors in Gabon that were active 1.8 billion years ago, therefore e-metal should be an excellent waste form for the immobilization of the long-lived isotopes 107Pd (6.5 ? 106 a) and 99Tc (2.13 ? 105 a), with 99Tc being the isotope of interest for repository performance. Therefore, the chemical durability of this potential waste form is assessed in this study. Typically, corrosion rates for metallic materials are measured electrochemically because they are quick, inexpensive, and can reveal the mechanism by which a metal corrodes, at least initially. However, in a repository the waste form would be subjected to slowly flowing water without an applied electrical potential over long time periods. Therefore, the corrosion rates of e-metal specimens were measured with both electrochemical tests and the single-pass flow-through test (SPFT). Potentiodynamic and potentiostatic polarization results suggest that a thin passive film exists on the alloy surface, which seems to be responsible for its high corrosion resistance. Additionally, X-Ray photoelectron spectroscopic results suggest that Pd oxides are significantly enriched in the passive film Results from the SPFT show that the dissolution rates were weakly dependent on pH. Only Mo and Re were found in solution and were used for the calculation of the dissolution rates. In general, the electrochemically determined corrosion rates agree reasonably well with the initial dissolution rate measured with the SPFT test, but they are about one or two orders of magnitude higher than the steady state rates. The causes for this discrepancy are discussed.

epsilon-metal, nuclear waste form, corrosion, sing↗

DOE OSTI2022

The dispersion and propagation of topological Langmuir-cyclotron waves in cold magnetized plasmas

The topological Langmuir-cyclotron wave (TLCW) is a recently identified topological surface excitation in magnetized plasmas. We show that TLCW originates from the topological phase transition at the Langmuir wave-cyclotron wave resonance. By isofrequency surface analysis and two- and three-dimensional time-dependent simulations, we demonstrate that the TLCW can propagate robustly along complex phase transition interfaces in a unidirectional manner and without scattering. Because of these desirable features, the TLCW could be explored as an effective mechanism to drive current and flow in magnetized plasmas. As a result, the analysis also establishes a close connection between the newly instituted topological phase classification of plasmas and the classical Clemmow-Mullaly-Allis (CMA) diagram of plasma waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Multiphysics Simulations of MSRE with NEAMS Thermal Hydraulics Tools

This report documents the benchmarks being developed and simulations performed using tools and codes developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing MSRE experimental data. In FY23, three main work scopes were investigated under the NEAMS MSR work package at ANL. The first scope investigated the Griffin-SAM coupling model for simulating the pump startup transient experiment of MSRE. The analyses start with a simple model (single-channel, single-lattice), gradually adding more details (multi-channel, full-core) into the model. The results show that the reactivity loss curve is very sensitive to the axial boundary conditions and the radial core discretization. The simple model can predict a similar reactivity trend as that of the more sophisticated model, which is likely due to error cancellation. Accurately modeling the axial boundary condition may further improve the reactivity trend but would require significant efforts to generate the mesh of the MSRE inlet and upper plenum. The core channel radial discretization for the Griffin-SAM coupled model also depends on the flow distribution. Given the complex geometry in the inlet plenum, the flow distribution needed to be calculated from CFD analysis, which was performed using the NekRS code. This analysis employed a MSRE CAD model developed by Copenhagen Atomics. The CAD model was disassembled to keep the inlet plenum region only, which was subsequently cleaned and modified so that the mesh generated is under the memory limit. The results are merged to a few radial regions to show that the flow rate is highest in the central region. This would be useful for future improvement of the Griffin-SAM coupling model of the MSRE core. The last task investigated is tritium transport modeling using the standalone SAM code. This task aimed to initiate the effort to demonstrate and validate the tritium transport model implemented in SAM. The preliminary investigation employed an MSRE model consisting of the primary loop. Three tritium transport pathways were examined including the retention in the graphite, the permeation through the HX tube wall, and the removal from the off-gas system. The results compare well with the MSRE data, but improvements are still needed on the initial conditions (i.e., the present state may not have reached equilibrium), the boundary conditions, the off-gas system modeling, and a better numerical strategy to reach the equilibrium state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗