Efficient Probabilistic Analysis Methods for Complex Numerical Models with Large Numbers of Random Variables, 209190 Printer Friendly VersionPrincipal Investigators Inclusive Dates: 04/01/00  Current Background  Physicsbased probabilistic reliability and risk analysis of engineered and natural system is emerging as a less expensive addition to field and laboratorytestbased methods. However, new challenges exist because highly complicated physicsbased models are computationally intensive and have a large number of variables required to conduct accurate systemlevel failure analyses. The TotalSystem Performance Assessment code developed at the Center for Nuclear Waste Regulatory Analyses (a division of the Southwest Research Institute) to evaluate a highlevel radioactive waste (HLW) disposal facility and the codes for reliabilitybased motor vehicle design involving hundreds of input variables are typical examples. These problems cannot be solved efficiently and accurately by the existing methods that have been successfully used for problems with a small number of variables. Approach  The objective of the project is to develop and demonstrate a fast and accurate method to tackle risk and reliability problems involving computationally intensive and numerically complex models with large numbers of parameters. The framework of the method features a hybrid approach that combines the sampling approach that explores the parameter space and the advanced reliability methods that focus the analysis at the tail of the probability distribution. The use of probabilistic sensitivity measures that can effectively identify the most important random variables is being developed. The method is being extended to combine the lowprobability highconsequence scenarios with the most likely scenarios to obtain a single performance measure. An efficient reliability analysis method and an efficient reliabilitybased design optimization method are being developed. The demonstration examples include nuclear waste management, oil and gas exploration and transportation, and automotive applications. Accomplishments  A formal and comprehensive study of the method was pursued by investigating two samplingbased sensitivity measures, (∂p/p)/(∂µ_{i }/σ_{i}) and (∂p/p)/(∂σ_{i }/σ_{i}), where p, µ, and σ are the CDF function, the mean, and the standard deviation, respectively. Analytical solutions were derived for linear response functions with Gaussian variables and numerical solutions for secondorder response surfaces with normal or lognormal variables. The analytical solutions suggest that, in general, the (∂p/p)/(∂µ_{i }/σ_{i}) measure has a better discriminating power in identifying influential variables than the (∂p/p)/(∂σ_{i }/σ_{i}) measure; however, both measures have good discriminating power at the tails of the distribution. The method and its application to an HLW problem were published recently in the Reliability Engineering and System Safety peerreviewed journal (volume 73, p 167176, 2001). Two new sensitivity measures developed since then are particularly relevant to HLW disposal regulatory criteria. These measures are referred to as performancemeanbased sensitivity measures, ∂µ_{Y}/∂µ_{Xi}_{ } and ∂µ_{Y}/∂σ_{Xi}_{ }, where Y is the model output, X_{i} are input parameters, and µ_{Y} is the peak expected dose. Based on ∂µ_{Y}/∂µ_{Xi}_{ }, eleven out of 246 variables are identified as having significant sensitivities at a 95percent acceptance limit. Also, a close relationship between (∂p/p)/(∂σ_{i }/σ_{i} ) and ∂µ_{Y}/∂µ_{Xi}_{ } has been discovered. Based on the calculated ∂µ_{Y}/∂µ_{Xi}_{ }, eleven variables are identified as having significant sensitivities at a 95percent acceptance limit. The top four variables are in the same rank as those identified by ∂µ_{Y}/∂µ_{Xi}_{ }. It appears that fewer than 1,000 samples may be sufficient to identify influential variables, which is very encouraging with respect to the computational burden typically associated with large problems. Application of the sensitivity analysis approach is currently underway for an oil and gas recovery problem.
