Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transient State”

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 469 records · Page 26

Employing Long-Range Inductive Effects to Modulate Metal-to-Ligand Charge Transfer Photoluminescence in Homoleptic Cu(I) Complexes

Four Cu(I) bis(phenanthroline) photosensitizers formulated from a new ligand structural motif (Cu1-Cu4) coded according to their 2,9-substituents, were synthesized, structurally characterized, and fully evaluated using steady-state and time-resolved absorption and photoluminescence (PL) measurements as well as electrochemistry. The 2,9-disubstituted-3,4,7,8-tetramethyl-1,10-phenanthroline ligands feature the following six-membered ring systems prepared through photochemical synthesis: 4,4-dimethylcyclohexyl (1), tetrahydro-2H-pyran-4-yl (2), tetrahydro-2H-thiopyran-4-yl (3), and 4,4-difluorocyclohexyl (4). Universally, these Cu(I) metal-to-ligand charge transfer (MLCT) chromophores display excited state lifetimes on the microseconds time scale at room temperature, including the three longest-lived homoleptic cuprous phenanthroline excited states measured to date in deaerated CH 2 Cl 2 , τ = 2.5 to 4.3 μs. This series of molecules also feature high PL quantum efficiencies (ΦPL = 5.3 to 12% in CH 2 Cl 2 ). Temperature-dependent PL lifetime experiments confirmed that all these molecules exhibit reverse intersystem crossing and display thermally activated delayed PL from a 1 MLCT excited state lying slightly above the 3 MLCT state, 1050 – 1490 cm -1 . Ultrafast and conventional transient absorption measurements confirmed that the PL originates from the MLCT excited state, which remains sterically arrested, preventing an excessive flattening distortion even when dissolved in Lewis basic CH 3 CN. Combined PL and electrochemical data provided evidence that Cu1-Cu4 are highly potent photoreductants (E ox * = -1.73 to -1.62 V vs. Fc +/o in CH 3 CN), whose potentials are altered solely based on which heteroatoms or substituents are resident on the 2,9-appended ring derivatives. It is proposed that long-range electronic inductive effects are responsible for the systematic modulation observed in the PL spectra, excited state lifetimes, as well as the ground state absorption spectra and redox potentials. Cu1-Cu4 quantitatively follow the energy gap law, correlating well with structurally related cuprous phenanthrolines and are also shown to triplet photosensitize the excited states of 9,10-diphenylanthracene with bimolecular rate constants ranging from 1.61 to 2.82 × 10 8 M –1 s –1 . Furthermore, the ability to tailor both photophysical and electrochemical properties using long-range inductive effects imposed by the 2,9-ring platforms advocates new directions for future MLCT chromophore discovery.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benchmark Modeling and Simulation of the FFTF LOFWOS Test #13 Using SAM

The Fast Flux Test Facility (FFTF) was a 400 MW thermal powered, oxide-fueled, liquid sodium cooled test reactor, built to assist development and testing of advanced fuels and materials for fast breeder reactors. In July 1986, a series of unprotected Loss of Flow Without Scram (LOFWOS) transients were performed in FFTF as part of the Passive Safety Testing (PST) program. The LOFWOS Test #13, which was initiated at 50% power and 100% flow with the pump pony motors left off, has been chosen as a benchmark case by IAEA to support collaborative efforts within international partnerships on the validation of simulation tools and models in the area of sodium fast reactor passive safety in an IAEA Coordinated Research Project (CRP), launched in October 2018. The System Analysis Module (SAM) is an advanced and modern system analysis tool under development at Argonne National Laboratory for advanced non-LWR safety analysis. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. To participate the IAEA CRP and enhance the SAM validation base for advanced reactor transient safety analysis, benchmark simulations of the FFTF LOFWOS Test #13 are performed using the SAM code. In this first phase of the validation effort, the thermal-hydraulic behavior of the reactor system is the focus and the reactor kinetics is not considered in the SAM FFTF model. Instead, the results of Argonne’s neutronics calculations are directly used, including the power shape of the active core region and the power history during the transient. The simulation results of FFTF at steady state agreed well with the measured data from the test. During the transient, reasonably good agreement were also obtained. Future work to improve the model will focus on introducing the reactivity predictions into the model, as well as better understanding or resolving the current discrepancies with the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

