We systematically study the dissipative anomaly in compressible magnetohydrodynamic (MHD) turbulence using direct numerical simulations, and show that the total dissipation remains finite as viscosity diminishes. The dimensionless dissipation rate $\mathcal {C}_{\varepsilon }$ fits well with the model $\mathcal {C}_{\varepsilon } = \mathcal {C}_{\varepsilon,\infty } + \mathcal {D}/R_L^-$ for all levels of flow compressibility considered here, where $R_L^-$ is the generalized large-scale Reynolds number. The asymptotic value $\mathcal {C}_{\varepsilon,\infty }$ describes the total energy transfer flux, and decreases with increase of the flow compressibility, indicating non-universality of the dimensionless dissipation rate in compressible MHD turbulence. After introducing an empirically modified dissipation rate, the data from compressible cases collapse to a form similar to the incompressible MHD case depending only on the modified Reynolds number.

Direct numerical simulations of temporally developing mixing layers have been performed to investigate the effects of compressibility on statistics and structures near the interfaces of high- and low-speed large-scale structures (LSSs), covering a range of convective Mach numbers from $M_c=0.2$ to $1.8$ and Taylor Reynolds numbers up to 290. The interfaces of LSSs are directly defined by the isosurface of zero fluctuating streamwise velocity. The characteristic velocity jump at the interfaces grows rapidly in the transition stage and then decreases until reaching a gradual self-similar stage. The evolution of velocity jump is evidently delayed as the convective Mach number increases. The interface layer is formed by the shearing motion of neighbouring LSSs. A small vortical motion characterized by the Kolmogorov scale is induced in the interface layer by shear-dominated outer regions. The strengths of outer shearing motion and central vortical motion are greater at the leading edge, but smaller at the trailing edge, which is also reflected in a larger velocity jump at the leading edge. As the convective Mach number increases, the small-scale vortical structure is obviously suppressed by compressibility. At high convective Mach number $M_c=1.8$, the compressive shear-dominated outer regions are linked with a sheet-like structure passing through the centre of the expansion region corresponding to a small-scale vortical structure. The compressibility and shearing strength near the interface are highly dependent on the interface orientation.

Based on the characteristics of large-scale plume currents in rotating plane Poiseuille flows (RPPF), an injection/suction control strategy is introduced to augment the intensity of plume currents and improve the turbulent transport of passive scalar. The control strategy maintains the conventional non-penetrative condition on the stable side, and applies a wall-normal velocity that varies in the spanwise direction on the unstable side. For comparison, the control with non-penetrative condition on the unstable side and injection/suction on the stable side is also examined. The RPPF at $Re_\tau$ ranging from $180$ to $300$ and $Ro_\tau$ ranging from $0$ to $30$ are studied. Injection/suction with fixed root-mean-squared wall-normal velocity (below $1\,\%$ of the bulk mean velocity) at the wall is considered. Direct numerical simulations reveal that the injection/suction on the stable side has minor effect for all cases studied, whereas the injection/suction with properly distributed slots on the unstable side can significantly enhance the large-scale plume currents and the turbulent transport efficiency at small or moderate rotation numbers. This is attributed to the fact that plumes are generated on the unstable side, and controlling the origin of plumes is more effective. A proper strategy for maximum enhancement of turbulent transport is limiting the distance between injection slots equal to or slightly larger than the intrinsic distance between plume currents. This strategy is valid for all physical and computational parameters currently considered, and is potentially applicable over a much wider parameter range.

The wall cooling effect on the spectra and structures of thermodynamic variables are investigated in hypersonic turbulent boundary layers. The density and temperature can be divided into the acoustic and entropic modes based on the Kovasznay decomposition. The intensities of the pressure and the acoustic modes of density and temperature attain the maximum values near the wall, while those of the entropy and the entropic modes of density and temperature achieve their primary peaks near the edge of the boundary layer. In the near-wall region, the pressure and the acoustic modes of density and temperature are significantly enhanced when the wall is strongly cooled, which can be attributed to the appearance of the travelling-wave-like alternating positive and negative structures. Moreover, the intensities of the entropy and the entropic modes of density and temperature become stronger near the wall as the wall temperature decreases, due to the appearance of the streaky entropic structures (SES). The SES are mainly caused by the advection effect of the strong positive wall-normal gradient of the mean temperature associated with ejection and sweep events. It is also found that the profiles of the intensities of the entropy, density and temperature are similar to each other far from the wall, which is mainly due to the reason that the entropic modes are dominant in the fluctuating density and temperature in the far-wall region. The acoustic modes of density and temperature only have significant contributions in the near-wall region.

