Resolvent analysis provides a framework to predict coherent spatio-temporal structures of the largest linear energy amplification, through a singular value decomposition (SVD) of the resolvent operator, obtained by linearising the Navier–Stokes equations about a known turbulent mean velocity profile. Resolvent analysis utilizes a Fourier decomposition in time, which has thus far limited its application to statistically stationary or time-periodic flows. This work develops a variant of resolvent analysis applicable to time-evolving flows, and proposes a variant that identifies spatio-temporally sparse structures, applicable to either stationary or time-varying mean velocity profiles. Spatio-temporal resolvent analysis is formulated through the incorporation of the temporal dimension to the numerical domain via a discrete time-differentiation operator. Sparsity (which manifests in localisation) is achieved through the addition of an $l_1$-norm penalisation term to the optimisation associated with the SVD. This modified optimisation problem can be formulated as a nonlinear eigenproblem and solved via an inverse power method. We first showcase the implementation of the sparse analysis on a statistically stationary turbulent channel flow, and demonstrate that the sparse variant can identify aspects of the physics not directly evident from standard resolvent analysis. This is followed by applying the sparse space–time formulation on systems that are time varying: a time-periodic turbulent Stokes boundary layer and then a turbulent channel flow with a sudden imposition of a lateral pressure gradient, with the original streamwise pressure gradient unchanged. We present results demonstrating how the sparsity-promoting variant can either change the quantitative structure of the leading space–time modes to increase their sparsity, or identify entirely different linear amplification mechanisms compared with non-sparse resolvent analysis.

We report direct numerical simulations results of the rough-wall channel, focusing on roughness with high $k_{rms}/k_a$ statistics but small to negative $Sk$ statistics, and we study the implications of this new dataset on rough-wall modelling. Here, $k_{rms}$ is the root mean square, $k_a$ is the first-order moment of roughness height, and $Sk$ is the skewness. The effects of packing density, skewness and arrangement of roughness elements on mean streamwise velocity, equivalent roughness height ($z_0$) and Reynolds and dispersive stresses have been studied. We demonstrate that two-point correlation lengths of roughness height statistics play an important role in characterizing rough surfaces with identical moments of roughness height but different arrangements of roughness elements. Analysis of the present as well as historical data suggests that the task of rough-wall modelling is to identify geometric parameters that distinguish the rough surfaces within the calibration dataset. We demonstrate a novel feature selection procedure to determine these parameters. Further, since there is no finite set of roughness statistics that distinguish between all rough surfaces, we argue that obtaining a universal rough-wall model for making equivalent sand-grain roughness ($k_s$) predictions would be challenging, and that each rough-wall model would have its applicable range. This motivates the development of group-based rough-wall models. The applicability of multi-variate polynomial regression and feedforward neural networks for building such group-based rough-wall models using the selected features has been shown.

In hypersonic flight the shock wave and turbulent boundary layer interaction (STBLI) sharply increases wall heat transfer that intensifies the aerodynamic heating problems. In this work the STBLI is modelled by compression ramp flow with a Mach number of 5, a Reynolds number based on momentum thickness of 4652 and a wall to recovery temperature ratio of 0.5. The aerodynamic heat generation and transport mechanisms are investigated in the interaction based on theoretical analysis and direct numerical simulation (DNS) that agrees with previous studies. A prediction correlation of wall heat flux in STBLI is deduced theoretically and validated by some representative data including the present DNS, which improves the prediction accuracy and can be applied to a wider $Ma$ range compared with the canonical Q-P theory. The correlation indicates that the sharp increase of wall heat transfer in the STBLI can be explained by the boundary layer compression and the convection transport enhancement. Based on the DNS results, the aerodynamic heat generation and transport mechanisms are revealed in the separation, recirculation and reattachment zones in the STBLI. From this perspective, the peak heat flux can be further explained by the enhancement of near-wall turbulent energy dissipation, compression aerodynamic heat generation and the near-wall turbulent transport. The generation and transport of compression aerodynamic heat reveal the underlying mechanism of the strong correlation between the peak heat flux ratios and the pressure ratios in STBLIs.

