The seminal Bolgiano–Obukhov (BO) theory established the fundamental framework for turbulent mixing and energy transfer in stably stratified fluids. However, the presence of BO scalings remains debatable despite their being observed in stably stratified atmospheric layers and convective turbulence. In this study, we performed precise temperature measurements with 51 high-resolution loggers above the seafloor for 46 h on the continental shelf of the northern South China Sea. The temperature observation exhibits three layers with increasing distance from the seafloor: the bottom mixed layer (BML), the mixing zone and the internal wave zone. A BO-like scaling $\alpha =-1.34\pm 0.10$ is observed in the temperature spectrum when the BML is in a weakly stable stratified ($N\sim 0.0018$ rad s$^{-1}$) and strongly sheared ($Ri\sim 0.0027$) condition, whereas in the unstably stratified convective turbulence of the BML, the scaling $\alpha =-1.76\pm 0.10$ clearly deviated from the BO theory but approached the classical $-$5/3 scaling in isotropic turbulence. This suggests that the convective turbulence is not the promise of BO scaling. In the mixing zone, where internal waves alternately interact with the BML, the scaling follows the Kolmogorov scaling. In the internal wave zone, the scaling $\alpha =-2.12 \pm 0.15$ is observed in the turbulence range and possible mechanisms are provided.

We report the unified constitutive law of vibroconvective turbulence in microgravity, i.e. $Nu \sim a^{-1} Re_{os}^\beta$ where the Nusselt number $Nu$ measures the global heat transport, $a$ is the dimensionless vibration amplitude, $Re_{os}$ is the oscillational Reynolds number and $\beta$ is the universal exponent. We find that the dynamics of boundary layers plays an essential role in vibroconvective heat transport and the $Nu$-scaling exponent $\beta$ is determined by the competition between the thermal boundary layer (TBL) and vibration-induced oscillating boundary layer (OBL). Then a physical model is proposed to explain the change of scaling exponent from $\beta =2$ in the TBL-dominant regime to $\beta = 4/3$ in the OBL-dominant regime. Our finding elucidates the emergence of universal constitutive laws in vibroconvective turbulence, and opens up a new avenue for generating a controllable effective heat transport under microgravity or even microfluidic environment in which the gravity effect is nearly absent.

Soft drink consumption has become a highly controversial public health issue. Given the pattern of consumption in China, sugar-sweetened beverage is the main type of soft drink consumed. Due to containing high levels of fructose, a soft drink may have a deleterious effect on handgrip strength (HGS) due to oxidative stress, inflammation and insulin resistance. However, few studies show an association between soft drink consumption and HGS in adults. We aimed to investigate the association between soft drink consumption and longitudinal changes in HGS among a Chinese adult population. A longitudinal population-based cohort study (5-year follow-up, median: 3·66 years) was conducted in Tianjin, China. A total of 11 125 participants (56·7 % men) were enrolled. HGS was measured using a handheld digital dynamometer. Soft drink consumption (mainly sugar-containing carbonated beverages) was measured at baseline using a validated FFQ. ANCOVA was used to evaluate the association between soft drink consumption and annual change in HGS or weight-adjusted HGS. After adjusting for multiple confounding factors, the least square means (95 % CI) of annual change in HGS across soft drink consumption frequencies were −0·70 (–2·49, 1·09) for rarely drinks, −0·82 (–2·62, 0·97) for < 1 cup/week and −0·86 (–2·66, 0·93) for ≥ 1 cup/week (Pfor trend < 0·05). Likewise, a similar association was observed between soft drink consumption and annual change in weight-adjusted HGS. The results indicate that higher soft drink consumption was associated with faster HGS decline in Chinese adults.