On the Recoverability of Reactor Dynamics from Point Kinetics Data using SINDYc

Data driven models for reactor dynamics tend to face a few notable challenges. First, the broad range of timescales present across the various feedback mechanisms. Next, high standards for safety and importance of model performance in all operational domains. Finally, the correlations between state variables observed in most transients leads to difficulties when methods attempt to attribute certain dynamic phenomena to a particular cause. The current paper seeks to support efforts towards incorporating physics knowledge into one particular data-driven method for finding reactor dynamics called ``Sparse Identification of Nonlinear Dynamics with Control (SINDYc). The incorporation of physical knowledge into SINDYc will help address the challenges listed above by giving the model an initial understanding of the system being modeled. The demonstration of an exact representation of a set of point reactor dynamics equations in SINDYc is provided. Then, this allows for further discussion with mathematical justification as to why SINDYc, and other data-driven methods, may be difficult to apply to nuclear reactor dynamics due to correlations in the state variables. A particularly useful result of this work is the set of candidate functions required in SINDYc to exactly represent the reactor dynamics under a point approximation. In the future, SINDYc can be integrated with point models for reactor dynamics before being applied to the physical system to yield higher accuracy.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗

Thermal and electric multidomain dynamic model for integration of power grid distribution with behind-the-meter devices

As renewable energy sources like solar and wind power become more integrated into the grid, coordinated control of behind-the-meter devices is crucial for enhancing grid flexibility and reliability and for meeting cost targets, with standardized models being developed to support this transition. The increasing flexibility and uncertainty of integrated renewable energy grids, along with interactions between various subsystems, make traditional steady-state modeling insufficient to capture transient and dynamic behaviors. Current models (e.g., composite load and battery equivalent models) focus on thermodynamic or electrical characteristics but overlook critical electromechanical interactions. This limits the ability to share performance information for grid services and hampers fast dynamic simulations. In addition, motor stalling is usually triggered by a fault event and attributed to the characteristics of the mechanical torque of the motor, resulting in absorption of a large amount of reactive power during the stalling period. Further, this significant withdrawal of reactive power will deteriorate the dynamic voltage stability of power grids and cause delayed voltage recovery. Therefore, an in-depth modeling of the thermodynamics or mechanical torque is essential to study the impacts of the realistic torque characteristics of those behind-the-meter devices on power system voltage stability. This study developed a dynamic multidomain model for building HVAC systems, such as air-source heat pumps, to simulate their thermal and electrical responses to grid transients. The model can accurately predict power metrics with a mean absolute percentage error of 10%, by validating against with power system computer-aided design performance data. Case studies demonstrate the model capability of capturing the transient response to sudden voltage changes, rapid load fluctuations, and system shutdowns respectively. During a sudden voltage drop (30% for 0.1s), a fully loaded heat pump’s motor speed dropped, continued declining, and shut down after 3.6s, with severe power oscillations and a torque spike. A partially loaded unit experienced temporary oscillations but stabilized. Under higher building loads, compressor speed increased from 64% to 100%, with power and torque rising before stabilizing. In safety-triggered shutdowns, power decreased after minor fluctuations, and torque briefly spiked before dropping to zero.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Thermal loading analysis of the ring injection dump for the Spallation Neutron Source facility

