Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model discrepancy”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Mapping TpPa-1 covalent organic framework (COF) molecular interactions in mixed solvents via atomistic modeling and experimental study

Complex solvent environments continue to limit the widespread adoption of organic solvent nanofiltration (OSN) in many chemical industry applications. In this paper we employ a commercially available covalent organic framework (COF), TpPa-1, and force field models to molecularly map separation performance of TpPa-1 membrane in mixed solvents. To minimize time and length scale mismatch between atomistic modeling and experiments, solvent permeance was normalized with water in modeling and experimental results to enable direct comparison. Model outputs, such as organic solvent permeance and solute rejection rate, matched well with filtration results. Since the atomistic models assume that all mass transfer is via through-pore transport, the discrepancies between modeling and experimental results provide insights on the effect of linear polymer defects, adsorption and interstitial mass transfer on polycrystalline COF membrane performance. In sum, force field models can serve as digital twins of COF membranes to simulate separation processes while capturing the effects of COF structure, chemistry, and crystallinity on membrane performance in complex organic solvent environments. Finally, this approach will provide insight into future COF design and synthesis for persisting separation challenges.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2016

Interpreting space-based trends in carbon monoxide with multiple models

Here, we use a series of chemical transport model and chemistry climate model simulations to investigate the observed negative trends in MOPITT CO over several regions of the world, and to examine the consistency of time-dependent emission inventories with observations. We also found that simulations driven by the MACCity inventory, used for the Chemistry Climate Modeling Initiative (CCMI), reproduce the negative trends in the CO column observed by MOPITT for 2000–2010 over the eastern United States and Europe. However, the simulations have positive trends over eastern China, in contrast to the negative trends observed by MOPITT. The model bias in CO, after applying MOPITT averaging kernels, contributes to the model–observation discrepancy in the trend over eastern China. This demonstrates that biases in a model's average concentrations can influence the interpretation of the temporal trend compared to satellite observations. The total ozone column plays a role in determining the simulated tropospheric CO trends. A large positive anomaly in the simulated total ozone column in 2010 leads to a negative anomaly in OH and hence a positive anomaly in CO, contributing to the positive trend in simulated CO. Our results demonstrate that accurately simulating variability in the ozone column is important for simulating and interpreting trends in CO.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

High‐Frequency Intermittency in Observed and Model‐Simulated Precipitation

Abstract A newly reprocessed, bias‐corrected version of hourly satellite observations that provides global coverage of precipitation at high space‐time resolution is evaluated and compared with climate model simulations. Irregular subdaily fluctuations are the dominant component around the world, greater than variance of daily mean precipitation, and much greater than variance associated with the mean diurnal cycle of precipitation. Irregular subdaily fluctuations of precipitation are severely underestimated by models, even after taking into account the observational “error bars” implied by different space‐time resolutions. Variance of daily mean precipitation is less severely underestimated. Although mean diurnal cycle amplitudes vary among the models, this component is but a small part of total precipitation variance. Therefore, the total precipitation variance is significantly underestimated by models in general. Further exploration of model‐data discrepancies in precipitation at high‐time frequency may lead to new and useful climate model diagnostics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Enhanced production of 60 Fe in massive stars

Massive stars are a major source of chemical elements in the cosmos, ejecting freshly produced nuclei through winds and core-collapse supernova explosions into the interstellar medium. Among the material ejected, long-lived radioisotopes, such as 60 Fe (iron) and 26 Al (aluminum), offer unique signs of active nucleosynthesis in our galaxy. There is a long-standing discrepancy between the observed 60 Fe/ 26 Al ratio by γ-ray telescopes and predictions from supernova models. This discrepancy has been attributed to uncertainties in the nuclear reaction networks producing 60 Fe, and one reaction in particular, the neutron-capture on 59 Fe. Here we present experimental results that provide a strong constraint on this reaction. We use these results to show that the production of 60 Fe in massive stars is higher than previously thought, further increasing the discrepancy between observed and predicted 60 Fe/ 26 Al ratios. The persisting discrepancy can therefore not be attributed to nuclear uncertainties, and points to issues in massive-star models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Interplay of Laser-Plasma Interactions and Inertial Fusion Hydrodynamics

