We present an experimental study on the effects of polymer additives on the turbulent/non-turbulent interface (TNTI) in a fully developed round water jet. The Reynolds number based on the jet diameter is fixed at $Re=7075$. The Weissenberg number $Wi$ ranges from 24 to 86. We employ time-resolved simultaneous particle image velocimetry and laser-induced fluorescence measurements to investigate the local entrainment and engulfment process along the TNTI in two regimes: entrainment transition and enhancement regimes. In polymer-laden jets, the TNTI fluctuates more intermittently in the radial direction and more ambient fluid can be engulfed into the turbulent region due to the augmented large scale motion. Though the contribution of engulfment to the total flux increases with $Wi$, engulfment is still not the major contribution to the entrainment in polymer-laden jets. We further show that the local entrainment velocity is increased in both regimes compared with the pure water jet, due to two contributions: polymer elastic stress and the more intermittent character of the TNTI. In the entrainment transition regime, we observe smaller fractal dimension and shorter length of TNTI compared with the Newtonian case, consistent with previous numerical simulations (Abreu et al. J. Fluid Mech. vol. 934, 2022, A36); whereas those in the enhancement regime remain largely unchanged. The difference between the two regimes results from the fact that the jet flow decays in the streamwise direction. In the entrainment transition regime, turbulence intensity is strong enough to significantly stretch the polymers, resulting in a smoother TNTI in the inertial range. However, in the entrainment enhancement regime, the polymer’s feedback is not strong enough to alter the fractal dimension due to the low elasticity. The above mentioned differences of entrainment velocity and TNTI in the entrainment reduction/transition and enhancement regimes also explain the reduced and enhanced spreading rate of the viscoelastic jet observed in previous numerical simulations and experiments (Guimarães et al. J. Fluid Mech. 2020,vol. 899, A11; Peng et al. Phys. Fluids, 2023, vol. 35, 045110).

Recent experiments and simulations have sparked growing interest in the study of Rayleigh–Bénard convection in very slender cells. One pivotal inquiry arising from this interest is the elucidation of the flow structure within these very slender cells. Here we employ tomographic particle image velocimetry, for the first time, to capture experimentally the full-field three-dimensional and three-component velocity field in a very slender cylindrical cell with aspect ratio $\Gamma =1/10$. The experiments cover a Rayleigh number range $5.0 \times 10^8 \leqslant Ra \leqslant 5.0 \times 10^9$ and Prandtl number 5.7. Our experiments reveal that the flow structure in the $\Gamma =1/10$ cell is neither in the multiple-roll form nor in the simple helical form; instead, the ascending and descending flows can intersect and cross each other, resulting in the crossing events. These crossing events separate the flow into segments; within each segment, the ascending and descending flows ascend or descend side by side vertically or in the twisting manner, and the twisting is not unidirectional, while the segments near the boundary can also be in the form of a donut like structure. By applying the mode decomposition analyses to the measured three-dimensional velocity fields, we identified the crossing events as well as the twisting events for each instantaneous flow field. Statistical analysis of the modes reveals that as $Ra$ increases, the average length of the segments becomes smaller, and the average number of segments increases from 2.5 to 3.9 in the $Ra$ range of our experiments.