Accurate prediction of the global properties of wall-bounded turbulence holds significant importance for both fundamental research and engineering applications. In atmospheric boundary layers, the relationship between friction drag and geostrophic wind is described by the geostrophic drag law (GDL). We use carefully designed large-eddy simulations to study nocturnal stable atmospheric boundary layers (NSBLs), which are characterized by a negative potential temperature flux at the surface and neutral stratification higher up. Our simulations explore a wider range of the Kazanski–Monin parameter, $\mu = L_f / L_s = [16.7, 193.3]$, with $L_f$ the Ekman length scale and $L_s$ the surface Obukhov length. We show collapse of the GDL coefficients onto single curves as functions of $\mu$, thereby validating the GDL's applicability to NSBLs over a very wide $\mu$ range. We show that the boundary-layer height $h$ scales with $\sqrt {L_f L_s}$, while both the streamwise and spanwise wind gradients scale with $u_*^2 / (h^2 f)$, where $u_*$ represents the friction velocity and $f$ the Coriolis parameter. Leveraging these insights, we developed new analytical expressions for the GDL coefficients, significantly enhancing our understanding of the GDL for turbulent boundary layers. These formulations facilitate the analytical prediction of the geostrophic drag coefficient and cross-isobaric angle.

This work introduces a formulation of resolvent analysis that uses wavelet transforms rather than Fourier transforms in time. Under this formulation, resolvent analysis may extend to turbulent flows with non-stationary mean states. The optimal resolvent modes are augmented with a temporal dimension and are able to encode the time-transient trajectories that are most amplified by the linearised Navier–Stokes equations. We first show that the wavelet- and Fourier-based resolvent analyses give equivalent results for statistically stationary flow by applying them to turbulent channel flow. We then use wavelet-based resolvent analysis to study the transient growth mechanism in the near-wall region of a turbulent channel flow by windowing the resolvent operator in time and frequency. The computed principal resolvent response mode, i.e. the velocity field optimally amplified by the linearised dynamics of the flow, exhibits characteristics of the Orr mechanism, which supports the claim that this mechanism is key to linear transient energy growth. We also apply this method to non-stationary parallel shear flows such as an oscillating boundary layer, and three-dimensional channel flow in which a sudden spanwise pressure gradient perturbs a fully developed turbulent channel flow. In both cases, wavelet-based resolvent analysis yields modes that are sensitive to the changing mean profile of the flow. For the oscillating boundary layer, wavelet-based resolvent analysis produces oscillating principal forcing and response modes that peak at times and wall-normal locations associated with high turbulent activity. For the turbulent channel flow under a sudden spanwise pressure gradient, the resolvent modes gradually realign themselves with the mean flow as the latter deviates. Wavelet-based resolvent analysis thus captures the changes in the transient linear growth mechanisms caused by a time-varying turbulent mean profile.

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the $-u^{\prime }w^{\prime }$ co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

This study investigates the correlation between the fluctuating wall heat flux, and the distribution and transport of Reynolds shear stress and turbulent heat flux in compressible boundary layers at Mach number 5.86 and friction Reynolds number 420, with a relatively weaker and a stronger wall cooling imposed. As illustrated from the probability density functions of the wall-heat-flux perturbations, with increasing wall cooling, the extreme wall heat flux is intensified and tends to be more negatively skewed. To examine the role of the extreme events in the transport of the momentum and heat, conditional analysis of the extreme positive and negative wall-heat-flux-perturbation events is conducted. In most regions of the boundary layer, the positive events are predominantly associated with an increase in Reynolds shear stress and a decrease in turbulent heat flux. Joint probability density functions of velocity and wall-heat-flux perturbations in the near-wall region indicate that the extreme positive events tend to be more correlated with ejections, which is particularly evident in the stronger wall-cooling case. To further shed light on the underlying mechanisms of the connections between wall heat flux and transport budgets, a transport equation for turbulent heat flux is derived, in a similar manner to that for Reynolds shear stress. The energy balance is inspected, with conditional analysis applied to budget terms and mean flow properties so as to quantify the correlation between wall-heat-flux fluctuations and energy evolution.