The effects of laser-plasma interactions (LPI) on the dynamics of inertial confinement fusion hohlraums are investigated in this work via a new approach that self-consistently couples reduced LPI models into radiation-hydrodynamics numerical codes. The interplay between hydrodynamics and LPI—specifically stimulated Raman scatter and crossed-beam energy transfer (CBET)—mostly occurs via momentum and energy deposition into Langmuir and ion acoustic waves. This spatially redistributes energy coupling to the target, which affects the background plasma conditions and thus, modifies laser propagation. In conclusion, this model shows reduced CBET and significant laser energy depletion by Langmuir waves, which reduce the discrepancy between modeling and data from hohlraum experiments on wall x-ray emission and capsule implosion shape.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Local structural variation with oxygen fugacity in Fe 2 SiO 4+x fayalitic iron silicate melts

Here, the structure of molten Fe 2 SiO 4+x has been studied using both high-energy X-ray diffraction and Fe K-edge X-ray absorption near-edge structure (XANES) spectroscopy, combined with aerodynamic levitation and laser beam heating. A wide range of Fe 3+ contents were accessed by varying the levitation and atmospheric gas composition. Diffraction measurements were made in the temperature ( T ) and oxygen partial pressure ranges 1624(21) < T < 2183(94) K (uncertainties in parentheses) and –5.6(3) < ΔFMQ < +2.8(5) log units (relative to the Fayalite-Magnetite-Quartz buffer). Iron K-edge XANES measurements covered the ranges 1557(33) < T < 1994(36) K and –2.1(3) < ΔFMQ < +4.4(3) log units. Fe 3+ contents, x = Fe 3+ /ΣFe, estimated directly from the pre-edge peaks of the XANES spectra varied between 0.15(1) and 0.40(2). While these agree in some cases with semi-empirical models, notable discrepancies are discussed in the context of the redox kinetics and the limitations in both the models and in the calibrations used to derive oxidation state from XANES spectra. XANES pre-edge peak areas imply average Fe–O coordination numbers, n FeO , close to 5 for all Fe 3+ /ΣFe. Diffraction measurements yielded values of 4.4(2) < n FeO < 4.7(1). There is limited evidence for a linear trend n FeO (x) = 4.46(3) + 0.4(1)x. Asymmetric Fe–O bond length distributions peak at around 1.96 Å and have a shoulder arising from longer interatomic distances. Mean r FeO lie close to 2.06 Å, consistent with n FeO close to 5. These observations suggest that Fe 2+ is less efficient at stabilizing tetrahedral Fe 3+ compared to large monovalent alkali cations. Comparison of in-situ XANES estimates of Fe 3+ /ΣFe in the melts to those of the quenched solids obtained from XANES as well as Mössbauer spectroscopy indicate rapid oxidation during cooling, enabled by stirring of the melt by the levitation gas flow. As such, the oxidation state of hot komatiitic and other highly fluid melts may not be retained, even during rapid cooling, as it is for cooler basaltic and more silicic magmas.

58 GEOSCIENCES↗
DOE OSTI2021

An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast—Part 2: Circulation, Boundary Layer, and Clouds

