We present an uncertainty quantification methodology for density estimation from Background Oriented Schlieren (BOS) measurements, in order to provide local, instantaneous, a-posteriori uncertainty bounds on each density measurement in the field of view. Displacement uncertainty quantification algorithms from cross-correlation based Particle Image Velocimetry (PIV) are used to estimate the uncertainty in the dot pattern displacements obtained from cross-correlation for BOS and assess their feasibility. In order to propagate the displacement uncertainty through the density integration procedure, we also develop a novel methodology via the Poisson solver using sparse linear operators. Testing the method using synthetic images of a Gaussian density field showed agreement between the propagated density uncertainties and the true uncertainty. Subsequently the methodology is experimentally demonstrated for supersonic flow over a wedge, showing that regions with sharp changes in density lead to an increase in density uncertainty throughout the field of view, even in regions without these sharp changes. The uncertainty propagation is influenced by the density integration scheme, and for the Poisson solver the density uncertainty increases monotonically on moving away from the regions where the Dirichlet boundary conditions are specified.

We present an image generation methodology based on ray tracing that can be used to render realistic images of particle image velocimetry (PIV) and background oriented schlieren (BOS) experiments in the presence of density/refractive index gradients. This methodology enables the simulation of aero-thermodynamics experiments for experiment design, error, and uncertainty analysis. Images are generated by emanating light rays from the particles or dot pattern, and propagating them through the density gradient field and the optical elements, up to the camera sensor. The rendered images are realistic, and can replicate the features of a given experimental setup, like optical aberrations and perspective effects, which can be deliberately introduced for error analysis. We demonstrate this methodology by simulating a BOS experiment with a known density field obtained from direct numerical simulations of homogeneous buoyancy driven turbulence, and comparing the light ray displacements from ray tracing to results from BOS theory. The light ray displacements show good agreement with the reference data. This methodology provides a framework for further development of simulation tools for use in experiment design and development of image analysis tools for PIV and BOS applications. An implementation of the proposed methodology in a Python-CUDA program is made available as an open source software for researchers.

We propose a dot-tracking methodology for processing Background Oriented Schlieren (BOS) images. The method significantly improves the accuracy, precision and spatial resolution compared to conventional cross-correlation algorithms. Our methodology utilizes the prior information about the dot pattern such as the location, size and number of dots to provide near 100% yield even for high dot densities (20 dots/32x32 pix.) and is robust to image noise. We also propose an improvement to the displacement estimation step in the tracking process, especially for noisy images, using a "correlation correction", whereby we combine the spatial resolution benefit of the tracking method and the smoothing property of the correlation method to increase the dynamic range of the overall measurement process. We evaluate the performance of the method with synthetic BOS images of buoyancy driven turbulence rendered using ray tracing simulations, and experimental images of flow in the exit plane of a converging-diverging nozzle.

No abstract

The cooling process associated with the flow induced by a spark plasma discharge generated between a pair of electrodes is measured using stereoscopic particle image velocimetry (S-PIV) and background oriented schlieren (BOS). Density measurements show that the hot gas kernel initially cools fast by convective cooling, followed by a slower cooling process. The cooling rates during the fast regime range from being 2 to 10 times those in the slower regime. An analytical model is developed to relate the cooling observed in the fast regime from BOS, to the total entrainment of cold ambient fluid per unit volume of the hot gas kernel, measured from S-PIV. The model calculates the cooling ratio to characterize the cooling process and shows that the cooling ratio estimated from the density measurements are in close agreement with those calculated from the entrainment. These measurements represent the first ever quantitative density and velocity measurements of the flow induced by a spark discharge and reveal the role of entrainment on the cooling of the hot gas kernel. These results underscore that convective cooling of the hot gas kernel, in the fast regime, leads to approximately 50% of the cooling and occurs within the first millisecond of the induced flow.

No abstract

We theoretically analyze the effect of density/refractive-index gradients on the measurement precision of Volumetric Particle Tracking Velocimetry (V-PTV) and Background Oriented Schlieren (BOS) experiments by deriving the Cramer-Rao lower bound (CRLB) for the 2D centroid estimation process. A model is derived for the diffraction limited image of a particle or dot viewed through a medium containing density gradients that includes the effect of various experimental parameters such as the particle depth, viewing direction and f-number. Using the model we show that non-linearities in the density gradient field lead to blurring of the particle/dot image. This blurring amplifies the effect of image noise on the centroid estimation process, leading to an increase in the CRLB and a decrease in the measurement precision. The ratio of position uncertainties of a dot in the reference and gradient images is a function of the ratio of the dot diameters and dot intensities. We term this parameter the Amplification Ratio (AR), and we propose a methodology for estimating the position uncertainties in tracking-based BOS measurements. The theoretical predictions of the particle/dot position estimation variance from the CRLB are compared to ray tracing simulations with good agreement. The uncertainty amplification is also demonstrated on experimental BOS images of flow induced by a spark discharge, where we show that regions of high amplification ratio correspond to regions of density gradients. This analysis elucidates the dependence of the position error on density and refractiveindex gradient induced distortion parameters, provides a methodology for accounting its effect on uncertainty quantification and provides a framework for optimizing experiment design.

We propose an improved density integration methodology for Background Oriented Schlieren (BOS) measurements that overcomes the noise sensitivity of the commonly used Poisson solver. The method employs a weighted least-squares (WLS) optimization of the 2D integration of the density gradient field by solving an over-determined system of equations. Weights are assigned to the grid points based on density gradient uncertainties to ensure that a less reliable measurement point has less effect on the integration procedure. Synthetic image analysis with a Gaussian density field shows that WLS constrains the propagation of random error and reduces it by 80% in comparison to Poisson for the highest noise level. Using WLS with experimental BOS measurements of flow induced by a spark plasma discharge show a 30% reduction in density uncertainty in comparison to Poisson, thereby increasing the overall precision of the BOS density measurements.