The ring injection dump (RID) is the largest beam dump in the Spallation Neutron Source facility at Oak Ridge National Laboratory and accepts a fraction of the beam from the ring that is not captured in the ring during injection. Thermo-fluid modeling of the RID components was performed during the initial design by using the software ANSYS CFX to characterize the power rating for the RID. From the calculations, the power rating was lowered from 200 kW to 150 kW, at 1.3 GeV beam energy due to concerns regarding the heating of the concrete structure. The Proton Power Upgrade project scope included reevaluating the power limit to the RID in anticipation that increasing it would provide more operational flexibility. The focus was on the shielding because it was the limiting factor, and data collected from thermocouples installed on the structure were used to benchmark the analysis. A three-stage process was adopted in the latest study. First, the already existing steady-state calculations were validated and extended with the latest Monte Carlo N-Particle source term calculations. Second, transient calculations were conducted to capture the dynamic state of the system in response to the energy dumped to the RID during the following operational period of the beam from 2005–2019 (14 years), and sensitivity analyses of crucial parameters were performed to benchmark the model with reasonable accuracy. Finally, an idealized transient cycle analysis was conducted with realistic duty factors to predict the temperature distribution at higher beam powers. Overall, the results confirmed the accuracy of the original steady state model but allowed the development of a new validated transient model with higher accuracy for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Advanced Fuels Campaign 2019 Accomplishments

The mission of the Advanced Fuels Campaign (AFC) is to perform or support research, development, and demonstration (RD&D) activities to identify and mature innovative fuels and cladding materials and associated technologies with the potential to improve the performance and safety of current and future reactors; increase the efficient utilization of nuclear energy resources; contribute to enhancing proliferation resistance of the nuclear fuel cycle; and address challenges related to waste management and ultimate disposal. AFC pursues its mission objectives using a goal-oriented, science-based approach that seeks to establish a fundamental understanding of fuel and cladding behaviors under conditions that arise during fabrication, normal steady-state irradiation, off-normal transient scenarios, and storage/disposal. This approach includes advancing the theoretical understanding of fuel behavior, conducting fundamental and integral experiments, and supporting the mechanistic, multi-scale modeling of nuclear fuels to inform and guide fuel development projects, advance the technological readiness of promising fuel candidates, and ultimately support fuel qualification and licensing initiatives. In the area of advanced tools for the modeling and simulation of nuclear fuels, the AFC works in close partnership with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, participating with NEAMS in developing mechanistic fuel behavior models and providing experimental data to inform and validate its most advanced tools. Specifically, AFC objectives in thecoming five year horizon include: 1. support the industry-led development of Accident Tolerant Fuel (ATF) technologies with improved reliability and performance under normal operations and enhanced tolerance during hypothetical accident scenarios, with implementation of lead test rods/assemblies of one or more ATF concepts in commercial reactor(s) by or before 2022 and batch reloads of near-term concept(s) in the 2023-2026 timeframe; 2. lead research and development on innovative fuel and cladding technologies for applications to future advanced reactors, especially fast-spectrum reactors, including reactors that utilize both once-through and recycle scenarios; 3. continue the development and demonstration of a multi-scale,science-based approach to fuel development and testing, andcontribute to the establishment of a state-of-the-art R&D infrastructure necessary to accelerate the development of new fuel concepts; and collaborate with NEAMS on the development and validation of multi-scale, multi-physics, and increasingly predictive fuel performance models and codes. 4. collaborate with NEAMS on the development and validation of multi-scale, multi-physics, and increasingly predictive fuel performance models and codes.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Elucidating Structural Disorder in Ultra‐Thin Bi‐Rich Bismuth Oxyhalide Photocatalysts