Turbulent flows over rough surfaces can be encountered in a wide range of engineering applications. Despite the progress made after several decades of studies, the prediction of drag and roughness function from the surface geometrical parameters remains an open question. Several methods have shown encouraging results. However, they lack generality and present some scatter in the data. In this paper we propose a new parameter, the effective distribution ($ED$), which lays foundation on the effective slope with some changes to take into account the sheltering effect of large roughness elements and the drag induced by pinnacles higher than the average roughness elements. To develop this new correlation between geometrical features of the wall and the drag, we performed a set of simulations of the turbulent flow over a rough surface made of triangular elements varying their height and spatial distribution. The $ED$ correlates quite well both with the drag and the roughness function for a wide range of cases having different mean roughness height, skewness and kurtosis. To further validate the $ED$, and assessing how it can be generalized to real rough wall, an irregular wall made from the superposition of random sinusoidal function was considered. Results were consistent with the correlation here presented.

Superhydrophobic (SHPo) surfaces can capture a thin layer of air called a plastron under water to reduce skin friction. Although a ~30 % drag reduction has been recently reported with longitudinal micro-trench SHPo surfaces under a boat and in a towing tank, the results lacked the consistency to establish a clear trend. Designed based on Yu et al. (J. Fluid Mech, vol. 962, 2023, A9), this work develops and tests a series of high-performance SHPo surface coupons that can sustain a pinned plastron underneath a passenger motorboat revamped to reach 14 knots. Importantly, plastrons in a pinned state, not just their existence, are confirmed during flow experiments for the first time. All the drag-reduction data measured on different longitudinal micro-trenches are found to collapse if plotted against slip length in wall units. In comparison, aligned posts and transverse trenches show less and little drag reduction, respectively, confirming the adverse effect of the spanwise slip in turbulent flows. This report not only verifies SHPo surfaces can provide a consistent drag reduction at high speeds in open sea but also shows that one may predict the amount of drag reduction in turbulent flows using the physical slip length obtained for Stokes flows.

In the present study, we propose a Reynolds analogy model for compressible wall turbulence. This model is demonstrated to be able to alleviate the defects of the generalized Reynolds analogy model (GRA) (Zhang et al., J. Fluid Mech., vol. 739, 2014, pp. 392–420), and maintains its success in describing the mean velocity–temperature relation. Furthermore, the present model is superior to the GRA in depicting the relationship between their fluctuating fields and also bridges the gap between the phenomenological model and the mathematical representation of the Reynolds analogy. The key points of the present model are validated by analysing the data of compressible wall-bounded turbulence with different Mach numbers, Reynolds numbers and wall thermal conditions.

Thwaites (Aeronaut. Q., vol. 1, 1949, pp. 245–280) developed an approximate method for determining the evolution of laminar boundary layers. The approximation follows from an assumption that the growth of a laminar boundary layer in the presence of pressure gradients could be parameterized solely as a function of the Holstein–Bohlen flow parameter, thus reducing the von Kármán momentum integral to a first-order ordinary differential equation. This method is useful for the analysis of laminar flows, and in computational potential flow solvers to account for the viscous effects. In this work, an approximate method for determining the momentum thickness of a two-dimensional, turbulent boundary layer is proposed following Thwaites’ work. It is shown that the method provides good estimates of the momentum thickness for multiple boundary layers, including both favourable and adverse pressure gradient effects, up to the point of separation. In the limit of high Reynolds numbers, it is possible to derive a criterion for the onset of separation from the proposed model, which is shown to be in agreement with prior empirical observations (Alber, 9th Aerospace Sciences Meeting, 1971). The sensitivity of the separation location with respect to upstream perturbations is also analysed through this model for the NASA/Boeing speed bump and the transonic Bachalo–Johnson bump.