We studied flow organization and heat transfer properties in mixed turbulent convection within Poiseuille–Rayleigh–Bénard channels subjected to temporally modulated sinusoidal wall temperatures. Three-dimensional direct numerical simulations were performed for Rayleigh numbers in the range $10^6 \leqslant Ra \leqslant 10^8$, a Prandtl number $Pr = 0.71$ and a bulk Reynolds number $Re_b \approx 5623$. We found that high-frequency wall temperature oscillations had minimal impact on flow structures, while low-frequency oscillations induced adaptive changes, forming stable stratified layers during cooling. Proper orthogonal decomposition (POD) analysis revealed a dominant streamwise unidirectional shear flow mode. Large-scale rolls oriented in the streamwise direction appeared as higher POD modes and were significantly influenced by lower-frequency wall temperature variations. Long-time-averaged statistics showed that the Nusselt number increased with decreasing frequency by up to 96 %, while the friction coefficient varied by less than 15 %. High-frequency modulation predominantly influenced near-wall regions, enhancing convective effects, whereas low frequencies reduced these effects via stable stratified layer formation. Phase-averaged statistics showed that high-frequency modulation resulted in phase-stable streamwise velocity and temperature profiles, while low-frequency modulation caused significant variations due to weakened turbulence. Turbulent kinetic energy (TKE) profiles remained high near the wall during both heating and cooling at high frequency, but decreased during cooling at low frequencies. A TKE budget analysis revealed that during heating, TKE production was dominated by shear near the wall and by buoyancy in the bulk region; while during cooling, the production, distribution and dissipation of TKE were all nearly zero.

We present an experimental study on the drag reduction by polymers in Taylor–Couette turbulence at Reynolds numbers ($Re$) ranging from $4\times 10^3$ to $2.5\times 10^4$. In this $Re$ regime, the Taylor vortex is present and accounts for more than 50 % of the total angular velocity flux. Polyacrylamide polymers with two different average molecular weights are used. It is found that the drag reduction rate increases with polymer concentration and approaches the maximum drag reduction (MDR) limit. At MDR, the friction factor follows the $-0.58$ scaling, i.e. $C_f \sim Re^{-0.58}$, similar to channel/pipe flows. However, the drag reduction rate is about $20\,\%$ at MDR, which is much lower than that in channel/pipe flows at comparable $Re$. We also find that the Reynolds shear stress does not vanish and the slope of the mean azimuthal velocity profile in the logarithmic layer remains unchanged at MDR. These behaviours are reminiscent of the low drag reduction regime reported in channel flow (Warholic et al., Exp. Fluids, vol. 27, no. 5, 1999, pp. 461–472). We reveal that the lower drag reduction rate originates from the fact that polymers strongly suppress the turbulent flow while only slightly weaken the mean Taylor vortex. We further show that polymers steady the velocity boundary layer and suppress the small-scale Görtler vortices in the near-wall region. The former effect reduces the emission rate of both intense fast and slow plumes detached from the boundary layer, resulting in less flux transport from the inner cylinder to the outer one and reduces energy input into the bulk turbulent flow. Our results suggest that in turbulent flows, where secondary flow structures are statistically persistent and dominate the global transport properties of the system, the drag reduction efficiency of polymer additives is significantly diminished.

We studied the transport and deposition behaviour of point particles in Rayleigh–Bénard convection cells subjected to Couette-type wall shear. Direct numerical simulations (DNSs) are performed for Rayleigh number ($Ra$) in the range $10^{7} \leq Ra \leq ~10^9$ with a fixed Prandtl number $Pr = 0.71$, while the wall-shear Reynolds number ($Re_w$) is in the range $0 \leq Re_w \leq ~12\,000$. With the increase of $Re_w$, the large-scale rolls expanded horizontally, evolving into zonal flow in two-dimensional simulations or streamwise-oriented rolls in three-dimensional simulations. We observed that, for particles with a small Stokes number ($St$), they either circulated within the large-scale rolls when buoyancy dominated or drifted near the walls when shear dominated. For medium $St$ particles, pronounced spatial inhomogeneity and preferential concentration were observed regardless of the prevailing flow state. For large $St$ particles, the turbulent flow structure had a minor influence on the particles’ motion; although clustering still occurred, wall shear had a negligible influence compared with that for medium $St$ particles. We then presented the settling curves to quantify the particle deposition ratio on the walls. Our DNS results aligned well with previous theoretical predictions, which state that small $St$ particles settle with an exponential deposition ratio and large $St$ particles settle with a linear deposition ratio. For medium $St$ particles, where complex particle–turbulence interaction emerges, we developed a new model describing the settling process with an initial linear stage followed by a nonlinear stage. Unknown parameters in our model can be determined either by fitting the settling curves or using empirical relations. Compared with DNS results, our model also accurately predicts the average residence time across a wide range of $St$ for various $Re_w$.