Abstract Advancing the field of photocatalysis requires the elucidation of structural properties that underpin the photocatalytic properties of promising materials. The focus of the present study is layered, Bi‐rich bismuth oxyhalides, which are widely studied for photocatalytic applications yet poorly structurally understood, due to high levels of disorder, nano‐sized domains, and the large number of structurally similar compounds. By connecting insights from multiple scattering techniques, utilizing electron‐, X‐ray‐ and neutron probes, the crystal phase of the synthesized materials is allocated as layered Bi 24 O 31 X 10 (X = Cl, Br), albeit with significant deviation from the reported 3D crystalline model. The materials comprise anisotropic platelet‐shaped crystalline domains, exhibiting significant in‐plane ordering in two dimensions but disorder and an ultra‐thin morphology in the layer stacking direction. Increased synthesis pH tailored larger, more ordered crystalline domains, leading to longer excited state lifetimes determined via femtosecond transient absorption spectroscopy (fs‐TAS). Although this likely contributes to improved photocatalytic properties, assessed via the photooxidation of benzylamine, increasing the overall surface area facilitated the most significant improvement in photocatalytic performance. This study, therefore, enabled both phase allocation and a nuanced discussion of the structure‐property relationship for complicated, ultra‐thin photocatalysts.

Chemistry↗

Thermal hydraulics modeling of On-line refueling for the Advanced High Temperature Reactor (AHTR)

As part of the effort to move forward the Fluoride salt-cooled High-temperature Reactor (FHR) technology, the Advanced High Temperature Reactor (AHTR) is being developed at Oak Ridge National Laboratory (ORNL) and several other institutions. Due to its plank fuel design and inherently low heavy metal loading, which challenges fuel utilization, online refueling is considered as an option for improving the economic viability of this reactor. This paper presents the thermal hydraulic modeling of the online refueling for the AHTR. Analyses at the single channel and fuel assembly level have been performed in order to develop a model of the reactor. The reactor modeling has been integrated with CFD studies for characterizing the steady state conditions and the refueling transient. The simulations of the transient have been complemented by the analysis of the stability of the refueling operation. This work demonstrates the viability of online refueling from a thermal hydraulic standpoint and develops a modeling approach for this type of operational transient.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Code validation of SAM using natural-circulation experimental data from the compact integral effects test (CIET) facility

The primary objective of this study is to validate the system analysis code, SAM, using experimental data from the Compact Integral Effects Test (CIET) experimental loop. SAM is a modern system analysis code being developed at Argonne National Laboratory for safety analysis of designs for advanced non-light water reactors (non-LWRs), such as sodium-cooled fast reactors, high-temperature gas-cooled reactors, and fluoride salt-cooled high-temperature reactors (FHRs). To support SAM code development for the wide range of non-LWR applications, it is of paramount importance to validate the code against experiments highly relevant to these reactor concepts. Additionally, the CIET facility, which was designed to reproduce the thermal-hydraulics response of FHRs under both forced- and natural-circulation conditions, has been identified and selected as one of the benchmark test facilities for SAM code validation. In this study, two sets of available CIET tests were selected for SAM code validation purposes, namely, power step change transient tests and steady-state natural-circulation tests. For all selected tests, SAM-predicted results show very good agreement with experimental data. The successful validation of SAM against these selected CIET experiments demonstrates that the computer code is well suited for thermal-hydraulics analysis of FHR designs.

42 ENGINEERING↗

Loss and Isotopic Fractionation of Alkali Elements during Diffusion-Limited Evaporation from Molten Silicate: Theory and Experiments