Direct numerical simulations of oscillatory flow over a bed made of ripples have been performed. Two oscillatory flow forcing mechanisms have been compared: (i) a sinusoidal external pressure gradient (pressure-driven flow); and (ii) a sinusoidal velocity boundary conditions on the rippled bed (shear-driven flow). In the second case, the oscillations of the bed are such that when observed from a reference frame fixed with the bed, the free stream follows the same harmonic oscillation as in the pressure-driven case. While the outer layers have the same dynamics in the two cases, close to the bed differences are observed during the cycle, mostly because the large form drag across the ripples cannot be reproduced in the shear-driven case. A comparison against experimental data from an oscillating tray apparatus provides a relatively good agreement for the phase-averaged flow when the same forcing is considered (i.e. a shear-driven flow). The pressure-driven case has a comparable error to the shear-driven numerical results over the crest of the ripples, whereas the discrepancy is larger at the troughs. The discrepancies between the two cases are more limited for time-averaged flow quantities, such as the mean flow pattern and the time-averaged Reynolds stress distribution. This suggests that numerical or experimental shear-driven configurations may capture well the net effects of coastal transport processes (which occur in pressure-driven oscillatory flow), but care should be exercised in interpreting phase-dependent dynamics near the troughs. More work is needed to fully assess the sensitivity to the forcing mechanisms in different flow regimes.

This study aims to quantify how turbulence in a channel flow mixes momentum in the mean sense. We applied the macroscopic forcing method (Mani & Park, Phys. Rev. Fluids, 2021, 054607) to direct numerical simulation (DNS) of a turbulent channel flow at $Re_\tau =180$ using two different forcing strategies that are designed to separately assess the anisotropy and non-locality of momentum mixing. In the first strategy, the leading term of the Kramers–Moyal expansion of the eddy viscosity is quantified, revealing all 81 tensorial coefficients that essentially characterise the local-limit eddy viscosity. The results indicate the following: (1) the eddy viscosity has significant anisotropy, (2) Reynolds stresses are generated by both the mean strain rate and mean rotation rate tensors associated with the momentum field and (3) the local-limit eddy viscosity generates asymmetric Reynolds stress tensors. In the second strategy, the eddy viscosity is quantified as an integration kernel revealing the non-local influence of the mean momentum gradient at each wall-normal coordinate on all nine components of the Reynolds stresses over the channel width. Our results indicate that while the shear component of the Reynolds stress is reasonably reproduced by the local mean gradients, other components of the Reynolds stress are highly non-local. These results provide an understanding of anisotropy and non-locality requirements for closure modelling of momentum transport in attached wall-bounded turbulent flows.

Turbulent shear flows driven by a combination of a pressure gradient and buoyancy forcing are investigated using direct numerical simulations. Specifically, we consider the set-up of a differentially heated vertical channel subject to a Poiseuille-like horizontal pressure gradient. We explore the response of the system to its three control parameters: the Grashof number $Gr$, the Prandtl number $Pr$, and the Reynolds number $Re$ of the pressure-driven flow. From these input parameters, the relative strength of buoyancy driving to the pressure gradient can be quantified by the Richardson number $Ri=Gr/Re^2$. We compare the response of the mixed vertical convection configuration to that of mixed Rayleigh–Bénard convection, and find a nearly identical behaviour, including an increase in wall friction at higher $Gr$, and a drop in the heat flux relative to natural convection for $Ri=O(1)$. This closely matched response is despite vastly different flow structures in the systems. No large-scale organisation is visible in visualisations of mixed vertical convection – an observation that is confirmed quantitatively by spectral analysis. This analysis, combined with a statistical description of the wall heat flux, highlights how moderate shear suppresses the growth of small-scale plumes and reduces the likelihood of extreme events in the local wall heat flux. Vice versa, starting from a pure shear flow, the addition of thermal driving enhances the drag due to the emission of thermal plumes.

We investigate a turbulent boundary layer (TBL) under uniform pressure force variations, focusing on understanding its response to local pressure force, local pressure force variation (local disequilibrating effect) and upstream history. The flow starts as a zero-pressure-gradient (ZPG) TBL, followed by a uniform increase in the ratio of pressure force to turbulent force in the outer region and concludes with a uniform decrease of the same magnitude. This last zone includes a subzone with a diminishing adverse pressure gradient (APG), followed by an increasing favourable pressure gradient (FPG). In both subzones, the impact remains the same: mean momentum gain and turbulence reduction. In the outer region, the mean flow responds to force balance changes with a considerable delay. The accumulated flow history leads to a FPG TBL at the domain's end with a momentum defect comparable to APG TBLs. Below $y^+=10$, the mean flow responds almost instantaneously to pressure force changes. In the overlap layer, velocity profiles deviate from the conventional logarithmic law of ZPG TBL. Outer-layer turbulence decays more slowly than it increases initially, the latter turbulence increase persisting even after the pressure force begins to decrease. As a result of the slow turbulence decay, the FPG TBL at the domain's end exhibits unusually high outer turbulence levels. Near the wall, turbulence responds with a delay to pressure force changes, partly due to the influence of large-scale turbulence. All these significant cumulative effects of continuous pressure force variation indicate that parameters based solely on local variables cannot fully describe the physics of non-equilibrium TBLs.

