Search

The evolution and release of fission gas impacts the performance of UO2 nuclear fuel. We have created a Bayesian framework to calibrate a novel model for fission gas transport that predicts diffusion rates of uranium and xenon in UO2 under both thermal equilibrium and irradiation conditions. Data sets are taken from historical diffusion, gas release, and thermodynamic experiments. These data sets consist invariably of summary statistics, including a measurement value with an associated uncertainty. Our calibration strategy uses synthetic data sets in order to estimate the parameters in the model, such that the resulting model predictions agree with the reported summary statistics. In doing so, the reported uncertainties are effectively reflected in the inferred uncertain parameters. Furthermore, to keep our approach computationally tractable, we replace the fission gas evolution model by a polynomial surrogate model with a reduced number of parameters, which are identified using global sensitivity analysis. We discuss the efficacy of our calibration strategy, and investigate how the contribution of the different data sets, taken from multiple sources in the literature, can be weighted in the likelihood function constructed as part of our Bayesian calibration setup, in order to account for the different number of data points in each set of data summaries. Our results indicate a good match between the calibrated diffusivity and non-stoichiometry predictions and the given data summaries. We demonstrate a good agreement between the calibrated xenon diffusivity and the established fit from Turnbull et al. (1982), indicating that the dominant uranium vacancy diffusion mechanism in the model is able to capture the trends in the data.

We construct a global sensitivity analysis framework for a coupled multiphysics model used to predict the changes in material properties and surface morphology of helium plasma-facing components in future fusion reactors. The model combines the particle dynamics simulator F-TRIDYN, that predicts the helium implantation profile, with the cluster dynamics simulator Xolotl, that predicts the growth and evolution of subsurface helium gas bubbles. In order to keep the sensitivity analysis tractable, we first construct a sparse, high-dimensional polynomial chaos expansion surrogate model for each output quantity of interest, which allows the efficient extraction of sensitivity information. The sensitivity analysis is performed for two problem settings: one for ITER-like conditions, and one that resembles conditions inside the PISCES-A linear plasma device. We present a systematic comparison of important parameters, for both F-TRIDYN and Xolotl in isolation as well as for the coupled model, and discuss the physical interpretation of these results.

Deep operator networks (DeepONets) are powerful architectures for fast and accurate emulation of complex dynamics. As their remarkable generalization capabilities are primarily enabled by their projection-based attribute, we investigate connections with low-rank techniques derived from the singular value decomposition (SVD). We demonstrate that some of the concepts behind proper orthogonal decomposition (POD)-neural networks can improve DeepONet's design and training phases. These ideas lead us to a methodology extension that we name SVD-DeepONet. Moreover, through multiple SVD analyses, we find that DeepONet inherits from its projection-based attribute strong inefficiencies in representing dynamics characterized by symmetries. Inspired by the work on shifted-POD, we develop flexDeepONet, an architecture enhancement that relies on a pre-transformation network for generating a moving reference frame and isolating the rigid components of the dynamics. In this way, the physics can be represented on a latent space free from rotations, translations, and stretches, and an accurate projection can be performed to a low-dimensional basis. In addition to flexibility and interpretability, the proposed perspectives increase DeepONet's generalization capabilities and computational efficiencies. For instance, we show flexDeepONet can accurately surrogate the dynamics of 19 variables in a combustion chemistry application by relying on 95% less trainable parameters than the ones of the vanilla architecture. We argue that DeepONet and SVD-based methods can reciprocally benefit from each other. In particular, the flexibility of the former in leveraging multiple data sources and multifidelity knowledge in the form of both unstructured data and physics-informed constraints has the potential to greatly extend the applicability of methodologies such as POD and PCA.

The i-V curve describes the current drawn from a flame as a function of the voltage difference applied across the reaction zone. Since combustion diagnostics and flame control strategies based on electric fields depend on the amount of current drawn from flames, there is significant interest in modeling and understanding i-V curves. We implement and apply a detailed model for the simulation of the production and transport of ions and electrons in one dimensional premixed flames. An analytical reduced model is developed based on the detailed one, and analytical expressions are used to gain insight into the characteristics of the i-V curve for various flame configurations. In order for the reduced model to capture the spatial distribution of the electric field accurately, the concept of a dead zone region, where voltage is constant, is introduced, and a suitable closure for the spatial extent of the dead zone is proposed and validated. The results from the reduced modeling framework are found to be in good agreement with those from the detailed simulations. The saturation voltage is found to depend significantly on the flame location relative to the electrodes, and on the sign of the voltage difference applied. Furthermore, at sub-saturation conditions, the current is shown to increase linearly or quadratically with the applied voltage, depending on the flame location. These limiting behaviors exhibited by the reduced model elucidate the features of i-V curves observed experimentally. The reduced model relies on the existence of a thin layer where charges are produced, corresponding to the reaction zone of a flame. Consequently, the analytical model we propose is not limited to the study of premixed flames, and may be applied easily to others configurations, e.g. nonpremixed counterflow flames., Comment: Proc. Combust. Inst., 36, 2016