Fragmentation Prediction for a Brittle Ceramic: a Two-Scale Model Approach

Glass and oxide-ceramic radioactive material forms (e.g., PuO₂) are susceptible to brittle fracture during handling and transportation. Fracture of these materials may result in the generation of fragments small enough to become airborne, and thus easily dispersed and inhaled. Entrainment and dispersal of particles with dimensions less than 10 μm is generally anticipated under normal aerodynamic conditions. The assessment of the risk to the public and the environment from these materials due to events resulting in their fracture is based on the determination of the source term for environmental release. This source term is typically quantified by an airborne release fraction (ARF) and a respirable fraction (RF) of material generated by a given event. While the recommended empirical bounding values prove satisfactory to comply with the regulatory process for handling and transportation of these materials, they are based on test data for only a few materials subjected to a limited range of loading conditions. Because the final state of fragmentation (in terms of the size and distribution of fragments) is strongly dependent not only on the properties of the material and the energetic loading conditions that lead to breakup, but also on the many intrinsic length scales associated with fragmentation, it is difficult to precisely predict fragmentation characteristics for other brittle solid wastes forms subjected to a wider range of loading conditions based on the test data alone. The purpose of this manuscript is to present a two-scale modeling approach that may be used to simulate the dynamic fragmentation of brittle solid radioactive material forms subjected to general loading and boundary conditions. The method provides the resulting fragment characteristics (average size and size distribution) across the entire fragment range of interest. At the macro-scale, the approach uses a gradient damage elastic continuum mechanics material model within a finite element framework to simulate dynamic fracture under general loading and boundary conditions. Results from the macroscopic calculations are passed to a one-dimensional (1-D) lower length scale model that incorporates the relevant fine scale effects to determine the micro-scale fragmentation. This 1-D micro-scale fragmentation model simulates the internal crack nucleation and opening processes to predict the average fragment size and the fragment size distribution as a function of the material properties, loading conditions, and microscopic characteristic length scale of the material. Results from the lower length scale model are then combined with the macro-scale model fragmentation information to determine the resulting fragment size distribution spanning both length scales. This two-scale model has been exercised on a representative problem drawn from experimental studies, consisting of an impact of a falling weight onto a brittle ceramic. (authors)