-
Unbalanced optimal transport for stochastic particle tracking
Authors:
Kairui Hao,
Atharva Hans,
Pavlos Vlachos,
Ilias Bilionis
Abstract:
Non-invasive flow measurement techniques, such as particle tracking velocimetry, resolve 3D velocity fields by pairing tracer particle positions in successive time steps. These trajectories are crucial for evaluating physical quantities like vorticity, shear stress, pressure, and coherent structures. Traditional approaches deterministically reconstruct particle positions and extract particle track…
▽ More
Non-invasive flow measurement techniques, such as particle tracking velocimetry, resolve 3D velocity fields by pairing tracer particle positions in successive time steps. These trajectories are crucial for evaluating physical quantities like vorticity, shear stress, pressure, and coherent structures. Traditional approaches deterministically reconstruct particle positions and extract particle tracks using tracking algorithms. However, reliable track estimation is challenging due to measurement noise caused by high particle density, particle image overlap, and falsely reconstructed 3D particle positions. To overcome this challenge, probabilistic approaches quantify the epistemic uncertainty in particle positions, typically using a Gaussian probability distribution. However, the standard deterministic tracking algorithms relying on nearest-neighbor search do not directly extend to the probabilistic setting. Moreover, such algorithms do not necessarily find globally consistent solutions robust to reconstruction errors. This paper aims to develop a globally consistent nearest-neighborhood algorithm that robustly extracts stochastic particle tracks from the reconstructed Gaussian particle distributions in all frames. Our tracking algorithm relies on the unbalanced optimal transport theory in the metric space of Gaussian measures. Specifically, we optimize a binary transport plan for efficiently moving the Gaussian distributions of reconstructed particle positions between time frames. We achieve this by computing the partial Wasserstein distance in the metric space of Gaussian measures. Our tracking algorithm is robust to position reconstruction errors since it automatically detects the number of particles that should be matched through hyperparameter optimization. Finally, we validate our method using an in vitro flow experiment using a 3D-printed cerebral aneurysm.
△ Less
Submitted 5 July, 2024;
originally announced July 2024.
-
Shock generated vorticity in spark discharges
Authors:
Bhavini Singh,
Lalit K. Rajendran,
Pavlos P. Vlachos,
Sally P. M. Bane
Abstract:
Spark discharges induce a complex flow field consisting of a shock wave at early times (~ 1 $μ$ s), a pair of vortex rings, and a hot gas kernel. The vortex rings entrain ambient gas into the hot gas kernel and control its cooling and expansion. In this work, we investigate the shock wave's contribution in producing the vortex ring vorticity. We analyze high-speed (700 kHz) schlieren images of the…
▽ More
Spark discharges induce a complex flow field consisting of a shock wave at early times (~ 1 $μ$ s), a pair of vortex rings, and a hot gas kernel. The vortex rings entrain ambient gas into the hot gas kernel and control its cooling and expansion. In this work, we investigate the shock wave's contribution in producing the vortex ring vorticity. We analyze high-speed (700 kHz) schlieren images of the shock wave for a range of electrical energies to measure the shock properties and estimate shock velocity and curvature. These measurements are combined with a model to calculate the vorticity generated. The measurements show that the highest vorticity is generated near the peak shock curvature location, and the shock curvature and strength increase with electrical energy deposited. A comparison of the vorticity estimated from the model to vorticity measurements from stereoscopic particle image velocimetry shows the results to be statistically equivalent. This suggests that the shock curvature and velocity contribute to the vortex rings induced by spark discharges.
△ Less
Submitted 19 March, 2021; v1 submitted 18 January, 2021;
originally announced January 2021.
-
Meta-Uncertainty for Particle Image Velocimetry
Authors:
Lalit K. Rajendran,
Sayantan Bhattacharya,
Sally P. M. Bane,
Pavlos P. Vlachos
Abstract:
Uncertainty quantification for Particle Image Velocimetry (PIV) is critical for comparing flow fields with Computational Fluid Dynamics (CFD) results, and model design and validation. However, PIV features a complex measurement chain with coupled, non-linear error sources, and quantifying the uncertainty is challenging. Multiple assessments show that none of the current methods can reliably measur…
▽ More
Uncertainty quantification for Particle Image Velocimetry (PIV) is critical for comparing flow fields with Computational Fluid Dynamics (CFD) results, and model design and validation. However, PIV features a complex measurement chain with coupled, non-linear error sources, and quantifying the uncertainty is challenging. Multiple assessments show that none of the current methods can reliably measure the actual uncertainty across a wide range of experiments. Because the current methods differ in assumptions regarding the measurement process and calculation procedures, it is not clear which method is best to use for an experiment. To address this issue, we propose a method to estimate an uncertainty method's sensitivity and reliability, termed the Meta-Uncertainty. The novel approach is automated, local, and instantaneous, and based on perturbation of the recorded particle images. We developed an image perturbation scheme based on adding random unmatched particles to the interrogation window pair considering the signal-to-noise (SNR) of the correlation plane. Each uncertainty scheme's response to several trials of random particle addition is used to estimate a reliability metric, defined as the rate of change of the inter-quartile range (IQR) of the uncertainties with increasing levels of particle addition. We also propose applying the meta-uncertainty as a weighting metric to combine uncertainty estimates from individual schemes, based on ideas from the consensus forecasting literature. We use PIV measurements across a range of canonical flows to assess the performance of the uncertainty schemes.The results show that the combined uncertainty method outperforms the individual methods, and establish the meta-uncertainty as a useful reliability assessment tool for PIV uncertainty quantification.
△ Less
Submitted 21 December, 2020;
originally announced December 2020.
-
Filamentary Surface Plasma Discharge Flow Length and Time Scales
Authors:
Lalit K. Rajendran,
Bhavini Singh,
Pavlos P. Vlachos,
Sally P. M. Bane
Abstract:
Nanosecond Surface Dielectric Barrier Discharge (ns-SDBDs) are a class of plasma actuators that utilize a high-voltage pulse of nanosecond duration between two surface-mounted electrodes to create an electrical breakdown of air, along with rapid heating. These actuators usually produce multiple filaments when operated at high pulse frequencies, and the rapid heating leads to the formation of shock…
▽ More
Nanosecond Surface Dielectric Barrier Discharge (ns-SDBDs) are a class of plasma actuators that utilize a high-voltage pulse of nanosecond duration between two surface-mounted electrodes to create an electrical breakdown of air, along with rapid heating. These actuators usually produce multiple filaments when operated at high pulse frequencies, and the rapid heating leads to the formation of shock waves and complex flow fields. In this work we replicate a single filament of the ns-SDBDs and characterize the induced flow using velocity measurements from particle image velocimetry and density measurements from background-oriented schlieren. The discharge is produced by a high voltage electrical pulse between two copper electrodes on an acrylic base. A hot gas kernel characterizes the flow field formed close to the electrodes that expands and cools over time and a vortex ring that propagates away from the surface while entraining cold ambient fluid. The gas density deficit inside the kernel displays a power-law decay over time. Based on the observations, we develop a simplified theoretical model based on vortex-driven cooling and perform a scaling analysis to obtain the induced flow length and time scales. The results show that the cooling process's time scales correspond to a circulation-based time scale of the vortex ring, and the length scale of the kernel corresponds to the vortex ring radius. These findings can guide the choice of optimal filament spacing and pulse frequencies in the design, deployment, and operation of nanosecond surface dielectric barrier discharges (ns-SDBDs) for flow control.
△ Less
Submitted 16 October, 2020;
originally announced October 2020.
-
Vortex rings drive entrainment and cooling in flow induced by a spark discharge
Authors:
Bhavini Singh,
Lalit K. Rajendran,
Jiacheng Zhang,
Pavlos P. Vlachos,
Sally P. M. Bane
Abstract:
Spark plasma discharges induce vortex rings and a hot gas kernel. We develop a model to describe the late stage of the spark induced flow and the role of the vortex rings in the entrainment of cold ambient gas and the cooling of the hot gas kernel. The model is tested in a plasma-induced flow, using density and velocity measurements obtained from simultaneous stereoscopic particle image velocimetr…
▽ More
Spark plasma discharges induce vortex rings and a hot gas kernel. We develop a model to describe the late stage of the spark induced flow and the role of the vortex rings in the entrainment of cold ambient gas and the cooling of the hot gas kernel. The model is tested in a plasma-induced flow, using density and velocity measurements obtained from simultaneous stereoscopic particle image velocimetry (S-PIV) and background oriented schlieren (BOS). We show that the spatial distribution of the hot kernel follows the motion of the vortex rings, whose radial expansion increases with the electrical energy deposited during the spark discharge. The vortex ring cooling model establishes that entrainment in the convective cooling regime is induced by the vortex rings and governs the cooling of the hot gas kernel, and the rate of cooling increases with the electrical energy deposited during the spark discharge.
△ Less
Submitted 21 July, 2020;
originally announced July 2020.
-
Weighted Least Squares (WLS) Density Integration for Background Oriented Schlieren (BOS)
Authors:
Lalit Rajendran,
Jiacheng Zhang,
Sally Bane,
Pavlos Vlachos
Abstract:
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 den…
▽ More
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.
△ Less
Submitted 2 April, 2020;
originally announced April 2020.
-
Direct estimation of global longitudinal strain from echocardiograms using a logarithm-scaled Fourier magnitude correlation
Authors:
Brett A. Meyers,
Melissa C. Brindise,
Vivek Jani,
Shelby Kutty,
Pavlos P. Vlachos
Abstract:
We present a new method for measuring global longitudinal strain and global longitudinal strain rate from 2D echocardiograms using a logarithmic-transform correlation (LTC) method. In contrast to traditional echocardiography strain analysis, our approach does not require cardiac chamber segmentation and regularization. The algorithm was benchmarked against two conventional strain analysis methods…
▽ More
We present a new method for measuring global longitudinal strain and global longitudinal strain rate from 2D echocardiograms using a logarithmic-transform correlation (LTC) method. In contrast to traditional echocardiography strain analysis, our approach does not require cardiac chamber segmentation and regularization. The algorithm was benchmarked against two conventional strain analysis methods using synthetic left ventricle ultrasound images. Measurement error was assessed as a function of contrast-to-noise ratio (CNR) using mean absolute error and root-mean-square error and showed better agreement to the ground truth for strain($R^2$ = 0.91) and strain rate($R^2$ = 0.85) as compared to conventional algorithms (strain($R^2$ = 0.7), strain rate($R^2$ = 0.7)). Also, our method was unaffected by CNR. A 200% increase in strain measurement accuracy was observed compared to the conventional algorithms. Subsequently, we tested the method using a 54-subject clinical cohort (20 subjects diseased with cardiomyopathy, 34 healthy controls). Our method distinguished between normal and abnormal left ventricular function with an AUC = 0.85, a 10% improvement over the conventional GLS algorithms.
△ Less
Submitted 24 May, 2021; v1 submitted 25 March, 2020;
originally announced March 2020.
-
To Seal or Not To Seal
Authors:
Javad Eshraghi,
Sunghwan Jung,
Pavlos P. Vlachos
Abstract:
When an object impacts the free surface of a liquid, it ejects a splash curtain upwards and creates an air cavity below the free surface. As the object descends into the liquid, the air cavity eventually closes under the action of hydrostatic pressure (deep seal). In contrast, the surface curtain may splash outwards or dome over and close, creating a surface seal. In this paper we experimentally i…
▽ More
When an object impacts the free surface of a liquid, it ejects a splash curtain upwards and creates an air cavity below the free surface. As the object descends into the liquid, the air cavity eventually closes under the action of hydrostatic pressure (deep seal). In contrast, the surface curtain may splash outwards or dome over and close, creating a surface seal. In this paper we experimentally investigate how the splash curtain dynamics are governed by the interplay of cavity pressure difference, gravity, and surface tension and how they control the occurrence, or not, of surface seal. Based on the experimental observations and measurements, we develop an analytical model to describe the trajectory and dynamics of the splash curtain. The model enables us to reveal the scaling relationship for the dimensionless surface seal time and discover the existence of a critical dimensionless number that predicts the occurrence of surface seal. This scaling indicates that the most significant parameter governing the occurrence of surface seal is the velocity of the airflow rushing into the cavity. This is in contrast to the current understanding which considers the impact velocity as the determinant parameter.
△ Less
Submitted 12 December, 2019;
originally announced December 2019.
-
Uncertainty-based pressure field reconstruction from PIV/PTV flow measurements with generalized least-squares
Authors:
Jiacheng Zhang,
Sayantan Bhattacharya,
Pavlos P. Vlachos
Abstract:
A novel uncertainty-based pressure reconstruction method is proposed to evaluate the instantaneous pressure fields from velocity fields measured using particle image velocimetry (PIV) or particle tracking velocimetry (PTV). First, the pressure gradient fields are calculated from velocity fields, while the local and instantaneous pressure gradient uncertainty is estimated from the velocity uncertai…
▽ More
A novel uncertainty-based pressure reconstruction method is proposed to evaluate the instantaneous pressure fields from velocity fields measured using particle image velocimetry (PIV) or particle tracking velocimetry (PTV). First, the pressure gradient fields are calculated from velocity fields, while the local and instantaneous pressure gradient uncertainty is estimated from the velocity uncertainty using a linear-transformation based algorithm. The pressure field is then reconstructed by solving an overdetermined linear system which involves the pressure gradients and boundary conditions. This linear system is solved with generalized least-squares (GLS) which incorporates the previously estimated variances and covariances of the pressure gradient errors as inverse weights to optimize the reconstructed pressure field. The method was validated with synthetic velocity fields of a 2D pulsatile flow and the results show significantly improved pressure accuracy with an error reduction of as much as 250% compared to the existing baseline method of solving the pressure Poisson equation (PPE). The GLS was more robust to the velocity errors and provides greater improvement with spatially correlated velocity errors. For experimental validation, the volumetric pressure fields were evaluated from a laminar pipe flow velocity field measured using 3D PTV. The GLS reduced the median absolute pressure errors by as much as 96%.
△ Less
Submitted 29 November, 2019;
originally announced December 2019.
-
Volumetric Particle Tracking Velocimetry (PTV) Uncertainty Quantification
Authors:
Sayantan Bhattacharya,
Pavlos P. Vlachos
Abstract:
We introduce the first comprehensive approach to determine the uncertainty in volumetric Particle Tracking Velocimetry (PTV) measurements. Volumetric PTV is a state-of-the-art non-invasive flow measurement technique, which measures the velocity field by recording successive snapshots of the tracer particle motion using a multi-camera set-up. The measurement chain involves reconstructing the three-…
▽ More
We introduce the first comprehensive approach to determine the uncertainty in volumetric Particle Tracking Velocimetry (PTV) measurements. Volumetric PTV is a state-of-the-art non-invasive flow measurement technique, which measures the velocity field by recording successive snapshots of the tracer particle motion using a multi-camera set-up. The measurement chain involves reconstructing the three-dimensional particle positions by a triangulation process using the calibrated camera mapping functions. The non-linear combination of the elemental error sources during the iterative self-calibration correction and particle reconstruction steps increases the complexity of the task. Here, we first estimate the uncertainty in the particle image location, which we model as a combination of the particle position estimation uncertainty and the reprojection error uncertainty. The latter is obtained by a gaussian fit to the histogram of disparity estimates within a sub-volume. Next, we determine the uncertainty in the camera calibration coefficients. As a final step the previous two uncertainties are combined using an uncertainty propagation through the volumetric reconstruction process. The uncertainty in the velocity vector is directly obtained as a function of the reconstructed particle position uncertainty. The framework is tested with synthetic vortex ring images. The results show good agreement between the predicted and the expected RMS uncertainty values. The prediction is consistent for seeding densities tested in the range of 0.01 to 0.1 particles per pixel. Finally, the methodology is also successfully validated for an experimental test case of laminar pipe flow velocity profile measurement where the predicted uncertainty is within 17% of the RMS error value.
△ Less
Submitted 27 November, 2019;
originally announced November 2019.
-
A Wavelet Approach for the Estimation of Left Ventricular Early Filling Wave Propagation Velocity from Color-M-Mode Echocardiograms
Authors:
Sreyashi Chakraborty,
Hiroyuki Iwano,
Michael E. Hall,
Pavlos P. Vlachos
Abstract:
Objective: This work evaluates a new approach for calculating the left-ventricular (LV) early filling propagation velocity (VP) from color M-Mode (CMM) echocardiograms using wavelet analysis. Unlike current approaches, the method requires no assumptions, user inputs, or heuristic conventions. Background: Current methods for measuring VP using CMM echocardiography do not account for the spatiotempo…
▽ More
Objective: This work evaluates a new approach for calculating the left-ventricular (LV) early filling propagation velocity (VP) from color M-Mode (CMM) echocardiograms using wavelet analysis. Unlike current approaches, the method requires no assumptions, user inputs, or heuristic conventions. Background: Current methods for measuring VP using CMM echocardiography do not account for the spatiotemporal variation of the filling wave propagation velocity. They are instead confined by empirical assumptions and user inputs that significantly hinder the accuracy of VP, subsequently limiting its clinical utility. Methods: We evaluated three methods for measuring LV early filling VP: conventional VP, the strength of propagation (VS), and VP determined from the most energetically important wave (Peak-VW), using 125 patients (Group A) with normal filling (n=50), impaired relaxation (n=25), pseudonormal filling (n=25) and restrictive filling (n=25), and in 69 patients (Group B) with normal (n=32), dilated (n=15), and hypertrophic ventricles (n=22). Results: Peak-VW most accurately distinguished normal left ventricles from diseased ones. Specifically, for Group A which spanned the stages of diastolic function, using receiver operator characteristics (ROC) and measuring their corresponding area under the curve (AUC), the AUC for Peak-VW was 0.92, versus 0.62 for conventional VP, 0.63 for VS and 0.58 for intraventricular pressure difference (IVPD). These correspond to a 50-70% improvement in classification ability. Similar improvements were found using Group B. Conclusion: A new determination of the LV early filling using wavelet analysis, Peak-VW, may provide a more accurate evaluation of diastolic function than the standard method of determining Vp and enable better diagnostic classification of patients suffering with diastolic dysfunction.
△ Less
Submitted 25 October, 2019;
originally announced October 2019.
-
Uncertainty amplification due to density/refractive-index gradients in volumetric PTV and BOS experiments
Authors:
Lalit K. Rajendran,
Sally P. M. Bane,
Pavlos P. Vlachos
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 containin…
▽ More
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 refractive-index gradient induced distortion parameters, provides a methodology for accounting its effect on uncertainty quantification and provides a framework for optimizing experiment design.
△ Less
Submitted 21 October, 2019;
originally announced October 2019.
-
Two Regime Cooling in Flow Induced by a Spark Discharge
Authors:
Bhavini Singh,
Lalit Rajendran,
Pavlos P Vlachos,
Sally Bane
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…
▽ More
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.
△ Less
Submitted 19 September, 2019;
originally announced September 2019.
-
Uncertainty Quantification in density estimation from Background Oriented Schlieren (BOS) measurements
Authors:
Lalit K. Rajendran,
Jiacheng Zhang,
Sayantan Bhattacharya,
Sally P . M. Bane,
Pavlos P. Vlachos
Abstract:
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 uncerta…
▽ More
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.
△ Less
Submitted 14 September, 2019;
originally announced September 2019.
-
4D-Flow MRI Pressure Estimation Using Velocity Measurement-Error based Weighted Least-Squares
Authors:
Jiacheng Zhang,
Melissa C. Brindise,
Sean Rothenberger,
Susanne Schnell,
Michael Markl,
David Saloner,
Vitaliy L. Rayz,
Pavlos P. Vlachos
Abstract:
This work introduces a 4D-flow magnetic resonance imaging (MRI) pressure reconstruction method which employs weighted least-squares (WLS) for pressure integration. Pressure gradients are calculated from the velocity fields, and velocity errors are estimated from the velocity divergence for incompressible flow. Pressure gradient errors are estimated by propagating the velocity errors through Navier…
▽ More
This work introduces a 4D-flow magnetic resonance imaging (MRI) pressure reconstruction method which employs weighted least-squares (WLS) for pressure integration. Pressure gradients are calculated from the velocity fields, and velocity errors are estimated from the velocity divergence for incompressible flow. Pressure gradient errors are estimated by propagating the velocity errors through Navier-Stokes momentum equation. A weight matrix is generated based on the pressure gradient errors, then employed for pressure reconstruction. The pressure reconstruction method was demonstrated and analyzed using synthetic velocity fields as well as Poiseuille flow measured using in vitro 4D-flow MRI. Performance of the proposed WLS method was compared to the method of solving the pressure Poisson equation which has been the primary method used in the previous studies. Error analysis indicated that the proposed method is more robust to velocity measurement errors. Improvement on pressure results was found to be more significant for the cases with spatially-varying velocity error level, with reductions in error ranging from 50% to over 200%. Finally, the method was applied to flow in a patient-specific cerebral aneurysm. Validation was performed with in vitro flow data collected using Particle Tracking Velocimetry (PTV) and Shake the Box (STB) method, and in vivo flow measurement obtained using 4D-flow MRI. Pressure calculated by WLS, as opposed to the Poisson equation, was more consistent with the flow structures and showed better agreement between the in vivo and in vitro data. These results suggest the utility of WLS method to obtain reliable pressure field from clinical flow measurement data.
△ Less
Submitted 30 April, 2019;
originally announced April 2019.
-
Flow field evolution and entrainment in a free surface plunging jet
Authors:
Syed Harris Hassan,
Tianqi Guo,
Pavlos P. Vlachos
Abstract:
We investigate ambient fluid entrainment and near-field flow characteristics of a free surface plunging jet for five Reynolds numbers ranging from 3000 to 10000 using time-resolved stereo particle image velocimetry (SPIV). We present time-averaged velocities, RMS velocity fluctuations, mean entrainment and unsteady flow features and compare them with previous studies on free jets. We find that plu…
▽ More
We investigate ambient fluid entrainment and near-field flow characteristics of a free surface plunging jet for five Reynolds numbers ranging from 3000 to 10000 using time-resolved stereo particle image velocimetry (SPIV). We present time-averaged velocities, RMS velocity fluctuations, mean entrainment and unsteady flow features and compare them with previous studies on free jets. We find that plunging jets have a smaller potential core length, and earlier decay of the mean centerline velocity. The peak RMS velocity fluctuations occur at a location significantly upstream compared to the free jets reported in the literature. Near-field ambient fluid entrainment of plunging jets is measured for the first time and is found to be considerably higher than free jets in the low Reynolds number range. For the plunging jet case at Re = 3000, faster jet decay, higher levels of turbulent intensity in the near-field, and augmented mass entrainment result from strong primary vortices that give the turbulent/non-turbulent interface (TNTI) its convoluted shape which facilitates both bulk entrapment of ambient fluid and small scale nibbling because of larger surface area. These primary vortices occur right below the free surface and disintegrate into secondary structures at axial locations that are upstream compared to those of free jets. At higher Reynolds numbers, primary vortices are smaller in size, weak in swirling strength, and disintegrate prematurely, resulting in suppressed mixing and reduced entrainment efficiency.
△ Less
Submitted 5 March, 2019;
originally announced March 2019.
-
Patient-Specific Cerebral Aneurysm Hemodynamics: Comparison of in vitro Volumetric Particle Velocimetry, Computational Fluid Dynamics (CFD), and in vivo 4D Flow MRI
Authors:
Melissa C. Brindise,
Sean Rothenberger,
Benjamin Dickerhoff,
Susanne Schnell,
Michael Markl,
David Saloner,
Vitaliy L. Rayz,
Pavlos P. Vlachos
Abstract:
Typical approaches to patient-specific hemodynamic studies of cerebral aneurysms use image based computational fluid dynamics (CFD) and seek to statistically correlate parameters such as wall shear stress (WSS) and oscillatory shear index (OSI) to risk of growth and rupture. However, such studies have reported contradictory results, emphasizing the need for in-depth comparisons of volumetric exper…
▽ More
Typical approaches to patient-specific hemodynamic studies of cerebral aneurysms use image based computational fluid dynamics (CFD) and seek to statistically correlate parameters such as wall shear stress (WSS) and oscillatory shear index (OSI) to risk of growth and rupture. However, such studies have reported contradictory results, emphasizing the need for in-depth comparisons of volumetric experiments and CFD. In this work, we conducted tomographic particle velocimetry experiments using two patient-specific cerebral aneurysm models (basilar tip and internal carotid artery) under pulsatile flow conditions and processed the particle images using Shake the Box (STB), a particle tracking method. The STB data was compared to that obtained from in vivo 4D flow MRI and CFD. Although qualitative agreement of flow pathlines across modalities was observed, each modality maintained notably unique spatiotemporal distributions of low normalized WSS regions. Analysis of time averaged WSS (TAWSS), OSI, and Relative Residence Time (RRT) demonstrated that non-dimensional parameters, such as OSI, may be more robust to the varying assumptions, limitations, and spatial resolutions of each subject and modality. These results suggest a need for further multi-modality analysis as well as development of non-dimensional hemodynamic parameters and correlation of such metrics to aneurysm risk of growth and rupture.
△ Less
Submitted 4 March, 2019;
originally announced March 2019.
-
Dot-Tracking Methodology for Background Oriented Schlieren (BOS)
Authors:
Lalit K. Rajendran,
Sally P. M. Bane,
Pavlos P. Vlachos
Abstract:
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 densiti…
▽ More
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 th exit plane of a converging-diverging nozzle.
△ Less
Submitted 27 December, 2018;
originally announced December 2018.
-
Color-Doppler Echocardiography Flow Field Velocity Reconstruction Using a Streamfunction -Vorticity Formulation
Authors:
Brett A Meyers,
Craig J Goergen,
Patrick Segers,
Pavlos P Vlachos
Abstract:
We introduce a new method, (Doppler Velocity Reconstruction or DoVeR), for reconstructing two-component velocity fields from color Doppler scans. DoVeR employs the streamfunction-vorticity equation, which satisfies mass conservation while accurately approximating the flow rate of rotation. We validated DoVeR using artificial color Doppler images generated from computational fluid dynamics models o…
▽ More
We introduce a new method, (Doppler Velocity Reconstruction or DoVeR), for reconstructing two-component velocity fields from color Doppler scans. DoVeR employs the streamfunction-vorticity equation, which satisfies mass conservation while accurately approximating the flow rate of rotation. We validated DoVeR using artificial color Doppler images generated from computational fluid dynamics models of left ventricle (LV) flow. We compare DoVeR against the conventional intraventricular Vector Flow Mapping (iVFM1D) and reformulated iVFM (iVFM2D). LV model error analysis showed DoVeR is more robust to noise and probe placement, with noise RMS errors (nRMSE) between 3.81%-6.67%, while the iVFM methods delivered 4.16%-24.17% for iVFM1D and 4.06%-400.21% for iVFM2D. We test the DoVeR and iVFM methods using in-vivo mouse-LV ultrasound scans. DoVeR yielded more hemodynamically accurate reconstructions, suggesting that it can provide a more reliable approach for robust quantification of cardiac flow.
△ Less
Submitted 21 October, 2020; v1 submitted 26 December, 2018;
originally announced December 2018.
-
PIV/BOS Synthetic Image Generation in Variable Density Environments for Error Analysis and Experiment Design
Authors:
Lalit K. Rajendran,
Sally P. M. Bane,
Pavlos P. Vlachos
Abstract:
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 generat…
▽ More
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 (DNS) 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.
△ Less
Submitted 12 December, 2018;
originally announced December 2018.
-
Particle Image Velocimetry (PIV) Uncertainty Quantification Using Moment of Correlation (MC) Plane
Authors:
Sayantan Bhattacharya,
John J. Charonko,
Pavlos P. Vlachos
Abstract:
We present a new uncertainty estimation method for Particle Image Velocimetry (PIV), that uses the correlation plane as a model for the probability density function (PDF) of displacements and calculates the second order moment of the correlation (MC). The cross-correlation between particle image patterns is the summation of all particle matches convolved with the apparent particle image diameter.…
▽ More
We present a new uncertainty estimation method for Particle Image Velocimetry (PIV), that uses the correlation plane as a model for the probability density function (PDF) of displacements and calculates the second order moment of the correlation (MC). The cross-correlation between particle image patterns is the summation of all particle matches convolved with the apparent particle image diameter. MC uses this property to estimate the PIV uncertainty from the shape of the cross-correlation plane. In this new approach, the Generalized Cross-Correlation (GCC) plane corresponding to a PIV measurement is obtained by removing the particle diameter contribution. The GCC primary peak represents a discretization of the displacement PDF, from which the standard uncertainty is obtained by convolving the GCC plane with a Gaussian function. Then a Gaussian least-squares-fit is applied to the peak region, accounting for the stretching and rotation of the peak, due to the local velocity gradients and the effect of the convolved Gaussian. The MC method was tested with simulated image sets and the predicted uncertainties show good sensitivity to the error sources and agreement with the expected RMS error. Subsequently, the method was demonstrated in three PIV challenge cases and two experimental datasets and was compared with the published image matching (IM) and correlation statistics (CS) techniques. Results show that the MC method has a better response to spatial variation in RMS error and the predicted uncertainty is in good agreement with the expected standard uncertainty. The uncertainty prediction was also explored as a function PIV interrogation window size, and the MC method outperforms the other uncertainty methods.
△ Less
Submitted 28 April, 2018;
originally announced April 2018.
-
Confocal laser scanning microscopy image correlation for nanoparticle flow velocimetry
Authors:
Brian Jun,
Matthew Giarra,
Haisheng Yang,
Russell Main,
Pavlos Vlachos
Abstract:
We present a new particle image correlation technique for resolving nanoparticle flow velocity using confocal laser scanning microscopy (CLSM). The two primary issues that complicate nanoparticle scanning laser image correlation (SLIC) based velocimetry are (1) the use of diffusion dominated nanoparticles as flow tracers, which introduce a random decorrelating error into the velocity estimate, and…
▽ More
We present a new particle image correlation technique for resolving nanoparticle flow velocity using confocal laser scanning microscopy (CLSM). The two primary issues that complicate nanoparticle scanning laser image correlation (SLIC) based velocimetry are (1) the use of diffusion dominated nanoparticles as flow tracers, which introduce a random decorrelating error into the velocity estimate, and (2) the effects of the scanning laser image acquisition, which introduces a bias error. To date, no study has quantified these errors or demonstrated a means to deal with them in SLIC velocimetry. In this work, we build upon the robust phase correlation (RPC) and existing methods of SLIC to quantify and mitigate these errors. First, we implement an ensemble RPC instead of using an ensemble standard cross correlation, and develop an SLIC optimal filter that maximizes the correlation strength in order to reliably and accurately detect the correlation peak representing the most probable average displacement of the nanoparticles. Secondly, we developed an analytical model of the SLIC measurement bias error due to image scanning of diffusion dominated tracer particles. We show that the bias error depends only on the ratio of the mean velocity of the tracer particles to that of the laser scanner and we use this model to correct the induced errors. We validated our technique using synthetic images and experimentally obtained SLIC images of nanoparticle flow through a micro-channel. Our technique reduced the error by up to a factor of ten compared to other SLIC algorithms for the images tested in this study. Moreover, our optimized RPC filter is reducing the number of image pairs required for the convergence of the ensemble correlation by two orders of magnitude compared to the standard cross correlation.
△ Less
Submitted 10 March, 2016;
originally announced March 2016.
-
Hemodynamics of Stent Implantation Procedures in Coronary Bifurcations: an in vitro study
Authors:
Melissa C. Brindise,
Claudio Chiastra,
Francesco Burzotta,
Francesco Migliavacca,
Pavlos P. Vlachos
Abstract:
Stent implantation in coronary bifurcations presents unique challenges and currently there is no universally accepted stent deployment approach. Despite clinical and computational studies, to date, the effect of each stent implantation method on the coronary artery hemodynamics is not well understood. In this study the hemodynamics of stented coronary bifurcations under pulsatile flow conditions w…
▽ More
Stent implantation in coronary bifurcations presents unique challenges and currently there is no universally accepted stent deployment approach. Despite clinical and computational studies, to date, the effect of each stent implantation method on the coronary artery hemodynamics is not well understood. In this study the hemodynamics of stented coronary bifurcations under pulsatile flow conditions were investigated experimentally. Three implantation methods, provisional side branch (PSB), culotte (CUL), and crush (CRU), were investigated using time-resolved particle image velocimetry (PIV) to measure the velocity fields. Subsequently, hemodynamic parameters including wall shear stress (WSS), oscillatory shear index (OSI), and relative residence time (RRT) were calculated and the pressure field through the vessel was non-invasively quantified. The effects of each stented case were evaluated and compared against an un-stented case. CRU provided the lowest compliance mismatch, but demonstrated detrimental stent interactions. PSB, the clinically preferred method, and CUL maintained many normal flow conditions. However, PSB provided about a 300% increase in both OSI and RRT. CUL yielded a 10% and 85% increase in OSI and RRT, respectively. The results of this study support the concept that different bifurcation stenting techniques result in hemodynamic environments that deviate from that of un-stented bifurcations, to varying degrees.
△ Less
Submitted 10 March, 2016;
originally announced March 2016.
-
Time-resolved particle image velocimetry measurements with wall shear stress and uncertainty quantification for the FDA benchmark nozzle model
Authors:
Jaime S. Raben,
Prasanna Hariharan,
Ronald Robinson,
Richard Malinauskas,
Pavlos P. Vlachos
Abstract:
We present validation of benchmark experimental data for computational fluid dynamics (CFD) analyses of medical devices using advanced Particle Image Velocimetry (PIV) processing and post-processing techniques. This work is an extension of a previous FDA-sponsored multi-laboratory study, which used a medical device mimicking geometry referred to as the FDA benchmark nozzle model. Time-resolved PIV…
▽ More
We present validation of benchmark experimental data for computational fluid dynamics (CFD) analyses of medical devices using advanced Particle Image Velocimetry (PIV) processing and post-processing techniques. This work is an extension of a previous FDA-sponsored multi-laboratory study, which used a medical device mimicking geometry referred to as the FDA benchmark nozzle model. Time-resolved PIV analysis was performed in five overlapping regions of the model for Reynolds numbers in the nozzle throat of 500, 2,000, 5,000, and 8,000. Images included a two-fold increase in spatial resolution in comparison to the previous study. Data was processed using ensemble correlation, dynamic range enhancement, and phase correlations to increase signal-to-noise ratios and measurement accuracy, and to resolve flow regions with large velocity ranges and gradients, which is typical of many blood-contacting medical devices. Parameters relevant to device safety, including shear stress at the wall and in bulk flow, were computed using radial basis functions (RBF) to improve accuracy. In-field spatially resolved pressure distributions, Reynolds stresses and energy dissipation rates were computed from PIV measurements. Velocity measurement uncertainty was estimated directly from the PIV correlation plane, and uncertainty analysis for wall shear stress at each measurement location was performed using a Monte Carlo model. (truncated abstract)
△ Less
Submitted 13 May, 2014;
originally announced May 2014.
-
Particle image velocimetry correlation signal-to-noise ratio metrics and measurement uncertainty quantification
Authors:
Zhenyu Xue,
John J. Charonko,
Pavlos P. Vlachos
Abstract:
In particle image velocimetry (PIV) the measurement signal is contained in the recorded intensity of the particle image pattern superimposed on a variety of noise sources. The signal-to-noise-ratio (SNR) strength governs the resulting PIV cross correlation and ultimately the accuracy and uncertainty of the resulting PIV measurement. Hence we posit that correlation SNR metrics calculated from the c…
▽ More
In particle image velocimetry (PIV) the measurement signal is contained in the recorded intensity of the particle image pattern superimposed on a variety of noise sources. The signal-to-noise-ratio (SNR) strength governs the resulting PIV cross correlation and ultimately the accuracy and uncertainty of the resulting PIV measurement. Hence we posit that correlation SNR metrics calculated from the correlation plane can be used to quantify the quality of the correlation and the resulting uncertainty of an individual measurement. In this paper we present a framework for evaluating the correlation SNR using a set of different metrics, which in turn are used to develop models for uncertainty estimation. The SNR metrics and corresponding models presented herein are expanded to be applicable to both standard and filtered correlations. In addition, the notion of a valid measurement is redefined with respect to the correlation peak width in order to be consistent with uncertainty quantification principles and distinct from an outlier measurement. Finally the type and significance of the error distribution function is investigated. These advancements lead to robust uncertainty estimation models, which are tested against both synthetic benchmark data as well as actual experimental measurements. For all cases considered here, standard uncertainties are estimated at the 68.5% confidence level while expanded uncertainties are estimated at 95% confidence level. For all cases the resulting calculated coverage factors approximate the expected theoretical values thus demonstrating the applicability of these new models for estimation of uncertainty for individual PIV measurements.
△ Less
Submitted 12 May, 2014;
originally announced May 2014.
-
Demonstration of Experimental Three Dimensional Finite Time Lyapunov Exponents with Inertial Particles
Authors:
Samuel G. Raben,
Shane D. Ross,
Pavlos P. Vlachos
Abstract:
This work provides an experimental method for simultaneously measuring finite time Lyapunov exponent fields for multiple particle groups, including non-flow tracers, in three-dimensional multiphase flows. From sequences of particle images, e.g., from experimental fluid imaging techniques, we can directly compute the flow map and coherent structures, with out performing the computationally costly n…
▽ More
This work provides an experimental method for simultaneously measuring finite time Lyapunov exponent fields for multiple particle groups, including non-flow tracers, in three-dimensional multiphase flows. From sequences of particle images, e.g., from experimental fluid imaging techniques, we can directly compute the flow map and coherent structures, with out performing the computationally costly numerical integration. This is particularly useful to find three-dimensional Lagrangian coherent structures for inertial particles, that do not follow the bulk fluid velocity, as we demonstrate for a grid turbulence experiment. The technique described may provide a new means for exploring the physics of experimental multi-phase flows.
△ Less
Submitted 4 October, 2013; v1 submitted 12 September, 2013;
originally announced September 2013.
-
Lift and wakes of flying snakes
Authors:
Anush Krishnan,
John J. Socha,
Pavlos P. Vlachos,
L. A. Barba
Abstract:
Flying snakes use a unique method of aerial locomotion: they jump from tree branches, flatten their bodies and undulate through the air to produce a glide. The shape of their body cross-section during the glide plays an important role in generating lift. This paper presents a computational investigation of the aerodynamics of the cross-sectional shape. Two-dimensional simulations of incompressible…
▽ More
Flying snakes use a unique method of aerial locomotion: they jump from tree branches, flatten their bodies and undulate through the air to produce a glide. The shape of their body cross-section during the glide plays an important role in generating lift. This paper presents a computational investigation of the aerodynamics of the cross-sectional shape. Two-dimensional simulations of incompressible flow past the anatomically correct cross-section of the species Chrysopelea paradisi show that a significant enhancement in lift appears at a 35-degrees angle of attack, above Reynolds numbers 2000. Previous experiments on physical models also obtained an increased lift, at the same angle of attack. The flow is inherently three-dimensional in physical experiments, due to fluid instabilities, and it is thus intriguing that the enhanced lift also appears in the two-dimensional simulations. The simulations point to the lift enhancement arising from the early separation of the boundary layer on the dorsal surface of the snake profile, without stall. The separated shear layer rolls up and interacts with secondary vorticity in the near-wake, inducing the primary vortex to remain closer to the body and thus cause enhanced suction, resulting in higher lift.
△ Less
Submitted 7 February, 2014; v1 submitted 11 September, 2013;
originally announced September 2013.
-
Straight-line and turning locomotion of Paramecia
Authors:
Saikat Jana,
Matthew Giarra,
Pavlos Vlachos,
Sunghwan Jung
Abstract:
In this fluid dynamics video we investigate the flow field around straight-line swimming and right and left turning Paramecia using micro- particle image velocimetry (μPIV). A Paramecium controls its ciliary beating to produce different fluid velocities on either side of its body. This phenomenon is visualized by applying μPIV to images in which Paramecia swim in a dilute suspension of 1 μm polyst…
▽ More
In this fluid dynamics video we investigate the flow field around straight-line swimming and right and left turning Paramecia using micro- particle image velocimetry (μPIV). A Paramecium controls its ciliary beating to produce different fluid velocities on either side of its body. This phenomenon is visualized by applying μPIV to images in which Paramecia swim in a dilute suspension of 1 μm polystyrene spheres. Paramecia that swim straight exhibit similar magnitudes of velocity on either side of their bodies. In contrast, right-turning Paramecia exhibit greater magnitudes of velocity on their right sides, while left-turning organisms show the opposite.
△ Less
Submitted 12 October, 2012;
originally announced October 2012.