The numerical investigation focuses on the flow patterns around a rectangular cylinder with three aspect ratios ($L/D=5$, $10$, $15$) at a Reynolds number of $1000$. The study delves into the dynamics of vortices, their associated frequencies, the evolution of the boundary layer and the decay of the wake. Kelvin–Helmholtz (KH) vortices originate from the leading edge (LE) shear layer and transform into hairpin vortices. Specifically, at $L/D=5$, three KH vortices merge into a single LE vortex. However, at $L/D=10$ and $15$, two KH vortices combine to form a LE vortex, with the rapid formation of hairpin vortex packets. A fractional harmonic arises due to feedback from the split LE shear layer moving upstream, triggering interaction with the reverse flow. Trailing edge (TE) vortices shed, creating a Kármán-like street in the wake. The intensity of wake oscillation at $L/D=5$ surpasses that in the other two cases. Boundary layer transition occurs after the saturation of disturbance energy for $L/D=10$ and $15$, but not for $L/D=5$. The low-frequency disturbances are selected to generate streaks inside the boundary layer. The TE vortex shedding induces the formation of a favourable pressure gradient, accelerating the flow and fostering boundary layer relaminarization. The self-similarity of the velocity defect is observed in all three wakes, accompanied by the decay of disturbance energy. Importantly, the decrease in the shedding frequency of LE (TE) vortices significantly contributes to the overall decay of disturbance energy. This comprehensive exploration provides insights into complex flow phenomena and their underlying dynamics.

This study investigates the effect of vibration on the flow structure transitions in thermal vibrational convection (TVC) systems, which occur when a fluid layer with a temperature gradient is excited by vibration. Direct numerical simulation (DNS) of TVC in a two-dimensional enclosed square box is performed over a range of dimensionless vibration amplitudes $0.001 \le a \le 0.3$ and angular frequencies $10^{2} \le \omega \le 10^{7}$, with a fixed Prandtl number of 4.38. The flow visualisation shows the transition behaviour of flow structure upon the varying frequency, characterising three distinct regimes, which are the periodic-circulation regime, columnar regime and columnar-broken regime. Different statistical properties are distinguished from the temperature and velocity fluctuations at the boundary layer and mid-height. Upon transition into the columnar regime, columnar thermal coherent structures are formed, in contrast to the periodic oscillating circulation. These columns are contributed by the merging of thermal plumes near the boundary layer, and the resultant thermal updrafts remain at almost fixed lateral position, leading to a decrease in fluctuations. We further find that the critical point of this transition can be described nicely by the vibrational Rayleigh number ${{Ra}}_{vib}$. As the frequency continues to increase, entering the so-called columnar-broken regime, the columnar structures are broken, and eventually the flow state becomes a large-scale circulation (LSC), characterised by a sudden increase in fluctuations. Finally, a phase diagram is constructed to summarise the flow structure transition over a wide range of vibration amplitude and frequency parameters.

Intertemporal choices involve tradeoffs between outcomes that occur at different times. Most of the research has used pure gains tasks and the discount rates yielding from those tasks to explain and predict real-world behaviors and consequences. However, real decisions are often more complex and involve mixed outcomes (e.g., sooner-gain and later-loss or sooner-loss and later-gain). No study has used mixed gain-loss intertemporal tradeoff tasks to explain and predict real-world behaviors and consequences, and studies involving such tasks are also scarce. Considering that tasks involving a combination of gains and losses may yield different discount rates and that existing pure gains tasks do not explain or predict real-world outcomes well, this study conducted two experiments to compare the discount rates of mixed gain-loss intertemporal tradeoffs with those of pure gains or pure losses (Experiment 1) and to examine whether these tasks predicted different real-world behaviors and consequences (Experiment 2). Experiment 1 suggests that the discount rate ordering of the four tasks was, from highest to lowest, pure gains, sooner-loss and later-gain, pure losses, and sooner-gain and later-loss. Experiment 2 indicates that the evidence supporting the claim that the discount rates of the four tasks were related to different real-world behaviors and consequences was insufficient.