High-fidelity simulations are conducted to investigate the turbulent boundary layers around a finite-span NACA0012 wing with a rounded wing-tip geometry at a chord-based Reynolds number of $Re_c=200\,000$ and at various angles of attack up to $10^\circ$. The study aims to discern the differences between the boundary layers on the finite-span wing and those on infinite-span wings at equivalent angles of attack. The finite-span boundary layers exhibit: (i) an altered streamwise and a non-zero spanwise pressure gradient as a result of the variable downwash induced by the wing-tip vortices (an inviscid effect typical of finite-span wings); (ii) differences in the flow history at different wall-normal distances, caused by the variable flow angle in the wall-normal direction (due to constant pressure gradients and variable momentum normal to the wall); (iii) laminar flow entrainment into the turbulent boundary layers near the wing tip (due to a laminar–turbulent interface); and (iv) variations in boundary layer thickness across the span, attributed to the variable wall-normal velocity in that direction (a primarily inviscid effect). These physical effects are then used to explain the differences in the Reynolds stress profiles and other boundary layer quantities, including the reduced near-wall peak of the streamwise Reynolds stress and the elevated Reynolds stress levels near the boundary layer edge, both observed in the finite-span wings. Other aspects of the flow, such as the downstream development of wing-tip vortices and their interactions with the surrounding flow, are reserved for future investigations.

The stirring and mixing of heat and momentum in the ocean surface boundary layer (OSBL) are dominated by 1 to 10 km fluid flows – too small to be resolved in global and regional ocean models. Instead, these processes are parametrized. Two main parametrizations include vertical mixing by surface-forced metre-scale turbulence and overturning by kilometre-scale submesoscale frontal flows and instabilities. In present models, these distinct parametrizations are implemented in tandem, yet ignore meaningful interactions between these two scales that may influence net turbulent fluxes. Using a large-eddy simulation of frontal spin down resolving processes at both scales, this work diagnoses submesoscale and surface-forced turbulence impacts that are the foundation of OSBL parametrizations, following a traditional understanding of these flows. It is shown that frontal circulations act to suppress the vertical buoyancy flux by surface forced turbulence, and that this suppression is not represented by traditional boundary layer turbulence theory. A main result of this work is that current OSBL parametrizations excessively mix buoyancy and overestimate turbulence dissipation rates in the presence of lateral flows. These interactions have a direct influence on the upper ocean potential vorticity and energy budgets with implications for global upper ocean budgets and circulation.

The present study investigates streamwise ($\overline {u^2}$) energy-transfer mechanisms in the inner and outer regions of turbulent boundary layers (TBLs). Particular focus is placed on the $\overline {u^2}$ production, its inter-component and wall-normal transport as well as dissipation, all of which become statistically significant in the outer region with increasing friction Reynolds number ($Re_{\tau }$). These properties are analysed using published data sets of zero, weak and moderately strong adverse-pressure-gradient (APG) TBLs across a decade of $Re_{\tau }$, revealing similarity in energy-transfer pathways for all these TBLs. It is found that both the inner and outer peaks of $\overline {u^2}$ are always associated with local maxima in the $\overline {u^2}$ production and its inter-component transport, and the regions below/above each of these peaks are always dominated by wall-ward/away-from-wall transport of $\overline {u^2}$, thereby classifying the $\overline {u^2}$ profiles into four distinct regimes. This classification reveals existence of phenomenologically similar energy-transfer mechanisms in the ‘inner’ and ‘outer’ regions of moderately strong APG TBLs, which meet at an intermediate location coinciding with the minimum in $\overline {u^2}$ profiles. Conditional averaging suggests existence of similar phenomena even in low $Re_{\tau }$ canonical and/or weak APG TBLs, albeit with the outer-region mechanisms weaker than those in the inner region. This explains the absence of their $\overline {u^2}$ outer peak and the dominance of $\overline {u^2}$ wall-normal transport away from the wall, which potentially originates from the inner region. Given that the wall-ward/away-from-wall transport of $\overline {u^2}$ is governed by the $Q_4$(sweeps)/$Q_2$(ejections) quadrants of the Reynolds shear stress, it is argued that the emergence of the $\overline {u^2}$ outer peak corresponds with the statistical dominance of $Q_4$ events in the outer region. Besides unravelling the dynamical significance of $Q_2$ and $Q_4$ events in the outer region of TBLs, the present analysis also proposes new phenomenological arguments for testing on canonical wall-turbulence data at very high $Re_{\tau }$.