Moderately volatile elements (MVEs) are variably depleted in planetary bodies, reflecting the imprints of nebular and planetary processes. Among MVEs, Na, K, and Rb are excellent tracers for unraveling the history of MVE depletion in planetary bodies because they have similar geochemical behaviors but can be chemically fractionated by evaporation and condensation processes. Furthermore, K and Rb are amenable to high-precision isotopic analyses, which can help constrain the conditions of evaporation and condensation. To quantitatively understand why Na, K, and Rb are depleted in planetary bodies, we have carried out vacuum evaporation experiments from basaltic melt at 1200 and 1400 °C to study their evaporation kinetics and isotopic fractionations. We chose this composition because it is relevant to evaporation from small differentiated planetesimals. The Rb isotopic compositions of the evaporation residues were measured by multicollector inductively coupled plasma mass spectrometry (MC-ICPMS), and the K isotopic compositions were measured along profiles across the residues by secondary ion mass spectrometry (SIMS). In the 1400 °C run products, we found that the concentrations of both K and Rb in the run products decreased from core to rim, which was accompanied by a heavy K isotope enrichment near the surface. This indicates that, in this run, evaporation was limited by diffusion. To use those data quantitatively, we derive analytical equations that describe the evaporation rate and isotopic fractionation associated with diffusion-limited evaporation from a sphere, slab, and cylinder in transient and quasi-steady state regimes. This model is used to tease out the roles that diffusive transport in the melt and evaporation at the melt/gas interface play in setting the elemental depletion and isotopic composition of the residue. Under our experimental conditions, volatility decreases in the order of Na, Rb, and K. Using our experimental results in a thermodynamic model, we have estimated the product γΓ of activity coefficients × evaporation coefficients of Na, Rb, and K. The measured isotopic compositions of the residues are well explained using Rayleigh distillations, whereby the relative volatilities of K and Rb isotopes are given by the square root of their masses. We use our results and previously published data to predict how K and Rb could have been lost as a function of temperature, melt composition, oxygen fugacity, and saturation degree relevant to Vesta’s building blocks. We find that the K and Rb depletions, K/Rb elemental fractionation, and δ 41 K and δ 87 Rb isotopic fractionations of Vesta (as sampled by howardite-eucrite-diogenite (HED) meteorites) are best explained by evaporation of submillimeter size objects for 0.1-10 years at moderate temperatures (~1050 °C) in a medium ~98.8% saturated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Maven: a multimodal foundation model for supernova science

Abstract A common setting in astronomy is the availability of a small number of high-quality observations, and larger amounts of either lower-quality observations or synthetic data from simplified models. Time-domain astrophysics is a canonical example of this imbalance, with the number of supernovae observed photometrically outpacing the number observed spectroscopically by multiple orders of magnitude. At the same time, no data-driven models exist to understand these photometric and spectroscopic observables in a common context. Contrastive learning objectives, which have grown in popularity for aligning distinct data modalities in a shared embedding space, provide a potential solution to extract information from these modalities. We present Maven, the first foundation model for supernova science. To construct Maven, we first pre-train our model to align photometry and spectroscopy from 0.5 M synthetic supernovae using a contrastive objective. We then fine-tune the model on 4702 observed supernovae from the Zwicky transient facility. Maven reaches state-of-the-art performance on both classification and redshift estimation, despite the embeddings not being explicitly optimized for these tasks. Through ablation studies, we show that pre-training with synthetic data improves overall performance. In the upcoming era of the Vera C. Rubin observatory, Maven will serve as a valuable tool for leveraging large, unlabeled and multimodal time-domain datasets.

Zhang, Gemma (ORCID:0000000280198082)↗

Giant photoinduced lattice distortion in oxygen vacancy ordered SrCoO 2.5 thin films

Despite the progress in understanding the formation and migration of oxygen vacancies in strongly correlated complex oxides, few studies focused on the dynamic behaviors of the oxygen vacancies under the transient and non-equilibrium states. Here we report a series of multi-timescale ultrafast X-ray diffraction experiments using picosecond synchrotron and femtosecond table-top X-ray sources to monitor the structural dynamics in oxygen-vacancy-ordered SrCoO 2.5 thin films excited by photons. A giant photoinduced structure distortion in c-axis direction (with strain $Δc/c$ > 1%) was observed, higher than any previously reported data in the other transition metal oxide films. The femtosecond X-ray diffraction reveals the formation and propagation of the coherent acoustic phonons, indicating an instantaneous and a much larger photoinduced stress within the first 100 ps. Density Function Theory reproduced the photostrictive responses and the strong dependence on excitation wavelength as observed in the experiment. The combined experimental and theoretic results demonstrate that the photoexcitation of bonding to antibonding state via charge transfer is the dominant mechanism in the SrCoO 2.5 thin films distinct from the depolarization effect by photoinduced carriers in the other perovskite oxides.