We investigate the energy transport and heat transfer efficiency in turbulent Rayleigh–Bénard (RB) convection laden with radiatively heated inertial particles. Direct numerical simulations combined with the Lagrangian point-particle mode were carried out in the range of density ratio $854.7\le \rho _p/\rho _0 \le 8547$ and radiation intensity $1\le \phi /\phi _{solar}\le 100$ for both two-dimensional (2-D) and three-dimensional (3-D) simulations. The Rayleigh number ranges from $2\times 10^6$ to $10^8$ for 2-D cases, and is $10^7$ for 3-D cases for $Pr=0.71$. It is found that particles with small density ratio that encounter strong radiation significantly alter the flow momentum transport and fluid heat transfer, so the fluid temperature of bulk is remarkably heated. We then derived the theoretical relation of the Nusselt number for interphase heat transfer in the heated particle-laden RB convection, which reveals that the heat transfer difference between the top and bottom plates stems from the interphase heat transfer. We further found that both the interphase heat transfer and the interphase thermal energy transport exhibit universal properties. They are both increased linearly with the reciprocal of the normalized density ratio. Additionally, both the interphase heat transfer and the interphase thermal energy transport increase linearly with the increase of radiation intensity. The growth rates exhibit specific scaling relations versus Rayleigh number and density ratio. Two different regimes distinguished by the critical density ratio, i.e. the exothermic particle regime and the endothermic particle regime, are observed. We further derived the power-law relation of the critical density ratios versus Rayleigh number and radiation intensity, i.e. $\rho _p/\rho _c \sim (\phi /\phi _{solar})^{1/2}\,Ra^{1/3}$, which is in remarkable agreement with the 3-D simulations.

The flow past a cylinder in proximity to a plane wall is investigated numerically for small gap ratios. Three vortex dynamic processes associated with different hairpin vortex generation mechanisms are identified for the first time, and the wake-induced turbulent transition is analysed. The vortex shedding is suppressed at $G/D = 0.1$, while the spanwise vortex is generated via a Kelvin–Helmholtz instability and evolves into hairpin vortices. For $G/D= 0.3$, the upper and lower rollers alternatively shedding from the cylinder, interact with the secondary vortex. The split secondary vortex merges with the upper roller and results in a new vortex downstream, which develops into hairpin vortices. When $G/D = 0.9$, the secondary vortex interacts with the lower roller and then evolves into hairpin vortices. A tertiary vortex induced by the secondary vortex is observed, rotating in the opposite direction to the secondary vortex the wake-induced transitions share the same route. The velocity fluctuations deviate from the optimal growth theory in the pre-transitional region. In the transitional region low-frequency disturbances penetrate the sheltering edge to generate streaks where the disturbance energy declines. In the turbulent region the logarithmic layer is formed, indicating that the turbulent equilibrium is established.

We report that vertical vibration with small amplitude and high frequency can tame convective heat transport in Rayleigh–Bénard convection in a turbulent regime. When vertical vibration is applied, a dynamically averaged ‘anti-gravity’ results that stabilizes the thermal boundary layer and inhibits the eruption of thermal plumes. This eventually leads to the attenuation of the intensity of large-scale mean flow and a significant suppression of turbulent heat transport. Accounting for both the thermally led buoyancy and the vibration-induced anti-gravitational effects, we propose an effective Rayleigh number that helps to extend the Grossmann–Lohse theory to thermal vibrational turbulence. The prediction of the reduction on both the Nusselt and Reynolds numbers obtained by the extended model is found to agree well with the numerical data. In addition, vibrational influences on the mean flow structure and the temporal evolution of Nusselt and Reynolds numbers are investigated. The non-uniform characteristic of vibration-induced ‘anti-gravity’ is discussed. The present findings provide a powerful basis for studying thermal vibrational turbulence and put forward a novel strategy for actively controlling thermal turbulence.