To date, a comprehensive understanding of the influence of the Prandtl number ($Pr$) on flow topology in turbulent Rayleigh–Bénard convection (RBC) remains elusive. In this study, we present an experimental investigation into the evolution of flow topology in quasi-two-dimensional turbulent RBC with $7.0 \leq Pr \leq 244.2$ and $2.03\times 10^{8} \leq Ra \leq 2.81\times 10^{9}$. Particle image velocimetry (PIV) measurements reveal the flow transitions from multiple-roll state to single-roll state with increasing $Ra$, and the transition is hindered with increasing $Pr$, i.e. the transitional Rayleigh number $Ra_t$ increases with $Pr$. We mapped out a phase diagram on the flow topology change on $Ra$ and $Pr$, and identified the scaling of $Ra_t$ on $Pr$: $Ra_t \sim Pr^{0.93}$ in the low $Pr$ range, and $Ra_t \sim Pr^{3.3}$ in the high $Pr$ range. The scaling in the low $Pr$ range is consistent with the model of balance of energy dissipation time and plume travel time that we proposed in our previous study, while the scaling in the high $Pr$ range implies a new governing mechanism. For the first time, the scaling of $Re$ on $Ra$ and $Pr$ is acquired through full-field PIV velocity measurement, $Re \sim Ra^{0.63}\,Pr^{-0.87}$. We also propose that increasing horizontal velocity promotes the formation of the large-scale circulation (LSC), especially for the high $Pr$ case. Our proposal was verified by achieving LSC through introducing horizontal driving force $Ra_H$ by tilting the convection cell with a small angle.

We report an experimental study about the effect of an obstructed centre on heat transport and flow reversal by inserting an adiabatic cylinder at the centre of a quasi-two-dimensional Rayleigh–Bénard convection cell. The experiments are carried out in a Rayleigh number ($Ra$) range of $2\times 10^7 \leq Ra \leq 2\times 10^9$ and at a Prandtl number ($Pr$) of $5.7$. It is found that for low $Ra$, the obstructed centre leads to a heat transfer enhancement of up to 21 $\%$, while as $Ra$ increases, the magnitude of the heat transfer enhancement decreases and the heat transfer efficiency ($Nu$) eventually converges to that of the unobstructed normal cell. Particle image velocimetry measurements show that the heat transfer enhancement originates from the change in flow topology due to the presence of the cylindrical obstruction. In the low-$Ra$ regime the presence of the obstruction promotes the transition of the flow topology from the four-roll state to the abnormal single-roll state then to the normal single-roll state with increasing obstruction size. While in the high-$Ra$ regime, the flow is always in the single-roll state regardless of the obstruction size, although the flow becomes more coherent with the size of the obstruction. We also found that in the presence of the cylindrical obstruction, the stability of the corner vortices is significantly reduced, leading to a large reduction in the frequency of flow reversals.

We studied the flow organization and heat transfer properties in two-dimensional and three-dimensional Rayleigh–Bénard cells that are imposed with different types of wall shear. The external wall shear is added with the motivation of manipulating flow mode to control heat transfer efficiency. We imposed three types of wall shear that may facilitate the single-roll, the horizontally stacked double-roll, and the vertically stacked double-roll flow modes, respectively. Direct numerical simulations are performed for fixed Rayleigh number $Ra = 10^{8}$ and fixed Prandtl number $Pr = 5.3$, while the wall-shear Reynolds number ($Re_{w}$) is in the range $60 \leqslant Re_{w} \leqslant 6000$. Generally, we found enhanced heat transfer efficiency and global flow strength with the increase of $Re_{w}$. However, even with the same magnitude of global flow strength, the heat transfer efficiency varies significantly when the cells are under different types of wall shear. An interesting finding is that by increasing the wall-shear strength, the thermal turbulence is relaminarized, and more surprisingly, the heat transfer efficiency in the laminar state is higher than that in the turbulent state. We found that the enhanced heat transfer efficiency at the laminar regime is due to the formation of more stable and stronger convection channels. We propose that the origin of thermal turbulence laminarization is the reduced amount of thermal plumes. Because plumes are mainly responsible for turbulent kinetic energy production, when the detached plumes are swept away by the wall shear, the reduced number of plumes leads to weaker turbulent kinetic energy production. We also quantify the efficiency of facilitating heat transport via external shearing, and find that for larger $Re_{w}$, the enhanced heat transfer efficiency comes at a price of a larger expenditure of mechanical energy.

