Combining experimental data and computer simulations, with an application to flyer plate experiments
Authors:
- Jim GattikerJGJim Gattiker
- Dave HigdonDHDave Higdon
- Sallie Keller-McNultySKSallie Keller-McNulty
- Michael McKayMMMichael McKay
- Leslie MooreLMLeslie Moore
- Brian WilliamsBWBrian Williams
Open Access
- 1 December 2006
- journal article
- Published by Institute of Mathematical Statistics in Bayesian Analysis
- Vol. 1 (4) , 765-792
Abstract
A flyer plate experiment involves forcing a plane shock wave through stationary test samples of material and measuring the free surface velocity of the target as a function of time. These experiments are conducted to learn about the behavior of materials subjected to high strain rate environments. Computer simulations of flyer plate experiments are conducted with a two-dimensional hydrodynamic code developed under the Advanced Strategic Computing (ASC) program at Los Alamos National Laboratory. This code incorporates physical models that contain parameters having uncertain values. The objectives of the analyses presented in this paper are to assess the sensitivity of free surface velocity to variations in the uncertain inputs, to constrain the values of these inputs to be consistent with experiment, and to predict free surface velocity based on the constrained inputs. We implement a Bayesian approach that combines detailed physics simulations with experimental data for the desired statistical inference (Kennedy and O'Hagan 2001; Higdon, Kennedy, Cavendish, Cafeo and Ryne 2004). The approach given here allows for: uncertainty regarding model inputs (i.e. calibration); accounting for uncertainty due to limitations on the number of simulations that can be carried out; discrepancy between the simulation code and the actual physical system; and uncertainty in the observation process that yields the actual field data on the true physical system. The resulting analysis accomplishes the objectives within a unified framework.Keywords
This publication has 19 references indexed in Scilit:
- Sensitivity analysis when model outputs are functionsReliability Engineering & System Safety, 2006
- Probabilistic Sensitivity Analysis of Complex Models: A Bayesian ApproachJournal of the Royal Statistical Society Series B: Statistical Methodology, 2004
- An Algorithm for Constructing Orthogonal and Nearly-Orthogonal Arrays With Mixed Levels and Small RunsTechnometrics, 2002
- Bayesian Calibration of Computer ModelsJournal of the Royal Statistical Society Series B: Statistical Methodology, 2001
- Bayesian Forecasting for Complex Systems Using Computer SimulatorsJournal of the American Statistical Association, 2001
- Spatiotemporal Hierarchical Bayesian Modeling Tropical Ocean Surface WindsJournal of the American Statistical Association, 2001
- Controlling Correlations in Latin Hypercube SamplesJournal of the American Statistical Association, 1994
- Orthogonal Array-Based Latin HypercubesJournal of the American Statistical Association, 1993
- Design and Analysis of Computer ExperimentsStatistical Science, 1989
- ACED: Algorithms for the Construction of Experimental DesignsThe American Statistician, 1985