In this paper, we report that reversals of the large-scale circulation in two-dimensional Rayleigh–Bénard (RB) convection can be suppressed by imposing sinusoidally distributed heating to the bottom plate. We examine how the frequency of flow reversals depends on the dimensionless wavenumber $k$ of the spatial temperature modulation with various modulation amplitude $A$. For sufficiently large $k$, the flow reversal frequency is close to that in the standard RB convection under uniform heating. However, when $k$ decreases, the frequency of flow reversal gradually becomes lower and can even be largely reduced. Furthermore, we examine the growth of the corner roll and the global flow structure based on Fourier mode decomposition, and reveal that the size of the corner roll diminishes as the wavenumber decreases. The reason is that the regions occupied by the cold phase can absorb heat from the hot plumes and thus lower their temperature, which reduces the corner roll's kinetic energy input provided by the buoyancy force, and weakens the feeding process of the corner rolls. This results in the locking of the corner roll into a smaller region near the corner, making it harder for a reversal to occur. Using the concept of horizontal convection caused by non-uniform heating, we find a relevant parameter $k/A$ to describe briefly how the reversal frequency depends on wavenumber and modulation amplitude. The present work provides a new way to control the flow reversals in RB convection through modifying temperature boundary conditions.

Let $\mathbb {N}$ be the set of all nonnegative integers. For $S\subseteq \mathbb {N}$ and $n\in \mathbb {N}$ , let $R_S(n)$ denote the number of solutions of the equation $n=s_1+s_2$ , $s_1,s_2\in S$ and $s_1<s_2$ . Let A be the set of all nonnegative integers which contain an even number of digits $1$ in their binary representations and $B=\mathbb {N}\setminus A$ . Put $A_l=A\cap [0,2^l-1]$ and $B_l=B\cap [0,2^l-1]$ . We prove that if $C \cup D=[0, m]\setminus \{r\}$ with $0<r<m$ , $C \cap D=\emptyset $ and $0 \in C$ , then $R_{C}(n)=R_{D}(n)$ for any nonnegative integer n if and only if there exists an integer $l \geq 1$ such that $m=2^{l}$ , $r=2^{l-1}$ , $C=A_{l-1} \cup (2^{l-1}+1+B_{l-1})$ and $D=B_{l-1} \cup (2^{l-1}+1+A_{l-1})$ . Kiss and Sándor [‘Partitions of the set of nonnegative integers with the same representation functions’, Discrete Math. 340 (2017), 1154–1161] proved an analogous result when $C\cup D=[0,m]$ , $0\in C$ and $C\cap D=\{r\}$ .

External modulation on thermal convection has been studied extensively to achieve the control of flow structures and heat-transfer efficiency. In this paper, we carry out direct numerical simulations on Rayleigh–Bénard convection accounting for both the modulation of wall shear and roughness over the Rayleigh number range $1.0 \times 10^6 \le Ra \le 1.0 \times 10^8$, the wall shear Reynolds number range $0 \le Re_w \le 5000$, the aspect-ratio range $2 \le \varGamma \le 4{\rm \pi}$, and the dimensionless roughness height range $0 \le h \le 0.2$ at fixed Prandtl number $Pr = 1$. Under the combined actions of wall shear and roughness, with increasing $Re_w$, the heat flux is initially enhanced in the buoyancy-dominant regime, then has an abrupt transition near the critical shear Reynolds number $Re_{w,cr}$, and finally enters the purely diffusion regime dominated by shear. Based on the crossover of the kinetic energy production between the buoyancy-dominant and shear-dominant regimes, a physical model is proposed to predict the transitional scaling behaviour between $Re_{w,cr}$ and $Ra$, i.e. $Re_{w,cr} \sim Ra^{9/14}$, which agrees well with our numerical results. The reason for the observed heat-transport enhancement in the buoyancy-dominant regime is further explained by the fact that the moving rough plates introduce an external shear to strengthen the large-scale circulation (LSC) in the vertical direction and serve as a conveyor belt to increase the chances of the interaction between the LSC and secondary flows within cavities, which triggers more thermal plumes, efficiently transports the trapped hot (cold) fluids outside cavities.