36 MATERIALS SCIENCE↗

Nanosecond Optically Induced Phase Transformation in Compressively Strained BiFeO 3 on LaAlO 3

Above-band-gap optical illumination of compressively strained BiFeO 3 induces a transient reversible transformation from a state of coexisting tilted tetragonal-like and rhombohedral-like phases to an untilted tetragonal-like phase. Time-resolved synchrotron x-ray diffraction reveals that the transformation is induced by an ultrafast optically induced lattice expansion that shifts the relative free energies of the tetragonal-like and rhombohedral-like phases. The transformation proceeds at interfaces between regions of the tetragonal-like phase and regions of a mixture of tilted phases, consistent with the motion of a phase boundary. The optically induced transformation demonstrates that there are new optically driven routes towards nanosecond-scale control of phase transformations in ferroelectrics and multiferroics.

36 MATERIALS SCIENCE↗

Development and Validation of SAM Multi-dimensional Flow Model for Thermal Mixing and Stratification Modeling

Thermal mixing and stratification in large pools or enclosures are very important phenomena that are critical to nuclear reactor safety. Because of the wide ranges of time and length scales associated with such phenomena, accurate modeling and simulation of thermal mixing and stratification remain as the key unresolved, challenging problems for reactor transient analyses. In traditional system analysis codes, simplified zero-dimensional (0-D) models are widely used for their high numerical efficiency, but they generally suffer from very limited prediction accuracies or range of applicability. Like traditional system analysis codes, the current version of the SAM code has implemented such simplified 0-D and one-dimensional mixing models. On the other hand, high-resolution Computational Fluid Dynamics (CFD) tools are often used to model complex thermal mixing and stratification phenomena. They are, however, generally numerically expensive, and they require large amounts of computational resources. It is therefore desirable to implement advanced and efficient thermal mixing and stratification modeling capabilities embedded in a system analysis code. This approach will improve the accuracy of reactor safety analyses when thermal mixing and stratification are involved, and also avoid using the large computational resources needed for high-resolution CFD analysis. Currently, with the support of the U.S. DOE Office of Nuclear Energy’s Nuclear Energy Advanced Modeling and Simulation program, an effort has been launched to develop and implement a multi-dimensional flow model in the system analysis code SAM, and demonstrate its applications to model thermal mixing and stratification phenomena in large enclosures. The main outcomes of this research and development activity are summarized in this report, which presents an attempt to include a built-in advanced multi-dimensional flow model in a system analysis code with the focus on overcoming the simulation challenges of thermal mixing and stratification phenomena. In this report, we start with the introduction of existing SAM code capabilities to simulate thermal mixing and stratification phenomena, which is followed by a short summary of the multi-dimensional model implemented in the SAM code, including both the physical model and the Finite Element Method code implementation. In this study, two options were implemented in the SAM code to model turbulent flows: a relatively simple built-in turbulence model and an interface to accept externally computed turbulent viscosities (e.g., from a high- resolution CFD simulation). Code validation studies on this newly added capability were then carried out to compare SAM simulation results with experimental data from the SUPERCAVNA facility, which was designed to study the complex flow recirculation and thermal stratification phenomena relevant to sodium fast reactor designs. In this study, one transient and two steady-state test cases were used for code validation. Different approaches have been used to model the complex turbulence flow fields in the SUPERCAVNA facility. A highly simplified zero-equation turbulence model was first used, but it was determined that it is too simple to capture the complex turbulence flow fields in these test cases. Subsequently, the code validation continued with the use of turbulent viscosity data from high-resolution STAR-CCM+ CFD simulations to improve the accuracy of the results. Using this approach, the SAM simulation results showed very good agreement with both the SUPERCAVNA experimental data and STAR-CCM+ simulation results. In this report, we demonstrate the development, implementation, and successful validation of a multi-dimensional flow model in the SAM code, which aims to improve the simulation accuracy for complex thermal mixing and stratification phenomena. Lessons have also been learned, including that in cases where the flow fields are not well predicted by the zero-equation model, the iinput of turbulent viscosities from an external source can enhance the overall predictive capabilities needed to accurately capture complex thermal-hydraulic phenomena. Therefore, future research will be needed to further improve the code’s capabilities, such as by developing a more efficient and robust approach to capture the turbulence effects in the SAM code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

