Sensitivity Studies on Fuel Burnup and Plutonium Buildup in Pebble Bed Reactors for Fuel Parameter Perturbations
Explore the source record for details and available documents.
Hartanto, Donny↗
SEARCH · Engineering Papers
Results for “Sensitivity Study”
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.
Sensitivity Study of Epithermal Resonance Up-scattering Using SCALE/Polaris-PARCS Code Procedure for LWR Physics Analysis
Explore the source record for details and available documents.
Kim, Kang Seog [ORNL] (ORCID:0000000186229321)↗
Sensitivity Study of Decomposing Carbon Fiber Epoxy Composites Model
Explore the source record for details and available documents.
Scott, Sarah Nicole↗
Thermal-Hydraulics Modeling and Simulations of Hot Pool Using the SAS-CFD Coupled Code
The main goal of this activity is to test the dynamic coupling of the SAS4A/SASSYS-1 (SAS) and CFD models, using a recently patched version of the SAS code intended to address an undocumented limitation that hindered the Versatile Test Reactor (VTR) simulation efforts in FY21. As described in previous VTR calculation reports, the undocumented limitation in SAS v5.4 does not allow the user to activate the CFD coupling option during restart calculations. Since the analysts were unaware of this limitation, prior SAS-CFD simulation results for the protected station blackout (PSBO) transient were erroneous. Root-cause analysis was performed to determine the cause of this undocumented limitation in SAS v5.4, the SAS software was updated in a new patch, and the SAS-CFD simulations were repeated with this patched software. The results of the SAS-CFD simulations documented in this report show that the software patch does address the cited issue, and that the patched software indeed supports the activation of the CFD coupling model in restart calculations. The SAS development team will determine the schedule for implementing the patch in an official software release. This report documents updated SAS-CFD simulations of the PSBO transient response in the VTR. The hot pool is modeled with the CFD code STAR-CCM+, which is coupled at the flow boundaries to the SAS model of the primary heat transport system. SAS computes the mass flow rate and temperature at each core subassembly outlet, the thermal insulation cavity bypass, and the IHX inlet windows. CFD in turn computes the absolute pressure and temperature at each of these boundaries. The SAS code will ignore the temperature data at flow boundaries where flow is directed into the hot pool, i.e., at the core subassembly outlets unless flow reversal occurs. Similarly, CFD will ignore temperature data at boundaries where the flow is directed out of the hot pool, i.e., at the IHX inlets except under flow reversal. The focus of this work is to ensure that the SAS software patch addresses the undocumented limitation described in prior VTR calculation reports, rather than the accurate assessment of thermal stratification in the VTR during protected transients. This motivates the development of a new, simplified CFD model with a coarser mesh to accelerate the testing process. The updated model, and simplifying assumptions, are documented in this report. In future work, the thermal stratification assessment should be performed in more detail. The simplified CFD model can be improved by performing grid convergence studies sensitivity studies of turbulence parameters (e.g., Prandtl number, turbulence production and dissipation parameters) on temperature distributions and thermal stratification.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Radiological Releases from Novel Fuel Forms in Advanced Reactors During Severe Accidents for Consequence Analyses
Various advanced reactor developers are exploring the potential for reductions in the size of physical security forces and emergency planning zones. These reductions are based on robust fuel forms and inherently safe reactor designs. However, such reductions in physical protection measures could increase the risk of sabotage. To assess the possibility of reducing these measures, sabotage-induced radiological consequence analyses were carried out. These analyses considered accident scenarios that were beyond design basis accidents and overly conservative (Shah, 2025a; Shah, 2025b; Shah and Hartanto, 2026), yielding very large release fractions. These fractions, which can be used to evaluate physical protection and emergency planning requirements, have been crudely determined and applied as demonstrations for a sodium-cooled fast reactor (SFR) (Shah and Hartanto, 2025a), a high-temperature gas-cooled reactor (HTGR) (Shah and Hartanto, 2025b), a heat pipe–cooled reactor (HPR) (Shah and Hartanto, 2025c), and a molten salt–cooled reactor (MSR) (Shah et al., 2026). A Sandia National Laboratories (SNL) team used MELCOR—a fully integrated severe accident analysis code—to demonstrate the code’s capability to analyze advanced (i.e., not light water–cooled) reactors (including a fluoride salt–cooled high-temperature reactor [FHR]) and calculate radiological releases to the environment during severe accidents (Wagner et al., 2022a, 2022b, 2022c, 2023a, and 2023b). Although the analyses were carried out to demonstrate MELCOR’s growing capability, the release source terms were estimated for advanced reactors, providing valuable insights into the accident progression and radiological releases. These findings from prior SNL studies, including estimated source terms and related sensitivity studies, were leveraged to derive source terms for postulated sabotage-induced accidents. Insights from these sensitivity studies informed the scaling of SNL’s estimated source terms for the defined accident scenarios. The derived release fractions for the severe accident scenarios for the respective reactor designs can be used to perform more nuanced dose consequence analyses to evaluate the reactors’ physical protection and emergency planning zone requirements. These analyses are in accordance with the risk-informed, performance-based approach proposed under 10 CFR Part 53. This study builds on the prior source term analyses and associated sensitivity studies by SNL to derive time-dependent and design-informed release fractions. Section 2 describes the diverse advanced reactor designs analyzed by the SNL team. Section 3 discusses the severe accident analyses, the release fractions calculated, and the limitations and assumptions of the demonstration project. Section 4 presents the release percentages derived for the hypothetical sabotage-induced severe accidents at the advanced reactors. Section 5 summarizes the study’s findings and conclusions.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Understanding asymmetries using integrated simulations of capsule implosions in low gas-fill hohlraums at the National Ignition Facility
Current capsule implosions at the National Ignition Facility (NIF) using high-density-carbon ablators and laser energies close to 2 MJ have shown neutron yields in excess of 50 kJ. Improving on this performance requires understanding of the different degradation mechanisms. For many NIF implosions, nuclear diagnostic signatures have inferred residual hot-spot velocities that correlate with fuel areal density variations consistent with a low-order mode-1 asymmetry. A current working hypothesis attributes these asymmetries to a combination of beam to beam variations in the laser delivery and possibly coupling to target features, such as the diagnostic holes needed for x-ray imaging. Recently, a new source has been identified, thickness variation in the ablator shell. To gain better understanding and eventually mitigate the causes of <ρR> asymmetries, 3D integrated simulations using the actual delivered laser powers are needed. To capture the effect of the diagnostic holes (DHs) using direct numerical simulation would require significantly large computational resources. Instead, our 3D simulations make use of a subgrid model developed using highly-resolved 2D simulations that include several details of the DH engineering complexity. Simulations of NIF shots using hohlraums without DHs, to isolate the effect of beam-to-beam variations, reproduce fairly well the observed nuclear diagnostic signatures. Similarly, reasonable agreement between data and simulations is also obtained in the presence of diagnostic holes. To account for the remaining discrepancies a sensitivity study of ablator thickness variation showed that 1% thickness asymmetries are comparable in effect to 1% peak drive mode-1 asymmetries. Furthermore, this study identified sensitivity to variations in the imbedded doped layer (needed to shield the DT ice from the high energy x-rays generated in the hohlraum) thickness even when the inner and outer surface of the ablator are perfectly spherical.
radiation hydrodynamics↗
304L Can Crush Validation Studies
Accurate prediction of ductile behavior of structural alloys up to and including failure is essential in component or system failure assessment, which is necessary for nuclear weapons alteration and life extensions programs of Sandia National Laboratories. Modeling such behavior requires computational capabilities to robustly capture strong nonlinearities (geometric and material), rate- dependent and temperature-dependent properties, and ductile failure mechanisms. This study's objective is to validate numerical simulations of a high-deformation crush of a stainless steel can. The process consists of identifying a suitable can geometry and loading conditions, conducting the laboratory testing, developing a high-quality Sierra/SM simulation, and then drawing comparisons between model and measurement to assess the fitness of the simulation in regards to material model (plasticity), finite element model construction, and failure model. Following previous material model calibration, a J 2 plasticity model with a microstructural BCJ failure model is employed to model the test specimen made of 304L stainless steel. Simulated results are verified and validated through mesh and mass-scaling convergence studies, parameter sensitivity studies, and a comparison to experimental data. The converged mesh and degree of mass-scaling are the mesh discretization with 140,372 elements, and a mass scaling with a target time increment of 1.0e-6 seconds and time step scale factor of 0.5, respectively. Results from the coupled thermal-mechanical explicit dynamic analysis are comparable to the experimental data. Simulated global force vs displacement (F/D) response predicts key points such as yield, ultimate, and kinks of the experimental F/D response. Furthermore, the final deformed shape of the can and field data predicted from the analysis are similar to that of the deformed can, as measured by 3D optical CMM scans and DIC data from the experiment.
36 MATERIALS SCIENCE↗
Overview of Base Model in Parametric Studies Specific to Performance of U-Mo Plates
This paper provides an overview of the base model specifically developed to perform parametric sensitivity studies on the U-10Mo monolithic fuel system. U-Mo monolithic fuels are being considered for the conversion of test reactors into high-performance research reactors that operate using proliferation-resistant, low-enriched uranium (LEU) fuels. These plate-type fuels contain a high-density, low-enrichment fuel sandwiched between zirconium diffusion barriers and encapsulated in aluminum claddings. All U.S. high-performance research reactors have released the designs of their LEU monolithic fuel reactor cores. These designs include nearly 50 distinct fuel plate geometries with different operational parameters. Consequently, a single generic plate geometry representing all the extreme points in this design matrix is unrealistic. To evaluate the performance for various parameters, a set of sensitivity studies was performed. These studies considered various input parameters (i.e., geometric, operational, and material property-related). The results revealed valuable information about plate performance and the sensitivity of this performance to various modeling inputs. To establish a reference state for comparing these result, base model featuring representative irradiation conditions was developed. To capture in-reactor behavior accurately, incorporation of representative constitutive models capable of evolving properties with respect to temperature, irradiation time, and burnup was needed. The behavioral models considered burnup-dependent properties, swelling, creep, and degradation. This paper introduces the base model created for the parametric sensitivity studies. The detailed description of the procedure includes the model geometry, model discretization, thermo-mechanical coupling, material properties and behavioral models. This paper also provides selected results and assesses the performance of the base model.
42 ENGINEERING↗
Initial ASME code rule analysis on Simple and Full assessment
ASME & ASTM codes present simple and full assessment methodologies to qualify nuclear grade graphite components. This report explains the theory behind the statistical portion of the codes and describes the issues discussed in a 2020 workshop in the simple assessment, as well as the full assessment. The methodology for the full assessment was developed by Hindley []. He used a validation methodology to match the experimental average failure load of test specimens from several geometries to the calculated 50% POF load to tune the grouping criteria parameters based on an RMSE penalty function. The grouping criteria consists of a minimum volume to satisfy the weakest link theory of the Weibull distribution and a minimum stress range parameter. Hindley’s work found a minimum link volume of 10 times the grain size and a stress range parameter of 7% [satisfy the validation requirements. Several studies found that the link volume of 10 times the grain size was not satisfactory for the volume grouping criteria for graphite grades with fine grains. ASME 2021 code adopted a new volume grouping criteria based on fracture toughness for calculation of the process zone volume. However, errors were found in the ASME 2021 code process zone volume equation and the question is now open as to what volume grouping criteria is satisfactory. This report presents results from sensitivity studies that were done to evaluate the volume grouping criteria effect on the full assessment POF. This report also presents results from sensitivity studies on other aspects of the code that affect the POF, the mesh size and the choice of the Weibull threshold parameter. All three sensitivity studies: the volume grouping criteria, the mesh size, and the threshold parameter choice are found to affect the component qualification decision for components of structural reliability class 1 (SRC-1). The sensitivity analysis results are presented in the body of the report for NBG-18, which is the same grade of graphite as Hindley used in his thesis. Results from NBG-17, IG-110, 2114, and PCEA are presented in the Appendix. For NBG-18, it was found that increasing the threshold increases the POF, finer meshes result in higher POFs, and that smaller link volume requirements lead to higher POFs (which is inconsistent with Hindley’s findings). This report proposes a new validation study similar to Hindley’s []. The purpose of the new validation study is to tune the threshold, mesh size and grouping criteria based on results from multiple grades of graphite, including finer grain graphites T220 (1-3 microns) and NG-CT-50 (5 microns). The test specimens will be broken at 3 labs, such that we can measure between lab measurement variability, with at least 5 test specimens per lab such that the within lab uncertainty can also be measured. Hindley’s validation method tried matching the average experimental load with the median. The future study will attempt to match failures at lower percentiles, which would be closer to the SRC POF limit. We are currently looking into what geometries make sense and the best method to measure multi-axial stress states. These design decisions are still being set.
36 MATERIALS SCIENCE↗
“Experimental investigation of the governing parameters of atmospheric ice nucleation using field-collected and laboratory generated aerosol particles and its application in cloud resolving models” (Final Report)
The objective of this research project is to improve our understanding of the role of aerosol particles acting as ice-nucleating particles (INPs) which in turn define the mixed-phase and cirrus cloud radiative properties and thus climate. The focus is placed on how the physicochemical particle population properties determine the particles’ ability to initiate ice nucleation. This research project combined micro-spectroscopic particle analysis, experimental ice nucleation studies, and model sensitivity studies to advance our predictive capability of the formation of mixed-phase and cirrus clouds. These project activities have led to new ice nucleation data from laboratory generated and ambient (authentic) aerosol particles furthering process-level understanding, insights in the role of organic aerosol in ice formation, and advancements in the interpretation and parameterization of ice nucleation.
54 ENVIRONMENTAL SCIENCES↗
New Cross Section Measurements and Evaluations of Zirconium Isotopes [Slides]
This presentation gives a look to new cross section measurements and evaluations of Zirconium (Zr) isotopes. This lecture includes highlights of accomplishments made in 2022. Additionally, it explains the motivation for investigation, new measurements, and re-evaluation of the Zr isotopes. This presentation concludes with a preliminary evaluation result from RPI 90 Zr evaluation and preliminary results from Eigenvalue sensitivity study on Zr sensitive critical benchmark experiments.
07 ISOTOPE AND RADIATION SOURCES↗
Development of a U-19Pu-10Zr fuel performance benchmark case based on the IFR-1 experiment
Metallic nuclear fuels are subject to research and development for use in advanced reactors. Robust, accurate metallic fuel performance models are important for the design, safety analysis, and licensing of these reactors. However, metallic fuel performance models require additional development; they are not as mature as UO 2 fuel performance models. A benchmark case based on the IFR-1 experiment was developed to better gauge the accuracy of existing models, identify models for high-priority development, and potentially quantify any future improvements made by further model development.This work collected publicly available information on the IFR-1 experiment and used it to develop the benchmark case. Fuel behavior during the IFR-1 irradiation was simulated by using the fuel performance code BISON, and the predicted results were compared with postirradiation examination data from the IFR-1 experiment. Furthermore, a sensitivity study and tuning studies were performed as a preliminary investigation into the causes of inaccurate temperature and dimensional change predictions.The benchmark predicted reasonably accurate values for the burnup and fission gas release. There was error in the predicted temperatures, which could be explained by uncertainty in the input parameters and legacy temperatures. BISON over-predicted dimensional changes in the fuel and cladding. The sensitivity study showed that the dimensional changes were most sensitive to the fuel swelling anisotropy and the cladding void swelling model. Future benchmark and model development should focus on cladding swelling behaviors to improve dimensional change predictions.
11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
Particle Swarm Optimization of Dynamic Load Model Parameters in Large Systems
This paper considers two dynamic load models that are widely used in industry to account for induction motor behavior: CMLD and CLOD. These models must be parametrized for the specific utility system in a general way so that they can be used in planning studies and provide a conservative but realistic representation of load behavior. This study considers a measurement-based approach to tuning both models. The load modeling study compares the response of the tuned models to generic candidate models using historical events. This study considers one area-based subsystem to simplify the modeling approach and reduce the number of models required for simulations. Additionally, because dynamic load models often produce similar results for different sets of parameters, a sensitivity study was conducted to assess the parameter impacts on the voltage response. The sensitivity study covers the parameters that are tuned using event measurements. The process to estimate the parameters uses the particle-swarm optimization algorithm. Overall, the performance of the tuned model more accurately captures recovery voltage, delayed recovery, and settling voltage than its predecessor models while not being overly tuned so that it remains general for peak summer conditions.
dynamic load modeling↗
The role of electron-scale turbulence in the JET tokamak: experiments and modelling
Dedicated electron heat transport experiments have been carried out in L- and H-mode Deuterium plasmas of the JET-ILW tokamak to identify the amount of electron heat carried by electron-scale electron temperature gradient (ETG) modes. Ion cyclotron resonance heating at different positions has been used to probe the response of the electron temperature inverse gradient length R/L Te to changes in electron heat flux qe, while different amounts of neutral beam heating allowed to scan the ratio of ion to electron temperature T e /T i , which is a key parameter for the onset of ETGs. Results indicate a steepening of the normalized qe vs R/L Te curve above R/LT e ~ 8 for T e /T i ≤ 1, suggestive of the ETG onset. Ion-scale gyro-kinetic (GK) simulations match the ion heat flux and the low-R/L Te part of the qe curve, but do not reproduce such steepening at high R/L Te . Multi-scale GK simulations covering both ion and electron scales and including one impurity bundling light and heavy species indicate an ETG contribution only for R/L Te values larger than the experimental ones. Sensitivity studies of such result are difficult to achieve due to limitation in numerical resources. The quasi-linear TGLF model has been used for sensitivity studies. With the same bundled impurity as the GK multi-scale, TGLF shows the qe steepening at much larger R/L Te values than in experiment, but when using the real mix of light impurities neglecting the heavy impurities, TGLF gets closer to the experimental results. Profile simulations with TGLF including both light and heavy impurities show over-prediction of T e profiles and in some cases also of density, but good T i predictions, confirming issues with the model electron stiffness for these plasmas.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Integrated Steady-State System Package for Nuclear Thermal Propulsion Analysis Using Multi-Dimensional Thermal Hydraulics and Dimensionless Turbopump Treatment
Nuclear thermal propulsion is an evolving technology that can be utilized for long-distance space travel. This technology yields the advantage of a high thrust and specific impulse, but requires an examination of the potential design adjustments necessary to enhance its feasibility. The development of nuclear thermal propulsion requires a comprehensive understanding of the system-level behavior during transient and steady-state operation. This paper extends our previous research by including the proper handling of turbomachinery with multi-channel thermal hydraulic simulations only for steady-state solutions. The system-level approach presented here enables the treatment of the turbopump components through non-dimensional analysis that eliminates the assumption of constant efficiencies. All the other components within the system (e.g., reflector and core) can be discretized to multiple channels and layers, in which the full thermal hydraulic solution is established. The approach chosen here enables the realistic modeling of the propellant flow within the expander cycle by capturing the pressure losses, mass flow rate splits, and enthalpy gain for various operational conditions. The verification of the package is completed through point comparisons of previous investigations into similar system designs. Furthermore, sensitivity studies are used to benchmark the capabilities of the package and investigate solution variations due to the perturbation of operational conditions and regimes. The sensitivity studies performed here are important to capture variation in flow characteristics (e.g., temperature, pressure, mass flow rates) for different design objectives such as the thrust and specific impulse. This work demonstrates that system-level simulations lacking multi-channel capability and proper turbomachinery treatment may yield higher uncertainties in understanding the engine’s response and characteristics to changing various requirements. This is extremely important when screening the design space of such propulsion systems and when transient simulations are required.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Cost Drivers for a Tokamak-Based Compact Pilot Plant
A comprehensive systems code that includes a range of physics and engineering considerations along with a simplified costing model has been utilized to evaluate the primary cost drivers for a compact tokamak pilot plant. The systems code has been benchmarked against several tokamak reactor designs and is utilized with sophisticated optimization algorithms to develop optimal solutions for a set of user-specified assumptions and design constraints. Here, in contrast to previous models that have focused on the cost of electricity as the key cost metric, this study uses the estimated capital cost of the facility. The analysis suggests that a pilot plant with the following features may offer potential for a cost-attractive pilot plant: A ~ 3, H 98y2 > 1.5, P net = 200 MW, 𝜏 𝑝𝑢𝑙𝑠𝑒 ~ 1 to 2 h utilizing rare-earth barium copper oxide (REBCO) magnet technology and a plug-bucked central solenoid/toroidal field (CS/TF) magnet support structure. While REBCO magnets offer some advantages relative to Nb 3 Sn magnets in all cases, the most gain is obtained when combined with the plug-bucked CS/TF bucking solution. Pulsed operation reduces capital cost requirements relative to steady-state operation, especially at low confinement. Cost sensitivity studies indicate that there are significant cost uncertainties associated with the achievable confinement quality, tritium breeding capability, attainable thermal efficiency, and achievable neutron wall loading, suggesting that these areas are the most critical areas in reducing the cost risk for a compact tokamak pilot plant. Further cost sensitivity studies indicate that the estimated cost is most sensitive to the underlying cost of the magnetic coils, providing further impetus to better establish cost-effective means for producing fusion magnets.
Compact pilot plant↗
Techno-economic comparison of sCO 2 cycles for particle-based CSP at design-point conditions
In this work, we compare the techno-economic performance of supercritical carbon dioxide power cycles integrated in a particle CSP system. We model four core cycle configurations: simple (with optional bypass), recompression (with optional bypass), partial cooling, and turbine split flow, which each demonstrate different benefits in a CSP system, such as high efficiency, low cost, or large HTF temperature differences. We parametrically sweep cycle design variables for each configuration. The set of power cycle performance results are then combined with a design point particle CSP system model which calculates the system specific cost. The simple cycle and turbine split flow cycles have the best performance in the baseline results, with system specific costs of 5,912 and 5,899 $\$$/kWe respectively. In addition to the baseline set of results, we also vary key parameters and costs in a sensitivity study. The cycle designs with the best system performance limit their efficiency to ~45 %, despite demonstrating higher maximum efficiencies, due to the rapid increase in cost of recuperation as efficiencies rise. The simple cycle has strong performance in the analysis and is on average only 1.4 % worse than the optimal configuration. Lowered turbine inlet temperatures from the sensitivity study improve performance by reducing the PHX and turbine cost. Decreasing the inlet temperature from 700 to 625°C results in an >8 % decrease in system specific cost. Future work should expand sensitivity analyses to colder turbine inlet temperatures and calculate system performance by simulating annual performance with off-design solar and cycle component models.
14 SOLAR ENERGY↗
Air Blast Meshing & Pressure Mapping [Slides]
Project Goals: Mesh sensitivity study – Determine at what mesh density results converge; Build a model that reflects theoretical solutions. Map experimental pressure data onto the surface of the cone – Track rigid body motion. Stresses in conic Body – Observe how the model reacts to stresses different than the ones tested in the mesh sensitivity study.
42 ENGINEERING↗