The Western North Atlantic Ocean (WNAO) is a complex land-ocean-atmosphere system that experiences a wide range of atmospheric phenomena, which in turn drive unique aerosol transport pathways, cloud morphologies, and boundary layer variability. This work, Part 2 of a 2-part paper series, provides an overview of the atmospheric circulation, boundary layer variability, three-dimensional cloud structure, and precipitation over the WNAO; the companion paper (Part 1) focused on chemical characterization of aerosol, gas, and wet deposition. Seasonal changes in the circulation and sea surface temperature explain a clear transition in cloud morphologies from small shallow cumulus clouds, convective clouds, and tropical storms in summer, to stratus/stratocumulus and multi-layer cloud systems associated with winter storms. Synoptic variability in cloud fields is estimated using satellite-based weather states, and the role of postfrontal conditions (cold-air outbreaks) in the development of stratiform clouds is further analyzed. Precipitation is persistent over the ocean, with a regional peak over the Gulf Stream path, where offshore sea surface temperature gradients are large and surface fluxes reach a regional peak. Satellite data show a clear annual cycle in cloud droplet number concentration along the coast in winter, and minima in summer, suggesting a marked annual cycle in aerosol-cloud interactions. The performance of four climate models are assessed against satellite cloud retrievals. The models qualitatively reproduce the annual cycle in cloud cover and liquid water path, but with large discrepancies across models, especially in the extra-tropics. The paper concludes with a summary of outstanding issues and recommendations for future work.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Poleward Migration of the Latitude of Maximum Tropical Cyclone Intensity—Forced or Natural?

Abstract Past studies have shown a significant observed poleward trend in the latitude at which tropical cyclones reach their lifetime maximum intensity (LMI), especially in the northwest Pacific basin. Given the brevity of the historical record, it remains difficult to separate the forced trend from internal variability of the climate system. A recently developed tropical cyclone downscaling model is used to downscale the Community Earth System Model, version 2 (CESM2), preindustrial control simulation. It is found that the observed trend in the latitude at which tropical cyclones reach their LMI in the northwest Pacific is very unlikely to be caused by internal variability. The same downscaling model is then used to downscale CESM2 simulations under historical forcing. The resulting trend distribution shows a significant poleward migration of tropical cyclone LMI even after regressing out both natural variability and the part of the forced warming pattern that projects onto natural variability. The results indicate that the observed poleward migration of the latitude at which tropical cyclones reach their LMI in the northwest Pacific basin is likely to be, at least in part, forced. However, the magnitude of the projected poleward trend in climate models can be significantly modulated by the simulated spatial pattern of ocean warming. This highlights how discrepancies between models and observations, with regard to projected changes to the equatorial zonal sea surface temperature gradient under anthropogenic forcing, can lead to large uncertainties in projected changes to the LMI latitude of tropical cyclones. Significance Statement Observations in the northwest Pacific basin show that the latitude at which tropical cyclones are at their most intense has been trending northward in the recent half century. These changes are important since tropical cyclones could bring hazardous weather to coastal areas that are poorly equipped to handle them. Here, we show that natural variations in Earth’s climate are very unlikely to explain the observed poleward trend in the latitude that tropical cyclone reach their maximum intensity. We find that it is much more likely that the observed trend is forced by human-related emissions, though the spatial pattern of warming in response to greenhouse emissions can have significant impacts on the magnitude of the trend.

Lin, Jonathan↗
DOE OSTI2026

A theoretical index for understanding distinct land relative humidity trends in observations, reanalyses, and models

Land surface relative humidity (RH) is a key variable in the coupled land–atmosphere system that profoundly influences terrestrial hydroclimate and ecosystems. Yet historical changes in land RH are not well understood due to limited observations, biased reanalyses, and the lack of a framework for interpreting RH changes under multiple influencing factors. Here, we show that the spatiotemporal variability of land RH and its distinct historical trends among observations, reanalyses, and Earth system models are captured by a simple index based on the ratio of precipitation (P) to a modified potential evapotranspiration formulated independently of RH (PET°). The index provides a physical calibration of biased land RH in reanalyses and a quantitative framework for interpreting land RH changes. Over 1973–2024, land RH has decreased substantially, owing to the intrinsic rise in PET° with temperature and little increase in land precipitation. Reanalyses overestimate the observed RH decrease, consistent with exaggerated surface warming and precipitation decline. The index captures this coherent bias and enables a calibration using observed precipitation and temperature. Models simulate a wide range of land RH trends, but nearly all runs underrepresent the historical drying. The index captures the model spread and discrepancy and attributes them to contributions of precipitation and PET°. Weaker land RH decreases in models arise mainly from weaker subtropical precipitation declines, linked to muted intensification of subtropical highs and biased subtropical climatology. The model–observation discrepancy is unlikely explained by internal variability, implying model underestimation of forced RH decrease and a drier land future than current projections.

