Abstract The influence of several potential error sources and non-ideal experimental effects on the accuracy of a wavelet-based optical flow velocimetry (wOFV) method when applied to tracer particle images is evaluated using data from a series of synthetic flows. Out-of-plane particle displacements, severe image noise, laser sheet thickness reduction, and image intensity non-uniformity are shown to decrease the accuracy of wOFV in a similar manner to correlation-based particle image velocimetry (PIV). For the error sources tested, wOFV displays a similar or slightly increased sensitivity compared to PIV, but the wOFV results are still more accurate than PIV when the magnitude of the non-ideal effects remain within expected experimental bounds. For the majority of test cases, the results are significantly improved by using image pre-processing filters and the magnitude of improvement is consistent between wOFV and PIV. Flow divergence does not appear to have an appreciable effect on the accuracy of wOFV velocity estimation, even though the underlying fluid transport equation on which wOFV is based implicitly assumes that the motion is divergence-free. This is a significant finding for the broader applicability of planar velocimetry measurements using wOFV. Finally, it is noted that the accuracy of wOFV is not reduced notably in regions of the image between tracer particles, as long as the overall seeding density is not too sparse i.e. below 0.02 particles per pixel. This explicitly demonstrates that wOFV (when applied to particle images) yields an accurate whole field measurement, and not only at or adjacent to the discrete particle locations.
more »
« less
This content will become publicly available on May 1, 2026
A Bayesian approach to locally varying regularization in optical flow velocimetry
Optical flow velocimetry (OFV) is a method for determining dense and accurate velocity fields from a pair of particle images by solving the classical optical flow problem. However, this is an ill-posed inverse problem, which generally entails minimizing a weighted sum of two terms–fidelity and regularization–and the weights in the sum are parameters that require manual tuning based on the properties of both the flow and the particle images. This manual tuning has historically been a consistent challenge that has limited the general applicability of OFV for experimental data, as the calculated velocity field is sensitive to the value of the weights. This work proposes a hierarchical model for the weighting parameters in the framework of a maximum a posteriori-based Bayesian optimization approach. The method replaces the classical Lagrange multiplier weighting parameter with a new, less-sensitive parameter that can be automatically predetermined from experimental images. The resulting method is tested on three different synthetic particle image velocimetry (PIV) datasets and on experimental particle images. The method is found to be capable of self-adjusting the local weights of the optimization process in real-time while simultaneously determining the velocity field, leading to an optimally regularized estimate of the velocity field without requiring any dataset-specific manual tuning of the parameters. The presented approach is the first truly general, parameter-free optical flow method for particle image velocimetry (PIV) images. The developed method is freely available as a part of the PIVlab package.
more »
« less
- Award ID(s):
- 2204618
- PAR ID:
- 10630131
- Publisher / Repository:
- AIP Publishing
- Date Published:
- Journal Name:
- Physics of Fluids
- Volume:
- 37
- Issue:
- 5
- ISSN:
- 1070-6631
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract High fidelity near-wall velocity measurements in wall bounded fluid flows continue to pose a challenge and the resulting limitations on available experimental data cloud our understanding of the near-wall velocity behavior in turbulent boundary layers. One of the challenges is the spatial averaging and limited spatial resolution inherent to cross-correlation-based particle image velocimetry (PIV) methods. To circumvent this difficulty, we implement an explicit no-slip boundary condition in a wavelet-based optical flow velocimetry (wOFV) method. It is found that the no-slip boundary condition on the velocity field can be implemented in wOFV by transforming the constraint to the wavelet domain through a series of algebraic linear transformations, which are formulated in terms of the known wavelet filter matrices, and then satisfying the resulting constraint on the wavelet coefficients using constrained optimization for the optical flow functional minimization. The developed method is then used to study the classical problem of a turbulent channel flow using synthetic data from a direct numerical simulation (DNS) and experimental particle image data from a zero pressure gradient, high Reynolds number turbulent boundary layer. The results obtained by successfully implementing the no-slip boundary condition are compared to velocity measurements from wOFV without the no-slip condition and to a commercial PIV code, using the velocity from the DNS as ground truth. It is found that wOFV with the no-slip condition successfully resolves the near-wall profile with enhanced accuracy compared to the other velocimetry methods, as well as other derived quantities such as wall shear and turbulent intensity, without sacrificing accuracy away from the wall, leading to state of the art measurements in the region of the turbulent boundary layer when applied to experimental particle images.more » « less
-
Optical distortion caused by changes in the refractive index of fluid flow is a common issue in flow visualization using techniques, such as particle image velocimetry (PIV). In thermally driven convection, this distortion can severely interfere with PIV results due to the ubiquitous density and, therefore, refractive index heterogeneity in the fluid. The distortion also varies spatially and temporally, adding to the challenge. We propose a composite filter, the shadow-affected PIV region filter, which combines a series of conventional image filters to address this issue, focusing on optical distortion of thermal plumes in laminar flow. We verify the effectiveness of the filter using both synthetic particle images created from ray tracing and real particle images from the laboratory. For the first time, we effectively mitigate the optical distortion from plumes while preserving the in-plane plume velocity and overall flow pattern, with the PIV data alone. Our filter is efficient and does not require additional measurements, expensive ray tracing, or a large dataset to begin with. It can be extended to separate the flow field and the effect of optical distortion in other fluid experiments when the two components are visually distinct. Additionally, this filter can serve as a baseline algorithm for comparison when developing more advanced methods like neural networks.more » « less
-
A new method for fluid–structure interaction (FSI) diagnostics to simultaneously capture time-resolved three-dimensional, three-component (3D3C) velocity fields and structural deformations using a single light field camera is presented. A light field camera encodes both spatial and angular information of light rays collected by a conventional imaging lens that allows for the 3D reconstruction of a scene from a single image. Building upon this capability, a light field fluid–structure interaction (LF FSI) methodology is developed with a focus on experimental scenarios with low optical access. Proper orthogonal decomposition (POD) is used to separate particle and surface information contained in the same image. A correlation-based depth estimation technique is introduced to reconstruct instantaneous surface positions from the disparity between angular perspectives and conventional particle image velocimetry (PIV) is used for flow field reconstruction. Validation of the methodology is achieved using synthetic images of simultaneously moving flat plates and a vortex ring with a small increase in uncertainty under ~0.5 microlenses observed in both flow and structure measurement compared to independent measurements. The method is experimentally verified using a flat plate translating along the camera’s optical axis in a flow field with varying particle concentrations. Finally, simultaneous reconstructions of the flow field and surface shape around a flexible membrane are presented, with the surface reconstruction further validated using simultaneously captured stereo images. The findings indicate that the LF FSI methodology provides a new capability to simultaneously measure large-scale flow characteristics and structural deformations using a single camera.more » « less
-
Three-dimensional and three-component (3D3C) velocity measurements have long been desired to resolve the 3D spatial structures of turbulent flows. Recent advancements have demonstrated tomographic particle image velocimetry (tomo-PIV) as a powerful technique to enable such measurements. The existing tomo-PIV technique obtains 3D3C velocity field by cross-correlating two frames of 3D tomographic reconstructions of the seeding particles. A most important issue in 3D3C velocity measurement involves uncertainty, as the derivatives of the measurements are usually of ultimate interest and uncertainties are amplified when calculating derivatives. To reduce the uncertainties of 3D3C velocity measurements, this work developed a regularized tomo-PIV method. The new method was demonstrated to enhance accuracy significantly by incorporating the conservation of mass into the tomo-PIV process. The new method was demonstrated and validated both experimentally and numerically. The results illustrated that the new method was able to enhance the accuracy of 3D3C velocity measurements by 40%–50% in terms of velocity magnitude and by 0.6°–1.1° in terms of velocity orientation, compared to the existing tomo-PIV technique. These improvements brought about by the new method are expected to expand the application of tomo-PIV techniques when accuracy and quantitative 3D flow properties are required.more » « less
An official website of the United States government