Data reconstruction of rotating turbulent snapshots is investigated utilizing data-driven tools. This problem is crucial for numerous geophysical applications and fundamental aspects, given the concurrent effects of direct and inverse energy cascades. Additionally, benchmarking of various reconstruction techniques is essential to assess the trade-off between quantitative supremacy, implementation complexity and explicability. In this study, we use linear and nonlinear tools based on the proper orthogonal decomposition (POD) and generative adversarial network (GAN) for reconstructing rotating turbulence snapshots with spatial damages (inpainting). We focus on accurately reproducing both statistical properties and instantaneous velocity fields. Different gap sizes and gap geometries are investigated in order to assess the importance of coherency and multi-scale properties of the missing information. Surprisingly enough, concerning point-wise reconstruction, the nonlinear GAN does not outperform one of the linear POD techniques. On the other hand, the supremacy of the GAN approach is shown when the statistical multi-scale properties are compared. Similarly, extreme events in the gap region are better predicted when using GAN. The balance between point-wise error and statistical properties is controlled by the adversarial ratio, which determines the relative importance of the generator and the discriminator in the GAN training.

The mean flow behaviour of a turbulent boundary layer over rough walls is expected to exhibit symmetries that govern the flow dynamics. In particular, when roughness elements are arranged in a spanwise symmetric manner, the mean flow above them should also exhibit spanwise symmetry. This symmetrical consideration has garnered substantial empirical support. We conduct direct numerical simulations (DNS) of flow over aligned cube arrays to test such symmetry considerations further. We vary the surface coverage density from 0.25 % to 6.25 %, and employ an averaging time of about 100 large-eddy turnover times, which is longer than the typical averaging time in prior DNS studies of rough-wall boundary layers. The results suggest the presence of spanwise asymmetry in the mean flow. Specifically, we observe the development of a prominent secondary vortex on one side of the cubical roughness, accompanied by a relatively smaller secondary vortex on the other side. This asymmetry becomes most pronounced when the surface coverage density is approximately 0.59 %, and diminishes as the coverage density approaches either a low or a high value. We also establish that this mean flow asymmetry is robust across variations in the domain size, the initial condition, and the placement of the cubes in the spanwise direction.

Theoretical models have ambiguous predictions on how workplace gender composition affects the incidence of marriage. Marital search theory suggests that having more opportunities for interactions between members of the opposite gender increases the likelihood of marriage. Yet, according to overload choice theory, people with more options could actually delay or forgo marriage if the increase in the number of choices makes it more difficult for them to make marriage decisions. I explore how changes in the gender composition within occupation and industry over the past 40 years affect marriage decisions. I find that a higher share of opposite gender coworkers within a person's occupation-industry is associated with a decreased likelihood of ever having been married.

The Reynolds number and wall cooling effects on correlations between the thermodynamic variables are systematically investigated in hypersonic turbulent boundary layers by direct numerical simulations. The correlations between the thermodynamic variables and the streamwise velocity are also analysed. The Kovasznay decomposition is introduced to decompose the fluctuating density and temperature into the acoustic and entropic modes. It is found that in the strongly cooled wall cases, the travelling-wave-like alternating positive and negative structures (TAPNSs) are found in the fluctuating pressure and the acoustic modes of density and temperature, and the streaky entropic structures (SESs) are identified in the fluctuating entropy and the entropic modes of density and temperature near the wall. Furthermore, both the acoustic and the entropic modes of density and temperature give significant contributions to the correlations involving density and temperature in the near-wall region, while these correlations are almost totally contributed by the entropic modes in the far-wall region. The entropic modes of the density and temperature are almost linearly correlated with the fluctuating entropy. Therefore, the fact that the fluctuating entropy is strongly correlated with the fluctuating density and temperature far from the wall is mainly due to the dominance of the entropic modes in the fluctuating density and temperature. Moreover, the fluctuating temperature is strongly positively correlated with the fluctuating streamwise velocity near the wall in strongly cooled wall cases, which can be ascribed to the appearance of the TAPNSs and SESs.

