The dynamics of air bubbles in turbulent Rayleigh–Bénard (RB) convection is described for the first time using laboratory experiments and complementary numerical simulations. We performed experiments at Ra =5.5x10^9 and 1.1x 10^10, where streams of 1 mm bubbles were released at various locations from the bottom of the tank along the path of the roll structure. Using three-dimensional particle tracking velocimetry, we simultaneously tracked a large number of bubbles to inspect the pair dispersion, R2(t), for a range of initial separations, r, spanning one order of magnitude, namely 25η < r < 225η; here η is the local Kolmogorov length scale. Pair dispersion, R2(t), of the bubbles within a quiescent medium was also determined to assess the effect of inhomogeneity and anisotropy induced by the RB convection. Results show that R2(t) underwent a transition phase similar to the ballistic-to-diffusive (t^2-to-t^1) regime in the vicinity of the cell centre; it approached a bulk behavior t3/2 in the diffusive regime as the distance away from the cell centre increased. At small r, R2(t) ~ t^1 is shown in the diffusive regime with a lower magnitude compared to the quiescent case, indicating that the convective turbulence reduced the amplitude of the bubble’s fluctuations. This phenomenon associated to the bubble path instability was further explored by the autocorrelation of the bubble’s horizontal velocity. At large initial separations, R2(t) ∝ t^2 was observed, showing the effect of the roll structure
more »
« less
On the dynamics of air bubbles in Rayleigh–Bénard convection
The dynamics of air bubbles in turbulent Rayleigh–Bénard (RB) convection is described for the first time using laboratory experiments and complementary numerical simulations. We performed experiments at $Ra=5.5\times 10^{9}$ and $1.1\times 10^{10}$ , where streams of 1 mm bubbles were released at various locations from the bottom of the tank along the path of the roll structure. Using three-dimensional particle tracking velocimetry, we simultaneously tracked a large number of bubbles to inspect the pair dispersion, $R^{2}(t)$ , for a range of initial separations, $r$ , spanning one order of magnitude, namely $25\unicode[STIX]{x1D702}\leqslant r\leqslant 225\unicode[STIX]{x1D702}$ ; here $\unicode[STIX]{x1D702}$ is the local Kolmogorov length scale. Pair dispersion, $R^{2}(t)$ , of the bubbles within a quiescent medium was also determined to assess the effect of inhomogeneity and anisotropy induced by the RB convection. Results show that $R^{2}(t)$ underwent a transition phase similar to the ballistic-to-diffusive ( $t^{2}$ -to- $t^{1}$ ) regime in the vicinity of the cell centre; it approached a bulk behavior $t^{3/2}$ in the diffusive regime as the distance away from the cell centre increased. At small $r$ , $R^{2}(t)\propto t^{1}$ is shown in the diffusive regime with a lower magnitude compared to the quiescent case, indicating that the convective turbulence reduced the amplitude of the bubble’s fluctuations. This phenomenon associated to the bubble path instability was further explored by the autocorrelation of the bubble’s horizontal velocity. At large initial separations, $R^{2}(t)\propto t^{2}$ was observed, showing the effect of the roll structure.
more »
« less
- Award ID(s):
- 1912824
- NSF-PAR ID:
- 10157234
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 891
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
For each $t\in \mathbb{R}$ , we define the entire function $$\begin{eqnarray}H_{t}(z):=\int _{0}^{\infty }e^{tu^{2}}\unicode[STIX]{x1D6F7}(u)\cos (zu)\,du,\end{eqnarray}$$ where $\unicode[STIX]{x1D6F7}$ is the super-exponentially decaying function $$\begin{eqnarray}\unicode[STIX]{x1D6F7}(u):=\mathop{\sum }_{n=1}^{\infty }(2\unicode[STIX]{x1D70B}^{2}n^{4}e^{9u}-3\unicode[STIX]{x1D70B}n^{2}e^{5u})\exp (-\unicode[STIX]{x1D70B}n^{2}e^{4u}).\end{eqnarray}$$ Newman showed that there exists a finite constant $\unicode[STIX]{x1D6EC}$ (the de Bruijn–Newman constant ) such that the zeros of $H_{t}$ are all real precisely when $t\geqslant \unicode[STIX]{x1D6EC}$ . The Riemann hypothesis is equivalent to the assertion $\unicode[STIX]{x1D6EC}\leqslant 0$ , and Newman conjectured the complementary bound $\unicode[STIX]{x1D6EC}\geqslant 0$ . In this paper, we establish Newman’s conjecture. The argument proceeds by assuming for contradiction that $\unicode[STIX]{x1D6EC}<0$ and then analyzing the dynamics of zeros of $H_{t}$ (building on the work of Csordas, Smith and Varga) to obtain increasingly strong control on the zeros of $H_{t}$ in the range $\unicode[STIX]{x1D6EC}
more » « less We mimic a flapping wing through a fluid–structure interaction (FSI) framework based upon a generalized lumped-torsional flexibility model. The developed fluid and structural solvers together determine the aerodynamic forces, wing deformation and self-propelled motion. A phenomenological solution to the linear single-spring structural dynamics equation is established to help offer insight and validate the computations under the limit of small deformation. The cruising velocity and power requirements are evaluated by varying the flapping Reynolds number ( $20\leqslant Re_{f}\leqslant 100$ ), stiffness (represented by frequency ratio, $1\lesssim \unicode[STIX]{x1D714}^{\ast }\leqslant 10$ ) and the ratio of aerodynamic to structural inertia forces (represented by a dimensionless parameter $\unicode[STIX]{x1D713}$ ( $0.1\leqslant \unicode[STIX]{x1D713}\leqslant 3$ )). For structural inertia dominated flows ( $\unicode[STIX]{x1D713}\ll 1$ ), pitching and plunging are shown to always remain in phase ( $\unicode[STIX]{x1D719}\approx 0$ ) with the maximum wing deformation occurring at the end of the stroke. When aerodynamics dominates ( $\unicode[STIX]{x1D713}>1$ ), a large phase difference is induced ( $\unicode[STIX]{x1D719}\approx \unicode[STIX]{x03C0}/2$ ) and the maximum deformation occurs at mid-stroke. Lattice Boltzmann simulations show that there is an optimal $\unicode[STIX]{x1D714}^{\ast }$ at which cruising velocity is maximized and the location of optimum shifts away from unit frequency ratio ( $\unicode[STIX]{x1D714}^{\ast }=1$ ) as $\unicode[STIX]{x1D713}$ increases. Furthermore, aerodynamics administered deformations exhibit better performance than those governed by structural inertia, quantified in terms of distance travelled per unit work input. Closer examination reveals that although maximum thrust transpires at unit frequency ratio, it is not transformed into the highest cruising velocity. Rather, the maximum velocity occurs at the condition when the relative tip displacement ${\approx}\,0.3$ .more » « lessLet $f\in C^{2}(\mathbb{T}^{2})$ have mean value 0 and consider $$\begin{eqnarray}\sup _{\unicode[STIX]{x1D6FE}\,\text{closed geodesic}}\frac{1}{|\unicode[STIX]{x1D6FE}|}\biggl|\int _{\unicode[STIX]{x1D6FE}}f\,d{\mathcal{H}}^{1}\biggr|,\end{eqnarray}$$ where $\unicode[STIX]{x1D6FE}$ ranges over all closed geodesics $\unicode[STIX]{x1D6FE}:\mathbb{S}^{1}\rightarrow \mathbb{T}^{2}$ and $|\unicode[STIX]{x1D6FE}|$ denotes its length. We prove that this supremum is always attained. Moreover, we can bound the length of the geodesic $\unicode[STIX]{x1D6FE}$ attaining the supremum in terms of the smoothness of the function: for all $s\geq 2$ , $$\begin{eqnarray}|\unicode[STIX]{x1D6FE}|^{s}{\lesssim}_{s}\biggl(\max _{|\unicode[STIX]{x1D6FC}|=s}\Vert \unicode[STIX]{x2202}_{\unicode[STIX]{x1D6FC}}f\Vert _{L^{1}(\mathbb{T}^{2})}\biggr)\Vert \unicode[STIX]{x1D6FB}f\Vert _{L^{2}}\Vert f\Vert _{L^{2}}^{-2}.\end{eqnarray}$$more » « lessLet $\unicode[STIX]{x1D6FC}\in \mathbb{R}\backslash \mathbb{Q}$ and $\unicode[STIX]{x1D6FD}(\unicode[STIX]{x1D6FC})=\limsup _{n\rightarrow \infty }(\ln q_{n+1})/q_{n}<\infty$ , where $p_{n}/q_{n}$ is the continued fraction approximation to $\unicode[STIX]{x1D6FC}$ . Let $(H_{\unicode[STIX]{x1D706},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D703}}u)(n)=u(n+1)+u(n-1)+2\unicode[STIX]{x1D706}\cos 2\unicode[STIX]{x1D70B}(\unicode[STIX]{x1D703}+n\unicode[STIX]{x1D6FC})u(n)$ be the almost Mathieu operator on $\ell ^{2}(\mathbb{Z})$ , where $\unicode[STIX]{x1D706},\unicode[STIX]{x1D703}\in \mathbb{R}$ . Avila and Jitomirskaya [The ten Martini problem. Ann. of Math. (2), 170 (1) (2009), 303–342] conjectured that, for $2\unicode[STIX]{x1D703}\in \unicode[STIX]{x1D6FC}\mathbb{Z}+\mathbb{Z}$ , $H_{\unicode[STIX]{x1D706},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D703}}$ satisfies Anderson localization if $|\unicode[STIX]{x1D706}|>e^{2\unicode[STIX]{x1D6FD}(\unicode[STIX]{x1D6FC})}$ . In this paper, we develop a method to treat simultaneous frequency and phase resonances and obtain that, for $2\unicode[STIX]{x1D703}\in \unicode[STIX]{x1D6FC}\mathbb{Z}+\mathbb{Z}$ , $H_{\unicode[STIX]{x1D706},\unicode[STIX]{x1D6FC},\unicode[STIX]{x1D703}}$ satisfies Anderson localization if $|\unicode[STIX]{x1D706}|>e^{3\unicode[STIX]{x1D6FD}(\unicode[STIX]{x1D6FC})}$ .more » « less