@tsinghua.org.cn
Department of Aerospace Engineering
Tsinghua University
Turbulent boundary layer;
Computational Fluid Mechanics
Scopus Publications
Ming Yu, Bo Li, QingQing Zhou, Dong Sun, and XianXu Yuan
Elsevier BV
Ming Yu, SiWei Dong, QiLong Guo, ZhiGong Tang, XianXu Yuan, and ChunXiao Xu
Cambridge University Press (CUP)
Very-large-scale motions are commonly observed in moderate- and high-Reynolds-number wall turbulence, constituting a considerable portion of the Reynolds stress and skin friction. This study aims to investigate the behaviour of these motions in high-speed and high-Reynolds-number turbulent boundary layers at varying Mach numbers. With the aid of high-precision numerical simulations, numerical experiments and theoretical analysis, it is demonstrated that the very-large-scale motions are weakened in high-Mach-number turbulence at the same friction Reynolds numbers, leading to the reduction in turbulent kinetic energy in the outer region. Conversely, the lower wall temperature enhances the very-large-scale motions but shortens the scale separation between the structures in the near-wall and outer regions.
Ming Yu, Davide Modesti, and Sergio Pirozzoli
Cambridge University Press (CUP)
We study turbulent flow in open channels with a free surface and rectangular cross-section, for various Reynolds numbers and duct aspect ratios. Direct numerical simulations are used to obtain accurate characterization of the secondary motions, which are found to be more intense than in closed ducts, and to scale with the bulk, rather than with the friction velocity. A notable feature of open-duct flows is the presence of a velocity dip, namely the peak velocity is achieved at some depth underneath the free surface. We find that the depth of the velocity peak increases with the Reynolds number, and correspondingly the flow becomes more symmetric with respect to the horizontal midplane. This is also confirmed from the change of the topology of the secondary motions, which exhibit a strong corner circulation at the free-surface/wall corners at low Reynolds number, which, however, weakens at higher $Re$ . The structure of the mean velocity field is such that the log law applies with good approximation in the direction normal to the nearest wall, which allows us to explain why predictive friction formulae based on the hydraulic diameter concept are successful. Additional analysis shows that the secondary motions account for a large fraction of the frictional drag (up to $15$ %).
Ming Yu, Qingqing Zhou, Hongmin Su, Qilong Guo, and Xianxu Yuan
Springer Science and Business Media LLC
Ming Yu, QingQing Zhou, SiWei Dong, XianXu Yuan, and ChunXiao Xu
Cambridge University Press (CUP)
In the present study, we investigate the compressibility effects in supersonic and hypersonic turbulent boundary layers under the influence of wall disturbances by exploiting direct numerical simulation databases at Mach numbers up to 6. Such wall disturbances enforce extra Reynolds shear stress on the wall and induce mean streamline curvature in rough wall turbulence that leads to the intensification of turbulent motions in the outer region. The turbulent and fluctuating Mach numbers, the density and the velocity divergence fluctuation intensities suggest that the compressibility effects are enhanced by the increment of the free-stream Mach number and the implementation of the wall disturbances. The differences between the Reynolds and Favre average due to the density fluctuations constitute approximately $9\\,\\%$ of the mean velocity close to the wall and $30\\,\\%$ of the Reynolds stress near the edge of the boundary layer, indicating their non-negligibility in turbulent modelling strategies. The comparatively strong compressive events behaving as eddy shocklets are observed at the free-stream Mach number of $6$ only in the cases with wall disturbances. By further splitting the velocity into the solenoidal and dilatational components with the Helmholtz decomposition, we found that the dilatational motions are organized as travelling wave packets in the wall-parallel planes close to the wall and as forward inclined structures in the form of radiated waves in the vertical planes. Despite their increased magnitudes and higher portion in the Reynolds normal and shear stresses, the dilatational motions show no tendency of contributing significantly to the skin friction and the production of turbulent kinetic energy due to their mitigation by the cross-correlation between the solenoidal and dilatational velocity components.
PengXin Liu, JunYang Li, HongMin Su, Dong Sun, Ming Yu, and XianXu Yuan
Elsevier BV
Ming Yu, Dong Sun, QingQing Zhou, PengXin Liu, and XianXu Yuan
AIP Publishing
In the present study, we investigate the evolution of turbulent statistics and coherent structures in hypersonic turbulent boundary layers at the Mach number of 5 impinged by oblique shock waves generated by the wedge with the angles of 14°, 10°, and 6°, inducing strong, mild, and incipient flow separation, by exploiting direct numerical simulation databases, for the purpose of revealing the underlying flow physics that are of significance to turbulent modeling. We found that the large-scale structures are amplified within the interaction zone, manifested in the form of large-scale low- and high-speed streaks with the spanwise length scale of boundary layer thickness, and gradually decay downstream, the process of which is extremely long. The abrupt variation in the characteristic length, time, and velocity scales as well as the incompatible viscous dissipation of the mean and turbulent kinetic energy results in the incorrect predictions by the Reynolds-Averaged Navier–Stokes (RANS) equation simulations, provided the models are established based on solving the transport equations of the turbulent kinetic equation and its viscous dissipation (k−ε or k−ω models, for instance). To amend this issue, we propose to refine the parameters in the model as the functions of wall pressure, the flow quantities related to multiple flow features. The RANS simulations with the k−ω SST model utilizing the proposed refinement improve greatly the accuracy of the skin friction, wall heat flux, and Reynolds shear stress downstream of the interaction zone, and the wall pressure distributions in hypersonic turbulence over compression ramp, suggesting its promising prospect in engineering applications.
Ming YU, Yalu FU, Zhigong TANG, Xianxu YUAN, and Chunxiao XU
Elsevier BV
Ming Yu, SiWei Dong, PengXin Liu, ZhiGong Tang, XianXu Yuan, and ChunXiao Xu
Cambridge University Press (CUP)
The oblique shock impinging on the supersonic turbulent boundary layer leads to a mixing layer and the emergence of large-scale coherent structures within the interaction zone which leave significant velocity defect and turbulence amplification downstream. In the present study, we investigate the turbulence recovery in the post-shock region by exploiting direct numerical simulation data of the oblique-shock/turbulent boundary layer interaction flow at the incoming Mach number of $2.28$ and the shock angle of $33.2^\\circ$ , with special attention paid to the contribution of the mixing layer and large-scale structures to flow dynamics. For that purpose, we propose to split the mean velocity, Reynolds stresses and spanwise spectra into a canonical portion that is constructed according to the statistics of canonical turbulent boundary layers, and a mixing-layer-induced portion. We found that the hidden mixing layer grows with the boundary layer thickness and that the induced mean shear and Reynolds stresses decay at different rates. The mean velocity recovers to the canonical profiles at a distance of 13 boundary layer thicknesses downstream where the mixing-layer-induced mean shear ceases to have strong impacts. The recovery of Reynolds stresses requires 10 boundary layer thicknesses in the near-wall region but a much longer streamwise extent in the outer region due to the slow decay of large-scale motions. These large-scale motions superpose on the near-wall turbulence, intensifying the turbulent fluctuations, yet having a trivial impact on the skin friction, for the contribution of the mixing-layer-induced mean shear and Reynolds shear stress are balanced by the advection term. We further establish a simple physical model capable of approximately predicting the streamwise evolution of mixing-layer-induced mean shear and turbulent kinetic energy. This model suggests that the complete recovery of turbulence in the outer region requires a streamwise extent of approximately 50 boundary layer thicknesses.
Ming Yu, PengXin Liu, ZhiGong Tang, XianXu Yuan, and ChunXiao Xu
AIP Publishing
In the present study, we perform direct numerical simulations to investigate the spatial development and basic flow statistics in the supersonic turbulent boundary layers at the free-stream Mach number of 2.0 over smooth and disturbed walls, the latter of which enforces extra Reynolds shear stress in the streamwise direction to emulate the drag increment and mean streamline curvature effects of rough walls. Such disturbances escalate the growth rate of turbulent boundary layer thickness and the shape factor. It is found that under the rescaled global coordinate, the mean velocity, Reynolds stress, and pressure fluctuation variance manifest outer-layer similarity, whereas the average and fluctuation variances of temperature and density do not share such a property. Compressibility effects are enhanced by the wall disturbances, yet not sufficiently strong to directly impact the turbulent kinetic energy transport under the presently considered flow parameters. The generalized Reynolds analogy that relates the mean velocity and temperature can be satisfied by incorporating the refinement in modifying the generalized recovery coefficient, and that associates the fluctuating velocity and temperature work reasonably well, indicating the passive transport of temperature fluctuations. The dispersive motions are dominant and decay exponentially below the equivalent sand grain roughness height ks, above which the wall disturbances are distorted to form unsteady motions responsible for the intensified density and pressure fluctuations in the free-stream traveling isentropically as the acoustic radiations.
YaLu Fu, QingQing Zhou, Ming Yu, HongMin Su, QiLong Guo, and XianXu Yuan
Informa UK Limited
HengYu Cai, Ming Yu, Dong Sun, ZhengYin Ye, PengXin Liu, and XianXu Yuan
AIP Publishing
In the present study, we investigate influences of shock intensity on wall pressure fluctuations by performing direct numerical simulations of supersonic turbulence boundary layers over compression ramps with different turning angles. We found that as the turning angle increases, low-frequency motions of the separation shock are enhanced, accompanied by enlarged energetic pressure structures with lower convection velocities. By inspecting wavenumber-frequency spectra under the assumption of streamwise homogeneity, we further identified two energetic modes convected at different velocities. The one with the lower convection velocity, namely, the “slow mode,” inherited from the upstream pressure fluctuations of the turbulent boundary layer, is decelerated when passing through the oblique shock, during which the “rapid mode” with pressure fluctuations convected at higher speeds are generated. The increasing turning angle decelerates the slow mode and intensifies the fast mode. The reconstruction of the flow field suggests that the rapid mode is associated with the shear layer generated adjacent to the interaction zone, while the slow mode is associated with the Görtler vortices on the ramp.
Ming Yu, MingXiang Zhao, ZhiGong Tang, XianXu Yuan, and ChunXiao Xu
Cambridge University Press (CUP)
Turbulence amplification and the large-scale coherent structures in shock wave/turbulent boundary layer interaction flows have been studied at length in previous research, while the direct association between these two flow features is still lacking. In the present study, the transport equation of turbulent kinetic energy spectra is derived and utilized to analyse the scale-by-scale energy budget across the interaction zone, enabling us to reveal the association between the genesis of the large-scale motions and the turbulence amplification. For the presently considered flow with incipient shock-induced separation, we identified in turbulent kinetic energy spectra distribution that the most energetic motions are converted from the near-wall small-scale motions to large-scale motions consisting of velocity streaks and cross-stream circulations as they go through the interaction zone. The amplification of streamwise velocity fluctuation is triggered first, resulting in the emergence of large-scale velocity streaks, which is attributed to the adverse pressure gradient, as indicated by the spectra of the production term. The energy carried by large-scale velocity streaks is transferred to other velocity components by the pressure-strain term, producing large-scale cross-stream circulations. When large-scale motions are convected downstream, their energy is transferred via turbulent cascade to smaller scales and dissipated by viscosity. The spanwise uniform fluctuations, reminiscent of the unsteadiness of the separation bubble, are contributed primarily by the inter-scale energy transfer from the finite spanwise scale motions.
Ming Yu, Yalu Fu, Pengxin Liu, Zhigong Tang, Xianxu Yuan, and Chunxiao Xu
Springer Science and Business Media LLC
JunYang Li, Ming Yu, Dong Sun, PengXin Liu, and XianXu Yuan
AIP Publishing
In this paper, we investigate the differences in wall heat transfer between the low- and high-enthalpy turbulent boundary layers by exploiting direct numerical simulation databases of hypersonic turbulent boundary layers at the free-stream Mach number of 4.5 and the friction Reynolds number of 800. For that purpose, we refine the integral formula of decomposing the wall heat flux proposed by Sun et al. [“A decomposition formula for the wall heat flux of a compressible boundary layer,” Adv. Aerodyn. 4, 1–13 (2022)], enabling us to scrutinize the contribution of different physical processes. Statistical results show that the mean wall heat transfer is primarily contributed by the heat conduction, the turbulent heat transfer, viscous dissipation of mean kinetic energy, and turbulent kinetic energy production. Among these processes, the contribution of the turbulent heat flux in the high-enthalpy case is 10% higher than that in the low-enthalpy case. Such discrepancy is caused by the turbulent–chemistry interaction consisting of velocity and species mass fraction fluctuations. Coherent structures in the conditionally averaged fields related to this process reveal that the sweep in the viscous sublayer and ejection in the logarithmic layer bringing the hot fluid downward and upward, respectively, significantly alter the distribution of the species mass fraction. The wall heat flux fluctuations are slightly enhanced in the high-enthalpy flows, which is ascribed to be the intensification of traveling wave packets.
Siwei Dong, Fulin Tong, Ming Yu, Jianqiang Chen, Xianxu Yuan, and Qian Wang
AIP Publishing
The negative and positive fluctuations of wall shear stress [Formula: see text] and wall heat flux [Formula: see text] can be related to the wall-attached paired up large-scale velocity and temperature streaks. It is justifiable to infer the spatially paired-up coexistence of those wall flow quantities. The present study aims at testifying this hypothesis. We establish such relations between the negative and positive wall shear stress by exploiting a direct numerical simulation database over heated and cooled walls at the friction Reynolds number of 800 and the Mach number of 2.25. The clustering method is adopted for the search of the in-pair structures. It is found that the τx- and qw-structures are less self-similar for flows over cold walls. As they become wider, the τx-structures are increasingly more streamwise stretched, while the trend is reversed for qw-structures. τx-structures of opposite signs are paired up and aligned in the spanwise directions as the wall-attached streamwise velocity, and are left behind by streamwise rollers. The relative position between qw-structures of opposite signs, on the other hand, is sensitive to the wall temperature. Scrutinizing the statistical structures, we elucidate that such spatial coherence is determined by the meandering of velocity streaks that yields strong streamwise gradients of the streamwise velocity.
Ming Yu, Alessandro Ceci, and Sergio Pirozzoli
Elsevier BV
Ming Yu, Peng Xin Liu, Ya Lu Fu, Zhi Gong Tang, and Xian Xu Yuan
AIP Publishing
Wall shear stress, pressure, and heat flux are of significant importance in engineering applications. In this two-part study, we investigate the compressibility effects on wall shear stress, pressure, and heat flux fluctuations in compressible wall-bounded turbulence by exploiting direct numerical simulation databases. In Paper I, we primarily deal with the one-point statistics, whereas in this second part, we report the effects of compressibility on the frequency spectra, wavenumber-frequency spectra of these flow quantities, and the two-point cross-correlations between them. It is found that the scaling laws of the spectra at low and high frequencies are retained as those of incompressible flows, whereas the spectra intensities at mid frequencies increase with the enhancement of compressibility effects, which is identified to be related to the ever-predominating traveling wave packets. These wave packets are convected downstream at the same velocity of [Formula: see text] as that of pressure fluctuations, higher than that of the streaky structures [Formula: see text] ( Ub the bulk velocity), and enhance the space and time cross correlation between wall shear stress, pressure, and heat flux fluctuations. By extracting the envelopes of the traveling wave packets and inspecting the time and space correlations between the envelopes and the streaky structures, we found that the emergence of traveling wave packets comes later than the streaky structures, both in time and space. Based on these observations, we provide a depiction of the physical processes regarding the formation and evolution of the traveling wave packets.
Ming Yu, PengXin Liu, YaLu Fu, ZhiGong Tang, and XianXu Yuan
AIP Publishing
This two-part study investigates the effects of Mach number and wall temperature on the statistics of wall shear stress, pressure, and heat flux fluctuations in compressible wall-bounded turbulence. In the first part, we focus on their one-point statistics, including the root mean square (r.m.s.), skewness factor (third-order moment), flatness factor (fourth-order moment), and their correlations. By exploiting the direct numerical simulation databases, we found that the r.m.s. of the streamwise wall shear stress and pressure, the skewness factor of all the flow quantities considered, and the flatness factor of streamwise wall shear stress monotonically vary with the friction Mach number ([Formula: see text]), while for the rest, the wall heat flux and global temperature parameters should be taken into account as well for a monotonic trend of variation. The correlation coefficients between wall shear stress, pressure, and heat flux fluctuations increase with the Mach number [Formula: see text], suggesting the underlying interactions between dynamic and thermodynamic processes. The distributions of spectra and probability density functions indicate that the increased correlation is induced by the highly intermittent traveling wave packets among the streaky structures, as reflected by the “double-peak” feature of the spectra that gradually emerges with the increasing compressibility effects. The probability density distribution also manifests the alteration of the occurrence of extreme events caused by these structures. By accordingly decomposing the fluctuations with cutoff filtering, it is found that the root mean squares of streamwise wall shear stress and heat flux fluctuations related to the streaky structures are Mach number-independent, while those related to the traveling wave packets monotonically increase with the friction Mach number.
Siwei Dong, Fulin Tong, Ming Yu, Jianqiang Chen, Xianxu Yuan, and Qian Wang
AIP Publishing
In the present study, we investigate two-point statistics of fluctuating streamwise wall shear stress [Formula: see text] and wall heat flux [Formula: see text] by exploiting a direct numerical simulation database of supersonic turbulent boundary layers over a heated wall and a cooled wall at the friction Reynolds number around 800. By separately investigating positive and negative families of [Formula: see text] and [Formula: see text] with the aid of the conditional correlation analysis, we identify the asymmetrical deformation of [Formula: see text] and [Formula: see text], reminiscent of and ascribed to the asymmetrical deformations of sweeps and ejections events. The degree of such asymmetry is alleviated by the lower wall temperature. The spatial orientation of [Formula: see text] is insensitive to the wall temperature, whereas the spanwise elongated [Formula: see text] that is closely related to the wall pressure is manifested merely in the cooled-wall case. The cross correlation between [Formula: see text] and the fluctuating streamwise velocity u′ reveals that low-speed streaks related to negative [Formula: see text] are more inclined to the wall than high-speed ones related to positive [Formula: see text] by [Formula: see text], and that the phase lag between negative [Formula: see text] and u′ is larger than that between positive [Formula: see text] and u′ except in the near-wall region. Such a difference is proportional to the wall distance and should be considered for models predicting near-wall and wall quantities using signals in the logarithmic layer.
Ming Yu, ChunXiao Xu, JianQiang Chen, PengXin Liu, YaLu Fu, and XianXu Yuan
American Physical Society (APS)
Ming Yu and Chunxiao Xu
Cambridge University Press (CUP)
Predictive models for near-wall velocity and temperature fluctuations in compressible wall-bounded turbulence are developed in the present study based on the model proposed by Marusic et al. (Science, vol. 329 (5988), 2010, pp. 193–196), which incorporates the superposition and amplitude modulation effects of the large-scale motions in the outer region on near-wall turbulence. The density variation is involved in the predictive model for velocity fluctuations to achieve Mach number independence. The predictive model for temperature fluctuations is derived to keep its consistency with the strong Reynolds analogy, in which the modulation effect is supposed to be cast as the quadratic function of the large-scale velocity fluctuations. An algebraic method is proposed to directly determine the modulation coefficients and extract the universal signals. A direct numerical simulation (DNS) of turbulent channel flow at the friction Reynolds number of $1170$ and bulk Mach number of $2.88$ is carried out for parameter calibration and validations. The variances and joint probability density functions of the predicted velocity and temperature fluctuations agree well with the DNS results.
Dandan Yang, Yanfeng Gao, Ming Yu, Xiaoping Wen, and Ming-Xiang Zhao
AIP Publishing
Analysis of drag reduction effects due to axial oscillation of an inner cylinder in a turbulent Taylor–Couette (TC) flow is performed in the present study. The frictional Reynolds number on the inner cylinder is 218, and the non-dimensional oscillating period is varied from 8 to 32. By examining turbulence statistics, we uncover different impacts of the long- and short-period oscillations on the circumferential ( θ) and radial ( r) velocity fluctuations in large ([Formula: see text]) and small ([Formula: see text]) scales. One of the most surprising findings is that the short-period oscillation increases the large-scale Reynolds shear stress [Formula: see text] by the strong intensification of [Formula: see text] exceeding the suppression of [Formula: see text]. To understand the phenomena, the spectra of each term in the transport equations of the Reynolds normal stresses [Formula: see text] and [Formula: see text] are analyzed. First, it is shown that the short-period oscillation weakens the productions of [Formula: see text], and [Formula: see text] while it enhances that of [Formula: see text]. In contrast, the long-period oscillation reduces the productions of [Formula: see text] and [Formula: see text] while it mainly intensifies that of [Formula: see text]. Second, the investigations of the pressure–strain terms indicate that the short-period oscillation mainly impedes the inter-component energy transfer originating from the small-scale background turbulence. However, the long-period oscillation benefits the small-scale inter-component energy communication while it hinders the large-scale one. In addition, the inverse energy transfer in the turbulent TC flow is confirmed by inspecting the inter-scale energy transfer terms. The hindrance of the inter-scale energy transfer by the inner-cylinder oscillation plays a non-negligible role in the reduction of the wall friction drag.