Earth system models↗
DOE OSTI2020

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

FY-25 Progress on Computational Modeling of the Water Based NSTF

This report summarizes the system level modeling using RELAP5-3D of the Natural Convection Shutdown Heat Removal Test Facility (NSTF) completed in FY25. This year’s work focuses on a new tank configuration where the inlet of the tank was lowered in elevation by 45”. The stability boundaries of the NSTF are thoroughly studied and stability maps are constructed based on the stability and the oscillation patterns of the system. Five distinct operational modes are identified, namely single-phase liquid, uniform double peak oscillations, uniform sinusoidal oscillations, stable two-phase flow, and non-uniform oscillations. Next, the riser inlet throttling case of experimental test Run-104 is simulated with the RELAP5 model where good agreement is obtained between the model and the experimental data. The simulation also highlights the effects of backflow of water from the tank to the upper region of the chimney. Additionally, the decay heat removal test of Run-99 is simulated with the RELAP5 model. Comparison is carried out between this run and a similar run with the mid-tank inlet of Run-74 performed in FY22. With the lower tank inlet, the RELAP5 model is able to predict the experimental data more accurately than the previous mid tank inlet configuration. The discrepancy in model prediction accuracy highlights the non-symmetrical spatial effects in the tank that would otherwise be more easily captured with higher fidelity models. Lastly, two exploratory studies are conducted to investigate the behaviors of the NSTF when 1) heating is provided to the downcomer and 2) a bypass channel is added between the horizontal chimney section to the downcomer.

42 ENGINEERING↗
DOE OSTI2023

An Assessment of Kinetic Models for Ammonia Flame Extinction

The growing interest in using ammonia as a carbon-neutral replacement fuel has prompted numerous recent research efforts towards understanding fundamental combustion characteristics and chemical kinetics of ammonia oxidation. Despite these recent efforts, there remains a large deficiency in extinction limit data for ammonia flames, particularly under elevated pressure and temperature conditions relevant to gas turbines. Even the literature reaction mechanisms developed specifically for ammonia have been untested under gas turbine relevant conditions, and existing kinetic models have shown large variations in laminar flame speed, extinction strain rate, and speciation predictions for NH3–H2 mixtures at ambient conditions. Recent work by Thomas and co-workers, for example, compared different ammonia kinetic models in the literature for predicting their extinction strain rate measurements for non-premixed NH3-H2 counterflow flames, and again have highlighted the discrepancies in model predictions. Here, our work has undertaken to further highlight important reaction pathways for ammonia oxidation that may be responsible for these discrepancies for extinction strain rate predictions, as well as elucidate sensitive reaction steps with large uncertainties that require further attention. This work points to a need for an improved understanding of ammonia combustion chemistry as well as new datasets at elevated pressure and temperature conditions to better constrain rate expressions for key reaction steps.

ammonia, flame, combustion, kinetic, kinetics↗
DOE OSTI2023

Reproducible benchmark for the SNAP 8 experimental reactor at dry conditions

Here, this work provides fully reproducible benchmark models of the Systems for Nuclear Auxiliary Power (SNAP) 8 Experimental Reactor. Validation benchmarks of the criticality configuration experiments under dry conditions with no coolant are presented. In addition, relevant reactivity worth experiments, and measurement of power distributions are validated and presented here. Discrepancies between modeled and experimental results are at most 110 pcm for fuel and poison worths, and critical configurations that do not manipulate control elements are within 50 pcm. Larger discrepancies found in reflector element worth experiments are due to experimental methodology which was noted as being simplified due to limited calculational capabilities at the time, while discrepancy in the control element manipulated critical configuration is due the absence of structural components for modeling simplification. Power distributions closely follow what was observed in experiment with expected peaking in reflecting elements. All models are well documented with cited references describing and justifying assumptions, simplifications, and adjustments to reproduce the models with any nuclear safety codes; model inputs and outputs are stored in the snapReactors GitHub repository.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

