Calibration of the CO + OH CO 2 + H rate constant with direct and indirect information
|
|
- Hortense Cobb
- 5 years ago
- Views:
Transcription
1 Paper # 070RK-0246 Topic: Reaction Kinetics 8 th US National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, Calibration of the CO + OH CO 2 + H rate constant with direct and indirect information S. Plessis PECOS, ICES, University of Austin, TX While the literature contains a vast array of kinetics parameterizations of chemical reactions, identical (or nearly identical) physical models have reported parameter values that can vary over several orders of magnitude. In this work, the investigation of the kinetics description of the reaction CO + OH CO 2 + H is presented. Multiple models have been calibrated and compared using bayesian techniques with data of two different types: measurements of (k, T ) and laminar flamespeeds. First, measurements of (k, T ) in the range of K were used. These data are taken from several different sources. Three different models were calibrated and their results were compared. All three model formulations performed equally well when evaluated in the bayesian context. These direct measurements of (k, T ) were not as informative as one might expect resulting in highly-correlated posterior distributions that vary over several orders of magnitude. These results can be used to explain the large discrepancies found in the literature. In addition to using the above data, a second set of calibrations were performed using indirect measurements: laminar flamespeed in a CO/H 2 mixture. By starting with the GRI-Mech 3.0 [1] mechanism and using very large prior distributions on the model parameters (signifying a lack of knowledge), the flamespeed data were used to inform the model parameters via a statistical inverse problem. Calibrating against this information did not lead to agreement with direct measurements, therefore yielding unrealistic descriptions of the rate constants. Thus, it is reasonable to conclude that the previous calibration alone is not sufficient for reliable flamespeed calculations. In conclusion, this work proposes a methodology to parameterize rate constants in low-pressure conditions that appears to be consistent with all available experimental data of (k, T ). However, rate constant measurements alone do not sufficiently constrain kinetics parameters to reliably predict laminar flamespeeds and additional information is required in order to parameterize a model that is both consistent with (k, T ) experimental data that also accurately predicts laminar flamespeeds. 1 Introduction The bimolecular CO + OH CO 2 + H reaction plays an important role in hydrogen combustion systems [2] and has two competing reactive channels [3, Fig. 3], which makes it difficult to represent by a general kinetics model. This reaction is known to depend heavily on pressure [4, 5]. In this work, we will only consider atmospheric pressure for combustion in order to be able to neglect this dependence [6]. Mechanisms and kinetics models used in the literature [6 9] usually derive from already compiled databases [10] and then update all or a subset of the kinetics parameters [9] by integrating indirect Corresponding author splessis@ices.utexas.edu 1
2 information. The GRI-Mech 3.0 mechanism [1] is built by calibrating all kinetics parameters of an extended chemical mechanism against numerous experiments (direct and indirect information) in order to obtain the most general model for gas-phase chemistry, while Davis et al. updated only the pre-exponential parameters (see Eq. 1) by using a large set of experimental measurements [9, Fig. 1 and Tab. 1]. One of the downsides of this approach is that it provides model-dependent kinetics characterization, which translates for the CO + OH CO 2 + H chemical reaction to orders of magnitude discrepancies on kinetics parameters (Tab. 1). A reaction is modeled by an elementary process or a duplicate process when the pressure conditions allow one to neglect pressure dependent models as in this work. These processes are represented by one or the sum of several Kooij equations, respectively (Eq. 1). n 1 ( ) βi ( T k(t ) = A i exp E ) ai T ref RT i=1 with k the rate constant, T the temperature, T ref a reference temperature (in this work T ref = 300 K), A the pre-exponential factor, β the power parameter, E a the activation energy and R the ideal gas constant. In this work, we use different kinetics models found in the literature, i.e. n = 1, 2 and 3, respectively called M 1, M 2 and M 3. The M 1 model is used for elementary processes while the M 2 and M 3 models add degrees of freedom to account for chemical processes not well approximated by an elementary process. The most common kinetics model for this reaction is the simple Kooij equation M 1 [1, 2, 7, 8, 10, 11], while Davis et al. [9] optimized the rate constant using a sum of two Kooij equations (M 2 ), and Sun et al. [6] used the M 3 form (Tab. 1) following the low pressure expression of Troe et al. [5]. See Tab. 1 for detailed values. In this work, we used bayesian calibration methods to inform the kinetics parameters for these rate constants using direct measurements values of rate constant at low pressures with respect to temperature and indirect information flamespeed measurements, thus obtaining a kinetics description relevant to a combustion model. To inform the model with indirect information, we used a laminar flamespeed as the Quantity of Interest (QoI) on a CO/H 2 combustion model, starting from the GRI-Mech 3.0 [1] reactive scheme. We used the CANTERA library 1 to perform the calculations. Combustion conditions are represented by the parameter Φ. This parameter is used to represent the type of combustion: Φ = 1 is stoichiometric combustion, a Φ < 1 is a lean combustion and a Φ > 1 is a rich combustion. We have reproduced the experimental conditions of McLean et al. [7] on the 5% H 2-95% CO mixture, Φ = and Φ = 3.895, to compare to their results and to the larger calibration of Davis et al. [9], who used a mechanism of 30 chemical reactions and optimized this mechanism over 35 different experimental measurements, including the two considered in this paper. 1 (1) 2
3 Source A β E a T range T ref (cm 3 mol 1 s 1 ) GRI-Mech 3.0 [1] (cal mol 1 ) > Baulch et al. [10] ( Ea in K) R Mueller et al. [8] (cal mol 1 ) Li et al. [2] ( ) 1.89( ) 583( ) ( Ea R (K) in K) Davis et al. [9] (cal mol 1 ) > Sun et al. [6] from Troe et al. [5], (cal mol 1 ) low-pressure limit Table 1: Values of OH + CO H + CO 2 rate constant parameters with some parameterizations available in the literature. The uncertainty reported from Li et al. [2] are the 95% confidence intervals. : in cm 3 molecule 1 s 1. 2 Bayesian calibration 2.1 Parameter calibration In a bayesian framework, parameters are represented by probability density functions (pdf) to account for uncertainties. Bayesian inference is based on Bayes theorem [12 14] (Eq. 2) which formalizes the concept of updating knowledge by integrating new information. p(x y) = p(x)p(y x) p(y) (2) A parameter x, which prior knowledge is represented by the pdf p(x) is updated by the information contained in the hypothesis y. The pdf p(x y) denotes the updated knowledge of x knowing y. For a vector θ of the kinetics parameters of a chemical model, Bayes theorem is expressed as in Eq. 3. π(θ; D) p post (θ D) = p prior (θ) pprior (θ)π(θ; D)dθ (3) with D the vector of information one wishes to use to inform the parameters; p prior the prior pdf, which denotes the state of knowledge of the parameters independently of the data; p post the posterior pdf, which is the state of knowledge after having integrated the information contained in the data D. π is the likelihood function, measure of the agreement between the data D and the model prevision for specific values of the parameters. 3
4 2.2 Model calibration This method is generalizable to models, and thus allow model comparison. Let s consider a set M = {M 1, M 2,..., M K } of K candidate models. Eq. 3 rewritten for a model class M i gives: π(θ i, M i ; D) p post (θ i D, M i ) = p prior (θ i M i ) (4) pprior (θ i M i )π(θ i, M i ; D)dθ i with θ i the parameters of the model class M i. To compare different competitive model classes, one needs to compare the posterior distributions. Using Bayes theorem (Eq. 3), we obtain as posterior pdf for model M i π evid (M i ; D) p post (M i D, M) = p prior (M i M) K k=1 p prior(m k M)π evid (M k ; D) The prior distributions are the probabilities of model classes M i, i [1,K] before integrating the data D. Usually, little information is available, therefore equiprobability is appropriate, i.e. uniform probability (p prior (M i M) = 1/K) to all model classes M i. π evid is the evidence function, analog to the likelihood function in (3). The evidence function measures the consistency of the concerned model and the data considering the entire parameter space. It is given by the following integral over parameter space. π evid (M i ; D) = p prior (θ i M i )π(θ i, M i ; D)dθ i (6) which corresponds to the denominator of (4). It can be rewritten as: log(π evid (M i ; D)) = log(π like (θ, M i ; D))p post (θ D, M i )dθ ( ) ppost (θ D, M i ) log p post (θ D, M i )dθ p prior (θ M i ) The first term is the posterior expectation of the log-likelihood which measures how well the model is able to fit the data, averaged over the posterior pdf for the parameters. The second term is the relative information entropy (or Kullback-Leibler divergence [15]) between the posterior and the prior pdfs. It measures the information about the parameters that is gained from the data. For two model classes that fit the data equally well, as measured by the first term in (7), the model that requires the least tuning, as measured by information gain, is preferred. For more details and discussion, see Jaynes [14] (Chapter 20) and Muto and Beck [16] (Section 4). (5) (7) 2.3 Monte Carlo method Eq. 3 and 4 are usually not solvable analytically, and their estimations are computationally very demanding. Computing the posterior pdfs or the evidence for instance requires solving highdimensional integrals. To evaluate such quantities, we use a Monte Carlo approach [17] in which the pdfs are represented by samples, which enables to approximate the integrals using those samples. In this work, the stochastic method referred as the adaptive multi-level algorithm [18] is used, implemented in the QUESO library [19]. 4
5 2.4 Calibration cycles The kinetics parameters are characterized through two different calibration cycles which enable to integrate different type of information. Calibrating over measurements of a reaction rate will inform directly the kinetics parameters and therefore will provide the most accurate representation of the rate constant considering available knowledge. This calibration will provide a model-agnostic description: an unbiased database, that is, containing only the kinetics measurements biases and uncertainties. Any model calculation using this description will contain all the model inadequacies. This cycle will be called the C I cycle thereafter. A calibration process using a combustion model and measured values of a chosen QoI will inform the kinetics parameters with both the measurements information and the model characteristics. This calibration contains a large amount of information as every chemical reaction has an impact on every other ones. Thus the chosen mechanism (the considered chemical reactions), the kinetics model chosen for every reaction, and, in this case, the transport model have all an impact on the posterior pdfs of the informed quantities. The effect of such calibration is a tuning of the informed parameters to the model and QoI, which will adapt the said-informed parameters to the model and measurements biases. Therefore this calibration will produce the most reliable model in terms of accuracy to the QoI. This cycle will be called the C II cycle thereafter. These two calibration processes are pictured in Fig. 1. Database on kinetics parameters of reaction: (k, T ) surface Reaction alone Kinetics models: M 1 : one Kooij equation, M 2 : two Kooij equations, M 3 : three Kooij equations. Full mechanism Bayes Calibrated kinetics parameters Combustion model Combustion data Figure 1: To obtain information on a reaction, we can infer from measurements on the reaction rate (black path); to obtain the most accurate model, we can infer from measurements on a model s output (red path). 5
6 2.5 Likelihood For the calibration with direct data, they are available in the form of points in the (k, T ) space, k being the rate constant (unit varies) and T the temperature (in K). The uncertainties are considered in the k and T directions. A rate constant is a positive parameter over which measurements disagreements can be large (> 30%). This is well represented by a lognormal distribution [20], leading to use a lognormal expression in the k direction (Eq. 8). The uncertainty on the temperatures is smaller, we chose a normal expression for the T dimension (Eq. 9). For a data point (k data, T data ) with uncertainty ( k data, T data ). The likelihood for a modeled data point, (k M, T M ), will therefore be expressed in the k and T directions. For the k direction, we define the lognormal deviation as F = 1 + k data /k data. π km k data,f = LogN (k M ; k data, F ) ( 1 = 2π kdata ln(f ) exp 1 ( ) ) 2 ln(kdata ) ln(k M ) 2 ln(f ) (8) The likelihood in the T dimension is of the form: π TM T data, T data = Norm(T M ; T data, T data ) ( 1 = exp 1 ( ) ) 2 TM T data 2π Tdata 2 T data (9) For each (k M, T M ) point, we thus define the likelihood as: π (km,t M ) data = 1 N data data i π km k i,f i π TM T i, Ti (10) When the data uncertainty is not provided, representative values are used: 10% uncertainty over k ( k data = 0.1 k data ), and 1% uncertainty over the temperature ( T data = 0.01 T data ). The rate constant is calculated along a interval of temperature. To estimate a likelihood, it is necessary to define a temperature grid on which to calculate the likelihood at each grid point along a k(t ) curve. The arithmetic average will provide the likelihood value at each sampled point. In the indirect information case, we have a measured flamespeed associated to its uncertainty. In such case, a gaussian likelihood is adapted (Eq. 11): π fm f meas, f meas = Norm(f M ; f meas, f meas ) ( 1 = exp 1 ( ) ) 2 fm f meas 2π fmeas 2 f meas (11) with f meas the measured flamespeed, f meas its associated uncertainty and f M the modeled flamespeed. 6
7 Parameters C I C II A U(0, 20) in log 10 (cm 3 mol 1 s 1 ) N (11, 3) in ln(m 3 kmol 1 s 1 ) β U(0, 2) N (1.1, 0.45) E a U( 50, 50) in kj mol 1 N ( 560, 410) in cal mol 1 Table 2: A priori pdfs for the calibration cycles. In the C II cycle, we add the constrains k 300 K and k 2500 K in cm 3 mol 1 s Priors The default priors were chosen as large as possible, in order to simulate at best a complete lack of information. In the C I calibration cycle, the prior considered knowledge was: the pre-exponential factor is positive and lower than cm 3 mol 1 s 1, the temperature dependence is positive but never above the quadratic, and the activation energy is either positive or negative, below an absolute value of 50 kj mol 1. Following the maximum entropy method [21, 22], we thus define a uniform pdf on the log 10 values of the pre-exponential factor, the power parameters and the activation energy (Tab. 2). It was not possible to use as large priors as previously defined for the C II cycle, as this description defines some physically unreasonable rate constants, and prevent the flamespeed calculations to converge in such cases. Thus the priors were changed considering the C I cycle posterior pdfs, and filters were added to ensure convergence, see Tab Results and comparisons 3.1 Calculations conditions The rate constant measurements are taken from seven different sources (Fig. 2), in atmospheric pressure conditions, in agreement to the flamespeed measurements used in cycle C II [7, 9]. The grid used is defined by T start = 100 K, T end = 2500 K and T step = 150 K. The combustion model used is based on the Gri-Mech 3.0 [1] mechanism, in which only the reaction CO + OH H + CO 2 was parametrized by probabilistic values. The transport model used in these calculations is a mixture-averaged model. The laminar flamespeed calculations are performed with the CANTERA chemistry program [23]. The flamespeed measurements and conditions of calculations are taken from McLean et al. [7] to enable comparison with Davis et al. [9]. Results (Tab. 3, 4 and Fig. 3) are interpreted in terms of posterior rate constants (Fig. 4) and flamespeeds (Fig. 5). The different models are compared in terms of performances within the C I cycle (log(ev) column of table 3a), and only the best one in terms of bayesian evidence ( 2.2), will be used in the C II cycle, to limit computational time. 7
8 10 12 k (cm 3 mol 1 s 1 ) Wooldridge et al., 1994 Wooldridge et al., 1996 Westenberg and de Haas, 1973 Jonah et al., 1984 Golden et al., 1998 Frost et al., 1991 Frost et al., ,000 1,500 2,000 T (K) Figure 2: Experimental data for the calibration of CO + OH CO 2 + H. Units have been harmonized to cm 3 mol 1 s Kinetics parameters posterior pdfs The C I inversion cycle provides highly-correlated posterior pdfs that spreads over several orders of magnitude, especially for the M 2 and M 3 models. The posterior pdfs of A have a lognormal deviation between 3 and 7.5, E a s sign is not even characterized. This can explain the widely different parameterizations found in the literature (Fig. 3), as these results agree with all considered comparative results. Yet these results show that direct measurements alone are not sufficient to characterize the kinetics parameters of CO+OH CO 2 +H. As all tested model classes perform equally well in term of bayesian evidence (Tab. 3a), the chosen model for the C II cycle was the simplest, M 1. Posterior pdfs obtained with the C II inversion cycle are more constrained: β is set to 0, which, in agreement to the high-correlation observed in the previous calibration, constrain A to a relatively small interval. E a is still not characterized. 3.3 Resulting rate constants While the C I cycle ensures realistic rate constants description, the posterior pdfs of the kinetics parameters define rate constants that are in blatant disagreement with the direct measurements (Fig. 4, red and green curves). It appears that the information contained in the laminar flamespeed is not accurate with respect to the rate constant. The information integrated in this inversion cycle contain the model inadequacies, which could partly explain the disagreements between calibrated and measured values. 8
9 Φ Model A(F ) β (σ) E a (σ) cm 3 /mol/s cal/mol log(ev) M (4.7) 1.1 (0.5) -525 (400) M 11 (7.5) 1.0 (0.5) 880 (3150) (4.5) 1.2 (0.5) 800 (3750) M (5.3) 1.0 (0.5) 1550 (3650) (3.1) 1.1 (0.5) 3700 (4250) (4.1) 1.0 (0.6) 3000 (4150) (a) C I calibration cycle results A (F ) β (σ) E a (σ) cm 3 /mol/s cal/mol Flamespeed (cm s 1 ) (1.3) 0 (0) 600 (530) 40 ± (1.3) 0 (0) 600 (550) 50 ± 10 (b) C II calibration cycle results Table 3: Results summarized for the bayesian inversion cycles. The uncertainties of inverse parameters are given as standard deviation (σ) or relative multiplicative factor (F = 1 + σ A /A). Above: direct information inverse on the temperature grid T start = 100 K, T end = 2500 K and T step = 150 K. Below: indirect information inverse on the GRI-Mech 3.0 mechanism [1], flamespeed calculations performed by the CANTERA software [23] on a 5% H 2-95% CO mixture. A β E a A β E a 1 (a) C I calibration cycle A β E a A β - - E a 1 (b) C II calibration cycle with Φ = A β E a A β - - E a 1 (c) C II calibration cycle with Φ = Table 4: Rank correlation of the results for the bayesian inversion for a temperature grid of T start = 100 K, T end = 2500 K and T step = 150 K. 9
10 A (cm 3 mol 1 s 1 ) 6,000 4,000 2,000 0 E a (J mol 1 ) β GRI 3.0-Mech Baulch et al, 1992 Mueller et al., 1999 Li et al., 2007 This study, dir. cal. This study, ind. cal. Φ = This study, ind. cal. Φ = Figure 3: Visual comparison of M 1 parametrizations used in the literature and this study results for T ref = 300 K. Li et al., 2007: 95% confidence intervals; this study: 90% confidence intervals. 3.4 Corresponding flamespeed Inversely, the C II cycle provides the most accurate flamespeeds, which is expected, as the flamespeed and model information are integrated in the kinetics model. The combustion model, using the posterior pdfs of cycle C I provides overestimated flamespeed, but not inconsistent, which can be ascribed to model inadequacies, notably the transport model. 3.5 Comparison with Davis et al. Davis et al. [9] used a chemical mechanism close to the Gri-Mech 3.0 [1] (14 species and 30 chemical reactions), a M 2 kinetics model on the CO+OH CO 2 +H reaction and a multi-component transport model. As seen in 3.2, the difference between the kinetics models is negligible: they both provides very similar rate constants (Fig. 4). The mechanisms are also quite similar, as the GRI 3.0 [1] was used as a base for their mechanism. The main expected source of discrepancies is therefore the transport model. The calibration was performed only on A parameters, thus acting like a multiplicative coefficient on the rate constant curves, which translates as a calibration on constrained kinetics models. These constrains enable to ensure physical consistence of the kinetics models (Fig. 4b) without a lost of the flamespeed accuracy (Fig. 5), but necessitate strong prior information on the activation energy and the power parameter. The combustion model used in this paper, while overestimating the flamespeed using the C I cycle results, do not provide inconsistent results. Measurements from McLean and al. [7] and results of Davis et al. [9] are within the flamespeed prediction uncertainty. 10
11 k (cm 3 mol 1 s 1 ) ,000 1,500 2,000 2,500 3,000 T (K) (a) Inversions results k (cm 3 mol 1 s 1 ) GRI-Mech 3.0 Baulch & al., 1992 Mueller & al., 1999 Li & al., 2007 Davis & al., 2005 Sun & al., ,000 1,500 2,000 2,500 3,000 T (K) (b) Models in the literature Figure 4: Above: generated rate constants. Blue: direct calibration; red: indirect calibration with Φ = 1.148; green: indirect calibration with Φ = Below: available models in the literature. 11
12 Davis et al., 2005 Davis et al., Flamespeed (m s 1 ) (a) Φ = Flamespeed (m s 1 ) (b) Φ = Figure 5: Histograms of flamespeed, N sample = 1000, the black solid curve is the target measured by McLean et al. [7]. Blue: direct calibration; red: indirect calibration. Davis et al. [9] is given for comparison. 4 Conclusion This paper proposes a bayesian methodology to calibrate, through a statistical inverse problem, the kinetics parameters for the reaction CO + OH CO 2 + H using two types of information: direct measurements of type (k, T ) and indirect measurements of laminar flamespeed. The first inversion did not constrain the parameters while the second provides unrealistic rate constants. It therefore seems not possible to obtain, in one characterization, realistic kinetics parameters and a reliable flamespeed description for a combustion model. Thus an accurate and reliable description would require more information than kinetics or flamespeed measurements alone. Combining the two calibration cycles into a double-calibrated kinetics model by (i) constraining the kinetics parameters to a physical domain by a C I calibration cycle, and (ii) improving the model s reliability by applying a C II calibration cycle on obtained pdfs by step (i) could constrain the kinetics parameters while ensuring physical consistency and flamespeed calculation accuracy. This is pictured in Fig. 1. Another improvement would be to integrate more information, either more specifically on the kinetics parameters, as done in Davis et al. [9], or on the rate constants by integrating pressure dependence, using measurements at higher temperature, etc. Other interesting possibilities reside in using quantum calculations to provide additional calibration data, which are not constrained by experimental limitations. 12
13 Database on kinetics parameters of reaction: (k, T ) surface Reaction alone Kinetics models: M 1 : one Kooij equation, M 2 : two Kooij equations, M 3 : three Kooij equations. Bayes Full mechanism Calibrated kinetics parameters Combustion data Combustion model Figure 6: Proposed scheme for estimating kinetics parameters. A first calibration cycle (black path) against direct information ensures a physical description of the rate constants. The second calibration cycle (red path) would then integrate the model information to ensure the best flamespeed predictive capabilities of the model. 13
14 References [1] Gregory P. Smith, David M. Golden, Michael Frenklach, Nigel W. Moriarty, Boris Eiteneer, Mikhail Goldenberg, C. Thomas Bowman, Ronald K. Hanson, Soonho Song, William C. Gardiner Jr., Vitali V. Lissianski, and Zhiwei Qin. GRI-Mech 3.0, [2] Juan Li, Zhenwei Zhao, Andrei Kazakov, Marcos Chaos, Frederick L. Dryer, and James J. Scire. Int. J. Chem. Kinet., 39 (2007) [3] R. S. Zhu, E. G. W. Diau, M. C. Lin, and A. M. Mebel. The Journal of Physical Chemistry A, 105 (2001) [4] D. Fulle, H. F. Hamann, H. Hippler, and J. Troe. The Journal of Chemical Physics, 105 (1996) [5] Jürgen Troe. Symposium (International) on Combustion, 27 (1998) Twenty-Seventh Sysposium (International) on Combustion Volume One. [6] Hongyan Sun, S.I. Yang, G. Jomaas, and C.K. Law. Proceedings of the Combustion Institute, 31 (2007) [7] Ian C. McLean, David B. Smith, and Simon C. Taylor. Symposium (International) on Combustion, 25 (1994) Twenty-Fifth Symposium (International) on Combustion. [8] M. A. Mueller, R. A. Yetter, and F. L. Dryer. International Journal of Chemical Kinetics, 31 (1999) [9] Scott G. Davis, Ameya V. Joshi, Hai Wang, and Fokion Egolfopoulos. Proceedings of the Combustion Institute, 30 (2005) [10] D. L. Baulch, C. J. Cobos, R. A. Cox, C. Esser, P. Frank, Th. Just, J. A. Kerr, M. J. Pilling, J. Troe, R. W. Walker,, and J. Warnatz. Journal of Physical and Chemical Reference Data, 21 (1992) [11] NIST. Nist kinetics database. [12] Thomas Bayes and Richard Price. Philosophical Transactions of the Royal Society of London, 55 (1763) [13] R. T. Cox. The Algebra of Probable Inference. Johns Hopkins University Press, [14] E. T. Jaynes. Probability Theory: The Logic of Science. Cambridge University Press, [15] S. Kullback and R. A. Leibler. Annals of Mathematical Statistics, 22 (1951) [16] M. Muto and J. L. Beck. J. Vib. Control, 14 (2008) [17] JCGM. Evaluation of measurement data - Supplement 1 to the GUM: propagation of distributions using a Monte Carlo method. Technical Report 101, BIPM, IEC, IFCC, ISO, IUPAP and OIML, [18] S. H. Cheung and J. L. Beck. New bayesian updating methodology for model validation and robust predictions based on data from hierarchical subsystem tests. Technical report, [19] Ernesto Prudencio and Karl W. Schulz. The Parallel C++ statistical Library QUESO: Quantification of Uncertainty for Estimation, Simulation and Optimization. In Workshop on Algorithms and Programming Tools for Next-Generation High-Performance Scientific Software HPSS, Bordeaux, France, [20] D.L. Smith. A demonstration of the lognormal distribution. Report ANL/NDM-156, Argonne National Laboratory, [21] Hang K. Ryu. Journal of Econometrics, 56 (1993) [22] Lihua Feng and Weihu Hong. Applied Mathematical Modelling, 33 (2009) [23] D. Goodwin. Cantera: An object-oriented software toolkit for chemical kinetics, thermodynamics, and transport processes. Caltech, PASADENA, Available: 14
Exploring The Fundamentals In Catalytic Partial Oxidation Of Methane: The Interaction Between Diffusion And Reaction In A Packed Bed Reactor
Exploring The Fundamentals In Catalytic Partial Oxidation Of Methane: The Interaction Between Diffusion And Reaction In A Packed Bed Reactor Songjun Liu; Ana Obradović; Joris W. Thybaut; Guy B. Marin Laboratory
More informationA REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL ENTROPY INEQUALITY
THE PUBLISHING HOUSE PROCEEDINGS OF THE ROMANIAN ACADEMY, Series A, OF THE ROMANIAN ACADEMY Special Issue/2018, pp. 285 290 A REDUCED-ORDER METHANE-AIR COMBUSTION MECHANISM THAT SATISFIES THE DIFFERENTIAL
More informationHierarchical approach
Chemical mechanisms Examine (i) ways in which mechanisms are constructed, (ii)their dependence on rate and thermodynamic data and (iii) their evaluation using experimental targets Copyright 2011 by Michael
More informationKINETIC MODELING OF OXIDATION METHANE CONVERSION IN REGIME OF FILTRATION COMBUSTION WITH SUPERADIABATIC HEATING
KINETIC MODELING OF OXIDATION METHANE CONVERSION IN REGIME OF FILTRATION COMBUSTION WITH SUPERADIABATIC HEATING Anna A. Karnaukh, Avigeya N. Ivanova, Svetlana S. Kostenko, George B. Manelis, and Eugene
More informationHigh-Pressure Kinetic Mechanisms for Hydrogen and Hydrogen Syngas
High-Pressure Kinetic Mechanisms for Hydrogen and Hydrogen Syngas 1 st International Workshop on Flame Chemistry Warsaw, Poland July 28, 212 Michael P. Burke Chemical Sciences and Engineering Division,
More informationCombustion Chemistry
Combustion Chemistry Hai Wang Stanford University 2015 Princeton-CEFRC Summer School On Combustion Course Length: 3 hrs June 22 26, 2015 Copyright 2015 by Hai Wang This material is not to be sold, reproduced
More informationUQ in Reacting Flows
UQ in Reacting Flows Planetary Entry Simulations High-Temperature Reactive Flow During descent in the atmosphere vehicles experience extreme heating loads The design of the thermal protection system (TPS)
More informationCOMBUSTION CHEMISTRY OF PROPANE: A CASE STUDY OF DETAILED REACTION MECHANISM OPTIMIZATION
Proceedings of the Combustion Institute, Volume 28, 2000/pp. 1663 1669 COMBUSTION CHEMISTRY OF PROPANE: A CASE STUDY OF DETAILED REACTION MECHANISM OPTIMIZATION ZHIWEI QIN, 1 VITALI V. LISSIANSKI, 1 HUIXING
More informationVALIDATION OF THE ASVDADD CONSTRAINT SELECTION ALGORITHM FOR EFFECTIVE RCCE MODELING OF NATURAL GAS IGNITION IN AIR
Proceedings of the ASME 2016 International Mechanical Engineering Congress & Exposition IMECE2016 November 11-17, 2016, Phoenix, Arizona, USA IMECE2016-65323 VALIDATION OF THE ASVDADD CONSTRAINT SELECTION
More informationApproximation of chemical reaction rates in turbulent combustion simulation
Approximation of chemical reaction rates in turbulent combustion simulation Lars Frank Große and Franz Joos * Helmut-Schmidt-University University of the Federal Armed Forces Hamburg - Laboratory of Turbo
More informationMEASUREMENT UNCERTAINTY AND SUMMARISING MONTE CARLO SAMPLES
XX IMEKO World Congress Metrology for Green Growth September 9 14, 212, Busan, Republic of Korea MEASUREMENT UNCERTAINTY AND SUMMARISING MONTE CARLO SAMPLES A B Forbes National Physical Laboratory, Teddington,
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY
CALIFORNIA INSTITUTE OF TECHNOLOGY EARTHQUAKE ENGINEERING RESEARCH LABORATORY NEW BAYESIAN UPDATING METHODOLOGY FOR MODEL VALIDATION AND ROBUST PREDICTIONS BASED ON DATA FROM HIERARCHICAL SUBSYSTEM TESTS
More informationREDIM reduced modeling of quenching at a cold inert wall with detailed transport and different mechanisms
26 th ICDERS July 3 th August 4 th, 217 Boston, MA, USA REDIM reduced modeling of quenching at a cold inert wall with detailed transport and different mechanisms Christina Strassacker, Viatcheslav Bykov,
More informationDevelopment and Optimization of H 2 /CO-H 2 Reaction Model using
Development and Optimization of H 2 /CO-H 2 Reaction Model using PrIMe A. Mirzayeva 1, N.A. Slavinskaya 1, U. Riedel 1, M. Frenklach 2, A. Packard 2, W. Li 2, J. Oreluk 2, A. Hedge 2 1 Institute of Combustion
More informationA comparison of the Bader Deuflhard and the Cash Karp Runge Kutta integrators for the GRI-MECH 3.0 model based on the chemical kinetics code Kintecus
1368 A comparison of the Bader Deuflhard and the Cash Karp Runge Kutta integrators for the GRI-MECH 3.0 model based on the chemical inetics code Kintecus James C. Ianni Vast Technologies Development, Inc.,
More informationRisk Estimation and Uncertainty Quantification by Markov Chain Monte Carlo Methods
Risk Estimation and Uncertainty Quantification by Markov Chain Monte Carlo Methods Konstantin Zuev Institute for Risk and Uncertainty University of Liverpool http://www.liv.ac.uk/risk-and-uncertainty/staff/k-zuev/
More informationDensity Estimation. Seungjin Choi
Density Estimation Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr http://mlg.postech.ac.kr/
More information2D Image Processing (Extended) Kalman and particle filter
2D Image Processing (Extended) Kalman and particle filter Prof. Didier Stricker Dr. Gabriele Bleser Kaiserlautern University http://ags.cs.uni-kl.de/ DFKI Deutsches Forschungszentrum für Künstliche Intelligenz
More informationDynamics of Excited Hydroxyl Radicals in Hydrogen Based Mixtures Behind Reflected Shock Waves. Supplemental material
Dynamics of Excited Hydroxyl Radicals in Hydrogen Based Mixtures Behind Reflected Shock Waves Proceedings of the Combustion Institute 34, 22 Supplemental material R. MÉVEL, S. PICHON, L. CATOIRE, N. CHAUMEIX,
More informationA BAYESIAN MATHEMATICAL STATISTICS PRIMER. José M. Bernardo Universitat de València, Spain
A BAYESIAN MATHEMATICAL STATISTICS PRIMER José M. Bernardo Universitat de València, Spain jose.m.bernardo@uv.es Bayesian Statistics is typically taught, if at all, after a prior exposure to frequentist
More informationSTAT 499/962 Topics in Statistics Bayesian Inference and Decision Theory Jan 2018, Handout 01
STAT 499/962 Topics in Statistics Bayesian Inference and Decision Theory Jan 2018, Handout 01 Nasser Sadeghkhani a.sadeghkhani@queensu.ca There are two main schools to statistical inference: 1-frequentist
More informationICES REPORT Model Misspecification and Plausibility
ICES REPORT 14-21 August 2014 Model Misspecification and Plausibility by Kathryn Farrell and J. Tinsley Odena The Institute for Computational Engineering and Sciences The University of Texas at Austin
More informationPremixed MILD Combustion of Propane in a Cylindrical. Furnace with a Single Jet Burner: Combustion and. Emission Characteristics
Premixed MILD Combustion of Propane in a Cylindrical Furnace with a Single Jet Burner: Combustion and Emission Characteristics Kin-Pang Cheong a, c, Guochang Wang a, Jianchun Mi a*, Bo Wang a, Rong Zhu
More informationUTSR Fellowship Presentation Gas Turbine Industrial Fellowship Program 2006
UTSR Fellowship Presentation Gas Turbine Industrial Fellowship Program 2006 Predicting Lean Blowout Using the Damkohler Number Matthew J. Bloxham, Brigham Young University Ingersoll Rand Energy Systems
More informationCSC 2541: Bayesian Methods for Machine Learning
CSC 2541: Bayesian Methods for Machine Learning Radford M. Neal, University of Toronto, 2011 Lecture 10 Alternatives to Monte Carlo Computation Since about 1990, Markov chain Monte Carlo has been the dominant
More informationResults of turbulent flame speed for H 2 -rich and syngas fuel mixtures measured
Results of turbulent flame speed for H 2 -rich and syngas fuel mixtures measured Deliverable 1.1.4 SEVENTH FRAMEWORK PROGRAMME FP7-ENERGY-2008-TREN-1 ENERGY-2008-6-CLEAN COAL TECHNOLOGIES Project Acronym:
More informationPlasma Assisted Reforming of Methane: Two Stage Perfectly Stirred Reactor (PSR) Simulation. L. Bromberg N. Alexeev.
PSFC/JA-05-12 Plasma Assisted Reforming of Methane: Two Stage Perfectly Stirred Reactor (PSR) Simulation L. Bromberg N. Alexeev August 25, 2005 Massachusetts Institute of Technology Plasma Science and
More informationExperimental study of the combustion properties of methane/hydrogen mixtures Gersen, Sander
University of Groningen Experimental study of the combustion properties of methane/hydrogen mixtures Gersen, Sander IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF)
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationFlamelet Analysis of Turbulent Combustion
Flamelet Analysis of Turbulent Combustion R.J.M. Bastiaans,2, S.M. Martin, H. Pitsch,J.A.vanOijen 2, and L.P.H. de Goey 2 Center for Turbulence Research, Stanford University, CA 9435, USA 2 Eindhoven University
More informationVariational Methods in Bayesian Deconvolution
PHYSTAT, SLAC, Stanford, California, September 8-, Variational Methods in Bayesian Deconvolution K. Zarb Adami Cavendish Laboratory, University of Cambridge, UK This paper gives an introduction to the
More informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More informationPredictive Engineering and Computational Sciences. Research Challenges in VUQ. Robert Moser. The University of Texas at Austin.
PECOS Predictive Engineering and Computational Sciences Research Challenges in VUQ Robert Moser The University of Texas at Austin March 2012 Robert Moser 1 / 9 Justifying Extrapolitive Predictions Need
More informationCALCULATION OF THE UPPER EXPLOSION LIMIT OF METHANE-AIR MIXTURES AT ELEVATED PRESSURES AND TEMPERATURES
CALCULATION OF THE UPPER EXPLOSION LIMIT OF METHANE-AIR MIXTURES AT ELEVATED PRESSURES AND TEMPERATURES F. Van den Schoor 1, F. Verplaetsen 2 and J. Berghmans 1 1 Katholieke Universiteit Leuven, Department
More informationExperimental study of the combustion properties of methane/hydrogen mixtures Gersen, Sander
University of Groningen Experimental study of the combustion properties of methane/hydrogen mixtures Gersen, Sander IMPORTANT NOTE: You are advised to consult the publisher's version (publisher's PDF)
More informationExpectation Propagation for Approximate Bayesian Inference
Expectation Propagation for Approximate Bayesian Inference José Miguel Hernández Lobato Universidad Autónoma de Madrid, Computer Science Department February 5, 2007 1/ 24 Bayesian Inference Inference Given
More informationA new Hierarchical Bayes approach to ensemble-variational data assimilation
A new Hierarchical Bayes approach to ensemble-variational data assimilation Michael Tsyrulnikov and Alexander Rakitko HydroMetCenter of Russia College Park, 20 Oct 2014 Michael Tsyrulnikov and Alexander
More informationStatistical Data Analysis Stat 3: p-values, parameter estimation
Statistical Data Analysis Stat 3: p-values, parameter estimation London Postgraduate Lectures on Particle Physics; University of London MSci course PH4515 Glen Cowan Physics Department Royal Holloway,
More informationLecture 2: From Linear Regression to Kalman Filter and Beyond
Lecture 2: From Linear Regression to Kalman Filter and Beyond January 18, 2017 Contents 1 Batch and Recursive Estimation 2 Towards Bayesian Filtering 3 Kalman Filter and Bayesian Filtering and Smoothing
More informationA wide range kinetic modelling study of laminar flame speeds of reference fuels and their mixtures
A wide range kinetic modelling study of laminar flame speeds of reference fuels and their mixtures A. Frassoldati, R. Grana, A. Cuoci, T. Faravelli, E. Ranzi Dipartimento di Chimica, Materiali e Ingegneria
More informationBayesian search for other Earths
Bayesian search for other Earths Low-mass planets orbiting nearby M dwarfs Mikko Tuomi University of Hertfordshire, Centre for Astrophysics Research Email: mikko.tuomi@utu.fi Presentation, 19.4.2013 1
More informationConfirmation of paper submission
www.flame-structure-2014.com Berlin Institute of Technology (TU Berlin) Prof. Dr. Frank Behrendt Fakultät III: Prozesswissenschaften, Institut für Energietechnik Chair Energy Process Engineering and Conversion
More informationParameter Estimation. William H. Jefferys University of Texas at Austin Parameter Estimation 7/26/05 1
Parameter Estimation William H. Jefferys University of Texas at Austin bill@bayesrules.net Parameter Estimation 7/26/05 1 Elements of Inference Inference problems contain two indispensable elements: Data
More informationIntroduction to Bayesian Statistics. James Swain University of Alabama in Huntsville ISEEM Department
Introduction to Bayesian Statistics James Swain University of Alabama in Huntsville ISEEM Department Author Introduction James J. Swain is Professor of Industrial and Systems Engineering Management at
More informationNon-Equilibrium Reaction Rates in Hydrogen Combustion
25 th ICDERS August 2 7, 25 Leeds, UK Non-Equilibrium Reaction Rates in Hydrogen Combustion Stephen J. Voelkel, Venkat Raman 2, Philip Varghese The University of Texas at Austin, Austin, TX 7872, USA 2
More informationSelf Adaptive Particle Filter
Self Adaptive Particle Filter Alvaro Soto Pontificia Universidad Catolica de Chile Department of Computer Science Vicuna Mackenna 4860 (143), Santiago 22, Chile asoto@ing.puc.cl Abstract The particle filter
More informationFrequentist-Bayesian Model Comparisons: A Simple Example
Frequentist-Bayesian Model Comparisons: A Simple Example Consider data that consist of a signal y with additive noise: Data vector (N elements): D = y + n The additive noise n has zero mean and diagonal
More informationPropagation of Uncertainties in Measurements: Generalized/ Fiducial Inference
Propagation of Uncertainties in Measurements: Generalized/ Fiducial Inference Jack Wang & Hari Iyer NIST, USA NMIJ-BIPM Workshop, AIST-Tsukuba, Japan, May 18-20, 2005 p. 1/31 Frameworks for Quantifying
More informationBayesian inference: what it means and why we care
Bayesian inference: what it means and why we care Robin J. Ryder Centre de Recherche en Mathématiques de la Décision Université Paris-Dauphine 6 November 2017 Mathematical Coffees Robin Ryder (Dauphine)
More informationIntroduction to Probabilistic Machine Learning
Introduction to Probabilistic Machine Learning Piyush Rai Dept. of CSE, IIT Kanpur (Mini-course 1) Nov 03, 2015 Piyush Rai (IIT Kanpur) Introduction to Probabilistic Machine Learning 1 Machine Learning
More informationStat 535 C - Statistical Computing & Monte Carlo Methods. Arnaud Doucet.
Stat 535 C - Statistical Computing & Monte Carlo Methods Arnaud Doucet Email: arnaud@cs.ubc.ca 1 Suggested Projects: www.cs.ubc.ca/~arnaud/projects.html First assignement on the web: capture/recapture.
More informationThermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study
Feb 12 th 2004 Thermal NO Predictions in Glass Furnaces: A Subgrid Scale Validation Study Padmabhushana R. Desam & Prof. Philip J. Smith CRSIM, University of Utah Salt lake city, UT-84112 18 th Annual
More informationApproximate Bayesian Computation: a simulation based approach to inference
Approximate Bayesian Computation: a simulation based approach to inference Richard Wilkinson Simon Tavaré 2 Department of Probability and Statistics University of Sheffield 2 Department of Applied Mathematics
More informationExpectation propagation for signal detection in flat-fading channels
Expectation propagation for signal detection in flat-fading channels Yuan Qi MIT Media Lab Cambridge, MA, 02139 USA yuanqi@media.mit.edu Thomas Minka CMU Statistics Department Pittsburgh, PA 15213 USA
More informationThe role of diffusion at shear layers in irregular detonations
The role of diffusion at shear layers in irregular detonations Marco Arienti 1 Joseph E. Shepherd 2 1 United Technologies Research Center, 411 Silver Lane, East Hartford, CT 06108 2 California Institute
More informationLecture 3. G. Cowan. Lecture 3 page 1. Lectures on Statistical Data Analysis
Lecture 3 1 Probability (90 min.) Definition, Bayes theorem, probability densities and their properties, catalogue of pdfs, Monte Carlo 2 Statistical tests (90 min.) general concepts, test statistics,
More informationRobotics. Lecture 4: Probabilistic Robotics. See course website for up to date information.
Robotics Lecture 4: Probabilistic Robotics See course website http://www.doc.ic.ac.uk/~ajd/robotics/ for up to date information. Andrew Davison Department of Computing Imperial College London Review: Sensors
More informationExpectation Propagation Algorithm
Expectation Propagation Algorithm 1 Shuang Wang School of Electrical and Computer Engineering University of Oklahoma, Tulsa, OK, 74135 Email: {shuangwang}@ou.edu This note contains three parts. First,
More informationBayesian System Identification based on Hierarchical Sparse Bayesian Learning and Gibbs Sampling with Application to Structural Damage Assessment
Bayesian System Identification based on Hierarchical Sparse Bayesian Learning and Gibbs Sampling with Application to Structural Damage Assessment Yong Huang a,b, James L. Beck b,* and Hui Li a a Key Lab
More informationAt A Glance. UQ16 Mobile App.
At A Glance UQ16 Mobile App Scan the QR code with any QR reader and download the TripBuilder EventMobile app to your iphone, ipad, itouch or Android mobile device. To access the app or the HTML 5 version,
More informationNumerical Simulation of Entropy Generation in Hydrogen Enriched Swirl Stabilized Combustion
Saqr & Wahid CFD Letters Vol. 5(1) 13 www.cfdl.issres.net Vol. 5 (1) March 13 Numerical Simulation of Entropy Generation in Hydrogen Enriched Swirl Stabilized Combustion Khalid M. Saqr 1,* and Mazlan A.
More informationParametric Models. Dr. Shuang LIANG. School of Software Engineering TongJi University Fall, 2012
Parametric Models Dr. Shuang LIANG School of Software Engineering TongJi University Fall, 2012 Today s Topics Maximum Likelihood Estimation Bayesian Density Estimation Today s Topics Maximum Likelihood
More informationIntroduction to Smoothing spline ANOVA models (metamodelling)
Introduction to Smoothing spline ANOVA models (metamodelling) M. Ratto DYNARE Summer School, Paris, June 215. Joint Research Centre www.jrc.ec.europa.eu Serving society Stimulating innovation Supporting
More informationBlack-box α-divergence Minimization
Black-box α-divergence Minimization José Miguel Hernández-Lobato, Yingzhen Li, Daniel Hernández-Lobato, Thang Bui, Richard Turner, Harvard University, University of Cambridge, Universidad Autónoma de Madrid.
More informationCorrelation Effects in the Uncertainty Estimation of Two-Pressure Humidity Generators
Journal of Physics: Conference Series PAPER OPEN ACCESS Correlation Effects in the Uncertainty Estimation of Two-Pressure Humidity Generators To cite this article: S A Carvajal and C A Sánchez 08 J. Phys.:
More informationReview of the role of uncertainties in room acoustics
Review of the role of uncertainties in room acoustics Ralph T. Muehleisen, Ph.D. PE, FASA, INCE Board Certified Principal Building Scientist and BEDTR Technical Lead Division of Decision and Information
More informationRevision of the Guide to the expression of uncertainty in measurement impact on national metrology institutes and industry
Revision of the Guide to the expression of uncertainty in measurement impact on national metrology institutes and industry Maurice Cox/Peter Harris National Physical Laboratory, Teddington, UK CCRI BIPM,
More informationBayesian Calibration of Simulators with Structured Discretization Uncertainty
Bayesian Calibration of Simulators with Structured Discretization Uncertainty Oksana A. Chkrebtii Department of Statistics, The Ohio State University Joint work with Matthew T. Pratola (Statistics, The
More informationStatistical Methods in Particle Physics
Statistical Methods in Particle Physics Lecture 11 January 7, 2013 Silvia Masciocchi, GSI Darmstadt s.masciocchi@gsi.de Winter Semester 2012 / 13 Outline How to communicate the statistical uncertainty
More informationModeling ion and electron profiles in methane-oxygen counterflow diffusion flames
Abstract 8 th U. S. National Combustion Meeting Organized by the Western States Section of the Combustion Institute and hosted by the University of Utah May 19-22, 2013 Modeling ion and electron profiles
More informationProfessor David B. Stephenson
rand Challenges in Probabilistic Climate Prediction Professor David B. Stephenson Isaac Newton Workshop on Probabilistic Climate Prediction D. Stephenson, M. Collins, J. Rougier, and R. Chandler University
More informationA Fresh Look at the Bayes Theorem from Information Theory
A Fresh Look at the Bayes Theorem from Information Theory Tan Bui-Thanh Computational Engineering and Optimization (CEO) Group Department of Aerospace Engineering and Engineering Mechanics Institute for
More informationCLASS NOTES Models, Algorithms and Data: Introduction to computing 2018
CLASS NOTES Models, Algorithms and Data: Introduction to computing 208 Petros Koumoutsakos, Jens Honore Walther (Last update: June, 208) IMPORTANT DISCLAIMERS. REFERENCES: Much of the material (ideas,
More informationApproximate Bayesian Computation by Subset Simulation for model selection in dynamical systems
Available online at www.sciencedirect.com ScienceDirect Procedia Engineering 199 (2017) 1056 1061 X International Conference on Structural Dynamics, EURODYN 2017 Approximate Bayesian Computation by Subset
More informationProbabilities for climate projections
Probabilities for climate projections Claudia Tebaldi, Reinhard Furrer Linda Mearns, Doug Nychka National Center for Atmospheric Research Richard Smith - UNC-Chapel Hill Steve Sain - CU-Denver Statistical
More informationParametric Techniques Lecture 3
Parametric Techniques Lecture 3 Jason Corso SUNY at Buffalo 22 January 2009 J. Corso (SUNY at Buffalo) Parametric Techniques Lecture 3 22 January 2009 1 / 39 Introduction In Lecture 2, we learned how to
More informationAnalysis of NO-Formation for Rich / Lean - Staged Combustion
1 Analysis of NO-Formation for Rich / Lean - Staged Combustion P.Frank (a), Y.Tan (a), P.Griebel (b), H.Nannen (b), H.Eickhoff (b) Deutsche Forschungsanstalt für Luft-und Raumfahrt : (a) Institut für Physikalische
More informationIntroduction: MLE, MAP, Bayesian reasoning (28/8/13)
STA561: Probabilistic machine learning Introduction: MLE, MAP, Bayesian reasoning (28/8/13) Lecturer: Barbara Engelhardt Scribes: K. Ulrich, J. Subramanian, N. Raval, J. O Hollaren 1 Classifiers In this
More informationDevelopment of One-Step Chemistry Models for Flame and Ignition Simulation
Development of One-Step Chemistry Models for Flame and Ignition Simulation S.P.M. Bane, J.L. Ziegler, and J.E. Shepherd Graduate Aerospace Laboratories California Institute of Technology Pasadena, CA 91125
More informationAnnex xxx General guidelines for the assessment of uncertainties in solar thermal systems performance testing
Annex xxx General guidelines for the assessment of uncertainties in solar thermal systems performance testing 1 Introduction The expected energy output presents the most important quantity, within those
More informationSkeletal Kinetic Mechanism of Methane Oxidation for High Pressures and Temperatures
7 TH EUROPEAN CONFERENCE FOR AERONAUTICS AND SPACE SCIENCES (EUCASS) Skeletal Kinetic Mechanism of Methane Oxidation for High Pressures and Temperatures Victor P. Zhukov and Alan F. Kong Institute of Space
More informationEfficient Likelihood-Free Inference
Efficient Likelihood-Free Inference Michael Gutmann http://homepages.inf.ed.ac.uk/mgutmann Institute for Adaptive and Neural Computation School of Informatics, University of Edinburgh 8th November 2017
More informationarxiv: v1 [physics.chem-ph] 26 Apr 2011
Assessment of existing H 2 /O 2 chemical reaction mechanisms at reheat gas turbine conditions Torleif Weydahl a,, Madhavan Poyyapakkam b, Morten Seljeskog a, Nils Erland L. Haugen a arxiv:1104.4925v1 [physics.chem-ph]
More informationPressure dependent mechanism for H/O/C(1) chemistry
Paper # P-22 Topic: Kinetics 5 th US ombustion Meeting rganized by the Western States Section of the ombustion Institute and osted by the University of alifornia at San Diego March 25-28, 2007. Pressure
More informationMeasurement Uncertainty of Dew-Point Temperature in a Two-Pressure Humidity Generator
Int J Thermophys (2012) 33:1568 1582 DOI 10.1007/s10765-011-1005-z Measurement Uncertainty of Dew-Point Temperature in a Two-Pressure Humidity Generator L. Lages Martins A. Silva Ribeiro J. Alves e Sousa
More informationIntroduction to Bayesian Data Analysis
Introduction to Bayesian Data Analysis Phil Gregory University of British Columbia March 2010 Hardback (ISBN-10: 052184150X ISBN-13: 9780521841504) Resources and solutions This title has free Mathematica
More informationHierarchical sparse Bayesian learning for structural health monitoring. with incomplete modal data
Hierarchical sparse Bayesian learning for structural health monitoring with incomplete modal data Yong Huang and James L. Beck* Division of Engineering and Applied Science, California Institute of Technology,
More informationVariational Principal Components
Variational Principal Components Christopher M. Bishop Microsoft Research 7 J. J. Thomson Avenue, Cambridge, CB3 0FB, U.K. cmbishop@microsoft.com http://research.microsoft.com/ cmbishop In Proceedings
More informationMCMC for big data. Geir Storvik. BigInsight lunch - May Geir Storvik MCMC for big data BigInsight lunch - May / 17
MCMC for big data Geir Storvik BigInsight lunch - May 2 2018 Geir Storvik MCMC for big data BigInsight lunch - May 2 2018 1 / 17 Outline Why ordinary MCMC is not scalable Different approaches for making
More information13: Variational inference II
10-708: Probabilistic Graphical Models, Spring 2015 13: Variational inference II Lecturer: Eric P. Xing Scribes: Ronghuo Zheng, Zhiting Hu, Yuntian Deng 1 Introduction We started to talk about variational
More informationEVALUATING SYMMETRIC INFORMATION GAP BETWEEN DYNAMICAL SYSTEMS USING PARTICLE FILTER
EVALUATING SYMMETRIC INFORMATION GAP BETWEEN DYNAMICAL SYSTEMS USING PARTICLE FILTER Zhen Zhen 1, Jun Young Lee 2, and Abdus Saboor 3 1 Mingde College, Guizhou University, China zhenz2000@21cn.com 2 Department
More informationThe Jackknife-Like Method for Assessing Uncertainty of Point Estimates for Bayesian Estimation in a Finite Gaussian Mixture Model
Thai Journal of Mathematics : 45 58 Special Issue: Annual Meeting in Mathematics 207 http://thaijmath.in.cmu.ac.th ISSN 686-0209 The Jackknife-Like Method for Assessing Uncertainty of Point Estimates for
More informationVCMC: Variational Consensus Monte Carlo
VCMC: Variational Consensus Monte Carlo Maxim Rabinovich, Elaine Angelino, Michael I. Jordan Berkeley Vision and Learning Center September 22, 2015 probabilistic models! sky fog bridge water grass object
More informationBAYESIAN ESTIMATION OF UNKNOWN PARAMETERS OVER NETWORKS
BAYESIAN ESTIMATION OF UNKNOWN PARAMETERS OVER NETWORKS Petar M. Djurić Dept. of Electrical & Computer Engineering Stony Brook University Stony Brook, NY 11794, USA e-mail: petar.djuric@stonybrook.edu
More informationBayesian Learning. HT2015: SC4 Statistical Data Mining and Machine Learning. Maximum Likelihood Principle. The Bayesian Learning Framework
HT5: SC4 Statistical Data Mining and Machine Learning Dino Sejdinovic Department of Statistics Oxford http://www.stats.ox.ac.uk/~sejdinov/sdmml.html Maximum Likelihood Principle A generative model for
More informationLecture 13 : Variational Inference: Mean Field Approximation
10-708: Probabilistic Graphical Models 10-708, Spring 2017 Lecture 13 : Variational Inference: Mean Field Approximation Lecturer: Willie Neiswanger Scribes: Xupeng Tong, Minxing Liu 1 Problem Setup 1.1
More informationChemical Kinetics of Ethanol Oxidation. Tonawanda, NY 14150, USA. Princeton University Princeton, NJ 08544, USA
Chemical Kinetics of Ethanol xidation Juan Li 1, Andrei Kazakov 2, Frederick L. Dryer 2* 1 Praxair, Inc. Tonawanda, NY 1415, USA 2 Department of Mechanical and Aerospace Engineering Princeton University
More informationIntroduction to Bayesian Methods
Introduction to Bayesian Methods Jessi Cisewski Department of Statistics Yale University Sagan Summer Workshop 2016 Our goal: introduction to Bayesian methods Likelihoods Priors: conjugate priors, non-informative
More informationAn introduction to Bayesian statistics and model calibration and a host of related topics
An introduction to Bayesian statistics and model calibration and a host of related topics Derek Bingham Statistics and Actuarial Science Simon Fraser University Cast of thousands have participated in the
More informationON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION
ON THE FAILURE RATE ESTIMATION OF THE INVERSE GAUSSIAN DISTRIBUTION ZHENLINYANGandRONNIET.C.LEE Department of Statistics and Applied Probability, National University of Singapore, 3 Science Drive 2, Singapore
More information