This study uses high-fidelity simulations (direct numerical simulation or large-eddy simulation) and experimental datasets to analyse the effect of non-equilibrium streamwise mean pressure gradients (adverse or favourable), including attached and separated flows, on the statistics of boundary-layer wall-pressure fluctuations. The datasets collected span a wide range of Reynolds numbers ($Re_\theta$ from 300 to 23 400) and pressure gradients (Clauser parameter from $-0.5$ to 200). The datasets are used to identify an optimal set of variables to scale the wall-pressure spectrum: edge velocity, boundary layer thickness and the peak magnitude of Reynolds shear stress. Using the present datasets, existing semi-empirical models of the wall-pressure spectrum are shown unable to capture effects of strong, non-equilibrium adverse pressure gradients, due to inappropriate scaling of the wall pressure using the wall shear stress, calibration with limited types of flows and dependency on model parameters based on the friction velocity, which reduces to zero at the detachment point. To address these shortcomings, a generalized wall-pressure spectral model is developed with parameters that characterize the extent of the logarithmic layer and the strength of the wake. Derived from the local mean streamwise velocity profile, these two parameters inherently carry the effect of the Reynolds number, as well as those of the non-equilibrium pressure gradient and its history. Comparison with existing models shows that the proposed model behaves well and is more accurate in strong-pressure-gradient flows and in separated-flow regions.

As shown by Wenzel et al. (J. Fluid Mech., vol. 930, 2022, A1), the Eckert number $Ec$ defined using the difference between recovery temperature $\bar{T}_r$ and wall temperature $\bar{T}_w$ can be understood as a meaningful quantity to compare heat-transfer effects inside compressible turbulent boundary layers (for a calorically perfect gas), no matter whether these are caused by different Mach-number or wall-temperature conditions. While the named study deduced this comparative behaviour of $Ec$ from an integral perspective in a strict sense, Cogo et al. (J. Fluid Mech., vol. 974, 2023, A10) performed a systematic parameter study based on the previous findings to look at wall-normal profiles. They have shown that the diabatic parameter $\varTheta$, being equivalent to $Ec$, is capable of categorizing heat-transfer effects for cases at different Mach numbers, even to some extent for some of the wall-normal profiles. Building on this progress, the present paper provides a comprehensive classification of both existing and newly computed super- and hypersonic direct numerical simulation data at various wall temperature conditions into heated cases, adiabatic cases or weakly/moderately/strongly/quasi-incompressibly cooled cases. Hereby, the classification is largely based on the wall-normal position of the temperature peak occurring in cooled boundary-layer cases, which is one of the determining factors for the topological characteristics of diabatic boundary-layer profiles. Integrating high-enthalpy data into the analysis allowed us to confirm the reliability of the proposed classification also in more complex scenarios, where the calorically perfect gas assumption no longer applies and additional heat-transfer mechanisms come into play. While the Eckert number is shown to well characterize heat-transfer effects on most important temperature-related quantities for a wide range of Mach numbers and $\bar {T}_w/\bar {T}_r$ conditions, also the local Reynolds number $Re_{\tau }$ is shown to notably affect the strength of heat-transfer effects. Since both $Ec$ and $Re_{\tau }$ can be determined in advance – or estimated to a reasonable extent – a key advantage of the classification scheme is to allow for an effective a priori estimation of the extent to which heat-transfer effects are to be expected for a given compressible turbulent boundary-layer configuration.