Understanding water and energy fluxes in the Amazonia: Lessons from an observation‐model intercomparison

Abstract Tropical forests are an important part of global water and energy cycles, but the mechanisms that drive seasonality of their land‐atmosphere exchanges have proven challenging to capture in models. Here, we (1) report the seasonality of fluxes of latent heat (LE), sensible heat ( H ), and outgoing short and longwave radiation at four diverse tropical forest sites across Amazonia—along the equator from the Caxiuanã and Tapajós National Forests in the eastern Amazon to a forest near Manaus, and from the equatorial zone to the southern forest in Reserva Jaru; (2) investigate how vegetation and climate influence these fluxes; and (3) evaluate land surface model performance by comparing simulations to observations. We found that previously identified failure of models to capture observed dry‐season increases in evapotranspiration (ET) was associated with model overestimations of (1) magnitude and seasonality of Bowen ratios (relative to aseasonal observations in which sensible was only 20%–30% of the latent heat flux) indicating model exaggerated water limitation, (2) canopy emissivity and reflectance (albedo was only 10%–15% of incoming solar radiation, compared to 0.15%–0.22% simulated), and (3) vegetation temperatures (due to underestimation of dry‐season ET and associated cooling). These partially compensating model‐observation discrepancies (e.g., higher temperatures expected from excess Bowen ratios were partially ameliorated by brighter leaves and more interception/evaporation) significantly biased seasonal model estimates of net radiation ( R n ), the key driver of water and energy fluxes (LE ~ 0.6 R n and H ~ 0.15 R n ), though these biases varied among sites and models. A better representation of energy‐related parameters associated with dynamic phenology (e.g., leaf optical properties, canopy interception, and skin temperature) could improve simulations and benchmarking of current vegetation–atmosphere exchange and reduce uncertainty of regional and global biogeochemical models.

Restrepo‐Coupe, Natalia↗
DOE OSTI2023

Thermal transport in warm dense matter revealed by refraction-enhanced x-ray radiography with a deep-neural-network analysis

Abstract Transport properties of high energy density matter affect the evolution of many systems, ranging from the geodynamo in the Earth’s core, to hydrodynamic instability growth in inertial confinement fusion capsules. Large uncertainties of these properties are present in the warm dense matter regime where both plasma models and condensed matter models become invalid. To overcome this limit, we devise an experimental platform based on x-ray differential heating and time-resolved refraction-enhanced radiography coupled to a deep neural network. We retrieve the first measurement of thermal conductivity of CH and Be in the warm dense matter regime and compare our measurement with the most commonly adopted models. The discrepancies observed are related to the estimation of a correction term from electron-electron collisions. The results necessitate improvement of transport models in the warm dense matter regime and could impact the understanding of the implosion performance for inertial confinement fusion.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Visualizing Wind Farm Wakes Using SCADA Data