In this paper, we propose artificial-neural-network-based (ANN-based) nonlinear algebraic models for the large-eddy simulation (LES) of compressible wall-bounded turbulence. An innovative modification is applied to the invariants and the tensor bases of the nonlinear algebraic models through using the local grid widths along each direction to normalise the corresponding gradients of the flow variables. Furthermore, the dimensionless model coefficients are determined by the ANN method. The modified ANN-based nonlinear algebraic model (MANA model) has much higher correlation coefficients and much lower relative errors than the dynamic Smagorinsky model (DSM), Vreman model and wall-adapting local eddy-viscosity model in the a priori test. The significantly more accurate estimations of the mean subgrid-scale (SGS) fluxes of the kinetic energy and temperature variance are also obtained by the MANA models in the a priori test. Furthermore, in the a posteriori test, the MANA model can give much more accurate predictions of the flow statistics and the mean SGS fluxes of the kinetic energy and the temperature variance than other traditional eddy-viscosity models in compressible turbulent channel flows with untrained Reynolds numbers, Mach numbers and grid resolutions. The MANA model has a better performance in predicting the flow statistics in supersonic turbulent boundary layer. The MANA model can well predict both direct and inverse transfer of the kinetic energy and temperature variance, which overcomes the inherent shortcoming that the traditional eddy-viscosity models cannot predict the inverse energy transfer. Moreover, the MANA model is computationally more efficient than the DSM.

Enstrophy production and flow topology are numerically investigated for statistically stationary compressible isotropic turbulence in vibrational non-equilibrium with a large-scale thermal forcing. The net enstrophy production term is decomposed into solenoidal, dilatational and isotropic dilatational terms based on the Helmholtz decomposition. From the full flow field perspective, the net enstrophy production mainly stems from the solenoidal term. For the dilatational and isotropic dilatational terms, although their local magnitudes can be considerable, the positive values in the compression region and the negative values in the expansion region cancel out on average. For the solenoidal component of the deviatoric strain-rate tensor, the statistical properties of its eigenvalues and alignments between vorticity and its eigenvectors are nearly independent of the local dilatation and vibrational relaxation. The solenoidal components of enstrophy production along three eigendirections are thus mainly affected by the vorticity. For the dilatational component of deviatoric strain-rate tensor, the statistical properties of its eigenvalues and alignments between vorticity and its eigenvectors closely relate to the local dilatation and vibrational relaxation. The dilatational components of enstrophy production along three eigendirections are therefore affected by the vorticity, eigenvalues and alignments between the vorticity and eigenvectors. The topological classification proposed by Chong et al. (Phys. Fluids, vol. 2, issue 5, 1990, pp. 765–777) is employed to decompose the flow field into various flow topologies. In the strong compression and strong expansion regions, the relaxation effects on the volume fractions of flow topologies and their relative contributions to the local enstrophy production are significant.

