An official website of the United States government
Based on a generalized local Kolmogorov–Hill equation expressing the evolution of kinetic energy integrated over spheres of size$$\ell$$in the inertial range of fluid turbulence, we examine a possible definition of entropy and entropy generation for turbulence. Its measurement from direct numerical simulations in isotropic turbulence leads to confirmation of the validity of the fluctuation relation (FR) from non-equilibrium thermodynamics in the inertial range of turbulent flows. Specifically, the ratio of probability densities of forward and inverse cascade at scale$$\ell$$is shown to follow exponential behaviour with the entropy generation rate if the latter is defined by including an appropriately defined notion of ‘temperature of turbulence’ proportional to the kinetic energy at scale$$\ell$$.
more »
« less
This content will become publicly available on November 25, 2025
Direct numerical simulation of bubble rising in turbulence
We describe the rising trajectory of bubbles in isotropic turbulence and quantify the slowdown of the mean rise velocity of bubbles with sizes within the inertial subrange. We perform direct numerical simulations of bubbles, for a wide range of turbulence intensity, bubble inertia and deformability, with systematic comparison with the corresponding quiescent case, with Reynolds number at the Taylor microscale from 38 to 77. Turbulent fluctuations randomise the rising trajectory and cause a reduction of the mean rise velocity$$\tilde {w}_b$$compared with the rise velocity in quiescent flow$$w_b$$. The decrease in mean rise velocity of bubbles$$\tilde {w}_b/w_b$$is shown to be primarily a function of the ratio of the turbulence intensity and the buoyancy forces, described by the Froude number$$Fr=u'/\sqrt {gd}$$, where$$u'$$is the root-mean-square velocity fluctuations,$$g$$is gravity and$$d$$is the bubble diameter. The bubble inertia, characterised by the ratio of inertial to viscous forces (Galileo number), and the bubble deformability, characterised by the ratio of buoyancy forces to surface tension (Bond number), modulate the rise trajectory and velocity in quiescent fluid. The slowdown of these bubbles in the inertial subrange is not due to preferential sampling, as is the case with sub-Kolmogorov bubbles. Instead, it is caused by the nonlinear drag–velocity relationship, where velocity fluctuations lead to an increased average drag. For$$Fr > 0.5$$, we confirm the scaling$$\tilde {w}_b / w_b \propto 1 / Fr$$, as proposed previously by Ruthet al.(J. Fluid Mech., vol. 924, 2021, p. A2), over a wide range of bubble inertia and deformability.
more »
« less
- PAR ID:
- 10581919
- Publisher / Repository:
- Cambridge University Press
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 999
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
An analysis is presented for the axisymmetric lubrication resistance between permeable spherical particles. Darcy's law is used to describe the flow in the permeable medium and a slip boundary condition is applied at the interface. The pressure in the near-contact region is governed by a non-local integral equation. The asymptotic limit$$K=k/a^{2} \ll 1$$is considered, where$$k$$is the arithmetic mean permeability, and$$a^{-1}=a^{-1}_{1}+a^{-1}_{2}$$is the reduced radius, and$$a_1$$and$$a_2$$are the particle radii. The formulation allows for particles with distinct particle radii, permeabilities and slip coefficients, including permeable and impermeable particles and spherical drops. Non-zero particle permeability qualitatively affects the axisymmetric near-contact motion, removing the classical lubrication singularity for impermeable particles, resulting in finite contact times under the action of a constant force. The lubrication resistance becomes independent of gap and attains a maximum value at contact$$F=6{\rm \pi} \mu a W K^{-2/5}\tilde {f}_c$$, where$$\mu$$is the fluid viscosity,$$W$$is the relative velocity and$$\tilde {f}_c$$depends on slip coefficients and weakly on permeabilities; for two permeable particles with no-slip boundary conditions,$$\tilde {f}_c=0.7507$$; for a permeable particle in near contact with a spherical drop,$$\tilde {f}_c$$is reduced by a factor of$$2^{-6/5}$$.more » « less
-
Abstract We give an explicit raising operator formula for the modified Macdonald polynomials$$\tilde {H}_{\mu }(X;q,t)$$, which follows from our recent formula for$$\nabla $$on an LLT polynomial and the Haglund-Haiman-Loehr formula expressing modified Macdonald polynomials as sums of LLT polynomials. Our method just as easily yields a formula for a family of symmetric functions$$\tilde {H}^{1,n}(X;q,t)$$that we call$$1,n$$-Macdonald polynomials, which reduce to a scalar multiple of$$\tilde {H}_{\mu }(X;q,t)$$when$$n=1$$. We conjecture that the coefficients of$$1,n$$-Macdonald polynomials in terms of Schur functions belong to$${\mathbb N}[q,t]$$, generalizing Macdonald positivity.more » « less
-
Free surface flows driven by boundary undulations are observed in many biological phenomena, including the feeding and locomotion of water snails. To simulate the feeding strategy of apple snails, we develop a centimetric robotic undulator that drives a thin viscous film of liquid with the wave speed$$V_w$$. Our experimental results demonstrate that the behaviour of the net fluid flux$$Q$$strongly depends on the Reynolds number$$Re$$. Specifically, in the limit of vanishing$$Re$$, we observe that$$Q$$varies non-monotonically with$$V_w$$, which has been successfully rationalised by Pandeyet al.(Nat. Commun., vol. 14, no. 1, 2023, p. 7735) with the lubrication model. By contrast, in the regime of finite inertia ($${Re} \sim O(1)$$), the fluid flux continues to increase with$$V_w$$and completely deviates from the prediction of lubrication theory. To explain the inertia-enhanced pumping rate, we build a thin-film, two-dimensional model via the asymptotic expansion in which we linearise the effects of inertia. Our model results match the experimental data with no fitting parameters and also show the connection to the corresponding free surface shapes$$h_2$$. Going beyond the experimental data, we derive analytical expressions of$$Q$$and$$h_2$$, which allow us to decouple the effects of inertia, gravity, viscosity and surface tension on free surface pumping over a wide range of parameter space.more » « less
-
We explored the settling dynamics of vertically aligned particles in a quiescent, stratified two-layer fluid using particle tracking velocimetry. Glass spheres of$$d=4\,{\rm mm}$$diameter were released at frequencies of 4, 6 and 8 Hz near the free surface, traversing through an upper ethanol layer ($$H_1$$), whereHis height or layer thickess, varying from$$10d$$to$$40d$$and a lower oil layer. Results reveal pronounced lateral particle motion in the ethanol layer, attributed to a higher Galileo number ($$Ga = 976$$, ratio of buoyancy–gravity to viscous effects), compared with the less active behaviour in the oil layer ($$Ga = 16$$). The ensemble vertical velocity of particles exhibited a minimum just past the density interface, becoming more pronounced with increasing$$H_1$$, and suggesting that enhanced entrainment from ethanol to oil resulted in an additional buoyancy force. This produced distinct patterns of particle acceleration near the density interface, which were marked by significant deceleration, indicating substantial resistance to particle motion. An increased drag coefficient occurred for$$H_1/d = 40$$compared with a single particle settling in oil; drag reduced as the particle-release frequency ($$\,f_p$$) increased, likely due to enhanced particle interactions at closer proximity. Particle pair dispersions, lateral ($$R^2_L$$) and vertical ($$R^2_z$$), were modulated by$$H_1$$, initial separation$$r_0$$and$$f_p$$. The$$R^2_L$$dispersion displayed ballistic scaling initially, Taylor scaling for$$r_0 < H_1$$and Richardson scaling for$$r_0 > H_1$$. In contrast,$$R^2_z$$followed a$$R^2_z \sim t^{5.5}$$scaling under$$r_0 < H_1$$. Both$$R^2_L$$and$$R^2_z$$plateaued at a distance from the interface, depending on$$H_1$$and$$f_p$$.more » « less