SAM Finite Volume Method Development Status Update: GCR Application, Restart, and MultiApp

The System Analysis Module (SAM) is being developed as a modern system analysis code for advanced non-light-water-reactor safety analysis under the U.S. DOE NEAMS program. Previous feasibility studies have demonstrated that a staggered-grid finite volume method (SG-FVM), implemented under the MOOSE framework, can deliver more than an order of magnitude speedup over the existing continuous Galerkin finite element method (CG-FEM) solver for liquid-cooled, incompressible but thermally expandable flow systems. This work extends the previous effort to compressible, gas-cooled reactor applications, where pressure couples directly into the mass equation adding additional nonlinearity into the equation system. New code capabilities are implemented for pebble bed high-temperature gas-cooled reactor (PB-HTGR) analysis, including a pebble bed CoreChannel component, built-in pebble bed effective thermal conductivity model and channel-to-channel crossflow model. The capabilities are tested, benchmarked, and demonstrated for problems with increased level of model and physical complexities, including the HTTU effective thermal conductivity test, the SANA passive cooling test, and a demonstration case using the GPBR200 reactor design covering steady-state operation, DLOFC and PLOFC transients. Across all cases, the SG-FVM solver demonstrated strong robustness and efficiency, and the solutions agree well with reference results and data. The finding of this work proves that SG-FVM is a viable and efficient solver pathway for compressible, gas-cooled reactor system analysis in SAM. In addition, work has been done to successfully support SAM-FVM recover/restart code feature that is essential to reactor safety analysis applications, and MultiApp code feature that is essential to multi-scale and multi-physics simulations. In summary, this work continued from previous feasibility studies, and further demonstrated that the SG-FVM will serve as a strong foundation for SAM’s advanced solver algorithm for future deployment.

Zou, Ling↗

Parametric Controls on Vegetation Responses to Biogeochemical Forcing in the CLM5

Future projections of land carbon uptake in Earth System Models are affected by land surface model responses to both CO 2 and nitrogen fertilization. The Community Land Model, Version 5 (CLM5), contains a suite of modifications to carbon and nitrogen cycle representation. Globally, the CLM5 has a larger CO 2 response and smaller nitrogen response than its predecessors. To improve our understanding of the controls over the fertilization responses of the new model, we assess sensitivity to eight parameters pertinent to the cycling of carbon and nitrogen by vegetation, both under present-day conditions and with CO 2 and nitrogen fertilization. The impact of fertilization varies with both model parameters and with the balance of limiting factors (water, temperature, nutrients, and light) in the pre-fertilization model state. The model parameters that impact the pre-fertilization state are in general not the same as those that control fertilization responses, meaning that goodness of fit to present-day conditions does not necessarily imply a constraint on future transient projections. Where pre-fertilization state has low leaf area, fertilization-induced increases in leaf production amplify the model response to the initial fertilization via further increases in photosynthesis. Model responses to CO 2 and N fertilization are strongly impacted by how much plant communities can increase their rates of nitrogen fixation and also directly affected by costs of N extraction from soil and stoichiometric flexibility. Illustration of how parametric flexibility impacts model outputs should help inform the interpretation of carbon-climate feedbacks estimated by in fully coupled Earth system model simulations.

54 ENVIRONMENTAL SCIENCES↗