As wind farms scale to include more and more turbines, questions about turbine wake interactions become increasingly important. Turbine wakes reduce wind speed and downwind turbines suffer decreased performance. The cumulative effect of the wakes throughout a wind farm will therefore decrease the performance of the entire farm. These interactions are dynamic and complicated, and it is difficult to quantify the overall effect of the wakes. This problem has attracted some attention in terms of computational modelling for siting turbines on new farms, but less attention in terms of empirical studies and performance validation of existing farms. In this report, Supervisory Control and Data Acquisition (SCADA) data from an existing wind farm is analyzed in order to explore methods for documenting wake interactions. Visualization techniques are proposed and used to analyze wakes in a 67 turbine farm. The visualizations are based on directional analysis using power measurements, and can be considered to be normalized capacity factors below rated power. Wind speed measurements are not used in the analysis except for data pre-processing. Four wake effects are observed; including wake deficit, channel speed up, and two potentially new effects, single and multiple shear point speed up. In addition, an attempt is made to quantify wake losses using the same SCADA data. Power losses for the specific wind farm investigated are relatively low, estimated to be in the range of 3-5%. Finally, a simple model based on the wind farm geometrical layout is proposed. Key parameters for the model have been estimated by comparing wake profiles at different ranges and making some ad hoc assumptions. A preliminary comparison of six selected profiles shows excellent agreement with the model. Where discrepancies are observed, reasonable explanations can be found in multi-turbine speedup effects and landscape features, which are yet to be modelled.

17 WIND ENERGY↗
DOE OSTI2019

Estimation of inflow uncertainties in laminar hypersonic double-cone experiments

We propose herein a probabilistic framework for assessing the consistency of an experimental dataset, i.e., whether the stated experimental conditions are consistent with the measurements provided. In case the dataset is inconsistent, our framework allows one to hypothesize and test sources of inconsistencies. This is crucial in model validation efforts. The framework relies on Bayesian inference to estimate experimental settings deemed uncertain, from measurements deemed accurate. The quality of the inferred variables is gauged by its ability to reproduce held-out experimental measurements. We test the correctness of the framework on three double-cone experiments conducted in the CUBRC Inc.'s LENS-I shock tunnel, which have also been numerically simulated successfully. Thereafter, we use the framework to investigate two double-cone experiments (executed in the LENS-XX shock tunnel) which have encountered difficulties when used in model validation exercises. We detect an inconsistency with one of the LENS-XX experiments. In addition, we hypothesize two causes for our inability to simulate LEXS-XX experiments accurately and test them using our framework. We find that there is no single cause that explains all the discrepancies between model predictions and experimental data, but different causes explain different discrepancies, to larger or smaller extent. We end by proposing that uncertainty quantification methods be used more widely to understand experiments and characterize facilities, and we cite three different methods to do so, the third of which we present in this paper.

Thermophysics↗
DOE OSTI2023

P2M Simulation Exercise on Past Fuel Melting Irradiation Experiments

This paper presents the results of the Power To Melt and Maneuverability (P2M) Simulation Exercise on past fuel melting irradiation experiments, organized within the Organisation for Economic Co-operation and Development/Nuclear Energy Agency Framework for IrraDiation ExperimentS (FIDES) framework by the Core Group (CEA, EDF, and SCK·CEN) and open to all FIDES members. The exercise consisted in simulating two past power ramps where fuel melting was detected: (1) the xM3 staircase power transient [ramp terminal level (RTL) 70 kW·m -1 , average burnup 27 GWd·tU -1 ], carried out in 2005 in the R2 reactor at Studsvik (Sweden), where the rodlet maintained its integrity, and (2) the HBC4 fast power transient (RTL 66 kW·m -1 , average burnup 48 GWd·tU -1 ), carried out in 1987 in the BR2 reactor at SCK·CEN (Belgium), where the cladding failed during the experiment. The exercise was joined by 13 organizations from 9 countries using 11 different fuel performance codes. In this paper, the main results of the Simulation Exercise are presented and compared to available postirradiation examinations (PIE) or on-line measurements during the power ramps (fuel and clad diameters, rod elongation, pellet-clad gap, and fission gas release). Since the focus of the Simulation Exercise is on fuel melting assessment, determination of the boundary between melted/nonmelted fuel and the consequent definition of a melting radius from PIE are first discussed. During the HBC4 ramp, fuel melting was predicted by most of the codes despite differences in the melting models. Higher discrepancies were observed for the xM3 rod that can be attributed partly to power uncertainty and partly to the limited capability of the models to describe partial melting of the fuel during this ramp. Finally, possible code developments to improve simulation results are presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