We present an experimental study of the large-scale vortex (or large-scale circulation, LSC) in turbulent Rayleigh–Bénard convection in a $\varGamma =\text {diameter}/\text {height}=2$ cylindrical cell. The working fluid is deionized water with Prandtl number ($Pr$) around 5.7, and the Rayleigh number ($Ra$) ranges from $7.64\times 10^7$ to $6.06\times 10^8$. We measured the velocity field in various vertical cross-sectional planes by using the planar particle image velocimetry technique. The velocity measurement in the LSC central plane shows that the flow is in the single-roll form, and the centre of the single-roll (vortex) does not always stay at the centre of the cell; instead, it orbits periodically in the direction opposite to the flow direction of the LSC, with its trajectory in the shape of an ellipse. The velocity measurements in the three vertical planes in parallel to the LSC central plane indicate that the flow is in the vortex tube form horizontally filling almost the whole cell, and the centre line of the vortex tube is consistent with the so-called ‘jump rope’ form proposed by a previous study that combined numerical simulation and local velocity measurements in the low $Pr$ case (Vogt et al., Proc. Natl Acad. Sci. USA, vol. 115, 2018, pp. 12674–12679). In addition, we found that the oscillation of the local velocity in $\varGamma =2$ cells originates from the periodical orbiting of the vortex centre. Our velocity measurements further indicate that the vortex centre orbiting is absent in $\varGamma =1$ cells, at least in the $Ra$ range of our experiments.

In this paper, we designed two different configurations with locally isothermal sidewalls, where the temperature is set to be the bulk temperature, to control the large-scale circulation in turbulent Rayleigh–Bénard convection, namely two-point control and four-point control. At fixed Rayleigh number $Ra=10^8$ and Prandtl number $Pr=2$, a series of direct numerical simulations are performed on both two-dimensional (2-D) and quasi-two-dimensional (quasi-2-D) cavities with both types of control, where the width of the control area is fixed at $\delta _c=0.05$ and the vertical distance from the cavity centre $h_c$ varies from 0 to 0.45 with an interval of 0.05. Our results show that the control effect depends on $h_c$, the control configurations as well as the flow dimensions. For 2-D cavities, both two-point control and four-point control suppress the flow reversal when $h_c \geq 0.05$, accompanied by the enhancement of vertical heat transfer and the strength of the large-scale circulation. For quasi-2-D cavities, the suppression of the flow reversals is obvious with two-point control and $h_c\geq 0.05$, while the effect is rather limited with four-point control. Further experiments with $Pr=5.7$ and $Ra$ up to $7.36\times10^8$ show that two-point control with $h_c=0.15$ can effectively suppress the flow reversal, while two-point control with $h_c=0$ can suppress the reversals at low $Ra=1.93\times 10^8$ and activate them at higher $Ra=7.36\times 10^8$, which agrees well with our numerical simulations.