We investigate the dynamic couplings between particles and fluid in turbulent Rayleigh–Bénard (RB) convection laden with isothermal inertial particles. Direct numerical simulations combined with the Lagrangian point-particle mode were carried out in the range of Rayleigh number $1\times 10^6 \le {Ra}\le 1 \times 10^8$ at Prandtl number ${Pr}=0.678$ for three Stokes numbers ${St_f}=1 \times 10^{-3}$, $8 \times 10^{-3}$ and $2.5 \times 10^{-2}$. It is found that the global heat transfer and the strength of turbulent momentum transfer are altered a small amount for the small Stokes number and large Stokes number as the coupling between the two phases is weak, whereas they are enhanced a large amount for the medium Stokes number due to strong coupling of the two phases. We then derived the exact relation of kinetic energy dissipation in the particle-laden RB convection to study the budget balance of induced and dissipated kinetic energy. The strength of the dynamic coupling can be clearly revealed from the percentage of particle-induced kinetic energy over the total induced kinetic energy. We further derived the power law relation of the averaged particles settling rate versus the Rayleigh number, i.e. $S_p/(d_p/H)^2{\sim} Ra^{1/2}$, which is in remarkable agreement with our simulation. We found that the settling and preferential concentration of particles are strongly correlated with the coupling mechanisms.

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 carry out direct numerical simulations of turbulent Rayleigh–Bénard convection in a square box with rough conducting plates over the Rayleigh number range $10^7\leqslant Ra\leqslant 10^9$ and the Prandtl number range $0.01\leqslant Pr\leqslant 100$. In Zhang et al. (J. Fluid Mech., vol. 836, 2018, R2), it was reported that while the measured Nusselt number $Nu$ is enhanced at large roughness height $h$, the global heat transport is reduced at small $h$. The division between the two regimes yields a critical roughness height $h_c$, and we now focus on the effects of the Prandtl number ($Pr$) on $h_c$. Based on the variations of $h_c$, we identify three regimes for $h_c(Pr)$. For low $Pr$, thermal boundary layers become thinner with increasing $Pr$. This makes the boundary layers easier to be disrupted by rough elements, leading to the decrease of $h_c$ with increasing $Pr$. For moderate $Pr$, the corner-flow rolls become much more pronounced and suppress the global heat transport via the competition between the corner-flow rolls and the large-scale circulation (LSC). As a consequence, $h_c$ increases with increasing $Pr$ due to the intensification of the corner–LSC competition. For high $Pr$, the convective flow transitions to the plume-controlled regime. As the rough elements trigger much stronger and more frequent plume emissions, $h_c$ again decreases with increasing $Pr$.

We aim to determine the correlation between parental rearing, personality traits, and obsessive–compulsive disorder (OCD) in different quantiles. In particular, we created an intermediary effect model in which parental rearing affects OCD through personality traits. All predictors were measured at the time of the survey, comprising parental rearing (paternal rearing and maternal rearing), demographics (grade and gender), and personality traits (neuroticism, extroversion, and psychoticism). These results suggest that (a) paternal emotional warmth was negatively correlated with OCD at the 0.40–0.80 quantile, while maternal emotional warmth was positively correlated with the OCD at the 0.45–0.69 quantile. (b) The correlation between negative parental rearing and OCD ranged from the 0.67 to 0.95 quantile for paternal punishment, 0.14–0.82 quantile for paternal overprotection, 0.05–0.36 and >0.50 quantile for maternal over-intervention and overprotection, and 0.08–0.88 quantile for maternal rejection. (c) Extroversion, neuroticism, and psychoticism were not only associated with OCD in a particular quantile but also mediated between parental rearing (namely parental emotional warmth, paternal punishment, paternal overprotection, maternal rejection, maternal over-intervention, and overprotection) and OCD. These findings provide targets for early interventions of OCD to improve the form of family education and personality traits and warrant validation.

