Paul J. Hadwin

Bayesian approach to the design of chemical species tomography experiments.

Reconstruction accuracy in chemical species tomography depends strongly on the arrangement of optical paths transecting the imaging domain. Optimizing the path arrangement requires a scheme that can predict the quality of a proposed arrangement prior to measurement. This paper presents a new Bayesian method for scoring path arrangements based on the estimated a posteriori covariance matrix. This technique focuses on defining an objective function that incorporates the same a priori information about the flow needed to carry out limited data tomography. Constrained and unconstrained path optimization studies verify the predictive capabilities of the objective function, and that superior reconstruction quality is obtained with optimized path arrangements.

General error model for analysis of laser-induced incandescence signals.

This paper presents a novel error model for TiRe-LII signals and illustrates how the model can be used to diagnose a detection system, quantify uncertainties in TiRe-LII, and characterize fluctuations in the measured process. Noise in a single TiRe-LII measurement shot obeys a Poisson-Gaussian noise model. Variation in the aerosol results in shot-to-shot fluctuations in the measured signals. These fluctuations induce a quadratic relationship between the mean and variance of a set of signals. We show how this model can elucidate aspects of the measurement system and fundamental properties of the aerosol, by comparing the noise model to four sets of experimental data.

Measurement-based meshing, basis selection, and prior assignment in chemical species tomography

Gas distributions imaged by chemical species tomography (CST) vary in quality due to the discretization scheme, arrangement of optical paths, errors in the measurement model, and prior information included in reconstruction. There is currently no mathematically-rigorous framework for comparing the finite bases available to discretize a CST domain. Following from the Bayesian formulation of tomographic inversion, we show that Bayesian model selection can identify the mesh density, mode of interpolation, and prior information best-suited to reconstruct a set of measurement data. We validate this procedure with a simulated CST experiment, and generate accurate reconstructions despite limited measurement information. The flow field is represented using the finite element method, and Bayesian model selection is used to choose between three forms of polynomial support for a range of mesh resolutions, as well as four priors. We show that the model likelihood of Bayesian model selection is a good predictor of reconstruction accuracy.

Bayesian inference of tissue properties from glottal area waveforms using a 2D finite element model

Bayesian inference has recently been demonstrated to be effective in estimating stationary and non-stationary reduced-order vocal fold model parameters, along with the associated levels of uncertainty, from simulated glottal area waveform measures (Hadwin et al., 2016, 2017). In these studies, the fitting model was a three mass body-cover model with two degrees of freedom in the cover layer. While demonstrative, restricting the fitting model to two degrees of freedom assumes a priori that this is sufficient to capture salient vocal fold dynamics, thus limiting future clinical applicability. To overcome this, we employ Bayesian inference to directly estimate tissue properties of a two-dimensional (2D) finite element vocal fold model from glottal area waveforms generated by both numerical simulations and recorded videos of synthetic vocal fold oscillations. We demonstrate that the 2D finite element model is not only capable of producing meaningful estimates with reasonable uncertainties, but is also capable...

Predicting the heat of vaporization of iron at high temperatures using time-resolved laser-induced incandescence and Bayesian model selection

Competing theories have been proposed to account for how the latent heat of vaporization of liquid iron varies with temperature, but experimental confirmation remains elusive, particularly at high temperatures. We propose time-resolved laser-induced incandescence measurements on iron nanoparticles combined with Bayesian model plausibility, as a novel method for evaluating these relationships. Our approach scores the explanatory power of candidate models, accounting for parameter uncertainty, model complexity, measurement noise, and goodness-of-fit. The approach is first validated with simulated data and then applied to experimental data for iron nanoparticles in argon. Our results justify the use of Romans equation to account for the temperature dependence of the latent heat of vaporization of liquid iron.

An extended Kalman filter approach to non-stationary Bayesian estimation of reduced-order vocal fold model parameters

The Bayesian framework for parameter inference provides a basis from which subject-specific reduced-order vocal fold models can be generated. Previously, it has been shown that a particle filter technique is capable of producing estimates and associated credibility intervals of time-varying reduced-order vocal fold model parameters. However, the particle filter approach is difficult to implement and has a high computational cost, which can be barriers to clinical adoption. This work presents an alternative estimation strategy based upon Kalman filtering aimed at reducing the computational cost of subject-specific model development. The robustness of this approach to Gaussian and non-Gaussian noise is discussed. The extended Kalman filter (EKF) approach is found to perform very well in comparison with the particle filter technique at dramatically lower computational cost. Based upon the test cases explored, the EKF is comparable in terms of accuracy to the particle filter technique when greater than 6000 particles are employed; if less particles are employed, the EKF actually performs better. For comparable levels of accuracy, the solution time is reduced by 2 orders of magnitude when employing the EKF. By virtue of the approximations used in the EKF, however, the credibility intervals tend to be slightly underpredicted.