The effects of compressibility on the statistics and coherent structures of a temporally developing mixing layer are studied using numerical simulations at convective Mach numbers ranging from $M_c=0.2$ to $1.8$ and at Taylor Reynolds numbers up to 290. As the convective Mach number increases, the streamwise dissipation becomes more effective to suppress the turbulent kinetic energy. At $M_c=1.8$, the streamwise dissipation increases much faster than the other two components in the transition region, even larger than pressure–strain redistribution, correlating with the streamwise elongated vortical structures at a higher level of compressibility. We confirm the existence of the large-scale high- and low-speed structures in the mixing layers, which accompany the spanwise Kelvin–Helmholtz rollers at low convective Mach number and dominate the mixing layer at higher convective Mach number. Conditional statistics demonstrate that the large-scale low-speed structures are lifted upwards by a pair of counter-rotating quasi-streamwise rollers flanking the structures. The small-scale vortical structures have an apparent preference for clustering into the top of the low-speed regions, which is directly associated with high-shearing motions on top of the low-speed structures. The high-speed structures statistically exhibit central symmetry with the low-speed structures. The statistics and dynamics of large-scale high- and low-speed structures in the compressible mixing layers resemble those in the outer region of the turbulent boundary layers, which reveals the universality of the large-scale structures in free shear and wall-bounded turbulence. A conceptual model is introduced for the large-scale high- and low-speed structures in turbulent mixing layers.

The decompositions of the skin-friction and heat-transfer coefficients based on the twofold repeated integration in hypersonic transitional and turbulent boundary layers are analysed to give some major reasons of the overshoot phenomena of the wall skin friction and heat transfer. It is shown that the overshoot of the skin-friction coefficient is mainly caused by the drastic change of the mean velocity profiles, especially the strong negative streamwise gradient of the mean streamwise velocity far from the wall; and the overshoot of the heat-transfer coefficient is primarily due to the viscous dissipation, especially the strong positive vertical gradient of the mean streamwise velocity near the wall. These observations are different from the previous observations that the Reynolds shear stress and Reynolds heat flux are the reasons, respectively. Further investigations show that the above observations are independent of the set-up of the wall blowing and suction parameters, which indicates the universality of the major reasons of the overshoot phenomena in our numerical simulations. In the hypersonic turbulent boundary layers, it is observed that the strongly cooled wall temperature and the high Mach number can slightly enhance the contribution of the Reynolds shear stress, and weaken the contribution of the mean convection, mainly due to the strong compressibility effect. Moreover, the magnitudes of the relative contributions of the mean convection, pressure dilatation, viscous dissipation and the Reynolds heat flux increase as the wall temperature increases.

We perform a theoretical investigation of the instability of a helical vortex filament beneath a free surface in a semi-infinite ideal fluid. The focus is on the leading-order free-surface boundary effect upon the equilibrium form and instability of the vortex. This effect is characterised by the Froude number $F_r = U(gh^*)^{-{1}/{2}}$ where $g$ is gravity, and $U = \varGamma /(2{\rm \pi} b^*)$ with $\varGamma$ being the strength, $2{\rm \pi} b^*$ the pitch and $h^*$ the centre submergence of the helical vortex. In the case of $F_r \rightarrow 0$ corresponding to the presence of a rigid boundary, a new approximate equilibrium form is found if the vortex possesses a non-zero rotational velocity. Compared with the infinite fluid case (Widnall, J. Fluid Mech., vol. 54, no. 4, 1972, pp. 641–663), the vortex is destabilised (or stabilised) to relatively short- (or long-)wavelength sub-harmonic perturbations, but remains stable to super-harmonic perturbations. The wall-boundary effect becomes stronger for smaller helix angle and could dominate over the self-induced flow effect depending on the submergence. In the case of $F_r > 0$, we obtain the surface wave solution induced by the vortex in the context of linearised potential-flow theory. The wave elevation is unbounded when the $m$th wave mode becomes resonant as $F_r$ approaches the critical Froude numbers ${\mathcal {F}} (m) = (C_0^*/U)^{-1} (mh^*/b^*)^{-{1}/{2}}$, $m=1, 2, \ldots,$ where $C_0^*$ is the induced wave speed. We find that the new approximate equilibrium of the vortex exists if and only if $F_r < {\mathcal {F}}(1)$. Compared with the infinite fluid and $F_r \rightarrow 0$ cases, the wave effect causes the vortex to be destabilised to super-harmonic and long-wavelength sub-harmonic perturbations with generally faster growth rate for greater $F_r$ and smaller helix angle.