Meiosis is a highly conserved process, and is responsible for the production of haploid gametes and generation of genetic diversity. We previously identified the transferrin receptor (TFRC) in the proteome profile of mice neonatal testes, indicating the involvement of the TFRC in meiosis. However, the exact molecular role of the TFRC in meiosis remains unclear. In this study, we aimed to determine the function of the TFRC in neonatal testicular development by TFRC knockdown using the testis culture platform. Our results showed high TFRC expression in 2-week testes, corresponding to the first meiotic division. Targeting TFRC using morpholino oligonucleotides resulted in clear spermatocyte apoptosis. In addition, we used the chromosomal spread technique to show that a deficiency of TFRC caused the accumulation of leptotene and zygotene spermatocytes, and a decrease of pachytene spermatocytes, indicating early meiotic arrest. Moreover, the chromosomes of TFRC-deficient pachytene spermatocytes displayed sustained γH2AX association, as well as SYCP1/SYCP3 dissociation beyond the sex body. Therefore, our results demonstrated that the TFRC is essential for the progression of spermatocyte meiosis, particularly for DNA double-stranded break repair and chromosomal synapsis.

Six acidic dykes were discovered surrounding the Laiziling pluton, Xianghualing area, in the western Cathaysia Block, South China. A number of captured zircons are found in two of these acidic dykes. By detailed U–Pb dating, Lu–Hf isotopes and trace-element analysis, we find that these zircons have ages clustered at c. 2.5 Ga. Two acidic dyke samples yielded upper intersection point 206U/238Pb ages of 2505 ± 42 Ma and 2533 ± 22 Ma, and weighted mean 207Pb/206Pb ages of 2500 ± 30 Ma and 2535 ± 16 Ma. The majority of these zircons have high (Sm/La)N, Th/U and low Ce/Ce* ratios, indicating a magmatic origin, but some grains were altered by later hydrothermal fluid. Additionally, the magmatic zircons have high Y, U, heavy rare earth element, Nb and Ta contents, indicating that their host rocks were mainly mafic rocks or trondhjemite–tonalite–granodiorite rock series. Equally, their moderate Y, Yb, Th, Gd and Er contents also indicate that a mafic source formed in a continental volcanic-arc environment. These zircons have positive ϵHf(t) values (2.5–6.9) close to zircons from the depleted mantle, with TDM (2565–2741 Ma) and TDM2 (2608–2864 Ma) ages close to their formation ages, indicating that these zircons originated directly from depleted mantle magma, or juvenile crust derived from the depleted mantle in a very short period. We therefore infer that the Cathaysia Block experienced a crustal growth event at c. 2.5 Ga.

In this paper, we report that reversals of large-scale circulation in two-dimensional Rayleigh–Bénard convection could be suppressed or enhanced by imposing local constant-temperature control on sidewalls. When the control area is away from the centre of the sidewalls, the control can successfully eliminate the flow reversal if the size of the control region is large enough. With a proper location, the width can be as small as 1 % of the system size. When the control region is located around the centre, the control may enhance the flow reversal. It may also stimulate the occurrence of a double-roll mode when the control is located in the centre. Explanations are also discussed based on the twofold effects of the control region on the nearby plumes and the concept of symmetry. The present work provides a new way to control the flow reversals in Rayleigh–Bénard convection through modifying sidewall boundary conditions.

Bloodstream infection (BSI), caused by Klebsiella pneumoniae, is associated with high morbidity and mortality, where the pks gene cluster plays a major role in their occurrence and prevalence. Information on the prevalence and characteristics of this gene cluster in K. pneumoniae is currently limited in mainland China. We therefore undertook a multicentre longitudinal study which revealed the prevalence, overall, community-onset and hospital-acquired BSI to be 20.5%, 28.3% and 13.0%, respectively. Compared to pks-negative, pks-positive isolates were significantly more susceptible to antimicrobial agents with a low incidence (5.1%) of multidrug-resistance and with infrequent extended-spectrum beta-lactamase (ESBL) production. Among pks-positive isolates, ST23 (78/117) and ST65 (20/117) were the dominant sequence types, and the majority harboured virulence genes. Community-onset BSI patients infected with pks-positive isolates had a higher proportion of liver abscesses and a lower proportion of biliary obstructions (P < 0.05). The pks-positive isolates were mostly sporadic in the phylogenetic tree, with a 65.8 and 47.0 average allele difference between Clade 1 and Clade 2, respectively. We concluded that although pks-positive K. pneumoniae were generally susceptible to antimicrobials, the high prevalence of such isolates in community cases and the genotoxicity, merits further investigation.