Peaking Factor Uncertainty of VERA-MPACT with the ENDF/B-VII.1 51-Group Library through Benchmark Calculations for Critical Experiments
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Using neutron spin-echo spectroscopy, we studied the microscopic structural relaxation of a prototypical network ionic liquid ZnCl 2 at the structure factor primary peak and pre-peak. The results show that the relaxation at the primary peak is faster than the pre-peak and the activation energy is ≈33% higher. Stretched exponential relaxation is observed even at temperatures well above the melting point T m . Surprisingly, the stretching exponent shows a rapid increase upon cooling, especially at the primary peak, where it changes from stretched exponential to simple exponential on approaching T m . Furthermore, these results suggest that the appearance of glassy dynamics typical of the supercooled state even in the equilibrium liquid state of ZnCl 2 as well as the difference of activation energy at the two investigated length scales are related to the formation of network structure on cooling.
Nuclear fuel with extended enrichment (235U enrichment within 5-8 wt%) is one of the evolutionary changes that have been pursued in recent years by commercial light water reactor operators and fuel vendors to improve the fuel cycle economy and operation performance of a nuclear plant. This work assesses the performance of the Polaris/PARCS two-step approach in core physics modeling of the pressurized water reactor cores with extended enrichment fuel, referred to as “LEU+” in this report. A representative LEU+ core with a 24-month fuel cycle developed by Southern Nuclear Company (SNC) was modeled using this two-step approach. A representative LEU core with an 18-month fuel cycle was also modeled to provide a reference for the LEU+ core. As expected, significantly more burnable poison absorbers were used in the LEU+ core to accomodate its higher fuel enrichment. Nine different fuel assembly types were modeled using Polaris for each core to generate the assembly cross sections, which were then processed by GenPMAX to prepare the cross-section data for PARCS. The average specific powers of each fuel batch in each core were derived from VERA results and higher specific powers in fresh assemblies were found in the LEU core due to its less total uranium loading included in the VERA LEU model, given that the total core power was assumed to be the same for both cores. PARCS models were developed to simulate the steady-state operations of both cores. PARCS results on the LEU+ core were first compared with the VERA results for verification purpose; good agreements were seen in soluble boron and burnup distribution results, indicating that the Polaris/PARCS modeling and simulation were correctly implemented. Core physics parameters calculated by PARCS, at zero power physics tests, beginning of cycle (BOC), and end of cycle conditions (EOC), were compared between the LEU+ core and the LEU core, including soluble boron concentration, burnup distributions, assembly and pin power peaking factors, fuel temperature reactivity coefficients, moderator temperature and density reactivity coefficients, control rod worth, and shut down margin. The main differences in PARCS results between the LEU+ and the LEU cores are summarized below: 1)The critical boron concentrations were found to be much higher in the LEU+ core than in the LEU core (1582 vs. 1335 ppm for peak values). 2)Higher assembly radial power peaking factors (1.4 vs. 1.3 for peak values), 2D pin peaking factors (1.53 vs. 1.42 for peak values), and 3D pin peaking factors (1.89 vs. 1.81 for peak values) were found in the LEU+ core than in the LEU core. 3)Significantly higher reactivity coefficients of moderator temperature (and density) were found in LEU+ than in LEU.4)Significantly lower control rod worth at EOC were found in LEU+ than LEU for all but one control banks.5)Significantly lower shut down margins were found in the LEU+ core than in the LEU core. The Polaris/GenPMAX/PARCS code suite was found to be capable of modeling the LEU+ PWR core for steady-state operations and no unexpected results in core physics parameters were observed, in spite of that a) several bugs in PARCS were identified and workarounds were used; b) several features were found lacking in the current version of PARCS that would be useful for core modeling. A list of requests for bug fixes and feature upgrades for PARCS originated from this work were transmitted to the code developers. The assessments on the performance of the Polaris/PARCS two-step approach in core modeling for boiling water reactor with LEU+ fuel is ongoing.
This paper describes conceptual design of a 112 MW thermal (50 MW electric) Fast Modular Reactor (FMR) system operating at 7 MPa with inlet/outlet temperatures of 509/800 ºC. The reactor system includes the reactor core, fuel assemblies, fuel rods, reactor internals, reflector, neutron control system, flow control, and structural components. The nuclear design and analysis were conducted to search for a baseline core with a cycle length greater than 8-year and the power peaking factor less than 1.5 by adjusting the fuel assembly, reactor core, and reflector configurations. The neutronics calculations of the baseline core showed that a refueling interval of 9-year is achieved with a total peaking factor of 1.47 and a fuel rod average linear power of 3.6 kW/m.
This paper describes conceptual design of a 112 MW thermal (50 MW electric) Fast Modular Reactor (FMR) system operating at 7 MPa with inlet/outlet temperatures of 509/800 Celsius degrees. The reactor system includes the reactor core, fuel assemblies, fuel rods, reactor internals, reflector, neutron control system, flow control, and structural components. The nuclear design and analysis were conducted to search for a baseline core with a cycle length greater than 8-year and the power peaking factor less than 1.5 by adjusting the fuel assembly, reactor core, and reflector configurations. The neutron transport calculations of the baseline core showed that a refueling interval of 9-year is achieved with a total peaking factor of 1.47 and a fuel rod average linear power of 3.6 kW/m. (authors)
Two homogenized models of the MYRRHA (Multi-purpose hybrid Research Reactor for High-tech Applications) reactor in its critical configuration are presented in this study. Both models will be verified by comparing the main reactor parameters (K{sub eff}, β{sub eff} and Λ) with those of the heterogeneous model and a sensitivity analysis will be performed. A sensitivity analysis of some important parameters related to reactivity, namely the vacuum coefficient, the Doppler and the power peaking factors, is also presented. The study of the vacuum coefficient has been carried out by analyzing two possible scenarios: 50% vacuum in the coolant (lead bismuth eutectic) and 100%, obtaining a list of the 10 reactions and nuclides that most affect the k{sub eff} value. With respect to the Doppler coefficient, it will be seen that a 300 degree increase in fuel temperature (starting at 800 K) results in a reduction of the total reactivity of the system by 134 pcm.. Finally, the map of the power peaking factors will be shown, as well as a sensitivity analysis. For making the sensitivity calculations Serpent-2 reactor physics Monte Carlo code has been used and JEFF-3.3 library has been chosen in order to continue with the OECD/NEA WPEC SG46 work evaluating the nuclear data of this library. (authors)
The nuclear industry is pursuing approval of an increase in the length of the pressurized water reactor (PWR) cycle from 18 months to 24 months to reduce reactor downtime and enhance the economic competitiveness of nuclear energy. Such an increase in reactor cycle length will require that the maximum rod average burnup exceeds the current regulatory limit of 62 GWd/MTU, and it could peak at approximately 75 GWd/MTU, posing potential reactor safety and performance concerns. One such concern is that fuel fragmentation, relocation, and dispersal (FFRD) could occur during a severe loss-of coolant accident (LOCA) in which a fuel rod balloons and bursts, and pulverized fuel fragments are dispersed throughout the reactor’s primary coolant system. Previous analyses have identified which reactor operating conditions leave the core more susceptible to FFRD and have shown that FFRD susceptibility is strongly linked to fuel rod burnup and linear heat rate (LHR) history. The work described in this report uses an optimization strategy known as parallel simulated annealing (PSA) and a coarse mesh Purdue Advanced Reactor Core Simulator (PARCS) reactor physics model to develop two core fuel loading patterns, each with a different optimization objective. One core optimization maximized the core’s cycle length while still respecting regulatory limits on the radial peaking factor and soluble boron concentration with a peak rod average burnup of 75 GWd/MTU. The second optimization was aimed at minimizing FFRD susceptibility while still targeting a 24-month cycle length and respecting regulatory limits. PARCS model predictions were verified using the high-fidelity Virtual Environment for Reactor Applications (VERA). The two core designs were compared to highlight core design strategies to minimize FFRD susceptibility and to maximize economic viability.
New fuel design and development currently requires 20 to 25 years to be qualified for use by the nuclear power industry. The thermal-hydraulics community has taken advantage of scaling theory to design reduced scale experiments that correctly preserve dominant key phenomena while quantifying distorted phenomena. These techniques can be leveraged in the design and analysis of fuel performance experiments to help reduce the timeline associated with fuel design and development. This study uses the Dynamical System Scaling (DSS) method to analyze cladding temperature data from the recent SETH-C experiment in the TREAT facility and accompanying BISON simulations to assess dynamic distortions occurring throughout the fast power excursion transient. The DSS analysis revealed that on the cool down from peak cladding temperature that the fuel radial power profile is the most sensitive modeling parameter with a heterogenous radial peaking factor corresponding to the lowest distortion compared to a uniform energy deposition. For the heat up to peak cladding temperature the heterogeneous radial power profile corresponded to the shortest process action. Finally, for the heat up to peak cladding temperature, the gap conductance model sensitivity was quantified using process action and shows that the default Light Water Reactor gap conductance model corresponded to the longest process action.
Molybdenum (Z = 42) density spatial profiles have been obtained for Alcator C-Mod L-, I-, and H-mode plasmas from x-ray observations of the neon-like Mo 32+ 4D transition. Profile shapes range from flat to modestly peaked in the core. The observed peaking factors for molybdenum are in good agreement with those calculated by TGLF across all confinement regimes, and there is little dependence seen with core density gradient, sawtooth amplitude, or rotation velocity. Predictions for tungsten peaking in SPARC are presented.
Infrared thermography of the first wall in DIII-D is used to show the toroidal peaking of plasma radiation during mitigated disruptions with shattered pellet injection (SPI) and with massive gas injection (MGI). During MGI, the radiation peak location is shown to be due to the n = 1 magnetohydrodynamic (MHD) activity and continuously variable in toroidal phase based on externally applied error fields, consistent with previous experimental results. Furthermore, the measured toroidal peaking factor (TPF) is in agreement with the previous estimates based on radiometry, with a total TPF of 1.07 ± 0.05 when integrated over the entire duration of the disruption, and a value of 1.3 ± 0.1 during the thermal quench (TQ). For SPI, the location of the radiation peak is found to be determined predominantly by the particle source location, and thus can be varied on DIII-D by injecting from either of two toroidally separated injectors. 3D non-linear MHD simulations with the NIMROD code support this picture. The localized SPI particle source leads to higher peaking than for MGI, with a total TPF of 1.3 ± 0.1 over the full disruption and a TQ TPF of 1.9 +0.5/−0.3, a value also consistent with the NIMROD simulations. This TQ value is close to or potentially exceeding the allowable limit for the previously planned Be first wall in ITER, pointing to the importance of radiation asymmetries depending on the choice of first wall material.
Tungsten peaking is predicted in the core of a WEST plasma with total-f gyrokinetic simulations, including both collisional and turbulent transport. This prediction is validated with a synthetic diagnostic of the bolometry. Although nitrogen impurities are shown to reduce the neoclassical peaking of tungsten on-axis, the overall tungsten peaking increases when nitrogen impurities are present, as they reduce the turbulence screening off-axis. This finding is important for the plasma current ramp-up phase of ITER, where light impurities seeding will be desirable to achieve low temperatures at the plasma-facing components and reduce tungsten sputtering. It provides further argument for applying early ECRH heating to maintain margins on the core power balance. The neoclassical peaking factor is cross-verified between XGC and FACIT. The heat flux at separatrix and the heat load width are modeled by XGC and compared to WEST data.
The ability to deploy new nuclear fuels for current or future reactor concepts requires carefully designed experiments to generate data to support fuel qualification. Ideally these experiments would include state of-the-art sensing to maximize the amount of in situ data that can be collected during operation. Furthermore, advanced reactor systems can take advantage of integrated in-core sensing technologies to maximize fuel utilization, reduce unnecessary conservativism in design margins, and improve operator’s understanding of limiting peaking factors. Before any novel sensing technologies can be readily adopted for nuclear applications, they must first demonstrate acceptable performance in test reactors. This report summarizes the preliminary design and analysis of the most highly instrumented irradiation experiment ever performed in the removable beryllium (RB) positions of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory (ORNL). The Wireless Instrumented RB Experiment 2021 (WIRE-21) will test a wide range of sensors including wireless sensors being developed by Westinghouse Electric Company (WEC) that could provide in situ measurements of peak fuel temperatures and fuel rod pressurization due to fission gas release. The ability to wirelessly transmit a signal through the fuel rod’s cladding is critical to improving fuel monitoring capabilities without requiring signal penetrations through the cladding pressure boundary, which would significantly impact fuel fabrication, handling, and operation. Other sensors that will be tested in WIRE-21 include an array of thermocouples, self-powered neutron detectors (SPNDs), and spatially distributed fiber-optic temperature sensors. More generally, WIRE-21 will establish a flexible irradiation vehicle design to allow accelerated, economical testing of advanced sensor technologies while leveraging the extremely high neutron flux that is available in HFIR. This report summarizes the mechanical design for WIRE-21, the experimental test matrix, initial neutronic and thermal design analyses, and the active monitoring and control system enhancements necessary to support testing of advanced sensor technologies. The containment for WIRE-21 is similar to previous RB irradiation vehicles but includes a few modifications, most notably the use of integrated compression seals to pass a larger number of sensor leads through the experiment’s pressure boundary. In addition to the sensor leads, inert gas lines are passed into the experiment to enable active temperature control and the ability to pneumatically actuate a bellows-driven pressure sensor. WIRE-21 is targeting temperatures (300–350°C) and neutron fluence levels (~10 22 n/cm 2 ) relevant to light water reactors (LWRs), but the flexible design of the experiment vehicle allows much higher operating temperatures (>1,100°C). Neutronic calculations determine the neutron flux conditions as well as the nuclear heating within the experiments. These results are used as inputs to detailed thermal finite element calculations, which are required to evaluate the complex, three-dimensional heat transfer that occurs within WEC’s wireless sensor enclosures. Initial results show that the temperatures of the sensors’ enclosures and the metal bellows can be operated near the temperature range of LWR coolants and cladding while simultaneously increasing the temperature of a surrogate fuel material to values in the range of 800–1200°C to simulate centerline fuel temperatures during LWR operation.
Nuclear fuel vendors around the world are pursuing approaches to sustain the existing nuclear reactor fleet consisting primarily of pressurized-water reactors (PWRs) and boiling-water reactors (BWRs). To support the industry's efforts, advanced modeling and simulation tools need to be capable of analyzing both legacy reactor concepts. BWR fuel rods are significantly different than those used in PWRs, which can affect fuel performance analysis. BWR fuel rods include: (1) an extensive use of gadolinia dopant as a burnable absorber, (2) an axial variation in fuel enrichment and gadolinia content, (3) the inclusion of a liner on the inner cladding surface to mitigate the impact of pellet-clad mechanical interaction (which impacts hydrogen and hydride distribution), (4) a lower initial fill gas pressure, (5) bottom-entry control rods, and (6) a lower coolant pressure that results in the two-phase flow boiling phenomenon. Although the primary focus of BISON has been in the area of PWR and advanced reactor fuel analyses, this paper presents the developments in BISON to support BWR fuel performance analysis. An overview of the models that account for the effects of gadolinia is highlighted. Internal mesh generation capabilities to include a liner is presented. Normal operating and transient (reactivity insertion accident) demonstration problems are presented to illustrate the impact of gadolinia, the hydrogen and hydride evolution due to the presence of the liner, and BISON's ability to simulate axially varying enrichments and dopant concentration. Bottom-entry control effects are captured by the axial power peaking factors present in the demonstration cases. Comparisons to integral experiments from the Halden IFA-681 experiments are discussed as initial validation. Reasonable comparisons are obtained for fuel centerline temperature and rod internal pressure as a function of time. In conclusion, simulations of additional experiments containing Gd 2 O 3 -bearing fuel are necessary to completely validate the code for BWR applications.
The structural relaxation in water solutions of kosmotrope (structure maker) and chaotrope (structure breaker) salts, namely sodium chloride, potassium chloride, and cesium chloride, were studied through quasielastic neutron scattering measurements. We found that the collective dynamics relaxation time at the structure factor peak obtained using heavy water solutions shows a distinctively different behavior in the kosmotrope as opposed to the chaotrope solutions, increasing with the salt concentration in the former and decreasing in the latter. In both cases the trends are proportional to the concentration dependence of the relative viscosity of the solutions. Furthermore, these results indicate that kosmotropes and chaotropes influence the solution’s viscosity by impacting in opposite ways the hydrogen bond network of water, strengthening it in one case and softening it in the other.
Donor-acceptor (D-A) conjugated polymers are high-performance organic electronic materials that exhibit complex aggregation behavior. Understanding the solution state conformation and aggregation of conjugated polymers is crucial for controlling morphology during thin-film deposition and the subsequent electronic performance. However, a precise multiscale structure of solution state aggregates is lacking. Here, we present an in-depth small-angle X-ray scattering (SAXS) analysis of the solution state structure of an isoindigo-bithiophene-based D-A polymer (PII-2T) in chlorobenzene and decane as our primary system. Modeling the system as a combination of hierarchical fibrillar aggregates mixed with dispersed polymers, we extract information about conformation and multiscale aggregation and also clarify the physical origin of features often observed but unaddressed or misinterpreted in small-angle scattering patterns of conjugated polymers. The persistence length of the D-A polymer extracted from SAXS agrees well with a theoretical model based on the dihedral potentials. Additionally, we show that the broad high q structure factor peak seen in scattering profiles can be attributed to lamellar stacking occurring within the fibril aggregates and that the low q aggregate scattering is strongly influenced by the polymer molecular weight. Overall, the SAXS profiles of D-A polymers in general exhibit a sensitive dependence on the co-existence of fibrillar aggregate and dispersed polymer chain populations. We corroborate our findings from SAXS with electron microscopy of freeze-dried samples for direct imaging of fibrillar aggregates. Finally, we demonstrate the generality of our approach by fitting the scattering profiles of a variety of D-A polymers based on thieno-isoindigo (PTII-2T), diketopyrrolopyrrole (DPP2T-TT, DPP-BTZ, PDPP2FT-C- 16 ), naphthalenediimide (P(NDI2OD-T2)), and a conjugated block copolymer P3HT-b-DPPT-T. The results presented here establish a picture of the D-A polymer solution state structure and provide a general method of interpreting and analyzing their scattering profiles.
The authors developed an artificial intelligence (AI)-based algorithm for the design and optimization of a nuclear reactor core based on a flexible geometry and demonstrated a 3× improvement in the selected performance metric: temperature peaking factor. The rapid development of advanced, and specifically, additive manufacturing (3-D printing) and its introduction into advanced nuclear core design through the Transformational Challenge Reactor program have presented the opportunity to explore the arbitrary geometry design of nuclear-heated structures. The primary challenge is that the arbitrary geometry design space is vast and requires the computational evaluation of many candidate designs, and the multiphysics simulation of nuclear systems is very time-intensive. Therefore, the authors developed a machine learning-based multiphysics emulator and evaluated thousands of candidate geometries on Summit, Oak Ridge National Laboratory’s leadership class supercomputer. The results presented in this work demonstrate temperature distribution smoothing in a nuclear reactor core through the manipulation of the geometry, which is traditionally achieved in light water reactors through variable assembly loading in the axial direction and fuel shuffling during refueling in the radial direction. The conclusions discuss the future implications for nuclear systems design with arbitrary geometry and the potential for AI-based autonomous design algorithms.
In a previous paper, a Self-Organizing Map had proven to be able to identify the regions of the plasma operative space characterizing the pre-disruptive phase at JET without relying on any a priori information. One of the strengths of this disruption predictor lies in its inherent self-organization capability. The Self-Organizing Map discovers non-trivial relationships and captures the complicated interplay of device diagnostics on the internal plasma states directly from the experimental data. Moreover, the provided model allows the visualization of high-dimensional plasma parameters and facilitates easy interrogation of the model to understand the reasons behind its correlations. In this paper, an additional step is taken towards the interpretability of models for predicting disruptions by training a Decision Tree to classify the plasma states according to the interpretation provided by the Self-Organizing Map (stable or at high risk of disruptions). The Decision tree provides a set of rules which describe the transition of the plasma towards the pre-disruptive phase as visualized in the Self-Organizing Map. The obtained rules for the database explored in the study identify four regions in the map, two of which are at risk of disruption. These regions correspond to partitions of a 3D space based on the peaking factors of the core and divertor radiation, as well as the Locked Mode. The agreement between the Self-Organizing Map answers and the rules supplied by the Decision Tree is confirmed by the comparison of the performance exhibited by the two models in the prediction of disruptions.
Here we present techniques for evaluating the radiated power in JET disruptions. Disrupting plasmas are shown to have non-axisymmetric radiation profiles, motivating the re-evaluation of the standard techniques for calculating the total radiated power at JET using bolometry. Four single-channel bolometers at different toroidal locations are exploited to quantify the radiation asymmetry. Toroidal radiation peaking factors integrated over the entire disruption of up to 1.5 have been observed when varying the quantity of neon in pellets used in disruptions mitigated by shattered pellet injection. Using synthetic bolometer diagnostics developed with the Cherab spectroscopy modeling framework, we can estimate the systematic error on total power calculations for relevant radiation profiles and improve estimates of the total radiated power. We show that the component of the systematic error on the total power due to the poloidal radiation profile can be reduced from 70% to 10% with suitable assumptions about the structure of the poloidal profile.