A Gaussian process emulator approach for rapid contaminant characterization with an integrated multizone-CFD model

Abstract This paper explores a Gaussian process emulator based approach for rapid Bayesian inference of contaminant source location and characteristics in an indoor environment. In the pre-event detection stage, the proposed approach represents transient contaminant fate and transport as a random function with multivariate Gaussian process prior. Hyper-parameters of the Gaussian process prior are inferred using a set of contaminant fate and transport simulation runs obtained at predefined source locations and characteristics. This paper uses an integrated multizone-CFD model to simulate contaminant fate and transport. Mean of the Gaussian process, conditional on the inferred hyper-parameters, is used as a computationally efficient statistical emulator of the multizone-CFD simulator. In the post event-detection stage, the Bayesian framework is used to infer the source location and characteristics using the contaminant concentration data obtained through a sensor network. The Gaussian process emulator of the contaminant fate and transport is used for Markov Chain Monte Carlo sampling to efficiently explore the posterior distribution of source location and characteristics. Efficacy of the proposed method is demonstrated for a hypothetical contaminant release through multiple sources in a single storey seven room building. The method is found to infer location and characteristics of the multiple sources accurately. The posterior distribution obtained using the proposed method is found to agree closely with the posterior distribution obtained by directly coupling the multizone-CFD simulator with the Markov Chain Monte Carlo sampling.

[1]  Qingyan Chen,et al.  Using CFD Capabilities of CONTAM 3.0 for Simulating Airflow and Contaminant Transport in and around Buildings , 2010 .

[2]  M. J. Box.  A New Method of Constrained Optimization and a Comparison With Other Methods , 1965, Comput. J..

[3]  Henry P. Wynn,et al.  Screening, predicting, and computer experiments , 1992 .

[4]  Dave Higdon,et al.  Combining Field Data and Computer Simulations for Calibration and Prediction , 2005, SIAM J. Sci. Comput..

[5]  J. Koenderink.  Q… , 2014, Les noms officiels des communes de Wallonie, de Bruxelles-Capitale et de la communaute germanophone.

[6]  Liangzhu Wang,et al.  Forecasting simulations of indoor environment using data assimilation via an Ensemble Kalman Filter , 2013 .

[7]  James Axley.  Indoor air quality modeling :: phase II report , 1987 .

[8]  Jin Wen.  Sensor system design for building indoor air protection , 2005, SPIE Optics East.

[9]  Liangzhu Wang.  Coupling of multizone and CFD programs for building airflow and contaminant transport simulations , 2007 .

[10]  A. O'Hagan,et al.  Bayesian calibration of computer models , 2001 .

[11]  W. K. Hastings,et al.  Monte Carlo Sampling Methods Using Markov Chains and Their Applications , 1970 .

[12]  Q Chen,et al.  Theoretical and numerical studies of coupling multizone and CFD models for building air distribution simulations. , 2007, Indoor air.

[13]  Russell W. Wiener,et al.  Short-term dispersion of indoor aerosols : can it be assumed the room is well mixed? , 2006 .

[14]  Clifford Federspiel,et al.  Estimating the inputs of gas transport processes in buildings , 1997, IEEE Trans. Control. Syst. Technol..

[15]  John J. Borkowski,et al.  SIMPLE LATIN SQUARE SAMPLING + 1 : A SPATIAL DESIGN USING QUADRATS , 1996 .

[16]  Richard J. Beckman,et al.  A Comparison of Three Methods for Selecting Values of Input Variables in the Analysis of Output From a Computer Code , 2000, Technometrics.

[17]  Zhiqiang John Zhai,et al.  Prompt tracking of indoor airborne contaminant source location with probability-based inverse multi-zone modeling , 2009 .

[18]  Helmut E. Feustel,et al.  COMIS—an international multizone air-flow and contaminant transport model , 1998 .

[19]  Marios M. Polycarpou,et al.  Security-oriented sensor placement in intelligent buildings , 2013 .