The analogy between rotating shear flow and thermal convection suggests the existence of plumes, inertial waves and plume currents in plane Poiseuille flow under spanwise rotation. The existence of these flow structures is examined with the results of three-dimensional and two-dimensional three-component direct numerical simulations. The dynamics of plumes near the unstable side is embodied in a truncated exponential distribution of turbulent fluctuations. For large rotation numbers, inertial waves are identified near the stable side, and these can be used to explain the abnormal flow statistics, such as the large root-mean-square of the streamwise velocity fluctuation and the nearly negligible Reynolds shear stress. For small or medium rotation numbers, plumes generated from the unstable side form large-scale plume currents and the patterns of the plume currents show different capabilities in scalar transport.

Kinetic energy flux (KEF) is an important physical quantity that characterizes cascades of kinetic energy in turbulent flows. In large-eddy simulation (LES), it is crucial for the subgrid-scale (SGS) model to accurately predict the KEF in turbulence. In this paper, we propose a new eddy-viscosity SGS model constrained by the properly modelled KEF for LES of compressible wall-bounded turbulence. The new methodology has the advantages of both accurate prediction of the KEF and strong numerical stability in LES. We can obtain an approximate KEF by the tensor-diffusivity model, which has a high correlation with the real value. Then, using the artificial neural network method, the local ratios between the real KEF and the approximate KEF are accurately modelled. Consequently, the SGS model can be improved by the product of that ratio and the approximate KEF. In LES of compressible turbulent channel flow, the new model can accurately predict mean velocity profile, turbulence intensities, Reynolds stress, temperature–velocity correlation, etc. Additionally, for the case of a compressible flat-plate boundary layer, the new model can accurately predict some key quantities, including the onset of transitions and transition peaks, the skin-friction coefficient, the mean velocity in the turbulence region, etc., and it can also predict the energy backscatters in turbulence. Furthermore, the proposed model also shows more advantages for coarser grids.

In this paper, we propose a series expansion of the baroclinic torque in low-Mach-number flows, so that the accuracy and universality of any buoyancy term could be examined analytically, and new types of buoyancy terms could be constructed and validated. We first demonstrate that the purpose of introducing a buoyancy term is to approximate the baroclinic torque, and straightforwardly the error of any buoyancy term could be defined with the deviation of its curl from the corresponding baroclinic torque. Then a regular perturbation method is introduced for the elliptic equation of the hydrodynamic pressure in low-Mach-number flows, resulting in a sequence of Poisson equations, whose solutions lead to the series representation of the baroclinic torque and the new types of buoyancy terms. It is found that the frame invariance of the momentum equation is maintained with one of the new types of buoyancy terms. With the error definition of buoyancy terms and the series representation of the baroclinic torque, the validity and accuracy of previous and new buoyancy terms are examined. Finally, numerical simulations confirm that, with a decreasing density variation or an increasing order of our new buoyancy term, the simplified equations can converge to the original low-Mach-number equations.

The effect of wall temperature on the transfer of kinetic energy in a hypersonic turbulent boundary layer for different Mach numbers and wall temperature ratios is studied by direct numerical simulation. A cold wall temperature can enhance the compressibility effect in the near-wall region through increasing the temperature gradient and wall heat flux. It is shown that the cold wall temperature enhances the local reverse transfer of kinetic energy from small scales to large scales, and suppresses the local direct transfer of kinetic energy from large scales to small scales. The average filtered spatial convection and average filtered viscous dissipation are dominant in the near-wall region, while the average subgrid-scale flux of kinetic energy achieves its peak value in the buffer layer. It is found that the wall can suppress the inter-scale transfer of kinetic energy, especially for the situation of a cold wall. A strong local reverse transfer of fluctuating kinetic energy is identified in the buffer layer in the inertial range. Helmholtz decomposition is applied to analyse the compressibility effect on the subgrid-scale flux of kinetic energy. A strong transfer of the solenoidal component of fluctuating kinetic energy is identified in the buffer layer, while a significant transfer of the dilatational component of fluctuating kinetic energy is observed in the near-wall region. It is also shown that compression motions have a major contribution to the direct transfer of fluctuating kinetic energy, while expansion motions play a marked role in the reverse transfer of fluctuating kinetic energy.