-
Effects of settling on inertial particle slip velocity statistics in wall bounded flows
Authors:
Andrew P. Grace,
David Richter,
Tim Berk,
Andrew D. Bragg
Abstract:
Developing reduced order models for the transport of solid particles in turbulence typically requires a statistical description of the particle-turbulence interactions. In this work, we utilize a statistical framework to derive continuum equations for the moments of the slip velocity of inertial settling Lagrangian particles in a turbulent boundary layer. Using coupled Eulerian-Lagrangian direct n…
▽ More
Developing reduced order models for the transport of solid particles in turbulence typically requires a statistical description of the particle-turbulence interactions. In this work, we utilize a statistical framework to derive continuum equations for the moments of the slip velocity of inertial settling Lagrangian particles in a turbulent boundary layer. Using coupled Eulerian-Lagrangian direct numerical simulations, we then identify the dominant mechanisms controlling the slip velocity variance, and find that for a range of St+, Sv+, and Re, the slip variance is primarily controlled by local differences between the "seen" variance and the particle velocity variance, while terms appearing due to the inhomogeneity of the turbulence are sub-leading until Sv+ becomes large. We also consider several comparative metrics to assess the relative magnitudes of the fluctuating slip velocity and the mean slip velocity, and we find that the vertical mean slip increases rapidly with Sv+, rendering the variance relatively small -- an effect found to be most substantial for Sv+>1. Finally, we compare the results to a model of the acceleration variance Berk and Coletti (2021) based the concept of a response function described in Csanady (1963), highlighting the role of the crossing trajectories mechanism. We find that while there is good agreement for low Sv+, systematic errors remain, possibly due to implicit non-local effects arising from rapid particle settling and inhomogeneous turbulence. We conclude with a discussion of the implications of this work for modeling the transport of coarse dust grains in the atmospheric surface layer.
△ Less
Submitted 19 June, 2024;
originally announced June 2024.
-
Asymptotic analysis of mixing in stratified turbulent flows, and the conditions for an inertial sub-range
Authors:
Andrew D. Bragg,
Stephen M. de Bruyn Kops
Abstract:
In an important study, Maffioli et al. (J. Fluid Mech., Vol. 794 , 2016) used a scaling analysis to predict that in the weakly stratified flow regime $Fr_h\gg1$ ($Fr_h$ is the horizontal Froude number), the mixing coefficient $Γ$ (defined as the ratio of the dissipation rates of potential to kinetic energy) scales as $Γ\sim O(Fr_h^{-2})$. Direct numerical simulations confirmed this result, and als…
▽ More
In an important study, Maffioli et al. (J. Fluid Mech., Vol. 794 , 2016) used a scaling analysis to predict that in the weakly stratified flow regime $Fr_h\gg1$ ($Fr_h$ is the horizontal Froude number), the mixing coefficient $Γ$ (defined as the ratio of the dissipation rates of potential to kinetic energy) scales as $Γ\sim O(Fr_h^{-2})$. Direct numerical simulations confirmed this result, and also indicated that for the strongly stratified regime $Fr_h\ll 1$, $Γ\sim O(1)$. Furthermore, the study argued that $Γ$ does not depend on the buoyancy Reynolds number $Re_b$, but only on $Fr_h$. We present an asymptotic analysis to predict theoretically how $Γ$ should behave for $Fr_h\ll1$ and $Fr_h\gg1$ in the limit $Re_b\to\infty$. To correctly handle the singular limit $Re_b\to\infty$ we perform the asymptotic analysis on the filtered Boussinesq-Navier-Stokes equations, and demonstrate the precise sense in which the inviscid scaling analysis of Billant \& Chomaz (Phys. Fluids, vol. 13, 1645-1651, 2001) applies to viscous flows with $Re_b\to\infty$. The analysis yields $Γ\sim O(Fr_h^{-2}(1+Fr_h^{-2}))$ for $Fr_h\gg1$ and $Γ\sim O(1+Fr_h^{2})$ for $Fr_h\ll 1$, providing a theoretical basis for the numerical observation made by Maffioli et al, as well as predicting the sub-leading behavior. Our analysis also shows that the Ozmidov scale $L_O$ does not describe the scale below which buoyancy forces are sub-leading, which is instead given by $O(Fr_h^{1/2} L_O)$, and that the condition for there to be an inertial sub-range when $Fr_h\ll 1$ is not $Re_b\gg1$, but the more restrictive condition $Re_b\gg Fr_h^{-4/3}$.
△ Less
Submitted 16 February, 2024;
originally announced February 2024.
-
Investigating the parametric dependence of the impact of two-way coupling on inertial particle settling in turbulence
Authors:
Soumak Bhattacharjee,
Josin Tom,
Maurizio Carbone,
Andrew D. Bragg
Abstract:
Tom et al.\ (J.\ Fluid Mech.\ 947, A7, 2022) investigated the impact of two-way coupling (2WC) on particle settling in turbulence. For the limited parameter choices explored, it was found that 2WC substantially enhances particle settling compared to the one-way coupled (1WC) case, even at low mass loading $Φ_m$. Moreover, contrary to previous claims, it was demonstrated that preferential sweeping…
▽ More
Tom et al.\ (J.\ Fluid Mech.\ 947, A7, 2022) investigated the impact of two-way coupling (2WC) on particle settling in turbulence. For the limited parameter choices explored, it was found that 2WC substantially enhances particle settling compared to the one-way coupled (1WC) case, even at low mass loading $Φ_m$. Moreover, contrary to previous claims, it was demonstrated that preferential sweeping remains the mechanism responsible for the particles settling faster than the Stokes settling velocity in 2WC flows. However, crucial questions remain: 1) how small must $Φ_m$ be for the effects of 2WC on particle settling to be negligible? 2) does the preferential sweeping mechanism remain relevant in 2WC flows as $Φ_m$ is increased? To answer these, we explore a much broader portion of the parameter space, and our simulations cover cases where the impact of 2WC on the global fluid statistics ranges from negligible to strong. We find that even for $Φ_m=7.5\times 10^{-3}$, 2WC can noticeably increase the settling for some choices of the Stokes and Froude numbers. We also demonstrate that even when $Φ_m$ is large enough for the global fluid statistics to be strongly affected by the particles, preferential sweeping is still the mechanism responsible for the enhanced particle settling. The difference between the 1WC and 2WC cases is that, in the latter the particles are not merely swept around the downward-moving side of vortices, but they also drag the fluid with them as they move down.
△ Less
Submitted 14 September, 2023;
originally announced September 2023.
-
Understanding the effect of Prandtl number on momentum and scalar mixing rates in neutral and stably stratified flows using gradient field dynamics
Authors:
Andrew D. Bragg,
Stephen M. de Bruyn Kops
Abstract:
Recently, direct numerical simulations (DNS) of stably stratified turbulence have shown that as the Prandtl number ($Pr$) is increased from 1 to 7, the mean turbulent potential energy dissipation rate (TPE-DR) drops dramatically, while the mean turbulent kinetic energy dissipation rate (TKE-DR) increases significantly. Through an analysis of the equations governing the fluctuating velocity and den…
▽ More
Recently, direct numerical simulations (DNS) of stably stratified turbulence have shown that as the Prandtl number ($Pr$) is increased from 1 to 7, the mean turbulent potential energy dissipation rate (TPE-DR) drops dramatically, while the mean turbulent kinetic energy dissipation rate (TKE-DR) increases significantly. Through an analysis of the equations governing the fluctuating velocity and density gradients we provide a mechanistic explanation for this surprising behavior and test the predictions using DNS. We show that the mean density gradient gives rise to a mechanism that opposes the production of fluctuating density gradients, and this is connected to the emergence of ramp-cliffs. The same term appears in the velocity gradient equation but with the opposite sign, and is the contribution from buoyancy. This term is ultimately the reason why the TPE-DR reduces while the TKE-DR increases with increasing $Pr$. Our analysis also predicts that the effects of buoyancy on the smallest scales of the flow become stronger as $Pr$ is increased, and this is confirmed by our DNS data. A consequence of this is that the standard buoyancy Reynolds number does not correctly estimate the impact of buoyancy at the smallest scales when $Pr$ deviates from 1, and we derive a suitable alternative parameter. Finally, an analysis of the filtered gradient equations reveals that the mean density gradient term changes sign at sufficiently large scales, such that buoyancy acts as a source for velocity gradients at small scales, but as a sink at large scales.
△ Less
Submitted 1 August, 2023;
originally announced August 2023.
-
Comparing phase-space and phenomenological modeling approaches for Lagrangian particles settling in a turbulent boundary layer
Authors:
Andrew P. Grace,
David H. Richter,
Andrew D. Bragg
Abstract:
Under the right circumstances, inertial particles (such as sand or dust) settling through the atmospheric boundary layer can experience a net enhancement in their average settling velocity due to their inertia. Since this enhancement arises due to their interactions with the surrounding turbulence it must be modelled at coarse scales. Models for the enhanced settling velocity (or deposition) of th…
▽ More
Under the right circumstances, inertial particles (such as sand or dust) settling through the atmospheric boundary layer can experience a net enhancement in their average settling velocity due to their inertia. Since this enhancement arises due to their interactions with the surrounding turbulence it must be modelled at coarse scales. Models for the enhanced settling velocity (or deposition) of the dispersed phase that find practical use in mesoscale weather models are often ad hoc or are built on phenomenological closure assumptions, meaning that the general deposition rate of particle is a key uncertainty. Instead of taking a phenomenological approach, exact phase space methods can be used to model the physical mechanisms responsible for the enhanced settling, and a more general parameterization of the enhanced settling of inertial particles can be built. In this work, we use direct numerical simulations (DNS) and phase space methods to evaluate the efficacy of phenomenological modelling approaches for the enhanced settling velocity of inertial particles with varying friction Stokes numbers and settling velocity parameters. We use the DNS data to estimate profiles of a drift-diffusion based parameterization of the fluid velocity sampled by the particles, which is key for determining the settling velocity behaviour of particles with low to moderate Stokes number. We find that by increasing the settling velocity parameter at moderate friction Stokes number, the magnitude of preferential sweeping is modified, and this behaviour is explained by the drift component. We then use these profiles to argue that the eddy-diffusivity-like closure used in phenomenological models is incomplete, relying on inadequate empirical corrections. Finally, we discuss opportunities for reconciling exact phase space approaches with simpler phenomenological approaches for use in coarse-scale weather models.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
Fate of bubble clusters rising in a quiescent liquid
Authors:
Tian Ma,
Hendrik Hessenkemper,
Dirk Luca,
Andrew D. Bragg
Abstract:
We use experiments to study the evolution of bubble clusters in a swarm of freely rising, deformable bubbles. A new machine learning-aided algorithm allows us to identify and track bubbles in clusters and measure the cluster lifetimes. The results indicate that contamination in the carrier liquid can enhance the formation of bubble clusters and prolong the cluster lifetimes. The mean bubble rise v…
▽ More
We use experiments to study the evolution of bubble clusters in a swarm of freely rising, deformable bubbles. A new machine learning-aided algorithm allows us to identify and track bubbles in clusters and measure the cluster lifetimes. The results indicate that contamination in the carrier liquid can enhance the formation of bubble clusters and prolong the cluster lifetimes. The mean bubble rise velocities conditioned on the bubble cluster size are also explored, and we find a positive correlation between the cluster size and the rise speed of the bubbles in the cluster, with clustered bubbles rising up to $20\%$ faster than unclustered bubbles.
△ Less
Submitted 7 June, 2023; v1 submitted 3 June, 2023;
originally announced June 2023.
-
Logarithmic scaling of higher-order temperature moments in the atmospheric surface layer
Authors:
Kelly Y. Huang,
Matt K. Fu,
Clayton P. Byers,
Andrew D. Bragg,
Gabriel G. Katul
Abstract:
A generalized logarithmic law for high-order moments of passive scalars is proposed for turbulent boundary layers. This law is analogous to the generalized log law that has been proposed for high-order moments of the turbulent longitudinal velocity and is derived by combining the random sweeping decorrelation hypothesis with a spectral model informed by the attached eddy hypothesis. The proposed t…
▽ More
A generalized logarithmic law for high-order moments of passive scalars is proposed for turbulent boundary layers. This law is analogous to the generalized log law that has been proposed for high-order moments of the turbulent longitudinal velocity and is derived by combining the random sweeping decorrelation hypothesis with a spectral model informed by the attached eddy hypothesis. The proposed theory predicts that the high-order moments of passive scalar fluctuations within the inertial sublayer will vary logarithmically with wall-normal distance ($z$). The proposed theory is evaluated using high frequency time-series measurements of temperature and streamwise velocity fluctuations obtained in the first meter of the atmospheric surface layer (ASL) under near-neutral thermal stratification. The logarithmic dependence with $z$ within the inertial sublayer is observed in both the air temperature and velocity moments, with good agreement to the predictions from the proposed theory. Surprisingly, the proposed theory appears to be as, if not more, valid for transported passive scalars than for the longitudinal velocity.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Effects of surfactants on bubble-induced turbulence
Authors:
Tian Ma,
Hendrik Hessenkemper,
Dirk Lucas,
Andrew D. Bragg
Abstract:
We use experiments to explore the effect of surfactants on bubble-induced turbulence (BIT) at different scales, considering how the bubbles affect the flow kinetic energy, anisotropy and extreme events. To this end, high-resolution Particle Shadow Velocimetry measurements are carried out in a bubble column in which the flow is generated by bubble swarms rising in water for two different bubble dia…
▽ More
We use experiments to explore the effect of surfactants on bubble-induced turbulence (BIT) at different scales, considering how the bubbles affect the flow kinetic energy, anisotropy and extreme events. To this end, high-resolution Particle Shadow Velocimetry measurements are carried out in a bubble column in which the flow is generated by bubble swarms rising in water for two different bubble diameters ($3$ mm $\&$ $4$ mm) and moderate gas volume fractions ($0.5\%\sim1.3\%$). To contaminate the flow, different amounts of 1-Pentanol were added to the flow, leading to different bubble shapes and surface boundary conditions. The results reveal that with increasing surfactant concentration, the BIT generated increases in strength, even though bubbles of a given size rise more slowly with surfactants. We also find that the level of anisotropy in the flow is enhanced with increasing surfactant concentration for bubbles of the same size, and that for the same surfactant concentration, smaller bubbles generate stronger anisotropy in the flow. Concerning the intermittency quantified by the normalized probability density functions of the fluid velocity increments, our results indicate that extreme values in the velocity increments become more probable with decreasing surfactant concentration for cases with smaller bubbles and low gas void fraction, while the effect of the surfactant is much weaker for cases with larger bubble and higher void fractions.
△ Less
Submitted 7 December, 2022;
originally announced December 2022.
-
Asymptotic closure model for inertial particle transport in turbulent boundary layers
Authors:
Yan Zhang,
Guiquan Wang,
Andrew D Bragg
Abstract:
Transport equations for heavy inertial particles in turbulent boundary layers may be derived from an underlying phase-space probability density function (PDF) equation. These equations, however, are unclosed, and the standard closure approach is to use a quasi-Normal approximation (QNA) in which the fourth moments are approximated as behaving as if the velocities were Normally distributed. Except…
▽ More
Transport equations for heavy inertial particles in turbulent boundary layers may be derived from an underlying phase-space probability density function (PDF) equation. These equations, however, are unclosed, and the standard closure approach is to use a quasi-Normal approximation (QNA) in which the fourth moments are approximated as behaving as if the velocities were Normally distributed. Except for particles with weak inertia, the QNA leads to large quantitative errors, and is not consistent with the known asymptotic predictions of Sikovsky (Flow Turbulence Combust, vol. 92, 2014, pp. 41-64) for the moments of the PDF in the viscous sublayer. We derive a new closure approximation based on an asymptotic solution to the transport equations in regions where the effect of particle inertia is significant. The new closure is consistent with the asymptotic predictions of Sikovsky, but applies even outside the viscous sublayer. Comparisons with direct numerical simulations (DNS) show that the new closure gives similar results to the QNA (with the QNA results in slightly better agreement with the DNS) when the viscous Stokes number is $St<10$, but for $St>10$ the new model is in far better agreement with the DNS than the QNA. While the predictions from the new model leave room for improvement, the results suggest that this new closure strategy is a very effective alternative to the traditional QNA approach, and the closure could be refined in future work.
△ Less
Submitted 9 November, 2022;
originally announced November 2022.
-
Lagrangian model for passive scalar gradients in turbulence
Authors:
Xiaolong Zhang,
Maurizio Carbone,
Andrew D. Bragg
Abstract:
The equation for the fluid velocity gradient along a Lagrangian trajectory immediately follows from the Navier-Stokes equation. However, such an equation involves two terms that cannot be determined from the velocity gradient along the chosen Lagrangian path: the pressure Hessian and the viscous Laplacian. A recent model handles these unclosed terms using a multi-level version of the recent deform…
▽ More
The equation for the fluid velocity gradient along a Lagrangian trajectory immediately follows from the Navier-Stokes equation. However, such an equation involves two terms that cannot be determined from the velocity gradient along the chosen Lagrangian path: the pressure Hessian and the viscous Laplacian. A recent model handles these unclosed terms using a multi-level version of the recent deformation of Gaussian fields (RDGF) closure (Johnson \& Meneveau, Phys.~Rev.~Fluids, 2017). This model is in remarkable agreement with DNS data and works for arbitrary Taylor Reynolds numbers $\Rey_λ$. Inspired by this, we develop a Lagrangian model for passive scalar gradients in isotropic turbulence. The equation for passive scalar gradients also involves an unclosed term in the Lagrangian frame, namely the scalar gradient diffusion term, which we model using the RDGF approach. However, comparisons of the statistics obtained from this model with direct numerical simulation (DNS) data reveal substantial errors due to erroneously large fluctuations generated by the model. We address this defect by incorporating into the closure approximation information regarding the scalar gradient production along the local trajectory history of the particle. This modified model makes predictions for the scalar gradients, their production rates, and alignments with the strain-rate eigenvectors that are in very good agreement with DNS data. However, while the model yields valid predictions up to around $\Rey_λ\approx 500$, beyond this, the model breaks down.
△ Less
Submitted 31 October, 2022;
originally announced October 2022.
-
A stochastic model for the residence time of solid particles in turbulent Rayleigh-Bénard Flow
Authors:
Colin J. Denzel,
Andrew D. Bragg,
David H. Richter
Abstract:
The Pi Chamber, located at Michigan Technological University, generates moist turbulent Rayleigh-Bénard flow in order to replicate steady-state cloud conditions. We take inspiration from this setup and consider a particle-laden, convectively-driven turbulent flow using direct numerical simulation (DNS). The aim of our study is to develop a simple stochastic model that can accurately describe the r…
▽ More
The Pi Chamber, located at Michigan Technological University, generates moist turbulent Rayleigh-Bénard flow in order to replicate steady-state cloud conditions. We take inspiration from this setup and consider a particle-laden, convectively-driven turbulent flow using direct numerical simulation (DNS). The aim of our study is to develop a simple stochastic model that can accurately describe the residence times of the particles in the flow, this time being determined by the complex competition between the gravitational settling of the particles, and the interaction of the particles with the turbulent structures in the flow. A simple conceptual picture underlies the stochastic model, namely that the particles take repeated trips between the top and bottom boundaries, driven by the convective cells that occur in Rayleigh-Bénard turbulence, and that their residence times are determined by the time it takes to complete one of these trips, which varies from one trip to another, and the probability of falling out to the bottom boundary after each trip. Despite the simplicity of the model, it yields quantitatively accurate predictions of the distribution of the particle residence times in the flow. We independently vary the Stokes numbers and settling velocities in order to shed light on the independent roles that gravity and inertia play in governing these residence times.
△ Less
Submitted 27 August, 2022;
originally announced August 2022.
-
The effect of tilt on turbulent thermal convection for a heated soap bubble
Authors:
Xiao-Qiu He,
Yong-Liang Xiong,
Andrew D. Bragg,
Patrick Fischer,
Hamid Kellay
Abstract:
We use direct numerical simulation (DNS) to explore the effect of tilt on two-dimensional turbulent thermal convection on a half-soap bubble that is heated at its equator.In the DNS, the bubble is tilted by an angle $δ\in[0^{\circ},90^{\circ}]$, the Rayleigh number is varied between $Ra\in[3\times10^6, 3\times10^9]$, and the Prandlt number is fixed at $Pr=7$.The DNS reveals two qualitatively diffe…
▽ More
We use direct numerical simulation (DNS) to explore the effect of tilt on two-dimensional turbulent thermal convection on a half-soap bubble that is heated at its equator.In the DNS, the bubble is tilted by an angle $δ\in[0^{\circ},90^{\circ}]$, the Rayleigh number is varied between $Ra\in[3\times10^6, 3\times10^9]$, and the Prandlt number is fixed at $Pr=7$.The DNS reveals two qualitatively different flow regimes: the dynamic plume regime (DPR) and the stable plume regime (SPR).In the DPR, small dynamic plumes constantly emerge from random locations on the equator and dissipate on the bubble.In the SPR, the flow is dominated by a single large and stable plume rising from the lower edge of the bubble.The scaling behaviour of the Nusselt number $Nu$ and Reynolds number $Re$ are different in these two regimes,with $Nu\propto Ra^{0.3}$ for the DPR and $Nu\propto Ra^{0.24}$ for the SPR. Concerning $Re$, the scaling in the DPR lies between $Re\propto Ra^{0.48}$ and $Re\propto Ra^{0.53}$ depending on $Ra$ and $δ$,while in the SPR, the scaling lies between $Re\propto Ra^{0.44}$ and $Re\propto Ra^{0.45}$ depending on $δ$.The turbulent thermal and kinetic energy dissipation rates ($ε_{T^{\prime}}$ and $ε_{u^{\prime}}$, respectively) are also very different in the DPR and SPR.The probability density functions (PDF) of the normalized $\logε_{T^{\prime}}$ and $\logε_{u^{\prime}}$ are close to a Gaussian PDF for small fluctuations, but deviate considerably from a Gaussian at large fluctuations in the DPR.
△ Less
Submitted 4 August, 2022;
originally announced August 2022.
-
How two-way coupling modifies the multiscale preferential sweeping mechanism
Authors:
Josin Tom,
Maurizio Carbone,
Andrew D. Bragg
Abstract:
For one-way coupled (1WC) flows, Tom & Bragg (J. Fluid Mech., 871, pp. 244-270, 2019) advanced the analysis of Maxey (J. Fluid Mech., 174, pp. 441-465, 1987), which applied to weakly inertial particles, to particles of arbitrary inertia, and the new theoretical result revealed the role that different scales play in the preferential sweeping mechanism that leads to enhanced particle settling in tur…
▽ More
For one-way coupled (1WC) flows, Tom & Bragg (J. Fluid Mech., 871, pp. 244-270, 2019) advanced the analysis of Maxey (J. Fluid Mech., 174, pp. 441-465, 1987), which applied to weakly inertial particles, to particles of arbitrary inertia, and the new theoretical result revealed the role that different scales play in the preferential sweeping mechanism that leads to enhanced particle settling in turbulent flows. Monchaux & Dejoan (Phys. Rev. Fluids, 2, 104302, 2017) showed using direct numerical simulations (DNS) that while for low particle loading the effect of two-way coupling (2WC) on the global flow statistics is weak, 2WC enables the particles to drag the fluid in their vicinity down with them, significantly enhancing their settling, and they argued that 2WC suppresses the preferential sweeping mechanism. We explore this further by considering the impact of 2WC on the contribution made by eddies of different sizes on the particle settling. In agreement with Monchaux & Dejoan, we show that even for low loading, 2WC strongly enhances particle settling, and we show how 2WC modifies the contribution from different flow scales. However, contrary to their study, we show that preferential sweeping remains important in 2WC flows. In particular, for both 1WC and 2WC flows, the settling enhancement due to turbulence is dominated by contributions from particles in straining regions of the flow, but for the 2WC case, the particles in these regions also drag the fluid down with them, leading to an enhancement of their settling compared to the 1WC case.
△ Less
Submitted 2 February, 2022;
originally announced February 2022.
-
A Physics-Informed Vector Quantized Autoencoder for Data Compression of Turbulent Flow
Authors:
Mohammadreza Momenifar,
Enmao Diao,
Vahid Tarokh,
Andrew D. Bragg
Abstract:
Analyzing large-scale data from simulations of turbulent flows is memory intensive, requiring significant resources. This major challenge highlights the need for data compression techniques. In this study, we apply a physics-informed Deep Learning technique based on vector quantization to generate a discrete, low-dimensional representation of data from simulations of three-dimensional turbulent fl…
▽ More
Analyzing large-scale data from simulations of turbulent flows is memory intensive, requiring significant resources. This major challenge highlights the need for data compression techniques. In this study, we apply a physics-informed Deep Learning technique based on vector quantization to generate a discrete, low-dimensional representation of data from simulations of three-dimensional turbulent flows. The deep learning framework is composed of convolutional layers and incorporates physical constraints on the flow, such as preserving incompressibility and global statistical characteristics of the velocity gradients. The accuracy of the model is assessed using statistical, comparison-based similarity and physics-based metrics. The training data set is produced from Direct Numerical Simulation of an incompressible, statistically stationary, isotropic turbulent flow. The performance of this lossy data compression scheme is evaluated not only with unseen data from the stationary, isotropic turbulent flow, but also with data from decaying isotropic turbulence, and a Taylor-Green vortex flow. Defining the compression ratio (CR) as the ratio of original data size to the compressed one, the results show that our model based on vector quantization can offer CR $=85$ with a mean square error (MSE) of $O(10^{-3})$, and predictions that faithfully reproduce the statistics of the flow, except at the very smallest scales where there is some loss. Compared to the recent study based on a conventional autoencoder where compression is performed in a continuous space, our model improves the CR by more than $30$ percent, and reduces the MSE by an order of magnitude. Our compression model is an attractive solution for situations where fast, high quality and low-overhead encoding and decoding of large data are required.
△ Less
Submitted 11 January, 2022; v1 submitted 10 January, 2022;
originally announced January 2022.
-
Emulating Spatio-Temporal Realizations of Three-Dimensional Isotropic Turbulence via Deep Sequence Learning Models
Authors:
Mohammadreza Momenifar,
Enmao Diao,
Vahid Tarokh,
Andrew D. Bragg
Abstract:
We use a data-driven approach to model a three-dimensional turbulent flow using cutting-edge Deep Learning techniques. The deep learning framework incorporates physical constraints on the flow, such as preserving incompressibility and global statistical invariants of velocity gradient tensor. The accuracy of the model is assessed using statistical and physics-based metrics. The data set comes from…
▽ More
We use a data-driven approach to model a three-dimensional turbulent flow using cutting-edge Deep Learning techniques. The deep learning framework incorporates physical constraints on the flow, such as preserving incompressibility and global statistical invariants of velocity gradient tensor. The accuracy of the model is assessed using statistical and physics-based metrics. The data set comes from Direct Numerical Simulation of an incompressible, statistically stationary, isotropic turbulent flow in a cubic box. Since the size of the dataset is memory intensive, we first generate a low-dimensional representation of the velocity data, and then pass it to a sequence prediction network that learns the spatial and temporal correlations of the underlying data. The dimensionality reduction is performed via extraction using Vector-Quantized Autoencoder (VQ-AE), which learns the discrete latent variables. For the sequence forecasting, the idea of Transformer architecture from natural language processing is used, and its performance compared against more standard Recurrent Networks (such as Convolutional LSTM). These architectures are designed and trained to perform a sequence to sequence multi-class classification task in which they take an input sequence with a fixed length (k) and predict a sequence with a fixed length (p), representing the future time instants of the flow. Our results for the short-term predictions show that the accuracy of results for both models deteriorates across predicted snapshots due to autoregressive nature of the predictions. Based on our diagnostics tests, the trained Conv-Transformer model outperforms the Conv-LSTM one and can accurately, both quantitatively and qualitatively, retain the large scales and capture well the inertial scales of flow but fails at recovering the small and intermittent fluid motions.
△ Less
Submitted 6 December, 2021;
originally announced December 2021.
-
Analysis of scale-dependent kinetic and potential energy in sheared, stably stratified turbulence
Authors:
Xiaolong Zhang,
Rohit Dhariwal,
Gavin Portwood,
Stephen M. de Bruyn Kops,
Andrew D. Bragg
Abstract:
Budgets of turbulent kinetic energy (TKE) and turbulent potential energy (TPE) at different scales $\ell$ in sheared, stably stratified turbulence are analyzed using a filtering approach. Competing effects in the flow are considered, along with the physical mechanisms governing the energy fluxes between scales, and the budgets are used to analyze data from direct numerical simulation (DNS) at buoy…
▽ More
Budgets of turbulent kinetic energy (TKE) and turbulent potential energy (TPE) at different scales $\ell$ in sheared, stably stratified turbulence are analyzed using a filtering approach. Competing effects in the flow are considered, along with the physical mechanisms governing the energy fluxes between scales, and the budgets are used to analyze data from direct numerical simulation (DNS) at buoyancy Reynolds number $Re_b=O(100)$. The mean TKE exceeds TPE by an order of magnitude at the large scales, with the difference reducing as $\ell$ is decreased. At larger scales, buoyancy is never observed to be positive, with buoyancy always converting TKE to TPE. As $\ell$ is decreased, the probability of locally convecting regions increases, though it remains small at scales down to the Ozmidov scale. The TKE and TPE fluxes between scales are both downscale on average and their instantaneous values are positively correlated, but not strongly so, and this occurs due to the different physical mechanisms that govern these fluxes. Moreover, the contribution to these fluxes arising from the sub-grid fields are shown to be significant, in addition to the filtered scale contributions associated with the processes of strain-self amplification, vortex stretching, and density gradient amplification. Probability density functions (PDFs) of the $Q,R$ invariants of the filtered velocity gradient are considered. Unlike isotropic turbulence, as $\ell$ increases the sheared-drop shape of the PDF disappears and the PDF becomes symmetric about $R=0$, meaning regions of vortex stretching and compression become equi-probable, as well as regions of strain amplification or suppression.
△ Less
Submitted 19 October, 2021;
originally announced October 2021.
-
An experimental study on the multiscale properties of turbulence in bubble-laden flows
Authors:
Tian Ma,
Hendrik Hessenkemper,
Dirk Lucas,
Andrew D. Bragg
Abstract:
The properties of bubble-laden turbulent flows at different scales are investigated experimentally, focusing on the flow kinetic energy, energy transfer, and extreme events. The experiments employed particle shadow velocimetry measurements to measure the flow in a column generated by a homogeneous bubble swarm rising in water, for two different bubble diameters ($2.7$ mm $\&$ $3.9$ mm) and moderat…
▽ More
The properties of bubble-laden turbulent flows at different scales are investigated experimentally, focusing on the flow kinetic energy, energy transfer, and extreme events. The experiments employed particle shadow velocimetry measurements to measure the flow in a column generated by a homogeneous bubble swarm rising in water, for two different bubble diameters ($2.7$ mm $\&$ $3.9$ mm) and moderate gas volume fractions ($0.26\%\sim1.31\%$). The two velocity components were measured at high-resolution, and used to construct structure functions up to twelfth order for separations spanning the small to large scales in the flow. Concerning the flow anisotropy, the velocity structure functions are found to differ for separations in the vertical and horizontal directions of the flow, and the cases with smaller bubbles are the most anisotropic, with a dependence on void fraction. The degree of anisotropy is shown to increase as the order of the structure functions is increased, showing that extreme events in the flow are the most anisotropic. Our results show that the average energy transfer with the horizontal velocity component is downscale, just as for the three-dimensional single-phase turbulence. However, the energy transfer associated with the vertical component of the fluid velocity is upscale. The probability density functions of the velocity increments reveal that extreme values become more probable with decreasing Reynolds number, the opposite of the behaviour in single-phase turbulence. We visualize those extreme events and find that regions of intense small scale velocity increments occur near the turbulent/non-turbulent interface at the boundary of the bubble wake.
△ Less
Submitted 23 September, 2021;
originally announced September 2021.
-
A co-spectral budget model links turbulent eddies to suspended sediment concentration in channel flows
Authors:
Shuolin Li,
Andrew D Bragg,
Gabriel Katul
Abstract:
The vertical distribution of suspended sediment concentration (SSC) remains a subject of active research given its relevance to a plethora of problems in hydraulics, hydrology, ecology, and water quality control. Much of the classical theories developed over the course of 90 years represent the effects of turbulence on suspended sediments (SS) using an effective mixing length or eddy diffusivity w…
▽ More
The vertical distribution of suspended sediment concentration (SSC) remains a subject of active research given its relevance to a plethora of problems in hydraulics, hydrology, ecology, and water quality control. Much of the classical theories developed over the course of 90 years represent the effects of turbulence on suspended sediments (SS) using an effective mixing length or eddy diffusivity without explicitly accounting for the energetics of turbulent eddies across scales. To address this gap, the turbulent flux of sediments is derived using a co-spectral budget (CSB) model that can be imminently used in SS and other fine particle transport models. The CSB closes the pressure-redistribution effect using a spectral linear Rotta scheme modified to include isotropoziation of production and interactions between turbulent eddies and sediment grains through a modified scale-dependent de-correlation time. The result is a formulation similar in complexity to the widely used Rouse's equation but with all characteristic scales, Reynolds number, and Schmidt number effects derived from well-established spectral shapes of the vertical velocity and accepted constants from turbulence models. Finally, the proposed CSB model can recover Prandtl's and Rouse's equations under restricted conditions.
△ Less
Submitted 23 August, 2021;
originally announced August 2021.
-
Hydrodynamic interactions and extreme particle clustering in turbulence
Authors:
Andrew D. Bragg,
Adam L. Hammond,
Rohit Dhariwal,
Hui Meng
Abstract:
From new detailed experimental data, we found that the Radial Distribution Function (RDF) of inertial particles in turbulence grows explosively with $r^{-6}$ scaling as the collision radius is approached. We corrected a theory by Yavuz et al. (Phys. Rev. Lett. 120, 244504 (2018)) based on hydrodynamic interactions between pairs of weakly inertial particles, and demonstrate that even this corrected…
▽ More
From new detailed experimental data, we found that the Radial Distribution Function (RDF) of inertial particles in turbulence grows explosively with $r^{-6}$ scaling as the collision radius is approached. We corrected a theory by Yavuz et al. (Phys. Rev. Lett. 120, 244504 (2018)) based on hydrodynamic interactions between pairs of weakly inertial particles, and demonstrate that even this corrected theory cannot explain the observed RDF behavior. We explore several alternative mechanisms for the discrepancy that were not included in the theory and show that none of them are likely the explanation, suggesting new, yet to be identified physical mechanisms are at play.
△ Less
Submitted 6 April, 2021;
originally announced April 2021.
-
Scale-dependent anisotropy, energy transfer and intermittency in bubble-laden turbulent flows
Authors:
Tian Ma,
Bernhard Ott,
Jochen Frohlich,
Andrew D. Bragg
Abstract:
Data from Direct Numerical Simulations of disperse bubbly flows in a vertical channel are used to study the effect of the bubbles on the carrier-phase turbulence. A new method is developed, based on the barycentric map approach, that allows to quantify the anisotropy and componentiality of the flow at any scale. Using this the bubbles are found to significantly enhance flow anisotropy at all scale…
▽ More
Data from Direct Numerical Simulations of disperse bubbly flows in a vertical channel are used to study the effect of the bubbles on the carrier-phase turbulence. A new method is developed, based on the barycentric map approach, that allows to quantify the anisotropy and componentiality of the flow at any scale. Using this the bubbles are found to significantly enhance flow anisotropy at all scales compared with the unladen case, and for some bubble cases, very strong anisotropy persists down to the smallest flow scales. The strongest anisotropy observed was for the cases involving small bubbles. Concerning the inter-scale energy transfer, our results indicate that for the bubble-laden cases, the energy transfer is from large to small scales, just as for the unladen case. However, there is evidence of an upscale transfer when considering the transfer of energy associated with particular components of the velocity field. Although the direction of the energy transfer is the same with and without the bubbles, the transfer is much stronger for the bubble-laden cases, suggesting that the bubbles play a strong role in enhancing the activity of the nonlinear term in the flow. The normalized forms of the fourth and sixth-order structure functions are also considered, and reveal that the introduction of bubbles into the flow strongly enhances intermittency in the dissipation range, but suppresses it at larger scales. This strong enhancement of the dissipation scale intermittency has significant implications for understanding how the bubbles might modify the mixing properties of turbulent flows.
△ Less
Submitted 1 April, 2021;
originally announced April 2021.
-
Dimension Reduced Turbulent Flow Data From Deep Vector Quantizers
Authors:
Mohammadreza Momenifar,
Enmao Diao,
Vahid Tarokh,
Andrew D. Bragg
Abstract:
Analyzing large-scale data from simulations of turbulent flows is memory intensive, requiring significant resources. This major challenge highlights the need for data compression techniques. In this study, we apply a physics-informed Deep Learning technique based on vector quantization to generate a discrete, low-dimensional representation of data from simulations of three-dimensional turbulent fl…
▽ More
Analyzing large-scale data from simulations of turbulent flows is memory intensive, requiring significant resources. This major challenge highlights the need for data compression techniques. In this study, we apply a physics-informed Deep Learning technique based on vector quantization to generate a discrete, low-dimensional representation of data from simulations of three-dimensional turbulent flows. The deep learning framework is composed of convolutional layers and incorporates physical constraints on the flow, such as preserving incompressibility and global statistical characteristics of the velocity gradients. The accuracy of the model is assessed using statistical, comparison-based similarity and physics-based metrics. The training data set is produced from Direct Numerical Simulation of an incompressible, statistically stationary, isotropic turbulent flow.
The performance of this lossy data compression scheme is evaluated not only with unseen data from the stationary, isotropic turbulent flow, but also with data from decaying isotropic turbulence, a Taylor-Green vortex flow, and a turbulent channel flow. Defining the compression ratio (CR) as the ratio of original data size to the compressed one, the results show that our model based on vector quantization can offer CR$=85$ with a mean square error (MSE) of $O(10^{-3})$, and predictions that faithfully reproduce the statistics of the flow, except at the very smallest scales where there is some loss. Compared to the recent study of Glaws. et. al. (Physical Review Fluids, 5(11):114602, 2020), which was based on a conventional autoencoder (where compression is performed in a continuous space), our model improves the CR by more than $30$ percent...
△ Less
Submitted 24 May, 2021; v1 submitted 1 March, 2021;
originally announced March 2021.
-
When is settling important for particle concentrations in wall-bounded turbulent flows?
Authors:
Andrew D. Bragg,
David H. Richter,
Guiquan Wang
Abstract:
We explore the role of gravitational settling on inertial particle concentrations in a wall-bounded turbulent flow. While it may be thought that settling can be ignored when the settling parameter $Sv\equiv v_s/u_τ$ is small ($v_s$ - Stokes settling velocity, $u_τ$ - fluid friction velocity), we show that even in this regime the settling may make a leading order contribution to the concentration p…
▽ More
We explore the role of gravitational settling on inertial particle concentrations in a wall-bounded turbulent flow. While it may be thought that settling can be ignored when the settling parameter $Sv\equiv v_s/u_τ$ is small ($v_s$ - Stokes settling velocity, $u_τ$ - fluid friction velocity), we show that even in this regime the settling may make a leading order contribution to the concentration profiles. This is because the importance of settling is determined, not by the size of $v_s$ compared with $u_τ$ or any other fluid velocity scale, but by the size of $v_s$ relative to the other mechanisms that control the vertical particle velocity and concentration profile. We explain this in the context of the particle mean-momentum equation, and show that in general, there always exists a region in the boundary layer where settling cannot be neglected, no matter how small $Sv$ is (provided it is finite). Direct numerical simulations confirm the arguments, and show that the near-wall concentration is highly dependent on $Sv$ even when $Sv\ll 1$, and can reduce by an order of magnitude when $Sv$ is increased from $O(10^{-4})$ and $O(10^{-2})$. The results also show that the preferential sampling of ejection events in the boundary layer by inertial particles when $Sv=0$ is profoundly altered as $Sv$ is increased, and is replaced by a preferential sampling of sweep events due to the onset of the preferential sweeping mechanism.
△ Less
Submitted 9 September, 2021; v1 submitted 12 January, 2021;
originally announced January 2021.
-
Turbulence and heat transfer on a rotating, heated half soap bubble
Authors:
X. Q. HE,
A. D. Bragg,
Y. L. Xiong,
P. Fischer
Abstract:
We use Direct Numerical Simulations to study the two-dimensional flow of a rotating, half soap bubble that is heated at its equator. The heating produces buoyancy and rotation generates a Coriolis forces in the fluid. However, due to the curved surface of the bubble, the buoyancy and Coriolis forces vary with latitude on the bubble, giving rise to rich flow behavior. We first explore the single-po…
▽ More
We use Direct Numerical Simulations to study the two-dimensional flow of a rotating, half soap bubble that is heated at its equator. The heating produces buoyancy and rotation generates a Coriolis forces in the fluid. However, due to the curved surface of the bubble, the buoyancy and Coriolis forces vary with latitude on the bubble, giving rise to rich flow behavior. We first explore the single-point properties of the flow, including the Reynolds and Nusselt numbers, mean fields, and Reynolds stresses, all as a function of latitude. For a given Rayleigh number, we observe a non-monotonic dependence on the Rossby number Ro, and large scale mean circulations that are strongly influenced by rotation. We then consider quantities that reveal the multiscale nature of the flow, including spectrums and spectral fluxes of kinetic and thermal energy, and enstrophy, and structure functions of velocity and temperature. The fluxes show that just a for nonbuoyant two-dimensional turbulence on a flat surface, there is an upscale flux of kinetic energy at larger scales (fed by buoyancy injection of turbulent kinetic energy at smaller scales), and a downscale flux of enstrophy at smaller scales. The kinetic energy spectrum and velocity structure functions are well described by Bolgiano-Obukhov (BO) scaling at scales where the effects of rotation are weak. The temperature structure functions do not, however, satisfy BO scaling in general, due to strong intermittency in the temperature field.
△ Less
Submitted 18 November, 2020;
originally announced November 2020.
-
Mechanisms governing the settling velocities and spatial distributions of inertial particles in wall-bounded turbulence
Authors:
Andrew D Bragg,
David H Richter,
Guiquan Wang
Abstract:
We use theory and Direct Numerical Simulations (DNS) to explore the average vertical velocities and spatial distributions of inertial particles settling in a wall-bounded turbulent flow. The theory is based on the exact phase-space equation for the Probability Density Function describing particle positions and velocities. This allowed us to identify the distinct physical mechanisms governing the p…
▽ More
We use theory and Direct Numerical Simulations (DNS) to explore the average vertical velocities and spatial distributions of inertial particles settling in a wall-bounded turbulent flow. The theory is based on the exact phase-space equation for the Probability Density Function describing particle positions and velocities. This allowed us to identify the distinct physical mechanisms governing the particle transport. We then examined the asymptotic behavior of the particle motion near the wall, revealing the fundamental differences to the near wall behavior that is produced when incorporating gravitational settling. When the average vertical particle mass flux is zero, the averaged vertical particle velocity is zero away from the wall due to the particles preferentially sampling regions where the fluid velocity is positive, which balances with the downward Stokes settling velocity. When the average mass flux is negative, the combined effects of turbulence and particle inertia lead to average vertical particle velocities that can significantly exceed the Stokes settling velocity, by as much as ten times. Sufficiently far from the wall, the enhanced vertical velocities are due to the preferential sweeping mechanism. However, as the particles approach the wall, the contribution from the preferential sweeping mechanism becomes small, and a downward contribution from the turbophoretic velocity dominates the behavior. Close to the wall, the particle concentration grows as a power-law, but the nature of this power law depends on the particle Stokes number. Finally, our results highlight how the Rouse model of particle concentration is to be modified for particles with finite inertia.
△ Less
Submitted 6 October, 2020;
originally announced October 2020.
-
Exploring the turbulent velocity gradients at different scales from the perspective of the strain-rate eigenframe
Authors:
Josin Tom,
Maurizio Carbone,
Andrew D Bragg
Abstract:
Expressing the evolution equations for the filtered velocity gradient tensor (FVGT) in the strain-rate eigenframe provides an insightful way to disentangle and understand various processes such as strain self-amplification, vortex stretching and tilting, and to consider their properties at different scales in the flow. Using data from Direct Numerical Simulation (DNS) of the forced Navier-Stokes e…
▽ More
Expressing the evolution equations for the filtered velocity gradient tensor (FVGT) in the strain-rate eigenframe provides an insightful way to disentangle and understand various processes such as strain self-amplification, vortex stretching and tilting, and to consider their properties at different scales in the flow. Using data from Direct Numerical Simulation (DNS) of the forced Navier-Stokes equation, we consider the relative importance of local and non-local terms in the FVGT eigenframe equations across the scales using statistical analysis. The analysis of the eigenframe rotation-rate, that drives vorticity tilting, shows that the anisotropic pressure Hessian plays a key role, with the sub-grid stress making an important contribution outside the dissipation range, and the local spinning due to vorticity making a much smaller contribution. The results also show the striking behavior that the vorticity tilting term remains highly intermittent even at relatively large scales. We derive a generalization of the Lumley triangle that allows us to show that the pressure Hessian has a preference for two-component axisymmetric configurations at small scales, with a transition to a more isotropic state at larger scales. Correlations between the sub-grid stress and other terms in the eigenframe equations are considered, highlighting the coupling between the sub-grid and nonlinear amplification terms, with the sub-grid term playing an important role in regularizing the system. These results provide useful guidelines for improving Lagrangian models of the FVGT, since current models fail to capture a number of subtle features observed in our results.
△ Less
Submitted 8 May, 2020;
originally announced May 2020.
-
Fluctuation-theorem and extended thermodynamics of turbulence
Authors:
Amilcare Porporato,
Milad Hooshyar,
Andrew D Bragg,
Gabriel Katul
Abstract:
Turbulent flows are out-of-equilibrium because the energy supply at large scales and its dissipation by viscosity at small scales create a net transfer of energy among all scales. Here, the energy cascade is approximated by a combined contribution of a forward drift and diffusion that recover accepted phenomenological theories of turbulence. The fluctuation theorem (FT) is then shown to describe t…
▽ More
Turbulent flows are out-of-equilibrium because the energy supply at large scales and its dissipation by viscosity at small scales create a net transfer of energy among all scales. Here, the energy cascade is approximated by a combined contribution of a forward drift and diffusion that recover accepted phenomenological theories of turbulence. The fluctuation theorem (FT) is then shown to describe the scale-wise statistics of forward and backward energy transfer and their connection to irreversibility and entropy production. The ensuing turbulence entropy may be used to formulate an extended turbulence thermodynamics.
△ Less
Submitted 4 May, 2020;
originally announced May 2020.
-
Gauge symmetry and dimensionality reduction of the anisotropic pressure Hessian
Authors:
Maurizio Carbone,
Michele Iovieno,
Andrew D Bragg
Abstract:
Analyzing the fluid velocity gradients in a Lagrangian reference frame provides an insightful way to study the small-scale dynamics of turbulent flows, and further insight is provided by considering the equations in the eigenframe of the strain-rate tensor. The dynamics of the velocity gradient tensor is governed in part by the anisotropic pressure Hessian, which is a non-local functional of the v…
▽ More
Analyzing the fluid velocity gradients in a Lagrangian reference frame provides an insightful way to study the small-scale dynamics of turbulent flows, and further insight is provided by considering the equations in the eigenframe of the strain-rate tensor. The dynamics of the velocity gradient tensor is governed in part by the anisotropic pressure Hessian, which is a non-local functional of the velocity gradient field. This anisotropic pressure Hessian plays a key role in the velocity gradient dynamics, for example in preventing finite-time singularities, but it is difficult to understand and model due to its non-locality and complexity. In this work a gauge symmetry for the pressure Hessian is introduced to the eigenframe equations of the velocity gradient, such that when the gauge is added to the original pressure Hessian, the dynamics of the eigenframe variables remain unchanged. We then exploit this gauge symmetry to perform a rank reduction on the three-dimensional anisotropic pressure Hessian, which, remarkably, is possible everywhere in the flow. The dynamical activity of the newly introduced rank-reduced anisotropic pressure Hessian is confined to two dimensional manifolds in the three dimensional flow, and exhibits striking alignment properties with respect to the strain-rate eigenframe and the vorticity vector. The dimensionality reduction, together with the strong preferential alignment properties, leads to new dynamical insights for understanding and modelling the role of the anisotropic pressure Hessian in three-dimensional flows.
△ Less
Submitted 19 November, 2019;
originally announced November 2019.
-
Local analysis of the clustering, velocities and accelerations of particles settling in turbulence
Authors:
Mohammadreza Momenifar,
Andrew D. Bragg
Abstract:
Using 3D Vorono\text{ï} analysis, we explore the local dynamics of small, settling, inertial particles in isotropic turbulence using Direct Numerical Simulations (DNS). We independently vary the Taylor Reynolds number $R_λ\in[90,398]$, Froude number $Fr\equiv a_η/g\in[0.052,\infty]$ (where $a_η$ is the Kolmogorov acceleration, and $g$ is the acceleration due to gravity), and Kolmogorov scale Stoke…
▽ More
Using 3D Vorono\text{ï} analysis, we explore the local dynamics of small, settling, inertial particles in isotropic turbulence using Direct Numerical Simulations (DNS). We independently vary the Taylor Reynolds number $R_λ\in[90,398]$, Froude number $Fr\equiv a_η/g\in[0.052,\infty]$ (where $a_η$ is the Kolmogorov acceleration, and $g$ is the acceleration due to gravity), and Kolmogorov scale Stokes number $St\equivτ_p/τ_η\in[0,3]$. In agreement with previous results using global measures of particle clustering, such as the Radial Distribution Function (RDF), we find that for small Vorono\text{ï} volumes (corresponding to the most clustered particles), the behavior is strongly dependent upon $St$ and $Fr$, but only weakly dependent upon $R_λ$, unless $St>1$. However, larger Vorono\text{ï} volumes (void regions) exhibit a much stronger dependence on $R_λ$, even when $St\leq 1$, and we show that this, rather than the behavior at small volumes, is the cause of the sensitivity of the standard deviation of the Vorono\text{ï} volumes that has been previously reported. We also show that the largest contribution to the particle settling velocities is associated with increasingly larger Vorono\text{ï} volumes as the settling parameter $Sv\equiv St/Fr$ is increased. Our local analysis of the acceleration statistics of settling inertial particles shows that clustered particles experience a net acceleration in the direction of gravity, while particles in void regions experience the opposite. The particle acceleration variance, however, is a convex function of the Vorono\text{ï} volumes, with or without gravity, which seems to indicate a non-trivial relationship between the Vorono\text{ï} volumes and the sizes of the turbulent flow scales. Results for the variance of the fluid acceleration at the inertial particle "..."
△ Less
Submitted 1 August, 2019;
originally announced August 2019.
-
Is vortex stretching the main cause of the turbulent energy cascade?
Authors:
Maurizio Carbone,
Andrew D. Bragg
Abstract:
In three dimensional turbulence there is on average a cascade of kinetic energy from the largest to the smallest scales of the flow. While the dominant idea is that the cascade occurs through the physical process of vortex stretching, evidence for this is debated. In the framework of the Karman-Howarth equation for the two point turbulent kinetic energy, we derive a new result for the average flux…
▽ More
In three dimensional turbulence there is on average a cascade of kinetic energy from the largest to the smallest scales of the flow. While the dominant idea is that the cascade occurs through the physical process of vortex stretching, evidence for this is debated. In the framework of the Karman-Howarth equation for the two point turbulent kinetic energy, we derive a new result for the average flux of kinetic energy between two points in the flow that reveals the role of vortex stretching. However, the result shows that vortex stretching is in fact not the main contributor to the average energy cascade; the main contributor is the self-amplification of the strain-rate field. We emphasize the need to correctly distinguish and not conflate the roles of vortex stretching and strain-self amplification in order to correctly understand the physics of the cascade, and also resolve a paradox regarding the differing role of vortex stretching on the mechanisms of the energy cascade and energy dissipation rate. Direct numerical simulations are used to confirm the results, as well as provide further results and insights on vortex stretching and strain-self amplification at different scales in the flow. Interestingly, the results imply that while vortex stretching plays a sub-leading role in the average cascade, it may play a leading order role during large fluctuations of the energy cascade about its average behavior.
△ Less
Submitted 17 June, 2019;
originally announced June 2019.
-
Multiscale preferential sweeping of particles settling in turbulence
Authors:
Josin Tom,
Andrew D Bragg
Abstract:
In a seminal article, \citet[J. Fluid Mech., 174:441-465]{maxey87} presented a theoretical analysis showing that enhanced particle settling speeds in turbulence occur through the preferential sweeping mechanism, which depends on the preferential sampling of the fluid velocity gradient field by the inertial particles. However, recent Direct Numerical Simulation (DNS) results in \citet[J. Fluid Mech…
▽ More
In a seminal article, \citet[J. Fluid Mech., 174:441-465]{maxey87} presented a theoretical analysis showing that enhanced particle settling speeds in turbulence occur through the preferential sweeping mechanism, which depends on the preferential sampling of the fluid velocity gradient field by the inertial particles. However, recent Direct Numerical Simulation (DNS) results in \citet[J. Fluid Mech., 796:659--711]{ireland16b} show that even in a portion of the parameter space where this preferential sampling is absent, the particles nevertheless exhibit enhanced settling velocities. Further, there are several outstanding questions concerning the role of different turbulent flow scales on the enhanced settling, and the role of the Taylor Reynolds number $R_λ$. The analysis of Maxey does not explain these issues, partly since it was restricted to particle Stokes numbers $St\ll1$. To address these issues, we have developed a new theoretical result, valid for arbitrary $St$, that reveals the multiscale nature of the mechanism generating the enhanced settling speeds. In particular, it shows how the range of scales at which the preferential sweeping mechanism operates depends on $St$. This analysis is complemented by results from DNS where we examine the role of different flow scales on the particle settling speeds by coarse-graining the underlying flow. The results show how the flow scales that contribute to the enhanced settling depend on $St$, and that contrary to previous claims, there can be no single turbulent velocity scale that characterizes the enhanced settling speed. The results explain the dependence of the particle settling speeds on $R_λ$, and show how the saturation of this dependence at sufficiently large $R_λ$ depends upon $St$. The results also show ...
△ Less
Submitted 20 April, 2019; v1 submitted 20 December, 2018;
originally announced December 2018.
-
Multiscale fluid--particle thermal interaction in isotropic turbulence
Authors:
Maurizio Carbone,
Andrew D. Bragg,
Michele Iovieno
Abstract:
We use direct numerical simulations to investigate the interaction between the temperature field of a fluid and the temperature of small particles suspended in the flow, employing both one and two-way thermal coupling, in a statistically stationary, isotropic turbulent flow. Using statistical analysis, we investigate this variegated interaction at the different scales of the flow. We find that the…
▽ More
We use direct numerical simulations to investigate the interaction between the temperature field of a fluid and the temperature of small particles suspended in the flow, employing both one and two-way thermal coupling, in a statistically stationary, isotropic turbulent flow. Using statistical analysis, we investigate this variegated interaction at the different scales of the flow. We find that the variance of the fluid temperature gradients decreases as the thermal response time of the suspended particles is increased. The probability density function (PDF) of the fluid temperature gradients scales with its variance, while the PDF of the rate of change of the particle temperature, whose variance is associated with the thermal dissipation due to the particles, does not scale in such a self-similar way. The modification of the fluid temperature field due to the particles is examined by computing the particle concentration and particle heat fluxes conditioned on the magnitude of the local fluid temperature gradient. These statistics highlight that the particles cluster on the fluid temperature fronts, and the important role played by the alignments of the particle velocity and the local fluid temperature gradient. The temperature structure functions, which characterize the temperature fluctuations across the scales of the flow, clearly show that the fluctuations of the fluid temperature increments are monotonically suppressed in the two-way coupled regime as the particle thermal response time is increased. Thermal caustics dominate the particle temperature increments at small scales, that is, particles that come into contact are likely to have very large differences in their temperature. This is caused by the nonlocal thermal dynamics of the particles...
△ Less
Submitted 22 October, 2018;
originally announced October 2018.
-
Enhanced and suppressed multiscale dispersion of bidisperse inertial particles due to gravity
Authors:
Rohit Dhariwal,
Andrew D. Bragg
Abstract:
Using Direct Numerical Simulations (DNS), we investigate how gravity modifies the multiscale dispersion of bidisperse inertial particles in isotropic turbulence. The DNS has a Taylor Reynolds number $R_λ=398$, and we simulate Stokes numbers (based on the Kolmogorov timescale) in the range $St\leq 3$ , and consider Froude numbers $Fr = 0.052$ and $\infty$, corresponding to strong gravity and no gra…
▽ More
Using Direct Numerical Simulations (DNS), we investigate how gravity modifies the multiscale dispersion of bidisperse inertial particles in isotropic turbulence. The DNS has a Taylor Reynolds number $R_λ=398$, and we simulate Stokes numbers (based on the Kolmogorov timescale) in the range $St\leq 3$ , and consider Froude numbers $Fr = 0.052$ and $\infty$, corresponding to strong gravity and no gravity, respectively. The degree of bidispersity is quantified by the difference in the Stokes number of the particles $|ΔSt|$. We first consider the mean-square separation of bidisperse particle-pairs and find that without gravity (i.e. $Fr = \infty$), bidispersity leads to an enhancement of the the mean-square separation over a significant range of scales. When $|ΔSt|\geq O(1)$, the relative dispersion is further enhanced by gravity due to the large difference in the settling velocities of the two particles. However, when $|ΔSt|\ll1$, gravity suppresses the relative dispersion as the settling velocity contribution is small, and gravity suppresses the non-local contribution to the particle dynamics. In order to gain further insights, we consider separately the relative dispersion in the vertical (parallel to gravity) and horizontal directions. As expected, the vertical relative dispersion can be strongly enhanced by gravity due to differences in the settling velocities of the two particles. However, a key finding of our study is that gravity can also significantly enhance the horizontal relative dispersion. This non-trivial effect occurs because fast settling particles experience rapid fluctuations in the fluid velocity field along their trajectory, leading to enhanced particle accelerations and relative velocities. For sufficiently large initial particle separations...
△ Less
Submitted 5 August, 2018;
originally announced August 2018.
-
The influence of Reynolds and Froude number on the motion of settling, bidisperse inertial particles in turbulence
Authors:
Mohammadreza Momenifar,
Rohit Dhariwal,
Andrew D. Bragg
Abstract:
Using Direct Numerical Simulations (DNS), we examine the effects of Taylor Reynolds number, $R_λ$, and Froude number, $Fr$, on the motion of settling, bidisperse inertial particles in isotropic turbulence. Particle accelerations play a key role in the relative motion of bidisperse particles, and we find that reducing $Fr$ leads to an enhancement of the accelerations, but a suppression of their int…
▽ More
Using Direct Numerical Simulations (DNS), we examine the effects of Taylor Reynolds number, $R_λ$, and Froude number, $Fr$, on the motion of settling, bidisperse inertial particles in isotropic turbulence. Particle accelerations play a key role in the relative motion of bidisperse particles, and we find that reducing $Fr$ leads to an enhancement of the accelerations, but a suppression of their intermittency. For Stokes numbers $St>1$, the effect of $R_λ$ on the accelerations is enhanced by gravity, since settling causes the particle accelerations to be affected by a larger range of flow scales. The results for the Probability Density Function (PDF) of the particle relative velocities show that for bidisperse particles, decreasing $Fr$ leads to an enhancement of their relative velocities in both the vertical (parallel to gravity) and horizontal directions. Importantly, our results show that even when the particles are settling very fast, turbulence continues to play a key role in their vertical relative velocities, and increasingly so as $R_λ$ is increased. This occurs because although the settling velocity may be much larger than typical velocities of the turbulence, due to intermittency, there are significant regions of the flow where the turbulence contribution to the particle motion is of the same order as that from gravitational settling...
△ Less
Submitted 4 August, 2018;
originally announced August 2018.
-
Particle-Pair Relative Velocity Measurement in High-Reynolds-Number Homogeneous and Isotropic Turbulence Using 4-Frame Particle Tracking Velocimetry
Authors:
Zhongwang Dou,
Peter J. Ireland,
Andrew D. Bragg,
Zach Liang,
Lance R. Collins,
Hui Meng
Abstract:
The radial relative velocity between particles suspended in turbulent flow plays a critical role in droplet collision and growth. We present a simple and accurate approach to RV measurement in isotropic turbulence - planar 4-frame particle tracking velocimetry - using routine PIV hardware. This study demonstrates the feasibility of accurately measuring RV using routine hardware and verifies, for t…
▽ More
The radial relative velocity between particles suspended in turbulent flow plays a critical role in droplet collision and growth. We present a simple and accurate approach to RV measurement in isotropic turbulence - planar 4-frame particle tracking velocimetry - using routine PIV hardware. This study demonstrates the feasibility of accurately measuring RV using routine hardware and verifies, for the first time, the path-history and inertial filtering effects on particle-pair RV at large particle separations experimentally.
△ Less
Submitted 20 December, 2017;
originally announced December 2017.
-
Do fluid particles separate exponentially in the dissipation range?
Authors:
Rohit Dhariwal,
Andrew D. Bragg
Abstract:
In this paper we consider how the statistical moments of the separation between two fluid particles grow with time when their separation lies in the dissipation range of turbulence. In this range the fluid velocity field varies smoothly, and the relative velocity of two fluid particles depends linearly upon their separation. While this may suggest that the rate at which fluid particles separate is…
▽ More
In this paper we consider how the statistical moments of the separation between two fluid particles grow with time when their separation lies in the dissipation range of turbulence. In this range the fluid velocity field varies smoothly, and the relative velocity of two fluid particles depends linearly upon their separation. While this may suggest that the rate at which fluid particles separate is exponential in time, this is not guaranteed because the strain-rate governing their separation is a strongly fluctuating quantity in turbulence. Indeed, the recent paper by Afik \& Steinberg (Nat. Commun. \textbf{8}: 468, 2017) argues that there is no convincing evidence that the moments of the separation between fluid particles grow exponentially with time in the dissipation range of turbulence. Motivated by this, we use Direct Numerical Simulations (DNS) to compute the moments of the particle separation over very long periods of time to see if we ever see evidence for exponential separation. Our results show that if the initial separation between the particles is infinitesimal, the moments of the particle separation first grow as power laws in time, but we then observe, for the first time, convincing evidence that at sufficiently long times the moments do grow exponentially. However, this exponential growth is only observed after extremely long times $\gtrsim 200τ_η$, where $τ_η$ is the Kolmogorov timescale. This is due to fluctuations in the strain-rate about its mean value measured along the particle trajectories, the effect of which on the moments of the particle separation persists for very long times. We also consider the Backward-in-Time (BIT) moments of the particle separation, and observe that they too grow exponentially in the long-time regime. However, a dramatic consequence of the exponential separation is that at long-times the difference between..
△ Less
Submitted 11 December, 2017;
originally announced December 2017.
-
Effects of Reynolds Number and Stokes Number on Particle-pair Relative Velocity in Isotropic Turbulence: A Systematic Experimental Study
Authors:
Zhongwang Dou,
Andrew D. Bragg,
Adam L. Hammond,
Zach Liang,
Lance R. Collins,
Hui Meng
Abstract:
The effects of Reynolds number and Stokes number on particle-pair relative velocity (RV) were investigated systematically using a recently developed planar four-frame particle tracking technique in a novel homogeneous and isotropic turbulence chamber.
The effects of Reynolds number and Stokes number on particle-pair relative velocity (RV) were investigated systematically using a recently developed planar four-frame particle tracking technique in a novel homogeneous and isotropic turbulence chamber.
△ Less
Submitted 6 November, 2017;
originally announced November 2017.
-
Extremes, intermittency and time reversibility of atmospheric turbulence at the cross-over from production to inertial scales
Authors:
E. Zorzetto,
A. D. Bragg,
G. Katul
Abstract:
The effects of mechanical generation of turbulent kinetic energy and buoyancy forces on the statistics of air temperature and velocity increments are experimentally investigated at the cross over from production to inertial range scales. The ratio of an approximated mechanical to buoyant production (or destruction) of turbulent kinetic energy can be used to form a dimensionless stability parameter…
▽ More
The effects of mechanical generation of turbulent kinetic energy and buoyancy forces on the statistics of air temperature and velocity increments are experimentally investigated at the cross over from production to inertial range scales. The ratio of an approximated mechanical to buoyant production (or destruction) of turbulent kinetic energy can be used to form a dimensionless stability parameter $ζ$ that classifies the state of the atmosphere as common in many atmospheric surface layer studies. Here, we assess how $ζ$ affects the scale-wise evolution of the probability of extreme air temperature excursions, their asymmetry and time reversibility. The analysis makes use of high frequency velocity and air temperature time series measurements collected at $z$=5 m above a grass surface at very large Reynolds numbers $Re=u_* z/ν> 1\times 10^5$ ($u_*$ is the friction velocity and $ν$ is the kinematic viscosity of air). Using conventional higher-order structure functions, temperature exhibits larger intermittency and wider multifractality when compared to the longitudinal velocity component, consistent with laboratory studies and simulations conducted at lower $Re$. Moreover, deviations from the classical Kolmogorov scaling for the longitudinal velocity are shown to be reasonably described by the She-Leveque vortex filament model that has no 'tunable' parameters and is independent of $ζ$. The work demonstrates that external boundary conditions, and in particular the magnitude and sign of the sensible heat flux, have a significant impact on temperature advection-diffusion dynamics within the inertial range. In particular, atmospheric stability affects both the buildup of intermittency and the persistent asymmetry and time irreversibility observed in the first two decades of inertial sub-range scales.
△ Less
Submitted 2 December, 2017; v1 submitted 5 October, 2017;
originally announced October 2017.
-
Irreversibility-inversions in 2 dimensional turbulence
Authors:
Andrew D Bragg,
Filippo De Lillo,
Guido Boffetta
Abstract:
In this paper we consider a recent theoretical prediction (Bragg \emph{et al.}, Phys. Fluids \textbf{28}, 013305 (2016)) that for inertial particles in 2D turbulence, the nature of the irreversibility of the particle-pair dispersion inverts when the particle inertia exceeds a certain value. In particular, when the particle Stokes number, ${\rm St}$, is below a certain value, the forward-in-time (F…
▽ More
In this paper we consider a recent theoretical prediction (Bragg \emph{et al.}, Phys. Fluids \textbf{28}, 013305 (2016)) that for inertial particles in 2D turbulence, the nature of the irreversibility of the particle-pair dispersion inverts when the particle inertia exceeds a certain value. In particular, when the particle Stokes number, ${\rm St}$, is below a certain value, the forward-in-time (FIT) dispersion should be faster than the backward-in-time (BIT) dispersion, but for ${\rm St}$ above this value, this should invert so that BIT becomes faster than FIT dispersion. This non-trivial behavior arises because of the competition between two physically distinct irreversibility mechanisms that operate in different regimes of ${\rm St}$. In 3D turbulence, both mechanisms act to produce faster BIT than FIT dispersion, but in 2D turbulence, the two mechanisms have opposite effects because of the flux of energy from the small to the large scales. We supplement the qualitative argument given by Bragg \emph{et al.} (Phys. Fluids \textbf{28}, 013305 (2016)) by deriving quantitative predictions of this effect in the short time limit. We confirm the theoretical predictions using results of inertial particle dispersion in a direct numerical simulation of 2D turbulence. A more general finding of this analysis is that in turbulent flows with an inverse energy flux, inertial particles may yet exhibit a net downscale flux of kinetic energy because of their non-local in-time dynamics.
△ Less
Submitted 26 September, 2017;
originally announced September 2017.
-
Small-scale dynamics of settling, bidisperse particles in turbulence
Authors:
Rohit Dhariwal,
Andrew D. Bragg
Abstract:
We use Direct Numerical Simulations (DNS) to investigate the dynamics of settling, bidisperse particles in isotropic turbulence. In agreement with previous studies, we find that without gravity (i.e. $Fr=\infty$, where $Fr$ is the Froude number), bidispersity leads to an enhancement of the relative velocities, and a suppression of their spatial clustering. For $Fr<1$, the relative velocities in th…
▽ More
We use Direct Numerical Simulations (DNS) to investigate the dynamics of settling, bidisperse particles in isotropic turbulence. In agreement with previous studies, we find that without gravity (i.e. $Fr=\infty$, where $Fr$ is the Froude number), bidispersity leads to an enhancement of the relative velocities, and a suppression of their spatial clustering. For $Fr<1$, the relative velocities in the direction of gravity are enhanced by the differential settling velocities of the bidisperse particles, as expected. However, we also find that gravity can strongly enhance the relative velocities in the "horizontal" directions (i.e. in the plane normal to gravity). This non-trivial behavior occurs because fast settling particles experience rapid fluctuations in the fluid velocity field along their trajectory, leading to enhanced particle accelerations and relative velocities. Indeed, the results show that even when $Fr\ll1$, turbulence can still play an important role, not only on the horizontal motion, but also on the vertical motion of the particles, with significant implications for understanding the mixing rates of settling bidisperse particles in turbulence. We also find that gravity drastically reduces the clustering of bidisperse particles. These results are strikingly different to the monodisperse case, for which recent results have shown that when $Fr<1$, gravity strongly suppresses the relative velocities in all directions, and can enhance clustering. We then consider the implications of these results for the collision rates of settling, bidisperse particles in turbulence. We find that for $Fr=0.052$, the collision kernel is almost perfectly predicted by the collision kernel for bidisperse particles settling in quiescent flow, such that the effect of turbulence may be ignored...
△ Less
Submitted 1 September, 2017;
originally announced September 2017.
-
Is the kinetic equation for turbulent gas-particle flows ill-posed?
Authors:
M. W. Reeks,
D. C. Swailes,
A. D. Bragg
Abstract:
This paper is about well-posedness and realizability of the kinetic equation for gas-particle flows and its relationship to the Generalized Langevin Model (GLM) PDF equation. Previous analyses claim that this kinetic equation is ill-posed, that in particular it has the properties of a backward heat equation and as a consequence, its solutions will in the course of time exhibit finite-time singular…
▽ More
This paper is about well-posedness and realizability of the kinetic equation for gas-particle flows and its relationship to the Generalized Langevin Model (GLM) PDF equation. Previous analyses claim that this kinetic equation is ill-posed, that in particular it has the properties of a backward heat equation and as a consequence, its solutions will in the course of time exhibit finite-time singularities. We show that the analysis leading to this conclusion is fundamentally incorrect because it ignores the coupling between the phase space variables in the kinetic equation and the time and particle inertia dependence of the phase space diffusion tensor. This contributes an extra $+ve$ diffusion that always outweighs the contribution from the$-ve$ diffusion associated with the dispersion along one of the principal axes of the phase space diffusion tensor. This is confirmed by a numerical evaluation of analytic solutions of these $+ve$ and $-ve$ contributions to the particle diffusion coefficient along this principal axis. We also examine other erroneous claims and assumptions made in previous studies that demonstrate the apparent superiority of the GLM PDF approach over the kinetic approach. In so doing we have drawn attention to the limitations of the GLM approach which these studies have ignored or not properly considered, to give a more balanced appraisal of the benefits of both PDF approaches.
△ Less
Submitted 31 July, 2017; v1 submitted 3 May, 2017;
originally announced May 2017.
-
Analysis of the forward and backward in time pair-separation PDFs for inertial particles in isotropic turbulence
Authors:
Andrew D. Bragg
Abstract:
In this paper we investigate, using theory and Direct Numerical Simulations (DNS), the Forward In Time (FIT) and Backward In Time (BIT) Probability Density Functions (PDFs) of the separation of inertial particle-pairs in isotropic turbulence. In agreement with our earlier study (Bragg et al., Phys. Fluids 28, 013305 (2016)), where we compared the FIT and BIT mean-square separations, we find that i…
▽ More
In this paper we investigate, using theory and Direct Numerical Simulations (DNS), the Forward In Time (FIT) and Backward In Time (BIT) Probability Density Functions (PDFs) of the separation of inertial particle-pairs in isotropic turbulence. In agreement with our earlier study (Bragg et al., Phys. Fluids 28, 013305 (2016)), where we compared the FIT and BIT mean-square separations, we find that inertial particles separate much faster BIT than FIT, with the strength of the irreversibility depending upon the final/initial separation of the particle-pair and their Stokes number $St$. However, we also find that the irreversibility shows up in subtle ways in the behavior of the full PDF that it does not in the mean-square separation. In the theory, we derive new predictions, including a prediction for the BIT/FIT PDF for ${St\geq O(1)}$, and for final/initial separations in the dissipation regime. The prediction shows how caustics in the particle relative velocities in the dissipation range affect the scaling of the pair-separation PDF, leading to a PDF with an algebraically decaying tail. The predicted functional behavior of the PDFs is universal, in that it does not depend upon the level of intermittency in the underlying turbulence. We also analyze the pair-separation PDFs for fluid particles at short-times, and construct theoretical predictions using the multifractal formalism to describe the fluid relative velocity distributions. The theoretical and numerical results both suggest that the extreme events in the inertial particle-pair dispersion at the small-scales are dominated by their non-local interaction with the turbulent velocity field, rather than due to the strong dissipation range intermittency of the turbulence itself...
△ Less
Submitted 14 August, 2017; v1 submitted 29 November, 2016;
originally announced November 2016.
-
Developments and difficulties in predicting the relative velocities of inertial particles at the small-scales of turbulence
Authors:
Andrew D. Bragg
Abstract:
In this paper, we use our recently developed theory for the backward-in-time (BIT) relative dispersion of inertial particles in turbulence (Bragg \emph{et al.}, Phys. Fluids 28, 013305, 2016) to develop the theoretical model by Pan \& Padoan (J. Fluid Mech. 661 73, 2010) for inertial particle relative velocities in isotropic turbulence. We focus on the most difficult regime to model, the dissipati…
▽ More
In this paper, we use our recently developed theory for the backward-in-time (BIT) relative dispersion of inertial particles in turbulence (Bragg \emph{et al.}, Phys. Fluids 28, 013305, 2016) to develop the theoretical model by Pan \& Padoan (J. Fluid Mech. 661 73, 2010) for inertial particle relative velocities in isotropic turbulence. We focus on the most difficult regime to model, the dissipation range, and find that the modified Pan \& Padoan model (that uses the BIT dispersion theory) can lead to significantly improved predictions for the relative velocities, when compared with Direct Numerical Simulation (DNS) data. However, when the particle separation distance, $r$, is less than the Kolmogorov length scale, $η$, the modified model overpredicts the DNS data. We explain how these overpredictions arise from two assumptions in the BIT dispersion theory, that are in general not satisfied when the final separation of the BIT dispersing particles is $<η$. We then demonstrate the failure of both the original and modified versions of the Pan \& Padoan model to predict the correct scale-invariant forms for the inertial particle relative velocity structure functions in the dissipation regime. It is shown how this failure, which is also present in other models, is associated with our present inability to correctly predict not only the quantitative, but also the qualitative behavior of the Radial Distribution Function in the dissipation range when $St=\mathcal{O}(1)$
△ Less
Submitted 13 September, 2016;
originally announced September 2016.
-
Spectral Model of Non-Stationary, Inhomogeneous Turbulence
Authors:
Andrew D. Bragg,
Susan Kurien,
Timothy T. Clark
Abstract:
We compare results from a spectral model for non-stationary, inhomogeneous turbulence (Besnard et al., Theor. Comp. Fluid. Dyn., vol. 8, pp 1-35, 1996) with Direct Numerical Simulation (DNS) data of a shear-free mixing layer (SFML) (Tordella et al., Phys. Rev. E, vol. 77, 016309, 2008). The SFML is used as a test case in which the efficacy of the model closure for the physical-space transport of t…
▽ More
We compare results from a spectral model for non-stationary, inhomogeneous turbulence (Besnard et al., Theor. Comp. Fluid. Dyn., vol. 8, pp 1-35, 1996) with Direct Numerical Simulation (DNS) data of a shear-free mixing layer (SFML) (Tordella et al., Phys. Rev. E, vol. 77, 016309, 2008). The SFML is used as a test case in which the efficacy of the model closure for the physical-space transport of the fluid velocity field can be tested in a flow with inhomogeneity, without the additional complexity of mean-flow coupling. The model is able to capture certain features of the SFML quite well for intermediate to long-times, including the evolution of the mixing-layer width and turbulent kinetic energy. At short-times, and for more sensitive statistics such as the generation of the velocity field anisotropy, the model is less accurate. We present arguments, supported by the DNS data, that a significant cause of the discrepancies is the local approximation to the intrinsically non-local pressure-transport in physical-space that was made in the model, the effects of which would be particularly strong at short-times when the inhomogeneity of the SFML is strongest.
△ Less
Submitted 18 November, 2015;
originally announced November 2015.
-
The effect of Reynolds number on inertial particle dynamics in isotropic turbulence. Part I: Simulations without gravitational effects
Authors:
Peter J. Ireland,
Andrew D. Bragg,
Lance R. Collins
Abstract:
In this study, we analyze the statistics of both individual inertial particles and inertial particle pairs in direct numerical simulations of homogeneous isotropic turbulence in the absence of gravity. The effect of the Taylor microscale Reynolds number $R_λ$ on the particle statistics is examined over the largest range to date (from $R_λ= 88-597$). We first explore the effect of preferential samp…
▽ More
In this study, we analyze the statistics of both individual inertial particles and inertial particle pairs in direct numerical simulations of homogeneous isotropic turbulence in the absence of gravity. The effect of the Taylor microscale Reynolds number $R_λ$ on the particle statistics is examined over the largest range to date (from $R_λ= 88-597$). We first explore the effect of preferential sampling on the single-particle statistics, and use our understanding of preferential sampling to provide a physical explanation for many of the trends in the particle velocity gradients, kinetic energies, and accelerations at low $St$. As $St$ increases, inertial filtering effects become more important, causing the particle kinetic energies and accelerations to decrease.
We then consider particle-pair statistics, and focus our attention on the relative velocities and radial distribution functions (RDFs) of the particles. The relative velocity statistics indicate that preferential-sampling effects are important for $St \lesssim 0.1$ and that path-history/non-local effects become increasingly important for $St \gtrsim 0.2$. The lower-order relative velocity statistics are only weakly sensitive to changes in Reynolds number at low $St$. We find that the RDFs peak near $St$ of order unity, that they exhibit power-law scaling for low and intermediate $St$, and that they are largely independent of Reynolds number for low and intermediate $St$. We also observe that at large $St$, changes in the RDF are related to changes the scaling exponents of the relative velocity variances. The particle collision kernel is found to be largely insensitive to the flow Reynolds number, suggesting that relatively low-Reynolds-number simulations may be able to capture much of the relevant physics of droplet collisions and growth in the adiabatic cores of atmospheric clouds.
△ Less
Submitted 24 July, 2015;
originally announced July 2015.
-
The effect of Reynolds number on inertial particle dynamics in isotropic turbulence. Part II: Simulations with gravitational effects
Authors:
Peter J. Ireland,
Andrew D. Bragg,
Lance R. Collins
Abstract:
In Part I of this study, we analyzed the motion of inertial particles in isotropic turbulence in the absence of gravity using direct numerical simulation (DNS). Here, in Part II, we introduce gravity and study its effect over a wide range of flow Reynolds numbers, Froude numbers, and particle Stokes numbers. We see that gravity causes particles to sample the flow more uniformly and reduces the tim…
▽ More
In Part I of this study, we analyzed the motion of inertial particles in isotropic turbulence in the absence of gravity using direct numerical simulation (DNS). Here, in Part II, we introduce gravity and study its effect over a wide range of flow Reynolds numbers, Froude numbers, and particle Stokes numbers. We see that gravity causes particles to sample the flow more uniformly and reduces the time particles can spend interacting with the underlying turbulence. We also find that gravity tends to increase inertial particle accelerations, and we introduce a model to explain that effect.
We then analyze the particle relative velocities and radial distribution functions (RDFs), which are generally seen to be independent of Reynolds number for low and moderate Kolmogorov-scale Stokes numbers $St$. We see that gravity causes particle relative velocities to decrease, and that the relative velocities have higher scaling exponents with gravity. We observe that gravity has a non-trivial effect on clustering, acting to decrease clustering at low $St$ and to increase clustering at high $St$. By considering the effect of gravity on the clustering mechanisms described in the theory of Zaichik & Alipchenkov (New J. Phys., 11:103018, 2009), we provide an explanation for this non-trivial effect of gravity. We also show that when the effects of gravity are accounted for in the theory of Zaichik & Alipchenkov, the results compare favorably with DNS. The relative velocities and RDFs exhibit considerable anisotropy at small separations, and this anisotropy is quantified using spherical harmonic functions. We use the relative velocities and the RDFs to compute the particle collision kernels, and find that the collision kernel remains as it was for the case without gravity, namely nearly independent of Reynolds number for low and moderate $St$.
△ Less
Submitted 24 July, 2015;
originally announced July 2015.
-
On the relationship between the non-local clustering mechanism and preferential concentration
Authors:
Andrew D. Bragg,
Peter J. Ireland,
Lance R. Collins
Abstract:
`Preferential concentration' (\emph{Phys. Fluids} \textbf{A3}:1169--78, 1991) refers to the clustering of inertial particles in the high-strain, low-rotation regions of turbulence. The `centrifuge mechanism' of Maxey (\emph{J. Fluid Mech.} \textbf{174}:441--65, 1987) appears to explain this phenomenon. In a recent paper, Bragg \& Collins (\emph{New J. Phys.} \textbf{16}:055013, 2014) showed that t…
▽ More
`Preferential concentration' (\emph{Phys. Fluids} \textbf{A3}:1169--78, 1991) refers to the clustering of inertial particles in the high-strain, low-rotation regions of turbulence. The `centrifuge mechanism' of Maxey (\emph{J. Fluid Mech.} \textbf{174}:441--65, 1987) appears to explain this phenomenon. In a recent paper, Bragg \& Collins (\emph{New J. Phys.} \textbf{16}:055013, 2014) showed that the centrifuge mechanism is dominant only in the regime ${St\ll1}$, where $St$ is the Stokes number based on the Kolmogorov time scale. Outside this regime, the centrifuge mechanism gives way to a non-local, path-history symmetry breaking mechanism. However, despite the change in the clustering mechanism, the instantaneous particle positions continue to correlate with high-strain, low-rotation regions of the turbulence. In this paper, we analyze the exact equation governing the radial distribution function and show how the non-local clustering mechanism is influenced by, but not dependent upon, the preferential sampling of the fluid velocity gradient tensor along the particle path-histories in such a way as to generate a bias for clustering in high-strain regions of the turbulence. We also show how the non-local mechanism still generates clustering, but without preferential concentration, in the limit where the timescales of the fluid velocity gradient tensor measured along the inertial particle trajectories approaches zero (such as white-noise flows or for particles in turbulence settling under strong gravity). Finally, we use data from a direct numerical simulation of inertial particles suspended in Navier-Stokes turbulence to validate the arguments presented in this study.
△ Less
Submitted 15 January, 2015;
originally announced January 2015.
-
Mechanisms for the clustering of inertial particles in the inertial range of isotropic turbulence
Authors:
Andrew D. Bragg,
Peter J. Ireland,
Lance R. Collins
Abstract:
In this paper, we consider the physical mechanism for the clustering of inertial particles in the inertial range of isotropic turbulence. We analyze the exact, but unclosed, equation governing the radial distribution function (RDF) and compare the mechanisms it describes for clustering in the dissipation and inertial ranges. We demonstrate that in the limit ${St_r\ll1}$, where $St_r$ is the Stokes…
▽ More
In this paper, we consider the physical mechanism for the clustering of inertial particles in the inertial range of isotropic turbulence. We analyze the exact, but unclosed, equation governing the radial distribution function (RDF) and compare the mechanisms it describes for clustering in the dissipation and inertial ranges. We demonstrate that in the limit ${St_r\ll1}$, where $St_r$ is the Stokes number based on the eddy turnover timescale at separation $r$, the clustering in the inertial range can be understood to be due to the preferential sampling of the coarse-grained fluid velocity gradient tensor at that scale. When ${St_r\gtrsim\mathcal{O}(1)}$ this mechanism gives way to a non-local clustering mechanism. These findings reveal that the clustering mechanisms in the inertial range are analogous to the mechanisms that we identified for the dissipation regime (see \emph{New J. Phys.} \textbf{16}:055013, 2014). Further, we discuss the similarities and differences between the clustering mechanisms we identify in the inertial range and the "sweep-stick" mechanism developed by Coleman \& Vassilicos (Phys. Fluids 21:113301, 2009). We argue that when ${St_r\ll1}$ the sweep-stick mechanism is equivalent to our mechanism in the inertial range if the particles are suspended in Navier-Stokes turbulence, but that the sweep-stick mechanism breaks down for ${St_r\gtrsim\mathcal{O}(1)}$. The argument also explains why the sweep-stick mechanism is unable to predict particle clustering in kinematic simulations. We then consider the closed, model equation for the RDF given in Zaichik \& Alipchenkov (Phys. Fluids. 19:113308, 2007) and use this, together with the results from our analysis, to predict the analytic form of the RDF in the inertial range for ${St_r\ll1}$, which, unlike that in the dissipation range, is not scale-invariant. The results are in good agreement with direct...
△ Less
Submitted 26 November, 2014;
originally announced November 2014.
-
Forward and backward in time dispersion of fluid and inertial particles in isotropic turbulence
Authors:
Andrew D. Bragg,
Peter J. Ireland,
Lance R. Collins
Abstract:
In this paper we investigate both theoretically and numerically the forward in time (FIT) and backward in time (BIT) dispersion of fluid and inertial particle pairs in isotropic turbulence. Fluid particles are known to separate faster BIT than FIT in three-dimensional turbulence, and we find that inertial particles do the same. However, we find that the irreversibility in the inertial particle dis…
▽ More
In this paper we investigate both theoretically and numerically the forward in time (FIT) and backward in time (BIT) dispersion of fluid and inertial particle pairs in isotropic turbulence. Fluid particles are known to separate faster BIT than FIT in three-dimensional turbulence, and we find that inertial particles do the same. However, we find that the irreversibility in the inertial particle dispersion is in general much stronger than that for fluid particles. For example, the ratio of the BIT to FIT mean-square separation can be up to an order of magnitude larger for inertial particles than for the fluid particles. We also find that for both the inertial and fluid particles the irreversibility becomes stronger as the scale of their separation decreases. Regarding the physical mechanism for the irreversibility, we argue that whereas the irreversibility of fluid particle-pair dispersion can be understood in terms of a directional bias arising from the energy transfer process in turbulence, inertial particles experience an additional source of irreversibility arising from the non-local contribution to their velocity dynamics, a contribution which vanishes in the limit ${St\to0}$, where $St$ is the particle Stokes number. For each given initial (final, in the backward in time case) separation $\bm{r}^0$ there is an optimum value of $St$ for which the dispersion irreversibility is strongest, as such particles are optimally affected by both sources of irreversibility. We derive analytical expressions for the BIT, mean-square separation of inertial particles and compare the predictions with...
△ Less
Submitted 22 July, 2015; v1 submitted 21 March, 2014;
originally announced March 2014.