[20]  P V Nielsen,et al.  Computational fluid dynamics and room air movement. , 2004, Indoor air.

[21]  N. Metropolis,et al.  Equation of State Calculations by Fast Computing Machines , 1953, Resonance.

[22]  T. J. Mitchell,et al.  Bayesian Prediction of Deterministic Functions, with Applications to the Design and Analysis of Computer Experiments , 1991 .

[23]  Ashok J. Gadgil,et al.  Towards improved characterization of high-risk releases using heterogeneous indoor sensor systems , 2011 .

[24]  Alex K. Jones,et al.  Indoor environmental quality in a dynamic life cycle assessment framework for whole buildings: Focus on human health chemical impacts , 2013 .

[25]  Y. Lisa Chen,et al.  The selection of the most appropriate airflow model for designing indoor air sensor systems , 2012 .

[26]  A. Bagtzoglou,et al.  State of the Art Report on Mathematical Methods for Groundwater Pollution Source Identification , 2001 .

[27]  J. Oakley.  Eliciting Gaussian process priors for complex computer codes , 2002 .

[28]  A. OHagan,et al.  Bayesian analysis of computer code outputs: A tutorial , 2006, Reliab. Eng. Syst. Saf..

[29]  Michael D. Sohn,et al.  Influence of indoor transport and mixing time scales on the performance of sensor systems for characterizing contaminant releases , 2007 .

[30]  A. O'Hagan,et al.  Bayesian emulation of complex multi-output and dynamic computer models , 2010 .

[31]  Etienne Wurtz,et al.  Two- and three-dimensional natural and mixed convection simulation using modular zonal models in buildings , 1999 .

[32]  Michael Goldstein,et al.  Probabilistic Formulations for Transferring Inferences from Mathematical Models to Physical Systems , 2005, SIAM J. Sci. Comput..

[33]  Bithin Datta,et al.  Optimal Monitoring Network and Ground-Water–Pollution Source Identification , 1997 .

[34]  Michael D. Sohn,et al.  A stiff, variable time step transport solver for CONTAM , 2013 .

[35]  Rhodri Hayward,et al.  Screening , 2008, The Lancet.

[36]  Zhiqiang Zhai,et al.  Application of CFD to Predict and Control Chemical and Biological Agent Dispersion in Buildings , 2003 .

[37]  Thomas J. Santner,et al.  Design and analysis of computer experiments , 1998 .

[38]  Priya Sreedharan,et al.  Bayesian based design of real-time sensor systems for high-risk indoor contaminants , 2007 .

[39]  Leon R. Glicksman,et al.  Application of integrating multi-zone model with CFD simulation to natural ventilation prediction , 2005 .

[40]  Michael D. Sohn,et al.  Evaluating sensor characteristics for real-time monitoring of high-risk indoor contaminant releases , 2006 .

[41]  Geoffrey E. Hinton,et al.  Evaluation of Gaussian processes and other methods for non-linear regression , 1997 .

[42]  Marc C. Kennedy,et al.  Case studies in Gaussian process modelling of computer codes , 2006, Reliab. Eng. Syst. Saf..

[43]  William W. Nazaroff,et al.  Mixing of a Point Source Pollutant by Natural Convection Flow within a Room , 1994 .

[44]  Thomas J. Santner,et al.  The Design and Analysis of Computer Experiments , 2003, Springer Series in Statistics.

[45]  R. Adler,et al.  Random Fields and Geometry , 2007 .

[46]  X Liu,et al.  Inverse modeling methods for indoor airborne pollutant tracking: literature review and fundamentals. , 2007, Indoor air.

[47]  Y. Lisa Chen,et al.  Comparison of sensor systems designed using multizone, zonal, and CFD data for protection of indoor environments , 2010 .

[48]  M. Sohn,et al.  Rapidly Locating and Characterizing Pollutant Releases in Buildings , 2000, Journal of the Air & Waste Management Association.