Fracture Probability Integral Applied to Reactor Vessel Life Estimate¹

The conventional method of fracture probability calculations such as that adopted by the NRC-sponsored PRAISE CODE and the FAVOR CODE developed in this laboratory are both based on Monte Carlo simulation. Heavy computations are required. A new method of fracture probability calculation is developed by direct probability integration. The preliminary version of the development was published in an earlier paper. More detailed development of the method is presented here. The present approach offers simple and expedient method to obtain numerical values of fracture probability. This method can be applied to problems as general as the method of Monte Carlo simulation. This approach also provides a clear physical picture on the meaning of the probability of fracture. Parametric studies are made in this paper to show the variation of the numerical values of the probabilities of fracture as a result of the change of the standard deviation of either fracture toughness or the radiation-induced temperature shift. Also, it is shown numerically that a limiting probability can be obtained if the standard deviation of the fracture toughness approaches zero that implies a deterministic fracture toughness. It confirms the theoretical proof shown in Eq. (¹¹). The limiting probability is the simplistic probability of crack count used by this author where both toughness and temperature shift are assumed to be deterministic values. The general probability of fracture developed here is simply a generalization of the crack count, except the crack count is selected with the appropriate fracture toughness in the toughness distribution. The toughness for the problem considered here is then multiplied by the appropriate temperature shift in the distribution function of the temperature shift. Although the present development is based on linear fracture mechanics assumption and applied to the radiated reactor vessel steel, there is no difficulty in viewing the present development as a general formulation that is capable of handling as many random variables as required by the fracture model. The multiplicity of the integration corresponds to the number of random variables. The probability integral is applied in this paper to calculate the probability of fracture for the high flux isotope reactor (HFIR) vessel that has been weakened due to the radiation embrittlement. The random variables used here are the crack length, the fracture toughness, and the radiation-induced temperature shift that is needed in the parametric representation of the radiated vessel steel.