BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Andrew D. Bragg

Andrew D. Bragg appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.flu-dynphysics.comp-ph

Catalog footprint

What is connected

15works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2025arXiv

Mechanism generating reverse buoyancy flux at the small scales of stably stratified turbulence

Previous studies have shown that at the small-scales of stably stratified turbulence, the scale-dependent buoyancy flux reverses sign, such that there is a conversion of turbulent potential energy (TPE) back into turbulent kinetic energy (TKE) at these scales. Moreover, the magnitude of the reverse flux becomes stronger with increasing Prandtl number $Pr$. Using a filtering analysis we demonstrate analytically how this flux reversal is connected to the mechanism identified in Bragg \& de Bruyn Kops (JFM 2024 Vol 991 A10) that is responsible for the surprising observation that the TKE dissipation rate increases while the TPE dissipation rate decreases with increasing $Pr$ in stratified turbulence. The mechanism identified by Bragg \& de Bruyn Kops, which is connected to the formation of ramp-cliff structures in the density field, is shown to give the scale-local contribution to the buoyancy flux. At the smallest-scales this local contribution dominates and explains the flux reversal, while at larger scales a non-local contribution is important. Direct numerical simulations (DNS) of 3D statistically stationary, stably stratified turbulence in the strongly stratified regime confirm the theoretical analysis, and indicate that while on average the local contribution only dominates the buoyancy flux at the smallest scales, it remains strongly correlated with the buoyancy flux at all scales. The results show that ramp-cliffs are not only connected to the reversal of the local buoyancy flux but also the non-local part. At the small scales (approximately below the Ozmidov scale), ramp structures contribute exclusively to reverse buoyancy flux events, whereas cliff structures contribute to both forward and reverse buoyancy flux events.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A stochastic model for the residence time of solid particles in turbulent Rayleigh-Bénard Flow

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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

How two-way coupling modifies the multiscale preferential sweeping mechanism

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An experimental study on the multiscale properties of turbulence in bubble-laden flows

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Analysis of scale-dependent kinetic and potential energy in sheared, stably stratified turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dimension Reduced Turbulent Flow Data From Deep Vector Quantizers

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...

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Is vortex stretching the main cause of the turbulent energy cascade?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Local analysis of the clustering, velocities and accelerations of particles settling in turbulence

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 "..."

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Multiscale fluid--particle thermal interaction in isotropic turbulence

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...

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Forward and backward in time dispersion of fluid and inertial particles in isotropic turbulence

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...

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the relationship between the non-local clustering mechanism and preferential concentration

`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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spectral Model of Non-Stationary, Inhomogeneous Turbulence

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of Reynolds number on inertial particle dynamics in isotropic turbulence. Part I: Simulations without gravitational effects

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effect of Reynolds number on inertial particle dynamics in isotropic turbulence. Part II: Simulations with gravitational effects

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$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mechanisms for the clustering of inertial particles in the inertial range of isotropic turbulence

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...

0 citations0 reviewsNo code linkedOpen access