We present a numerical study of the flow states and reversals of the large-scale circulation (LSC) in a two-dimensional circular Rayleigh–Bénard cell. Long-time direct numerical simulations are carried out in the Rayleigh number ($Ra$) range $10^{7} \leqslant Ra \leqslant 10^{8}$ and Prandtl number ($Pr$) range $2.0 \leqslant Pr \leqslant 20.0$. We found that a new, long-lived, chaotic flow state exists, in addition to the commonly observed circulation states (the LSC in the clockwise and counterclockwise directions). The circulation states consist of one primary roll in the middle and two secondary rolls near the top and bottom circular walls. The primary roll becomes stronger and larger, while the two secondary rolls diminish, with increasing $Ra$. Our results suggest that the reversal of the LSC is accompanied by the secondary rolls growing, breaking the primary roll and then connecting to form a new primary roll with reversed direction. We mapped out the phase diagram of the existence of the LSC and the reversal in the $Ra\text{--}Pr$ space, which reveals that the flow is in the circulation states when $Ra$ is large and $Pr$ is small. The reversal of the LSC can only occur in a limited $Pr$ range. The phase diagram can be understood in terms of competition between the thermal and viscous diffusions. We also found that the internal flow states manifested themselves into global properties such as Nusselt and Reynolds numbers.

We present an experimental study on controlling the number of vortices and the torque in a Taylor–Couette flow of water for Reynolds numbers from 660 to 1320. Different flow states are achieved in the annulus of width $d$ between the inner rotating and outer stationary cylinders through manipulating the initial height of the water annulus. We show that the torque exerted on the inner cylinder of the Taylor–Couette system can be reduced by up to 20 % by controlling the flow at a state where fewer than the nominal number of vortices develop between the cylinders. This flow state is achieved by starting the system with an initial water annulus height $h_0$ (which nominally corresponds to $h_0/d$ vortices), then gradually adding water into the annulus while the inner cylinder keeps rotating. During this filling process the flow topology is so persistent that the number of vortices does not increase; instead, the vortices are greatly stretched in the axial (vertical) direction. We show that this state with stretched vortices is sustainable until the vortices are stretched to around 2.05 times their nominal size. Our experiments reveal that by manipulating the initial height of the liquid annulus we are able to generate different flow states and demonstrate how the different flow states manifest themselves in global momentum transport.

Blockchain technology has rapidly emerged as a decentralized trusted network to replace the traditional centralized intermediator. Especially, the smart contracts that are based on blockchain allow users to define the agreed behaviour among them, the execution of which will be enforced by the smart contracts. Based on this, we propose a decentralized booking system that uses the blockchain as the intermediator between hoteliers and travellers. The system enjoys the trustworthiness of blockchain, improves efficiency and reduces the cost of the traditional booking agencies. The design of the system has been formally modelled using the CSP# language and verified using the model checker Process Analysis Toolkit. We have implemented a prototype decentralized booking system based on the Ethereum ecosystem.

We report a comparative experimental study of the reversal of the large-scale circulation in turbulent Rayleigh–Bénard convection in a quasi-two-dimensional corner-less cell where the corner vortices are absent and in a quasi-two-dimensional normal cell where the corner vortices are present. It is found that in the corner-less cell the reversal frequency exhibits a slow decrease followed by a fast decrease with increasing Rayleigh number $Ra$, separated by a transitional $Ra$ ($Ra_{t,r}$). The transition is similar to that in the normal cell, and $Ra_{t,r}$ is almost the same for both cells. Despite the similarities, the reversal frequency is greatly reduced in the corner-less cell. The reduction of the reversal frequency is more significant, in terms of both the amplitude and the scaling exponent, in the high-$Ra$ regime. In addition, we classified the reversals into main-vortex-led and corner-vortex-led, and found that both types exist in the normal cell while only the former exists in the corner-less cell. The frequency of main-vortex-led reversal in the normal cell is found to be in excellent agreement with the frequency of reversals in the corner-less cell. Our results reveal for the first time the quantitative role of the corner vortices in the occurrence of the reversals of the large-scale circulation.

The head-on collision between two dust-acoustic solitary waves in a non-magnetized, collisionless and strongly coupled dust plasma has been studied. The application scope of the analytical solution of the head-on collision is given in the present paper by using the particle-in-cell simulation method. It is noted that the analytical results are valid if the amplitudes of both of the solitary waves are small enough. The effects of the coupling parameters on both the head-on collision and the waveform are also studied in the present paper.